BackgroundMixtures 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.ResultsCellular 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.ConclusionThe 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.