- Research article
- Open Access
Cross-species toxicogenomic analyses and phenotypic anchoring in response to groundwater low-level pollution
BMC Genomics volume 15, Article number: 1067 (2014)
Comparison of toxicogenomic data facilitates the identification of deregulated gene patterns and maximizes health risk prediction in human.
Here, we performed phenotypic anchoring on the effects of acute exposure to low-grade polluted groundwater using mouse and zebrafish. Also, we evaluated two windows of chronic exposure in mouse, starting in utero and at the end of lactation. Bioinformatic analysis of livers microarray data showed that the number of deregulated biofunctions and pathways is higher after acute exposure, compared to the chronic one. It also revealed specific profiles of altered gene expression in all treatments, pointing to stress response/mitochondrial pathways as major players of environmental toxicity. Of note, dysfunction of steroid hormones was also predicted by bioinformatic analysis and verified in both models by traditional approaches, serum estrogens measurement and vitellogenin mRNA determination in mice and zebrafish, respectively.
In our report, phenotypic anchoring in two vertebrate model organisms highlights the toxicity of low-grade pollution, with varying susceptibility based on exposure window. The overlay of zebrafish and mice deregulated pathways, more than single genes, is useful in risk identification from chemicals implicated in the observed effects.
In vitro assays for compound toxicity are commonly based on the assumption that toxicants exposure results in changes in gene expression, a biological phenomenon predictive of successive morphological abnormalities [1–3]. Toxicogenomics, defined as changes in genome function that occur with toxicant interaction , is a sensitive, informative and measurable assay to complement traditional toxicological endpoints [5–7]. These advantages prompted the use of toxicogenomics to test the effect of single molecules or simple chemical mixtures [8, 9]. The objectives of transcriptomics in environmental studies (ecotoxicogenomics) are the accomplishment of classical toxicological and new molecular endpoints in the identification of exposure-related alterations, and proper consideration of the complex nature of anthropogenic pollution and bioaccumulation events [10–17]. Besides environments are frequently contaminated with multiple classes of compounds, only a limited number of toxicological studies have recently addressed this problem by using omics approaches to fish species, in environmental field [11, 18–20].
Ecotoxicogenomics is faced with determination of specific patterns of gene expression elicited by environmental samples with known or potential toxicity . Transcriptome analysis has been successfully applied in testing low doses of environmental stressors in biological systems, thus leading to the identification of biomarkers that are easily detectable and related to the observed phenotype, the so called phenotypic anchoring [21, 22]. In this process, the integration of toxicogenomics data from different models is pivotal to validate deregulated patterns, to challenge the low signal to noise ratio and to predict potential risks for human health [23, 24]. Mouse and zebrafish studies indicate that gene expression profile approaches are successful in identifying chemical-specific patterns of altered gene expression [2, 25–27]; for this reason, and for their genetics and biology, these models are widely accepted by the scientific community for environmental toxicology studies [10, 28].
In populations living near waste dumpsites, the correlation between the exposure to chemical mixtures and health disorders has been monitored with different results [29–32]. Typically, low-level exposure to pollutant mixtures is frequently unappreciated and little is known about the consequences of chronic exposure in infants. Among people exposed to contaminants, infants and foetuses are thought to be more susceptible to insults from toxic chemicals because of the period of rapid development [33, 34]. This is an important issue since the adverse effects of a long-term corollary of foetal/neonatal exposure to different pollutants can remain undetected till diseases develop in the adulthood. Several studies have investigated the leachate composition [35–37] and related cytotoxicity/mutagenicity in eukaryotic systems, suggesting the potential of leachate to cause harmful effects to public health through seepage into groundwater. Poorly concentrated pollutants remain undetected while they are transformed and enter the food chain. Moreover, their toxicity is underestimated if cocktail effect and bioaccumulation over long-term exposure is not considered.
In the present study, we investigate the effects of exposure to environmental low-level polluted water for distinct exposure time and developmental windows, with a focus on liver toxicity in two model systems, mouse and zebrafish. Methodologically, we correlate microarray data with phenotypic and chemical parameters after short-term exposure of mice and zebrafish, and long-term exposure of mice, to environmental low-grade polluted water. Our findings are a proof that toxicogenomics applied to environmental toxicology studies permits new biomarker identification and risk assessment in the common situation of low-grade pollution and different exposure timing.
Toxicogenomic evaluation of "acute" exposure in mice
To investigate the effects of exposure to environmental polluted waters, samples were collected from dumps located upstream (U) and downstream (D) a sanitary landfill, both insist on the same aquifer. Analyses aimed at the definition of aromatic and heavy metal contents in water samples were carried out by GC-MS and ICP-MS techniques. Additional file 1 reports results representative of different samplings, underlining the similarity between U and D waters compared with control water.
The impact of acute exposure to sampled groundwater was investigated by treating 21 PND CD1 mice for 3 months. Outbred CD1 mice were chosen to avoid the influence of genetic background on any phenotypic aspect. During the treatment, no differences were recorded in water/food consumption as well as in mortality and body weight between the two treatment groups compared with control (data not shown).
Since the liver is highly sensitive to toxicant exposure, gene expression profiling analysis of the "acute response" was performed on RNAs obtained from 9 livers/group and using Affymetrix mouse whole genome. Different probe sets were retrieved in U- and D-treated mice, as shown in Volcano Plots and related tables (Additional file 2). The Heatmap and the Venn Diagram highlight many common and unique Differentially Expressed Genes highlight many common and unique Differentially Expressed Genes (DEGs) in U and D group mice (Figure 1A, B). The observation of U- and D-specific DEGs suggests that these waters were similar but not identical (Additional file 1). Common DEGs were further investigated as true markers of groundwater exposure. Using the Ingenuity Pathway Analysis program (IPA, see Materials and Methods section), we identified several biofunctions relevant to the data set, including Hepatic System Disease (p-value = 2,73E-03) and Renal and Urological Disease (p-value = 1,37E-03). The IPA analysis of canonical pathways suggested that transcriptomics-based evidence of alterations in liver activity may be related to perturbation of pathways involved in liver stress responses (Eif2 signalling), or/and to mitochondrial dysfunction (Figure 1C). In the deregulated mitochondrial pathway, we found alteration of several genes involved in the respiratory chain, mostly codifying for components of the NADH:ubiquinone oxidoreductase (complex I, Ndufaf1, Ndufs5, Ndufab1, Ndufb6, Ndufa3). The liver alteration is also supported by toxfunction prediction performed within the IPA analysis (Figure 1D). qRT-PCR analysis of genes related to the two top impinged toxfunctions confirmed microarray data and other IPA tools bioinformatics observations (Figure 1E).We analysed hepatic and renal activity in acutely treated mice by measuring serum alanine and aspartate transaminases (ALT, AST), alkaline phosphatase (AP), and urea levels. Values of urea, AP and AST were significantly increased in U and D groups (Figure 2A-C), in agreement with bioinformatics prediction of liver suffering. Conversely, ALT levels were higher only in U animals (Figure 2D). However, hematoxylin/eosin staining of liver sections showed no macroscopic signs of tissue alteration (data not shown).
Comparing mouse and zebrafish transcriptomes
The identification of deregulated genes in common between mammals and fish may elicit the identification of reliable biomarkers of environmental exposure. Due to zebrafish susceptibility to groundwater treatment, 3 months old fish were enrolled in each treatment group. After 3 months exposure, the treatment was stopped to avoid suffering since animal weight and length were reduced relative to control (Additional file 3).
Transcriptome analysis revealed more changes in gene expression in fish treated with D water (100 genes, 86 down-regulated) than in those exposed to U water (24, 20). Almost all U-specific altered probes were included in the D list (Figure 3A). Volcano plots for these experiments are reported in Additional file 3. Validation by qRT-PCR of two common DEGs is showed in Figure 3B. Gene Ontology analysis has been performed to compare the two acute exposure profilings in mice and zebrafish. 81 gene ontology terms were significantly deregulated in mice with a p-value of 0.05 (Additional file 4). In zebrafish only 3 gene ontology terms were significantly deregulated (Additional file 5). Among them the mitochondrion was shared by the two profilings. Among the 3 genes enrinched in GO mitochondrion term in zebrafish, we validated the reduction of ctp2 transcript, codifying for a protein involved in the β-oxidation of long-chain fatty acids in the mitochondrial inner membrane. Besides the unchanged expression of respiratory chain genes, lower ctp2 levels can result in a decrease in the production of ATP. Reduced cpt2expression, a PPARα-target gene, was implicated in the development of hepatic steatosis and toxicity in the livers of both zebrafish exposed to the environmental contaminant perfluorononanoic acid and high fat diet-fed mice treated with PCB 15 [38, 39]. IPA analysis of common DEGs highlighted steroid hormone biosynthesis among deregulated canonical pathways (Figure 3C) but could not predict any liver function alteration (data not shown).
Comparison of data sets from acutely exposed mouse and zebrafish revealed two molecular responses in common, i.e. the altered activity of steroid hormones synthesis/signalling biofunctions, and the prediction of aromatase (Cyp19A) as upstream regulator (p-value = 4,95E-05 in mouse, p-value = 4,68E-05 in zebrafish), in agreement with the detection of trace chemicals affecting fertility . To verify IPA prediction, estradiol synthesis was examined in both models. Expression of the vitellogenin (vtg) gene, typically estrogen-dependent, was reduced in treated zebrafish, in line with estrogen findings. Likewise, plasmatic estradiol levels in acutely exposed mice were in line with IPA-based prediction (Figure 4B). Thus, data suggest that the comparison of toxicogenomics data sets in zebrafish and mouse allows recognition of altered biofunctions.
Toxicogenomic evaluation of chronic exposure started in utero or at 21 PND
In order to analyse the effects of chronic exposure (12 months) under different exposure windows, mice treatment started at 21 PND or in utero (F1 from four months old exposed CD1 mice) (Figure 5).Transcriptome analyses conducted on livers at the end of treatment showed extreme reduction of deregulated probes in chronic exposed parents (P) (Figure 6A) and offspring (F1) (Figure 7A) compared to the acute treatment. This reduction could be interpreted as indication of a compensation effect in chronically exposed mice.
Probe sets deregulated in P mice watered with U and D water shared several DEGs (Figure 6A). Furthermore, IPA analysis recognized numerous biofunctions related to impinged liver function (e.g. delay in initiation of liver repair; data not shown). Two canonical pathways, involving the complement system, were significantly altered (Figure 6B). They were not highlighted by acute exposure data, suggesting the specificity of chronic exposure effects. For this reason, elements of these complement system canonical pathways were used to validate microarray data by qRT-PCR (Additional file 6).
Bioinformatics analysis indicated seleno-methyl-selenocysteine (p-value = 7.12E-04) and phenol derivative (p-value = 3.44E-03) as upstream regulators of chemical toxicity. Remarkably, traces of selenium and phenolic compounds were retrieved in the sampled waters (Additional file 1). Analysis by toxfunctions (liver damage, p-value = 3.7E-02; liver cholestasis, p-value = 3.7E-02; liver necrosis/cell death, p-value = 3.7E-02) evidenced predictable signs of hepatic toxicity (Additional file 6), prompting us to deepen the assessment of liver/renal functionality in P mice by assaying urea, AP, ALT and AST serum levels. We found that levels of urea and AP were significantly increased in the treated mice versus control (Additional file 6). Histopathology evidenced the alteration of hepatic tissue, with the presence of typical signs of micro-vescicular steatosis in treated mice (Additional file 6).A high number of U/D-shared deregulated probes was observed also in chronically treated F1animals (Figure 7A). IPA analysis of common DEGs identified many biofunctions related to cell cycle alteration and cancer (Figure 7B) that were not detected in the P mice analysis. Thus, they could be considered related to the foetal exposure.
In order to identify exposure-specific signatures, we looked for common DEGs deregulated by U and D treatment in both chronic and acute exposure conditions in mice. Among 30 DEGs common to both data sets, the majority codifies for mitochondrial tRNA (Table 1), in agreement with mitochondrial alteration in liver stressing conditions . 15 out of the 22 mitochondrial tRNAs were strongly up-regulated in both acute and chronic exposure; the reason for their increase is not clear. However it has been reported that hypoxia induces a similar effect in cardiomyocytes .
To identify hallmarks of groundwater chronic toxicity, we looked for genes deregulated in both P and F1 upon 12 months treatment. As described in Material and Methods and depicted in Figure 5, P and F1 animals received the same waters for most of the treatment, the only difference being the exposure window. However, only 9 deregulated probes were present in both lists, of which only two DEGs were similarly regulated (Usmg5, Nd6) (Table 2). The reported data suggest that the effects of acute vs chronic exposures as well as the effects of chronic exposure windows were mainly specific, likely as a consequence of adaptive response.
Toxicogenomics aims to identify genes whose altered expression is associated to observed phenotypes or unpredicted outcomes. Matching high-throughput transcriptional approaches to traditional toxicological criteria is considered a powerful method for testing the impact of low-dose pollutant mixtures, particularly in natural settings. Here, we evaluated low-level contaminated groundwater toxicity by both approaches in vertebrate mouse and zebrafish. Determination of polluted groundwater toxicity in the experimentally controlled animal facility reduces disturbing factors encountered with environmental sampling of wild organisms.
Bioinformatics analyses of liver transcriptomics profiles in the acute exposure treatments revealed the effect of pollutants, whose presence in groundwater was supported by chemical testing (e.g. dibutyl phthalate, p-value = 1.46E-4) (Additional file 1). The alteration of pathways involved in stress response (EIF2-, eIF4-, p70S6K- signalling pathways) suggested their relevance in the evaluation of environmental exposure. Phenotypic and molecular data obtained in the analysis of acutely exposed mice point at liver toxicity, and advise a particular role for the stress response and mitochondria related pathways in its development. Indeed, the IPA analysis on acutely treated mice highlighted, among the canonical pathways, the mitochondrial dysfunction including genes for the mitochondrial respiratory chain (Ndufaf1, Cox6b10, Ndufs5, Cox6a1, Uqcr11, Ndufab1, Ndufb6, Ndufa3). All of them are nuclear genes down-regulated in the liver of acutely exposed mice. This condition could lead to the impairment in the electron transport and ATP synthesis and, overall, could compromise the health of liver cells. Inhibition of the respiratory and impairment of complexes I as well as of mitochondrial β-oxidation has been frequently associated to acute exposure to drugs and hepatocytes toxicity . Furthermore, down-regulation of electron transport complex genes as well as mitochondrial alterations has been associated with several diseases, being chronic liver disease among them [44–49].
Among the genes whose expression level was altered in the acutely exposed mice, we found downregulation of Cdkn1a and Nfe2l2. Cdkn1a is the gene codifying for p21, which is a p53-dependent key regulator of cell fate, as it triggers cell cycle arrest in the G1 phase under various stress conditions . It has been recently suggested that reduced p21 expression can support liver fibrosis through reduced hepatic stellate cells senescence . It is relevant to note that in our experimental condition, hepatic stellate cells activation was one of the disturbed canonical pathways (p-value = 4,3E-001). As stated in the Results section, Nfe2l2 is one of the few genes deregulated also in the chronic exposed mice (Table 1), where we found evidence of micro-vescicular steatosis (Additional file 6). The transcription factor Nfe2l2 plays a protective role in hepatic cells as a key regulator for induction of detoxifying enzymes, antioxidative stress genes and several other enzymes involved in cellular protection. Its deficiency leads to exacerbation of chemical hepatotoxicity and to a considerable increase in micro and macro-vescicular steatosis .
Indeed, chronically exposed mice showed progressive alteration of liver parenchyma and increased level of alkaline phosphatase, as reported in human non-alcoholic steatohepatitis where it often evolves in hepatocellular carcinoma, a disease pathway whose alteration was seen in the present research. A compensation effect was evidenced at molecular level by the reduction of DEGs. Sets of deregulated biofunctions were almost overlapping in acute and chronic exposure (data not shown). One exception was the specific alteration of the complement system. C1qa, c1qb and c1qc are members of the complement cascade, a fundamental component in the innate immunity helping in the clearance of pathogens. Recent studies demonstrated that complement activation may contribute to cancer growth by facilitating the dysregulation of mitogenic signaling pathways . The three complement system genes were up-regulated in the liver of chronically exposed mice; in particular, overexpression of C1q complement system components may be envisaged as indicative of hetapotoxicity since its transcription in liver hematopoietic cells was induced in several types of damage [54, 55].
The possibility that exposure during foetal life could be more harmful than in the adulthood has been explored [33, 34, 56]. Here, our results confirm the selective impact of the exposure window in DEG changes, as only few genes (Usmg5 and Nd6) were shared among parents and offspring. Usmg5 was first recognized as a gene whose mRNA level increased during skeletal muscle growth in rats . USMG5 protein was shown to be associated with ATP synthase in the mitochondria  and to contribute to the maintenance of ATP synthase population, an indication of its importance in cellular energy metabolism . ND6 is one of the 7 mitochondrial DNA (mtDNA) encoded subunits of respiratory Complex I, again suggesting the alteration of mitochondria biology in the exposed animals also in chronic exposure conditions.
DEGs involved in cancer development (cell cycle and DNA repair) were specifically observed in offspring. Even if cancers were not observed in the collected organs, genotype context or longer exposure time could play a major role in their development. The identified genes could represent a predictive biomarker of exposure and effects whose identification is pivotal in monitoring. Interestingly, increased incidence of cancers has been reported in people leaving in landfill areas, where they are potentially exposed to polluted groundwater used for agriculture or other human activities. Further work is necessary to explore this hypothesis.
In toxicogenomics, studies of gene expression alteration in phylogenetically distant species are supposed to assist in the identification of altered pathways that impinge on human health. To this aim, we investigated the effects of polluted groundwater exposure on mice and zebrafish liver by toxicogenomics. Even if no common DEGs were retrieved in both vertebrates, bioinformatics analysis and phenotypic assays indicate aromatase and estrogen biosynthesis changes. These findings agreed with the suggestion that functional pathways could be better and reliable markers than single genes. This observation is consistent with the documented effect of several water pollutants (metals, plasticizers and others) on hormonal axes linked to reproductive health outcomes [60–64], as seen also in our experimental settings (manuscript in preparation). Here, toxicogenomics showed its predictive strength in identifying biomarkers, steroid hormones, that can be easily detected in blood, suitable for monitoring the impact of anthropogenic activities on human and identifying risks for human populations even if not indicative of specific pollutants.
In all treatment conditions, observed phenotypic effects can correlate to alterations in the stress response pathways and, mainly, in mitochondria activities, as supported by the ability of leachates to induce oxidative stress in organs [65, 66]. These pathways could be effectors of low-dose mixture toxicity, playing a role in the alteration of liver function and in the regulation of hormones biosynthesis, as previously shown .
The reported data suggest that toxicogenomics analyses of different animal models and exposure conditions (timing and window) have the potential to disclose unpredicted outcomes and, most importantly, pathways useful for human and environmental health risk assessment also in conditions of low-level exposure. Our study points to the importance of considering pathways more than single gene alteration in toxicological assessment. It also allows the identification of specific pathways, with mitochondria as key factors of toxicity response, and sexual hormones as biomarker in environmental assessment.
Water sampling was performed in a landfill located in southern Italy (Lat 41° 12′ 23,54′ ‘N; Long 15° 12′ 27,63′ ‘E, datum WGS84), with the permission of the municipal local authority. The present study did not involve endangered or protected species.
Water samples were collected from five 15 meters deep piezometers that are located in the saturated portion of an aquifer consisting of poorly permeable lithologies. Piezometers are situated transverse to the axis of groundwater flow: two upstream and three downstream the landfill (referring to the water flow direction). The sampling occurred every 15 days under dynamic conditions by means of a micro pump.
Aliquots of water solution (1 ml) from each sample were directly analysed by Inductively Coupled Plasma Mass Spectometry (ICP-MS) with an Agilent 7700 ICP-MS (Agilent Technologies), as described in detail in Additional file 1. Upon liquid-liquid extraction with chloroform, hexane and dichloromethane (1:1 v/v, performed 3 times), organic substances were collected, dried and dissolved in 200 μl of hexane, before GC-MS analysis performed on a 5390 MSD quadrupole mass spectrometer (Agilent Technologies), as detailed in Additional file 7.
Animals and treatments
Groundwater used for animal treatment was filtered through glass and 0.22 μm (Millipore) filters and stored at 4°C.
Animal experiments were performed in compliance with the European Council Directive 86/609/EEC and the Italian Legislation on Animal Experimentation (D.Lvo 116/92) and procedures were approved (ID number 21-2009) by the Ethical committee named CESA (Committee for the Ethics of the Experimentations on Animals) of the Biogem Institute of Genetics Research "Gaetano Salvatore" (IRGS). The project has been communicated to the employed office of the Ministry of Health following the rules of the D.Lvo 116/92.
Mice were kept under standard facility conditions (22 ± 2°C, 55 + 10% humidity, 12:12 h light-dark cycle) in a specific pathogen-free facility. Animals received water and standard diet (4RF21 form Mucedola) "ad libitum"; type II EU cages, in polysulfone, had space to allow motility and parental care. Mice were sacrificed by carbon dioxide inhalation.Fifty 21 Post Natal Days (PND) CD1 mice (outbred, 30 male and 30 female) were randomly recruited for each of the three F0 treatment groups that received drinking tap-water (control, Ctrl), waters sampled from dumps upstream (U) or downstream (D) the landfill. Sampling from two different locations around the landfill area were undertaken to perform a robust toxicogenomic analyses, based on different (up and down) data sets whose shared elements were used for bioinformatics analyses. We chose this approach in order to eliminate from the analyses genes whose deregulation could be related to other reasons (noise). After 3 months of treatment (acute exposure, P generation), twenty mice/group were sacrificed, blood and several organs collected. Remaining animals were mated for F1 generation and continuously treated till sacrifice (12 months, chronic exposure, P generation) for blood and organ sampling. Thirty F1 mice per group were treated for 12 months before sacrifice. Schematic representation of the mice study design is depicted in Figure 5.
Fish were maintained in 50 l sterile glass tanks with 20 animals each, under standard laboratory conditions (28 ± 0.5°C, 14:10 h light: dark cycle). The water was continuously aerated, filtered and 1/3 of its volume was manually renewed every three days. Ammonia, nitrate, nitrite, pH and water hardness was monitored twice a week using commercial kit (Tetra GmbH). Fish were daily fed freshly hatched Artemia nauplii once and granular food (Special Diets Services) twice. Sixty 3 months old zebrafish (AB strain, 30 males and 30 females) were randomly recruited from each Ctrl, U and D treatment group, distributed in three replicate tanks for condition. After 3 months treatment, fish were sacrificed by overdose of tricaine methane sulfonate (MS222, 200-300 mg/l) by prolonged immersion and their lengths and weight recorded before organ dissection.
Serum prepared from collected blood samples was -80°C frozen until assayed by Estradiol EIA kit (Cayman) following manufacturer’s instructions. Each sample was analyzed in triplicates. p-value was calculated by t-student test.
RNA extraction and real-time RT-PCR
RNA was isolated from mouse and zebrafish liver tissues using Trizol reagent (Invitrogen) and purified by RNeasy mini kit (Qiagen). cDNA synthesis and qRT-PCR analysis were set for 8 animals, each carried out in triplicate as previously reported . Data obtained were normalized on the relative expression of reference gene Tubulin in mouse and rpl13a in zebrafish and reported as ratio between U or D vs Ctrl expression values. Primer sequences are reported in Table 3. p-value was calculated by t-student test.
Microarray and bioinformatic analysis
Hepatic RNA for microarray was extracted and cleaned as reported above. Nine animals per group were randomly divided into three equivalent sets, and a constant amount of RNA from animals in the same set was pooled into one single sample in order to eliminate individual differences within group. Three arrays were used for each group. cRNA was generated by using the Affymetrix One-Cycle Target Labeling and Control Reagent kit (AffymetrixInc), following the manufacturer’s protocol, starting from 5 μg of total RNA. Biotinylated cRNA was hybridized to the GeneChip Mouse Gene 1.0 ST Array [(MoGene-1_0-st-v1, Affymetrix) and to GeneChip Zebrafish Gene 1.0 ST Arrays (Affymetrix). Chips were washed and scanned on the Affymetrix Complete GeneChip System, generating digitized image data (DAT) files.
The datasets obtained were analyzed with GeneSpring GX 12 Software (Agilent Technologies). Robust multichip average (RMA) algorithm [67, 68] was used for summarization and normalization. Hybridization quality was assessed by spiked-in controls. Principal Component Analysis (PCA) was performed to check data quality that resulted adequate for all samples. Transcripts were filtered by their signal intensity values, selecting transcripts with intensity values between 20 and 100 percentile in at least 1 out of each set samples for differential analysis.
Differentially expressed transcripts between exposed livers vs controls were filtered for absolute fold change ≥ 1.5 and corrected p-value ≤ 0.05. Statistical analysis was performed using Oneway ANOVA adjusted for multiple comparison by the Benjamini-Hochberg method.
Functional annotation for differentially expressed transcripts was performed using Ingenuity Pathway Analysis (IPA; http://www.ingenuity.com), a web-based tool for the identification of biological functions, canonical pathways, transcription factors as well as toxofunctions that are most significant to the dataset. Fisher Exact test was used to calculate the p-value determining the likelihood that the association between the set of focus genes in the dataset and a given process or pathway is due to chance alone. Corrected p-value calculation (based on the Benjamini-Hochberg method) controlled the error rate in analysis results and focus in on the most significant biological functions associated with DEG. Gene ontology (GO) analysis was performed using David software (http://david.abcc.ncifcrf.gov) [69, 70]. The listed GO terms included three or more DEGs p-value <0.05.
Availability of supporting data
Microarray data are available in the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress) under accession numbers E-MTAB-2905 (for mouse arrays) and E-MTAB-2906 (for zebrafish arrays).
Quantitative real time PCR
Deregulated expressed genes
Ingenuity pathway analysis.
Wei X, Ai J, Deng Y, Guan X, Johnson DR, Ang CY, Zhang C, Perkins EJ: Identification of biomarkers that distinguish chemical contaminants based on gene expression profiles. BMC Genomics. 2014, 15: 248-10.1186/1471-2164-15-248.
Hamadeh HK, Bushel PR, Jayadev S, Martin K, DiSorbo O, Sieber S, Bennett L, Tennant R, Stoll R, Barrett JC, Blanchard K, Paules RS, Afshari CA: Gene expression analysis reveals chemical-specific profiles. Toxicol Sci. 2002, 67 (2): 219-231. 10.1093/toxsci/67.2.219.
Shukla SJ, Huang R, Austin CP, Xia M: The future of toxicity testing: a focus on in vitro methods using a quantitative high-throughput screening platform. Drug Discov Today. 2010, 15 (23–24): 997-1007.
Aardema MJ, MacGregor JT: Toxicology and genetic toxicology in the new era of "toxicogenomics": impact of "-omics" technologies. Mutat Res. 2002, 499 (1): 13-25. 10.1016/S0027-5107(01)00292-5.
Suter L, Babiss LE, Wheeldon EB: Toxicogenomics in predictive toxicology in drug development. Chem Biol. 2004, 11 (2): 161-171. 10.1016/S1074-5521(04)00031-6.
Afshari CA, Nuwaysir EF, Barrett JC: Application of complementary DNA microarray technology to carcinogen identification, toxicology, and drug safety evaluation. Cancer Res. 1999, 59 (19): 4759-4760.
Hu YF, Kaplow J, He Y: From traditional biomarkers to transcriptome analysis in drug development. Curr Mol Med. 2005, 5 (1): 29-38. 10.2174/1566524053152915.
Hook SE, Skillman AD, Gopalan B, Small JA, Schultz IR: Gene expression profiles in rainbow trout, Onchorynchus mykiss, exposed to a simple chemical mixture. Toxicol Sci. 2008, 102 (1): 42-60.
Hook SE, Skillman AD, Small JA, Schultz IR: Temporal changes in gene expression in rainbow trout exposed to ethynyl estradiol. Comp Biochem Physiol C Toxicol Pharmacol. 2007, 145 (1): 73-85. 10.1016/j.cbpc.2006.10.011.
Zhang Y, Huang K, Deng Y, Zhao Y, Wu B, Xu K, Ren H: Evaluation of the toxic effects of municipal wastewater effluent on mice using omic approaches. Environ Sci Technol. 2013, 47 (16): 9470-9477. 10.1021/es401615y.
Kosmehl T, Otte JC, Yang L, Legradi J, Bluhm K, Zinsmeister C, Keiter SH, Reifferscheid G, Manz W, Braunbeck T, Strähle U, Hollert H: A combined DNA-microarray and mechanism-specific toxicity approach with zebrafish embryos to investigate the pollution of river sediments. Reprod Toxicol. 2012, 33 (2): 245-253. 10.1016/j.reprotox.2012.01.005.
Lettieri T: Recent applications of DNA microarray technology to toxicology and ecotoxicology. Environ Health Perspect. 2006, 114 (1): 4-9.
Williams TD, Diab A, Ortega F, Sabine VS, Godfrey RE, Falciani F, Chipman JK, George SG: Transcriptomic responses of European flounder (Platichthys flesus) to model toxicants. Aquat Toxicol. 2008, 90 (2): 83-91. 10.1016/j.aquatox.2008.07.019.
Williams TD, Diab AM, Gubbins M, Collins C, Matejusova I, Kerr R, Chipman JK, Kuiper R, Vethaak AD, George SG: Transcriptomic responses of European flounder (Platichthys flesus) liver to a brominated flame retardant mixture. Aquat Toxicol. 2013, 142-143C: 45-52.
Falciani F, Diab AM, Sabine V, Williams TD, Ortega F, George SG, Chipman JK: Hepatic transcriptomic profiles of European flounder (Platichthys flesus) from field sites and computational approaches to predict site from stress gene responses following exposure to model toxicants. Aquat Toxicol. 2008, 90 (2): 92-101. 10.1016/j.aquatox.2008.07.020.
Williams TD, Davies IM, Wu H, Diab AM, Webster L, Viant MR, Chipman JK, Leaver MJ, George SG, Moffat CF, Robinson CD: Molecular responses of European flounder (Platichthys flesus) chronically exposed to contaminated estuarine sediments. Chemosphere. 2014, 108: 152-158.
Richter CA, Martyniuk CJ, Annis ML, Brumbaugh WG, Chasar LC, Denslow ND, Tillitt DE: Methylmercury-induced changes in gene transcription associated with neuroendocrine disruption in largemouth bass (Micropterus salmoides). Gen Comp Endocrinol. 2014, 203: 215-224.
Williams TD, Turan N, Diab AM, Wu H, Mackenzie C, Bartie KL, Hrydziuszko O, Lyons BP, Stentiford GD, Herbert JM, Abraham JK, Katsiadaki I, Leaver MJ, Taggart JB, George SG, Viant MR, Chipman KJ, Falciani F: Towards a system level understanding of non-model organisms sampled from the environment: a network biology approach. PLoS Comput Biol. 2011, 7 (8): e1002126-10.1371/journal.pcbi.1002126.
Christiansen HE, Mehinto AC, Yu F, Perry RW, Denslow ND, Maule AG, Mesa MG: Correlation of gene expression and contaminant concentrations in wild largescale suckers: a field-based study. Sci Total Environ. 2014, 484: 379-389.
Baker ME, Vidal-Dorsch DE, Ribecco C, Sprague LJ, Angert M, Lekmine N, Ludka C, Martella A, Ricciardelli E, Bay SM, Gully JR, Kelley KM, Schlenk D, Carnevali O, Šášik R, Hardiman G: Molecular analysis of endocrine disruption in hornyhead turbot at wastewater outfalls in southern california using a second generation multi-species microarray. PLoS One. 2013, 8 (9): e75553-10.1371/journal.pone.0075553.
Xu H, Lam SH, Shen Y, Gong Z: Genome-wide identification of molecular pathways and biomarkers in response to arsenic exposure in zebrafish liver. PLoS One. 2013, 8 (7): e68737-10.1371/journal.pone.0068737.
Antherieu S, Le Guillou D, Coulouarn C, Begriche K, Trak-Smayra V, Martinais S, Porceddu M, Robin MA, Fromenty B: Chronic exposure to low doses of pharmaceuticals disturbs the hepatic expression of circadian genes in lean and obese mice. Toxicol Appl Pharmacol. 2014, 276 (1): 63-72. 10.1016/j.taap.2014.01.019.
Iguchi T, Watanabe H, Katsu Y: Application of ecotoxicogenomics for studying endocrine disruption in vertebrates and invertebrates. Environ Health Perspect. 2006, 114 (Suppl 1): 101-105.
Mattingly CJ, Rosenstein MC, Colby GT, Forrest JN, Boyer JL: The Comparative Toxicogenomics Database (CTD): a resource for comparative toxicological studies. J Exp Zool A Comp Exp Biol. 2006, 305 (9): 689-692.
Hu T, Gibson DP, Carr GJ, Torontali SM, Tiesman JP, Chaney JG, Aardema MJ: Identification of a gene expression profile that discriminates indirect-acting genotoxins from direct-acting genotoxins. Mutat Res. 2004, 549 (1–2): 5-27.
Yang L, Kemadjou JR, Zinsmeister C, Bauer M, Legradi J, Muller F, Pankratz M, Jakel J, Strahle U: Transcriptional profiling reveals barcode-like toxicogenomic responses in the zebrafish embryo. Genome Biol. 2007, 8 (10): R227-10.1186/gb-2007-8-10-r227.
Lam SH, Mathavan S, Tong Y, Li H, Karuturi RK, Wu Y, Vega VB, Liu ET, Gong Z: Zebrafish whole-adult-organism chemogenomics for large-scale predictive and discovery chemical biology. PLoS Genet. 2008, 4 (7): e1000121-10.1371/journal.pgen.1000121.
Dai YJ, Jia YF, Chen N, Bian WP, Li QK, Ma YB, Chen YL, Pei DS: Zebrafish as a model system to study toxicology. Environ Toxicol Chem. 2014, 33 (1): 11-17. 10.1002/etc.2406.
Vrijheid M: Health effects of residence near hazardous waste landfill sites: a review of epidemiologic literature. Environ Health Perspect. 2000, 108 (Suppl 1): 101-112.
Palmer SR, Dunstan FD, Fielder H, Fone DL, Higgs G, Senior ML: Risk of congenital anomalies after the opening of landfill sites. Environ Health Perspect. 2005, 113 (10): 1362-1365. 10.1289/ehp.7487.
Mattiello A, Chiodini P, Bianco E, Forgione N, Flammia I, Gallo C, Pizzuti R, Panico S: Health effects associated with the disposal of solid waste in landfills and incinerators in populations living in surrounding areas: a systematic review. Int J Public Health. 2014, 58 (5): 725-735.
Porta D, Milani S, Lazzarino AI, Perucci CA, Forastiere F: Systematic review of epidemiological studies on health effects associated with management of solid waste. Environ Health. 2009, 8: 60-10.1186/1476-069X-8-60.
Makri A, Goveia M, Balbus J, Parkin R: Children’s susceptibility to chemicals: a review by developmental stage. J Toxicol Environ Health B Crit Rev. 2004, 7 (6): 417-435. 10.1080/10937400490512465.
Kishi R, Kobayashi S, Ikeno T, Araki A, Miyashita C, Itoh S, Sasaki S, Okada E, Kashino I, Itoh K, Nakajima S: Ten years of progress in the Hokkaido birth cohort study on environment and children’s health: cohort profile–updated 2013. Environ Health Prev Med. 2013, 18 (6): 429-450. 10.1007/s12199-013-0357-3.
Bernard C, Guido P, Colin J, Le Du-Delepierre A: Estimation of the hazard of landfills through toxicity testing of leachates—I: determination of leachate toxicity with a battery of acute tests. Chemosphere. 1996, 33 (11): 2303-2320. 10.1016/0045-6535(96)00322-0.
Baun A, Ledin A, Reitzel LA, Bjerg PL, Christensen TH: Xenobiotic organic compounds in leachates from ten Danish MSW landfills–chemical analysis and toxicity tests. Water Res. 2004, 38 (18): 3845-3858. 10.1016/j.watres.2004.07.006.
Kjeldsen P, Barlaz M, Rooker A, Baun A, Ledin A, Christensen CH: Present and long-term composition of MSW landfill leachate: a review. Crit Rev Environ Sci Technol. 2002, 32: 297-336. 10.1080/10643380290813462.
Zhang W, Liu Y, Zhang H, Dai J: Proteomic analysis of male zebrafish livers chronically exposed to perfluorononanoic acid. Environ Int. 2012, 42: 20-30.
Wahlang B, Falkner KC, Gregory B, Ansert D, Young D, Conklin DJ, Bhatnagar A, McClain CJ, Cave M: Polychlorinated biphenyl 153 is a diet-dependent obesogen that worsens nonalcoholic fatty liver disease in male C57BL6/J mice. J Nutr Biochem. 2013, 24 (9): 1587-1595. 10.1016/j.jnutbio.2013.01.009.
Tan SW, Meiller JC, Mahaffey KR: The endocrine effects of mercury in humans and wildlife. Crit Rev Toxicol. 2009, 39 (3): 228-269. 10.1080/10408440802233259.
Jaeschke H, McGill MR, Ramachandran A: Oxidant stress, mitochondria, and cell death mechanisms in drug-induced liver injury: lessons learned from acetaminophen hepatotoxicity. Drug Metab Rev. 2012, 44 (1): 88-106. 10.3109/03602532.2011.602688.
Jian B, Wang D, Chen D, Voss J, Chaudry I, Raju R: Hypoxia-induced alteration of mitochondrial genes in cardiomyocytes: role of Bnip3 and Pdk1. Shock. 2010, 34 (2): 169-175. 10.1097/SHK.0b013e3181cffe7d.
Pessayre D, Fromenty B, Berson A, Robin MA, Letteron P, Moreau R, Mansouri A: Central role of mitochondria in drug-induced liver injury. Drug Metab Rev. 2012, 44 (1): 34-87. 10.3109/03602532.2011.604086.
Grattagliano I, Russmann S, Diogo C, Bonfrate L, Oliveira PJ, Wang DQ, Portincasa P: Mitochondria in chronic liver disease. Curr Drug Targets. 2011, 12 (6): 879-893. 10.2174/138945011795528877.
Sun X, Wang JF, Tseng M, Young LT: Downregulation in components of the mitochondrial electron transport chain in the postmortem frontal cortex of subjects with bipolar disorder. J Psychiatry Neurosci. 2006, 31 (3): 189-196.
Anitha A, Nakamura K, Thanseem I, Matsuzaki H, Miyachi T, Tsujii M, Iwata Y, Suzuki K, Sugiyama T, Mori N: Downregulation of the expression of mitochondrial electron transport complex genes in autism brains. Brain Pathol. 2012, 23 (3): 294-302.
Dahlman I, Forsgren M, Sjogren A, Nordstrom EA, Kaaman M, Naslund E, Attersand A, Arner P: Downregulation of electron transport chain genes in visceral adipose tissue in type 2 diabetes independent of obesity and possibly involving tumor necrosis factor-alpha. Diabetes. 2006, 55 (6): 1792-1799. 10.2337/db05-1421.
Soronen J, Laurila PP, Naukkarinen J, Surakka I, Ripatti S, Jauhiainen M, Olkkonen VM, Yki-Jarvinen H: Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects. BMC Med Genomics. 2012, 5: 9-10.1186/1755-8794-5-9.
Olsson AH, Yang BT, Hall E, Taneera J, Salehi A, Nitert MD, Ling C: Decreased expression of genes involved in oxidative phosphorylation in human pancreatic islets from patients with type 2 diabetes. Eur J Endocrinol. 2011, 165 (4): 589-595. 10.1530/EJE-11-0282.
Warfel NA, El-Deiry WS: p21WAF1 and tumourigenesis: 20 years after. Curr Opin Oncol. 2012, 25 (1): 52-58.
Valenti L, Dongiovanni P: CDKN1A: a double-edged sword in fatty liver?. Cell Cycle. 2014, 13 (9): 1371-1372. 10.4161/cc.28754.
Shin SM, Yang JH, Ki SH: Role of the Nrf2-ARE pathway in liver diseases. Oxid Med Cell Longev. 2013, 2013: 763257-
Rutkowski MJ, Sughrue ME, Kane AJ, Mills SA, Parsa AT: Cancer and the complement cascade. Mol Cancer Res. 2010, 8 (11): 1453-1465. 10.1158/1541-7786.MCR-10-0225.
Rensen SS, Slaats Y, Driessen A, Peutz-Kootstra CJ, Nijhuis J, Steffensen R, Greve JW, Buurman WA: Activation of the complement system in human nonalcoholic fatty liver disease. Hepatology. 2009, 50 (6): 1809-1817. 10.1002/hep.23228.
Cohen JI, Roychowdhury S, McMullen MR, Stavitsky AB, Nagy LE: Complement and alcoholic liver disease: role of C1q in the pathogenesis of ethanol-induced liver injury in mice. Gastroenterology. 2010, 139 (2): 664-674. 10.1053/j.gastro.2010.04.041. 674 e661
Barouki R, Gluckman PD, Grandjean P, Hanson M, Heindel JJ: Developmental origins of non-communicable disease: implications for research and public health. Environ Health. 2012, 11: 42-10.1186/1476-069X-11-42.
Paivarinne H, Kainulainen H: DAPIT, a novel protein down-regulated in insulin-sensitive tissues in streptozotocin-induced diabetes. Acta Diabetol. 2001, 38 (2): 83-86. 10.1007/s005920170018.
Chen R, Runswick MJ, Carroll J, Fearnley IM, Walker JE: Association of two proteolipids of unknown function with ATP synthase from bovine heart mitochondria. FEBS Lett. 2007, 581 (17): 3145-3148. 10.1016/j.febslet.2007.05.079.
Ohsakaya S, Fujikawa M, Hisabori T, Yoshida M: Knockdown of DAPIT (diabetes-associated protein in insulin-sensitive tissue) results in loss of ATP synthase in mitochondria. J Biol Chem. 2011, 286 (23): 20292-20296. 10.1074/jbc.M110.198523.
Krogenaes AK, Ropstad E, Gutleb AC, Hardnes N, Berg V, Dahl E, Fowler PA: In utero exposure to environmentally relevant concentrations of PCB 153 and PCB 118 disrupts fetal testis development in sheep. J Toxicol Environ Health A. 2014, 77 (9–11): 628-649.
Palus J, Rydzynski K, Dziubaltowska E, Wyszynska K, Natarajan AT, Nilsson R: Genotoxic effects of occupational exposure to lead and cadmium. Mutat Res. 2003, 540 (1): 19-28. 10.1016/S1383-5718(03)00167-0.
Yoon K, Kwack SJ, Kim HS, Lee BM: Estrogenic endocrine-disrupting chemicals: molecular mechanisms of actions on putative human diseases. J Toxicol Environ Health B Crit Rev. 2014, 17 (3): 127-174. 10.1080/10937404.2014.882194.
Ye T, Kang M, Huang Q, Fang C, Chen Y, Shen H, Dong S: Exposure to DEHP and MEHP from hatching to adulthood causes reproductive dysfunction and endocrine disruption in marine medaka (Oryzias melastigma). Aquat Toxicol. 2014, 146: 115-126.
Vested A, Ramlau-Hansen CH, Olsen SF, Bonde JP, Kristensen SL, Halldorsson TI, Becher G, Haug LS, Ernst EH, Toft G: Associations of in utero exposure to perfluorinated alkyl acids with human semen quality and reproductive hormones in adult men. Environ Health Perspect. 2013, 121 (4): 453-458. 458e451-455
Bertoldi K, Spindler C, Dos Santos MF, Vanzella C, Lovatel GA, Elsner VR, Rodrigues MA, Siqueira IR: Effect of landfill leachate on oxidative stress of brain structures and liver from rodents: modulation by photoelectrooxidation process. Ecotoxicol Environ Saf. 2012, 84: 319-324.
Bakare AA, Patel S, Pandey AK, Bajpayee M, Dhawan A: DNA and oxidative damage induced in somatic organs and tissues of mouse by municipal sludge leachate. Toxicol Ind Health. 2012, 28 (7): 614-623. 10.1177/0748233711420466.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.
Irizarry RA, Ooi SL, Wu Z, Boeke JD: Use of mixture models in a microarray-based screening procedure for detecting differentially represented yeast mutants. Stat Appl Genet Mol Biol. 2003, 2: Article1-
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4 (1): 44-57.
Huang DW, Sherman BT, Lempicki RA: Bionformatics enrichment tools: paths toward the comprehensive fuctional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.
We thank G. Mercadante and animal facility personnel for technical support. Finally we thank IRGS, Biogem for the support. We thank Dr. Rachele Zair for the careful reading and editing of the manuscript. This work was funded by Italian Workers’ Compensation Authority grant (E95E10000510005) and by IRGS Biogem Institute. PS was funded by MIUR (PONa3_0239 and FIRB RBFR12QW4I). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors have declared that they have no competing interests.
PS, MDF, MM and CA conceived and designed the experiments. IP, DG and EC carried out organ sampling, molecular and histological analysis. IP performed the fish exposure and organ sampling, and drafted the manuscript. An Am and FP performed the water analysis. An Af and MS participated in the mouse exposure. LE and FMG carried out the geochemical evaluation of the area. AN and SM performed the water sampling. PDL performed the microarray experiments. FDA and MC performed the bioinformatics analysis. PS, MDF, CA, MM and FMG provided reagents. CA, PS and MM prepared the final version of the manuscript. All authors read and approved the final manuscript.
Mario De Felice and Concetta Ambrosino contributed equally to this work.
Electronic supplementary material
Additional file 2: Transcriptomic analyses in liver of acutely exposed mice. Volcano plots of microarray data in U (A) and D (B) compared to Ctrl treated animals. The y-axis value is the negative logarithm base 10 of the corrected p-value. A green horizontal line on the plot represents the user-defined significant threshold for p-value. The x-axis is shown as the logarithm base 2 of the fold change in expression level between treated and control livers. The vertical green lines on the plot represent the user-defined thresholds for fold change. Red dots are up-regulated probes, green dots down-regulated probes. The number of down/up-regulated probes for each Volcano plot is reported in the underlying table. (TIFF 2 MB)
Additional file 3: Zebrafish phenotypic and molecular changes induced by groundwater treatment. (A) Body weight and (B) body length of adult zebrafish exposed for 3 months to Ctrl, U or D waters. Data are reported separately for female and male fish and, for each point, 20 animals were recorded. Volcano plots of microarray data in U (C) and D (D) compared to Ctrl treated animals. For plot description see the caption to the Additional file 1. The number of down/up-regulated probes for each Volcano plot is reported in the underlying table. (TIFF 3 MB)
Additional file 6: Liver molecular and phenotypic characterization of chronic exposed mice. (A) qRT-PCR validation of selected genes (C1qa, C1qb, C1qc) deregulated in chronically exposed livers. Data are reported as fold change value calculated as ratio between average expression in U/D and in Ctrl exposed animals. 8 animals were analysed for each group. (B) Toxfunctions deregulated in U and D chronically exposed animal livers. Serum Urea (C), AP (D), AST (E) and ALT (F) levels measured in mice treated with Ctrl (black), U (orange) or D (blue). The analysis was conducted on 20 animals per treatment groups. Each sign is a single mouse. For each treatment group mean and standard deviation is reported. (G) Haematoxylin/eosin staining of liver sections of Ctrl, U and D chronically exposed mice. *p-value ≤0,05, **p-value ≤0,01. (TIFF 10 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Porreca, I., D’Angelo, F., Gentilcore, D. et al. Cross-species toxicogenomic analyses and phenotypic anchoring in response to groundwater low-level pollution. BMC Genomics 15, 1067 (2014). https://doi.org/10.1186/1471-2164-15-1067
- Chemical mixture
- Groundwater pollution
- Animal models