- Research article
- Open Access
Liver proteome response of pre-harvest Atlantic salmon following exposure to elevated temperature
BMC Genomicsvolume 19, Article number: 133 (2018)
Atlantic salmon production in Tasmania (Southern Australia) occurs near the upper limits of the species thermal tolerance. Summer water temperatures can average over 19 °C over several weeks and have negative effects on performance and health. Liver tissue exerts important metabolic functions in thermal adaptation. With the aim of identifying mechanisms underlying liver plasticity in response to chronic elevated temperature in Atlantic salmon, label-free shotgun proteomics was used to explore quantitative protein changes after 43 days of exposure to elevated temperature.
A total of 276 proteins were differentially (adjusted p-value < 0.05) expressed between the control (15 °C) and elevated (21 °C) temperature treatments. As identified by Ingenuity Pathway Analysis (IPA), transcription and translation mechanisms, protein degradation via the proteasome, and cytoskeletal components were down-regulated at elevated temperature. In contrast, an up-regulated response was identified for NRF2-mediated oxidative stress, endoplasmic reticulum stress, and amino acid degradation. The proteome response was paralleled by reduced fish condition factor and hepato-somatic index at elevated temperature.
The present study provides new evidence of the interplay among different cellular machineries in a scenario of heat-induced energy deficit and oxidative stress, and refines present understanding of how Atlantic salmon cope with chronic exposure to temperature near the upper limits of thermal tolerance.
A major impact of climate change on fisheries and aquaculture is the increased seawater temperature she . Aquatic ectotherm species in the temperate areas are particularly vulnerable given the high warming rate of this area . In Tasmania (Southern Australia), seawater Atlantic salmon is produced near the upper limits of thermal tolerance . Post-smolt Atlantic salmon are generally considered to have an optimum temperature range of 13–15 °C for growth and an upper critical range of 22–33 °C. , while the sea surface temperature in Tasmania averages over 19 °C for several weeks during the summer period (IMOS-OceanCurrent, unpublished data) and further increases are projected . Sea-caged Atlantic salmon cannot escape the increased surface water temperatures by vertical migration and thus risk long-term exposure to heat stress and physiological challenges that both impair production efficiency and raise welfare issues .
Heat stress co-occurs with hypoxia, and although oxygen limited thermal tolerance is of major importance in the fish physiological response , the temperature-dependent response occurs independently of water oxygen saturation [8, 9]. Any temperature rise increases the metabolic rate and consequently the maintenance of energy requirements, leading to a state of metabolic remodeling to compensate for increased energy demand . This response is dependent on the species’ thermal tolerance (eurythermals vs. stenothermals) and the exposure regime (acute vs. short-term vs. long-term) . Individual tissues show highly divergent responses to thermal stress that are likely related to their physiological role in the body . The liver is the central hub for the regulation of nutrient metabolism and detoxification , and thus constitutes an excellent target to characterize mechanisms of acclimation to chronically elevated temperature. To our knowledge, the molecular plasticity of liver of Atlantic salmon exposed to chronically elevated temperature is currently limited to analysis at the transcriptome level . The most notable findings of this study were reduced protein synthesis and increased xenobiotic metabolism in fish held at 19 °C. Shotgun proteomics can provide an unprecedented view of the liver response of this species to chronic heat stress, and offers potential to unravel adaptive physiological mechanisms that can be only postulated by the transcriptome. In addition, biological interpretation of data generated from Atlantic salmon is greatly benefited from increased genomic resources such the recently published reference genome , associated protein database improvements (e.g. UniProt), and the continuous improvement of the bioinformatics tools used for gene ontology (GO) enrichment (e.g. Ingenuity Pathway Analysis, DAVID). In relation to previous fish (i.e. non-salmon species) proteomics research, the liver response to thermal stress has been only examined using gel-based approach [15, 16]. Shotgun proteomics can extend the range of quantifiable proteins and provide a more detailed characterization of the affected mechanisms that are common to different species and thermal stresses.
Here we explored the long-term exposure to elevated temperature on the molecular response in liver tissue of pre-harvest Tasmanian Atlantic salmon using label-free shotgun proteomics. We compared two water temperatures, 15 °C and 21 °C, under the same oxygen saturation levels and after 43 days of exposure. Quantitative changes in the protein expression pattern can contribute to the mechanistic understanding of how Atlantic salmon specifically, and fish generally, cope with chronically elevated temperatures. This study expands fundamental and applied information about fish nutrition under limiting environmental conditions previously conducted in Tasmania [17,18,19,20]. These findings may be also of value for the aquaculture industry in further development of dietary formulations for Atlantic salmon specific to the summer period.
Growth trial and sampling
The trial was conducted at the Experimental Aquaculture Facility (EAF) of the Institute for Marine and Antarctic Studies, University of Tasmania (Taroona, Tasmania, Australia) in accordance with University of Tasmania Animal Ethics (Investigation A0015208). All female Atlantic salmon post-smolt (average weight ± STD: 961 ± 172 g) from a single cohort were sourced from a commercial hatchery (Huon Aquaculture, Tasmania, Australia), haphazardly allocated amongst 6 × 2500 L circular tanks at an initial stocking density of 19 fish tank− 1 and acclimated for 38 days. Each tank had an independent recirculation system equipped with a heat-exchanger, protein skimmer, drum filter, UV filter and biological filter. Seawater was continually supplied and progressively replaced 1.5 times hour− 1 with 10% fresh water exchange day− 1. Water temperature was maintained at 15 °C during acclimation. Photoperiod was maintained at 12 h light:12 h dark. Water quality parameters (dissolved oxygen, pH, nitrate and nitrite) were recorded daily and maintained within limits for Atlantic salmon ; dissolved oxygen was maintained at 101 ± 0.8% and pH at 7.8 ± 0.0. Fish were fed a commercial diet (Optiline 8 mm, Skretting, Tasmania, Australia) in excess four times day− 1 with automatic feeders, and uneaten pellets were collected after termination of each meal to calculate daily feed intake. The feed contained 42% crude protein, 29% lipid, and 21.4 MJ kg− 1 digestible energy.
In order to examine the response to elevated chronic temperature, fish were exposed to two temperatures, 15 ± 0.0 °C and 21 ± 0.1 °C. Temperature in triplicate tanks was steadily increased (0.5 °C day− 1) over 13 days until it reached 21 °C, and then maintained for 43 days until the end of the experiment. At the beginning of the experiment, fish were anaesthetized (Aqui-S® 50 mg L− 1)  and benchmarked (wet weight, fork length and skin/fin condition). At the end of the experiment fish were euthanized (Aqui-S® 500 mg L− 1) and re-measured for wet weight and fork length. No signs of pathogenic infection to skin or gill were found in any of the treatments. Livers were dissected and weighed, immediately frozen in liquid nitrogen, and stored at − 80 °C for proteomic analysis. Growth rate was correlated with condition factor (Pearson’s r = 0.69, p < 0.01) and livers of fish with k < 1.2 (11% of fish at the end of the trial) were discarded from proteomic analysis in order to avoid confounding effects of poor fish condition on protein expression. The effect of elevated temperature on growth biometry and liver proteome was therefore assessed on an individual basis in three fish tank− 1 (nine treatment− 1) randomly selected among those with k > 1.2.
Liver preparation for proteomic analysis
Livers (~ 60 mg of frontal lobe tissue) from each treatment group (nine treatment− 1) were individually homogenized for 5 s in Eppendorf tubes containing lysis buffer (7 M urea, 2 M thiourea, 50 mM pH 8 Tris) and protease inhibitor cocktail (Roche, NSW, Australia) using Tissue-Tearor homogenator (Biospec Products, OK, USA). Each extraction was performed for 18–24 h at 4 °C with overnight rotation. After removal of insoluble material by centrifugation (13,000 rpm, 15 min at 4 °C), an aliquot was precipitated with 100% ethanol (9:1, v/v) overnight. Protein pellets were washed twice in 70% ethanol and re-suspended in lysis buffer. Protein concentrations were estimated with Bradford Protein Assay (Bio-Rad, NSW, Australia) using plate reader (Synergy TMHT, BioTek, QLD, Australia) and the volumes were adjusted with lysis buffer to achieve a concentration of 1 μg μL− 1 for each extract.
Nano-liquid chromatography and tandem mass spectrometry (LTQ-Orbitrap XL)
Protein samples were trypsin-digested using standard procedures  and analyzed by nanoLC-MS/MS using an LTQ-Orbitrap XL and Ultimate 3000 RSLCnano HPLC system (ThermoFisher Scientific, MA, USA). Tryptic peptides (~ 1 μg) were loaded onto a 20 mm × 75 μm PepMap 100 trapping column (3 μm C18) at 5 μl/min, using 98% water, 2% acetonitrile and 0.05% TFA. Peptides were separated at 0.3 μl/min on a 250 mm × 75 μm PepMap 100 RSLC column (2 μm C18) held at 40 °C, using a stepped gradient from 97% mobile phase A (0.1% formic acid in water) to 50% mobile phase B (0.08% formic acid in 80% acetonitrile and 20% water) comprising 3–10% B over 10 min, 10–40% B over 180 min, 40–50% B over 10 min, holding at 95% B for 10 min then re-equilibration in 3% B for 15 min. The LTQ-Orbitrap XL was controlled using Xcalibur 2.1 software in data-dependent mode and MS/MS spectra were acquired as described .
Database searching and criteria for protein identification
RAW files from the LTQ-Orbitrap were imported into MaxQuant software version 184.108.40.206 for peptide matching to MS/MS spectra and label-free protein quantification using the max LFQ algorithm . MS/MS spectra were searched against the Salmonidae database (http://uniprot.org/taxonomy/8030; 17,795 entries) using the Andromeda search engine. Default settings for protein identification were used, including a maximum of two missed cleavages, mass error tolerances of 20 ppm then 4.5 ppm for initial and main peptide searches, respectively, 0.5 Da tolerance for fragment ions, variable oxidation of methionine and fixed carbamidomethylation of cysteine. The false discovery rates (FDR) for peptide-spectrum matches and protein identification were both set to 0.01. MaxQuant output files of the complete peptide and protein-level mass spectrometry are provided in Additional files 1 and 2, respectively.
Calculations and statistical analysis
Standard formulae were used to assess growth biometrical data. Fulton’s condition factor was calculated as k = W / FL3, where W is fish wet weight (g) and FL is fork length (cm). Hepato-somatic index was determined as HSI = (LW / W) × 100, where LW is liver weight (g) and W is fish wet weight (g).
Statistical analyses of biometrical data was performed using R software . Individual fish data was analysed using the Generalized Estimating Equations (GEE) model to control the cluster (tank) correlation derived from the sampling of individuals from different tanks within each treatment . Tank and fish nested within tank were considered random variables. The Wald test was used to detect significant differences (p < 0.05) between treatments and results were expressed as mean ± standard error (SEM) (n = 9).
For statistical analysis of LTQ-Orbitrap mass spectrometry, the “ProteinGroups” output file generated by MaxQuant analysis of liver extracts was analysed in R  using the limma package . Proteins identified on the basis of a single matching peptide were excluded and only proteins detected in at least six out of nine biological replicates in any one treatment group were considered. The effect of temperature was investigated by fitting a linear model with log2 protein group intensity as the response and including tank and fish nested to tank as explanatory random variables. Prior to model fitting, intensity values were normalized using cyclic loess normalization  and the method of empirical array quality weights  was used to calculate sample reproducibility and down-weight less reproducible samples. After initial model fitting, empirical Bayes  was used to calculate moderated test statistics and Benjamini Hochberg correction was applied to adjust p-values for multiple testing. Missing values for all remaining proteins were excluded from the analysis with degrees of freedom adjusted accordingly. To gain further insight into the potential mechanisms of the effects of elevated temperature, differentially (adjusted p-value < 0.05) expressed proteins were selected for Ingenuity Pathway Analysis (IPA, https://qiagenbioinformatics.com/products/ingenuity-pathway-analysis/). The salmonidae genes were first mapped to human orthologues using PANTHER , and then gene symbols and the corresponding fold change in protein expression were submitted to IPA for identification of canonical and toxicity pathways and mapping of interaction networks.
Elevated temperature had a significant and negative effect on k and HSI, although this effect was not significantly reflected in final weight (Table 1). Feed intake (measured by tank) was not significantly affected by elevated temperature.
A total of 842 proteins were identified on the basis of two or more unique matching peptide sequences and presence in at least six of the nine biological replicates in either treatment group (Additional file 3). Peptide length ranged from seven to 43 residues and averaged 15 residues. Multidimensional scaling (MDS) showed maximum separation between temperature treatments (Figure 1a). On the basis of an adjusted p-value < 0.05, comparison of intensity values identified differences in 276 proteins, which are shown on the volcano plot in Fig. 1b. Of these proteins, 89 and 187 proteins were up-regulated and down-regulated in 21 °C livers relative to 15 °C livers, respectively, with fold changes ranging from 1.2 to 5.4. Differentially abundant proteins showing fold changes > 2 are listed in Table 2. One protein, ferritin (FTL), showed a significant (adjusted p-value < 0.05) degree of correlation (Pearson’s r = 0.92) with the measured HSI (Fig. 2) in addition to being significantly upregulated at the higher temperature.
Biological pathway analysis of regulated proteins using IPA software revealed 43 significant canonical pathways. The significance of the association between the data set and the pathway was determined based on the p-value, which determines the probability that the association between the data set and the pathway is explained only by chance, and on the ratio value, representing the number of proteins from the data found in each pathway over the total number of proteins in that pathway. The top five significant canonical pathways included “EIF2 Signaling” (p-value = 6.3 × 10− 30; ratio = 0.17), “Protein Ubiquitination Pathway” (p-value = 1.6 × 10− 19; ratio = 0.10), “Regulation of eIF4 and p70S6K Signaling” (p-value = 4.9 × 10− 10; ratio = 0.09), “mTOR Signaling” (p-value = 1.3 × 10− 8; ratio = 0.07) and “Mitochondrial Dysfunction” (p-value = 4.0 × 10− 6; ratio = 0.06), while those with the higher ratio included “Spliceosomal Cycle” (p-value = 2.0 × 10− 2; ratio = 0.5), “Fatty acid β-oxidation” (p-value = 1.5 × 10− 3; ratio = 0.33), “Pentose Phosphate Pathway” (p-value = 2.8 × 10− 4; ratio = 0.18), “β-alanine degradation” (p-value = 4.5 × 10− 3; ratio = 0.20), “Superoxide radicals degradation” (p-value = 4.5 × 10− 3; ratio = 0.20), “Endoplasmic Reticulum Stress Pathway” (p-value = 5.6 × 10− 5; ratio = 0.19) and “Valine Degradation” (p-value = 2.7 × 10− 5; ratio = 0.14). The activity pattern of the canonical pathway can be predicted by the activation z-score, a statistical measure of the match between expected relationship direction and observed expression . Predicted activity was found to be increased for “NRF2-mediated Oxidative Stress Response” (p-value = 1.7 × 10− 5; ratio = 0.06; z-score = 0.45). The full list of canonical pathways as determined by temperature-regulated proteins is shown in Additional file 4.
Toxicity pathways are canonical pathways that are significantly associated with toxicity lists. These are functional groupings based on critical biological processes and key toxicological responses, and describe adaptive, defensive, or reparative responses to xenobiotic insults. The significance of the association is also defined by a p-value and a ratio value. A total of 15 significant toxicity pathways were mined from the temperature-regulated proteins. The five most significant pathways also showed the highest ratios and included “Fatty acid metabolism” (−log p-value = 1.6 × 10− 25; ratio = 0.09), “Increases Transmembrane Potential of Mitochondria and Mitochondrial Membrane” (p-value = 4.0 × 10− 21; ratio = 0.14), “Mitochondrial Dysfunction” (p-value = 1.3 × 10− 19; ratio = 0.06), “Glutathione Depletion-CYP Induction and Reactive Metabolites” (p-value = 2.5 × 10− 18; ratio = 0.33) and “Oxidative Stress” (−log p-value = 5.0 × 10− 16; ratio = 0.11). The components of the top five toxicity pathways are shown in Additional file 5.
IPA software identified 15 top networks and ranked them by a score that considers the number of focus proteins and the size of the network to approximate the relevance of the network to the original list of focus proteins . Networks with scores ≥2 have at least a 99% confidence level of not being generated by random chance alone. Top five networks revealed links with “RNA Post-Transcriptional Modification, Cellular Assembly and Organization” (score 52), “Cell-To-Cell Signaling and Interaction”, “Cancer, Cell Death and Survival, Organismal Injury and Abnormalities” (score 44), “Protein Synthesis, Gene Expression, Developmental Disorder” (score 39), “Protein Trafficking, Molecular Transport, Cellular Compromise” (score 37), and “Lipid Metabolism, Small Molecule Biochemistry, Molecular Transport” (score 35). The components of the top four networks are shown in Fig. 3.
We have used shotgun proteomics to characterize changes in the liver proteome of Atlantic salmon following exposure to chronically elevated temperature (21 °C). To our knowledge, this is the first study to provide proteome-level evidence of such response in Atlantic salmon. The metabolic role and plasticity of liver tissue in the thermal adaptive process was reflected in the clear separation between temperature treatments by MDS and in the high number of proteins that were differentially regulated at 21 °C. The general pattern of down-regulation is in line with the liver transcriptome response of similar size Atlantic salmon upon chronically elevated temperature . The molecular dynamics of adaptation to thermal stress are dependent on tissue type examined, type of stress, and thermal exposure regime . Specifically in liver tissue, metabolic down-regulation and small magnitude response of stress indicators appears to be associated with chronic exposure to elevated temperature [13, 33], while an up-regulated response and larger fold changes are elicited in response to more acute thermal stresses [15, 16].
Physiological adjustments to maintain tissue homeostasis were paralleled by constraints in performance. Metabolic rate increases exponentially as temperature increases and, at any given temperature, the difference between feed intake and metabolic rate will determine the energy available for tissue development and fish growth . Given the similar feed intake between temperature treatments, the lower biometric indices (HSI, k) in the 21 °C fish suggested that the energy deficit caused by the increased energy demand was not met by feed intake. The proteomic data described here correlate with the biometric data and provide new insights into the main biological processes involved in the liver adaptive response of Atlantic salmon to chronically elevated temperature. We linked the pathways and networks generated by IPA analysis with expression patterns of protein groups and greatest fold changes to identify prevailing processes affected by elevated temperature and describe them in the following four sections.
Protein synthesis and degradation
Protein synthesis is a major energy consuming process that can account for up to 42% of total energy expenditure in fish . Of all tissues, the liver tissue has one of the highest rates of protein synthesis , and so is considered a highly sensitive physiological indicator of the impact of elevated temperature [35, 36]. Suppression of protein synthesis observed in the current study is consistent with previous liver transcriptional data in Atlantic salmon , and indicates a possible compensatory response to the energy deficit induced under chronically elevated sub-lethal temperature. In both studies, two independent mechanisms: eukaryotic initiation factor 2 (EIF2) signaling and the mechanistic target of the rapamycin (mTOR) signaling pathway were implicated. These regulatory mechanisms respond to energy status and stress within the cell via the activity of serine/threonine kinases and phosphatases through phosphorylation/dephosphorylation of EIF2α  and downstream effectors such as eiF4 and EEF2 . Accordingly, reduced expression levels of ribosomal subunits (25 proteins), translation initiation (eight proteins) and elongation factors (five proteins), including EIF2, eiF4 and EEF2 subunits, were accompanied by reduced serine/threonine phosphatase activity (PP1CC and PPM1G).
Global repression of the translational machinery was also reflected by the inhibition of pre-translational regulatory mechanisms, and by inhibition of cytoplasmic chaperones that assist in the folding and trafficking of nascent proteins. Down-regulation of heterogeneous nuclear ribonucleoproteins (HNRNPK, HNRNPH1, HNRNPAB, HNRNPA0, SYNCRIP, PCBP2), responsible for packing and stabilizing freshly transcribed pre-mRNA , mirrored a reduced demand for mRNA export out of the nucleus to translation active sites in the cytoplasm. Noteworthy was the down-regulation of chromatin regulators, including high mobility group proteins (HMGB2 and HMGB3), the intracellular hyaluran-binding protein (HABP4) and nucleophosmin (NMP1). Chromatin remodelling appears to be a critical process in compensating for elevated temperature effects , as also suggested in the liver transcriptome analysis of heat-stressed eurythermal annual killifish (Austrofundulus limnaeus) . In cytoplasm, reduced chaperone demand was indicated by the 2.5-fold down-regulation of the heat shock cognate 71 kDa protein (HSPA8), which contributes to the overall cytoplasmic folding by binding approximately 20% of newly translated proteins . The chaperonin-containing TCP1 complex (TCP1, CCT2, CCT3, CCT4, CCT8) was also down-regulated. While the TCP1 complex is involved in the folding of only ~ 1% of newly synthesized proteins , these mostly include microtubule proteins, which were also regulated as a consequence of thermal stress, and are discussed below in the context of cytoskeletal integrity.
An important finding of this study was the inhibitory effect of chronically elevated temperature on protein degradation via the ATP-dependent ubiquitin-proteasome pathway. Proteins involved in ubiquitin conjugation (UBE2V1, UBE2D2, UBQLN4), ubiquitin ligase activation (SKP, TCEB2) and de-ubiquitination (UCLH5, ADRM1), and in the subsequent proteasome degradation (five PSMCs, five PSMDs, five PSMAs and one PSME), were down-regulated. This was additionally paralleled by the down-regulation of cathepsins (CTSD, CATM) involved in lysosomal protein degradation pathway. Specifically, CTSD has been previously measured as an indicator of protein degradation in response to temperature stress in fish [20, 44]. Signs of suppressed protein degradation were however not detected in the liver transcriptome of similar size Atlantic salmon exposed to chronically elevated temperature, in that case at 19 °C . This observation corroborates the impact that the extra energy deficit generated by an additional rise of 2 °C has on protein metabolism; under restricted feeding, such a difference in rainbow trout reduced the liver rates of protein synthesis and degradation by ~ 35–55% . Reduced protein degradation is therefore a compensatory mechanism that follows protein synthesis inhibition in order to maintain growth and conserve energy under a temperature-induced energy deficit not compensated by feed intake. This concept is well established at temperatures within the upper critical range , however, it also seems to apply to the chronic exposure to temperatures considered to be within the range of optimal thermal tolerance, albeit at the upper limit.
Energy and lipid metabolism
The liver is a central hub for the storage and conversion of high-energy substrates. Under optimal temperature and feeding conditions, fish use amino acids rather than glucose as preferential energy source, while the contribution of fatty acid oxidation to energy production is correlated to dietary lipid levels . At elevated temperature, energy metabolism is remodeled to compensate for the consequences of increased metabolic rate [10, 13]. Consistent with the trend observed in the plasma metabolome of Atlantic salmon , our proteomic data suggested an increased dependence on amino acids rather than glucose and fatty acids for energy production upon chronically elevated temperature. This was reflected in the up-regulation of mitochondrial enzymes involved in the degradation of valine (HIBADH, HIBC), tryptophan (AFMID, GCDH) and leucine (MCCC1, MCCC2, IVD) towards the formation of citric acid (TCA) cycle intermediates, and supported by the increased expression of TCA enzymes such as glutamate dehydrogenase (GLUD1) and aspartate aminotransferase (GOT2). The apparent increased TCA flux from amino acid catabolism was paralleled by signs of reduced glucose availability, including up-regulation of glycogen catabolic enzymes (PYGB, MTAP), down-regulation of several proteins involved in the pentose phosphate pathway (PGD, RPE, TKT, TALDO1), and the down-regulation of pyruvate kinase (KPYK). Mitochondrial β-oxidation of fatty acids also showed signs of suppression (ACAD11, ECI2). This may be linked to the fact that, under unrestricted feeding, TCA dependence on fatty acid oxidation diminishes with increasing exposure time to elevated temperature .
Increasing temperature increases the fluidity of cell membrane, leading to alterations in lipid metabolism to maintain and stabilize fluidity [3, 48]. In this study, several proteins associated with fatty acid metabolism and lipid transport were down-regulated at elevated temperature and linked to mechanisms involved in modulating membrane fluidity. The peroxisomal acyl-coenzyme A oxidase 3 (ACOX3), the rate-limiting enzyme in the oxidative breakdown of methyl-branched fatty acids , showed the largest fold change in abundance (× 5.4). Since methyl-branched fatty acids have been suggested to enhance the fluidity of the membrane lipid bilayers over different environmental conditions , ACOX3 down-regulation may be part of a mechanism to increase membrane fluidity at elevated temperature. Membrane fluidity is stabilized by the incorporation of cholesterol into lipid bilayers, with increased cholesterol levels associated to increased temperatures . An important mechanism in cholesterol homeostasis is reverse cholesterol transport (RCT), whereby cholesterol is transported in high-density lipoproteins (HDL) from peripheral tissues back to the liver to be further eliminated in the bile . Elevated temperature appeared to inhibit RCT, as reflected by the parallel down-regulation of apolipoprotein A1 (APOA1), the major component of HDL, and lipase (LIPA), involved in HDL uptake and a well-known rate-limiting enzyme in RCT . Subsequent signs of altered cholesterol metabolism included down-regulation of epididymal secretory protein E1 (NCP2), an intracellular cholesterol transporter that regulates cholesterol biliary secretion, and of 3-oxo-5-beta-steroid 4-dehydrogenase (AKR1D1), involved in cholesterol breakdown towards the synthesis of bile acids. These observations re-emphasize the importance of cholesterol in the acclimatization to elevated temperature , and further implicate suppression of RCT and cholesterol catabolism to increase peripheral retention under conditions of thermal stress.
Stress-induced perturbation of transcription and translational mechanisms can block cell growth and proliferation via effects on cytoskeletal integrity . Thus, the microtubule network of α- and β-tubulins forming the cell cytoskeleton is recognized to play a role in maintenance of cell homeostasis and execution of a variety of cell stress responses . In line with a previous observation in the liver transcriptome of Antarctic fish , the liver proteome of Atlantic salmon reflected suppression of the microtubule dynamics in response to chronically elevated temperature. Down-regulation of tubulins (TUBB, TUBB1) was paralleled by down-regulation of the chaperonin-containing TCP1 complex (TCP1, CCT2, CCT3, CCT4, CCT8), which is particularly implicated in the folding and assembly of tubulins in an ATP-dependent manner . Notably, microtubule-stabilizing proteins were also down-regulated, including tubulin-folding factor (TBCB) , translationally-controlled tumor protein (TPT1) , N-myc downstream regulated gene 1 (NDRG1) , and stathmin (STMN1) . We highlight the larger than two-fold changes in TPT1 and NDRG1, which are established hallmarks and mediators of cell proliferation through microtubule stabilization. A novel outcome of the heat-induced cytoskeleton remodeling was the striking (4.4-fold) down-regulation of the canopy homologue 2 (CNPY2). CNPY2 is a transmembrane protein that regulates myosin regulatory light (MRLC) protein levels, a protein that links cytoskeleton to membrane proteins and stimulates cell growth . Reduced cytoskeletal integrity at elevated temperature has been attributed to maintenance of a less dense subcellular structure resulting from alterations in cytosol solubility and viscosity , or to an increase in the pool of soluble tubulin as a result of oxidative stress-induced microtubule depolymerization . In contrast to our results, acute heat stress induced up-regulation of genes encoding cytoskeleton components in gill tissue of Pacific salmon . This difference between studies highlights the tissue- and exposure regime- specificity of the heat stress response across different protein groups and biological processes .
Oxidative stress and endoplasmic reticulum (ER) stress
Exposing ectotherms to elevated temperatures challenges the cellular redox balance and leads to increased production of reactive oxygen species (ROS) . Mitochondria are well-known as major sites of ROS production and consumption and are implicated as the main source of thermally-induced oxidative stress . ROS imbalance occurs due to the uncoupling of mitochondrial respiration; a proton gradient generated by complex I, III and IV, but not coupled with consumption by complex V, will leak, consequently increasing the mitochondrial membrane potential and ROS formation [64, 65]. This mechanism, previously proposed as a cause of ROS formation in polar and temperate fish upon long-term warm acclimation , is consistent with our findings here in Atlantic salmon, with imbalances in the electron transport complex IV (COXAI2, COX6B1, COX7A2L, MT-CO2) and complex V (ATPL5L) proteins and in other mitochondrial proteins (CAT, SOD2, CASP3, MSRB2) indicative of mitochondrial dysfunction and increased mitochondrial transmembrane potential. The ROS-induced cellular antioxidant response was predicted to be mediated via nuclear factor erythroid 2-related factor 2 (Nrf2) through the regulation of different mechanisms, which as previously described  included induction of catabolism of the ROS superoxide through mitochondrial superoxide dismutase (SOD2) and the peroxiredoxin system, and metal chelation by ferritin (FTL). The members of the peroxiredoxin family exhibited opposing patterns of regulation (e.g. reduced PRDX 1 and PRDX6 and increased PRDX3 and PRDX4) that are characteristic of other disorders related to oxidative stress  and attributed to the increased exposure to the superoxide derivative hydrogen peroxide [69, 70]. Specifically, the observed co-regulation of PRDX3 and PRDX4 has been proposed as an outcome of antioxidant response and oxidative damage in cancer . The 3-fold up-regulation of FTL was similar to the increased transcript and protein abundance in Channa liver upon chronically elevated temperature . Iron interacts with hydrogen peroxide via the Fenton reaction leading to the production of hydroxyl radical ROS, a very reactive initiator of lipid peroxidation . As FTL maintains iron homeostasis, overexpression is proposed as an attempt to prevent Fenton type reactions and ROS accumulation . The Nrf2-mediated antioxidant response was further supported by up-regulation of metallothionein B (mtb), a cysteine-rich metal binding protein involved in metal homeostasis and ROS scavenging . Induction of the ROS scavenging system was also paralleled by the up-regulation of cytosolic epoxide hydrolases (EPHX1, EPHX2) and aldo-keto reductases (AKR7A2, AKR1B1), respectively involved in the detoxification of fatty acid epoxides and aldehydes [75, 76], and suggesting some of degree of lipid peroxidation in 21 °C livers, as previously detected in heat-stressed fish [77, 78]. A contribution of our study towards biomarker discovery was the significant and negative correlation between FTL expression and HSI across fish of both 15 °C and 21 °C groups (Figure 2). FTL is used as a reliable blood marker of liver disorders in human medicine , and it would be interesting to further explore the FTL expression levels in blood as a surrogate and less invasive measurement of liver condition and oxidative stress in fish.
Oxidative stress is a condition of imbalance between the formation of ROS and the biological system’s ability to detoxify the reactive intermediates, thus dysfunction of the antioxidant and detoxifying activity is also reflective of oxidative damage . This concept is consistent with previous thermal stress studies in fish examining varying exposure regimes [11, 62] and was supported here by the down-regulation of catalase (CAT), glutathione transferases (GSTs: GSTT1 and GSTP1), and betaine-homocysteine S-methyltransferase 1 (BHMT) at elevated temperature. CAT is a mitochondrial and peroxisomal consumer of hydrogen peroxide, though down-regulation in fish has been attributed to its sensitivity to the fluctuation of superoxide radicals . Glutathione plays a central role in the cellular defence against lipid peroxidation . While GSTs catalyse the metabolism of lipid peroxides by conjugation with glutathione and NADPH , BHMT maintains steady levels of the glutathione precursor s-adenosyl-methionine (SAM) (i.e. SAM is the amino acid methionine bound to an ATP molecule) . Co- regulation of GSTs and BHMT was thus indicative of glutathione depletion and reduced hepato-protection, as previously detected and associated to increased lipid peroxidation in liver of fish exposed to thermal stress [15, 77, 78]. These results collectively target the methionine cycle as a possible route to mitigate the impaired enzymatic hepato-protection associated to thermal stress. Enhanced protection was recently suggested in seabream fed a winter-specific diet supplemented with methionine . A comparison of our proteomic data with the liver transcriptome of similar size Atlantic salmon upon chronic exposure to 19 °C  implicates variation in the oxidative stress response that attributes to temperature-specific effects and the consequent ROS production [11, 63]. This highlights the risk of drawing general conclusions and also the importance of other experimental factors such as dietary formulation (e.g. antioxidant supplementation) and methodology used (e.g. proteomics vs. transcriptomics, pooling vs. individual-based analysis) in identifying thermal stressors across studies.
There is accumulating evidence for an intrinsic link between cellular oxidative stress and ER stress . One proposed mechanism for this interrelation is that mitochondrial ROS formation promotes the calcium release from the ER, leading to the accumulation of unfolded proteins in the ER lumen and further contributing to ROS formation due to the excessive calcium influx into the mitochondria . The ER unfolded protein response (UPR) is subsequently activated and mediated through the induction of molecular chaperones to restore proteostasis and avoid apoptosis. At elevated temperature, activation of the UPR response was indicated by the up-regulation of chaperones and other proteins required for protein folding and stabilization of pre- and un-folded proteins, specifically calreticulin (CALR), disulfide isomerases (PDIA3, PDIA6), 78 kDa glucose-regulated protein (HSPA5/Bip or GRP78) and ER glycosyltransferases (UGT2A1, RPN1). Our data confirm previous detections of ER stress and UPR in liver of thermal-stressed fish [15, 16, 33, 84] and demonstrate its specific importance in the adaptive response to chronically elevated temperature. Another indicator of the crosstalk between cellular oxidative stress and ER stress was the 3-fold up-regulation of the collagen-specific chaperone Serpin H1 (SERPINH1 or HSP47) which is also induced by stress induced-lipid peroxidation [87, 88]. In line with our study, increased transcript levels of SERPINH1 were also found in gill tissue of Pacific salmon under different heat exposure regimes [61, 89, 90].
This study increases our understanding of the molecular mechanisms occurring in the liver of pre-harvest Atlantic salmon that are important for coping with chronically elevated temperature, which is of increasing importance in temperate production areas such as Tasmania. Suppression of protein synthesis and degradation pathways appears to be the main energy saving mechanism for the increased metabolic demand, which is also reflected in an increased dependence of amino acid catabolism towards energy production. Other chronic heat-stress related mechanisms included the reverse transport and catabolism of cholesterol, cytoskeletal dynamics, the Nrf2-mediated antioxidant response and the endoplasmic reticulum UPR. Many of the proteins regulated here are common to the thermal stress response across fish species and tissues. These included protein groups (TCP1 complex proteins, high mobility group proteins, disulfide isomerases, peroxiredoxins and glutathione transferases) and individual proteins (CTSD, SOD2, FTL, BHMT, HSPA5, CALR and SERPINH1). Other stressors, particularly several cytoskeletal-related proteins (STMN1, TPT1, NDRG1 and CNPY2) were for first time here reported in response to thermal stress. Finally, opportunities for new research towards the development of salmon feed formulations specific to the summer period are raised here. Further understanding of the potential for methionine supplementation to improve the liver detoxifying capacity, or the evaluation of ingredients and additives to compensate the energy deficit and to spare amino acid degradation towards energy production are warranted.
De Silva SS, Soto D. Climate change and aquaculture: potential impacts, adaptation and mitigation. In: Cochrane K, De Young C, Soto D, Bahr T, editors. Climate change implications for fisheries and aquaculture: overview of current scientific knowledge, vol. 530. Rome: FAO Fisheries and Aquaculture; 2009. p. 151–212.
Neuheimer AB, Thresher RE, Lyle JM, Semmens JM. Tolerance limit for fish growth exceeded by warming waters. Nat Clim Chang. 2011;1(2):110–3.
Miller MR, Nichols PD, Barnes J, Davies NW, Peacock EJ, Carter CG. Regiospecificity profiles of storage and membrane lipids from the gill and muscle tissue of Atlantic salmon (Salmo salar L.) grown at elevated temperature. Lipids. 2006;41(9):865–76.
Jonsson B, Jonsson N. A review of the likely effects of climate change on anadromous Atlantic salmon Salmo salar and brown trout Salmo trutta, with particular reference to water temperature and flow. J Fish Biol. 2009;75(10):2381–447.
Oliver ECJ, Wotherspoon SJ, Chamberlain MA, Holbrook NJ. Projected Tasman Sea extremes in sea surface temperature through the twenty-first century. J Clim. 2014;27(5):1980–98.
Battaglene S, Carter C, Hobday A, Lyne V, Nowak BF. Scoping study into adaptation of the Tasmanian salmonid aquaculture industry to potential impacts of climate change. Hobart, Australia: Tasmanian Aquaculture and Fisheries Industry, University of Tasmania; 2008.
Portner HO, Knust R. Climate change affects marine fishes through the oxygen limitation of thermal tolerance. Science. 2007;315(5808):95–7.
Kullgren A, Jutfelt F, Fontanillas R, Sundell K, Samuelsson L, Wiklander K, Kling P, Koppe W, Larsson DG, Bjornsson BT, et al. The impact of temperature on the metabolome and endocrine metabolic signals in Atlantic salmon (Salmo salar). Comp Biochem Physiol A Mol Integr Physiol. 2013;164(1):44–53.
Hevroy EM, Waagbo R, Torstensen BE, Takle H, Stubhaug I, Jorgensen SM, Torgersen T, Tvenning L, Susort S, Breck O, et al. Ghrelin is involved in voluntary anorexia in Atlantic salmon raised at elevated sea temperatures. Gen Comp Endocrinol. 2012;175(1):118–34.
Jobling M. In: Jobling M, editor. Part 2: physiological energetics: feeding, metabolism and growth. In: Fish bioenergetics. London, UK: Chapman and Hall; 1994.
Logan CA, Buckley BA. Transcriptomic responses to environmental temperature in eurythermal and stenothermal fishes. J Exp Biol. 2015;218(Pt 12):1915–24.
Johnson LL, Anulacion BA, Arkoosh MR, Burrows DG, Da Silva DAM, Dietrich JP, Myers MS, Spormberg J, Ylitalo GM. Effects of legacy persistent organic pollutants (POS) in fish - current and future challenges. In: Tierney K, Farrell A, Brauner CJ, editors. Fish physiology: organic chemical toxicology of fishes, vol. 33. San Diego, USA: Elsevier Science; 2013.
Olsvik PA, Vikeså V, Lie KK, Hevrøy EM. Transcriptional responses to temperature and low oxygen stress in Atlantic salmon studied with next-generation sequencing technology. BMC Genomics. 2013;14(1):1–21.
Lien S, Koop BF, Sandve SR, Miller JR, Kent MP, Nome T, Hvidsten TR, Leong JS, Minkley DR, Zimin A, et al. The Atlantic salmon genome provides insights into rediploidization. Nature. 2016;533(7602):200–5.
Ibarz A, Martín-Pérez M, Blasco J, Bellido D, de Oliveira E, Fernández-Borràs J. Gilthead sea bream liver proteome altered at low temperatures by oxidative stress. Proteomics. 2010;10:963–75.
Mahanty A, Purohit GK, Banerjee S, Karunakaran D, Mohanty S, Mohanty BP. Proteomic changes in the liver of Channa striatus in response to high temperature stress. Electrophoresis. 2016;37(12):1704–17.
Katersky RS, Carter CG. A preliminary study on growth and protein synthesis of juvenile barramundi, Lates calcarifer at different temperatures. Aquaculture. 2007;267(1–4):157–64.
Katersky RS, Carter C. Growth efficiency of juvenile barramundi, Lates calcarifer (Bloch), at high temperatures. Aquaculture. 2005;250:775–80.
Katersky RS, Carter CG. High growth efficiency occurs over a wide temperature range for juvenile barramundi Lates calcarifer fed a balanced diet. Aquaculture. 2007;272(1–4):444–50.
Carter C, Katersky RS, Barnes J, Bridle A, Hauler RC. Assessment of fish growth peformance under limiting evironmental conditions: aquaculture nutrition subprogram. In., vol. project no. 2004/237. Deakin, Australia: Fisheries Research and Development Corporation; 2008.
Wedemeyer GA. Physiology of fish in intensive culture systems. In: Wedemeyer GA, editor. Interactions with water quality conditions. New York: Chapman and Hall; 1996.
Javahery S, Nekoubin H, Moradlu AH. Effect of anaesthesia with clove oil in fish (review). Fish Physiol Biochem. 2012;38:1545–52.
Wilson R, Golub SB, Rowley L, Angelucci C, Karpievitch YV, Bateman JF, Fosang AJ. Novel elements of the chondrocyte stress response identified using an in vitro model of mouse cartilage degradation. J Proteome Res. 2016;4(15):1033–50.
Cox J, Hein MY, Luber CA, Paron I, Nagaraj N, Mann M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol Cell Proteomics. 2014;13(9):2513–26.
Core Team R. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016.
Højsgaard S, Halekoh U, Yan J. The R package geepack for generalized estimating equations. 2005;15(2):1–10.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19(2):185–93.
Ritchie ME, Diyagama D, Neilson J, van Laar R, Dobrovic A, Holloway A, Smyth GK. Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006;7:261.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.
Mi H, Muruganujan A, Thomas PD. PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 2013;41(Database issue):D377–86.
Krämer A, Green J, Pollard J, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523–30.
Mininni AN, Milan M, Ferraresso S, Petochi T, Di Marco P, Marino G, Livi S, Romualdi C, Bargelloni L, Patarnello T. Liver transcriptome analysis in gilthead sea bream upon exposure to low temperature. BMC Genomics. 2014;15(1):1–12.
Carter C, Houlihan DF. Protein synthesis. In: Fish physiology: nitrogen excretion, vol. 20. London: Wright and Andersen; 2001. p. 31–75.
Morgan IJ, McDonald DG, Wood CM. The cost of living for freshwater fish in a warmer, more polluted world. Glob Chang Biol. 2001;7:345–55.
Katersky RS, Carter CG. The effect of temperature on post-prandial protein synthesis in juvenile barramundi, Lates calcarifer. Comp Biochem Physiol A Mol Integr Physiol. 2010;156(4):529–36.
Stolboushkina EA, Garber MB. Eukaryotic type translation initiation factor 2: structure–functional aspects. Biochem Mosc. 2011;76(3):283–94.
Showkat M, Beigh MA, Andrabi KI. mTOR signaling in protein translation regulation: implications in cancer genesis and therapeutic interventions. Mol Biol Int. 2014;2014:14.
Geuens T, Bouhy D, Timmerman V. The hnRNP family: insights into their role in health and disease. Hum Genet. 2016;135(8):851–67.
Portner HO, Peck LS, Somero G. Mechanisms defining thermal limits and adpatation in marine ectotherms: an integrate view. In: Rogers AD, Nadine MJ, Murphy EJ, West Sussex CA, editors. Antarctic ecosystems an extreme environment in a changing world. UK: Blackwell Publishing Ltd; 2012.
Podrabsky JE, Somero GN. Changes in gene expression associated with acclimation to constant temperatures and fluctuating daily temperatures in an annual killifish Austrofundulus Limnaeus. J Exp Biol. 2004;207(Pt 13):2237–54.
Stricher F, Macri C, Ruff M, Muller S. HSPA8/HSC70 chaperone protein: structure, function, and chemical targeting. Autophagy. 2013;9(12):1937–54.
Dekker C, Stirling PC, McCormack EA, Filmore H, Paul A, Brost RL, Costanzo M, Boone C, Leroux MR, Willison KR. The interaction network of the chaperonin CCT. EMBO J. 2008;27(13):1827–39.
Mommsen TP. Salmon spawning migration and muscle protein metabolism: the august Krogh principle at work. Comp Biochem Physiol B Biochem Mol Biol. 2004;139(3):383–400.
Morgan IJ, D'Cruz LM, Dockray JJ, Linton TK, Wood CM. The effects of elevated summer temperature and sublethal pollutants (ammonia, low pH) on protein turnover in the gill and liver of rainbow trout (Oncorhynchus mykiss) on a limited food ration. Comp Biochem Physiol A Mol Integr Physiol. 1999;123(1):43–53.
Reid SD, Dockray JJ, Linton TK, McDonald DG, Wood CM. Effects of chronic environmental acidification and a summer global warming scenario: protein synthesis in juvenile rainbow trout (Oncorhynchus mykiss). Can J Fish Aquat Sci. 1997;54:2014–24.
Médale F, Guillaume J. Nutritional energetics. In: Guillaume J, Kaushik S, Bergot P, Chichester MR, editors. Nutrition and feeding of fish and crustaceans. UK: Springer; 2001. p. 59–79.
Los DA, Murata N. Membrane fluidity and its roles in the perception of environmental signals. Biochim Biophys Acta. 2004;1666(1–2):142–57.
Hunt MC, Siponen MI, Alexson SEH. The emerging role of acyl-CoA thioesterases and acyltransferases in regulating peroxisomal lipid metabolism. Biochim Biophys Acta Mol basis Dis. 2012;1822(9):1397–410.
Poger D, Caron B, Mark AE. Effect of methyl-branched fatty acids on the structure of lipid Bilayers. J Phys Chem B. 2014;118(48):13838–48.
Robertson JC, Hazel JR. Membrane constraints to physiological function at different temperatures: does cholesterol stabilize membranes at elevated temperatures? In: Wood CM, DG MD, editors. Global warming: implications for freshwater and marine fish, vol. 61. Cambridge: Cambridge University Press; 1997.
Ghosh S. Early steps in reverse cholesterol transport: cholesteryl ester hydrolase and other hydrolases. Curr Opin Endocrinol Diabetes Obes. 2012;19(2):136–41.
Kim S, Coulombe PA. Emerging role for the cytoskeleton as an organizer and regulator of translation. Nat Rev Mol Cell Biol. 2010;11(1):75–81.
Parker AL, Kavallaris M, McCarroll JA. Microtubules and their role in cellular stress in cancer. Front Oncol. 2014;4
Windisch HS, Frickenhaus S, John U, Knust R, Pörtner HO, Lucassen M. Stress response or beneficial temperature acclimation: transcriptomic signatures in Antarctic fish (Pachycara Brachycephalum). Mol Ecol. 2014;23(14):3469–82.
Szymanski D. Tubulin folding cofactors: half a dozen for a Dimer. Curr Biol. 2002;12(22):R767–9.
Chan THM, Chen L, Guan X-Y. Role of Translationally controlled tumor protein in cancer progression. Biochem Res Int. 2012;2012:5.
Melotte V, Qu X, Ongenaert M, van Criekinge W, de Bruine AP, Baldwin HS, van Engeland M. The N-myc downstream regulated gene (NDRG) family: diverse functions, multiple applications. FASEB J. 2010;24(11):4153–66.
Rubin CI, Atweh GF. The role of stathmin in the regulation of the cell cycle. J Cell Biochem. 2004;93(2):242–50.
Hatta K, Guo J, Ludke A, Dhingra S, Singh K, Huang ML, Weisel RD, Li RK. Expression of CNPY2 in mouse tissues: quantification and localization. PLoS One. 2014;9(11):e111370.
Tomalty KM, Meek MH, Stephens MR, Rincon G, Fangue NA, May BP, Baerwald MR. Transcriptional response to acute thermal exposure in juvenile chinook salmon determined by RNAseq. G3 (Bethesda, Md). 2015;5(7):1335–49.
Tomanek L. Proteomic responses to environmentally induced oxidative stress. J Exp Biol. 2015;218(12):1867–79.
Banh S, Wiens L, Sotiri E, Treberg JR. Mitochondrial reactive oxygen species production by fish muscle mitochondria: potential role in acute heat-induced oxidative stress. Comp Biochem Physiol B Biochem Mol Biol. 2016;191:99–107.
Rousset S, Alves-Guerra MC, Mozo J, Miroux B, Cassard-Doulcier AM, Bouillaud F, Ricquier D. The biology of mitochondrial uncoupling proteins. Diabetes. 2004;53(Suppl 1):S130–5.
Srinivasan S, Avadhani NG. Cytochrome c oxidase dysfunction in oxidative stress. Free Radic Biol Med. 2012;53(6):1252–63.
Mark FC, Lucassen M, Portner HO. Thermal sensitivity of uncoupling protein expression in polar and temperate fish. Comp Biochem Physiol D Genomics Proteomics. 2006;1(3):365–74.
Ma Q. Role of nrf2 in oxidative stress and toxicity. Annu Rev Pharmacol Toxicol. 2013;53:401–26.
Pagano G, Aiello Talamanca A, Castello G, Cordero MD, Ischia M, Gadaleta MN, Pallard FV, Petrovi S, Tiano L, Zatterale A. Oxidative stress and mitochondrial dysfunction across broad-ranging pathologies: toward mitochondria-targeted clinical strategies. Oxidative Med Cell Longev. 2014;2014:27.
Rhee SG, Kang SW, Jeong W, Chang TS, Yang KS, Woo HA. Intracellular messenger function of hydrogen peroxide and its regulation by peroxiredoxins. Curr Opin Cell Biol. 2005;17(2):183–9.
Cox AG, Pullar JM, Hughes G, Ledgerwood EC, Hampton MB. Oxidation of mitochondrial peroxiredoxin 3 during the initiation of receptor-mediated apoptosis. Free Radic Biol Med. 2008;44(6):1001–9.
Basu A, Banerjee H, Rojas H, Martinez SR, Roy S, Jia Z, Lilly MB, De Leon M, Casiano CA. Differential expression of peroxiredoxins in prostate cancer: consistent upregulation of PRDX3 and PRDX4. Prostate. 2011;71(7):755–65.
Chen X, Guo C, Kong J. Oxidative stress in neurodegenerative diseases. Neural Regener Res. 2012;7(5):376–85.
Orino K, Lehman L, Tsuji Y, Ayaki H, Torti SV, Torti FM. Ferritin and the response to oxidative stress. Biochem J. 2001;357(Pt 1):241–7.
Ruttkay-Nedecky B, Nejdl L, Gumulec J, Zitka O, Masarik M, Eckschlager T, Stiborova M, Adam V, Kizek R. The role of metallothionein in oxidative stress. Int J Mol Sci. 2013;14(3):6044–66.
Fretland AJ, Omiecinski CJ. Epoxide hydrolases: biochemistry and molecular biology. Chem Biol Interact. 2000;129(1–2):41–59.
Penning TM. The aldo-keto reductases (AKRs): overview. Chem Biol Interact. 2015;234:236–46.
Lushchak VI, Bagnyukova TV. Temperature increase results in oxidative stress in goldfish tissues. 1. Indices of oxidative stress. Comp Biochem Physiol C Toxicol Pharmacol. 2006;143(1):30–5.
Kaur M, Atif F, Ali M, Rehman H, Raisuddin S. Heat stress-induced alterations of antioxidants in the freshwater fish Channa Punctata Bloch. J Fish Biol. 2005;67(6):1653–65.
Adams P. Management of elevated serum ferritin levels. Gastroenterology & Hepatol. 2008;4(5):333–4.
Wilhelm Filho D. Fish antioxidant defenses - a comparative approach. Braz J Med Biol Res. 1996;29(12):1735–42.
Ayala A, Muñoz MF, Argüelles S. Lipid peroxidation: production, metabolism, and signaling mechanisms of malondialdehyde and 4-hydroxy-2-nonenal. Oxidative Med Cell Longev. 2014;2014:31.
Xiao T, Shoeb M, Siddiqui MS, Zhang M, Ramana KV, Srivastava SK, Vasiliou V, Ansari NH. Molecular cloning and oxidative modification of human lens ALDH1A1: implication in impaired detoxification of lipid aldehydes. J Toxicol Environ Health A. 2009;72(9):577–84.
Obeid R. The metabolic burden of methyl donor deficiency with focus on the Betaine Homocysteine Methyltransferase pathway. Nutrients. 2013;5(9):3481–95.
Richard N, Silva TS, Wulff T, Schrama D, Dias JP, Rodrigues PML, Conceição LEC. Nutritional mitigation of winter thermal stress in gilthead seabream: associated metabolic pathways and potential indicators of nutritional state. J Proteome. 2016;142:1–14.
Malhotra JD, Kaufman RJ. Endoplasmic reticulum stress and oxidative stress: a vicious cycle or a double-edged sword? Antioxid Redox Signal. 2007;9(12):2277–93.
Cribb AE, Peyrou M, Muruganandan S, Schneider L. The endoplasmic reticulum in xenobiotic toxicity. Drug Metab Rev. 2005;37(3):405–42.
Cichoż-Lach H, Michalak A. Oxidative stress as a crucial factor in liver diseases. World J Gastroenterol. 2014;20(25):8082–91.
Wu D, Cederbaum AI. Oxidative stress and alcoholic liver disease. Semin Liver Dis. 2009;29(2):141–54.
Jeffries KM, Hinch SG, Sierocinski T, Pavlidis P, Miller KM. Transcriptomic responses to high water temperature in two species of Pacific salmon. Evol Appl. 2014;7(2):286–300.
Jeffries KM, Hinch SG, Sierocinski T, Clark TD, Eliason EJ, Donaldson MR, Li S, Pavlidis P, Miller KM. Consequences of high temperatures and premature mortality on the transcriptome and blood physiology of wild adult sockeye salmon (Oncorhynchus Nerka). Ecol Evol. 2012;2(7):1747–64.
The authors would like to thank Polly Hilder, Basseer Codabaccus and Rafael León, for technical assistance, and to Lynn Rowley and John Bateman (Murdoch Childrens Research Institute, Melbourne) for the IPA analysis.
This research was funded by Nuseed Pty Ltd. and CSIRO Ocean and Atmosphere, formerly as the Food Futures Flagship. WGN was supported by an International Postgraduate Research Scholarship funded by the Australian government.
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 PXD004985.
WGN, CGC and PDN conceived and designed the study and sampled the experimental animals. WGN performed the laboratory work. RW performed MS analysis. IRC performed statistical analysis of proteomics data. WGN interpreted the data and wrote the manuscript. All listed authors edited and approved the final manuscript.
Ethics approval and consent to participate
All procedures implemented during this experiment were approved by the University of Tasmania Animal Ethics Committee (Investigation A0015208).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
MaxQuant output files of the complete peptide-level mass spectrometry. A total of nine biological replicates (fish) per temperature treatment were analysed. (XLSX 5578 kb)
MaxQuant output files of the complete protein-level mass spectrometry. A total of nine biological replicates (fish) per temperature treatment were analysed. (XLSX 2354 kb)
Total number of quantifiable proteins. Proteins identified on the basis of two or more unique matching peptide sequences and presence in at least six of the nine biological replicates in either treatment group (XLSX 160 kb)
Full list of canonical pathways as determined by temperature-regulated proteins and predicted by IPA analysis. The significance of the association between the data set and the pathway is based on the p-value, which determines the probability that the association between the data set genes and the pathway is explained only by chance, and on the ratio value, representing the number of genes from the data found in each pathway over the total number of genes in that pathway. (XLSX 56 kb)
Components of the top five toxicity pathways as determined by temperature-regulated proteins and predicted by IPA analysis. Exp fold change and Exp p-value correspond with the fold change value and the adjusted p-value (using Benjamini Hochberg correction), respectively, reported in Additional file 3. (XLSX 52 kb)