Next Article in Journal
Commercial Vaccines Do Not Confer Protection against Two Genogroups of Piscirickettsia salmonis, LF-89 and EM-90, in Atlantic Salmon
Next Article in Special Issue
Transcriptome and Genome Analyses Applied to Aquaculture Research
Previous Article in Journal
Growth and Nutrient Utilization in Basil Plant as Affected by Applied Nutrient Quantity in Nutrient Solution and Light Spectrum
Previous Article in Special Issue
Whole-Transcriptome Analysis Identifies Gender Dimorphic Expressions of Mrnas and Non-Coding Rnas in Chinese Soft-Shell Turtle (Pelodiscus sinensis)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 11
Issue 7
10.3390/biology11070990
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Temperature Stress Effect on the Red Cusk-Eel (Geypterus chilensis) Liver: Transcriptional Modulation and Oxidative Stress Damage

by
Phillip Dettleff
1,*,
Rodrigo Zuloaga
2,3,
Marcia Fuentes
2,3,
Pamela Gonzalez
2,3,
Jorge Aedo
2,3,
Juan Manuel Estrada
4,
Alfredo Molina
2,3 and
Juan Antonio Valdés
2,3,*
1
Escuela de Medicina Veterinaria, Facultad de Agronomía e Ingeniería Forestal, Facultad de Ciencias Biológicas y Facultad de Medicina, Pontificia Universidad Católica de Chile, Santiago 7820436, Chile
2
Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago 8370186, Chile
3
Interdisciplinary Center for Aquaculture Research (INCAR), Concepción 4030000, Chile
4
Centro de Investigación Marina Quintay (CIMARQ), Universidad Andrés Bello, Quintay 2340000, Chile
*
Authors to whom correspondence should be addressed.
Biology 2022, 11(7), 990; https://doi.org/10.3390/biology11070990
Submission received: 10 May 2022 / Revised: 22 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022
(This article belongs to the Special Issue Transcriptome and Genome Analyses Applied to Aquaculture Research)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The red cusk-eel (Genypterus chilensis) is a native Chilean species important for aquaculture diversification in Chile. The effect of high-temperature stress on the liver, a key organ for fish metabolism, is unknown. In this study we determined for the first time the effects of high-temperature stress on the liver of red cusk-eel. The results showed that high-temperature stress increased hepatic enzyme activity in the plasma of stressed fish. Additionally, this stressor generated oxidative damage in liver, and generated a transcriptional response with 1239 down-regulated and 1339 up-regulated transcripts associated with several processes, including unfolded protein response, heat shock response and oxidative stress, among others. Together, these results indicate that high-temperature stress generates a relevant impact on liver, with should be considered for the aquaculture and fisheries industry of this species under a climate change scenario.

Abstract

Environmental stressors, such as temperature, are relevant factors that could generate a negative effect on several tissues in fish. A key fish species for Chilean aquaculture diversification is the red cusk-eel (Genypterus chilensis), a native fish for which knowledge on environmental stressors effects is limited. This study evaluated the effects of high-temperature stress on the liver of red cusk-eel in control (14 °C) and high-temperature (19 °C) groups using multiple approaches: determination of plasmatic hepatic enzymes (ALT, AST, and AP), oxidative damage evaluation (AP sites, lipid peroxidation, and carbonylated proteins), and RNA-seq analysis. High-temperature stress generated a significant increase in hepatic enzyme activity in plasma. In the liver, a transcriptional regulation was observed, with 1239 down-regulated and 1339 up-regulated transcripts. Additionally, high-temperature stress generated oxidative stress in the liver, with oxidative damage and transcriptional modulation of the antioxidant response. Furthermore, an unfolded protein response was observed, with several pathways enriched, as well as a heat shock response, with several heat shock proteins up regulated, suggesting candidate biomarkers (i.e., serpinh1) for thermal stress evaluation in this species. The present study shows that high-temperature stress generated a major effect on the liver of red cusk-eel, knowledge to consider for the aquaculture and fisheries of this species.

1. Introduction

Environmental factors are important for the physiology of fish, particularly those associated with water conditions. Among these, temperature, pH, and dissolved oxygen (DO) are crucial for the homeostasis of the fish in marine environments, and changes in these factors could lead to generating a stressful status for the animal [1]. It has been reported that an increase in water temperature could lead to stress and negative effects on marine fish, including salmonids [2,3], catfish [4], Atlantic cod [5], and gilthead seabream [6]. Understanding the effects of water temperature on fish is crucial in the actual scenario of climate change, considering the average rising sea temperature per decade [7]. Additionally, the effect of climate change on sea temperatures is expected to influence relevant phenomena of the Pacific coast, such as El Niño–Southern Oscillation (ENSO), increasing their intensity and frequency [8], which is relevant for marine species of the South Pacific coast, considering the increase in water temperatures associated with this phenomena. Therefore, it is important to understand the effect of an increase in temperature and how it affects the stress response in marine fish with aquaculture potential.
Chile is a relevant country in seafood production, with an important fishery and aquaculture industry. Additionally, the Chilean aquaculture sector is recognized for its important salmon and mussel industries [9]. However, in the last decade, an important effort has been made by the public private association to diversify Chilean aquaculture with native fish of commercial value. One of these species is the red cusk-eel (Genypterus chilensis), part of the Genypterus genus, endemic to the South Pacific coast, an economically relevant fish for fisheries and, recently, part of the Chilean aquaculture diversification program [10,11], with recent elucidation of the complete production cycle [12]. The red cusk-eel is a demersal fish with a carnivorous diet; it lives on rocky bottoms [13] and is characterized reproductively as a multiple spawner [11]. The high value of its flesh makes this species an attractive product; however, tons of fisheries have presented variable levels, with a decreasing tendency in the last decade [14]. In this sense, it is important to understand how environmental factors could affect this species and the mechanisms involved.
Stress in fish can be characterized as an adaptative response to danger, which generates physiological changes to prepare the fish to respond and survive threats. This stress response in fish is mediated by the hypothalamic–pituitary–interrenal (HPI) axis though several key chemical mediators, including corticotropin releasing factor, adrenocorticotropic hormone, α-melanocyte-stimulating hormone, adrenaline, and cortisol, a key hormone that increases in plasma under stress [15]. If the stressor is maintained for prolongated periods, a chronic stress status is generated in the fish, leading to several negative physiological effects, including decreased growth, reproductive problems, behavior modifications, and immune response [16,17]. Moreover, at a cellular level, the stress can lead to an increase in reactive oxygen species (ROS), which could lead to an oxidative stress status, as previously observed for several teleost fish species [18,19,20,21,22], an effect also observed for red cusk-eel in response to several stressors [23,24,25]. Red cusk-eel has shown low tolerance to intensive farming stressors, with limited information related to the stress response capacity in this species [13]. However, the specific stress response varies according to species, as well as the effect on each tissue. Our previous studies on G. chilensis showed a variable tissue response under handling stress, with altered metabolic status in the liver, a modulation of the immune response in the head kidney, and an induction of atrophy in skeletal muscle through coding and noncoding regulation [26]. Additionally, it has been observed that thermal stress could induce muscle atrophy in this species [24], as well as oxidative damage in eggs, with a minor effect on the ovary [23]. One of the most important organs for fish metabolism is the liver, which is directly involved in stress response by metabolizing and liberating stored energy to respond to stress [17]. Nevertheless, the liver response to thermal stress in Genypterus species has not been previously studied, nor has the impact of this stressor on the oxidative status of this tissue. It is important to consider that sea temperatures will increase through sea heat waves due to the effect of ENSO under a climate change scenario, which will affect the Chilean coast associated with the geographic range of G. chilensis. Therefore, the objective of this study was to evaluate the effect of high-temperature stress on the liver of G. chilensis in terms of the transcriptomic and oxidative stress status to determine the negative impact of this type of stressor on liver.

2. Materials and Methods

2.1. Ethics Statement

All procedures with the red cusk-eel individuals and all scientific activities adhered to animal welfare procedures and were approved by the bioethical committees of the Universidad Andres Bello (007/2018) and the National Commission for Scientific and Technological Research (CONICYT) of the Chilean government.

2.2. Fish Sampling and Experimental Design

In this study, we used reproductively immature red cusk-eel juveniles (G. chilensis) of 12 months of age (average weight of 665 ± 52.7 g; average length of 60 ± 4.8 cm), collected from the Centro de Investigación Marina de Quintay (CIMARQ), maintained under natural photoperiod conditions (L:D 12:12), and controlled temperature (14 ± 1 °C), and fed daily with commercial pellet food. Fish were separated into control and stress groups, with the stress group subjected to a standardized thermal stress protocol proven to generate stress in red cusk-eel [24]. Briefly, this protocol consists of increasing the temperature over 24 h at a rate of 1 °C in 5 h. This protocol maintains the thermal stress temperature (19 ± 1 °C) for 5 days. This high temperature protocol was selected considering heat waves observed in the summer season on the Chilean coast in recent years [27]. The control group was maintained at the control temperature (14 ± 1 °C) during the assay. At the end of the experiment, six individuals per group (two tanks per group, with three animals sampled per tank, total of N = 12 sampled fish) were netted and sampled. Blood samples were collected via caudal puncture using heparinized tubes, immediately centrifuged at 5000× g for 10 min for serum obtention and stored at −80 °C until analysis. After blood sampling, fish were euthanized (overdose of anesthetic 3-aminobenzoic acid ethyl ester, 300 mg/L). Fish livers were collected and stored for RNA extraction in RNAsave solution (Biological Industries, Cromwell, CT, USA) or immediately frozen in liquid nitrogen and stored at −80 °C until analysis for oxidative damage evaluation.

2.3. AST, ALT and AP Evaluation

The plasmatic activities of aspartate aminotransferase (AST), alanine aminotransferase (ALT), and alkaline phosphatase (AP) were determined using commercially available kits from Valtek (Santiago, Chile) following the manufacturer’s instruction. Briefly, these kits determined the enzymatic activity (IU/L) via the generation of colorimetric products from glutamate (colorimetric product: 450 nm), pyruvate (535 nm), and p-nitrophenol (405 nm) for AST, ALT and AP activity, respectively.

2.4. Oxidative Stress Assays in Liver

To determine the oxidative damage in the liver of red cusk-eel in response to high-temperature stress, DNA oxidative damage, protein carbonylation, and lipid peroxidation were evaluated using commercially available kits. The DNA oxidative damage assay was performed with genomic DNA (gDNA) purified from 25 mg of the liver using the DNAzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol and quantified with the Epoch Spectrophotometer System (BioTek, Winooski, VT, USA). Then, the apurinic/apyrimidinic sites (AP sites) were determined in the gDNA using the OxiSelect Oxidative DNA Damage Quantification Kit (Cell Biolabs, CA, USA) following the manufacturer’s protocol. The protein carbonylation assay was performed using total protein extracted from 100 mg of the liver in 1 mL of lysis buffer containing 50 mM Tris-HCl pH 7.4, 150 mM NaCl, 1 mM EDTA, and 1% NP-40, solubilized at 4 °C after 12,000× g centrifugation. Proteins were quantified using the Pierce BCA Protein Assay Kit (Thermo Scientific, Batavia, IL, USA). Then, carbonylated protein content was quantified using the OxiSelect Protein Carbonyl Spectrophotometric Assay (Cell Biolabs, San Diego, CA, USA) following the manufacturer’s protocol. The lipid peroxidation determination in the liver was performed using the OxiSelect HNE Adduct Competitive ELISA Kit (Cell Biolabs, San Diego, CA, USA), following the manufacturer’s protocol. This kit determines lipid peroxidation through the quantification of hydroxynonenal (HNE) protein adducts in the extracted proteins of the liver.

2.5. Liver RNA Extraction, Library Preparation and Illumina Sequencing

Total RNA was extracted from the liver stored in RNAsave solution (Biological Industries, Cromwell, CT, USA) using the TRIzol® reagent (Invitrogen, Carlsbad, CA, USA) protocol. Total RNA was quantified by fluorometry with the Qubit® RNA quantitation assay (Invitrogen, Carlsbad, CA, USA) and purity was determined according to 260/280 ratio using the Epoch Spectrophotometer System (BioTek, Winooski, VT, USA). The RNA integrity was measured according to RNA Quality Measurement Number (RQN) through a Fragment Analyzer with the Standard Sensitivity RNA Analysis kit (Advanced Analytical Technologies, Fiorenzuola, Italy), selecting samples with RQN ≥ 8. The cDNA libraries construction were performed with 1 µg of total RNA per sample using the Illumina® TruSeq RNA Sample Prep Kit v2 (Illumina®, San Diego, CA, USA), following the manufacturer’s protocol. The sizes of the mRNA libraries were determined through a Fragment Analyzer using the NGS Fragment Analysis kit (Advanced Analytical Technologies) and quantified by qPCR using the Kapa Library Quantification kit (Roche, Little Falls, NJ, USA). Paired-end sequencing (2 × 100 bp) was performed on a Hiseq 2500 (Illumina®) platform in Macrogen Inc. (Seul, South Korea).

2.6. Reads Filtering, Differential Expression, and GO Enrichment Analysis

The raw reads obtained from Illumina sequencing were trimmed to remove the remaining Illumina adapter, low-quality sequences, and short sequences (<50 bp), using the CLC Genomics Workbench v.7.0.3 software. The filtered reads were mapped to a G. chilensis reference transcriptome previously annotated by our group ([28], NCBI accession number SRS614525) using the CLC Genomics Workbench v.7.0.3 software with the following parameters: mismatch cost = 2, insertion cost = 3, deletion cost = 3, length fraction = 0.8, and a similarity fraction = 0.8. The expression values were used for clustering and heatmap chart generation with R. Differential expression analysis was performed with the R package DESeq2 (version 1.2.10) [29] to determine differentially expressed transcripts in the liver between the control and stress groups. Transcripts presenting an adjusted p-value of <0.05 and an absolute log2 fold change of >1 were considered as differentially expressed between the groups. Enrichment analysis was performed on the list of differentially expressed transcripts and GO terms to determine the overrepresented processes in response to high-temperature stress in the liver, considering up and down-expressed transcripts. This analysis was performed using the enrichment analysis tool implemented in the Blast2GO software [30]. Additionally, the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway database was used to build the represented pathways through the KEGG Automatic Annotation Server (KAAS) [31] using the KEGG Orthology (KO) identifiers of the differentially expressed list.

2.7. RNA-seq Validation by qPCR

Total RNA previously extracted was DNAse-treated to remove residual gDNA, and 1 μg of RNA was reverse transcribed into cDNA using the QuantiTect® Reverse Transcription kit (Qiagen, Germantown, MD, USA), following the manufacturer’s protocol. A total of thirteen differentially expressed transcripts were selected for qPCR validation, corresponding to: hsp60, hsp70, gpx7, ddit4, leptin, msh2, msh3, c1ql1, ccl20, atg12, atg4b, casp3 and c3. Primers were designed with Primer 3 software v0.4.0 (http://frodo.wi.mit.edu/primer3/, accessed on 6 May 2021) using the reference transcriptome previously described. The primer sequences, amplicon size, Tm, and efficiency are presented in Table 1. The qPCR was performed in a Stratagene MX3000P qPCR system (Stratagene, La Jolla, CA, USA). All qPCR assays were performed in triplicates, using no-template and no-RT controls, in compliance with the MIQE guidelines [32]. The qPCR reaction mixture contained 7.5 μL of 2× Brilliant® II SYBR® Green master mix (Agilent Technologies, CA, USA), 100 ng of cDNA per reaction and 200 nM of each primer in a 15 μL final volume. Thermal cycling conditions were an initial activation of 2 min at 95 °C, followed by 40 cycles of 30 s at 95 °C, 30 s at 62 °C, and 30 s at 72 °C. A melt curve analysis was performed to confirm a single qPCR product and a standard curve using two-fold series dilutions was used to estimate the efficiency of each primer set. The expression of target genes was normalized using the geometric means of two reference genes (actb and taf12) previously validated for red cusk-eel in the liver [26] and following the methodology described by [33].

2.8. Statistical Analysis

Significant differences between means of the control and stress groups for DNA damage (AP sites), enzymatic activity (AST, ALT, and AP) and differential expression of the qPCR validated genes were determined using a t-test with a significance threshold of p < 0.05. All statistical analyses were performed using GraphPad Prism, v.5.00 (GraphPad Software, San Diego, CA, USA).

3. Results

3.1. Hepatic Enzyme Activity and Oxidative Stress Response to High-Temperature Stress

We previously reported that thermal stress for 5 days significantly increased plasmatic levels of cortisol and glucose in G. chilensis [24]. To understand how thermal stress affects the metabolism in this species, we studied the effect of high temperature on the liver. To evaluate the effect on the hepatic function of this type of stressor, we measured plasmatic markers of liver damage, i.e., ALT, AST, and AP enzymatic activity. High-temperature stress significantly increased the plasmatic activity of ALT (Figure 1A), AST (Figure 1B) and AP (Figure 1C) in the stress group, evidencing altered hepatic function in the liver of stressed fish. To determine the oxidative damage in the liver generated by high-temperature stress, the DNA oxidative damage, protein carbonylation and lipid peroxidation were determined in the liver. High temperatures generated DNA damage, evidenced by the significant increase in apurinic/apyrimidinic sites in the stressed group (Figure 2A). Oxidative damage was also observed in lipid and proteins, determined by a significant increase in protein carbonylation (Figure 2B) and lipid peroxidation (Figure 2C) in response to high temperature in the stress group.

3.2. Differentially Expressed Transcripts in Hepatic Response to High-Temperature Stress

To understand the complexity of the stress response in the liver of G. chilensis, we performed RNA-seq analysis on the liver of each experimental group. The sequencing generated a total of 754,678,455 paired-end reads, with an average of 58,052,189 raw paired-end reads per library. The raw data are available from NCBI under BioProject PRJNA835467 with BioSamples accession number SAMN28102858, SAMN28102859, SAMN28102860, SAMN28102861, SAMN28102862, SAMN28102863, SAMN28102864, SAMN28102865, SAMN28102866, and SAMN28102867. After trimming by quality, adapters, and size, we obtained an average of 58,018,122 high-quality filtered paired-end reads per library (Table 2). These reads were mapped to the G. chilensis reference transcriptome (NCBI accession number SRS614525), obtaining an average of 85.2% mapped reads (Table 2). Expression values were used for normalization and differential expression analysis with the R package DESeq2 (version 1.2.10) [29], obtaining a total of 2578 differentially expressed transcripts between control and stressed groups. Of these transcripts, 1239 were down-regulated (Table S1) and 1339 were up-regulated (Table S2) in the stressed group (Figure 3 and Figure S1), evidencing five clusters of transcripts with different patterns of expression for the control and thermal stress groups (Figure S2).

3.3. GO Enrichment and Pathway Analysis in the Liver

The differentially expressed transcripts were used for an enrichment analysis using GO terms to determine the overrepresented processes in response to high-temperature stress in the up- and down-regulated transcripts of the liver. The enrichment analysis showed 38 enriched processes for up-regulated transcripts (Table 3), including 18 biological processes (BP), 16 cellular components (CC), and 4 molecular functions (MF). Several of the enriched terms were related to protein metabolism, including protein folding (GO:0006457), protein transport (GO:0015031), protein localization (GO:0008104), protein retention in ER lumen (GO:0006621), and unfolded protein binding (GO:0051082). Additionally, the BP term response to heat (GO:0009408) was also enriched in the up-regulated transcripts in the high-temperature stress group. No enriched GO terms were identified in the down-regulated transcripts (FDR < 0.05).
In the KEGG pathway analysis, we identified several pathways represented in the down-regulated and up-regulated transcripts, including apoptosis, cell cycle, MAPK signaling pathway, autophagy, protein processing in endoplasmic reticulum, and ubiquitin-mediated proteolysis. Among these pathways, protein processing in endoplasmic reticulum was one of the most represented on the differentially expressed transcripts (Figure 4), with 31 up-regulated and three down-regulated transcripts represented in this pathway, including several heat shock proteins; heat shock protein 40 (hsp40, also known as dnaJ homolog subfamily B member 1), heat shock protein 70 (hsp70) and heat shock protein 90 (hsp90).

3.4. qPCR Validation of Differentially Expressed Genes

The validation of the RNA-seq results was performed by selecting 13 differentially expressed transcripts between the control and stressed groups, including heat shock protein 60 (hsp60), hsp70, DNA damage-inducible transcript 4 (ddit4), glutathione peroxidase 7-like (gpx7), leptin, DNA mismatch repair Msh2-like (msh2), DNA mismatch repair Msh3-like (msh3), complement C1q 2-like (clql1), C-C motif chemokine 20-like (ccl20), ubiquitin ATG12-like (atg12), cysteine protease ATG4B-like (atg4b), caspase 3 (casp3) and complement C3-like (c3) (Table 1). All of the selected differentially expressed transcripts on RNA-seq were validated by qPCR, confirming the results of the RNA-seq analysis (Figure 5), with a Pearson’s correlation coefficient of r = 0.847, confirming the significant differential expression of all tested transcripts by qPCR between the control and high-temperature group (Figure S3).

4. Discussion

Stress in fish is a relevant issue for aquaculture species and native populations of fish [15]. The environmental stress factor associated with water parameters represents a key issue in fish, especially considering the variations in the environment in the short- and long-term associated with global warming, climate change, ENSO, and pollution of the oceans, which could lead to modification of sea temperature, pH, and DO level, as well as increased levels of pollutants, including microplastics and toxic compounds, which could generate stress in teleost fish [7,8,34,35,36]. These environmental stressors could also affect the red cusk-eel, a relevant species for Chilean fisheries and aquaculture diversification. However, studies aimed at understanding the effect of stressors in Genypterus species are limited [25,26,37,38], with no information about the impact of thermal stress on the liver for Genypterus species. In this sense, thermal stress has been previously studied in other tissues of Genypterus species, including the skeletal muscle, ovary, and post-ovulatory eggs [23,24]. Therefore, to understand the effect of high-temperature stress at the hepatic level, we evaluated the hepatic enzymes, oxidative stress response, and transcriptional regulation in red cusk-eel.

4.1. High-Temperature Effect on Hepatic Enzymes

In a previous study, we determined that high temperature (19 °C) could generate a stress response in red cusk-eel with an increase in the plasmatic level of cortisol and glucose [24]. In the present study, we determined that a high temperature increased the plasmatic activity of ALT, AST, and AP enzymes, evidencing the effect at the hepatic level. This effect was previously observed in other teleost fish under stress conditions, including the plasmatic activity of ALT and AST under high densities and metal pollution in rohu (Labeo rohita) [39,40], metal toxicity and pesticides in Nile tilapia (Oreochromis niloticus) [41,42] and spotted snakehead (Channa punctatus) [43], and elevated ALT, AST and AP under handling stress in red cusk-eel [25]. In terms of high-temperature stress, different results have been observed, with an increase in the plasmatic activity of AST and ALT in pufferfish (Takifugu obscurus) [44], which is consistent with our results for red cusk-eel, while no effect for AST and ALT was observed in Turbot (Scophthalmus maximus) [45] under high-temperature stress, evidencing that the hepatic response to thermal stress could vary according to fish species. The increased levels on ALT, AST an AP enzymes could be indicative of liver dysfunction, reflecting hepatocyte damage due to thermal stress, which is concordant with the results observed in mammals [46].

4.2. Oxidative Stress under High-Temperature Stress

Oxidative stress corresponds to a disturbance between the production of ROS, which can accumulate in cells, and the antioxidant defenses generated by the cellular systems to detoxify these ROS [47]. OS participate in different normal cellular functions, acting as a second messenger in signal transduction [48], but they can also generate cellular damage, including oxidative damage generating lipid peroxidation, as well as damage to proteins and DNA [49]. Environmental stress, including thermal stress, can lead to oxidative stress status in marine animals [3]. It was previously observed that thermal stress could modulate the oxidative stress status in several tissues of teleost fish [50,51]. One of the relevant organs in which oxidative status could be affected by thermal stress is the liver, as it was previously observed that thermal stress by low temperature could increase antioxidant enzymes in milkfish (Chanos Chanos) [52], as well as generate oxidative damage, leading to lipid peroxidation in Pacu (Piaractus mesopotamicus) under low-temperature stress [53,54]. High-temperature stress could also generate a relevant impact on the oxidative status of fish, with increased antioxidant enzyme activities, as observed for Senegalese sole (Solea senegalensis) [55]. Additionally, oxidative damage was observed in the liver of rohu (Labeo rohita), with lipid peroxidation and DNA fragmentation [56]. This is consistent with the findings of our study, evidencing that thermal stress generates an important oxidative effect on the liver of red cusk-eel. This response was also observed at the transcriptional level, with the up-regulation of GPx genes in several tissues of teleost fish under thermal stress, including black porgy (Acanthopagrus schlegeli) [57] and pufferfish [58], as well as in red cusk-eel eggs under thermal stress [23], where gpx1 was increased, an effect not observed for this species in skeletal muscle [24]. Additionally, the effect of high-temperature stress on the liver related to DNA damage was also observed at the transcriptional level, with the up-regulation of genes involved in DNA mismatch repair (msh2 and msh3), concordant with the previously reported effect of thermal stress on zebrafish (Danio rerio) [59] and American lobster [60]. However, the thermal stress response associated with oxidative stress in teleost fish could vary according to species and tissues, as observed in sheepshead minnow (Cyprinodon variegatus), where a limited effect on antioxidant enzymes and no lipid peroxidation were present [61], in contrast to the variable lipid peroxidation observed in Senegalese sole [55]. We previously observed in red cusk-eel that the impact of high-temperature stress on oxidative damage could vary according to tissue, with lipid peroxidation and DNA damage observed for skeletal muscle [24], but no oxidative damage was observed in the ovaries under high-temperature stress [23], showing that the liver is a sensitive organ under thermal stress in red cusk-eel under an ENSO temperature increase scenario.

4.3. High-Temperature Stress in Hepatic Protein Processing and Folding

Temperature can modulate several cellular processes through gene expression regulation [2]. Here, we used RNA-seq analysis to evaluate, for the first time, the effect of thermal stress on the liver of red cusk-eel, considering that no studies have previously evaluated the temperature effect on the liver in any species of the Genypterus genus. We observed a higher hepatic transcriptional response to high-temperature stress compared to other types of stressors, such as handling stress, previously observed for red cusk-eel (4.6 times the differentially expressed transcripts) [25], evidencing that a high temperature could have a higher transcriptional impact on the liver than others stressors in this species. At the cellular level, the function of endoplasmic reticulum is key to protein synthesis, folding, and exporting [62]. However, external processes such as thermal stress could generate alterations in homeostasis, affecting normal protein folding, which leads to endoplasmic reticulum (ER) stress. To alleviate this stress, the unfolded protein response molecular mechanism is activated in cells [63]. This ER stress and subsequent unfolded protein response is concordant with our results in the liver, where we observed enriched processes associated with this type of response in the up-regulated genes in the stress group, including protein folding, protein transport, protein localization, protein retention in ER lumen, and unfolded protein binding. In this sense, we found several genes associated with protein processing and folding pathways, including transcripts involved in: protein export, such as signal recognition particle 14 and 19 kDa (srp14 and srp19) and signal peptidase complex subunit 2-like (spcs2) [64]; transcripts associated with protein processing in endoplasmic reticulum (Figure 4), with part of this pathway down-regulated, including TNF receptor-associated factor 2-like (traf2), and mainly an up-regulation of this pathway, including ubiquilin-4 (ubqln4), cytoskeleton-associated protein 4 (ckap4); ribosome-binding protein 1 (rrbp1, also known as p180); and glucosidase 2 subunit beta (prkcsh). Others include transcripts associated with the unfolded protein response, including eukaryotic translation initiation factor 2 alpha kinase 1 (eik2ak1, also known as hri), 3 (eik2ak3, also known as perk) and 5 (eik2ak5), cyclic AMP-dependent transcription factor ATF-4-like (atf4) and DNA damage-inducible transcript 3 -like (ddit4, also known as chop). The unfolded protein response could be initiated by EIK2AK3 kinase activation of eIF2α, leading to ribosome inhibition and attenuating protein synthesis [65]. Additionally, the ATF4 gene was activated, which would lead to DDIT4 gene regulation and an antioxidant response in cells [66], evidencing that thermal stress in the liver of red cusk-eel modulates the protein processing and generates an unfolded protein response that is not able to control the oxidative stress and damage in this tissue. This is concordant with the results observed in previous studies where unfolded protein response genes were activated under thermal stress in mammals [67,68] and fish, e.g., in the liver of Tambaqui (Colossoma macropomum) [69] and gilthead sea bream (Sparus aurata) [70].

4.4. Heat Shock Protein as Thermal Stress Biomarkers

The cellular response to thermal stress is a key process to preserve protein integrity; this is known as the heat shock response, which includes heat shock proteins (Hsps) to re-fold the proteins damaged by temperature [71]. This response was observed in the liver of red cusk-eel with the enrichment of the term response to heat (GO:0009408), highlighting several heat shock proteins differentially expressed in response to high-temperature stress, including hsp40, hsp60, hsp70, hsp90, and serpinh1. The up-regulation of these Hsps was previously observed in other fish under thermal stress, with an increase in the liver of three-spined stickleback (Gasterosteus aculeatus) (hsp60, hsp70, and hsp90) [72], Atlantic salmon (Salmo salar) (hsp70) [73,74], Atlantic cod (Gadus morhua) (hsp70 and hsp90) [74], and Wuchang bream (Megalobrama amblycephala) (hsp60, hsp70, and hsp90) [75]. The hsp60 and hsp70 expression levels were also regulated under thermal stress in the postovulatory eggs and skeletal muscle of red cusk-eel [23,24], showing that these genes could be valuable biomarkers of thermal stress in this species. However, the most up-regulated gene in the liver under high-temperature stress was serpinh1 (5.7log2 Fold Change, Table S2). SerpinH1 (also known as Hsp47) is a chaperone involved in the biosynthesis of collagen at the ER [76], with a key role in the restoration of homeostasis during high-temperature stress and oxidative stress [22]. It was shown to be a good biomarker for thermal stress in several salmonid species, including rainbow trout (Oncorhynchus mykiss), sockeye salmon (Oncorhynchus nerka), and Chinook salmon (Oncorhynchus tshawytscha) [3,77], as well as zebrafish [78], evidencing that sherpinh1 is one of the most suitable transcriptional biomarkers of high-temperature stress in fish. This was also the case for red cusk-eel in this study, representing a useful tool to evaluate thermal stress status in this species under a climate change scenario.

5. Conclusions

The present study evaluated for the first time the effects of high-temperature stress in the liver of G. chilensis, using a multiple approach of plasmatic hepatic enzymes, oxidative damage evaluation and RNA-seq analysis. We showed that high-temperature stress under heatwaves ENSO-associated scenario generated a major effect in the liver, affecting hepatic enzymes, generating oxidative damage in this tissue, as well as generating an unfolded protein response at the molecular level in several associated pathways, including a heat shock response, evidencing the affection of red cusk-eel under this type of stressor. This study contributes to knowledge about thermal stress under a climate change scenario, generating candidate biomarkers for thermal stress evaluation in this species, information that should be relevant for the aquaculture and fisheries industry of red cusk-eel.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11070990/s1, Figure S1: Differentially expressed transcripts in response to high-temperature stress in G. chilensis expressed as average transcription per group; Figure S2: Differentially expressed transcripts and clustering groups in response to high-temperature stress in G. chilensis; Figure S3: Validation of selected differentially expressed transcripts by qPCR on liver of G. chilensis on response to high-temperature stress; Table S1: Down-regulated differentially expressed transcripts between control and stress groups of G. chilensis; Table S2: Up-regulated differentially expressed transcripts between control and stress groups of G. chilensis.

Author Contributions

Conceptualization and methodology, P.D. and J.A.V.; bioinformatic analysis, P.D.; animal experiment and sampling, P.D., R.Z., M.F., P.G., J.A. and J.M.E.; resources and funding acquisition, P.D., J.M.E., J.A.V. and A.M., writing-draft preparation, P.D. and J.A.V.; writing—review and editing, P.D., J.A.V., R.Z. and J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by CONICYT/FONDAP (grant number 15110027) awarded to Juan Antonio Valdés and Alfredo Molina, CONICYT FONDECYT Postdoctorado (grant number 3180283) awarded to Phillip Dettleff, and CONICYT FONDECYT Regular (grant number 1201498).

Institutional Review Board Statement

All procedures with the red cusk-eel individuals in this study were conducted according to the guidelines of the Declaration of Helsinki and approved by the Bioethics Committee of Universidad Andres Bello and the National Commission for Scientific and Technological Research (CONICYT) of the Chilean government (N 007/2018).

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw read sequences obtained from liver sequencing of red cusk-eel were deposited in NCBI under BioProject PRJNA835467, with BioSample accession number SAMN28102858, SAMN28102859, SAMN28102860, SAMN28102861, SAMN28102862, SAMN28102863, SAMN28102864, SAMN28102865, SAMN28102866, and SAMN28102867. The datasets generated and/or analyzed in the present study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank the staff at CIMARQ, Universidad Andres Bello, and Italo Barrios for their assistance with experiments and fish management.

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

References

  1. Somero, G.N. The physiology of climate change: How potentials for acclimatization and genetic adaptation will determine ‘winners’ and ‘losers’. J. Exp. Biol. 2010, 213, 912–920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Huang, J.; Li, Y.; Liu, Z.; Kang, Y.; Wang, J. Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney. Fish Shellfish Immunol. 2018, 82, 32–40. [Google Scholar] [CrossRef] [PubMed]
  3. Akbarzadeh, A.; Gunther, O.P.; Houde, A.L.; Li, S.; Ming, T.J.; Jeffries, K.M.; Hinch, S.G.; Miller, K.M. Developing specific molecular biomarkers for thermal stress in salmonids. BMC Genom. 2018, 19, 749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Dalvi, R.S.; Pal, A.K.; Tiwari, L.R.; Baruah, K. Influence of acclimation temperature on the induction of heat-shock protein 70 in the catfish Horabagrus brachysoma (Gunther). Fish Physiol. Biochem. 2012, 38, 919–927. [Google Scholar] [CrossRef]
  5. Aursnes, I.A.; Rishovd, A.L.; Karlsen, H.E.; Gjoen, T. Validation of reference genes for quantitative RT-qPCR studies of gene expression in Atlantic cod (Gadus morhua L.) during temperature stress. BMC Res. Notes 2011, 4, 104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Madeira, D.; Vinagre, C.; Costa, P.M.; Diniz, M.S. Histopathological alterations, physiological limits, and molecular changes of juvenile Sparus aurata in response to thermal stress. Mar. Ecol. Prog. Ser. 2014, 505, 256–266. [Google Scholar] [CrossRef]
  7. Dobretsov, S.; Coutinho, R.; Rittschof, D.; Salta, M.; Ragazzola, F.; Hellio, C. The oceans are changing: Impact of ocean warming and acidification on biofouling communities. Biofouling 2019, 35, 585–595. [Google Scholar] [CrossRef]
  8. Wang, B.; Luo, X.; Yang, Y.M.; Sun, W.; Cane, M.A.; Cai, W.; Yeh, S.W.; Liu, J. Historical change of El Nino properties sheds light on future changes of extreme El Nino. Proc. Natl. Acad. Sci. USA 2019, 116, 22512–22517. [Google Scholar] [CrossRef] [Green Version]
  9. Flores-Kossack, C.; Montero, R.; Kollner, B.; Maisey, K. Chilean aquaculture and the new challenges: Pathogens, immune response, vaccination and fish diversification. Fish Shellfish Immunol. 2020, 98, 52–67. [Google Scholar] [CrossRef]
  10. Vega, R.; Pradenas, M.; Estrada, J.M.; Ramirez, D.; Valdebenito, I.; Mardones, A.; Dantagnan, P.; Alfaro, D.; Encina, F.; Pichara, C. Evaluation and comparison of the efficiency of two incubation systems for Genypterus chilensis (Guichenot, 1848) eggs. Lat. Am. J. Aquat. Res. 2012, 40, 187–200. [Google Scholar] [CrossRef]
  11. Chong, J.; Gonzalez, P. Reproductive cycle and maturity mean size of the red cusk eel, Genypterus chilensis (Guichenot, 1881) in the coast off Talcahuano, Chile. Rev. Biol. Mar. Oceanogr. 2009, 44, 257–262. [Google Scholar] [CrossRef]
  12. Gonzalez, P.; Dettleff, P.; Valenzuela, C.; Estrada, J.M.; Valdes, J.A.; Meneses, C.; Molina, A. Evaluating the genetic structure of wild and commercial red cusk-eel (Genypterus chilensis) populations through the development of novel microsatellite markers from a reference transcriptome. Mol. Biol. Rep. 2019, 46, 5875–5882. [Google Scholar] [CrossRef] [PubMed]
  13. Vega, R.; Estrada, J.M.; Ramirez, D.; Flores, C.; Zamorano, J.; Encina, F.; Mardones, A.; Valdebenito, I.; Dantagnan, P. Growth of cusk eel Genypterus chilensis juveniles in culture conditions. Lat. Am. J. Aquat. Res. 2015, 43, 344–350. [Google Scholar] [CrossRef]
  14. SERNAPESCA. Anuario Estadístico de Pesca y Acuicultura 2020; Ministeriode Economía: Santiago, Chile, 2021. [Google Scholar]
  15. Aluru, N.; Vijayan, M.M. Stress transcriptomics in fish: A role for genomic cortisol signaling. Gen. Comp. Endocrinol. 2009, 164, 142–150. [Google Scholar] [CrossRef]
  16. Ellis, T.; Yildiz, H.Y.; Lopez-Olmeda, J.; Spedicato, M.T.; Tort, L.; Overli, O.; Martins, C.I.M. Cortisol and finfish welfare. Fish Physiol. Biochem. 2012, 38, 163–188. [Google Scholar] [CrossRef]
  17. Faught, E.; Aluru, N.; Vijayan, M. The molecular stress response. In Biology of Stress in Fish; Schreck, C.B., Tort, L., Farrell, A.P., Brauner, C.J., Eds.; Elsevier: Amsterdam, The Netherlands, 2016; Volume 35, pp. 114–149. [Google Scholar]
  18. Logan, C.A.; Somero, G.N. Effects of thermal acclimation on transcriptional responses to acute heat stress in the eurythermal fish Gillichthys mirabilis (Cooper). Am. J. Physiol.-Regul. Integr. Comp. Physiol. 2011, 300, R1373–R1383. [Google Scholar] [CrossRef] [Green Version]
  19. Wu, S.M.; Liu, J.H.; Shu, L.H.; Chen, C.H. Anti-oxidative responses of zebrafish (Danio rerio) gill, liver and brain tissues upon acute cold shock. Comp. Biochem. Phys. A 2015, 187, 202–213. [Google Scholar] [CrossRef]
  20. Yang, C.G.; Jiang, M.; Wen, H.; Tian, J.; Liu, W.; Wu, F.; Gou, G.W. Analysis of differential gene expression under low-temperature stress in Nile tilapia (Oreochromis niloticus) using digital gene expression. Gene 2015, 564, 134–140. [Google Scholar] [CrossRef]
  21. Nakano, T.; Kameda, M.; Shoji, Y.; Hayashi, S.; Yamaguchi, T.; Sato, M. Effect of severe environmental thermal stress on redox state in salmon. Redox Biol. 2014, 2, 772–776. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, Y.; Liu, Z.; Li, Z.; Shi, H.; Kang, Y.; Wang, J.; Huang, J.; Jiang, L. Effects of heat stress on respiratory burst, oxidative damage and SERPINH1 (HSP47) mRNA expression in rainbow trout Oncorhynchus mykiss. Fish Physiol. Biochem. 2016, 42, 701–710. [Google Scholar] [CrossRef]
  23. Dettleff, P.; Zuloaga, R.; Fuentes, M.; Estrada, J.M.; Molina, A.; Valdes, J.A. Temperature effect on oxidative stress and egg quality-related genes on post-ovulatory eggs and ovary of red cusk-eel (Genypterus chilensis). J. Fish Biol. 2021, 98, 1475–1480. [Google Scholar] [CrossRef] [PubMed]
  24. Dettleff, P.; Zuloaga, R.; Fuentes, M.; González, J.; Aedo, J.; Estrada, J.M.; Molina, A.; Valdés, J.A. Physiological and molecular responses to thermal stress in red cusk-eel (Genypterus chilensis) juveniles reveals atrophy and oxidative damage in skeletal muscle. J. Therm. Biol. 2020, 94, 102750. [Google Scholar] [CrossRef] [PubMed]
  25. Naour, S.; Espinoza, B.M.; Aedo, J.E.; Zuloaga, R.; Maldonado, J.; Bastias-Molina, M.; Silva, H.; Meneses, C.; Gallardo-Escarate, C.; Molina, A.; et al. Transcriptomic analysis of the hepatic response to stress in the red cusk-eel (Genypterus chilensis): Insights into lipid metabolism, oxidative stress and liver steatosis. PLoS ONE 2017, 12, e0176447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Dettleff, P.; Hormazabal, E.; Aedo, J.; Fuentes, M.; Meneses, C.; Molina, A.; Valdes, J.A. Identification and Evaluation of Long Noncoding RNAs in Response to Handling Stress in Red Cusk-Eel (Genypterus chilensis) via RNA-seq. Mar. Biotechnol. 2020, 22, 94–108. [Google Scholar] [CrossRef] [PubMed]
  27. SHOA. Servicio Hidrográfico y Oceanográfico de la Armada de Chile 2019; Servicio Hidrografico y Oceanografico de la Armada: Valparaiso, Chile, 2019. [Google Scholar]
  28. Aedo, J.E.; Maldonado, J.; Estrada, J.M.; Fuentes, E.N.; Silva, H.; Gallardo-Escarate, C.; Molina, A.; Valdes, J.A. Sequencing and de novo assembly of the red cusk-eel (Genypterus chilensis) transcriptome. Mar. Genom. 2014, 18, 105–107. [Google Scholar] [CrossRef]
  29. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  30. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [Green Version]
  31. Moriya, Y.; Itoh, M.; Okuda, S.; Yoshizawa, A.C.; Kanehisa, M. KAAS: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35, W182–W185. [Google Scholar] [CrossRef] [Green Version]
  32. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [Green Version]
  33. Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, research0034.1. [Google Scholar] [CrossRef] [Green Version]
  34. Doney, S.C.; Fabry, V.J.; Feely, R.A.; Kleypas, J.A. Ocean acidification: The other CO2 problem. Annu. Rev. Mar. Sci. 2009, 1, 169–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Alimi, O.S.; Farner Budarz, J.; Hernandez, L.M.; Tufenkji, N. Microplastics and Nanoplastics in Aquatic Environments: Aggregation, Deposition and Enhanced Contaminant Transport. Environ. Sci. Technol. 2018, 52, 1704–1724. [Google Scholar] [CrossRef] [PubMed]
  36. Bhagat, J.; Nishimura, N.; Shimada, Y. Toxicological interactions of microplastics/nanoplastics and environmental contaminants: Current knowledge and future perspectives. J. Hazard. Mater. 2021, 405, 123913. [Google Scholar] [CrossRef] [PubMed]
  37. Aballai, V.; Aedo, J.E.; Maldonado, J.; Bastias-Molina, M.; Silva, H.; Meneses, C.; Boltana, S.; Reyes, A.; Molina, A.; Valdes, J.A. RNA-seq analysis of the head-kidney transcriptome response to handling-stress in the red cusk-eel (Genypterus chilensis). Comp. Biochem. Phys. D 2017, 24, 111–117. [Google Scholar] [CrossRef]
  38. Aedo, J.E.; Maldonado, J.; Aballai, V.; Estrada, J.M.; Bastias-Molina, M.; Meneses, C.; Gallardo-Escarate, C.; Silva, H.; Molina, A.; Valdes, J.A. mRNA-seq reveals skeletal muscle atrophy in response to handling stress in a marine teleost, the red cusk-eel (Genypterus chilensis). BMC Genom. 2015, 16, 1024. [Google Scholar] [CrossRef] [Green Version]
  39. Pakhira, C.; Nagesh, T.S.; Abraham, T.J.; Dash, G.; Behera, S. Stress responses in rohu, Labeo rohita transported at different densities. Aquac. Rep. 2015, 2, 39–45. [Google Scholar] [CrossRef] [Green Version]
  40. Kumar, N.; Krishnani, K.K.; Singh, N.P. Oxidative and Cellular Metabolic Stress of Fish: An Appealing Tool for Biomonitoring of Metal Contamination in the Kolkata Wetland, a Ramsar Site. Arch. Environ. Contam. Toxicol. 2019, 76, 469–482. [Google Scholar] [CrossRef]
  41. Veisi, S.; Sarkheil, M.; Johari, S.A.; Safari, O. Dietary supplementation with melatonin: Influence on growth performance, oxidative stress status, and amelioration of silver nanoparticles-induced toxicity in Nile tilapia (Oreochromis niloticus). Trop. Anim. Health Prod. 2021, 53, 314. [Google Scholar] [CrossRef]
  42. Firat, O.; Tutus, R. Comparative Acute Toxicity Assessment of Organophosphate and Avermectin Insecticides on a Freshwater Fish Oreochromis niloticus. Bull. Environ. Contam. Toxicol. 2020, 105, 582–587. [Google Scholar] [CrossRef]
  43. Bharti, S.; Rasool, F. Analysis of the biochemical and histopathological impact of a mild dose of commercial malathion on Channa punctatus (Bloch) fish. Toxicol. Rep. 2021, 8, 443–455. [Google Scholar] [CrossRef]
  44. Cheng, C.H.; Guo, Z.X.; Luo, S.W.; Wang, A.L. Effects of high temperature on biochemical parameters, oxidative stress, DNA damage and apoptosis of pufferfish (Takifugu obscurus). Ecotoxicol. Environ. Saf. 2018, 150, 190–198. [Google Scholar] [CrossRef] [PubMed]
  45. Yang, S.; Zhao, T.; Ma, A.; Huang, Z.; Liu, Z.; Cui, W.; Zhang, J.; Zhu, C.; Guo, X.; Yuan, C. Metabolic responses in Scophthalmus maximus kidney subjected to thermal stress. Fish Shellfish Immunol. 2020, 103, 37–46. [Google Scholar] [CrossRef] [PubMed]
  46. Das, A. Heat stress-induced hepatotoxicity and its prevention by resveratrol in rats. Toxicol. Mech. Methods 2011, 21, 393–399. [Google Scholar] [CrossRef] [PubMed]
  47. Pizzino, G.; Irrera, N.; Cucinotta, M.; Pallio, G.; Mannino, F.; Arcoraci, V.; Squadrito, F.; Altavilla, D.; Bitto, A. Oxidative Stress: Harms and Benefits for Human Health. Oxid. Med. Cell. Longev. 2017, 2017, 8416763. [Google Scholar] [CrossRef]
  48. Betteridge, D.J. What is oxidative stress? Metab. Clin. Exp. 2000, 49, 3–8. [Google Scholar] [CrossRef]
  49. Lesser, M.P. Oxidative stress in marine environments: Biochemistry and physiological ecology. Annu. Rev. Physiol. 2006, 68, 253–278. [Google Scholar] [CrossRef] [Green Version]
  50. Madeira, D.; Narciso, L.; Cabral, H.N.; Vinagre, C.; Diniz, M.S. Influence of temperature in thermal and oxidative stress responses in estuarine fish. Comparative Biochemistry and Physiology. Part A Mol. Integr. Physiol. 2013, 166, 237–243. [Google Scholar] [CrossRef]
  51. Vinagre, C.; Madeira, D.; Narciso, L.; Cabral, H.N.; Diniz, M.S. Effect of temperature on oxidative stress in fish: Lipid peroxidation and catalase activity in the muscle of juvenile seabass, Dicentrarchus labrax. Ecol. Indic. 2012, 23, 274–279. [Google Scholar] [CrossRef]
  52. Chang, C.H.; Mayer, M.; Rivera-Ingraham, G.; Blondeau-Bidet, E.; Wu, W.Y.; Lorin-Nebel, C.; Lee, T.H. Effects of temperature and salinity on antioxidant responses in livers of temperate (Dicentrarchus labrax) and tropical (Chanos Chanos) marine euryhaline fish. J. Therm. Biol. 2021, 99, 103016. [Google Scholar] [CrossRef]
  53. Ale, A.; Bacchetta, C.; Rossi, A.S.; Scarabotti, P.A.; Cazenave, J. Low temperature stress in a cultured fish (Piaractus mesopotamicus) fed with Pyropia columbina red seaweed-supplemented diet. Fish Physiol. Biochem. 2021, 47, 829–839. [Google Scholar] [CrossRef]
  54. Bacchetta, C.; Ale, A.; Rossi, A.S.; Karakachoff, M.; Cazenave, J. Effects of cold stress on juvenile Piaractus mesopotamicus and the mitigation by beta-carotene. J. Therm. Biol. 2020, 88, 102497. [Google Scholar] [CrossRef] [PubMed]
  55. Castro, C.; Perez-Jimenez, A.; Guerreiro, I.; Peres, H.; Castro-Cunha, M.; Oliva-Teles, A. Effects of temperature and dietary protein level on hepatic oxidative status of Senegalese sole juveniles (Solea senegalensis). Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2012, 163, 372–378. [Google Scholar] [CrossRef] [PubMed]
  56. Roychowdhury, P.; Aftabuddin, M.; Pati, M.K. Thermal stress-induced oxidative damages in the liver and associated death in fish, Labeo rohita. Fish Physiol. Biochem. 2021, 47, 21–32. [Google Scholar] [CrossRef] [PubMed]
  57. Lopez-Anido, R.N.; Harrington, A.M.; Hamlin, H.J. Coping with stress in a warming Gulf: The postlarval American lobster’s cellular stress response under future warming scenarios. Cell Stress Chaperones 2021, 26, 721–734. [Google Scholar] [CrossRef]
  58. An, K.W.; Kim, N.N.; Shin, H.S.; Kil, G.S.; Choi, C.Y. Profiles of antioxidant gene expression and physiological changes by thermal and hypoosmotic stresses in black porgy (Acanthopagrus schlegeli). Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2010, 156, 262–268. [Google Scholar] [CrossRef]
  59. Cheng, C.H.; Yang, F.F.; Liao, S.A.; Miao, Y.T.; Ye, C.X.; Wang, A.L.; Tan, J.W.; Chen, X.Y. High temperature induces apoptosis and oxidative stress in pufferfish (Takifugu obscurus) blood cells. J. Therm. Biol. 2015, 53, 172–179. [Google Scholar] [CrossRef]
  60. Ho, T.N.; Paul, G.V.; Chen, Y.H.; Hsu, T. Heat stress upregulates G-T mismatch binding activities in zebrafish (Danio rerio) embryos preexposed and nonexposed to a sublethal level of cadmium (Cd). Chemosphere 2019, 218, 179–188. [Google Scholar] [CrossRef]
  61. Baker, B.P.; Van Wie, I.; Braun, E.; Jimenez, A.G. Thermal stability vs. variability: Insights in oxidative stress from a eurytolerant fish. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2020, 249, 110767. [Google Scholar] [CrossRef]
  62. Anelli, T.; Sitia, R. Protein quality control in the early secretory pathway. EMBO J. 2008, 27, 315–327. [Google Scholar] [CrossRef] [Green Version]
  63. Malini, M.K.; Lekshmy, V.S.; Pal, M.; Chinnusamy, V.; Nagaraj Kumar, M. Unfolded protein response (UPR) mediated under heat stress in plants. Plant Physiol. Rep. 2020, 25, 569–582. [Google Scholar] [CrossRef]
  64. Faoro, C.; Ataide, S.F. Noncanonical Functions and Cellular Dynamics of the Mammalian Signal Recognition Particle Components. Front. Mol. Biosci. 2021, 8, 679584. [Google Scholar] [CrossRef] [PubMed]
  65. Harding, H.P.; Zhang, Y.H.; Ron, D. Protein translation and folding are coupled by an endoplasmic-reticulum-resident kinase. Nature 1999, 397, 271–274. [Google Scholar] [CrossRef] [PubMed]
  66. Vattem, K.M.; Wek, R.C. Reinitiation involving upstream ORFs regulates ATF4 mRNA translation in mammalian cells. Proc. Natl. Acad. Sci. USA 2004, 101, 11269–11274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Adachi, M.; Liu, Y.; Fujii, K.; Calderwood, S.K.; Nakai, A.; Imai, K.; Shinomura, Y. Oxidative Stress Impairs the Heat Stress Response and Delays Unfolded Protein Recovery. PLoS ONE 2009, 4, A57–A66. [Google Scholar] [CrossRef] [Green Version]
  68. Ahmed, K.; Zaidi, S.F.; Mati-ur-Rehman; Rehman, R.; Kondo, T. Hyperthermia and protein homeostasis: Cytoprotection and cell death. J. Therm. Biol. 2020, 91, 102615. [Google Scholar] [CrossRef]
  69. Fe-Goncalves, L.M.; Araujo, J.D.A.; dos Santos, C.H.D.; Val, A.L.; de Almeida-Val, V.M.F. How will farmed populations of freshwater fish deal with the extreme climate scenario in 2100? Transcriptional responses of Colossoma macropomum from two Brazilian climate regions. J. Therm. Biol. 2020, 89, 102487. [Google Scholar] [CrossRef]
  70. Mininni, A.N.; 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 Genom. 2014, 15, 765. [Google Scholar] [CrossRef] [Green Version]
  71. Feder, M.E.; Hofmann, G.E. Heat-shock proteins, molecular chaperones and the stress response: Evolutionary and ecological physiology. Annu. Rev. Physiol. 1999, 61, 243–282. [Google Scholar] [CrossRef] [Green Version]
  72. Dammark, K.B.; Ferchaud, A.L.; Hansen, M.M.; Sorensen, J.G. Heat tolerance and gene expression responses to heat stress in threespine sticklebacks from ecologically divergent environments. J. Therm. Biol. 2018, 75, 88–96. [Google Scholar] [CrossRef]
  73. Corey, E.; Linnansaari, T.; Cunjak, R.A.; Currie, S. Physiological effects of environmentally relevant, multi-day thermal stress on wild juvenile Atlantic salmon (Salmo salar). Conserv. Physiol. 2017, 5, cox014. [Google Scholar] [CrossRef]
  74. Hori, T.S.; Gamperl, A.K.; Afonso, L.O.; Johnson, S.C.; Hubert, S.; Kimball, J.; Bowman, S.; Rise, M.L. Heat-shock responsive genes identified and validated in Atlantic cod (Gadus morhua) liver, head kidney and skeletal muscle using genomic techniques. BMC Genom. 2010, 11, 72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Liu, B.; Xu, P.; Brown, P.B.; Xie, J.; Ge, X.P.; Miao, L.H.; Zhou, Q.L.; Ren, M.C.; Pan, L.K. The effect of hyperthermia on liver histology, oxidative stress and disease resistance of the Wuchang bream, Megalobrama amblycephala. Fish Shellfish Immunol. 2016, 52, 317–324. [Google Scholar] [CrossRef] [PubMed]
  76. Nakai, A.; Satoh, M.; Hirayoshi, K.; Nagata, K. Involvement of the stress protein HSP47 in procollagen processing in the endoplasmic reticulum. J. Cell Biol. 1992, 117, 903–914. [Google Scholar] [CrossRef] [PubMed]
  77. Pandey, A.; Rajesh, M.; Baral, P.; Sarma, D.; Tripathi, P.H.; Akhtar, M.S.; Ciji, A.; Dubey, M.K.; Pande, V.; Sharma, P.; et al. Concurrent changes in thermal tolerance thresholds and cellular heat stress response reveals novel molecular signatures and markers of high temperature acclimation in rainbow trout. J. Therm. Biol. 2021, 102, 103124. [Google Scholar] [CrossRef]
  78. Long, Y.; Li, L.; Li, Q.; He, X.; Cui, Z. Transcriptomic characterization of temperature stress responses in larval zebrafish. PLoS ONE 2012, 7, e37209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Biology 11 00990 g001 550
Figure 1. Plasmatic activity of alanine aminotransferase (ALT), aspartate aminotransferase (AST), and alkaline phosphatase (AP) in control and high-temperature stress groups of G. chilensis. Levels of ALT (A), AST (B), and AP (C) activity are expressed in IU/L. Bars represent the mean ± SEM. Significant differences between the control and stress groups are indicated by asterisks; * (p < 0.05) and ** (p < 0.01).
Figure 1. Plasmatic activity of alanine aminotransferase (ALT), aspartate aminotransferase (AST), and alkaline phosphatase (AP) in control and high-temperature stress groups of G. chilensis. Levels of ALT (A), AST (B), and AP (C) activity are expressed in IU/L. Bars represent the mean ± SEM. Significant differences between the control and stress groups are indicated by asterisks; * (p < 0.05) and ** (p < 0.01).
Biology 11 00990 g001
Biology 11 00990 g002 550
Figure 2. Evaluation of oxidative damage in the liver of G. chilensis in control and high-temperature stress groups. Level of DNA damage in terms of AP sites (A), protein carbonylation (B), and lipid peroxidation in terms of HNE adducts (C) in control and stress groups. Bars represent the mean ± SEM. Significant differences between the control and stress groups are indicated by asterisks; * (p < 0.05), ** (p < 0.01) and **** (p < 0.0001).
Figure 2. Evaluation of oxidative damage in the liver of G. chilensis in control and high-temperature stress groups. Level of DNA damage in terms of AP sites (A), protein carbonylation (B), and lipid peroxidation in terms of HNE adducts (C) in control and stress groups. Bars represent the mean ± SEM. Significant differences between the control and stress groups are indicated by asterisks; * (p < 0.05), ** (p < 0.01) and **** (p < 0.0001).
Biology 11 00990 g002
Biology 11 00990 g003 550
Figure 3. Differentially expressed transcripts in response to high-temperature stress in G. chilensis. The heatmap presents the differentially expressed transcripts between control (C) and stress groups (HT). The numbers 1 to 6 indicate the sample library of control and stress fish.
Figure 3. Differentially expressed transcripts in response to high-temperature stress in G. chilensis. The heatmap presents the differentially expressed transcripts between control (C) and stress groups (HT). The numbers 1 to 6 indicate the sample library of control and stress fish.
Biology 11 00990 g003
Biology 11 00990 g004 550
Figure 4. Illustration of transcriptional changes in liver of G. chilensis in response to high-temperature stress using modified KEGG pathway map of protein processing in the endoplasmic reticulum. The red rectangles indicate down-regulated transcripts, while green rectangles indicates up-regulated transcripts between the control and high-temperature groups.
Figure 4. Illustration of transcriptional changes in liver of G. chilensis in response to high-temperature stress using modified KEGG pathway map of protein processing in the endoplasmic reticulum. The red rectangles indicate down-regulated transcripts, while green rectangles indicates up-regulated transcripts between the control and high-temperature groups.
Biology 11 00990 g004
Biology 11 00990 g005 550
Figure 5. qPCR validation of selected differentially expressed transcripts in the liver of G. chilensis in response to high-temperature stress. The transcript expression levels were normalized with the geometric means of actb and taf12. The differential expression levels according to qPCR (black bars) and RNA-seq (gray bars) for these selected genes are expressed as log2 fold changes. The log2FC represents the expression change in the stress group compared with the control group. Results are expressed as the mean ± standard error. Significant differences in the validated qPCR data between control and stress groups are indicated by asterisks in the log2 Fold qPCR bars; (* p-value < 0.05).
Figure 5. qPCR validation of selected differentially expressed transcripts in the liver of G. chilensis in response to high-temperature stress. The transcript expression levels were normalized with the geometric means of actb and taf12. The differential expression levels according to qPCR (black bars) and RNA-seq (gray bars) for these selected genes are expressed as log2 fold changes. The log2FC represents the expression change in the stress group compared with the control group. Results are expressed as the mean ± standard error. Significant differences in the validated qPCR data between control and stress groups are indicated by asterisks in the log2 Fold qPCR bars; (* p-value < 0.05).
Biology 11 00990 g005
Table
Table 1. Primer sequence, amplicon size, Tm and efficiency of validated genes by qPCR.
Table 1. Primer sequence, amplicon size, Tm and efficiency of validated genes by qPCR.
AbbreviationGene NameForward SequenceReverse SequenceAmplicon Size (bp)TmEfficiency (%)
taf12Transcription initiation factor TFIID subunit 12GATCTGTAACGACGACGAAGAACAAATCAGAGGGACGTCATGTA9262101
actbBeta actinTGTCCCTGTATGCTTCTGGTCCCCTCTCAGTCAGAATCTTCAT17262104
hsp60Heat shock protein 60GACGGTTCCAATCTCTACATCTCCGCTCTCCAAACCAGTTACA866299
hsp70Heat shock protein 70AAGATCAGCGACGACGATAAGCTGGTGCTCATACTCCTCTTTC1056295
ddit4DNA damage-inducible transcript 4GGGAATGAGGAGTTTGGTACATGAAGGAAGTGGTGGACCTTATT886296
gpx7Glutathione peroxidase 7-likeTCTCCTTCCCTCTGTTCAGTAAGAAATTCCAGTCGGGCTCTT1046299
leptinLeptinCGAAGAGACTTCCTGCTTCACCTGATGATCTGGGTGGACTTTC11062104
msh2DNA mismatch repair Msh2-likeGCCCGTTCCCAGATATTTGATCGACCGCAATGACTACTACAC1006298
msh3DNA mismatch repair Msh3-likeCGACTTCTTCAGGGACTTTGGTGGCTCTCTGAGTGTCTGT7862104
clql1Complement C1q 2-likeGATGTTTGTGGCGACGTATTTGGTTGCTTTCTCAGCCTCTGTA9962104
ccl20C-C motif chemokine 20-likeCAGCCGTGTGTTAGGGAATACAGTTGTCTCGTGTCTCTCTATC1236298
atg12Ubiquitin ATG12-likeGCCCTCACCAGATCAAGAAGAGAGAGTCAGAGTGGAGTTAGAG13362102
atg4bCysteine protease ATG4B-likeATCTGGGCGATCTGATGAATGCGGAGGGCAGAAACAAAGA9662102
casp3Caspase 3GCTCCAATTCTTTCCCGTATTTCAGATTTCCTCTACGCCTACTC12361103
c3Complement C3-likeCTGCTTCTGGTGACCTGTTTACTTCGTGTCCTCTCCATCTTTC9962103
Table
Table 2. Sequencing and mapping statistics of liver libraries from control and thermal stress groups of G. chilensis.
Table 2. Sequencing and mapping statistics of liver libraries from control and thermal stress groups of G. chilensis.
Experimental SamplesNumber of ReadsAverage Length Number of ReadsNumber of Reads after TrimmingAverage Length after TrimmingPercentage of Mapped Reads
Control 151,220,46810151,103,24895.085.1
Control 243,384,32810143,371,54794.986.0
Control 350,445,70810150,425,18394.985.4
Control 470,505,07410170,462,85894.984.7
Control 565,838,80210165,800,59294.985.7
Control 657,706,00410157,677,64694.985.3
Stress 161,733,30810161,711,17594.986.8
Stress 250,630,80610150,612,51894.984.1
Stress 358,086,70410158,053,69994.984.5
Stress 465,071,75610165,045,36594.985.5
Stress 562,370,66410162,343,21294.884.0
Stress 659,632,64410159,610,42194.985.8
Average58,052,18910158,018,12294.985.2
Table
Table 3. Enriched GO terms represented in up-regulated transcripts in response to high-temperature stress in the liver of G. chilensis.
Table 3. Enriched GO terms represented in up-regulated transcripts in response to high-temperature stress in the liver of G. chilensis.
GO NameGO CategoryGO IDFDRN° of Transcripts
Protein foldingBPGO:00064572.96 × 10−1128
Establishment of protein localizationBPGO:00451846.62 × 10−330
Protein transportBPGO:00150311.36 × 10−229
Nitrogen compound transportBPGO:00717051.65 × 10−233
Protein localization to endoplasmic reticulumBPGO:00709723.41 × 10−26
Macromolecule localizationBPGO:00330363.41 × 10−236
TransportBPGO:00068103.41 × 10−296
Cellular protein localizationBPGO:00346133.41 × 10−225
Organic substance transportBPGO:00717023.41 × 10−238
Protein localizationBPGO:00081043.41 × 10−231
Response to heatBPGO:00094083.41 × 10−25
Cellular macromolecule localizationBPGO:00707273.41 × 10−225
Establishment of localizationBPGO:00512343.80 × 10−296
Golgi organizationBPGO:00070303.80 × 10−24
Protein retention in ER lumenBPGO:00066214.09 × 10−23
Maintenance of protein localization in organelleBPGO:00725954.09 × 10−23
Maintenance of protein localization in endoplasmic reticulumBPGO:00354374.09 × 10−23
Intracellular transportBPGO:00469074.09 × 10−228
Endoplasmic reticulumCCGO:00057833.20 × 10−633
CytoplasmCCGO:00057379.42 × 10−6133
Nuclear outer membrane-endoplasmic reticulum membrane networkCCGO:00421751.31 × 10−524
Endoplasmic reticulum sub-compartmentCCGO:00988273.68 × 10−523
Endoplasmic reticulum membraneCCGO:00057893.68 × 10−523
Endomembrane systemCCGO:00125056.46 × 10−550
Organelle sub-compartmentCCGO:00319847.63 × 10−429
Organelle membraneCCGO:00310902.14 × 10−344
SarcomereCCGO:00300173.00 × 10−27
MyofibrilCCGO:00300163.00 × 10−27
Coated membraneCCGO:00484753.41 × 10−210
MembraneCCGO:00160203.41 × 10−2225
Membrane coatCCGO:00301173.41 × 10−210
Coated vesicle membraneCCGO:00306623.41 × 10−27
Cytoplasmic vesicle membraneCCGO:00306593.41 × 10−28
Contractile fiberCCGO:00432923.41 × 10−27
Unfolded protein bindingMFGO:00510821.67 × 10−919
Signal sequence bindingMFGO:00050482.70 × 10−35
ER retention sequence bindingMFGO:00469234.09 × 10−23
Phosphofructokinase activityMFGO:00084434.69 × 10−25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dettleff, P.; Zuloaga, R.; Fuentes, M.; Gonzalez, P.; Aedo, J.; Estrada, J.M.; Molina, A.; Valdés, J.A. High-Temperature Stress Effect on the Red Cusk-Eel (Geypterus chilensis) Liver: Transcriptional Modulation and Oxidative Stress Damage. Biology 2022, 11, 990. https://doi.org/10.3390/biology11070990

AMA Style

Dettleff P, Zuloaga R, Fuentes M, Gonzalez P, Aedo J, Estrada JM, Molina A, Valdés JA. High-Temperature Stress Effect on the Red Cusk-Eel (Geypterus chilensis) Liver: Transcriptional Modulation and Oxidative Stress Damage. Biology. 2022; 11(7):990. https://doi.org/10.3390/biology11070990

Chicago/Turabian Style

Dettleff, Phillip, Rodrigo Zuloaga, Marcia Fuentes, Pamela Gonzalez, Jorge Aedo, Juan Manuel Estrada, Alfredo Molina, and Juan Antonio Valdés. 2022. "High-Temperature Stress Effect on the Red Cusk-Eel (Geypterus chilensis) Liver: Transcriptional Modulation and Oxidative Stress Damage" Biology 11, no. 7: 990. https://doi.org/10.3390/biology11070990

APA Style

Dettleff, P., Zuloaga, R., Fuentes, M., Gonzalez, P., Aedo, J., Estrada, J. M., Molina, A., & Valdés, J. A. (2022). High-Temperature Stress Effect on the Red Cusk-Eel (Geypterus chilensis) Liver: Transcriptional Modulation and Oxidative Stress Damage. Biology, 11(7), 990. https://doi.org/10.3390/biology11070990

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop