Quantitative analyses of the hepatic proteome of methylmercury-exposed Atlantic cod (Gadus morhua) suggest oxidative stress-mediated effects on cellular energy metabolism
© The Author(s). 2016
Received: 20 February 2016
Accepted: 30 June 2016
Published: 5 August 2016
Methylmecury (MeHg) is a widely distributed environmental pollutant with considerable risk to both human health and wildlife. To gain better insight into the underlying mechanisms of MeHg-mediated toxicity, we have used label-free quantitative mass spectrometry to analyze the liver proteome of Atlantic cod (Gadus morhua) exposed in vivo to MeHg (0, 0.5, 2 mg/kg body weight) for 2 weeks.
Out of a toltal of 1143 proteins quantified, 125 proteins were differentially regulated between MeHg-treated samples and controls. Using various bioinformatics tools, we performed gene ontology, pathway and network enrichment analysis, which indicated that proteins and pathways mainly related to energy metabolism, antioxidant defense, cytoskeleton remodeling, and protein synthesis were regulated in the hepatic proteome after MeHg exposure. Comparison with previous gene expression data strengthened these results, and further supported that MeHg predominantly affects many energy metabolism pathways, presumably through its strong induction of oxidative stress. Some enzymes known to have functionally important oxidation-sensitive cysteine residues in other animals are among the differentially regulated proteins, suggesting their modulations by MeHg-induced oxidative stress. Integrated analysis of the proteomics dataset combined with previous gene expression dataset showed a more pronounced effect of MeHg on amino acid, glucose and fatty acid metabolic pathways, and suggested possible interactions of the cellular energy metabolism and antioxidant defense pathways.
MeHg disrupts mainly redox homeostasis and energy generating metabolic pathways in cod liver. The energy pathways appear to be modulated through MeHg-induced oxidative stress, possibly mediated by oxidation sensitive enzymes.
Mercury is released into the environment from both natural and anthropogenic sources, and occurs in different chemical forms in the environment, such as elemental mercury, inorganic mercury, and also as more toxic organic mercury compounds. The different forms of mercury can be interconverted in a global mercury cycle between the aquatic systems and the atmosphere . Methylmercury (MeHg), which is an organic compound, is formed from inorganic mercury in the aquatic environment by anaerobic bacteria . In aquatic environments, MeHg tends to bioconcentrate and bioaccumulate in organisms, and it has an extremely high potential for biomagnification up the food-chain . MeHg can therefore affect the health of both aquatic animals and humans, as humans can be exposed to MeHg through consumption of fish and other seafood. Because of safety concerns, food safety authorities often issue guidelines for the public to avoid or limit the consumption of some fish species from contaminated environments or with high muscle mercury levels .
MeHg is highly toxic to animal tissues and organs, especially the brain [1, 4, 5]. Its neurotoxic properties are partly due to its easy transport throughout the body, including the blood–brain and blood-placental barriers, as a complex with cysteine . The MeHg–L-cysteine complex has a similar structure as methionine and can be transported across cellular membranes by L-type amino acid transporters . Its neurotoxic effects have been well documented from tragic human poisoning incidents in Japan (Minamata Bay) and Iraq (rural areas), and also from controlled laboratory experiments [5, 7, 8]. Despite the focus of many studies on the brain where most serious effects are elicited, MeHg affects other animal organs and tissues as well [1, 4, 5]. For example, in fish, effects on other tissues such as the liver and kidney have been observed [9, 10]. Further, a recent epidemiological study has reported associations between blood mercury levels and changes in liver enzymes . Among the main cellular effects caused by MeHg are oxidative damage and mitochondrial dysfunction, disruption of Ca2+ homeostasis, apoptosis, inhibition of microtubule assembly, and inhibition of protein synthesis [5, 12–16].
The underlying molecular mechanism of MeHg toxicity seems to be related to its ability to bind thiol-groups present in biological molecules, including proteins. This may lead to disruption of enzymatic functions, and particularly inactivation of antioxidant systems, thereby promoting generation of reactive oxygen species (ROS) [12, 13]. Reactive thiol-groups of antioxidants such as glutathione and proteins in the antioxidant defense systems are known targets of MeHg-induced oxidation . Indeed, MeHg has been shown to inhibit components in the two major antioxidant defense pathways, the thioredoxin and glutaredoxin enzyme systems [17–19]. Compromised antioxidant defense can potentially result in severe oxidative damage to sensitive proteins, and thus affecting related pathways.
Analyses of gene expression and proteomic responses are increasingly useful in studying the effects and toxicity mechanisms of environmental contaminants . Quantitative proteomics approaches for mapping molecular targets of toxicants may lead to a better understanding of the molecular mechanisms of toxicity and discovery of new biomarkers . Recent studies using proteomics have reported that MeHg causes alterations of proteins and pathways related to antioxidant defense mechanisms, energy metabolism, and the cytoskeleton in various organs and cells in both mammals [22, 23] and fish [24–27]. These studies suggest that MeHg has similar molecular mechanisms of toxicity in various tissues of different animal species. Although the liver is the most important organ in metabolism and detoxification of environmental chemicals, it has not been well investigated as a target organ of MeHg toxicity, particularly by large-scale proteomics approaches.
The Atlantic cod (Gadus morhua) is a commercially important species, but increasingly susceptible to marine pollution from sources such as shipping activities, urban effluent discharges, and discharges from offshore petroleum industries. Importantly, Atlantic cod is increasingly becoming a useful model organism to investigate effects of contaminants, and as a sentinel species used in environmental monitoring of marine pollution [28–34]. It is one of few fish species with a sequenced and annotated genome , thus opening the possibility of using omics methods in molecular toxicological studies [25, 36, 37].
The aim of this study was to investigate the molecular targets and mechanisms of toxicity of MeHg on the liver proteome of Atlantic cod. In our previous work, global hepatic gene expression analysis with microarrays was used to study the effect of MeHg on gene expression using the same samples . As part of an integrated approach using both transcriptomics and proteomics, we have in this study performed label-free mass spectrometry analysis of the liver proteomes of cod exposed in vivo to 0, 0.5 and 2 mg/kg body weight MeHg. In the present study, 125 differentially regulated proteins were identified, and functional enrichment analysis showed that oxidative-stress responses and energy pathways were among the most affected by the MeHg treatment. These findings are in line with our previous gene expression data, and further suggest that MeHg may cause vulnerabilities of cellular energy pathways to oxidative stress.
Results and discussion
Fish exposure, survival and condition factor
To examine dose–response effects of MeHg on Atlantic cod, an in vivo exposure experiment with vehicle control and three different doses of MeHg (0.5, 2 and 8 mg/kg BW) were performed (n = 10/group). All fish in the highest (8 mg/kg BW) dose group, and 1 fish from each of the 0.5 and 2 mg dose groups, died by the end of the 14 days exposure time . The lowest dose of MeHg (0.5 mg/kg BW) was chosen based on approximates of mercury levels in fish reported in some field studies. For example, mean mercury concentrations of up to 0.46 mg/kg wet-weight were detected in fish caught in the Northwest Atlantic Ocean, and up to 2.07 mg/kg wet-weight in fish from the Mediterranean Sea have been reported . The maximum limits of Hg allowed in fish and fish products for human consumption is 0.5 or 1 mg/kg wet-weight, depending on the species . Previous chemical analysis of tissue concentrations of Hg in the fish studied here, showed accumulation in the liver with average levels of 0.8 and 2.8 mg Hg/kg wet weight for the 0.5 and 2 mg/kg BW MeHg doses, respectively . Higher Hg levels were detected in the muscle with 0.9 and 4.7 mg/kg wet weight for 0.5 and 2 mg/kg BW doses, respectively , consistent with a known tendency of MeHg to accumulate in muscle tissue . No significant changes were observed in behavioural parameters such as overall activity and balance during the 2 weeks of MeHg exposure period . There were small significant decreases in body weight at the end of the experiment in all groups, including the control (Additional file 1: Figure S1A). Condition factor decreased only slightly (not significant) in the treated groups (Additional file 1: Figure S1B) and there were no significant differences in liver-somatic indices between control and treated groups (Additional file 1: Figure S1C). Although the slight decrease in condition factor suggests some effects of the treatment, the observed decrease in body weight does not seem to be the main effect of MeHg treatment as it was also observed in the control group (Additional file 1: Figure S1A).
Differentially regulated proteins
Among the differentially regulated proteins, many were previously shown to have redox-sensitive Cys-residues in other studies. These include the thioredoxin enzyme system proteins , as well as proteins that are part of the fatty acid metabolism pathways, such as mitochondrial acyl-CoA dehydrogenases . In addition, chaperons, tubulin and hemoglobin subunits were among the differentially regulated proteins between MeHg-treated and control groups (Additional file 2: Table S1), which is also in line with previous two-dimensional gel electrophoresis analyses of the thiol-proteome obtained from the same liver samples .
Some of the proteins and enzymes in the cod liver that responded to MeHg exposure have previously been associated with different human diseases. For instance, one of the proteins that increased the most after MeHg exposure, T23O (Additional file 2: Table S1, Fig. 2a) was reported to have increased abundance in the brain of patients with schizophrenia and Alzheimer’s disease (AD) [45, 46]. T23O is a rate limiting-enzyme in the tryptophan degradation pathway, and regulates the levels of the amino acid tryptophan in the body. Apolipoprotein E (ApoE) variants and other cholesterol-related genes have been associated with the pathogenesis of AD . The protein levels of APOB, APOE, APOC1 and APOA4 decreased in the high dose MeHg-treated fish liver (Additional file 2: Table S1). It has been shown that genetic ablation of Apoa4 gene in a mouse model accelerated the development of AD . The differential regulation of the products of the genes previously associated with the AD appears to support the hypothesis that environmental chemicals such as MeHg may contribute to its etiology ).
The top significantly enriched GO BP, MF and MetaCore “Process Networks” annotationsa
A) GO processes
Small molecule metabolic process
Response to inorganic substance
Macromolecular complex remodeling
Plasma lipoprotein particle remodeling
Protein-lipid complex remodeling
Glycerolipid catabolic process
Single-organism metabolic process
Plasma lipoprotein particle organization
B) GO Molecular functions
O-Acyltransferase activator activity
Small molecule binding
Cholesterol transporter activity
Sterol transporter activity
Poly(A) RNA binding
C) Process Networks
Response to hypoxia and oxidative stress
Cell adhesion_Cell junctions
Selected reaction monitoring (SRM)
To verify the label-free quantitative proteomics data, six proteins (T23O, GLNA EPS8L2, APOA4, RAP1B, CATZ) among the top differentially regulated were selected and quantified with the more specific SRM method. At least one peptide originating from each of the six selected proteins was quantified using SRM. Peptides from five of the selected proteins (T23O, GLNA EPS8L2, APOA4, RAP1B) showed significant changes in the SRM analyses in agreement with the label-free quantification (Fig. 2a-e, Additional file 2: Table S1). Three other SRM-quantified peptides selected for the APOA4 protein also showed significant decreases (Additional file 1: Figure S3A-C), similar to the 4th peptide shown in Fig. 2d. The peptide for cathepsin Z (CATZ), which was found to have increased levels in MeHg-treated samples with the label-free method (Additional file 2: Table S1), did not change significantly using the SRM method (Additional file 1: Figure S3D). Thus, the SRM data showed largely good concordance with the label-free quantitation, confirming reliability of the quantitative data generated by the latter method.
Quantitative proteomics studies may help to identify proteins that can be used as biomarkers in mechanistic toxicological studies and environmental pollution monitoring. In this regard, using label-free quantitative proteomics in the discovery phase for identifying altered proteins, followed by verification with SRM, can be a useful method in identification and quantification of potential biomarker proteins. Verification of differentially regulated candidate biomarkers by the SRM method is particularly useful when antibodies are unavailable for quantification of the proteins of interest by immunological methods, such as Western blotting and ELISA.
A more promising biomarker candidate of particular interest is the enzyme uroporphyrinogen decarboxylase (DCUP/UROD), found down regulated here (Additional file 2: Table S1). UROD is an enzyme in the heme biosynthetic pathway and decreased activity of DCUP has been associated with the human diseases porphyria and liver damage . It has been shown that MeHg inhibits DCUP in rat kidney, leading to increased porphyrin excretion . The urinary porphyrin excretion profile was specific to mercury compounds, and thus it was suggested as a specific biomarker of mercury exposure . Therefore, down regulation of UROD in MeHg treated cod makes it a strong candidate as a biomarker of MeHg exposure. Liver and blood UROD activity and porphyrin profile might be evaluated as biomarkers for MeHg exposure in fish. Down regulation of other proteins involved in heme metabolism such as heme-binding protein 2 (HEBP2) and biliverdin reductase B (BLVRB) (Additional file 2: Table S1) also suggest perturbation of heme homeostasis in MeHg treated cod liver. Consistent with our results, lower liver UROD activities and increased porphyrin levels were detected in pike from River Rhine polluted with halogenated hydrocarbons and heavy metals including mercury . However, UROD activity and porphyrin profiles are not commonly used as biomarkers in fish , and our results encourage further potential for development. Aryl hydrocarbon receptor agonist halogenated hydrocarbons have been shown to induce prorphyrin accumulation in the fish hepatoma cell line (PLHC-1) . Thus, in further experiments, changes in expression and activity of UROD combined with porphyrin profile may be performed to evaluate the specificity of responses to mercury exposure in hepatoma cell line  and in ex vivo systems such as precision-cut liver slices .
Enrichment of the energy pathways supports the effects of MeHg on these cellular processes, particularly mitochondrial metabolic pathways as observed previously by gene expression analysis of the same samples . In the previous gene expression analysis, transcripts of many mitochondrial enzymes, such as acyl-CoA dehydrogenases were particularly affected, suggesting specific effects of MeHg-induced oxidative stress. Recently, Cys-residues of many mitochondrial enzymes have been shown to be oxidized upon Cadmium-induced oxidative stress in mouse liver . Notably, Cys residues of some of the proteins differentially regulated here (ACADS, GLNA, IDHP/IDH2, ERP44, PECR and TRFE) (Additional file 2: Table S1), were found to be differentially oxidized in the Cd-treated mouse liver samples . Among these, ACADS, IDHP/IDH2 and PECR are enzymes that are part of cellular energy pathways. Importantly, ERP44 is a multifunctional chaperone with a Cys-residue that is critical for its function in quality control of proteins in the early secretory pathway, redox homeostasis, and regulation of calcium flux . In our previous transcriptome analysis of the same samples, we also observed differential regulation of several genes participating in different energy pathways, particularly the fatty acid and amino acid metabolic pathways . The similar responses to cadmium chloride and MeHg suggest that differential regulation of the genes or proteins are mediated by induction of oxidative stress by these compounds.
Mitochondrial energy generating pathways (electron-transport chain and fatty acid beta-oxidation) were found to be affected by HgCl2 in zebrafish liver , which is also in agreement with our results. Similarly, proteins that are part of the antioxidant defense and energy pathways were modulated in Atlantic salmon fed with a MeHg-containing diet . Moreover, energy and redox pathways were also affected in the brain, demonstrated by two-dimensional gel electrophoresis analysis of the brain proteome of Atlantic cod  as well as transcriptomics studies of zebrafish . Recently, by proteomics analysis, it was also reported that chronic exposure to low doses of MeHg resulted in down-regulation of energy pathway proteins and metabolic deficit in the somatosensory cortex of the rat brain . In another recent proteomics study of the brain of common marmosets, proteins involved in the metabolism of carbohydrate derivatives were found to be down regulated, suggesting negative effects of MeHg on the carbohydrate metabolism . Taken together, these results strongly suggest that oxidative stress mediated effects on energy pathways is a general mechanism of MeHg-mediated toxicity.
Comparison of proteomics and gene expression data
Comparison of significantly enriched (FDR < 0.05) top 50 annotations (GO BP, MF and MetaCore process networks) between the trancriptomics  and proteomics data obtained from the same cod liver samples showed many shared pathways (Additional file 1: Figure S5). As expected, a higher number of genes (650 genes) were differentially regulated in the microarray analysis compared to the identified number of differentially regulated proteins (125 proteins). However, twenty-six annotations were shared between the two datasets. The main processes previously known to be affected by MeHg, such as energy metabolism, oxidoreductase activities, protein translation, and cytoskeleton remodelling, were enriched in both datasets, showing agreement between the two methods.
Although analysis of transcriptomics and proteomics data revealed many shared pathways, only 11 genes were present in both lists of differentially regulated genes and proteins. Of these 11 genes, 6 showed concordance in expression changes at mRNA and protein levels (genes encoding ACADS, METK2, SCOT1, EPS8L2 and FUBP2 increased, while CATZ decreased). In contrast, genes encoding IDHP/IDH2, GSTT1, RL19 and TRI16 increased at the mRNA levels and decreased at the protein levels. The gene encoding LPP60 decreased at mRNA level and increased at protein level. Thus, changes in 5 of the 11 proteins were not consistent with changes in mRNA levels. The overall correlation for 32 common transcripts/proteins detected at more relaxed differential regulation cut-offs for transcripts (FDR < 5 % and minimum 1.2 fold change) and proteins (p < 0.05) was low (r = 0.16), with only 14 of them changing in the same direction (not shown). The lack of strong correlation between transcript and protein changes is not surprising. Numerous studies in various organisms reported low to modest correlation between mRNA and protein levels [69–71]. In mouse liver, not more than 50 % of the genes showed a significant correlation of 0.27 between transcript and protein levels . In a study comparing transcript and protein expression responses to a toxicant (2A-DNT) treatment, poor correlation was observed, although higher order pathways affected were similar . Similarly, in our study proteomics and transcriptomics datasets shared many GO annotations despite low correlation at the individual transcript and protein expression levels (Additional file 1: Figure S5). Factors that may contribute to the low correlation may be post-translational modifications, localization (ex. secreted proteins), and different dynamics of mRNA and protein turnover. As discussed above, degradation of some proteins triggered by MeHg-induced oxidation could be one factor contributing to low correlation between protein and mRNA levels.
To further investigate integrated effects of MeHg in the cod liver, we combined the differentially regulated 650 genes and the 125 proteins and performed network analysis and enrichment of GO (BP and MF), KEGG, Wikipathways, and REACTOME pathways using the CluePedia application in Cytoscape. The generated network highlights major processes and pathways affected, including amino-acid metabolism, fatty-acid metabolism, glucose metabolism and protein synthesis (Additional file 1: Figure S6). The energy pathways dominate in the combined dataset due to the high number of energy pathway genes, particularly amino acid and fatty acid metabolism genes represented in the larger transcriptomics dataset . Processes and pathways mostly related to oxidative stress such as sulphur compound metabolism, glutathione metabolism and oxidation and reduction processes, and co-factor metabolism, are also highlighted (Additional file 1: Figure S6). Notably, the pentose phosphate pathway (PPP) that generates NADPH required as electron donor in the antioxidant systems is significantly enriched in the combined dataset (Additional file 1: Figure S6).
Finally, comparison of GO BP and MF affected in MeHg treated cod as determined by liver transcriptomics analysis , liver proteomics (this study) and brain proteomics (same fish as used in this study)  shows that processes related to energy metabolism and oxidative responses are shared among the three studies (Additional file 1: Figure S7). For example, among the 12 significantly enriched GO terms shared between the datasets, “small molecule metabolic process” and “oxidoreductase activity” are related to energy and redox pathways (Additional file 1: Figure S7).
Other enriched pathways include pathways related to translation, cytoskeleton, vesicle transport and cholesterol transport (Table 1; Additional file 4: Table S3, Additional file 1: Figure S4). Effects on protein synthesis and the cytoskeleton were also observed in both transcriptomics and proteomics analyses (Additional file 4: Table S3, Additional file 1: Figure S4). It is not clear if some of the enriched pathways are activated or suppressed from our analysis as the differential regulation of component proteins largely changed in both directions (Additional file 1: Figure S4). For example, some proteins involved in poly(A) RNA binding related to translation increased and others decreased, and the outcome on protein synthesis is not clear. However, previous studies have reported that MeHg inhibits protein synthesis . The enrichment of the cytoskeletal assembly pathway is possibly related to the known negative effect of MeHg on microtubule assembly . Tubulin alpha and beta chains were also among the polypeptides with reactive thiol groups identified in the liver thiol proteome analysis of the same cod liver samples .
Network analyses suggest possible interactions of oxidative stress and energy pathways
Network analysis of the differentially regulated proteins was performed in Cytoscape (http://www.cytoscape.org/) using the Cluepedia plug-in . Similar to the GO terms and process networks enriched in the MetaCore analysis (Table 1, Additional file 4: Table S3), the generated networks highlighted energy pathways (e.g. TCA cycle and fatty acid beta oxidation) and antioxidant responses (e.g. phase II conjugation) (Fig. 5). Other pathways related to the significantly enriched GO terms “vesicle-mediated transport” and “blood coagulation” and “cadherin binding” (Additional file 4: Table S3) can also be visualized in the Cytoscape network (Fig. 5).
Furthermore, protein-protein interaction network analysis of the differentially regulated proteins was performed using the STRING database and web-tool (http://string-db.org). STRING also enables enrichment analysis and visualization of KEGG, GO and protein domains in protein networks. Annotations consisting of 6 KEGG pathways, 9 GO BP terms, 2 GO MF terms, 16 GO CC terms and 1 Interpro protein domain were significantly enriched at FDR p-value < 0.05 (Additional file 1: Figure S4, Additional file 5: Table S4).
The networks suggest interesting interactions of proteins and pathways of energy metabolism and antioxidant systems. Notably, the networks show separate clustering of energy pathway (TCA cycle) proteins and antioxidant response proteins (disulfide oxido-reductase activity, phase II conjugation) linked by the mitochondrial isocitrate dehydrogenase enzyme (IDPH/IDH2) node (Fig. 5). IDH2, that was down regulated in MeHg exposed samples (Additional file 2: Table S1), is an energy pathway enzyme that is also a key component of the antioxidant defense, being important in the NADPH generation in the mitochondria . NADPH is important in providing reducing potential to the glutaredoxin and thioredoxin antioxidant systems. Although the consequence of IDH2 down regulation on energy pathways is unclear from our results, mitochondrial IDH2 function is known to be influenced by ROS , and it has been suggested that its activity is specifically modulated through oxidation of an important Cys-residue upon oxidative stress . Consistent with these previous reports, our network analysis suggests that IDH2 could be one of the links for interactions between antioxidant systems and energy pathways (Fig. 5). Indeed, recent studies have shown that ROS levels can regulate the activities of enzymes in the mitochondrial energy generating pathways through Cys thiol modifications .
Further, STRING protein-protein interaction networks and enrichment analysis also suggest crosstalk between energy pathways and redox homeostasis (Additional file 1: Figure S4). For example, ATP-citrate lyase (ACLY) that is known to be redox-regulated , forms an important hub in STRING protein-protein interaction network linked with many energy pathway and antioxidant proteins (Additional file 1: Figure S4). Another hub is formed by GAPDH, which is a glycolytic enzyme known to be involved in antioxidant responses . ACLY, GAPDH and many energy pathway enzymes (e.g. TCA cycle and fructose and mannose metabolism) were down regulated in MeHg treated samples, suggesting suppression of the cellular energy pathways (Additional file 1: Figure S4, Additional file 2: Table S1). Suppression of the energy pathways is consistent with presence of reactive Cys residues in many of the energy pathway enzymes discussed above that might make them targets of MeHg-indued oxidation. Importantly, it has been reported that oxidation of GAPDH enables the re-routing of metabolic flux from glycolysis to the Pentose Phosphate Pathway (PPP) in order to generate NADPH required for the antioxidant systems . This is another example of redox modulation of an energy pathway. The PPP is one of the significantly enriched carbohydrate metabolism pathways in the combined transcriptomics and proteomics dataset (Additional file 1: Figure S6), suggesting its activation in response to increased oxidative stress. These results provide some mechanistic insights into how MeHg-induced oxidative stress might be disrupting cellular redox balance and energy metabolism.
Recent studies have identified oxidation sensitive Cys residues in many mitochondrial enzymes and suggest modulatory effects of oxidative stress on the energy metabolism [40, 79], which is consistent with our results. It has been suggested that the interaction of redox and metabolic pathways are coordinated by the transcription factor nuclear factor erythroid 2-related factor (Nrf2), which regulates transcription of genes in phase II and phase III detoxification pathways, the antioxidant glutaredoxin and thioredoxin systems, and energy pathways, including carbohydrate metabolism and NADPH generation . In fact, many Nrf2 regulated genes were regulated in our MeHg treated samples, and the Nrf2 regulatory pathway was significantly enriched in the transcriptomics analysis . Long established effects of MeHg, such as mitochondrial dysfunction, altered calcium homeostasis, and apoptosis, are thought to be consequences of its oxidative stress promoting characteristics [12, 13]. Although other studies have reported effects of MeHg on energy pathways [22, 23, 27, 81], our integrated transcriptomics and proteomics data indicate more comprehensive effects of MeHg on major energy pathways, in particular the mitochondrial fatty acid metabolism. The apparently more coordinated effects observed on energy metabolic pathways in our studies, perhaps reflect the extent of the data (integrated transcriptomics and proteomics) and the fact that the main physiological role of the liver is metabolism. The overall results indicate that oxidative stress may lead to perturbations in energy pathways in cod liver, further supporting induction of oxidative stress as the underlying mechanism of cytotoxic effects of MeHg. Although oxidation of proteins is often reversible, persistent exposure to MeHg may lead to irreversible oxidation of thiol groups causing enzyme inactivation, mitochondrial dysfunction and other cellular damages.
Extensive proteomics analysis of liver samples from cod treated with MeHg indicated oxidative stress responses and perturbation of cellular energy pathways, confirming our earlier findings by transcriptomics analysis. The integrated analyses with both transcriptomics and proteomics data, together with previous studies, further support that MeHg-induced oxidative stress modulates cellular energy pathways, contributing to well-known toxic effects of this compound, such as mitochondrial dysfunction and other cellular damages.
Fish exposure and sampling
Fish exposure and sampling has been described previously . Briefly, juvenile Atlantic cod (G. morhua) of mixed gender (body weight of 260–530 g) were divided into four 500 L-tanks for control, 0.5, 2 and 8 mg/kg groups (n = 10/group), and acclimated for 6 days in continuously running seawater (temperature 10 ° C and salinity 34‰) at the facility in ILAB, Bergen, Norway. The fish were fed daily with commercial pellets (Ewos marin 6 mm, EWOS, Bergen, Norway). Methylmercury chloride (CH3 HgCl) (Strem Chemicals, Newburyport, USA) was dissolved in the vehicle (20 % acetone and 80 % soybean oil). The vehicle (control), 0.5, 2 or 8 mg/kg body weight methylmercury chloride was injected intraperitoneally by inserting the needle carefully into the abdominal cavity, not to puncture internal organs. The indicated doses were divided into two and the first half injected on day 1, and the second half after one week. After the second injection all the fish in the highest (8 mg) dose group and 1 fish from each of the 0.5 and 2 mg dose groups died. The fish were observed and behavioural parameters such as overall activity and balance was registered during the exposure period. After 2 weeks (14 days) the fish were sacrificed by decapitation and tissue samples collected. Liver samples were dissected, weighed and snap frozen in liquid nitrogen and stored at −80 °C. The National Animal Research Authority approved the exposure experiment. Fish length (in cm) and weight (in g) were recorded at the start and end of the experiment. Condition factor (100 * weight/length3) and liver-somatic index (100 * liver weight/BW) were calculated, and Student's t-test was used to calculate p-value between sample groups.
Homogenization of liver samples has been described previously . Briefly, frozen individual cod liver sample obtained from each of control (n = 10), 0.5 mg/kg BW MeHg (n = 9) and 2 mg/kg BW MeHg (n = 9) groups (from 28 fish in total) was thawed on ice, washed and homogenized in phosphate-buffered saline (PBS, pH 7.4). The homogenate was centrifuged at 13,000 × g for 10 min and the supernatant containing hepatic proteins was digested with trypsin, and the resulting peptides purified as described previously .
Label-free proteomics using LC-MS
The peptide sample (1 μg) from each of control (n = 10), 0.5 mg/kg BW MeHg (n = 9) and 2 mg/kg BW MeHg (n = 9) groups was dissolved in 1 % aqueous formic acid (FA) and injected into an Ultimate 3000 RSLC system (Thermo Scientific) connected online to an LTQ-Orbitrap Velos Pro mass spectrometer (Thermo Scientific) equipped with a nanospray Flex ion source (Thermo Scientific). The peptides were separated using a 90 min HPLC gradient with increasing acetonitrile concentration in 0.1 % FA and the acquisition time for LC-MS data was 80 min as previously described . The LTQ-Orbitrap Velos Pro settings were the same as previously described  except that the normalized collision energy was set to 35 % and the dynamic exclusion to 30s.
Analysis of the label-free data using Progenesis LC-MS
Progenesis LC-MS® v4.1 (Nonlinear Dynamics Ltd) was used for label-free quantification of automatically aligned LC-MS data. Only peptide features with charges between +2 to +5 containing the top 10 ranked MSMS spectra (with precursor intensities more than 0) and where MSMS was executed on the highest precursor isotopes (less than 6) were accepted. The MSMS spectra list was sorted by rank and the fragment ions in the MSMS spectra were limited to 200 allowing deisotoping and charge deconvolution prior to export as an mgf file for identification. The mgf file was searched using SearchGUI v1.18.1  using the option “Delete duplicate spectra titles”. The search engines X!Tandem, MS Amanda and MS-GF+ Beta were used for searches in the Ensembl Gadus morhua protein database. Protein identification was conducted against a concatenated target/decoy version of a complement of the database and the decoy sequences were created by reversing the target sequences in SearchGUI. The identification settings were: trypsin with a maximum of 1 missed cleavages; 10 ppm as MS1 and 0.5 Da as MS2 tolerances; carbamidomethyl Cys as fixed modification; oxidation of Met as variable modification.
Peptides and proteins were inferred from the spectrum identification results using PeptideShaker v0.28 . Peptide Spectrum Matches (PSMs), peptides and proteins were validated at a 1 % False Discovery Rate (FDR) estimated using the decoy hit distribution. The results were exported as validated PSMs and imported into Progenesis. The sum of the normalized abundances of 1143 quantified unique proteins was exported from Progenesis and analyzed further using statistical analyses.
Proteomics data analysis
For functional annotation and pathway analysis, the 1143 unique cod proteins quantified (Additional file 3: Table S2) were mapped to orthologs of human proteins in Swissprot database using the BioMart tool in Ensembl database (Ensembl.org) or BLAST searches in Swissprot database using the cod protein sequences. There are 22154 genes in the genome of Atlantic cod in the Ensembl database (Ensembl.org). Differential expression analysis was performed using one-way ANOVA followed by Dunnett´s test (JMP, SAS Institute, Cary NC, USA). In addition, we performed differential expression analysis using SAM (Significance Analysis of Microarrays) implemented in J-Express  (http://www.molmine.com/) that enables estimation of False Discovery Rates (FDR) in pairwise comparisons. Proteins with p value < 0.05 (Dunnett´s test) and FDR < 0.1 (SAM) satisfying a fold change (treated/ control) greater than 1.2 or less than 0.8 were considered differentially regulated. (Additional file 2: Table S1). Although minimum fold change of 1.3-2 is commonly used as cut-off for differential regulation in quantitative proteomics , a minimum fold change of 1.2 was used here, considering a relatively large number of biological replicates analyzed (n = 9–10 per group). Following these criteria, a total of 125 proteins were identified as significantly differentially regulated between the 2 mg/ kg BW MeHg dose and control groups, of which 44 were up-regulated and 81 were down regulated between controls and 2 mg/kg BW MeHg treated samples (Additional file 2: Table S1). Among the 125 proteins, 7 were differentially regulated in both 0.5 mg and 2 mg/kg BW MeHg dose groups (Additional file 2: Table S1) using the statistical criteria stated above. No proteins were differentially regulated exclusively in the low dose (0.5 mg/ kg BW) group. Hierarchical clustering (complete-linkage) and principal component analysis (PCA) of proteins differentially regulated by MeHg was performed based on log2-transformed abundance ratio values and arcSinh normalised protein abundances, respectively (JMP 11, SAS Institute, Cary NC, USA). Discriminant analyses were performed on arcSinh normalised protein abundances using the Linear method and stepwise variable selection in JMP 11.
Selected reacting monitoring (SRM)
SRM was performed for six proteins (T23O, GLNA EPS8L2, APOA4, RAP1B, CATZ) among the top differentially regulated. As certain criteria, such as peptide length, protease miscleavage, and amino acid composition must be considered when selecting peptides for SRM analyses, it was not possible to choose freely from the list of proteins/peptides quantified with the label-free approach. Stable isotope –labelled synthetic peptides (SIS) were purchased in crude quality from Thermo Scientific. The C-terminal lysine or arginine in the SIS peptides were replaced by isotope-labelled lysine (13C6, 15N2) or arginine (13C6, 15N4), resulting in a mass difference of 8 Da and 10 Da, respectively, to the corresponding endogenous peptide. The SIS peptides were spiked in equal amounts into the digested protein samples, at approximately the same level as the endogens peptide, prior to desalting with Oasis HLB μElution Plate (Waters). A Q-Trap 5500 (AB SCIEX) connected to a Dionex Ultimate NCR-3500RS LC system was used for the LC SRM analyses. The protein digest were dissolved in 2 % ACN, 0.1%FA (Buffer A) and 1 μg peptides was trapped on the pre-column (Dionex, Acclaim PepMap 100, 2 cm x 75 μm i.d, 3 μm C18 beads) in buffer A at a flowrate of 5 μl/min for 5 min before separation by reverse phase chromatography (Dionex, Acclaim PepMap 100, 15 cm x 75 μm i.d., 3 μm C18 beads) at a flow of 250 nL/min. The nano LC gradient consisted of a linear gradient starting at 5 % of 90 % ACN, 0.1 % FA (buffer B) and ramping to 40 % buffer B over 45 min (0–45 min). In half a minute the gradient was ramped to 90 % buffer B (45–45.5 min) and held for 6 min (45.5-51.5 min.) followed by ramping to 5 % buffer B for 3.5 min (51.5-55 min) and equilibration of the column in 15 min (55–70 min). The Collision energy (CE) for each peptide was calculated according to the following formulas; CE = 0.044 x m/z +10 for doubly charged precursors and CE = 0.05 x m/z + 9 for triply precursors. The SRM data was analyzed using Skyline v1.4  using the most abundant transition for quantification. Student's t-test was used to calculate the p-value between sample groups.
Glutathione-S-transferase (GST) assay
GST activity was determined according to a previously published method . GST activity per min per mg protein was calculated for each sample in control (n = 10), 0.5 mg/kg BW MeHg (n = 9) and 2 mg/kg BW MeHg (n = 9) groups. Treated and control groups were compared using Student’s t-test.
Annotation, pathway and gene ontology enrichment analysis
Enrichment analyses for MetaCore “Process Networks” and GO (Biological Processes and Molecular Functions) were performed in MetaCore™ (MetaCore™; GeneGo-Thomson Reuters, St. Joseph, MI). A false discovery rate (FDR) < 0.05 was used as the threshold for significant enrichment.
Network analysis, enrichment and visualization
Gene ontology pathway enrichment analysis and network visualization was performed using Cytoscape v3.2.1 (http://www.cytoscape.org/) with the plug-ins ClueGo v2.1.7  and CluPedia v1.1.7 , at FDR < 0.05. GO (BP and MF) and KEGG, REACTOME and WikiPathways databases were used. Protein-protein interaction network and KEGG pathway, GO and protein domain enrichment analyses were performed in STRING v10 .
ACADS, Short-chain specific acyl-CoA dehydrogenase, mitochondrial; ACLY, ATP citrate lyase; APOA4, Apolipoprotein A4; AT1B1, ATPase, Na+/K+ transporting, beta 1 polypeptide; BP, biological process; BW, Body weight; CATZ, Cathepsin Z; CC, cellular component; DHSO, sorbitol dehydrogenase; EPS8L2, Epidermal growth factor receptor kinase substrate 8-like protein 2; ERP44, Endoplasmic reticulum resident protein 44; FDR, false discovery rate; FUBP2, Far upstream element-binding protein 2; GAPDH (G3P), glyceraldehyde-3-phosphate dehydrogenase; GAPDH; GSTO1, glutathione S-transferase omega-1: GLNA, glutamine synthetase; GLNA, glutamine synthetase; GLRX, Glutaredoxin-1; GO, Gene Ontology; GRP78, 78 kDa glucose-regulated protein; GSTT1, glutathione S-transferase theta-1; GSTT2, glutathione S-transferase theta-2B; HXK2, Hexokinase-2; IDHP (IDH2), Isocitrate dehydrogenase [NADP], mitochondrial; KEGG, Kyoto Encyclopedia of Genes and Genomes; LPP60, 60 kDa lysophospholipase; MeHg, methylmercury; METK2:, methionine adenosyltransferase II, alpha; MF; molecular function; NADPH, Nicotinamide adenine dinucleotide phosphate; Nrf2, nuclear factor erythroid 2-related factor; PECR, Peroxisomal trans-2-enoyl-CoA reductase; PPP, Pentose Phosphate Pathway; PRDX2, Peroxiredoxin-2; RL19, 60S Ribosomal Protein L19; RPIA, Ribose-5-phosphate isomerase; SCOT1, Succinyl-CoA; SRM, Selected reacting monitoring; T23O (TDO2), tryptophan 2,3 dioxygenase; TRFE, Serotransferrin; TRI16, Tripartite motif-containing protein 16; TXNL1, Thioredoxin-like protein 1; UPR, Unfolded Protein Response.
The authors are grateful to the Norwegian Cod Genome Sequencing Consortium for providing access to the cod EST and genome assemblies ahead of the publication of the cod genome sequence.
This work was funded by the Norwegian Research Council grants iCOD (192441/I30) and iCOD 2.0 (244564).
Availability of data and materials
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE  partner repository with the dataset identifier PXD004311. The transcriptomics dataset was previously deposited in NCBI’s Gene Expression Omnibus (GEO) database (GEO accession GSE38746).
FY performed data analysis and wrote the manuscript. SB participated in sample preparation and assays. HKG, EO, FB performed proteomics data acquisition and reviewed the manuscript. AG conceived the study, participated in its design, coordinated the study, and reviewed the manuscript. OAK coordinated the project, participated in proteomics experiment design, data analysis, writing and review of the manuscript. All authors read and approved the final manuscript.
The authors declare no competing interests in relation to this work.
Consent for publication
Ethics approval and consent to participate
The Atlantic cod MeHg exposure experiment was conducted in compliance with the European Directive 2010/63/EU, and approved by the National Animal Research Authority (NARA) in Norway.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Clarkson TW. The toxicology of mercury. Crit Rev Clin Lab Sci. 1997;34(4):369–403.View ArticlePubMedGoogle Scholar
- Jensen S, Jernelov A. Biological methylation of mercury in aquatic organisms. Nature. 1969;223(5207):753–4.View ArticlePubMedGoogle Scholar
- Birch RJ, Bigler J, Rogers JW, Zhuang Y, Clickner RP. Trends in blood mercury concentrations and fish consumption among US women of reproductive age, NHANES, 1999–2010. Environ Res. 2014;133:431–8.View ArticlePubMedGoogle Scholar
- Clarkson TW, Magos L. The toxicology of mercury and its chemical compounds. Crit Rev Toxicol. 2006;36(8):609–62.View ArticlePubMedGoogle Scholar
- Syversen T, Kaur P. The toxicology of mercury and its compounds. J Trace Elem Med Biol. 2012;26(4):215–26.View ArticlePubMedGoogle Scholar
- Simmons-Willis TA, Koh AS, Clarkson TW, Ballatori N. Transport of a neurotoxicant by molecular mimicry: the methylmercury-L-cysteine complex is a substrate for human L-type large neutral amino acid transporter (LAT) 1 and LAT2. Biochem J. 2002;367:239–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Bakir F, Damluji SF, Aminzaki L, Murtadha M, Khalidi A, Alrawi NY, Tikriti S, Dhahir HI, Clarkson TW, Smith JC, et al. Methylmercury poisoning in Iraq –Interuniversity report. Science. 1973;181(4096):230–41.View ArticlePubMedGoogle Scholar
- Harada M. Minamata disease: methylmercury poisoning in Japan caused by environmental pollution. Crit Rev Toxicol. 1995;25(1):1–24.View ArticlePubMedGoogle Scholar
- Drevnick PE, Roberts AP, Otter RR, Hammerschmidt CR, Klaper R, Oris JT. Mercury toxicity in livers of northern pike (Esox lucius) from Isle Royale, USA. Comp Biochem Physiol C-Toxicol Pharmacol. 2008;147(3):331–8.View ArticlePubMedGoogle Scholar
- Mela M, Randi MAF, Ventura DF, Carvalho CEV, Pelletier E, Ribeiro CAO. Effects of dietary methylmercury on liver and kidney histology in the neotropical fish Hoplias malabaricus. Ecotoxicol Environ Saf. 2007;68(3):426–35.View ArticlePubMedGoogle Scholar
- Lee H, Kim Y, Sim CS, Ham JO, Kim NS, Lee BK. Associations between blood mercury levels and subclinical changes in liver enzymes among South Korean general adults: Analysis of 2008–2012 Korean national health and nutrition examination survey data. Environ Res. 2014;130:14–9.View ArticlePubMedGoogle Scholar
- Ceccatelli S, Dare E, Moors M. Methylmercury-induced neurotoxicity and apoptosis. Chem Biol Interact. 2010;188(2):301–8.View ArticlePubMedGoogle Scholar
- Farina M, Aschner M, Rocha JBT. Oxidative stress in MeHg-induced neurotoxicity. Toxicol Appl Pharmacol. 2011;256(3):405–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Limke TL, Bearss JJ, Atchison WD. Acute exposure to methylmercury causes Ca2+ dysregulation and neuronal death in rat cerebellar granule cells through an M3 muscarinic receptor-linked pathway. Toxicol Sci. 2004;80(1):60–8.View ArticlePubMedGoogle Scholar
- Vogel DG, Margolis RL, Mottet NK. The effects of methyl mercury binding to microtubules. Toxicol Appl Pharmacol. 1985;80(3):473–86.View ArticlePubMedGoogle Scholar
- Yin ZB, Milatovic D, Aschner JL, Syversen T, Rocha JBT, Souza DO, Sidoryk M, Albrecht J, Aschner M. Methylmercury induces oxidative injury, alterations in permeability and glutamine transport in cultured astrocytes. Brain Res. 2007;1131(1):1–10.View ArticlePubMedGoogle Scholar
- Branco V, Canario J, Holmgren A, Carvalho C. Inhibition of the thioredoxin system in the brain and liver of zebra-seabreams exposed to waterborne methylmercury. Toxicol Appl Pharmacol. 2011;251(2):95–103.View ArticlePubMedGoogle Scholar
- Carvalho CML, Chew E-H, Hashemy SI, Lu J, Holmgren A. Inhibition of the human thioredoxin system - A molecular mechanism of mercury toxicity. J Biol Chem. 2008;283(18):11913–23.View ArticlePubMedGoogle Scholar
- Franco JL, Posser T, Dunkley PR, Dickson PW, Mattos JJ, Martins R, Bainy ACD, Marques MR, Dafre AL, Farina M. Methylmercury neurotoxicity is associated with inhibition of the antioxidant enzyme glutathione peroxidase. Free Radic Biol Med. 2009;47(4):449–57.View ArticlePubMedGoogle Scholar
- Van Aggelen G, Ankley GT, Baldwin WS, Bearden DW, Benson WH, Chipman JK, Collette TW, Craft JA, Denslow ND, Embry MR, et al. Integrating Omic Technologies into Aquatic Ecological Risk Assessment and Environmental Monitoring: Hurdles, Achievements, and Future Outlook. Environ Health Perspect. 2010;118(1):1–5.PubMedGoogle Scholar
- Martyniuk CJ, Alvarez S, Denslow ND. DIGE and iTRAQ as biomarker discovery tools in aquatic toxicology. Ecotoxicol Environ Saf. 2012;76:3–10.View ArticlePubMedGoogle Scholar
- Kong HK, Wong MH, Chan HM, Lo SCL. Chronic Exposure of Adult Rats to Low Doses of Methylmercury Induced a State of Metabolic Deficit in the Somatosensory Cortex. J Proteome Res. 2013;12(11):5233–45.View ArticlePubMedGoogle Scholar
- Shao YT, Yamamoto M, Figeys D, Ning ZB, Chan HM. Proteomic Analysis of Cerebellum in Common Marmoset Exposed to Methylmercury. Toxicol Sci. 2015;146(1):43–51.View ArticlePubMedGoogle Scholar
- Berg K, Puntervoll P, Valdersnes S, Goksøyr A. Responses in the brain proteome of Atlantic cod (Gadus morhua) exposed to methylmercury. Aquat Toxicol. 2010;100(1):51–65.View ArticlePubMedGoogle Scholar
- Karlsen OA, Sheehan D, Goksøyr A. Alterations in the Atlantic COD (Gadus morhua) Hepatic Thiol-Proteome After Methylmercury Exposure. J Toxicol Environ Health Part A Curr Issues. 2014;77(9–11):650–62.View ArticleGoogle Scholar
- Nostbakken OJ, Martin SAM, Cash P, Torstensen BE, Amlund H, Olsvik PA. Dietary methylmercury alters the proteome in Atlantic salmon (Salmo salar) kidney. Aquat Toxicol. 2012;108:70–7.View ArticlePubMedGoogle Scholar
- Richter CA, Garcia-Reyero N, Martyniuk C, Knoebl I, Pope M, Wright-Osment MK, Denslow ND, Tillitt DE. Gene expression changes in female zebrafish (Danio rerio) brain in response to acute exposure tomethylmercury. Environ Toxicol Chem. 2011;30(2):301–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Amlund H, Lundebye A-K, Berntssen MHG. Accumulation and elimination of methylmercury in Atlantic cod (Gadus morhua L.) following dietary exposure. Aquat Toxicol. 2007;83(4):323–30.View ArticlePubMedGoogle Scholar
- Balk L, Hylland K, Hansson T, Berntssen MHG, Beyer J, Jonsson G, Melbye A, Grung M, Torstensen BE, Børseth JF, et al. Biomarkers in Natural Fish Populations Indicate Adverse Biological Effects of Offshore Oil Production. Plos One. 2011;6(5):e19735.View ArticlePubMedPubMed CentralGoogle Scholar
- Chesman BS, O'Hara S, Burt GR, Langston WT. Hepatic metallothionein and total oxyradical scavenging capacity in Atlantic cod Gadus morhua caged in open sea contamination gradients. Aquat Toxicol. 2007;84(3):310–20.View ArticlePubMedGoogle Scholar
- Goksøyr A, Andersson T, Hansson T, Klungsøyr J, Zhang Y, Förlin L. Species characteristics of the hepatic xenobiotic and steroid biotransformation systems of two teleost fish, Atlantic cod (Gadus morhua) and rainbow trout (Salmo gairdneri). Toxicol Appl Pharmacol. 1987;89(3):347–60.View ArticlePubMedGoogle Scholar
- Holth TF, Beylich BA, Skarphedinsdottir H, Liewenborg B, Grung M, Hylland K. Genotoxicity of Environmentally Relevant Concentrations of Water-Soluble Oil Components in Cod (Gadus morhua). Environ Sci Technol. 2009;43(9):3329–34.View ArticlePubMedGoogle Scholar
- Olsvik PA, Lie KK, Goksøyr A, Midtun T, Frantzen S, Maage A. Are Atlantic Cod in Store Lungegrdsvann, a Seawater Recipient in Bergen, Affected by Environmental Contaminants? A qRT-PCR Survey. J Toxicol Environ Health Part A Curr Issues. 2009;72(3–4):140–54.View ArticleGoogle Scholar
- Tollefsen KE, Sundt RC, Beyer J, Meier S, Hylland K. Endocrine Modulation in Atlantic Cod (Gadus morhua L.) Exposed to Alkylphenols, Polyaromatic Hydrocarbons, Produced Water, and Dispersed Oil. J Toxicol Environ Health Part A Curr Issues. 2011;74(7–9):529–42.View ArticleGoogle Scholar
- Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrom M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, et al. The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011;477(7363):207–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Yadetie F, Karlsen OA, Eide M, Hogstrand C, Goksøyr A. Liver transcriptome analysis of Atlantic cod (Gadus morhua) exposed to PCB 153 indicates effects on cell cycle regulation and lipid metabolism. BMC Genomics. 2014;15:481.View ArticlePubMedPubMed CentralGoogle Scholar
- Yadetie F, Karlsen OA, Lanzen A, Berg K, Olsvik P, Hogstrand C, Goksøyr A. Global transcriptome analysis of Atlantic cod (Gadus morhua) liver after in vivo methylmercury exposure suggests effects on energy metabolism pathways. Aquat Toxicol. 2013;126:314–25.View ArticlePubMedGoogle Scholar
- St Gelais AT, Costa-Pierce BA. Mercury concentrations in Northwest Atlantic winter-caught, male spiny dogfish (Squalus acanthias): A geographic mercury comparison and risk-reward framework for human consumption. Mar Pollut Bull. 2016;102(1):199–205.View ArticlePubMedGoogle Scholar
- EU: Setting maximum levels for certain contaminants in foodstuffs. Off J Eur Union 2006, L 364/5Google Scholar
- Go YM, Roede JR, Orr M, Liang YL, Jones DP. Integrated Redox Proteomics and Metabolomics of Mitochondria to Identify Mechanisms of Cd Toxicity. Toxicol Sci. 2014;139(1):59–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu W, Xu ZF, Deng Y, Xu B, Yang HB, Wei YG, Feng S. Excitotoxicity and Oxidative Damages Induced by Methylmercury in Rat Cerebral Cortex and the Protective Effects of Tea Polyphenols. Environ Toxicol. 2014;29(3):269–83.View ArticlePubMedGoogle Scholar
- Kaliman PA, Nikitchenko IV, Sokol OA, Strel'chenko EV. Regulation of heme oxygenase activity in rat liver during oxidative stress induced by cobalt chloride and mercury chloride. Biochem Moscow. 2001;66(1):77–82.View ArticleGoogle Scholar
- Panayiotidis MI, Stabler SP, Ahmad A, Pappa A, Legros LH, Hernandez-Saavedra D, Schneider BK, Allen RH, Vasiliou V, McCord JM, et al. Activation of a novel isoform of methionine adenosyl transferase 2A and increased S-adenosylmethionine turnover in lung epithelial cells exposed to hyperoxia. Free Radic Biol Med. 2006;40(2):348–58.View ArticlePubMedGoogle Scholar
- Pajares MA, Alvarez L, Perez-Sala D. How are mammalian methionine adenosyltransferases regulated in the liver? A focus on redox stress. Febs Letters. 2013;587(12):1711–6.View ArticlePubMedGoogle Scholar
- Miller CL, Llenos IC, Dulay JR, Barillo MM, Yolken RH, Weis S. Expression of the kynurenine pathway enzyme tryptophan 2,3-dioxygenase is increased in the frontal cortex of individuals with schizophrenia. Neurobiol Dis. 2004;15(3):618–29.View ArticlePubMedGoogle Scholar
- Wu W, Nicolazzo JA, Wen L, Chung R, Stankovic R, Bao SSS, Lim CK, Brew BJ, Cullen KM, Guillemin GJ. Expression of Tryptophan 2,3-Dioxygenase and Production of Kynurenine Pathway Metabolites in Triple Transgenic Mice and Human Alzheimer's Disease Brain. Plos One. 2013;8(4):11.Google Scholar
- Wollmer MA. Cholesterol-related genes in Alzheimer's disease. Biochimica Biophys Acta Mol Cell Biol Lipids. 2010;1801(8):762–73.View ArticleGoogle Scholar
- Cui YJ, Huang MW, He YB, Zhang SY, Luo YZ. Genetic Ablation of Apolipoprotein A-IV Accelerates Alzheimer's Disease Pathogenesis in a Mouse Model. Am J Pathol. 2011;178(3):1298–308.View ArticlePubMedPubMed CentralGoogle Scholar
- Yegambaram M, Manivannan B, Beach TG, Halden RU. Role of Environmental Contaminants in the Etiology of Alzheimer's Disease: A Review. Curr Alzheimer Res. 2015;12(2):116–46.View ArticlePubMedGoogle Scholar
- Moon JC, Hah YS, Kim WY, Jung BG, Jang HH, Lee JR, Kim SY, Lee YM, Jeon MG, Kim CW, et al. Oxidative stress-dependent structural and functional switching of a human 2-Cys peroxiredoxin Isotype II that enhances HeLa cell resistance to H2O2-induced cell death. J Biol Chem. 2005;280(31):28775–84.View ArticlePubMedGoogle Scholar
- Larose C, Canuel R, Lucotte M, Di Giulio RT. Toxicological effects of methylmercury on walleye (Sander vitreus) and perch (Perca flavescens) from lakes of the boreal forest. Comp Biochem Physiol C Toxicol Pharmacol. 2008;147(2):139–49.View ArticlePubMedGoogle Scholar
- Mieiro CL, Ahmad I, Pereira ME, Duarte AC, Pacheco M. Antioxidant system breakdown in brain of feral golden grey mullet (Liza aurata) as an effect of mercury exposure. Ecotoxicology. 2010;19(6):1034–45.View ArticlePubMedGoogle Scholar
- Grune T, Reinheckel T, Davies KJA. Degradation of oxidized proteins in mammalian cells. Faseb Journal. 1997;11(7):526–34.PubMedGoogle Scholar
- Thamsen M, Kumsta C, Li F, Jakob U. Is Overoxidation of Peroxiredoxin Physiologically Significant? Antioxid Redox Signal. 2011;14(4):725–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanda H, Shinkai Y, Kumagai Y. S-Mercuration of cellular proteins by methylmercury and its toxicological implications. J Toxicol Sci. 2014;39(5):687–700.View ArticlePubMedGoogle Scholar
- Kanda H, Toyama T, Shinohara-Kanda A, Iwamatsu A, Shinkai Y, Kaji T, Kikushima M, Kumagai Y. S-Mercuration of rat sorbitol dehydrogenase by methylmercury causes its aggregation and the release of the zinc ion from the active site. Arch Toxicol. 2012;86(11):1693–702.View ArticlePubMedGoogle Scholar
- Martyniuk CJ, Fang B, Koomen JM, Gavin T, Zhang L, Barber DS, LoPachin RM. Molecular Mechanism of Glyceraldehyde-3-phosphate Dehydrogenase Inactivation by alpha, beta-Unsaturated Carbonyl Derivatives. Chem Res Toxicol. 2011;24(12):2302–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakajima H, Amano W, Kubo T, Fukuhara A, Ihara H, Azuma YT, Tajima H, Inui T, Sawa A, Takeuchi T. Glyceraldehyde-3-phosphate Dehydrogenase Aggregate Formation Participates in Oxidative Stress-induced Cell Death. J Biol Chem. 2009;284(49):34331–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Rutkowski DT, Kaufman RJ. A trip to the ER: coping with stress. Trends Cell Biol. 2004;14(1):20–8.View ArticlePubMedGoogle Scholar
- Berntssen MHG, Aatland A, Handy RD. Chronic dietary mercury exposure causes oxidative stress, brain lesions, and altered behaviour in Atlantic salmon (Salmo salar) parr. Aquat Toxicol. 2003;65(1):55–72.View ArticlePubMedGoogle Scholar
- Kushner JP, Barbuto AJ, Lee GR. An inherited enzymatic defect in porphyria cutanea tarda: decreased uroporphyrinogen decarboxylase activity. J Clin Investig. 1976;58(5):1089–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Woods JS, Bowers MA, Davis HA. Urinary porphyrin profiles as biomarkers of trace metal exposure and toxicity: studies on urinary porphyrin excretion patterns in rats during prolonged exposure to methyl mercury. Toxicol Appl Pharmacol. 1991;110(3):464–76.View ArticlePubMedGoogle Scholar
- Koss G, Schuler E, Arndt B, Seidel J, Seubert S, Seubert A. A comparative toxicological study on pike (Esox lucius L.) from the River Rhine and River Lahn. Aquat Toxicol. 1986;8(1):1–9.View ArticleGoogle Scholar
- Schlenk D, Handy R, Steinert S, Depledge MH, Benson W. Biomarkers. Boca Raton: Crc Press-Taylor & Francis Group; 2008.View ArticleGoogle Scholar
- Hahn ME, Chandran K. Uroporphyrin accumulation associated with cytochrome P4501A induction in fish hepatoma cells exposed to aryl hydrocarbon receptor agonists, including 2,3,7,8-tetrachlorodibenzo-p-dioxin and planar chlorobiphenyls. Arch Biochem Biophys. 1996;329(2):163–74.View ArticlePubMedGoogle Scholar
- Eide M, Karlsen OA, Kryvi H, Olsvik PA, Goksøyr A. Precision-cut liver slices of Atlantic cod (Gadus morhua): An in vitro system for studying the effects of environmental contaminants. Aquat Toxicol. 2014;153:110–5.View ArticlePubMedGoogle Scholar
- Anelli T, Sannino S, Sitia R. Proteostasis and "redoxtasis" in the secretory pathway: Tales of tails from ERp44 and immunoglobulins. Free Radic Biol Med. 2015;83:323–30.View ArticlePubMedGoogle Scholar
- Ung CY, Lam SH, Hlaing MM, Winata CL, Korzh S, Mathavan S, Gong ZY. Mercury-induced hepatotoxicity in zebrafish: in vivo mechanistic insights from transcriptome analysis, phenotype anchoring and targeted gene expression validation. BMC Genomics. 2010;11:14.View ArticleGoogle Scholar
- Ghazalpour A, Bennett B, Petyuk VA, Orozco L, Hagopian R, Mungrue IN, Farber CR, Sinsheimer J, Kang HM, Furlotte N, et al. Comparative Analysis of Proteome and Transcriptome Variation in Mouse. Plos Genetics. 2011;7(6):17.View ArticleGoogle Scholar
- Gust KA, Nanduri B, Rawat A, Wilbanks MS, Ang CY, Johnson DR, Pendarvis K, Chen XF, Quinn MJ, Johnson MS, et al. Systems toxicology identifies mechanistic impacts of 2-amino-4,6-dinitrotoluene (2A-DNT) exposure in Northern Bobwhite. BMC Genomics. 2015;16:17.View ArticleGoogle Scholar
- Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13(4):227–32.PubMedPubMed CentralGoogle Scholar
- Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29(5):661–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Aon MA, Stanley BA, Sivakumaran V, Kembro JM, O'Rourke B, Paolocci N, Cortassa S. Glutathione/thioredoxin systems modulate mitochondrial H2O2 emission: An experimental-computational study. J Gen Physiol. 2012;139(6):479–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Jo SH, Son MK, Koh HJ, Lee SM, Song IH, Kim YO, Lee YS, Jeong KS, Kim WB, Park JW, et al. Control of mitochondrial redox balance and cellular defense against oxidative damage by mitochondrial NADP(+)-dependent isocitrate dehydrogenase. J Biol Chem. 2001;276(19):16168–76.View ArticlePubMedGoogle Scholar
- Kil IS, Park JW. Regulation of mitochondrial NADP(+)-dependent isocitrate dehydrogenase activity by glutathionylation. J Biol Chem. 2005;280(11):10846–54.View ArticlePubMedGoogle Scholar
- Mailloux RJ, Jin XL, Willmore WG. Redox regulation of mitochondrial function with emphasis on cysteine oxidation reactions. Redox Biol. 2014;2:123–39.View ArticlePubMedGoogle Scholar
- Wells TNC, Saxty BA. Redox control of catalysis in ATP-citrate lysate from rat liver. Eur J Biochem. 1992;204(1):249–55.View ArticlePubMedGoogle Scholar
- Hildebrandt T, Knuesting J, Berndt C, Morgan B, Scheibe R. Cytosolic thiol switches regulating basic cellular functions: GAPDH as an information hub? Biol Chem. 2015;396(5):523–37.View ArticlePubMedGoogle Scholar
- Hurd TR, Prime TA, Harbour ME, Lilley KS, Murphy MP. Detection of reactive oxygen species-sensitive thiol proteins by redox difference gel electrophoresis - Implications for mitochondrial redox signaling. J Biol Chem. 2007;282(30):22040–51.View ArticlePubMedGoogle Scholar
- Hayes JD, Dinkova-Kostova AT. The Nrf2 regulatory network provides an interface between redox and intermediary metabolism. Trends Biochem Sci. 2014;39(4):199–218.View ArticlePubMedGoogle Scholar
- Cambier S, Benard G, Mesmer-Dudons N, Gonzalez P, Rossignol R, Brethes D, Bourdineaud JP. At environmental doses, dietary methylmercury inhibits mitochondrial energy metabolism in skeletal muscles of the zebra fish (Danio rerio). Int J Biochem Cell Biol. 2009;41(4):791–9.View ArticlePubMedGoogle Scholar
- Berle M, Kroksveen AC, Garberg H, Aarhus M, Haaland OA, Wester K, Ulvik RJ, Helland C, Berven F. Quantitative proteomics comparison of arachnoid cyst fluid and cerebrospinal fluid collected perioperatively from arachnoid cyst patients. Fluids Barriers CNS. 2013;10(1):17.View ArticlePubMedPubMed CentralGoogle Scholar
- Vaudel M, Barsnes H, Berven FS, Sickmann A, Martens L. SearchGUI: An open-source graphical user interface for simultaneous OMSSA and X!Tandem searches. Proteomics. 2011;11(5):996–9.View ArticlePubMedGoogle Scholar
- Vaudel M, Burkhart JM, Zahedi RP, Oveland E, Berven FS, Sickmann A, Martens L, Barsnes H. PeptideShaker enables reanalysis of MS-derived proteomics data sets. Nat Biotechnol. 2015;33(1):22–4.View ArticlePubMedGoogle Scholar
- Dysvik B, Jonassen I. J-Express: exploring gene expression data using Java. Bioinformatics. 2001;17(4):369–70.View ArticlePubMedGoogle Scholar
- Mann M, Kelleher NL. Precision proteomics: The case for high resolution and high mass accuracy. Proc Natl Acad Sci U S A. 2008;105(47):18132–8.View ArticlePubMedPubMed CentralGoogle Scholar
- MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010;26(7):966–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Habig WH, Pabst MJ, Jakoby WB. Glutathione S-transferases. The first enzymatic step in mercapturic acid formation. J Biol Chem. 1974;249(22):7130–9.PubMedGoogle Scholar
- Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pages F, Trajanoski Z, Galon J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(D1):D447–52.View ArticlePubMedGoogle Scholar
- Vizcaino JA, Csordas A, del-Toro N, Dianes JA, Griss J, Lavidas I, Mayer G, Perez-Riverol Y, Reisinger F, Ternent T, et al. 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 2016;44(D1):D447–56.View ArticlePubMedGoogle Scholar