Characterization and Comparison of Ocular Surface Microbiome in Newborns.

The ocular microbiome is of fundamental importance for immune eye homeostasis, and its alteration would lead to an impairment of ocular functionality. Little evidence is reported on the composition of the ocular microbiota of term infants and on the impact of antibiotic prophylaxis. METHODS: A total of 20 conjunctival swabs were collected from newborns at birth and after antibiotic treatment. Samples were subjected to 16S rRNA sequencing via system MiSeq Illumina. The data were processed with the MicrobAT software and statistical analysis were performed using two-way ANOVA. RESULTS: Antibiotic prophylaxis with gentamicin altered the composition of the microbiota. In detail, a 1.5- and 2.01-fold reduction was recorded for Cutibacterium acnes (C. acnes) and Massilia timonae (M. timonae), respectively, whereas an increase in Staphylococcus spp. of 6.5 times occurred after antibiotic exposure. CONCLUSIONS: Antibiotic prophylaxis altered the ocular microbiota whose understanding could avoid adverse effects on eye health.

Native AMF Communities in an Italian Vineyard at Two Different Phenological Stages of Vitis vinifera.

Arbuscular mycorrhizal fungi (AMF) are beneficial soil microorganisms that can establish symbiotic associations with Vitis vinifera roots, resulting in positive effects on grapevine performance, both in terms of water use efficiency, nutrient uptake, and replant success. Grapevine is an important perennial crop cultivated worldwide, especially in Mediterranean countries. In Italy, Piedmont is one of the regions with the longest winemaking tradition. In the present study, we characterized the AMF communities of the soil associated or not with the roots of V. vinifera cv. Pinot Noir cultivated in a vineyard subjected to conventional management using 454 Roche sequencing technology. Samplings were performed at two plant phenological stages (flowering and early fruit development). The AMF community was dominated by members of the family Glomeraceae, with a prevalence of the genus Glomus and the species Rhizophagus intraradices and Rhizophagus irregularis. On the contrary, the genus Archaeospora was the only one belonging to the family Archaeosporaceae. Since different AMF communities occur in the two considered soils, independently from the plant phenological stage, a probable role of V. vinifera in determining the AMF populations associated to its roots has been highlighted.

Climatic Zone and Soil Properties Determine the Biodiversity of the Soil Bacterial Communities Associated to Native Plants from Desert Areas of North-Central Algeria.

Algeria is the largest country in Africa characterized by semi-arid and arid sites, located in the North, and hypersaline zones in the center and South of the country. Several autochthonous plants are well known as medicinal plants, having in common tolerance to aridity, drought and salinity. In their natural environment, they live with a great amount of microbial species that altogether are indicated as plant microbiota, while the plants are now viewed as a "holobiont". In this work, the microbiota of the soil associated to the roots of fourteen economically relevant autochthonous plants from Algeria have been characterized by an innovative metagenomic approach with a dual purpose: (i) to deepen the knowledge of the arid and semi-arid environment and (ii) to characterize the composition of bacterial communities associated with indigenous plants with a strong economic/commercial interest, in order to make possible the improvement of their cultivation. The results presented in this work highlighted specific signatures which are mainly determined by climatic zone and soil properties more than by the plant species.

Microbiota thrombus colonization may influence athero-thrombosis in hyperglycemic patients with ST segment elevation myocardialinfarction (STEMI). Marianella study.

OBJECTIVES: We examined the association of the coronary thrombus microbiota and relative metabolites with major adverse cardiovascular events (MACE) in hyperglycemic patients with ST segment elevation myocardial infarction (STEMI). BACKGROUND: Hyperglycemia during STEMI may affect both development and progression of coronary thrombus via gut and thrombus microbiota modifications. METHODS: We undertook an observational cohort study of 146 first STEMI patients treated with primary percutaneous coronary intervention (PPCI) and thrombus-aspiration (TA). Patients were clustered, based on admission blood glucose levels, in hyperglycemic (≥140 mg/dl) and normoglycemic (<140 mg/dl). We analyzed gut and thrombus microbiota in all patients. Moreover, we assessed TMAO, CD40L and von Willebrand Factor (vWF) in coronary thrombi. Cox regressions were used for the association between Prevotellaspp. and TMAO terziles and MACE. MACE endpoint at 1 year included death, re-infarction, unstable angina. RESULTS: In fecal and thrombus samples, we observed a significantly different prevalence of both Prevotellaspp. and Alistipesspp. between patients with hyperglycemia (n = 56) and those with normal glucose levels (n = 90). The abundance of Prevotella increased in hyperglycemic vs normoglycemic patients whereas the contrary was observed for Alistipes. Interestingly, in coronary thrombus, the content of Prevotella was associated with admission blood glucose levels (p < 0.01), thrombus dimensions (p < 0.01), TMAO, CDL40 (p < 0.01) and vWF (p < 0.01) coronary thrombus contents. Multivariate Cox-analysis disclosed a reduced survival in patients with high levels of Prevotella and TMAO in coronary thrombus as compared to patients with low levels of Prevotella and TMAO, after 1-year follow up. CONCLUSIONS: Hyperglycemia during STEMI may increase coronary thrombus burden via gut and thrombus microbiota dysbiosis characterized by an increase of Prevotella and TMAO content in thrombi. CLINICAL TRIAL REGISTRATION: NCT03439592. September 30, 2016. Ethic Committee Vanvitelli University: 268/2016.

Controlled ovarian stimulation and progesterone supplementation affect vaginal and endometrial microbiota in IVF cycles: a pilot study.

PURPOSE: Does controlled ovarian stimulation (COS) and progesterone (P) luteal supplementation modify the vaginal and endometrial microbiota of women undergoing in vitro fertilization? METHODS: Fifteen women underwent microbiota analysis at two time points: during a mock transfer performed in the luteal phase of the cycle preceding COS, and at the time of fresh embryo transfer (ET). A vaginal swab and the distal extremity of the ET catheter tip were analyzed using next-generation 16SrRNA gene sequencing. Heterogeneity of the bacterial microbiota was assessed according to both the Bray-Curtis similarity index and the Shannon diversity index. RESULTS: Lactobacillus was the most prevalent genus in the vaginal samples, although its relative proportion was reduced by COS plus P supplementation (71.5 ± 40.6% vs. 61.1 ± 44.2%). In the vagina, an increase in pathogenic species was observed, involving Prevotella (3.5 ± 8.9% vs. 12.0 ± 19.4%), and Escherichia coli-Shigella spp. (1.4 ± 5.6% vs. 2.0 ± 7.8%). In the endometrium, the proportion of Lactobacilli slightly decreased (27.4 ± 34.5% vs. 25.0 ± 29.9%); differently, both Prevotella and Atopobium increased (3.4 ± 9.5% vs. 4.7 ± 7.4% and 0.7 ± 1.5% vs. 5.8 ± 12.0%). In both sites, biodiversity was greater after COS (p < 0.05), particularly in the endometrial microbiota, as confirmed by Bray-Curtis analysis of the phylogenetic distance among bacteria genera. Bray-Curtis analysis confirmed significant differences also for the paired endometrium-vagina samples at each time point. CONCLUSIONS: Our findings suggest that COS and P supplementation significantly change the composition of vaginal and endometrial microbiota. The greater instability could affect both endometrial receptivity and placentation. If our findings are confirmed, they may provide a further reason to encourage the freeze-all strategy.

Discovering the bacteriome of Vitis vinifera cv. Pinot Noir in a conventionally managed vineyard.

The structure of the bacteriome associated with grapevine roots can affect plant development, health and grape quality. We previously investigated the bacterial biodiversity of the Vitis vinifera cv. Pinot Noir rhizosphere in a vineyard subjected to integrated pest management. The aim of this work is to characterize the bacteriome of V. vinifera cv. Pinot Noir in a conventionally managed vineyard using a metabarcoding approach. Comparisons between the microbial community structure in bulk soil and rhizosphere (variable space) were performed and shifts of bacteriome according to two sampling times (variable time) were characterized. Bacterial biodiversity was higher at the second than at the first sampling and did not differ according to the variable space. Actinobacteria was the dominant class, with Gaiella as the most represented genus in all the samples. Among Proteobacteria, the most represented classes were Alpha, Beta and Gamma-Proteobacteria, with higher abundance at the second than at the first sampling time. Bradyrhizobium was the most frequent genus among Alpha-Proteobacteria, while Burkholderia was the predominant Beta-Proteobacteria. Among Firmicutes, the frequency of Staphylococcus was higher than 60% in bulk soil and rhizosphere. Finally, the sampling time can be considered as one of the drivers responsible for the bacteriome variations assessed.

Cellulomonas fimi secretomes: In vivo and in silico approaches for the lignocellulose bioconversion.

Lignocellulose degradation is a challenging step for value added products and biofuels production. Cellulomonas fimi secretes complex mixtures of carbohydrate active enzymes (CAZymes) which synergistically degrade cellulose and hemicelluloses. Their characterization may provide new insights for enzymatic cocktails implementation. Bioinformatic analysis highlighted 1127 secreted proteins, constituting the in silico secretome, graphically represented in a 2DE map. According to Blast2GO functional annotation, many of these are involved in carbohydrates metabolism. In vivo secretomes were obtained, growing C. fimi on glucose, CMC or wheat straw for 24 h. Zymography revealed degradative activity on carbohydrates and proteomic analysis identified some CAZymes, only in secretomes obtained with CMC and wheat straw. An interaction between cellobiohydrolases is proposed as a strategy adopted by soluble multimodular cellulases. Such approach can be crucial for a better characterization and industrial exploitation of the synergism among C. fimi enzymes.

The Rhizosphere Bacterial Microbiota of Vitis vinifera cv. Pinot Noir in an Integrated Pest Management Vineyard.

Microorganisms associated with Vitis vinifera (grapevine) can affect its growth, health and grape quality. The aim of this study was to unravel the biodiversity of the bacterial rhizosphere microbiota of grapevine in an integrated pest management vineyard located in Piedmont, Italy. Comparison between the microbial community structure in the bulk and rhizosphere soil (variable: space) were performed. Moreover, the possible shifts of the bulk and rhizosphere soil microbiota according to two phenological stages such as flowering and early fruit development (variable: time) were characterized. The grapevine microbiota was identified using metagenomics and next-generation sequencing. Biodiversity was higher in the rhizosphere than in the bulk soil, independent of the phenological stage. Actinobacteria were the dominant class with frequencies ≥ 50% in all the soil samples, followed by Proteobacteria, Gemmatimonadetes, and Bacteroidetes. While Actinobacteria and Proteobacteria are well-known as being dominant in soil, this is the first time the presence of Gemmatimonadetes has been observed in vineyard soils. Gaiella was the dominant genus of Actinobacteria in all the samples. Finally, the microbiota associated with grapevine differed from the bulk soil microbiota and these variations were independent of the phenological stage of the plant.

Toxic effects of mercury on the cell nucleus of Dictyostelium discoideum.

Governmental agencies (www.epa.gov/mercury) and the scientific community have reported on the high toxicity due to mercury. Indeed, exposure to mercury can cause severe injury to the central nervous system and kidney in humans. Beyond its recognized toxicity, little is known regarding the molecular mechanisms involved in the actions of this heavy metal. Mercury has been also observed to form insoluble fibrous protein aggregates in the cell nucleus. We used D. discoideum to evaluate micronuclei formation and, since mercury is able to induce oxidative stress that could bring to protein aggregation, we assessed nuclear protein carbonylation by Western Blot. We observed a significant increase in micronuclei formation and 14 carbonylated proteins were identified. Moreover, we used isotope-coded protein label (ICPL) and mass spectrometry analysis of proteins obtained by lysis of purified nuclei, before of tryptic digestion to quantify nuclear proteins affected by mercury. In particular, we examined the effects of mercury that associate a classical genotoxic assay to proteomic effects into the nucleus. The data present direct evidences for mercury genotoxicity, nuclear protein carbonylation, quantitative change in core histones, and the involvement of pseudouridine synthase in mercury toxicity. © 2016 Wiley Periodicals, Inc. Environ Toxicol 32: 417-425, 2017.

Arbuscular mycorrhizal symbiosis affects the grain proteome of Zea mays: a field study.

Maize is one of the most important crops worldwide and is strongly dependent on arbuscular mycorrhiza (AM) fungi, organisms that form a mutualistic association with land plants. In maize, AM symbiosis enhances spike dry weight, spike length, spike circumference, and the dry weight and dimensions of the grain. Notwithstanding its ubiquitous nature, the detailed relationship between AM fungal colonization and plant development is not completely understood. To facilitate a better understanding of the effects of AM fungi on plants, the work reported here assessed the effects of a consortium of AM fungi on the kernel proteome of maize, cultivated in open-field conditions. To our knowledge, this is the first report of the modulation of a plant seed proteome following AM fungal inoculation in the field. Here, it was found that AM fungi modify the maize seed proteome by up-regulating enzymes involved in energetic metabolism, embryo development, nucleotide metabolism, seed storage and stress responses.

NGS-Trex: an automatic analysis workflow for RNA-Seq data.

RNA-Seq technology allows the rapid analysis of whole transcriptomes taking advantage of next-generation sequencing platforms. Moreover with the constant decrease of the cost of NGS analysis RNA-Seq is becoming very popular and widespread. Unfortunately data analysis is quite demanding in terms of bioinformatic skills and infrastructures required, thus limiting the potential users of this method. Here we describe the complete analysis of sample data from raw sequences to data mining of results by using NGS-Trex platform, a low user interaction, fully automatic analysis workflow. Used through a web interface, NGS-Trex processes data and profiles the transcriptome of the samples identifying expressed genes, transcripts, and new and known splice variants. It also detects differentially expressed genes and transcripts across different experiments.

NGS-Trex: Next Generation Sequencing Transcriptome profile explorer.

BACKGROUND: Next-Generation Sequencing (NGS) technology has exceptionally increased the ability to sequence DNA in a massively parallel and cost-effective manner. Nevertheless, NGS data analysis requires bioinformatics skills and computational resources well beyond the possibilities of many "wet biology" laboratories. Moreover, most of projects only require few sequencing cycles and standard tools or workflows to carry out suitable analyses for the identification and annotation of genes, transcripts and splice variants found in the biological samples under investigation. These projects can take benefits from the availability of easy to use systems to automatically analyse sequences and to mine data without the preventive need of strong bioinformatics background and hardware infrastructure. RESULTS: To address this issue we developed an automatic system targeted to the analysis of NGS data obtained from large-scale transcriptome studies. This system, we named NGS-Trex (NGS Transcriptome profile explorer) is available through a simple web interface http://www.ngs-trex.org and allows the user to upload raw sequences and easily obtain an accurate characterization of the transcriptome profile after the setting of few parameters required to tune the analysis procedure. The system is also able to assess differential expression at both gene and transcript level (i.e. splicing isoforms) by comparing the expression profile of different samples.By using simple query forms the user can obtain list of genes, transcripts, splice sites ranked and filtered according to several criteria. Data can be viewed as tables, text files or through a simple genome browser which helps the visual inspection of the data. CONCLUSIONS: NGS-Trex is a simple tool for RNA-Seq data analysis mainly targeted to "wet biology" researchers with limited bioinformatics skills. It offers simple data mining tools to explore transcriptome profiles of samples investigated taking advantage of NGS technologies.

Dictyostelium discoideum developmental cycle (DDDC) assay: a tool for Hg toxicity assessment and soil health screening.

The social amoeba Dictyostelium discoideum has been proposed for assessing stress responses to pollutants in soil and it has already been successfully employed in the aquatic environment. Presently, we developed the DDDC assay (D. discoideum developmental cycle assay) for both soil toxicity assessment and soil health screening. The DDDC assay is primarily aimed at determining the capacity of D. discoideum to undergo its developmental programme forming a fruiting body, measured in terms of fruiting body formation inhibition and fruiting body size factor, which may be considered an indication of its ecological fitness (potential for spore dispersal). A second objective of the solid phase DDDC assay is to identify potential mechanisms of toxic action on the developmental cycle, for which three checkpoints are examined: (a) aggregation arrest, (b) migration arrest, and (c) culmination arrest. Presently, conditions for the DDDC assay such as soil texture, soil water content, soil pH, food availability and incubation time were investigated and optimized. In addition, both solid and liquid phase variants of the DDDC assay were applied to assess the toxicity of Hg, at regulatory concentrations. The developmental cycle and ecological fitness were affected from the exposure to 0.3 mg Hg/kg dry-wt soil onwards. The DDDC assay has been shown to be a high sensitivity test.

Effects of nickel, chlorpyrifos and their mixture on the Dictyostelium discoideum proteome.

Mixtures of chemicals can have additive, synergistic or antagonistic interactions. We investigated the effects of the exposure to nickel, the organophosphate insecticide chlorpyrifos at effect concentrations (EC) of 25% and 50% and their binary mixture (Ec25 + EC25) on Dictyostelium discoideum amoebae based on lysosomal membrane stability (LMS). We treated D. discoideum with these compounds under controlled laboratory conditions and evaluated the changes in protein levels using a two-dimensional gel electrophoresis (2DE) proteomic approach. Nickel treatment at EC25 induced changes in 14 protein spots, 12 of which were down-regulated. Treatment with nickel at EC50 resulted in changes in 15 spots, 10 of which were down-regulated. Treatment with chlorpyrifos at EC25 induced changes in six spots, all of which were down-regulated; treatment with chlorpyrifos at EC50 induced changes in 13 spots, five of which were down-regulated. The mixture corresponding to EC25 of each compound induced changes in 19 spots, 13 of which were down-regulated. The data together reveal that a different protein expression signature exists for each treatment, and that only a few proteins are modulated in multiple different treatments. For a simple binary mixture, the proteomic response does not allow for the identification of each toxicant. The protein spots that showed significant differences were identified by mass spectrometry, which revealed modulations of proteins involved in metal detoxification, stress adaptation, the oxidative stress response and other cellular processes.

Interactions of a pesticide/heavy metal mixture in marine bivalves: a transcriptomic assessment.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

Effects of mercury on Dictyostelium discoideum: proteomics reveals the molecular mechanisms of physiological adaptation and toxicity.

Dictyostelium discoideum amoebae were exposed to Hg 2 microM corresponding to a sublethal concentration and Hg 10 microM with the first effects on mortality and replication rate. A total of 900 spots were visualized by 2-DE electrophoresis. Two-hundred fifty single proteins were identified by mass spectrometry. Low Hg concentration (2 microM) treatment induced up-regulation of 13 spots, mainly involved in oxidative stress response/detoxification, oxidoreductase activity, and metabolic processes. High Hg concentration (10 microM) treatment showed a different PES with 12 proteins downregulated and only two up-regulated, mainly involved in cellular metabolic processes, metal ion binding, and transferase activity. The analyses for the carbonylation show no changes after 2 microM Hg(2+) treatment and 13 differentially carbonylated proteins after 10 microM Hg(2+) involved in a broad range of cellular processes. Our findings provide insight into the mechanisms of physiological adaptation and toxicity to a low and an high mercury concentration, respectively, of Dictyostelium amoebae.

Transcriptomic and proteomic effects of a neonicotinoid insecticide mixture in the marine mussel (Mytilus galloprovincialis, Lam.).

BACKGROUND: Imidacloprid and Thiacloprid are two neonicotinoid insecticides whose use have been raising exponentially. Both act selectively as agonist of the insect nicotinic-Acetylcholine Receptor (nAChR) and therefore, by definition, they hold the same mode of action. Notwithstanding the growing attention to the ecotoxicological effects of neonicotinoids, there is a lack of information on their toxicodynamics and their mixture effects, in particular, in aquatic organisms. OBJECTIVES: The main objectives of this work were: (i) assess sublethal effects of two neonicotinoids-Imidacloprid and Thiacloprid-in the tissues of the marine mussel (Mytilus galloprovincialis Lam.); (ii) identify the molecular dynamics elicited by the two chemicals through gene/protein expression profiling and a functional genomics approach; (iii) assess the effects of a neonicotinoid binary mixture. METHODS: Sublethal effects were measured by means of digestive gland lysosomal membrane stability (LMS) and gill acetylcholinesterase (AChE) activity. Gene expression profiles were evaluated in the digestive gland using a 1.7K cDNA microarray and quantitative-PCR (Q-PCR). Proteome profiling was performed by means of two-dimensional electrophoresis of digestive gland cytosolic proteins. Functional genomics was based on the over-representation of Gene Ontology (GO) terms. RESULTS AND DISCUSSION: Our results showed that (i) biomarkers responded in the micromolar range; (ii) Imidacloprid and Thiacloprid elicited distinct toxicodynamics as depicted by the different transcriptomic and proteomic profiles and the opposite trend of AChE activity; (iii) at biomarkers level, the joint effects of the two chemicals appeared to fulfill the principle of independence, but this was less evident at molecular level where a novel specific molecular signature took place. CONCLUSIONS: These findings imply that different toxicodynamics may occur also as a response of chemicals with the same mode of action. Our results unveil also the incongruousness of the actual concept of pesticide mode of action in the context of ecological risk assessment analysis of chemical mixtures.

Use of highly sensitive sublethal stress responses in the social amoeba Dictyostelium discoideum for an assessment of freshwater quality.

In this work, the sensitivity of a battery of tests on the social amoeba Dictyostelium discoideum has been assessed within a freshwater toxicity study. The results obtained from the evaluation of survival and replication rate of D. discoideum were compared to those derived with a series of widely used tests for freshwater toxicity assessment, i. e. bioassays using Vibrio fischeri, Daphnia magna and Pseudokirchneriella subcapitata. The effects on sublethal endpoints, i.e. lysosomal membrane stability (LMS) and endocytotic rate, were analysed in conjunction with high-level endpoints to verify the potential to make a typical bioassay more sensitive. The field ecotoxicological investigation employing D. discoideum is part of a monitoring study assessing environmental quality of the Bormida River (Italy), subjected until recently to a chronic industrial pollution. The survey was carried out at several stations (upstream and downstream of a chemical factory outlet) in two different periods. In 2002, the results of chemical analyses performed on river water indicated no contamination. The ecotoxicological data obtained in this period showed that no evidence of biological effects was observed using V. fischeri and D. magna bioassays. In spite of the previous classical acute toxicity tests, significant differences in cell viability of D. discoideum were found. By analysing the effects measured on LMS and endocytotic rate, more relevant changes were observed for these sublethal stress biomarkers compared to survival. The chronic toxicity data showed significant changes in cell growth both of P. subcapitata and D. discoideum. Nevertheless, more sensitive and rapid responses were obtained when assessing the effects of exposure on D. discoideum. The chemical and ecotoxicological data obtained in 2006 indicated a full recovery of the quality of the river water (neither contamination nor toxicity found). Altogether, the results reported in this study underline that the use of a battery of biomarkers in conjunction with high-level endpoints may help follow the pollutant-induced stress syndrome in the organisms from early sublethal effects to starting mortality.