Interactions of a pesticide/heavy metal mixture in marine bivalves: a transcriptomic assessment
© Dondero et al; licensee BioMed Central Ltd. 2011
Received: 27 April 2010
Accepted: 16 April 2011
Published: 16 April 2011
Skip to main content
© Dondero et al; licensee BioMed Central Ltd. 2011
Received: 27 April 2010
Accepted: 16 April 2011
Published: 16 April 2011
Mixtures of chemicals present in aquatic environments may elicit toxicity due to additive or synergistic effects among the constituents or, vice versa, the adverse outcome may be reduced by antagonistic interactions. Deviations from additivity should be explained either by the perturbations of toxicokinetic parameters and/or chemical toxicodynamics. We addressed this important question in marine mussels exposed subchronically to a binary mixture made of two wide-spread pollutants: the heavy metal nickel and the organic phosphorus pesticide Chlorpyrifos. To this aim, we carried out in tissues of Mytius galloprovincialis (Lam) a systems approach based on the evaluation and integration of different disciplines, i.e. high throughput gene expression profiling, functional genomics, stress biomakers and toxicokinetics.
Cellular and tissue biomarkers, viz. digestive gland lysosomal membrane stability, lysosomal/cytosol volume ratio, neutral lipid content and gill acetylcholinesterase activity were, in general, altered by either the exposure to nickel and Chlorpyrifos. However, their joint action rendered (i) an overall decrease of the stress syndrome level, as evaluated through an expert system integrating biomarkers and (ii) statistically significant antagonistic deviations from the reference model systems to predict mixture toxicity. While toxicokinetic modeling did not explain mixture interactions, gene expression profiling and further Gene Ontology-based functional genomics analysis provided clues that the decrement of toxicity may arise from the development of specific toxicodynamics. Multivariate statistics of microarray data (238 genes in total, representing about 14% of the whole microarray catalogue) showed two separate patterns for the single chemicals: the one belonging to the heavy metal -135 differentially expressed genes (DEGs) was characterized by the modulation of transcript levels involved in nucleic acid metabolism, cell proliferation and lipid metabolic processes. Chlorpyrifos exposure (43 DEGs) yielded a molecular signature which was biased towards carbohydrate catabolism (indeed, chitin metabolism) and developmental processes. The exposure to the mixture (103 DEGs) elicited a composite complex profile which encompassed the core properties of the pesticide but also a relevant set of unique features. Finally, the relative mRNA abundance of twelve genes was followed by Q-PCR to either confirm or complement microarray data. These results, in general, were compatible with those from arrays and indeed confirmed the association of the relative abundance of two GM-2 ganglioside activator genes in the development of the hyperlipidosis syndrome observed in digestive gland lysosomes of single chemical exposed mussels.
The transcriptomic assessment fitted with biological data to indicate the occurrence of different toxicodynamic events and, in general, a decrease of toxicity, driven by the mitigation or even abolition of lysosomal responses. Furthermore, our results emphasized the importance of the application of mechanistic approaches and the power of systems assessment to study toxicological responses in ecologically relevant organisms.
Pollutants are present in aquatic environments in the form of complex mixtures whose single compounds may be not toxic, per se, at the considered concentrations [1, 2]. Although scientists generally have a good understanding of the toxicity of individual chemical pollutants, there is a great need to bridge the gap between our understanding of the toxic effects of exposure to individual xenobiotics and those effects from exposure to mixtures of such chemicals . The toxic effects of mixtures are usually predicted from reference models based on non-interaction among the single chemicals, i.e. concentration addition  and response addition  have been established. Several studies dealing in particular with classical ecotoxicological endpoints such as mortality and reproduction, have suggested the validity of both models as a first screening to analyze the toxicity of mixtures . However a considerable part of mixture data (20-40%) is not describable by such models, indicating a significant interaction among chemicals, leading to synergistic or antagonistic outcomes . Even if limited to binary mixture assessments, recently improved mathematical models have been proposed to describe such deviations . However, the great challenge still remains of predicting such occurrences and providing mechanistic explanations on mixture toxicity.
Mussels have been extensively used in biomonitoring projects through the application of a battery of physiological and cellular biological responses able to prove the occurrence of a stress syndrome and the biological risk associated with polluted environments [8–12]. Recently, furthermore, transcriptomics approaches have been successfully applied to these species to unveil the molecular mechanisms of adaptation to both natural and chemical stressors [13–15]. With the advent of the post-genomic and second generation sequencing era, ecotoxicologists have enthusiastically embraced trancriptomic profiling as a tool to assess exposure to environmental stressors . However, gene expression may not represent a direct marker of functional responses as only gene products and metabolites do represent the final cellular effectors Recently, systems biology -a combination of high-throughput molecular disciplines- was proposed to improve the landscape of protective and non-protective responses occurring in cell/tissues of a given bioindicator species. This specific application is known as systems toxicology and is defined also as the study of perturbations of biological systems by chemicals and stressors as well as monitoring changes in molecular expression and conventional toxicological parameters . A successful example is represented by the study of the non-model organism L. rubellus, from which was obtained a strong inference about the mechanistic effects of copper on earthworms .
A first attempt on marine bivalves was proposed by  who showed the benefit of the integration of biomarkers with gene expression changes to interpret the physiological status of specimens collected along a copper pollution gradient. Notwithstanding, that work was limited by the assessment of a few stress genes. More recently, our research group has developed a high density (1.7 K) cDNA array and used it to assess distinct molecular fingerprints in the tissues of mussels exposed either to model contaminants in laboratory conditions or collected from field sites challenged with different pollutant levels .
We wondered whether transcriptomics and, in particular, a systems toxicology assessment may effectively represent a valuable approach for studying the biological effects of a mixture of pollutants and specifically to predict the interactions among its members, as recently proposed . Here we present results based on transcriptomics, physiological biomarkers (lysosomal responses) and toxicokinetics measurements obtained in the digestive gland tissue of the marine bivalve Mytilus galloprovincialis exposed to a binary mixture of chemicals for four days. The toxic agents were selected as dissimilarly-acting toxicants from a panel of priority environmental pollutants within the framework of the European Project NoMiRACLE : Chlorpyrifos, a broad-spectrum organophosphate insecticide and nickel (Ni), an important heavy metal for its civil and industrial applications. Chlorpyrifos represents one of the most utilized phyto-pharmacological products in the world for both crop protection and pest control  It's present in the marine environment, even including biota, of the American coast since the early 1990's . Like other organophosphorus compounds, its mode of action is mainly based on the inhibition of the acetylcholinesterase system, both in vertebrate and invertebrate models [22–24]. In addition, several studies identified putative neurodevelopmental mechanisms that are independent of cholinesterase inhibition [25–28]. CHP has been shown to interfere with different components of cell signalling [29, 30] and to affect oxidative stress parameters in the developing brain, leading to shifts in expression and function of antioxidant genes [31, 32]. On the other hand, nickel is considered a dangerous pollutant, in particular for the recognized carcinogenic activity probably related to the production of oxidative DNA damage and the inhibition of DNA repair activity. Furthermore, nickel can produce oxidative stress that depletes glutathione, activates Ap1, NF-kB and other oxidatively sensitive transcription factors. Nevertheless, exact knowledge of the molecular mechanisms of nickel toxicity and carcinogenicity is still limited . Recently, a few studies reporting the ecotoxicological effects due to the joint Ni and Chlorpyrifos exposure in different model species were published as results of the NoMiRACLE project [34–37]. Here we report our original findings on mixture toxicity analysis using a systems toxicology approach in the marine mussel Mytilus galloprovincialis.
EC12.5 (0.25 TU)
EC25 (0.5 TU)
EC50 (1.0 TU)
Statistical testing of antagonistic interactions
Fit of CA
CA vs A
A vs DL
Fit of IA
IA vs A
A vs DL
Real time Q-PCR analysis: comparison with microarray data
metallothionein isoform mt-10b
metallothionein isoform mt-20-IV
gm2 ganglioside activator protein
gm2 ganglioside activator protein
Match (absolute; relative)
Common gene frequency
Unique Mix including opposite
Real-time quantitative PCR was carried out to confirm microarray data and refine the relative expression levels of a selection of paralogue genes such as four chitinases, two metallothionein genes (mt10, mt20) and two ganglioside GM2 activator proteins, the latter involved in lysosomal lipid metabolism. Moreover, another two genes implicated in lipid metabolic processes were selected to confirm the relevance of this pathway in exposed bivalves: hexosaminidase and apolipophorin precursor. An actin variant and the p53-like protein gene were also included in this survey. Microarray and Q-PCR analysis showed consistent outcomes in 27 out 33 (82%) comparisons made on the three classes of toxic treatments, but chitinase and GM2-AP genes, in some cases, could not be confirmed (Table 4). Considering the high degree of identity of such sequences (data not shown), it is likely that microarray probes could not provide a reliable assessment and therefore Q-PCR outcome was further considered for the discussion of data. In mixture exposed tissues, Q-PCR analysis could also confirm the up-regulation of the hexosaminidase gene which from microarray showed a positive, but not statistically significant, expression level value.
The basic concept for the description of a toxicological action of components of a mixture is based on the principle of non-interaction, which means that chemicals in the mixture do not affect the toxicity potency of one another or each other's probability to exert a toxic effect . Recent works focusing on the cumulative toxicity of Ni and CHP reported the occurrence of interactive effects at different levels of biological organization. An antagonistic interaction between the two chemicals on the locomotor activity of the zebrafish larvae was reported , while authors in  suggested a synergistic effect on the survival of the ground beetle Pterostichus oblongopunctatus indicating species-specific responses. In mussels, metabolomic profiling of digestive gland tissues obtained from animals exposed to Ni and CHP was compatible with a reduction in toxicity , arguing, therefore, antagonistic interactions. In the current study, we present for the first time in mussels, a systems biology assessment integrating the toxicokinetics and toxicodynamics of these chemicals and their mixture, thereby obtaining clues on the molecular mechanisms as the basis of pollutants interactions in mussel tissues.
The exposure of marine bivalves to two different chemicals and their mixture determined biological responses with discrete quantitative levels through which it was possible to rank the health status of the organisms according to their stress syndrome. According to the biomarker expert system outcome, animals exposed to the mixture exhibited a better health status compared to those exposed to single chemicals at the same nominal toxic dose (Figure 3). A particular evidence it was the suppressive effect of Chlorpyrifos over the hyperlipidosis syndrome diagnosed in digestive gland lysosomes of Ni exposed mussels (Figure 2). However, the occurrence of deviations from the common mixture toxicity reference models based on non-interaction was confirmed for all biomarkers, testing the significance through maximum likelihood analysis (Table 2).
To get clues on interactions between heavy metal and pesticide toxicity, we first looked at the toxicokinetics but it did not provide an explanation for the attenuation in toxicity observed in mixture-exposed animals. In fact, looking at the parameters obtained by fitting the data into a kinetic model, the elimination constant (k) of CHP in the absence/presence of the heavy metal is almost identical, indicating no changes in the metabolism of the biocide. Moreover, the elimination constant of Ni indicated an even lower degree of detoxification when CHP was jointly administered to mussels (Table 3). We conclude, therefore, that the decrease in toxicity observed at a biological level could not arise from interaction at a toxicokenetic level as reported by  for the soil-dwelling collembolum Folsomia candida.
As it is well known that the mussel digestive gland is the most active metabolic organ  suitable for ecotoxicogenomic profiling [19, 15], we, therefore, focused our attention to this tissue. Starting from nominal equitoxic effects on lysosomal membrane stability (Table 1), we carried out transcriptomic profiling using the 1.7 K cDNA chip (see Additional file 2; Figure 4) and further gene ontology-based functional genomics analysis (Figure 5-6, see Additional file 3). These analyses provided a picture of the biological processes and toxicodynamics implicated in the response to the two pollutants and their mixture. In the case of single exposures, our results pointed out two distinct gene expression profiles (Figure 4). Divergent gene expression profiles have already been reported in another model species, i.e. differentiating PC12 mouse brain cells challenged with the same chemicals [29, 45, 46]. In the case of the mixture, a complex pattern was observed which accounted for some genes inherited from the single chemicals -in particular Chlorpyrifos- showing the same relative expression trend; a set of unique sequences and a fraction of genes showing an opposed relative expression trend with respect to that observed in single chemical exposures (Table 5).
Ni exposure increased the expression level of the metallothionein mt10 gene, while negligible effects were observed for the cognate sequence mt20 (Table 4). In addition, transcriptomic profiling, by means of microarray analysis, allowed the detection of several genes which appear to be involved in epigenetic processes, as already reported for different heavy metals and also other environmental perturbations [19, 47]. In fact, among the significant molecular features observed exclusively in Ni-treated organisms there are those putatively involved in spliceosome assembly and the establishment of chromatin architecture, as outlined by the Blast2GO analysis (Figure 5; see Additional file 3). Ni elicited the specific down-regulation of four different, small, nuclear, ribonucleoprotein which are homologous with components of the Sm core complex  and a poly-A-binding cytoplasmic 4 like-protein, involved in the half-life regulation of labile mRNAs . With the same trend of down-regulation, we found also histone h3, histone h2a, histone aminotransferase 1 and a heterochromatin protein 1 family member which, in different organisms, is reported to bind histone H3 methylated at Lys 9, thus leading to gene silencing through heterochromatin promotion .
Moreover, Ni modulated several genes involved in the regulation of biological processes such as proliferation, growth and apoptosis (Figure 6; see Additional file 4). Examples are given by a b-cell translocation gene 1 homologue (upregulated) with putative anti-proliferative functions ; prohibitin (downregulated) which is an ubiquitous protein with a controversial role in cell proliferation processes  and biogenesis of mitochondria ; a (downregulated) defender-against-cell-death (dad-1)-like gene, identified as a negative regulator of programmed cell death in mammalian cells and with a role in tissue regeneration in the marine scallop Argopecten irradians . Finally, Ni exposure elicited the modulation of genes involved in lipid metabolism with two genes whose expression was up-regulated: an apolipophorin putatively involved in the transport of lipids to hemolymph in insects  -whose differential expression was also confirmed by Q-PCR (Table 4)- and a fatty acid-binding protein [GeneBank:AJ624395] with a role in the proliferation process . Still concerned with lipid metabolism are two non-allelic variants of the GM2 ganglioside activator protein (GM2-AP) [GeneBank:AJ624405 and GeneBank:AJ624495] whose expression levels were dramatically up- and down-regulated, respectively (Table 4). GM2-AP are ubiquitous lysosomal proteins  which, in mammals, act as substrate-specific co-factors of β-hexosaminidase-A for the hydrolysis of GM2 ganglioside, a class of glycolipids positive to oil-red staining. Defects on both lysosomal β-hexosaminidase A and GM2-AP result in a fatal syndrome characterized by a hyper-accumulation of glycolipids in neuronal cells . These two genes could represent the candidates for the strong increase of lipid accumulation in the digestive gland cells of mussels (Figure 1, panel C; Figure 2). The presence of two differentially expressed GM2-AP variants might represent a transcriptional mechanism to regulate lysosomal hexosaminidase activity, according to the physiological requirements of the organism. To this aim, it is important to point out that different GM2-AP expression patterns corresponded to different lipid disorder levels (Figure 1, panel C; Figure 2). The huge relative expression levels revealed by Q-PCR would suggest a typical switch-on/off transcriptional regulation.
Mussels exposed to Chlorpyrifos exhibited a marked decrease in acetylcholinesterase activity in the gills, independently of the concentration tested (Figure 1, panel D) and a similar effect was also previously demonstrated in the digestive gland by . Acetylcholinesterase inhibition is a well-known biomarker of exposure to organophosphate pesticides in several organisms. However, in marine invertebrates, heavy metals may affect such enzymatic activity too , and, in fact, we herein report the suppressive effect of Ni (Figure 1, panel D). In what concerns the transcriptomic assessment, the most relevant biological process involved in response to Clorpyrifos was represented by carbohydrate catabolic process, in particular that related to chitinase activity (Figure 5). In fact, up to 5 different highly homologous chitinase genes were found hugely up-regulated from microarray analysis. TaqMan assays set up for four of these sequences confirmed the relative expression levels (Table 4). It is difficult to speculate on the biological meaning of such findings; however, in mussels and other marine invertebrates, chitinases play a role in digestion  and moreover, in hemocytes, they participate in the innate immune response [60, 61]. Finally, GO enrichment analysis identified a series of statistically over-represented developmental and differentiation processes, even if these were related to the down-regulation of a single myosinase-like-gene.
Very little is known about the relationships between mixture toxicity and gene expression changes. Some studies reported that the transcriptional patterns found in mixture-exposed samples were largely inherited from the single chemicals, suggesting additive or weak interactive effects at a molecular level [62, 63]. Conversely, other works reported much more important interactions so that exposure to a mixture determines -along with overlapping genes- a relevant number of sequences whose expression is exclusively modulated in this condition. These findings were observed in different biological models -rat, fish and crustacean- and for different types of chemical mixtures, i.e. binary, ternary, similarly- and dissimilarly-acting chemicals [64–66]. In our study, multivariate analysis of gene expression patterns showed that the molecular signature observed in mixture exposed samples is biased towards the effects of the pesticide (Figure 4). A similar evidence emerged also from the comparison of functional gene annotations in which GO terms such as carbohydrate catabolic process, multicellular organismal development, anatomical structure morphogenesis, signal transduction were shared with the pesticide. Nevertheless, the exposure to the mixture elicited the modulation of several unique genes giving rise to a molecular fingerprint which appears to be characterized by at least two original features: catalytic and DNA-binding (including transcription factor) activity (Figure 5, see Additional file 3). The direct comparison of GO term distributions demonstrated that these attributes were driven neither by Ni nor Chlorpyrifos (Figure 5-6). This kind of composite signature bearing common and unique features is consistent with that found by our research group in mussels exposed to a binary mixture of neonicotinoids insecticides . Also Vandenbrouck and coworkers , studying transcriptomic effects of binary mixtures (nickel with other heavy metals) in the freshwater cladoceron Daphnia magna, found affected genes and pathways which were exclusive of the mixture exposure. More recently, another research group investigated in the worm Caenorhabditis elegans the joint effects of Chlorpyrifos and another organic phosphorus compound, diazinon, by means of gene expression techniques. They still concluded their work stating that the effect of a mix of low doses of the two biocides is not a summed effect of the single components, but at the same time, the similarities in the evoked pathways indicate the regulation of similar responses . In general, these findings are in accordance with ours and therefore the presence of common and original responses seem to represent a common rule. We, however, identified one additional feature. In the present study, the three toxic treatments modulated the expression of some genes involved in lipid metabolism. This process was even significantly over-represented in Ni and mixture exposed tissues (Figure 5-6, see Additional file 3 and 4). In compliance with Vinuela and coworkers , not all transcripts of this group were identical among the three treatments (see Additional files 3, 4), except the two GM2-AP. These two transcripts, however, showed diametrically opposed expression levels and very huge fold change which cannot be explained simply by means of an additivity model (Table 4). This finding seriously poses for the implication of similar pathways but with a different biological meaning and in fact the outcome of lysosomal biomarkers was markedly different (Figure 1). Lipid metabolism was not a single case because these two molecular signatures were also characterized by a certain number of (dissimilar) genes with a possible role in apoptosis and cell proliferation (for the mixture: gadd45, caspase7, translationally controlled tumor protein (tctp) [GeneBank:AJ624761]; for Ni: b-cell translocation gene-1 homologue, prohibitin, dad-1), which considering their putative role in other model species, might be in contrast at a functional level, i.e. pro-apoptotic and pro-proliferative, respectively. However, this hypothesis remains open and requires further investigations, as we did not evaluated specific functional assays for such processes.
We presented an analysis and comparison of the biological responses elicited by the exposure of marine mussels to two toxicants and a combination of both using a systems toxicology approach. Our results demonstrated the occurrence of interactive effects giving rise to an unpredicted response and finally to a decrease of toxicity. We identified and confirmed the core features of mixtures' gene expression profiles and highlighted novel developments. Furthermore, our findings underlined the fallacy of the "non-interaction" criterion usually applied to mixture toxicity prediction. This was particularly evident at a molecular level such as gene expression, thus putting a serious concern in the adoption of conventional mixture toxicity reference models in ecotoxicological surveys and risk assessment procedures.
This work demonstrated also that an integrated approach made of of transcriptomics, functional genomics, cellular and histological biomarkers can provide clues on complex biological responses, making links between different levels of biological organization, as we obtained in the case of lipid metabolism genes and lysosomal hyperlipidosis in the digestive gland.
Future implementation and development of massive sequencing applications and high density microarrays will fill the gap of genomic/transcriptomic information which is actually the major limitation in the use of Mytilus galloprovincialis as a model organism, thus providing more robust assessments into an ecologically relevant species ubiquitously present along the coastal marine environments.
Nickel was used in the form of chloride salt, obtained from Sigma Aldrich. Chlorpyrifos-ethyl was obtained from India Industrie Chimiche SPA (Padova, Italy). All other reagents were of analytical grade or "Chromasolv" grade for chemical analyses if not otherwise stated.
For the range-finding tests, specimens of Mytilus galloprovincialis (5-6 cm length) were taken from a mussel farm in Cesenatico, North-East-Italy during March 2005, and transferred to aquaria with re-circulating aerated seawater, at 16°C collected offshore, at a density of 1 animal/L. After an acclimatizing period of 6 days, mussels were divided in groups of 5 vessels per 30 individuals each (150 mussels per condition) and further used as control reference or subjected to semi-static exposure to chemicals. Each vessel represented a parallel replicated experiment. Chemicals were administered every day, together with a commercial algal preparation (Liquifry, Interpret Ltd., Dorking, Surrey, UK) and seawater renewed every two days. Nickel was administered as the chloride salt (NiCl2) from a concentrated stock solution (5000 X), while Chlorpyrifos was diluted in dimethylsulphoxide (DMSO) and added at the desired concentration from concentrated stock solutions (5000 X). DMSO was also added to a control seawater and Ni-exposed mussels at the concentration of 0.02%. In all experiments female individuals, screened by microscopic inspection of gonads, were used for subsequent analyses.
Mussels (40 individuals per condition) were exposed for 4 days to nickel and Chlorpyrifos respectively in the range of 0.01-15 mg/L and 0.1-10 mg/L. Digestive gland cryostat sections were further scored for lysosomal membrane stability (LSM). Control groups of vehicle-treated mussels were kept in the same conditions as the ones exposed to the chemical.
For mixture toxicity analysis, a fixed design was used, encompassing the nominal effective concentrations (ECs) for LSM obtained from the dose-range finding experiments. The following endpoints were selected: for single chemicals, EC10, EC25, EC50 (Table 1); mixtures, NiEC10+CHPEC10 and NiEC25+CHPEC25. Mussels were exposed for 4 days as previously indicated. After treatments, digestive glands from female sex specimens were rapidly removed and washed in artificial seawater buffered with 20 mM Hepes pH 7.4 and stored according to further analysis. For transcriptomics, the tissue was kept at -20°C in a RNA preserving solution (RNA Later, Sigma-Aldrich); for histochemistry the digestive gland of 5 animals sampled from each of the five vessels was mounted on aluminum chucks and frozen in super-cooled n-hexan as previously described . Two different chucks were prepared per group of mussels. The rest of the tissue, once excised, was snap-frozen in liquid nitrogen and stored at -80°C until analysis.
For toxicokinetics assessment, mussels (30 individuals per condition) were exposed for 4 days to nickel or Chlorpyrifos at the nominal EC50 for effects on LMS and to NiEC25+CHPEC25 in the case of the mixture. After the intoxication period, animals were transferred to clean seawater for 6 days. Sampling of test animals was performed at 3, 6, 12, 24 hours, 3 and 6 days. One vessel per condition, each containing 30 mussels was set up.
Cryostat sections (10 μm) were obtained through a Leica cryostat apparatus at -27°C. LSM was evaluated using the method described by . Staining intensity of lysosomes was obtained by means of an inverted Axiovert microscope (Zeiss) at 400×magnification, connected to a digital camera (Axiocam, Zeiss). Digital image analysis was carried out using the Scion Image software package (Scion Corp. Inc.) from 8-bit gray scale images.
The lipid content was assessed by staining tissue sections with the oil-soluble dye, Oil Red-O (ORO) as indicated by  and quantified by digital image analysis as described above.
Acetylcholinesterase (AChE) activity in gill extracts was determined using the Acetylcholinesterase Reagent Kit (Ikzus) using 2 mM acetylthiocholine as substrate, essentially as described in . For each condition, pools of gill portions obtained from 6 animals sampled from each of the five vessels were used for the analysis.
Dual color competitive hybridizations were carried out following a common reference design in which each experimental condition was hybridized against the same reference condition, i.e. digestive gland tissue from vehicle-treated animals. Five different biological replicates were obtained from each of the five vessels and further used to analyze each condition, with the exception of CHP for which seven samples were analyzed in an attempt to obtain more differentially expressed genes (negligible differences were obtained using five or seven arrays, data not shown). One replicate per array was used. Total RNA was extracted from pools of 6 digestive gland pieces by the acid phenol-chloroform procedure according to , using the TRI-Reagent (Sigma-Aldrich). RNA was further purified by precipitation in the presence of 1.5 M LiCl2. The quality of each RNA preparation was verified both by UV spectroscopy and TBE agarose gel electrophoresis, in the presence of formamide as described by . Competitive dual-color microarray hybridization analysis was performed using the Mytarray V1.0 and V1.1 platform [15, 67]. This array encompasses 3' cDNA probes representing 1748 independent mussel sequences obtained from unbiased M. galloprovincialis tissues-specific normalized cDNA libraries. cDNA fluorescence-labeled probes were obtained by the direct labeling procedure in the presence of modified cy-3 and cy5 dCTP (Perkin Elmer) . The procedure was carried out from 15 μg total RNA essentially as described by  with the exception that first strand synthesis was carried out with 0.5 μg of an anchored oligo dT(19)VN instead of random examers. Microarray slides pre-hybridised with the formamide-based buffer Northern Max (Ambion) for 1-2 h at 42°C were further hybridized overnight at 42°C with cDNA probes resuspended in 20 μl of the same buffer. After hybridisation, slides were washed to remove excess probes and unspecific binding. Three washing steps were carried out: first washing in 1× SSC, 0.2% SDS for 5 min; second washing in 0.1× SSC, 0.2% SDS for 3 min; third washing in 0.1×SSC for 5 min, the last repeated once for 3 min. All washing steps were performed with gentle shaking at room temperature. Laser scanning of microarrays was performed using a ChipReader at 5 micron resolution (Bio-Rad Laboratories, CA, USA). 16 bit TIFF images were analyzed by means of Genepix 6.0 (Axon) to get raw fluorescence data from each spot. Pre-processing and differentially expressed genes were obtained by means of an R-based package LIMMA  through the implementation of empirical Bayes statistics. B > 0 was used, where B-statistics represents the log-odds that a particular gene is differentially expressed. Microarray experiments were deposited in the Gene Expression Omnibus (GEO) database with the Series record [GSE21229].
Q-PCR analysis was carried out from 25 ng RNA reverse-transcribed cDNA obtained from the same pools used for microarray hybridization. Reverse transcription was performed from 1 μg total RNA according to . Relative expression levels of the following genes, actin [GeneBank:AJ625116], GM2-AP [GeneBank:AJ624495, GeneBank:AJ624405], hexosaminidase [GeneBank:AJ623463], apolipophorin precursor [GeneBank:AJ625863], chitinase [GeneBank:AJ624093, GeneBank:AJ625569, GeneBank:AJ624637, GeneBank:AJ625051] were expressed as group mean and geometrically normalized on 18S rRNA (L33452) and an invariant alkaline phosphatase gene [Geneank:AJ626187]. Genes of interest were amplified into a CFX384 quantitative thermal cycler (Bio-Rad) in triplex TaqMan assay in the presence either of the "No Rox Multiplex Quantitect" (Qiagen) or "iQTM Multiplex Powermix" (Bio-Rad) master mix. In all cases except one, it was followed the protocol for triplex mode with an annealing temperature of 60°C. In the case of the chitinase gene [GeneBank:AJ625051] the annealing temperature was selected through a preliminary gradient run and further set for the quantitative analysis at 47.7°C. Reference genes were run in a duplex assay. To this aim, 0.25 ng RNA reverse-transcribed to cDNA was amplified in the presence of 0.1 μM each dual labelled probe (Texas Red/BH2 for 18S rRNA; HEX/BH1 for alkaline phosphatase), 0.1 μM and 0.4 μM each forward and reverse primer pairs, respectively for 18S rRNA and alkaline phosphatase. The thermal protocol using the iQTM Multiplex Powermix" (Bio-Rad) was as follows: 30 sec at 95°C, followed by 40 cycles (10 s at 95°C, 20 s at 60°C) into a CFX 384 Bio-Rad PCR apparatus. All Q-PCR reactions were performed on four biological replicates and three technical replicates. Primers and probes -designed by means of Beacon Designer V. 3.0 software (Premier Biosoft International, Inc.) - are reported in Additional file 5. Metallothionein (mt10 and mt20) mRNA levels were evaluated as previously described . The p53-like mRNA level [Genebank:AJ625243] was evaluated as previously reported . All statistical computation and analysis of Q-PCR data were carried out using the REST and REST-mcs software .
The functional characterization of mussel genes present in the array was based on Gene Ontology annotation and was carried out by means of the universal platform Blast2GO (B2GO)  using default parameters. GO term enrichment analysis was carried out through the implementation of a hypergeometric statistic (p < 0.05) .
Nickel was determined in mussel tissues (about 0.5 g of a 1:1 homogenate in double distilled water) by inductively coupled plasma-mass spectrometry (ICP-MS). Samples were added of 5 ml of concentrated 65% nitric acid and introduced into a microwave oven for the mineralization. The sample was then filtered on a nitrocellulose membrane (0.45 μm) and the Ni quantified using a VG Plasma Quad 3 (VG Elemental) Inductively Couple Plasma (ICP) Mass Spectrometer. Procedure validation was performed using the Std CRM 145 R reference material containing added known amounts of metal . Chlorpyrifos-ethyl residual analysis was performed in the tissue of mussels (10 g) by means of homogenization in the presence of anhydrous NaSO4 (10 g), dichloromethane Soxhlet extraction (5 h, with reflux rate 3-5 min), concentration in rotavapor (T = 55°), and further Gas Chromatography Mass Spectrometry, equipped with an Electron Capture Detector (ECD). GC conditions were as follow: Splitless injector (splitless time: 1 min), T = 250°C; Carrier: The constant flux at 1 ml/min; Temperature gradient: 100°C 1 min; 5°C/min up to 250°C for 1 min; 30°C/min up to 300°C for 10 min. Detector settings were the following: transfer line, T = 300°C; Source, T = 200°C. Chlorpyrifo-ethyl ion was detected at 314 m/z. Internal and external calibration procedures were performed. The analytical measurements were carried out in triplicate from pools of 4-5 animals.
Biomarkers data where analyzed with the Mann Whitney U-test (n = 10). Log-logistic regression curves describing dose-dependent effects for LMS were obtained using the software package Sigma Plot 9 (Systat Inc.). The health status of mussels has been determined applying an expert system for classification able to rank the stress syndrome evolution by integrating the results from a battery of biomarkers . Results on pollutant-exposed organisms are compared with those obtained from control animals, applying the non-parametric statistical tests Mann-Whitney U-test. Significant changes (p < 0.05) are converted into alteration levels (ALs) by comparison with specific thresholds. Finally data are integrated into a health status index (HSI), ranging from A (healthy) to E (pathologically stressed). HSI levels are calculated integrating ALs with a classification algorithm based on rules in the "if...then..." form: synthetically, HSI depends on the number of altered biomarkers and on the level of biological organization affected by pollutants (i.e. cell, tissue, organism) .
Statistically significant deviations of biomarkers responses from mixture toxicity reference models were analyzed using the MIXTOX spreadsheet model developed by . Briefly, the MIXTOX software allows deviations from either the Concentration Addition  or Independent Action  model to be significance tested and analyzed in detail using a response surface analysis framework designed to identify biologically relevant response patterns like overall synergy, antagonism as well as more complex issues of ration and effect level dependent deviations.
The toxicokinetic models were computed using Sigma SYSTAT 10 from the following equation .
Ct = C0 · e -kt
Where C t is the pollutant concentration in an animal at time t, C 0 is the pollutant concentration at the beginning of the detoxification phase and k is the elimination rate constant.
This research was supported by a grant from the UE 6th FW Program, NoMiracle, IP Contract number, 003956. The MytArray slides were developed and printed by CRIBI BIOTECHNOLOGY CENTER - University of Padova, Via Ugo Bassi, No.58/B, 35121 Padova, Italy. Authors acknowledge Dr. Flavia Capri for technical assistance in cytochemical analysis. All auhors read and approved the final manuscript.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.