Mercury-induced hepatotoxicity in zebrafish: in vivo mechanistic insights from transcriptome analysis, phenotype anchoring and targeted gene expression validation
- Choong Yong Ung†1,
- Siew Hong Lam†1Email author,
- Mya Myintzu Hlaing1,
- Cecilia Lanny Winata1,
- Svetlana Korzh1,
- Sinnakaruppan Mathavan2 and
- Zhiyuan Gong1Email author
© Ung et al; licensee BioMed Central Ltd. 2010
Received: 10 September 2009
Accepted: 30 March 2010
Published: 30 March 2010
Mercury is a prominent environmental contaminant that causes detrimental effects to human health. Although the liver has been known to be a main target organ, there is limited information on in vivo molecular mechanism of mercury-induced toxicity in the liver. By using transcriptome analysis, phenotypic anchoring and validation of targeted gene expression in zebrafish, mercury-induced hepatotoxicity was investigated and a number of perturbed cellular processes were identified and compared with those captured in the in vitro human cell line studies.
Hepato-transcriptome analysis of mercury-exposed zebrafish revealed that the earliest deregulated genes were associated with electron transport chain, mitochondrial fatty acid beta-oxidation, nuclear receptor signaling and apoptotic pathway, followed by complement system and proteasome pathway, and thereafter DNA damage, hypoxia, Wnt signaling, fatty acid synthesis, gluconeogenesis, cell cycle and motility. Comparative meta-analysis of microarray data between zebrafish liver and human HepG2 cells exposed to mercury identified some common toxicological effects of mercury-induced hepatotoxicity in both models. Histological analyses of liver from mercury-exposed fish revealed morphological changes of liver parenchyma, decreased nucleated cell count, increased lipid vesicles, glycogen and apoptotic bodies, thus providing phenotypic evidence for anchoring of the transcriptome analysis. Validation of targeted gene expression confirmed deregulated gene-pathways from enrichment analysis. Some of these genes responding to low concentrations of mercury may serve as toxicogenomic-based markers for detection and health risk assessment of environmental mercury contaminations.
Mercury-induced hepatotoxicity was triggered by oxidative stresses, intrinsic apoptotic pathway, deregulation of nuclear receptor and kinase activities including Gsk3 that deregulates Wnt signaling pathway, gluconeogenesis, and adipogenesis, leading to mitochondrial dysfunction, endocrine disruption and metabolic disorders. This study provides important mechanistic insights into mercury-induced liver toxicity in a whole-animal physiology context, which will help in understanding the syndromes caused by mercury poisoning. The molecular conservation of mercury-induced hepatotoxicity between zebrafish and human cell line reveals the feasibility of using zebrafish to model molecular toxicity in human for toxicant risk assessments.
Mercury is a highly hazardous pollutant with an estimated global natural mercury emission of 1,800-5,800 tons per annum  and global anthropogenic mercury emission to the atmosphere was estimated to be 2,190 tons in 2000 . Since mercury is ubiquitous in the environment, it is nearly impossible for most humans and animals to avoid exposure to some forms of mercury, be it elementary, organic or inorganic. Mercury, through biotransformation and bioaccumulation, has found its way through the food chain to humans. Anthropogenic activities and industrialization are also sources of mercury pollution that had resulted in several catastrophe of mercury poisoning in Japan , the Amazon basin  and Iraq , and recently many high-risk sites have been identified in Asia . Moreover, mercury pollution and poisoning have imposed a huge economic cost on environmental remediation and public health [6, 7].
All forms of mercury cause toxic effects in a number of tissues and organs depending on the chemical form of mercury as well as the level, duration and the route of exposure [8–11]. Exposure to mercury compounds typically occurs by inhalation or ingestion. Ingested mercury is absorbed in the gastrointestinal tract and it is distributed to all tissues in about 30 h while inhaled through mercury vapor accumulates in red blood cells and is carried to all tissues in the body in less than 24 h . Mercury undergoes extensive biliary-hepatic cycling . It is secreted into bile and partly reabsorbed into the portal circulation and thereby returned to the liver. The high mobility of mercury in the body is attributed to the formation of water-soluble mercury complexes that are mainly, if not exclusively, attached to the sulfur atom of thiol groups such as glutathione . The glutathione moiety is degraded in the bile duct and gall bladder to a dipeptide and finally to an L-cysteine mercury complex before entering the circulatory system . It is in the form of the L-cysteine complex that mercury enters the endothelial cells of the blood-brain and placental barriers. Mercury can gradually accumulate in the central nervous system and kidney, thus causing damage to these organs. Hence, there have been several molecular mechanistic studies on neurological  and renal  toxicities induced by mercury. However, despite of the fact that the extensive biliary-hepatic cycling of mercury  and of some evidence suggesting that liver plays a role in renal tubular uptake of mercury , little is known with regard to the mechanism of mercury-induced hepatotoxicity. Based on the role of liver in mercury biotransforming and cycling, as well as its central roles in the control and synthesis of critical blood constituents that affect whole body physiology, we hypothesized that mercury exposure can cause liver injury that may lead to other syndromes. In view of the prevalence of mercury exposure in humans, understanding of the mechanism of mercury-induced hepatotoxicity is therefore important for elucidating its impact on liver health.
The advent of omics technologies such as genomics and transcriptomics has enabled the simultaneous assessment of expression profiles of thousands of genes that respond to a toxic compound within a particular cell type, tissue, or an organism . The zebrafish is an increasingly popular model not only for vertebrate development  but also for understanding human diseases  and toxicology . Many physiological and molecular similarities in xenobiotic metabolism and adaptive response to toxicant insults have been found between zebrafish and mammals, making the zebrafish an ideal toxicology model in toxicology [19, 20]. Although zebrafish has become popular as an in vivo system to model human diseases and toxicities, toxicogenomic researches using zebrafish seem not to have been extensively conducted. We and others have shown that zebrafish responded biologically to chemicals, such as small molecules, drugs and environmental toxicants, in a similar manner as mammals [21–24]. The availability of zebrafish in large numbers, its small size and easy husbandry makes zebrafish a cost-effective model when compared to the rodent model for toxicological studies. In addition, the amenability of the zebrafish system to various molecular techniques and the vast genomic resources, including the near-completed zebrafish genome project and available zebrafish microarrays, make it a highly versatile system for toxicogenomic studies. Moreover, fish have long been used as sentinels for biomonitoring of aquatic environmental pollutants and are good indicators of mercury pollution [19, 20].
To our knowledge, there is no research report to elucidate the in vivo hepatotoxicity mechanism induced by mercury as well as its induced syndromes. Recent work by Gonzalez et al had indicated that the liver of the zebrafish is one of the main organs of mercury accumulation during the first 7 days of exposure to methylmercury . In this study, we aimed to investigate the in vivo mechanism of mercury-induced hepatotoxicity using the zebrafish as a model and compared it with other published data. To date, most mechanistic studies of mercury-induced toxicity in human/mammals have been performed either in vitro or in vivo using non-liver tissue. However, mechanistic insight obtained from the in vitro systems can be quite different from those obtained from the in vivo environment because of the lack of coordination of biological processes in whole organism. As the liver is the main organ that performs detoxification processes as well as regulation of metabolic pathways, it is important to understand toxicity effects and malfunctions that may implicate metabolic disorders induced by some toxicants in liver in vivo.
Here, using the DNA microarray technology, we captured the kinetics of hepatic transcriptome changes in zebrafish upon mercury exposure in order to generate a global transcriptome view of mercury-induced hepatotoxicity. Using gene set enrichment analysis (GSEA) and knowledge-based data mining, we compared our findings with those reported in human liver cell line and/or in mammalian model. We performed phenotype anchoring by linking histological analyses to the transcriptome data and subsequently performed gene expression validation using real-time PCR on a separate set of samples from seperate mercury exposure experiments. Taken together, we proposed a plausible in vivo mechanistic model of mercury-induced hepatotoxicity that could lead to syndromes related to mitochondrial dysfunction, endocrine disruption and metabolic disorders.
Results and discussion
Characterization of mercury-induced hepatotoxicity responses in zebrafish from transcriptome analysis
Preliminary acute toxicity test was first performed to determine the appropriate concentrations of mercuric chloride (HgCl2) for DNA microarray experiments. Fishes were treated with different concentrations of HgCl2 ranging from 100 to 600 μg/L for 7 days. Acute toxicity and motility test for HgCl2 exposure treatment revealed that 600 μg/L HgCl2 caused all the fish to die within 2 days, while treatment at concentrations of 200 μg/L produced 90% of survival rate up to 7 days (See Additional file 1). All surviving fish treated with 200 μg/L displayed normal swimming behavior without obvious adverse effect. Therefore, 200 μg/L of HgCl2 was chosen for the subsequent DNA microarray experiments as the concentration was sufficient to induce acute toxicity with minimal mortality.
Bioaccumulation of HgCl2 was not determined in the present study. In a similar study recently published , 0.15 mg/L HgCl2 was used for exposure in another tropical freshwater fish, Brycon amazonicus, and 10.46 mg Hg per kg liver from treated fish was detected in fish liver at 96 h and thus a bioaccumulation factor of about 100 in the liver was estimated.
Based on the transcriptome data, genes that were significantly regulated by mercury were identified at different time points (8 h, 24 h, 48 h and 96 h) with p-value < 0.05 and mapped to human homologous genes (See Additional file 2). Compared to 8 h time point, the number of genes significantly deregulated (both up- and down-regulated) at 24 h increased dramatically and a further increase was observed at 96 h time point. The increasing numbers of differentially expressed genes may serves as a rough indication that more biological processes were deregulated with the increase of duration of mercury exposure up to 96 h.
Gene Set Enrichment Analysis (GSEA) of liver transcriptome of HgCl2-treated zebrafish.
Electron transport chain
Mitochondrial fatty acid beta-oxidation
Nuclear receptor signaling pathway
Extrinsic apoptotic pathway
Mitochondrial fatty acid beta-oxidation
Complement activation (classical)
Mitochondrial fatty acid beta-oxidation
Electron transport chain
Nuclear receptor signaling pathway
Non-substrate GSK3 interacting proteins
Complement activation (classical)
DNA damage signaling pathway
Wnt signaling pathway
Fatty acid synthesis
Intrinsic apoptotic pathway
Hypoxia-induced factor pathway
Electron transport chain
The transcriptome analysis using GSEA revealed that genes associated with electron transport chain, mitochondrial fatty acid beta-oxidation, nuclear receptor signaling and apoptotic pathway were down-regulated as early as 8 h of HgCl2 exposure. This was followed by up-regulation of genes associated with the complement system and proteasome pathway, suggesting increased activity of these processes at 24 h. Major processes that were affected throughout the latter three time points of HgCl2 treatments (24 h, 48 h and 96 h) were proteasome pathway, electron transport chain, and mitochondrial fatty acid beta-oxidation (Table 1). Interestingly, electron transport chain that was initially down-regulated at 8 h were subsequently up-regulated at 48 h and 96 h, implying an increased need of intracellular energy production presumably to compensate the early down-regulation and to meet the energy requirement of the adaptive response to HgCl2 exposure. At 96 h of HgCl2 exposure, more biological pathways were significantly affected such as up-regulation of actin pathway, fatty acid synthesis, DNA damage signaling pathway, classical complement activation, intrinsic apoptotic pathway, cell motility, hypoxia-induced factor pathway, and gluconeogenesis, as well as down-regulation of cell cycle, non-substrate Gsk3 interacting proteins, and Wnt signaling pathway.
Taken together, the transcriptome analysis revealed that early targets of mercury-induced hepatotoxicity involved mitochondrial processes (electron transport chain and fatty acid beta-oxidation), apoptotic pathway and nuclear receptor signaling. The following up-regulation of complement system and proteasome pathway suggest activation of acute phase response and increased degradation of damaged proteins within 24 h of HgCl2 exposure. By 96 h, several critical processes associated with cellular maintenance, stress, survival and metabolism were clearly deregulated in the liver. Earlier work by Gonzalez et al using the zebrafish to elucidate effects of methylmercury on gene expression in liver, skeletal muscle, and brain indicated genes involved in mitochondrial metabolism, detoxification process, DNA repair and apoptosis were generally deregulated. Hence, there are similar observed toxicity effects between inorganic mercury and organic mercury such as methylmercury. Thus, it is important to further understand the mechanism of their in vivo toxicity in future works.
Comparative analysis revealed some common mercury-induced processes in both zebrafish liver and human liver cell line
Gene Set Enrichment Analysis (GSEA) of transcriptome of mercury-treated human HepG2 cells.
DNA damage signaling*
Electron transport chain*
Nuclear receptor signaling pathway*
Mitochondrial fatty acid beta-oxidation*
However, unlike zebrafish, HepG2 cells showed up-regulation of cell cycle, pyrimidine metabolism and DNA repair owing to the proliferative state of HepG2 cells (Table 2). This may indicate differential toxicity responses between proliferative cell cultures and highly differentiated in vivo models. More importantly, being an in vivo model, the zebrafish liver data has an advantage to capture deregulated biological processes that involves interactions of various cell types and whole-organism physiological metabolic processes such as fatty acid synthesis, gluconeogenesis, and other processes involving non-substrate Gsk3 interacting proteins and complement activation. Thus, using an in vivo model such as zebrafish may help researchers to characterize pathways that are perturbed by toxicants under in vivo environments.
HgCl2 affected liver cell morphology and cell adhesion
In order to further elucidate the mechanism of mercury-induced hepatotoxicity, we performed phenotypic anchoring [28, 29] by linking histological evidence to the mercury-induced transcriptome changes. This was further validated by real-time PCR for targeted genes associated with mercury-induced phenotype/pathology and transcriptome changes in the liver. In these experiments, we treated a new batch of zebrafish with different concentrations of HgCl2, from environmental-relevant (10 and 50 μg/L) to sub-lethal concentration (200 μg/L). Real-time PCR validation of targeted genes were performed on liver samples from zebrafish treated with 10, 50 and 200 μg/L for 24 h while phenotype anchoring were performed using liver from zebrafish treated with 50 and/or 200 μg/L for 96 h.
HgCl2 decreased hepatocyte cell number via induction of apoptosis
HgCl2 induced fatty liver condition
HgCl2 induced glycogen accumulation in liver
Comparative analysis revealed similarity of mercury- and arsenic-induced processes in zebrafish liver
We have previously reported that arsenic exposure in adult zebrafish caused decreased number of hepatic parenchyma cells as well as altered cell adhesion and cytoskeleton . Due to similarity of the toxicity effects between these two compounds, we compared transcriptome profiles between mercury and arsenic in the zebrafish liver using GSEA. Similar statistical criteria were used to select enriched pathways induced by arsenic (See Additional file 3). Arsenic induced proteasome pathway and intrinsic apoptotic pathway as early as 8 h of arsenic exposure. Unlike mercury, extrinsic apoptotic pathway is not enriched. Also, DNA damage signaling pathway was significantly up-regulated at 96 h of arsenic treatment although positive NES values were observed at an earlier time point. Like the case of mercury, Wnt signaling was also significantly down-regulated in Asenic-treated livers. However, unlike mercury, gluconeogenesis and fatty acid synthesis were down-regulated by arsenic (See Additional file 3) implying that there are different toxicity modes between these two compounds.
Validation of targeted genes as potential biomarkers of mercury-induced hepatotoxicity
Mechanistic insights into in vivo mercury-induced hepatotoxicity
Depletion of thiol-containing antioxidants will cause direct (red fonts in Figure 7) and indirect (blue fonts in Figure 7) oxidative stress-mediated liver injury. Direct oxidative stress induces DNA damage, mitochondrial dysfunction such as increased oxidative phosphorylation or electron transport chain activity and reduced activity of fatty acid β-oxidation. These will subsequently cause induction of mitochondria-mediated intrinsic apoptotic pathway to release Caspase-8 to trigger apoptosis and cell death leading to liver injury. Depletion of thiol-antioxidants results in the overall intracellular redox balance to shift towards increased oxidative stress, especially in the mitochondria where active oxidative phosphorylation occurs. This is evident from the transcriptome analysis indicating the deregulation of electron transport chain and mitochondrial fatty acid beta-oxidation in at least 3 time-points (8 h, 24 h and 48 h; Table 1), suggesting mitochondrial dysfunction occurring as early as 8 h of mercury exposure. Moreover a shift from down-regulation of extrinsic apoptosis at 8 h time-point to up-regulation of mitochondrial-mediated intrinsic apoptosis signaling at 96 h time-point (Table 1) strongly suggests increased mitochondrial dysfunction and intracellular stress [38, 39]. Increased apoptosis was confirmed through real-time PCR of up-regulated dffa (encoding a substrate of Caspase 3 that triggers DNA fragmentation during apoptosis) as showed in Figure 6 and immunohistochemical-staining of apoptotic DNA fragments in Figure 3. Some of the observed histological phenotypes in the liver parenchyma, such as irregular cell-shape, loose contact with neighboring cells and significantly reduced nucleated-cell count, may be attributed to morphological features of apoptosis  as well as to mercury interaction with cytoskeletal proteins . E-cadherin and Keratin-8 staining further revealed the extent of the disrupted cell adhesion and cytoskeletal organization corroborating with the deregulated actin pathway and cell motility process captured by the transcriptome data. Unmitigated oxidative stress would cause cellular damage involving DNA and proteins and lead to tissue injury as well as cell lysis/death, as suggested by the up-regulation of DNA damage signaling, proteasome pathway, several acute phase protein genes (fga, agt) including complement activation. This further explains the down-regulation of cell cycle indicated in the transcriptome analysis which is necessary for DNA repair and assessment of cellular damage, and if necessary, to induce cell death instead of cell proliferation.
In addition to direct mercury-induced oxidative cellular injury, our study has also captured indirect mercury-induced oxidative stress-mediated gene expression/regulation (blue fonts in Figure 7). Depletion of cellular thiols and intracellular redox balance can alter activity in many transcription regulators and nuclear receptors [41, 42] as indicated by the early deregulation nuclear receptor signaling pathway (Table 1) that include genes such as esr1 (estrogen receptor 1) and nr2f1l, further confirming the endocrine-disrupting effect of mercury . Our findings suggest that the liver is a mercury-targeted organ for endocrine disruption. In addition to hormonal receptors, the up-regulation of transcription regulators such as junb, β-c/ebp, δ-c/ebp, and keap1 suggested modulation of protective/defense and survival mechanisms involving genes associated with antioxidant response element (ARE) (keap1, junb, junc, gstp1, gclm, sqstm1) and acute phase response (β-c/ebp, δ-c/ebp, cfb, c3, c9, fga, fgb, hpx) (See Additional file 4). Up-regulated glucocorticoid receptor signaling together with β-c/ebp could lead to gluconeogenesis and glycogen synthesis as indicated in the transcriptome analysis and PAS staining of glycogen as shown in Figure 5. Moreover, up-regulation β-c/ebp and δ-c/ebp showed in Figure 4 together with the disruption of mitochondrial fatty acid beta-oxidation could lead to fat synthesis and accumulation. All these incidents can eventually lead to metabolic disorders that may cause further liver injury.
Apart from transcription regulators, it was also observed that expression of genes encoding several kinases (cdc2, gsk3a, mapkapk2, ptk2 and sgk) were also affected suggesting that signal transduction via phosphorylation was altered by mercury-induce oxidative stress. Modulation of kinases involved in MAPK (Mitogen-Activated Protein Kinase) pathways could result in up-regulation of hypoxia signaling . Furthermore, deregulated kinases such as Gsk3 will trigger hypoxia-inducible pathways and intrinsic apoptotic pathway. Gsk3 had been shown to activate intrinsic but not extrinsic apoptotic pathway  and deregulate Wnt signaling pathway [46, 47] as indicated in the transcriptome analysis. Likewise, Wnt signaling can act as negative regulator for adipogenesis [31, 32] as its down-regulation may also contribute to adipogenesis observed in the study.
Finally, it should be pointed out that the analysis is based on the nominal concentration of mercury used in the exposure experiment and does not reflect the amount of mercury bioaccumulated in the tissue. Future studies should include bioaccumulation of mercury in target tissues in order to relate the observed effects to its actual concentration in the tissues. This is important for assessing and comparing the actual degree of mercury exposure in the tissue with its corresponding toxicological responses from different studies.
Our study revealed in vivo mechanistic insights of mercury-induced hepatotoxicity that may imply complications associated with metabolic disorders. To our knowledge, this is the first report of an in vivo mechanistic model for mercury-induced toxicity in the liver. Comparative analyses from arsenic-induced hepatotoxicity in zebrafish and human HepG2 cell line treated with HgCl2 showed some similar effects with mercury-induced hepatotoxicity in zebrafish. In addition, transcriptome analysis coupled with phenotypic anchoring and targeted gene validation had revealed that mercury can cause liver injury via broad processes associated to oxidative stress and cell death as well as deregulation of kinases such as Gsk3 in regulating the activities of Wnt signaling pathway, gluconeogenesis and adipogenesis that may eventually lead to syndromes such as mitochondrial dysfunction, endocrine disruption and metabolic disorders. Moreover, targeted gene validation also identified genes whose expressions are sensitive to mercury exposure suggesting their potential use for detecting or diagnosing mercury-induced toxicity. Thus, the understanding of the in vivo mechanism of hepatotoxicity using zebrafish could facilitate researches in toxicology and toxicogenomics, in particular in vivo heavy metal toxicity. These findings provide strong evidence of the feasibility of using zebrafish to model mercury-induced hepatotoxicity in human for risk assessment of mercury.
Adult zebrafish (6 months - 1 year old) were obtained from a local fish farm. The fishes were allowed to acclimatize in aquaria for several days before transferred into smaller containers for mercury exposure. For all three types of experiments (i.e. microarray, histology and PCR biomarker validation), zebrafish were exposed to mercury at different concentrations for 96 h at a density of one fish/200 ml at 27 ± 2°C. De-chlorinated water and chemical were renewed daily. All experimental with animals were conducted following the guidelines of Institutional Animal Care and Use Committee (IACUC) and the procedures were approved by IACUC.
HgCl2 or mercury (II) chloride [99.5% A.C.S reagent; Sigma, Aldrich] was chosen in consideration of its solubility in water that greatly facilitates daily preparation of chemical stock. Concentrations for treatment were chosen based on pre-LC50 and LC50 results. Concentration used for treatment is from 0.01 - 0.2 mg/L. Prior to obtaining the transcriptomic profiles of zebrafish, the appropriate dosage of HgCl2 needed was determined from acute toxicity tests. Survival analyses for zebrafish treated with various concentrations of HgCl2 were performed. The survival rate of sub-lethal 200 μg/L of HgCl2 was 90% at day 7 was chosen for subsequent microarray and toxicological studies. Treatment for microarray analyses was carried out in 4 days, using one concentration (200 μg/L) of HgCl2 with liver RNA sampling at 8, 24, 48, and 96 hours, respectively. Treatment was conducted using triplicate groups of four zebrafishes.
DNA Microarray Experiments
Total RNAs from triplicates (each replicate consisted of liver pooled from four fishes) at the respective sampling time-points were isolated by Trizol reagent from Gibco-BRL. Reference RNA was obtained by pooling equal amount of total RNA extracted from male and female wildtype fish. Reference RNA was co-hybridized with RNA samples either from treated or control fish on a glass array spotted with 16,416 zebrafish oligo probes (16,244 LEADS clusters and 172 copies of β-actin controls) . Both reference and sample RNAs were reverse-transcribed and labeled differently, using fluorescent dyes Cy-3 and Cy-5, respectively. After hybridization at 42°C for 16 hours in hybridization chambers (Gene Machines), the microarray slides were washed in a series of washing solutions (2× SSC with 0.1% SDS; 1× SSC with 0.1% SDS; 0.2× SSC and 0.05× SSC; 30 seconds each), dried by low-speed centrifugation and scanned for fluorescence detection using the GenePix 4000B microarray scanner (Axon Instruments). Further detailed protocols for microarray experiment and data acquisition have been recently published by us [49, 50].
Microarray data processing and transcriptome analysis
The raw microarray data was normalized using Lowess method in the R package http://www.braju.com/R/. The data had been submitted to Gene Expression Omnibus (GEO) with accession number of GSE18861. One-way ANOVA was performed on the normalized data followed by Student's t-test to identified genes whose expressions are significantly altered by HgCl2. Gene set enrichment analysis (GSEA)  was performed to characterize pathways or processes that are affected by HgCl2. Fishes were retreated with HgCl2 and quantitative real-time PCR was used to validate genes whose expressions were significantly altered as well as some regulatory genes in relevant pathways affected by HgCl2.
Analysis of Transcriptomics Profile Using Gene Set Enrichment Analysis (GSEA)
Gene Set Enrichment Analyses (GSEA) as described in detail by  was used to determine biological pathways perturbed by HgCl2. The zebrafish genes were mapped to human homologous genes as previously described . The human homologs of zebrafish genes from the HgCl2-induced transcriptome profiles were ranked based on their p-values using Student t-test. The "GSEAPreranked" option of GSEA was used. The ranking metric used was log10 (1/P) where P is the p-value of a gene from microarray data. Up-regulated genes will carry positive values of log10 (1/P) whereas down-regulated genes will carry negative values of log10 (1/P). The genes were then ranked in descending order based on values of log10 (1/P). The ranking procedure was conducted for array at each time point of HgCl2 treatment. The ranked list of genes for each time-point treatment will then be compared to 639 pre-defined gene sets or signatures of canonical pathways. These pre-defined gene sets can be downloaded from the GSEA website.
An enrichment score (ES) that reflects the degree to which a pre-defined gene set is overrepresented at the top or bottom rank of the ranked whole transcriptome profile was calculated by walking down the ranked profile. The ES is the maximum deviation from zero encountered in the random walk corresponds to a weighted Kolmogorov-Smirnov-like statistic. The statistical significant (nominal p-value, NP) of the ES was estimated by using an empirical phenotype-based permutation test procedure. The phenotype labels were permuted and the ES of the gene set for the permuted data were recomputed that generate a null distribution for the ES. The empirical, nominal p-value of the observed ES was then calculated relative to this null distribution. The estimated significance level was adjusted with multiple hypothesis testing. The ES for each gene set was first normalized to the size of the set yielding a normalized enrichment score (NES) with the following relation:
NES = actual ES/mean (ESs against all permutations of the dataset)
The number of permutation used was 1000. The proportion of false positive was then controlled by calculating the false discovery rate (FDR) corresponding to each NES.
Pathways with false discovery rate (FDR) < 0.3 or nominal p-value (NP) < 0.1 were considered statistically significant. Positive and negative values of normalized enrichment scores (NES) indicated the activities of pathways that were up- and down-regulated, respectively.
Biomarker Validation with Real Time Quantitative PCR
The quantitative Real-Time PCR reaction was performed from synthesized FirstStrand cDNA using LightCycler® 480 SYBR Green I Master kit according to the manufacturer's protocol (Roche). Six biological replicates were performed for all real-time PCR experiments. Analysis of the transcript levels were performed by using relative quantification between the PCR signal of the target transcript in a treatment group and that of an untreated control group after normalization with the transcript level of beta-actin. The detail information of primers used is given in Additional file 5.
Histological Sectioning and Sample Preparation
We repeated the treatment of zebrafish with HgCl2 for histological analysis. Adult zebrafish were treated with different concentrations (50 μg/L, 200 μg/L) of HgCl2 [99.5% A.C.S reagent; Sigma-Aldrich] for 96 hours at a density of 1 fish/200 ml at 27 ± 2°C in a static condition. Control fish were kept in water under similar conditions, and a total of 4 fish were used for each group. Treated solution and water were changed daily. Fish were slit open ventrally from heart to the anus to expose the digestive organs. The whole gut (alimentary tract including the liver) were removed and fixed in Formalin solution 10%, Neutral Buffered (Sigma-Aldrich) for 2-3 days at room temperature. Fixed gut samples were washed several times with 70% ethanol, followed by dehydration in a grade series of ethanol (70% - 100%) before clearing in Histoclear and embedding in paraffin. The paraffin-embedded tissues were sectioned serially at 5 μm thickness. The sections were stained with hematoxylin and eosin (H&E) for qualitative and quantitative histological analysis.
Periodic acid-Schiff (PAS) staining
PAS staining is mainly used for staining structures containing a high proportion of carbohydrate macromolecules (glycogen, glycoprotein, proteoglycans) typically found in connective t/s, mucus and basal lamina). Staining for glycogen was performed using Alcian Blue PAS stain kit without diastase according to the manufacturer's protocol (BioGenex).
Oil red O staining
Staining for adipocytes (Lipid and fat globules) was performed using Cryostat Sectioning from fresh frozen liver samples. Staining procedure was performed using Oil Red O stain according to the manufacturer's protocol (Sigma-Aldrich).
DNA fragmentations associated with ultrastructural changes in cellular morphology in apoptosis were detected using ApopTag®Plus Fluorescein In Situ Apoptosis Detection Kit according to the manufacturer's protocol (Chemicon). The 3'-OH ends of double-stranded or single-stranded DNA were labeled with the digoxigenin-nucleotide and then allowed to bind an anti-digoxigenin antibody (Anti-DIG) that is coujugated to alkaline phosphatase. For The DNA fragmentations localized in apoptotic bodies were detected enzymatically using 5-Bromo-4-chloro-3-indolyl phosphate (BCIP)/Nitroblue tetrazolium (NBT) substrate.
Histopathological analysis was performed to obtain a tissue-level image of HgCl2 effects on zebrafish at 96 hours. This also serves to further understand the gene expression profiles resulted from the microarray analyses. Imaging of liver sections was carried out using a compound microscope, Axioskop 2, Zeiss® equipped with an imaging system and were analyzed with a computer-assisted image analyzer program (Axiovision, Zeiss). Liver sections from treated fish were compared with control fish, both quantitatively (i.e. cell size and cell density) and qualitatively (i.e. changes observed in liver tissue) from hematoxylin and eosin-stained section. Density of nucleated cells (no. of nucleated cells/5,100 μm2) and the size of nucleus were determined with an image analyzer in both the HgCl2-treated group and control group. Three images from each portion (anterior, middle and posterior portions) of liver sections (1,000× magnification) for each liver from four control and four HgCl2-treated fish (n = 4 liver samples) were used to determine the density of nucleated cells and to measure the size of the nuclei for each image. A two-tailed t-test was performed to analyze the data for statistical significance (P < 0.05).
Images of the PAS stained sections (1,000× magnifications) were captured using Axioskop 2 for each liver from four control and four HgCl2-treated fishes. Comparison of the glycogen content in control and treated liver samples were performed. Images of the Oil-red O stained sections (1,000× magnifications) were also captured using Axioskop 2. Comparison of the adipocytes (lipid globules) stained with Oil-red O from each image of control and treated liver sections were performed.
This work is supported by Environment and Water Industry (EWI) Development Council of Singapore and BioMedical Research Council (BMRC) of Singapore.
- Li P, Feng XB, Qiu GL, Shang LH, Li ZG: Mercury pollution in Asia: A review of the contaminated sites. J Hazard Mater. 2009, 168: 591-601. 10.1016/j.jhazmat.2009.03.031.PubMedView ArticleGoogle Scholar
- Pacyna EG, Pacyna JM, Steenhuisen F, Wilson S: Global anthropogenic mercury emission inventory for 2000. Atmos Environ. 2006, 40: 4048-4063. 10.1016/j.atmosenv.2006.03.041.View ArticleGoogle Scholar
- Kudo A, Fujikawa Y, Miyahara S, Zheng J, Takigami H, Sugahara M, Muramatsu T: Lessons from minamata mercury pollution. Japan-After a continuous 22 years of observation. Water Sci Technol. 1998, 38: 187-193.View ArticleGoogle Scholar
- Pfeiffer WC, Lacerda LD: Mercury inputs in the Amazon region. Environ Technol Lett. 1988, 9: 325-330. 10.1080/09593338809384573.View ArticleGoogle Scholar
- Bakir F, Damluji SF, Amin-Zaki L, Murtadha M, Khalidi A, al-Rawi NY, Tikriti S, Dahahir HI, Clarkson TW, Smith JC, Doherty RA: Methylmercury poisoning in Iraq. Science. 1973, 181: 230-241. 10.1126/science.181.4096.230.PubMedView ArticleGoogle Scholar
- Trasande L, Landrigan PJ, Schechter C: Public health and economic consequences of methyl mercury toxicity to the developing brain. Environ Health Persp. 2005, 113: 590-596.View ArticleGoogle Scholar
- Hylande LD, Goodsite ME: Environmental costs of mercury pollution. SciTotal Environ. 2006, 368: 352-370. 10.1016/j.scitotenv.2005.11.029.Google Scholar
- Agency for Toxic Substances and Disease Registry (ATSDR): Toxicological Profile for Mercury. 1999, Atlanta, GA: U.S. Department of Health and Human Services, Public Health ServiceGoogle Scholar
- EPA US: Mercury Report to Congress Office of Air Quality and Standards. 1997, U.S. Environmental Protection Agency, Washington, DC, EPA 452/R-97-0003Google Scholar
- EPA US: Water Quality Criterion for the Protection of Human Health: Methyl Mercury. 2001, U.S. Environmental Protection Agency, Washington, DC, EPA 0823-R-01-001Google Scholar
- Clarkson TW: The Three Modern Faces of Mercury. Environ Health Perspect. 2002, 110 (suppl 1): 11-23.PubMed CentralPubMedView ArticleGoogle Scholar
- Dutczak WJ, Ballatori N: Transport of the glutathionemethylmercury complex across liver canalicular membranes on reduced glutathione carriers. J Biol Chem. 1994, 269: 9746-9751.PubMedGoogle Scholar
- Castoldi AF, Coccini T, Ceccatelli S, Manzo L: Neurotoxicity and molecular effects of methylmercury. Review Brain Res Bull. 2001, 55: 197-203. 10.1016/S0361-9230(01)00458-0.View ArticleGoogle Scholar
- Zalups RK: Molecular interactions with mercury in the kidney. Pharmacol Rev. 2000, 52: 113-143.PubMedGoogle Scholar
- Suter L, Babiss LE, Wheeldon EB: Toxicogenomics in predictive toxicology in drug development. Chem Biol. 2004, 11: 161-171.PubMedView ArticleGoogle Scholar
- Grunwald DJ, Eisen JS: Headwaters of the zebrafish-emergence of a new model vertebrate. Nat Rev Genet. 2002, 3: 717-723. 10.1038/nrg892.PubMedView ArticleGoogle Scholar
- Shin JT, Fishman MC: From zebrafish to human: modular medical models. Annu Rev Genomics Human Genet. 2002, 3: 311-340. 10.1146/annurev.genom.3.031402.131506.View ArticleGoogle Scholar
- Spitsbergen JM, Kent ML: The state of the art of the zebrafish model for toxicology and toxicologic pathology research-advantages and current limitation. Toxicol Pathol. 2003, 31: 62-87.PubMed CentralPubMedGoogle Scholar
- Carvan MJ, Dalton TP, Stuart GW, Nebert DW: Transgenic zebrafish as sentinels for aquatic pollution. Acad Sci. 2000, 919: 133-147.View ArticleGoogle Scholar
- Carvan MJ, Solis WA, Gedamu L, Nebert DW: Activation of transcription factors in zebrafish cell cultures by environmental pollutants. Arch Biochem Biophys. 2000, 376: 320-327. 10.1006/abbi.2000.1727.PubMedView ArticleGoogle Scholar
- Hill AJ, Teraoka H, Heideman W, Peterson RE: Zebrafish as a model vertebrate for investigating chemical toxicity. Toxicol Sci. 2005, 86: 6-19. 10.1093/toxsci/kfi110.PubMedView ArticleGoogle Scholar
- Lam SH, Wu YL, Vega VB, Miller LD, Spitsbergen J, Tong Y, Zhan H, Govindarajan KR, Lee S, Mathavan S, Murthy KR, Buhler DR, Liu ET, Gong Z: Conservation of gene expression signatures between zebrafish and human liver tumors and tumor progression. Nat Biotechnol. 2006, 24: 73-75. 10.1038/nbt1169.PubMedView ArticleGoogle Scholar
- Lam SH, Winata CL, Tong Y, Korzh S, Lim WS, Korzh V, Spitsbergen J, Mathavan S, Miller LD, Liu ET, Gong Z: Transcriptome kinetics of arsenic-induced adaptive response in zebrafish liver. Physiol Genomics. 2006, 27: 351-361. 10.1152/physiolgenomics.00201.2005.PubMedView ArticleGoogle Scholar
- 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: e1000121-10.1371/journal.pgen.1000121.PubMed CentralPubMedView ArticleGoogle Scholar
- Gonzalez P, Dominique Y, Massabuau JC, Boudou A, Bourdineaud JP: Comparative effects of dietary methylmercury on gene expression in liver, skeletal muscle, and brain of the zebrafish (Danio rerio). Environ Sci Technol. 2005, 39: 3972-3980. 10.1021/es0483490.PubMedView ArticleGoogle Scholar
- Monteiro DA, Rantin FT, Kalinin AL: Inorganic mercury exposure: toxicological effects, oxidative stress biomarkers and bioaccumulation in the tropical freshwater fish matrinxã, Brycon amazonicus (Spix and Agassiz, 1829). Ecotoxicology. 2010, 19: 105-123. 10.1007/s10646-009-0395-1.PubMedView ArticleGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMed CentralPubMedView ArticleGoogle Scholar
- Paules R: Phenotypic anchoring: linking cause and effect. Environ Health Perspect. 2003, 111: A338-A339.PubMed CentralPubMedView ArticleGoogle Scholar
- Moggs JG, Tinwell H, Spurway T, Chang HS, Pate I, Lim FL, Moore DJ, Soames A, Stuckey R, Currie R, Zhu T, Kimber I, Ashby J, Orphanides G: Phenotypic anchoring of gene expression changes during estrogen-induced uterine growth. Environ Health Perspect. 2004, 112: 1589-1606.PubMed CentralPubMedView ArticleGoogle Scholar
- Araragi S, Kondoh M, Kawase M, Saito S, Higashimoto M, Sato M: Mercuric chloride induces apoptosis via a mitochondrial-dependent pathway in human leukemia cells. Toxicology. 2003, 184: 1-9. 10.1016/S0300-483X(02)00443-2.PubMedView ArticleGoogle Scholar
- Rangwala SM, Lazar MA: Transcriptional control of adipogenesis. Annu Rev Nutr. 2000, 20: 535-559. 10.1146/annurev.nutr.20.1.535.PubMedView ArticleGoogle Scholar
- Rosen ED, Walkey CJ, Puigserver P, Spiegelman BM: Transcriptional regulation of adipogenesis. Genes Dev. 2000, 14: 1293-1307.PubMedGoogle Scholar
- Rao MV, Sharma PS: Protective effect of vitamin E against mercuric chloride reproductive toxicity in male mice. Reprod Toxicol. 2001, 15: 705-712. 10.1016/S0890-6238(01)00183-6.PubMedView ArticleGoogle Scholar
- Anh LL, Loi VD, Mon PG, Phu NL, Ngoc DM, Tuy P: Determination of lead and mercury in clinical samples during work and accidental exposure in Vietnam. Anal Sci. 2001, 17: a37-a39. 10.2116/analsci.17.37.View ArticleGoogle Scholar
- Butler WJ, Houseman J, Seddon L, McMullen E, Tofflemire K, Mills C, Corriveau A, Weber JP, LeBlanc A, Walker M, Donaldson SG, Van Oostdam J: Maternal and umbilical cord blood levels of mercury, lead, cadmium, and essential trace elements in Arctic Canada. Environ Res. 2006, 100: 295-318. 10.1016/j.envres.2005.05.006.View ArticleGoogle Scholar
- Counter SA, Buchanan LH: Mercury exposure in children: a review. Toxicol Appl Pharmacol. 2004, 198: 209-230. 10.1016/j.taap.2003.11.032.PubMedView ArticleGoogle Scholar
- Yoon S, Han SS, Rana SV: Molecular markers of heavy metal toxicity--a new paradigm for health risk assessment. J Environ Biol. 2008, 29: 1-14.PubMedGoogle Scholar
- Kroemer G, Reed JC: Mitochondrial control of cell death. Nat Med. 2000, 6: 513-519. 10.1038/74994.PubMedView ArticleGoogle Scholar
- Lockshin RA, Zakeri Z: Programmed cell death and apoptosis: origins of the theory. Nat Rev Mol Cell Biol. 2001, 2: 545-550. 10.1038/35080097.PubMedView ArticleGoogle Scholar
- Saraste A, Pulkki K: Morphologic and biochemical hallmarks of apoptosis. Cardiovasc Res. 2000, 45: 528-537. 10.1016/S0008-6363(99)00384-3.PubMedView ArticleGoogle Scholar
- Arrigo AP: Gene expression and the thiol redox state.Free Radic. Boil Med. 1999, 27: 936-944.Google Scholar
- Sen CK: Cellular thiols and redox-regulated signal transduction. Curr Top Cell. 2000, 36: 1-30. full_text.View ArticleGoogle Scholar
- Tan SW, Meiller JC, Mahaffey KR: The endocrine effects of mercury in humans and wildlife. Crit Rev Toxicol. 2009, 39: 228-269. 10.1080/10408440802233259.PubMedView ArticleGoogle Scholar
- Galanis A, Karapetsas A, Sandaltzopoulos R: Metal-induced carcinogenesis, oxidative stress and hypoxia signaling. Mutat Res. 2009, 674: 31-35.PubMedView ArticleGoogle Scholar
- Beurel E, Jope RS: The paradoxical pro- and anti-apoptotic actions of GSK3 in the intrinsic and extrinsic apoptosis signaling pathways. Prog Neurobiol. 2006, 79: 173-189. 10.1016/j.pneurobio.2006.07.006.PubMed CentralPubMedView ArticleGoogle Scholar
- Cohen P, Frame S: The renaissance of GSK3. Nat Rev Mol Cell Biol. 2001, 2: 769-776. 10.1038/35096075.PubMedView ArticleGoogle Scholar
- Jope RS, Johnson GV: The glamour and gloom of glycogen synthase kinase-3. Trends Biochem Sci. 2004, 29: 95-102. 10.1016/j.tibs.2003.12.004.PubMedView ArticleGoogle Scholar
- Mathavan S, Lee SG, Mak A, Miller LD, Murthy KR, Govindarajan KR, Tong Y, Wu YL, Lam SH, Yang H, Ruan Y, Korzh V, Gong Z, Liu ET, Lufkin T: Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005, 1: 260-276. 10.1371/journal.pgen.0010029.PubMedView ArticleGoogle Scholar
- Lam SH, Mathavan S, Gong Z: Zebrafish spotted-microarray for genome-wide expression profiling experiments. Part I: array printing and hybridization. Methods Mol Biol. 2009, 546: 175-195. full_text.PubMedView ArticleGoogle Scholar
- Lam SH, Krishna Murthy Karuturi R, Gong Z: Zebrafish spotted-microarray for genome-wide expression profiling experiments: data acquisition and analysis. Methods Mol Biol. 2009, 546: 197-226. full_text.PubMedView ArticleGoogle Scholar