Development of mussel mRNA profiling: Can gene expression trends reveal coastal water pollution?

Hypoxia in the Blue Mussel Mytilus chilensis Induces a Transcriptome Shift Associated with Endoplasmic Reticulum Stress, Metabolism, and Immune Response

The increase in hypoxia events, a result of climate change in coastal and fjord ecosystems, impacts the health and survival of mussels. These organisms deploy physiological and molecular responses as an adaptive mechanism to maintain cellular homeostasis under environmental stress. However, the specific effects of hypoxia on mussels of socioeconomic interest, such as Mytilus chilensis, are unknown. Using RNA-seq, we investigated the transcriptomic profiles of the gills, digestive gland, and adductor muscle of M. chilensis under hypoxia (10 days at 2 mg L-1) and reoxygenation (10 days at 6 mg L-1). There were 15,056 differentially expressed transcripts identified in gills, 11,864 in the digestive gland, and 9862 in the adductor muscle. The response varied among tissues, showing chromosomal changes in Chr1, Chr9, and Chr10 during hypoxia. Hypoxia regulated signaling genes in the Toll-like, mTOR, citrate cycle, and apoptosis pathways in gills, indicating metabolic and immunological alterations. These changes suggest that hypoxia induced a metabolic shift in mussels, reducing reliance on aerobic respiration and increasing reliance on anaerobic metabolism. Furthermore, hypoxia appeared to suppress the immune response, potentially increasing disease susceptibility, with negative implications for the mussel culture industry and natural bed populations. This study provides pivotal insights into metabolic and immunological adaptations to hypoxia in M. chilensis, offering candidate genes for adaptive traits.

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.

Antidepressants and their metabolites primarily affect lysosomal functions in the marine mussel, Mytilus galloprovincialis

Antidepressants widely occur as emerging contaminants in marine coastal waters, with concentrations reported in the low ng/L range. Although at relatively lower levels with respect to other pharmaceuticals, antidepressants - fluoxetine (FLX) in particular - have attracted attention because of their striking effects exerted at low doses on marine invertebrates. In this study, the effects of four antidepressants including FLX, sertraline (SER), and citalopram, as members of the selective serotonin reuptake inhibitor (SSRI) class, and venlafaxine (VEN) as a member of the serotonin and norepinephrine reuptake inhibitor (SNRI) class, were evaluated in the mussel Mytilus galloprovincialis. In addition, the effects of two main metabolites of FLX and VEN, i.e., norfluoxetine (NFL) and O-desmethylvenlafaxine (ODV) respectively, were compared to those of the parent compounds. Eight concentrations of each drug (0.5-500 ng/L range) were tested on the early life stage endpoints of gamete fertilization and larval development at 48 h post fertilization (hpf). Egg fertilization was reduced by all compounds, except for VEN. Larval development at 48 hpf was affected by all SSRIs, but not by SNRIs. The above effects were significant but never exceeded 20 % of control values. Adult mussels were exposed in vivo for 7 days to environmental concentrations of the drugs (0.5, 5, and 10 ng/L) and a battery of eight biomarkers was assessed. Antidepressants primarily targeted lysosomal functions, decreasing haemocyte lysosome membrane stability (up to 70 % reduction) and increasing of the lysosome/cytosol ratio (up to 220 %), neutral lipid (up to 230 %), and lipofuscin (up to 440 %) accumulation in digestive gland. Only SER and NFL significantly affected catalase and glutathione-S-transferase activities in gills and digestive gland. NFL and ODV, were effective and sometimes more active than the parent compounds. All compounds impaired mussel health status, as indicated by the low to high stress levels assigned using the Mussel Expert System.

Effect of Different Species of Prorocentrum Genus on the Japanese Oyster Crassostrea gigas Proteomic Profile

This paper assesses the effects of exposure to toxic concentrations (1200 to 6000 cells/mL) of the dinoflagellates Prorocentrum lima, Prorocentrum minimum, and Prorocentrum rhathymum and several concentrations of aqueous and organic extracts obtained from the same species (0 to 20 parts per thousand) on the Crassostrea gigas (5-7 mm) proteomic profile. Through comparative proteomic map analyses, several protein spots were detected with different expression levels, of which eight were selected to be identified by liquid chromatography-mass spectrometry (LC-MS/MS) analyses. The proteomic response suggests that, after 72 h of exposure to whole cells, the biological functions of C. gigas affected proteins in the immune system, stress response, contractile systems and cytoskeletal activities. The exposure to organic and aqueous extracts mainly showed effects on protein expressions in muscle contraction and cytoskeleton morphology. These results enrich the knowledge on early bivalve developmental stages. Therefore, they may be considered a solid base for new bioassays and/or generation of specific analytical tools that allow for some of the main effects of algal proliferation phenomena on bivalve mollusk development to be monitored, characterized and elucidated.

Single-molecule long-read (SMRT) transcriptome sequencing of Mercenaria mercenaria reveals a powerful anti-apoptotic system critical for air exposure endurance

Mercenaria mercenaria is an economically important clam species and exhibits an outstanding resistance to multiple environmental stressors. However, our understanding of their stress adaptability is limited due to a lack of genomic information, such as transcriptome resources. In this study, single-molecule long-read (SMRT) mRNA sequencing was performed to obtain the full-length gill transcriptome reference sequences of M. mercenaria under air exposure stress. In all, 14.5 G subreads were obtained and assembled into 64,603 unigenes, among which 50,613 were successfully annotated. Additionally, 56,295 SSRs, 1457 transcription factors, and 5924 lncRNAs were identified in M. mercenaria transcriptome. Furthermore, numerous apoptosis-related transcripts were identified according to Swiss-Prot annotation and their numbers were counted. We also found that most apoptosis-related transcripts exhibited typical domains of a certain protein family through conserved domain prediction. Additionally, eight typical sequences related to apoptosis pathway were detected by RT-PCR, with the aim to show the sequential variation of gene expression levels under air exposure. These results implied that the complicated apoptosis system, especially the powerful anti-apoptotic system was critical for M. mercenaria to endure air exposure.

Effects of Sudden Drop in Salinity on Osmotic Pressure Regulation and Antioxidant Defense Mechanism of Scapharca subcrenata

Salinity is an important ecological factor that impacts the growth and survival of aquatic organisms. The salinity of seawater in coastal and estuarine areas is often subject to dynamic changes because of seasonal rainfall and continental runoff. Thus, the current study investigated the effects of sudden changes in salinity on the survival rate and osmotic pressure regulation mechanisms of bottom-sowing seedlings of the economically important ark shell, Scapharca subcrenata. By simulating the sudden changes that occur in seawater salinity after rainstorms, the results showed that the osmotic pressure of the hemolymph and Na⁺, K⁺, Ca²⁺, and Cl^– concentrations first decreased and then increased. When the salinity decreased from 30 to 14‰, hemoglobin, soluble total protein, taurine, and total free amino acid gradually increased; maximum levels of hemoglobin, soluble total protein, and taurine occurred once the salinity increased to 22‰ at 96 h. After 96 h, the total free amino acid content increased until 144 h. The reactive oxygen species (ROS) content and total antioxidant capacity (T-AOC) peaked at 96 h, whereas the expression levels of Mn-superoxide dismutase (MnSOD) and catalase (CAT) increased earlier, indicating that, with continuous ROS generation, antioxidant defense mechanisms were activated to avoid oxidative damage. Expression levels of cathepsin C (CTSC), cathepsin D (CTSD), heat shock protein 20 (HSP20), and heat shock protein 70 (HSP70) were significantly higher than in the control group at 48 h (salinity level 14‰); the expression levels of HSP20, heat shock protein 90 (HSP90), MnSOD, and glutathione peroxidase (GP_x) remained high, indicating that they were still required for osmotic pressure regulation to maintain the dynamic balance between the generation and removal of ROS as the salinity level increased. These results not only add to our basic understanding of the aquatic ecology of S. subcrenata, but also provide a theoretical ground for improving the survival rate of bottom-sowing, propagation, and release of S. subcrenata seedlings.

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.

Dietary nano cerium oxide promotes growth, relieves ammonia nitrogen stress, and improves immunity in crab (Eriocheir sinensis)

Oxidative stress plays a crucial role in ammonia nitrogen toxicity. In this study, the beneficial effects of dietary nano cerium oxide (nano CeO₂) as a potent antioxidant were examined in the Chinese mitten crab (Eriocheir sinensis). Crabs were fed a diet supplemented with 0, 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, or 12.8 mg/kg nano CeO₂ for 60 d. The optimum supplementation level of nano CeO₂ that significantly increased weight gain rate and decreased feed coefficient was 0.8 mg/kg. This level also offered immune protection when crabs were kept under ammonia nitrogen stress and/or exposed to pathogen infection (Aeromonas hydrophila). Supplementation with 0.8 mg/kg of CeO₂ (i) relieved pathological damage to the hepatopancreas; (ii) increased hemocyte counts, including total number of hemocytes, granulocytes, and hyalinocytes; (iii) decreased malondialdehyde content and increased antioxidant enzyme activities of superoxide dismutase and catalase in the hemolymph; (iv) increased the activities of lysozyme, acid phosphatase, and alkaline phosphatase in the hemolymph; and (v) increased gene and protein expression of cathepsin L in the hepatopancreas. Mortality increased when crabs were injected with bacteria under ammonia nitrogen stress, but dietary supplementation with 0.8 mg/kg nano CeO₂ decreased the mortality rate. Thus, the results of this study suggested that dietary supplementation with nano CeO₂ in crabs promoted growth and up-regulated immunity to bacterial infection under ammonia nitrogen stress.

Mixed evidence for adaptation to environmental pollution

Adaptation to pollution has been studied since the first observations of heavy metal tolerance in plants decades ago. To document micro-evolutionary responses to pollution, researchers have used phenotypic, molecular genetics, and demographic approaches. We reviewed 258 articles and evaluated the evidence for adaptive responses following exposure to a wide range of pollutants, across multiple taxonomic groups. We also conducted a meta-analysis to calculate the magnitude of phenotypic change in invertebrates in response to metal pollution. The majority of studies that reported differences in responses to pollution were focused on phenotypic responses at the individual level. Most of the studies that used demographic assays in their investigations found that negative effects induced by pollution often worsened over multiple generations. Our meta-analysis did not reveal a significant relationship between metal pollution intensity and changes in the traits studied, and this was probably due to differences in coping responses among different species, the broad array of abiotic and biotic factors, and the weak statistical power of the analysis. We found it difficult to make broad statements about how likely or how common adaptation is in the presence of environmental contamination. Ecological and evolutionary responses to contamination are complex, and difficult to interpret in the context of taxonomic, and methodological biases, and the inconsistent set of approaches that have been used to study adaptation to pollution in the laboratory and in the field. This review emphasizes the need for: (a) long-term monitoring programs on exposed populations that link demography to phenotypic, genetic, and selection assays; (b) the use of standardized protocols across studies especially for similar taxa. Approaches that combine field and laboratory studies offer the greatest opportunity to reveal the complex eco-evolutionary feedback that can occur under selection imposed by pollution.

Genomics and immunity of the Mediterranean mussel Mytilus galloprovincialis in a changing environment

The Mediterranean mussel (Mytilus galloprovincialis) is a marine invasive species cultured all over the world. Mussels are an appreciated resource in local aquaculture enterprises because of their robust production and resilience that translates into a reliable economic value. So far, no massive mortalities have been reported in natural or cultured populations of this species. In the last years, the knowledge about its immune system has greatly improved but there are still many questions to be answered. One of them is why mussels, with their high filtering activity, are able to be exposed to a high number of potential pathogens without getting infected and without developing an elevated inflammatory response. The sequencing of the mussel genome has revealed a very complex organization with high heterozygosity, abundance of repetitive sequences and extreme intraspecific sequence diversity among individuals, mainly in immune related genes. Among those genes, antimicrobial peptides are the most expressed gene families in mussels, highly polymorphic and with antimicrobial effect against molluscs pathogens, but also against pathogens of lower vertebrates and humans. The combination of a complex genome with the adaptation of mussel immune system to a changing environment could explain this high variability, not only in immune-related genes, but also in the functional response among individuals sampled in the same location and date.

Revealing Mytilus galloprovincialis transcriptomic profiles during ontogeny

Mediterranean mussels are a worldwide spread bivalve species with extraordinary biological success. One of the reasons of this success could be the reproduction strategy of bivalves, characterized by the presence of trochophore larvae. Larval development in bivalves has been a topic of raising interest in the scientific community but it deserves much more attention. The principal objective of this work was to study the transcriptomic profile of the ontogeny of Mytilus galloprovincialis analyzing the gene expression in different developmental stages, from oocytes to juveniles. For this purpose, after conducting a 454 sequencing of the transcriptomes of mussel hemocytes, adult tissues and larvae, a new DNA microarray was designed and developed. The studied developmental stages: unfertilized oocytes, veliger, pediveliger, settled larvae and juveniles, showed very different transcriptomic profiles and clustered in groups defining their characteristic gene expression along ontogeny. Our results show that oocytes present a distinct and characteristic transcriptome. After metamorphosis, both settled larvae and juveniles showed a very similar transcriptome, with no enriched GO terms found between these two stages. This suggests: 1.- the progressive loss of RNA of maternal origin through larval development and 2.- the stabilization of the gene expression after settlement. On the other hand during metamorphosis a specific profile of differentially expressed genes was found. These genes were related to processes such as differentiation and biosynthesis. Processes related to the immune response were strongly down regulated. These suggest a development commitment at the expense of other non-essential functions, which are temporary set aside. Immune genes such as antimicrobial peptides suffer a decreased expression during metamorphosis. In fact, we found that the oocytes which express a higher quantity of genes such as myticins are more likely to reach success of the offspring, compared to oocytes poor in such mRNAs, whose progeny died before reaching metamorphosis.

Synthesis methods influence characteristics, behaviour and toxicity of bare CuO NPs compared to bulk CuO and ionic Cu after in vitro exposure of Ruditapes philippinarum hemocytes

Copper oxide (CuO) nanoparticles (NPs) are increasingly investigated, developed and produced for a wide range of industrial and consumer products. Notwithstanding their promising novel applications, concern has been raised that their increased use and disposal could consequently increase their release into marine systems and potentially affect species within. To date the understanding of factors and mechanisms of CuO (nano-) toxicity to marine invertebrates is still limited. Hence, we studied the characteristics and behaviour of two commercially available CuO NPs of similar size, but produced employing distinct synthesis methods, under various environmentally and experimentally relevant conditions. In addition, cell viability and DNA damage, as well as gene expression of detoxification, oxidative stress, inflammatory response, DNA damage repair and cell death mediator markers were studied in primary cultures of hemocytes from the marine clam Ruditapes philippinarum and, where applicable, compared to bulk CuO and ionic Cu (as CuSO₄) behaviour and effects. We found that the synthesis method can influence particle characteristics and behaviour, as well as the toxicity of CuO NPs to Ruditapes philippinarum hemocytes. Our results further indicate that under the tested conditions aggregating behaviour influences the toxicity of CuO NPs by influencing their rate of extra- and intracellular dissolution. In addition, gene expression analysis identified similar transcriptional de-regulation for all tested copper treatments for the here measured suite of genes. Finally, our work highlights various differences in the aggregation and dissolution kinetics of CuO particles under environmental (marine) and cell culture exposure conditions that need consideration when extrapolating in vitro findings.

Early ecotoxic effects of ZnO nanoparticle chronic exposure in Mytilus galloprovincialis revealed by transcription of apoptosis and antioxidant-related genes

Recently, China became one of the largest nanomaterial markets in the world. The wide use of ZnO nanoparticles in a number of products implies an increasing release in marine environment and consequently the evaluation of the potential effects upon marine organisms largely cultured in China for commercial purposes, such as invertebrate bivalves is a current need. To this aim, survival, bioaccumulation, and transcription pattern of key genes, p53, PDRP, SOD, CAT, and GST, involved in DNA damage/repair and antioxidation, in Mytilus galloprovincialis digestive gland, exposed to ZnO NPs (<100 nm) and ZnO bulk (150-200 nm) for 4 weeks, were evaluated. ZnSO₄ was also assessed to appraise the role of zinc ions. Starting from 72 h, increasing mortality values along the exposure time were observed for all ZnO compounds. The highest difference was evident after 28 d when NPs resulted three times more toxic than bulk, (LC₅₀) = 0.78 mg Zn/L (confidence limits: 0.64, 1.00) and 2.62 mg Zn/L (confidence limits: 1.00, 4.00), respectively. For ZnSO₄ the (LC₅₀) was always the lowest reaching the minimum value at 28 d 0.25 mg Zn/L (confidence limits: 0.10-0.40). Digestive gland showed higher uptake rate of ionic Zn respect to ZnO NPs and bulk during the first three days of exposure. In particular at the end of the exposure time (28 d) at 1 mg Zn/L the rank of Zn uptake rate was Zinc ion > ZnO NPs > ZnO bulk. The relative expression of investigated genes evidenced that distinct actions of apoptosis and antioxidation occurred in M. galloprovincialis exposed to ZnO NPs with a peculiar pattern dependent on exposure time and concentration. Application of the qRT-PCR technique revealed evidence of sensitivity to the nanomaterial since the first time of exposure.

Interactive effects of nutrition, reproductive state and pollution on molecular stress responses of mussels, Mytilus galloprovincialis Lamarck, 1819

Marine bivalves including mussels Mytilus galloprovincialis are commonly used as sentinels for pollution monitoring and ecosystem health assessment in the coastal zones. Use of biomarkers to assess the pollution effects assumes that the effects of pollutants on the biomarkers exceed the natural background variability; yet this assumption has rarely been tested. We exposed mussels at different reproductive stages and nutritive states to two concentrations of a polycyclic aromatic hydrocarbon (fluoranthene, 3 and 60 μg L^-1) for three weeks. Expression levels of the molecular biomarkers related to the detoxification and general stress response [cytochrome P450 oxidase (CYP450), glutathione S-transferases (GST-α; GST-S1; GST-S2), the multixenobiotic resistance protein P-glycoprotein (PgP), metallothioneins (MT10 and MT20), heat shock proteins (HSP22, HSP70-2; HSP70-3; HSP70-4), as well as mRNA expression of two reproduction-related genes, vitellogenin (Vitel) and vitelline coat lysin M7 (VCLM7)] were measured. The mussels' nutrition and reproductive state affected the baseline mRNA levels of molecular biomarkers and modulated the transcriptional responses of biomarker genes to the pollutant exposure. Thus, mussel physiological state could act as a confounding factor in the evaluation of the response of pollution through molecular biomarkers. The biomarker baseline levels must be determined across a range of physiological states to enable the use of biomarkers in monitoring programs.

Functional and molecular responses of the blue mussel Mytilus edulis' hemocytes exposed to cadmium - An in vitro model and transcriptomic approach

The bivalve mollusk, Mytilus edulis, is used as a sentinel species in several monitoring programs due to its ability to bio-accumulate contaminants. Its immune system consists of hemocytes and humoral components, which constitute the main part of the hemolymph. The present study is aimed at understanding the effects of Cd on the differentially expressed genes involved in the phagocytosis of M. edulis' hemocytes. Our approach focuses on an in vitro model by exposing hemocytes to different concentrations of Cd ranging from 10^-9 M to 10^-3 M. Phagocytosis and cell viability as functional markers were measured using flow cytometry. The molecular mechanisms regulated by Cd were investigated using RNA-seq and DGE analysis. Results showed that viability and phagocytosis of hemocytes exposed to 10^-3 M of Cd were significantly decreased after 21 h of exposure. RNA sequencing data showed that 1112 transcripts (out of 352,976 contigs) were differentially regulated by the highest concentration of Cd. Among these identified transcripts, 1028 and 84 were up and down-regulated respectively. The induction of super oxide dismutase (SOD), glutathion-s-transferase (GST), cytochrome P450 2C8 (CYP2C8), multidrug resistance protein (MRP1) and heat shock protein 70 (HSP70) suggests that Cd can regulate key molecular mechanisms. In addition, several toll-like receptors (TLR) as well as genes involved in phagocytosis (actin and CDC42) and apoptosis (caspase 8 and XIAP/IAP) were induced by Cd. Thus, our model highlights the effect of Cd on the phagocytic function of M. edulis' hemocytes along with the regulation of gene expression involved in innate immunity, detoxification and apoptosis. Further investigations need to be pursued to unravel the effects of Cd on the molecular mechanisms identified in this study.

Assessing the impact of Benzo[a]pyrene on Marine Mussels: Application of a novel targeted low density microarray complementing classical biomarker responses

Despite the increasing use of mussels in environmental monitoring and ecotoxicological studies, their genomes and gene functions have not been thoroughly explored. Several cDNA microarrays were recently proposed for Mytilus spp., but putatively identified partial transcripts have rendered the generation of robust transcriptional responses difficult in terms of pathway identification. We developed a new low density oligonucleotide microarray with 465 probes covering the same number of genes. Target genes were selected to cover most of the well-known biological processes in the stress response documented over the last decade in bivalve species at the cellular and tissue levels. Our new ‘STressREsponse Microarray’ (STREM) platform consists of eight sub-arrays with three replicates for each target in each sub-array. To assess the potential use of the new array, we tested the effect of the ubiquitous environmental pollutant benzo[a]pyrene (B[a]P) at 5, 50, and 100 μg/L on two target tissues, the gills and digestive gland, of Mytilus galloprovincialis exposed invivo for three days. Bioaccumulation of B[a]P was also determined demonstrating exposure in both tissues. In addition to the well-known effects of B[a]P on DNA metabolism and oxidative stress, the new array data provided clues about the implication of other biological processes, such as cytoskeleton, immune response, adhesion to substrate, and mitochondrial activities. Transcriptional data were confirmed using qRT-PCR. We further investigated cellular functions and possible alterations related to biological processes highlighted by the microarray data using oxidative stress biomarkers (Lipofuscin content) and the assessment of genotoxicity. DNA damage, as measured by the alkaline comet assay, increased as a function of dose.DNA adducts measurements using ³²P-postlabeling method also showed the presence of bulky DNA adducts (i.e. dG-N²-BPDE). Lipofiscin content increased significantly in B[a]P exposed mussels. Immunohistochemical analysis of tubulin and actin showed changes in cytoskeleton organisation. Our results adopting an integrated approach confirmed that the combination of newly developed transcriptomic approcah, classical biomarkers along with chemical analysis of water and tissue samples should be considered for environmental bioimonitoring and ecotoxicological studies to obtain holistic information to assess the impact of contaminants on the biota.

Transcriptomic effects of the non-steroidal anti-inflammatory drug Ibuprofen in the marine bivalve Mytilus galloprovincialis Lam

The transcriptomic effects of Ibuprofen (IBU) in the digestive gland tissue of Mytilus galloprovincialis Lam. specimens exposed at low environmental concentrations (250 ng L(-1)) are presented. Using a 1.7 K feature cDNA microarray along with linear models and empirical Bayes statistical methods 225 differentially expressed genes were identified in mussels treated with IBU across a 15-day period. Transcriptional dynamics were typical of an adaptive response with a peak of gene expression change at day-7 (177 features, representing about 11% of sequences available for analysis) and an almost full recovery at the end of the exposure period. Functional genomics by means of Gene Ontology term analysis unraveled typical mussel stress responses i.e. aminoglycan (chitin) metabolic processes but also more specific effects such as the regulation of NF-κB transcription factor activity.

Combined effects of n-TiO2 and 2,3,7,8-TCDD in Mytilus galloprovincialis digestive gland: A transcriptomic and immunohistochemical study

Despite the growing concern over the potential biological impact of nanoparticles (NPs) in the aquatic environment, little is known about their interactions with other pollutants. In the marine mussel Mytilus galloprovincialis, exposure to nanosized titanium dioxide (n-TiO2), one of the most widespread type of NPs in use, in combination with and 2,3,7,8-tetrachlorodibenzo-p-dioxins (TCDD), chosen as model organic xenobiotic, was shown to induce significant changes in different biomarkers in hemocytes, gills and digestive gland, with distinct effects depending on cell/tissue and type of response measured. In this work, the interactive effects of n-TiO2 and TCDD at the tissue level were further investigated in mussel digestive gland using an integrated approach transcriptomics/immunohistochemistry. Mussels were exposed to n-TiO2 (100μgL(-1)) and TCDD (0.25μgL(-1)), alone and in combination, for 96h. Transcriptomic analysis identified 48-, 49- and 62 Differentially Expressed Genes (DEGs) in response to n-TiO2, TCDD and n-TiO2/TCDD, respectively. Gene Ontology (GO) term analysis revealed distinct biological processes affected in different experimental conditions. n-TiO2 mainly up-regulated cytoskeletal genes, while TCDD up-regulated endocrine and signal transduction related processes. Co-exposure induced transcriptional changes common to individual treatments, and identified a newly generated process, response to chemical stimulus. Transcription of selected genes was verified by qPCR. Moreover, expression of tubulin, as an example of target protein of interest identified by gene transcription data, was evaluated in tissue sections by immunolabelling. Tissue TCDD accumulation was evaluated by immunofluorescence with an anti-dioxins antibody. The results demonstrate both distinct and interactive effects of n-TiO2 and TCDD in mussel digestive gland at the molecular and tissue level, identify the main molecular targets involved, and underline how exposure to the n-TiO2/TCDD mixture does not result in increased TCDD accumulation and overall stressful conditions in the tissue. These represent the first data on transcriptional responses of marine invertebrates to exposure not only to n-TiO2 as a model of NP, but also to a legacy contaminant like TCDD.

Cathepsin L is an immune-related protein in Pacific abalone (Haliotis discus hannai)--Purification and characterization

Cathepsin L, an immune-related protein, was purified from the hepatopancreas of Pacific abalone (Haliotis discus hannai) by ammonium sulfate precipitation and column chromatographies of SP-Sepharose and Sephacryl S-200 HR. Purified cathepsin L appeared as two bands with molecular masses of 28.0 and 28.5 kDa (namely cathepsin La and Lb) on SDS-PAGE under reducing conditions, suggesting that it is a glycoprotein. Peptide mass fingerprinting (PMF) analysis revealed that peptide fragments of 95 amino acid residues was high similarity to cathepsin L of pearl oyster (Pinctada fucata). The optimal temperature and pH of cathepsin L were 35 °C and pH 5.5. Cathepsin L was particularly inhibited by cysteine proteinase inhibitors of E-64 and leupeptin, while it was activated by metalloproteinase inhibitors EDTA and EGTA. The full-length cathepsin L cDNA was further cloned from the hepatopancreas by rapid PCR amplification of cDNA ends (RACE). The open reading frame of the enzyme was 981 bp, encoding 327 amino acid residues, with a conserved catalytic triad (Cys134, His273 and Asn293), a potential N-glycosylation site and conserved ERFNIN, GNYD, and GCGG motifs, which are characteristics of cathepsin L. Western blot and proteinase activity analysis revealed that the expression and enzyme activity of cathepsin L were significantly up-regulated in hepatopancreas at 8 h following Vibrio parahaemolyticus infection, demonstrating that cathepsin L is involved in the innate immune system of abalone. Our present study for the first time reported the purification, characterization, molecular cloning, and tissue expression of cathepsin L in abalone.

A multibiomarker approach to explore interactive effects of propranolol and fluoxetine in marine mussels

A multi-biomarker approach, including several lysosomal parameters, activity and mRNA expression of antioxidant enzymes, and DNA damage, was employed to investigate the nominal effects of 0.3 ng/L fluoxetine (FX) and 0.3 ng/L propranolol (PROP) alone or in combination (0.3 ng/L FX + 0.3 ng/L PROP) on Mediterranean mussels after a 7 day treatment. FX co-exposure appears to facilitate PROP bioaccumulation because PROP only accumulated in digestive gland of FX + PROP treated mussels. Lysosomal parameters were significantly impaired by FX + PROP treatment, while no clear antioxidant responses at the catalytic and transcriptional levels were observed. Biomarker responses led to a "medium stress level" diagnosis in FX + PROP treated mussels, according to the Expert System, whereas 0.3 ng/L PROP or FX alone did not induce consistent stress conditions. These findings suggest vulnerability of coastal marine mussels to FX and PROP contamination at environmentally relevant levels.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Does zebra mussel (Dreissena polymorpha) represent the freshwater counterpart of Mytilus in ecotoxicological studies? A critical review

One of the fundamentals in the ecotoxicological studies is the need of data comparison, which can be easily reached with the help of a standardized biological model. In this context, any biological model has been still proposed for the biomonitoring and risk evaluation of freshwaters until now. The aim of this review is to illustrate the ecotoxicological studies carried out with the zebra mussel Dreissena polymorpha in order to suggest this bivalve species as possible reference organism for inland waters. In detail,we showed its application in biomonitoring, as well as for the evaluation of adverse effects induced by several pollutants, using both in vitro and in vivo experiments. We discussed the advantages by the use of D. polymorpha for ecotoxicological studies, but also the possible limitations due to its invasive nature.

Effect of seawater desalination and oil pollution on the lipid composition of blue mussels Mytilus edulis L. from the White Sea

A study on the effect oil pollution under normal and reduced salinity had on blue mussels Mytilus edulis L. from the White Sea in an aquarium-based experiment and in the natural habitat revealed a change in gill total lipids as a compensatory response. The cholesterol concentration and the cholesterol/phospholipids ratio in gills were found to reflect the impact of the environmental factors (oil pollution and desalination), and evidence adaptive changes in the cell membrane structure. An elevated content of storage lipids (chiefly triacylglycerols) in the mussels in the aquarium experiment indicates, first of all, the uptake and accumulation of oil products in gill cells under both normal and reduced seawater salinity, while high triacylglycerols level in gill littoral mussels from 'control' biotope in the Gulf of Kandalaksha is primarily associated with the mussel׳s pre-spawning period.

Seasonal variations of gene expression biomarkers in Mytilus galloprovincialis cultured populations: temperature, oxidative stress and reproductive cycle as major modulators

The blue mussel Mytilus galloprovincialis has been used as monitoring organism in many biomonitoring programs because of its broad distribution in South European sea waters and its physiological characteristics. Different pollution-stress biomarkers, including gene expression biomarkers, have been developed to determine its physiological response to the presence of different pollutants. However, the existing information about basal expression profiles is very limited, as very few biomarker-based studies were designed to reflect the natural seasonal variations. In the present study, we analyzed the natural expression patterns of several genes commonly used in biomonitoring, namely ferritin, metallothionein, cytochrome P450, glutathione S-transferase, heat shock protein and the kinase responsive to stress KRS, during an annual life cycle. Analysis of mantle-gonad samples of cultured populations of M. galloprovincialis from the Delta del Ebro (North East Spain) showed natural seasonal variability of these biomarkers, pointing to temperature and oxidative stress as major abiotic modulators. In turn, the reproductive cycle, a process that can be tracked by VCLM7 expression, and known to be influenced by temperature, seems to be the major biotic factor involved in seasonality. Our results illustrate the influence of environmental factors in the physiology of mussels through their annual cycle, a crucial information for the correct interpretation of responses under stress conditions.

Molecular cloning and characterization of cathepsin L from freshwater mussel, Cristaria plicata

Cathepsin L is one of the crucial enzyme superfamilies and involved in the immune responses. The Cathepsin L cDNA and genome of Cristaria plicata（CpCL） was cloned from the hemocytes using degenerate primers by the rapid amplification of cDNA ends (RACE) PCR. The genomic DNA was 9353 bp long and had a total of six introns and seven exons. The full-length cDNA of CpCL was 1144 bp, the cDNA contained a 5' untranslated region (UTR) of 34 nucleotides, the 3' UTR of 108 bp with a canonical polyadenylation signal sequence AATAAA and a polyA tail, and an open reading frame (ORF) of 1002 bp, encoding 333 amino acid residues with 37.65 kDa predicted molecular weight. The theoretical isoelectric point was 8.61. The prepro-cathepsin L was consisted of a typical signal peptide (Met1-Gly20), a pro-region peptide (Leu21-Glu116) and a mature peptide (Tyr117-Val333). Many members of the papain family possessed of a proline residue at position 2 in the mature enzymem, this was also observed in CpCL. The preproprotein included an oxyanion hole (Gln 135), the active center formed by Cys141, His280 and Asn 300, the potential N-glycosylation site (Asn38, Asn 113 and Asn 272) and the conserved GCXGG motifs, which was characteristic of cathepsin, the conserved ERWNIN and GNFD motifs, which were characteristic for cathepsin L. Homology analysis revealed that the CpCL shared 49-87% identity to other known cathepsin L sequences. The phylogenetic tree showed that the CpCL clustered with the invertebrate cathepsin L cysteine proteases, and was closely related to the cathepsin L of Hyriopsis cumingii. The expression of CpCL mRNA was detected in hepatopancreas, hemocytes, mantle, gills and adductor muscle, and the higher expression level was in hepatopancreas. After A. hydrophila stimulation, the expression of the CpCL mRNA was up-regulated in hemocytes and hepatopancreas, and the expression level was significantly lower in gill than one after PBS challenge group.

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.

Proteomic response of mussels Mytilus galloprovincialis exposed to CuO NPs and Cu²⁺: an exploratory biomarker discovery

CuO NPs are one of the most used metal nanomaterials nowadays with several industrial and other commercial applications. Nevertheless, less is known about the mechanisms by which these NPs inflict toxicity in mussels and to what extent it differs from Cu(2+). The aim of this study was to investigate changes in protein expression profiles in mussels Mytilus galloprovincialis exposed for 15 days to CuO NPs and Cu(2+) (10 μg L(-1)) using a proteomic approach. Results demonstrate that CuO NPs and Cu(2+) induced major changes in protein expression in mussels' showing several tissue and metal-dependent responses. CuO NPs showed a higher tendency to up-regulate proteins in the gills and down-regulate in the digestive gland, while Cu(2+) showed the opposite tendency. Distinctive sets of differentially expressed proteins were found, either common or specific to each Cu form and tissue, reflecting different mechanisms involved in their toxicity. Fifteen of the differentially expressed proteins from both tissues were identified by MALDI-TOF-TOF. Identified proteins indicate common response mechanisms induced by CuO NPs and Cu(2+), namely in cytoskeleton and cell structure (actin, α-tubulin, paramyosin), stress response (heat shock cognate 71, putative C1q domain containing protein), transcription regulation (zinc-finger BED domain-containing protein 1, nuclear receptor subfamily 1G) and energy metabolism (ATP synthase F0 subunit 6). CuO NPs alone also had a marked effect on other biological processes, namely oxidative stress (GST), proteolysis (cathepsin L) and apoptosis (caspase 3/7-1). On the other hand, Cu(2+) affected a protein associated with adhesion and mobility, precollagen-D that is associated with the detoxification mechanism of Cu(2+). Protein identification clearly showed that the toxicity of CuO NPs is not solely due to Cu(2+) dissolution and can result in mitochondrial and nucleus stress-induced cell signalling cascades that can lead to apoptosis. While the absence of the mussel genome precluded the identification of other proteins relevant to clarify the effects of CuO NPs in mussels' tissues, proteomics analysis provided additional knowledge of their potential effects at the protein level that after confirmation and validation can be used as putative new biomarkers in nanotoxicology.

Identification, expression, and responses to bacterial challenge of the cathepsin C gene from the razor clam Sinonovacula constricta

Cathepsin C (dipeptidyl-peptidaseI, DPPI) is a lysosomal cysteine proteinase that belongs to the papain superfamily, and it is involved in protein degradation and proenzyme activation. However, very little is known about the function of cathepsin C in bivalves. In the present study, we identified the cathepsin C gene in the razor clam Sinonovacula constricta (Sc-CTSC). The full-length Sc-CTSC cDNA contained a complete open reading frame (ORF) of 1371 nt encoding 456 amino acids, a 98 bp 5' UTR, and a 1043 bp 3' UTR. The ORF of Sc-CTSC consisted of a putative signal peptide of 22 aa, a propeptide of 229 aa, and a mature peptide of 205 aa containing the active site triad of Cys, His, and Asn. The Sc-CTSC transcript was expressed in a wide range of tissues but exhibited the greatest level of expression in the digestive gland. During the early developmental stages, the transcript was detected widely. Upon injection with Vibrio anguillarum, the Sc-CTSC transcript was significantly up-regulated in digestive gland, mantle, and gill tissues. The results provided important information for further exploring the roles of cathepsin C in the innate immune responses.

A selection of reference genes and early-warning mRNA biomarkers for environmental monitoring using Mytilus spp. as sentinel species

mRNA biomarkers are promising tools for environmental health assessment and reference genes are needed to perform relevant qPCR analyses in tissue samples of sentinel species. In the present study, potential reference genes and mRNA biomarkers were tested in the gills and digestive glands of native and caged mussels (Mytilus spp.) exposed to harbor pollution. Results highlighted the difficulty to find stable reference genes in wild, non-model species and suggested the use of normalization indices instead of single genes as they exhibit a higher stability. Several target genes were found differentially expressed between mussel groups, especially in gills where cyp32, π-gst and CuZn-sod mRNA levels could be biomarker candidates. Multivariate analyses confirmed the ability of mRNA levels to highlight site-effects and suggested the use of several combined markers instead of individual ones. These findings support the use of qPCR technology and mRNA levels as early-warning biomarkers in marine monitoring programs.

The role of environmental biotechnology in exploring, exploiting, monitoring, preserving, protecting and decontaminating the marine environment

In light of the Marine Strategy Framework Directive (MSFD) and the EU Thematic Strategy on the Sustainable Use of Natural Resources, environmental biotechnology could make significant contributions in the exploitation of marine resources and addressing key marine environmental problems. In this paper 14 propositions are presented focusing on (i) the contamination of the marine environment, and more particularly how to optimize the use of biotechnology-related tools and strategies for predicting and monitoring contamination and developing mitigation measures; (ii) the exploitation of the marine biological and genetic resources to progress with the sustainable, eco-compatible use of the maritime space (issues are very diversified and include, for example, waste treatment and recycling, anti-biofouling agents; bio-plastics); (iii) environmental/marine biotechnology as a driver for a sustainable economic growth.

Transcriptomic responses to heat stress and nickel in the mussel Mytilus galloprovincialis

The exposure of marine organisms to stressing agents may affect the level and pattern of gene expression. Although many studies have examined the ecological effects of heat stress on mussels, little is known about the physiological mechanisms that maybe affected by co-exposure to heat stress and environmental contaminants such as nickel (Ni). In the present work, we investigated the effects of simultaneous changes in temperature and Ni supply on lysosomal membrane stability (LMS) and malondialdehyde accumulation (MDA) in the digestive gland (DG) of the blue mussel Mytilus galloprovincialis (Lam.). To elucidate how the molecular response to environmental stressors is modulated, we employed a cDNA microarray with 1673 sequences to measure relative transcript abundances in the DG of mussels exposed to Ni along with a temperature increase. A two-way ANOVA revealed that temperature and Ni rendered additive effects on LMS and MDA accumulation, increasing the toxic effects of metal cations. Ni loads in the DG were also affected by co-exposure to 26°C. In animals exposed only to heat stress, functional genomics analysis of the microarray data (171 differentially expressed genes (DEGs)) highlighted seven biological processes, largely dominated by the up-regulation of folding protein-related genes and the down-regulation of genes involved in cell migration and cellular component assembly. Exposure to Ni at 18°C and 26°C yielded 188 and 262 DEGs, respectively, exhibiting distinct patterns in terms of biological processes. In particular, the response of mussels exposed to Ni at 26°C was characterized by the up-regulation of proteolysis, ribosome biogenesis, response to unfolded proteins, and catabolic-related genes, as well as the down-regulation of genes encoding cellular metabolic processes. Our data provide new insights into the transcriptomic response in mussels experiencing temperature increases and Ni exposure; these data should be carefully considered in view of the biological effects of heat stress, particularly in polluted areas.

Development of a decision support system to manage contamination in marine ecosystems

In recent years, contamination and its interaction with climate-change variables have been recognized as critical stressors in coastal areas, emphasizing the need for a standardized framework encompassing chemical and biological data into risk indices to support decision-making. We therefore developed an innovative, expert decision support system (Exp-DSS) for the management of contamination in marine coastal ecosystems. The Exp-DSS has two main applications: (i) to determine environmental risk and biological vulnerability in contaminated sites; and (ii) to support the management of waters and sediments by assessing the risk due to the exposure of biota to these matrices. The Exp-DSS evaluates chemical data, both as single compounds and as total toxic pressure of the mixture, to compare concentrations to effect-based thresholds (TELs and PELs). Sites are then placed into three categories of contamination: uncontaminated, mildly contaminated, and highly contaminated. In highly contaminated sites, effects on high-level ecotoxicological endpoints (i.e. survival and reproduction) are used to determine risk at the organism-population level, while ecological parameters (i.e. alterations in community structure and ecosystem functions) are considered for assessing effects on biodiversity. Changes in sublethal biomarkers are utilized to assess the stress level of the organisms in mildly contaminated sites. In Triad studies, chemical concentrations, ecotoxicological high-level effects, and ecological data are combined to determine the level of environmental risk in highly contaminated sites; chemical concentration and ecotoxicological sublethal effects are evaluated to determine biological vulnerability in mildly contaminated sites. The Exp-DSS was applied to data from the literature about sediment quality in estuarine areas of Spain, and ranked risks related to exposure to contaminated sediments from high risk (Huelva estuary) to mild risk (Guadalquivir estuary and Bay of Cadiz). A spreadsheet-based version of the Exp-DSS is available at the MEECE and DiSIT web sites (www.meece.eu and www.disit.unipmn.it).

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.

Gene discovery in the tropical scallop Nodipecten nodosus: construction and sequencing of a normalized cDNA library

We report the construction and characterization of a normalized cDNA library from the digestive gland of the marine bivalve Nodipecten nodosus, a commercially valuable tropical scallop. A total of 288 clones were sequenced, and 250 unique sequences were obtained. The cDNA library showed a small sequence redundancy (2.3%) and high numbers of recombinant (99.9%) and independent clones (2.0 × 10(6) cfu), indicating that the cDNA library generated in this study is a profitable resource for efficient gene discovery for N. nodosus. EST functional annotation by Gene Ontology term assignment revealed the identification of sequences potentially involved in aquaculture and ecotoxicology relevant processes such as apoptosis, growth, lipid metabolism, reproduction, development, response to stress and immunity.

Molecular characterization and expression analysis of four cathepsin L genes in the razor clam, Sinonovacula constricta

Cathepsin L (CTSL) is a lysosomal cysteine protease involved in immune responses in vertebrates. However, few studies exist regarding the role of cathepsin L in bivalves. In this study, we isolated and characterized four cathepsin L genes from the razor clam Sinonovacula constricta, referred to as CTSL1, CTSL2, CTSL3 and CTSL4. These four genes contained typical papain-like cysteine protease structure and enzyme activity sites with ERWNIN-like and GNFD-like motifs in the proregion domain and an oxyanion hole (Gln) and a catalytic triad (Cys, His and Asn) in the mature domain. Expression analysis of the four transcripts revealed a tissue-specific pattern with high expression of CTSL1 and CTSL3 in liver and gonad tissues and high expression of CTSL2 and CTSL4 in liver and gill tissues. During the developmental stages, the four transcripts showed the highest expression in the juvenile stage; however, CTSL3 had a much higher expression level than the other three transcripts during embryogenesis. The four transcripts showed significant changes in expression as early as 4 h or 8 h after infection with Vibrio anguillarum. The fact that bacterial infection can induce expression of the four CTSL transcripts suggests that these transcripts are important components of the innate immunity system of the clam.

Transcriptional response of the mussel Mytilus galloprovincialis (Lam.) following exposure to heat stress and copper

Global warming is a major factor that may affect biological organization, especially in marine ecosystems and in coastal areas that are particularly subject to anthropogenic pollution. We evaluated the effects of simultaneous changes in temperature and copper concentrations on lysosomal membrane stability (N-acetyl-hexosaminidase activity) and malondialdehyde accumulation (MDA) in the gill of the blue mussel Mytilus galloprovincialis (Lam.). Temperature and copper exerted additive effects on lysosomal membrane stability, exacerbating the toxic effects of metal cations present in non-physiological concentrations. Mussel lysosomal membrane stability is known to be positively related to scope for growth, indicating possible effects of increasing temperature on mussel populations in metal-polluted areas. To clarify the molecular response to environmental stressors, we used a cDNA microarray with 1,673 sequences to measure the relative transcript abundances in the gills of mussels exposed to copper (40 µg/L) and a temperature gradient (16°C, 20°C, and 24°C). In animals exposed only to heat stress, hierarchical clustering of the microarray data revealed three main clusters, which were largely dominated by down-regulation of translation-related differentially expressed genes, drastic up-regulation of protein folding related genes, and genes involved in chitin metabolism. The response of mussels exposed to copper at 24°C was characterized by an opposite pattern of the genes involved in translation, most of which were up-regulated, as well as the down-regulation of genes encoding heat shock proteins and "microtubule-based movement" proteins. Our data provide novel information on the transcriptomic modulations in mussels facing temperature increases and high copper concentrations; these data highlight the risk of marine life exposed to toxic chemicals in the presence of temperature increases due to climate change.

Transcriptomic seasonal variations in a natural population of zebra mussel (Dreissena polymorpha)

The zebra mussel Dreissena polymorpha is a Caspian Sea bivalve that colonized freshwater bodies worldwide during the XX century. To analyze the impact of seasonal and environmental variations on the physiology and metabolism of this invasive species, we developed a custom microarray using 4057 publicly available DNA sequences from Dreissena and other related genera. Transcriptome profiles were analyzed using half-body samples from a relatively clean site (Riba-Roja, low Ebro River, N.E. Spain), at three different stages of the annual cycle: Pre-spawning (February), spawning (June), and gonad resorption (September). Transcripts from a total of 745 unique sequences showed significant changes among these three groups of samples. Functional characterization of these transcripts based on their closest known homologues showed that genes involved in stress defense (oxidative and infection) were overrepresented in September, whereas genes related to reproductive functions were overrepresented in the spawning and pre-spawning periods. This transcriptomic information can help to identify developmental stages at which the organism is more vulnerable for future control strategies. These data will also contribute to the implementation of gene expression-based assays for pollution monitoring in water bodies harboring stable zebra mussel populations.

Genes of the mitochondrial apoptotic pathway in Mytilus galloprovincialis

Bivalves play vital roles in marine, brackish, freshwater and terrestrial habitats. In recent years, these ecosystems have become affected through anthropogenic activities. The ecological success of marine bivalves is based on the ability to modify their physiological functions in response to environmental changes. One of the most important mechanisms involved in adaptive responses to environmental and biological stresses is apoptosis, which has been scarcely studied in mollusks, although the final consequence of this process, DNA fragmentation, has been frequently used for pollution monitoring. Environmental stressors induce apoptosis in molluscan cells via an intrinsic pathway. Many of the proteins involved in vertebrate apoptosis have been recognized in model invertebrates; however, this process might not be universally conserved. Mytilus galloprovincialis is presented here as a new model to study the linkage between molecular mechanisms that mediate apoptosis and marine bivalve ecological adaptations. Therefore, it is strictly necessary to identify the key elements involved in bivalve apoptosis. In the present study, six mitochondrial apoptotic-related genes were characterized, and their gene expression profiles following UV irradiation were evaluated. This is the first step for the development of potential biomarkers to assess the biological responses of marine organisms to stress. The results confirmed that apoptosis and, more specifically, the expression of the genes involved in this process can be used to assess the biological responses of marine organisms to stress.

Identification and expression of multiple CYP1-like and CYP3-like genes in the bivalve mollusk Mytilus edulis

Various sequencing projects over the last several years have aided the discovery of previously uncharacterized invertebrate sequences, including new cytochrome P450 genes (CYPs). Here we present data on the identification and characterization of two CYP1-like and three CYP3-like genes from the bivalve mollusk Mytilus edulis, and assess their potential as biomarkers based on their responses to several known vertebrate aryl hydrocarbon receptor (AHR) agonists. Quantitative real-time PCR was used to measure CYP transcript levels in digestive gland, labial palps, adductor muscle, gill, foot, and different regions of the mantle. Levels of both CYP1-like genes were highest in digestive gland, whereas labial palps had the highest expression levels of the three CYP3-like genes followed by digestive gland and outer margin of the mantle. Mussels were exposed by injection to the AHR agonists, β-naphthoflavone (BNF; 25 μg g(-1)), 3,3',4,4',5-polychlorinated biphenyl (PCB126; 2 μg g(-1)), or 6-formylindolo[3,2-b]carbazole (FICZ; 0.1 μg g(-1)), or to Aroclor 1254 (a mixture of PCBs; 50 μg g(-1)) for 24 h, followed by CYP expression analysis. There was no statistically significant change in expression of either of the CYP1-like genes after exposure to the various AHR agonists. The CYP3-like-1 gene was significantly up-regulated by BNF in gill tissues and the CYP3-like-2 gene was up-regulated in digestive gland by PCB126 and in gill tissue by BNF. These results suggest that distinct mechanisms of CYP gene activation could be present in M. edulis, although the importance of the CYP1-like and CYP3-like genes for xenobiotic and endogenous lipids biotransformation requires additional investigation.

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.

Transcriptomic analysis of Ruditapes philippinarum hemocytes reveals cytoskeleton disruption after in vitro Vibrio tapetis challenge

The Manila clam, Ruditapes philippinarum, is an economically-important, commercial shellfish; harvests are diminished in some European waters by a pathogenic bacterium, Vibrio tapetis, that causes Brown Ring disease. To identify molecular characteristics associated with susceptibility or resistance to Brown Ring disease, Suppression Subtractive Hybridization (SSH) analyzes were performed to construct cDNA libraries enriched in up- or down-regulated transcripts from clam immune cells, hemocytes, after a 3-h in vitro challenge with cultured V. tapetis. Nine hundred and ninety eight sequences from the two libraries were sequenced, and an in silico analysis identified 235 unique genes. BLAST and "Gene ontology" classification analyzes revealed that 60.4% of the Expressed Sequence Tags (ESTs) have high similarities with genes involved in various physiological functions, such as immunity, apoptosis and cytoskeleton organization; whereas, 39.6% remain unidentified. From the 235 unique genes, we selected 22 candidates based upon physiological function and redundancy in the libraries. Then, Real-Time PCR analysis identified 3 genes related to cytoskeleton organization showing significant variation in expression attributable to V. tapetis exposure. Disruption in regulation of these genes is consistent with the etiologic agent of Brown Ring disease in Manila clams.

In-depth proteomic analysis of a mollusc shell: acid-soluble and acid-insoluble matrix of the limpet Lottia gigantea

BACKGROUND

Invertebrate biominerals are characterized by their extraordinary functionality and physical properties, such as strength, stiffness and toughness that by far exceed those of the pure mineral component of such composites. This is attributed to the organic matrix, secreted by specialized cells, which pervades and envelops the mineral crystals. Despite the obvious importance of the protein fraction of the organic matrix, only few in-depth proteomic studies have been performed due to the lack of comprehensive protein sequence databases. The recent public release of the gastropod Lottia gigantea genome sequence and the associated protein sequence database provides for the first time the opportunity to do a state-of-the-art proteomic in-depth analysis of the organic matrix of a mollusc shell.

RESULTS

Using three different sodium hypochlorite washing protocols before shell demineralization, a total of 569 proteins were identified in Lottia gigantea shell matrix. Of these, 311 were assembled in a consensus proteome comprising identifications contained in all proteomes irrespective of shell cleaning procedure. Some of these proteins were similar in amino acid sequence, amino acid composition, or domain structure to proteins identified previously in different bivalve or gastropod shells, such as BMSP, dermatopontin, nacrein, perlustrin, perlucin, or Pif. In addition there were dozens of previously uncharacterized proteins, many containing repeated short linear motifs or homorepeats. Such proteins may play a role in shell matrix construction or control of mineralization processes.

CONCLUSIONS

The organic matrix of Lottia gigantea shells is a complex mixture of proteins comprising possible homologs of some previously characterized mollusc shell proteins, but also many novel proteins with a possible function in biomineralization as framework building blocks or as regulatory components. We hope that this data set, the most comprehensive available at present, will provide a platform for the further exploration of biomineralization processes in molluscs.

Suppressive subtractive hybridization libraries prepared from the digestive gland of the oyster Crassostrea brasiliana exposed to a diesel fuel water-accommodated fraction

Diesel fuel can cause adverse effects in marine invertebrates by mechanisms that are not clearly understood. The authors used suppressive subtractive hybridization to identify genes up- and downregulated in Crassostrea brasiliana exposed to diesel fuel. Genes putatively involved in protein regulation, innate immune, and stress response, were altered by diesel challenge. Three genes regulated by diesel were validated by quantitative real-time polymerase chain reaction. This study sheds light on transcriptomic responses of oysters to diesel pollution.

Proteomic discovery of biomarkers of metal contamination in Sydney Rock oysters (Saccostrea glomerata)

In the current study we examined the effects of metal contamination on the protein complement of Sydney Rock oysters. Saccostrea glomerata were exposed for 4 days to three environmentally relevant concentrations (100 μg/l, 50 μg/l and 5 μg/l) of cadmium, copper, lead and zinc. Protein abundances in oyster haemolymph from metal-exposed oysters were compared to those from non-exposed controls using two-dimensional electrophoresis to display differentially expressed proteins. Differentially expressed proteins were subsequently identified using tandem mass spectrometry (LC-MS/MS), to assign their putative biological functions. Unique sets of differentially expressed proteins were affected by each metal, in addition to proteins that were affected by more than one metal. The proteins identified included some that are commonly associated with environmental monitoring, such as HSP 70, and other novel proteins not previously considered as candidates for molecular biomonitoring. The most common biological functions of proteins were associated with stress response, cytoskeletal activity and protein synthesis.