Skip to main content

Transgenerational exposure to ocean acidification impacts the hepatic transcriptome of European sea bass (Dicentrarchus labrax)


Physiological effects of ocean acidification associated with elevated CO2 concentrations in seawater is the subject of numerous studies in teleost fish. While the short time within-generation impact of ocean acidification (OA) on acid-base exchange and energy metabolism is relatively well described, the effects associated with transgenerational exposure to OA are much less known. Yet, the impacts of OA can vary in time with the potential for acclimation or adaptation of a species. Previous studies in our lab demonstrated that transgenerational exposure to OA had extensive effects on the transcriptome of the olfactory epithelium of European sea bass (Dicentrarchus labrax), especially on genes related to ion balance, energy metabolism, immune system, synaptic plasticity, neuron excitability and wiring. In the present study, we complete the previous work by investigating the effect of transgenerational exposure to OA on the hepatic transcriptome of European sea bass. Differential gene expression analysis was performed by RNAseq technology on RNA extracted from the liver of two groups of 18 months F2 juveniles that had been exposed since spawning to the same AO conditions as their parents (F1) to either actual pH or end-of-century predicted pH levels (IPCC RCP8.5), respectively. Here we show that transgenerational exposure to OA significantly impacts the expression of 236 hepatic transcripts including genes mainly involved in inflammatory/immune responses but also in carbohydrate metabolism and cellular homeostasis. Even if this transcriptomic impact is relatively limited compared to what was shown in the olfactory system, this work confirmed that fish transgenerationally exposed to OA exhibit molecular regulation of processes related to metabolism and inflammation. Also, our data expand the up-regulation of a key gene involved in different physiological pathways including calcium homeostasis (i.e. pthr1), which we already observed in the olfactory epithelium, to the liver. Even if our experimental design does not allow to discriminate direct within F2 generation effects from transgenerational plasticity, these results offer the perspective of more functional analyses to determine the potential physiological impact of OA exposure on fish physiology with ecological relevance.

Peer Review reports


Most fish are effective acid-base regulators so that when exposed to ocean acidification (OA), they can regulate extracellular pH. However, pH homeostasis induces elevated PCO2 and HCO3- (hypercapnia) in their extracellular fluids that may have physiological effects [1]. The most documented effects are probably those related to the alteration of GABAergic system, due to an inversion of the HCO3- channels across neuronal membranes, and related consequences on sensory system and behaviour [2, 3]. The large amount of energy allocated for the ion regulation process to maintain the acid-base balance in the reduced-pH environment is not associated with a clear modification of the standard and maximal metabolic rates in most of the fish species studied [4]. However, several studies indicated OA-induced regulation in energetic requirement on isolated tissues suggesting trade-off in aerobic and anaerobic metabolism within an individual [1, 5]. For example, hypercapnia caused a shift from aerobic to anaerobic pathways of carbohydrate oxidation in the muscular tissue of Gilthead seabream (Sparus aurata) associated with the tampering of hepatic oxygen consumption and key mitochondrial enzyme activities in Atlantic Cod (Gadus morhua) and Marbled rockcod (Notothenia rossii) [5,6,7,8]. The regulation of energy metabolism is also reflected in the regulation of plasma or tissue metabolite levels, as well as in the regulation of the activity and/or the gene expression of key enzymes of aerobic and anaerobic metabolic pathways. Araújo and collaborators (2018) showed that hypercapnia-induced metabolic reprogramming in the muscle of Sparus aurata involved the regulation of the expression of glycolysis related enzymes such as triosephosphate isomerase B (tpisb), glyceraldehyde-3-phosphate dehydrogenase (g3p) and glycogen debranching enzyme (gde) as well as enzymes involved in lipid metabolism (apo1). Such metabolic regulation can differ among tissues depending on their function according to OA-induced energy reallocation from fitness-related traits (e.g. growth, reproduction, immune defence) [9, 10]. Moreover, OA-induced metabolic effects can vary depending on the duration of OA exposure and on the capacity of fish for OA acclimation [11]. Thus, OA induced adverse effects on metabolic rates in juvenile anemonefish, Amphiprion melanopu, but these effects did not occur when their parents experienced the simimar OA, revealing transgenerational acclimation [12].

It has been demonstrated that the environment in which teleost exist may generate long-term changes on its immune system [13]. While studies on shellfish demonstrated that hypercapnia has a significant negative impact on the immune system (e.g. hemocyte reactive oxygen species production, hemocyte count, lysosomal content), the scarce information available for fish indicated several OA effects on physiological parameters including immune function which, may be beneficial for the tolerance of fish to pathogens [14, 15]. For example, OA stimulated the innate immunity by increasing activity of lysozyme and complement system in Atlantic halibut (Hippoglossus hippoglossus) [16]. Proteomic analysis in Sparus aurata indicated that hypercapnia induced a regulation in the expression of factors related to cell adhesion molecule production that is relevant for immune system [7]. More recently, transcriptomic effects on the olfactory epithelium of European seabass (Dicentrarchus labrax), including stimulation of genes involved in innate antiviral immunity (pathogen recognition receptors and interferon-stimulated genes) were observed in F2 fish transgenerationally exposed to OA [17]. These regulations were associated with a better resistance in fish challenged with betanodavirus. The modulation of the expression of genes involved in immune response was associated with the extensive regulation of genes involved in odour transduction, ionic homeostasis and in energy metabolism. This long-term OA-induced complex regulation among biological processes within the olfactory epithelium of European sea bass is likely a part of a broader adaptive physiological response that would require further attention in other tissues.

Data regarding the effects of acidification on the hepatic biological processes of fish are very scarce. Since the liver is the primary organ for metabolism, is involved in immune regulation [18] and is sensitive to a great variety of environmental factors [19, 20], we hypothesized that hepatic transcriptome could also be regulated by transgenerational exposure to OA which could have physiological consequences for fish. While the liver has an extraordinary capacity to regenerate upon various acute injuries, chronic stress such as long-term OA exposure may induce hepatic tissue disorders such as fibrosis that may interfere with normal tissue function [21, 22]. In the present study, we characterized the effects of OA on biological processes by targeting differentially expressed genes with a special focus on genes related to immune responses, metabolism and cellular homeostasis that have been revealed as regulated at transcriptomic level in metazoans [23]. This study was carried out on the same batch of animals that were subjected to transcriptomic analysis in the olfactory rosette, allowing comparison of the global transcriptomic response induced by OA between the two tissues [17].

Materials and methods

Animal husbandry and experimental setup

Experiments were conducted at the French Research Institute for Exploitation of the Sea (IFREMER) in Plouzané, within the facilities of the PHYTNESS Laboratory (Agreement number: B29-212-05) following the European Commission recommendation 2007/526/EC and Directive 2010/63/EU for the accommodation and care of animals used for experimental and other scientific purposes. The experiment of OA conditioning was the subject of a specific authorization issued by a French Ethics Committee for animal testing [CEEA – 074: Comité d’éthique finistérien en expérimentation animale (CEFEA): Authorization APAFIS #2,018,032,209,421,223].

Animals used in this study were previously described [17]. Livers were sampled from 18 months-old F2 juveniles of European sea bass exposed from larval stage to control conditions (~ pH8.0) or to OA conditions (~ pH7.6) (Fig. 1). F2 fish were obtained from an in vitro fertilisation of 4 years old F1 parents (20 males X 6 females) exposed from larval to adult stages to control conditions (~ pH8.0) or to OA conditions (~ pH7.6). F2 juveniles were reared in water at the same pH as the brood stock F1 they originated from. No gene expression was performed on the immature gonads of the F2 fish and the absence of visual dimorphism did not allow us to sex them. Rearing conditions during larval and juvenile stages were similar to those described in [24]. 201 fish per treatment were distributed in triplicates in culture tanks (400 L, six tanks = three replicas per treatment) that were part of an open-circuit system. To guarantee high quality, seawater pumped 500 m off the coastline at a depth of 20 m passed through a sand filter, a tungsten heater, a degassing column packed with plastic rings, a 2-µm filter membrane, and a UV lamp. Seawater for the control treatment was then poured into each of the three replicas tanks. For each triplicate per treatment, temperature and salinity followed seasonality of the Bay of Brest. For the OA condition, CO2 was injected at constant flow in a header tank equipped with a degassing CO2 column to favour mixing and adjusted by a flow-control unit. From there, seawater flow poured hydrostatically into each of the 3 replicate tanks. pH and temperature were measured daily with a pH meter (WTW 3110; Xylem Analytics Germany, Weilheim, Germany, NBS scale). Total alkalinity was measured once a week following the adapted protocol of [25]. Summary data containing physico-chemical parameters of the rearing seawater is available in the SEANOE repository

Fish from both conditions were fed ad libitum and no significant difference was observed in the mean body weights between the two groups [17].

Fig. 1
figure 1

Experimental design summary. Rearing times applied on European sea bass parental linage (F1) and their offspring (F2) exposed to either the control (pH8.0, blue) or the acidified (pH7.6, orange) treatment

RNA extraction

Fish sampling was spread over two days during the same time period (9:00–11:00 a.m.) to avoid circadian variations. Fish from both groups were sampled each day to avoid the overlapping of group and day factors. Fish were fasted for 24 h before sampling. Fish were anesthetized (20 mg L-1), and then euthanized with a lethal dose (200 mg L-1) of tricaine methane sulfonate 222 (MS222, Pharmaq, Fordingbridge, Hampshire, UK). Full livers were collected from 15 individuals per group (5 fish per triplicate), weighed then flash frozen in liquid nitrogen and preserved at -80 °C. Total RNA was extracted from 7 individual powdered livers (ground under liquid nitrogen) per group (2 or 3 fish per triplicate). RNA extraction was performed using Extract-All reagent (Eurobio, Courtaboeuf, Essonne, France) combined with Nucleospin RNA column according to the manufacturer’s instructions (Macherey–Nagel, Düren, Germany). RNA concentration and purity were verified (260/280 ratio > 2) using an ND-1000 NanoDrop® spectrophotometer (Thermo Scientific Inc., Waltham, MA, USA) and the integrity of RNA was verified using an Agilent Bioanalyzer 2100 (Agilent Technologies Inc., Santa Clara, CA, USA). All samples showed an RNA integrity (RIN) score > 8. RNA samples were stored at − 80 °C for further RNA sequencing.

RNA-Sequencing analysis

The 14 total RNA samples extracted from liver (7 per group) were sent to the sequencing platform GenomiX (MGX, Montpellier, France) for transcriptome analysis through RNA sequencing (RNA-Seq). RNA-seq libraries were prepared using the Truseq stranded mRNA sample prep kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Briefly, polyadenylated RNAs were selected from 1 µg of total RNA using oligo-dT magnetic beads. The polyA + RNAs were fragmented using divalent cations at elevated temperature and reverse transcribed using random hexamers, Super Script II (Thermo Fisher Scientific, Carlsbad, CA) and actinomycin D. Deoxy-TTP was replaced by dUTP during the second strand synthesis to prevent its amplification by PCR. Double stranded cDNAs were adenylated at their 3’ ends and ligated to Illumina’s adapters containing unique dual indexes (UDI). Ligated cDNAs were PCR amplified for 15 cycles and the PCR products were purified using AMPure XP Beads (Beckman Coulter Genomics, Brea, CA, USA). The size distribution of the resulting libraries was monitored using a Fragment Analyzer (Agilent Technologies, Santa Clara, CA, USA) and the libraries were quantified using the KAPA Library quantification kit (Roche, Basel, Switzerland).

The libraries were denatured with NaOH, neutralized with Tris-HCl, and diluted to 150 pM. Clustering and sequencing were performed on a NovaSeq 6000 (Illumina, San Diego, CA, USA) using the single read 100 nt protocol on 2 lanes of a flow cell SP. Image analyses and base calling were performed using the NovaSeq Control Software and the Real-Time Analysis component (Illumina). Demultiplexing was performed using Illumina’s conversion software (bcl2fastq The quality of the raw data was assessed using FastQC (v 0.11.9) from the Babraham Institute, multiQC v1.9 [26] and the Illumina software SAV (Sequencing Analysis Viewer). FastqScreen (v 0.15.1) was used to identify potential contamination.

Transcriptomic and Gene Ontology (GO) analysis

Raw reads were processed using Trim Galore v0.6.7 [27] to perform quality trimming (Q = 15) and to remove the T overhang at 5’ end of reads, according to Illumina recommendations for stranded mRNA Prep protocols [28].

Trimmed reads were mapped to the European sea bass genome v1 [29] using STAR v2.7.9a [30] with standard parameters. STAR provided matrices of raw counts which were used to perform differential analysis.

To minimize the false-positive rate, the count matrices were filtered for low expressed transcripts (minimum counts > 10 in at least half of the samples). Differentially expressed genes (DEGs) between acidified and control samples were identified using DESeq2 v2.1.34 [31] with R v4.1.2 ( and a collection of graphics and statistics packages : rafalib [32] ; pheatmap [33]; Sartools [34] ; gplots [35] and ggplot2 [36]. DESeq2 method internally corrects for library size and uses negative binomial generalized linear models to test for differential expression. DEGs were annotated using available annotations published with GCA_000689215.1 genome release (Tine et al., 2014).

Statistical analysis

Liver mass and HSI

Shapiro-Wilk normality test were performed to ensure the normality of the data related to liver mass and hepatosomatic index. Data were log transformed to allow analysis of variance (ANOVA). Fisher’s test was used to test if the variances were equal. Two-way ANOVA was performed to determine the potential effect of pH condition and sampling day on the liver mass and the hepatosomatic index, with significant p-value fixed at 0.05.

RNAseq: In this study the statistical design considered two factors: pH condition and sampling day. The statistical model was built using ‘counts ~ sampling_day + group’ design formula, where the qualitative variable ‘group’ indicates sample group (acidified/control). Fold changes were adjusted using Bayesian shrinkage estimators for effect sizes with apeglm package [37]. All features with adjusted p-value smaller than 0.01 (Benjamini-Hochberg method) were reported as differentially expressed (DEGs). Over- and under-representation of Gene Ontology (GO) terms associated with DEGs were tested using goatools v0.8.12 [38]. Goatools process is based on Fisher’s exact test and GO database from geneontology (go-basic.obo release 2017-07-07,


The two way ANOVAs did not reveal any significant effect of pH condition and sampling day on the mean liver mass and hepatosomatic index of the 18-months-old juveniles (Table 1).

Among the 16.816 genes analysed, we detected 236 genes (1.4% of total) as differentially expressed in the liver of acidified (A) vs. control (C) fish groups with p-value ≤ 0.01, 111 and 125 genes exhibiting higher expression levels in the C and A groups, respectively (additional Tables 1 and 2). GO enrichment analysis did not reveal either biological process, molecular pathway or cell component over-represented among the differentially expressed genes. A set of 13 genes related to carbohydrate metabolism was differentially expressed between the A and C groups, 4 genes (akr1a1, enosf1, cpt1a and fggy) being up-regulated, whereas 9 genes exhibited down regulation (pfkfb4, irs2, ppp1cc, g6pc2, ppp1r3c, mpi, adpgk, uxs1 and pck1) in the A group (Table 2).

We found 15 differentially expressed genes related to cellular homeostasis. Among the down-regulated genes in the A group we found slc12a9, sox4 and irs2 while pth1r, txnrd1 while sodc were found up-regulated (Table 3). 46 of the 236 differentially expressed genes are involved in a biological process related to immune system process and/or inflammation (innate immune response, regulation of immune system process, inflammatory response), 17 and 29 genes being up-regulated in the A and C groups, respectively (Table 4). Several nlcr3-like genes exhibited significantly differential expression, 1 and 5 being up and down-regulated in the A group, respectively.

8 and 11 genes involved in cell adhesion were down- and up-regulated in group A, respectively (Table 5). Among them, fibulin-2 gene was highly down-regulated (x0.11), while plexin C-1 was highly up-regulated (x14).

Among the 236 genes differentially expressed, 90 (38% of the differentially expressed genes) were also previously found regulated in the olfactory epithelium of fish from same experimental group [17] and 58 were regulated in the same way. Among them, pvrl1 and pth1r genes exhibited acidification-induced higher expression in both tissues, while the opposite was observed for irs2, ppp1r3c and ppp1c genes (additional Table 3).

Only 20 genes involved in various biological processes were shown differentially expressed between the two sampling days (additional Table 4). No interaction between sampling day and pH condition factors was found for any gene.

Table 1 ANOVA testing for differences in Liver mass Hepatosomatic Index (HSI), in response to pH condition and day of sampling (fixed factors) and interaction effects. Means ± SD in brackets
Table 2 Genes involved in carbohydrate metabolism regulated (p-value ≤ 0.01) by OA in the liver of D. labrax
Table 3 Genes involved in the process of cellular homeostasis regulated (p-value ≤ 0.01) by OA in the liver of D. labrax
Table 4 Genes involved in the immune system and inflammation related processes regulated (p-value ≤ 0.01) by OA in the liver of D. labrax
Table 5 Genes involved in the process of cell adhesion regulated (p-value ≤ 0.01) by OA in the liver of D. labrax


This study reports the effects of transgenerational OA exposure on the hepatic transcriptome of European sea bass. We assume that the experimental protocol used in the present study did not allow to investigate the transgenerational plasticity to OA since hepatic transcriptome was not parallelly analysed in F2 juveniles whose parents had not been previously exposed to OA. Some studies revealed positive transgenerational acclimation of fish to OA since OA-induced decrease in the growth rate or survival of juvenile was compensated when parents also experienced similar conditions [12, 39]. In contrary, negative transgenerational effects of OA exposure were observed in terms of survival and body size in marine three-spined stickleback (Gasterosteus aculeatus) when parents were acclimated to the high-CO2 environment [40]. Even if our results do not allow to distinguish direct effect of OA on F2 hepatic transcriptome from transgenerational effects of parents’ previous exposure, they reveal the impact of OA on hepatic transcriptome considering the transgenerational exposure of fish to a realistic scenario of future ocean acidification, which is of ecological relevance.

In order to limit false positive differentially expressed genes, only RNAs with adjusted p-value below 0.01 were retained as differentially expressed, which ensures robustness of our data. Only 20 genes involved in a wide range of biological processes were shown differentially expressed between the sampling days. No relevant scientific information can be derived from the regulation of so few genes and these data will not be discussed.

Our results displayed 236 genes significantly differentially expressed by OA. This number of differentially expressed genes is relatively low compared to the 9112 found differentially expressed with the same adjusted p-value ≤ 0.01 in the olfactory rosette of the same fish batch under the same conditions [17]. This difference can be explained, at least in part, by the fact that the olfactory rosette is in direct contact with the surrounding environment, therefore it continuously activates mechanisms for regulating ion exchange that interfere with other biological processes related to metabolism and neural activity. Accordingly, we only found one gene with potential function in inorganic ion homeostasis (slc12a9) regulated by OA in the liver [41]. This data confirms the efficiency of fish to compensate for a hypercapnic acid-base disturbance through regulation of ion exchange occurring chiefly in water exposed epithelia [1, 42]. In the same way, while we cannot rule out post transcriptional regulation of genes and associated differential activity of metabolic pathways in the liver, our data indicated limited transcriptomic effects of transgenerational exposure to OA on the hepatic metabolic pathways in European sea bass. Since the liver is a metabolic overachiever organ which governs energy balances in the body and that tissues were sampled in unstimulated and fasted individuals, this data suggests limited long-term effects of OA on the global energy metabolism of resting European sea bass. This result agrees with a previous study indicating that 1.5 years exposure to OA did not affect fish standard metabolic rate (SMR) in European sea bass [43]. Moreover, in this latter study, authors revealed that OA induced higher aerobic capacities (maximal metabolite rates) and concluded by asking about the potential impact of OA on anaerobic metabolic capacity. Anaerobic capacity of fish depends, among other, on their ability to use glucose as substrate through anaerobic glycolysis. Previous studies performed in gilthead seabream revealed that the intermediary metabolism, especially the balance between aerobic and anaerobic metabolism pathways, was reorganized in the liver of fish exposed five weeks to hypercapnia in order to maintain systemic ion homeostasis (Ruiz-Jarabo et al., 2021). In the present study, 13 genes (5% of the differentially expressed genes) involved in carbohydrate metabolism were found differentially expressed. Among them, some genes promoting glycolysis (pfkfb4, irs2 and adpgk) and regulating glycogen metabolism (ppp1cc, ppp1r3) were found down-regulated, suggesting an OA-induced long term discrete fine-tuning regulation of glucose homeostasis in the liver [44,45,46,47]. According to that, we found the gene coding for the PTHrP receptor (pth1r) significantly up-regulated in the liver of fish exposed to hypercapnia. Interestingly, PTHrP has a role in the regulation of hepatic gluconeogenesis in European sea bass where it is supposed to promote hepatic glucose export to peripheral tissues [48]. It is noteworthy that part of the genes (e.g. pth1r, ppp1cc, ppp1r3) related to carbohydrate metabolism we found regulated in the liver was similarly regulated in the olfactory rosette, suggesting comparable OA-induced regulation of cellular pathways among different tissues of fish. In line with the potential regulation of glucose export from the liver, we observed an up-regulation (p value < 0.05, additional Table 1) of the mRNA coding for glucose transporter 2 (glut2) in the fish exposed to OA [49]. The liver is a relevant tissue for glycogen storage and has a major gluconeogenic function regulating blood glucose levels. Therefore the up and down-regulation of carbohydrate metabolic pathway genes we observed may be associated with regulation of glycemia and anaerobic capacity. Since the hepatic glycogen and plasmatic glucose levels were not analysed in the present study, we are not able to determine whether the transcriptomic regulation was associated with a perturbation of the blood glucose content. Previous study in European sea bass indicated that while higher circulating glucose levels were observed after short term exposure to hypercapnia, no difference was observed after several weeks of acclimation [50]. Potential regulation of anaerobic capacity may be revealed through analysis of fish performance in environment with O2 concentration below the critical oxygen level (the oxygen threshold below which SMR is no longer sustainable aerobically).

The regulation of pth1r gene that is essential for calcium mobilization, vitellogenesis or cortisol production in other tissues [51] may have consequence on fish physiology including bone mineralization, sexual maturation or stress response. According to that, Di Santo (2019) found changes in mineralization in the skeleton of L. erinacea [52]. Moreover, a recent study performed in European sea bass indicates that OA alters the physiological response to acute stress via the neuroendocrine regulation of the corticotropic axis [53]. Also, recent data revealed that OA regulates the reproductive function both in terms of sexual maturation and gamete quality [54]. Interestingly, we observed in the present study a down-regulation in the liver of fish exposed to OA of the gene coding estradiol receptor alpha that are responsible for the liver production of vitellogenins. We are fully aware that the regulation of the pth1r gene cannot be formally associated with the regulation of these physiological functions in other tissues of fish long-term exposed to OA, but our result offer possibilities for further investigation that also call for a comparison of metabolomes.

The very few genes (i.e. cpt1a) involved in lipid metabolism exhibiting differential expression in the present study confirmed previous data obtained in cod showing the resilience of lipid metabolism for fish when exposed for long time to elevated pCO2 [55].

Our data revealed also several differentially expressed genes involved in cellular homeostasis suggesting regulation in cellular mechanisms for maintaining normal liver function in fish transgenerationally exposed to OA. Especially, two genes involved in redox homeostasis (txnrd1 and sodc) were found up-regulated in the liver of fish exposed to OA which is consistent with the associated increase in akr1a1 gene in yellow catfish, whose expressions has been associated with the regulation of oxidative status [56]. While the present data suggest that transgenerational exposure to OA may be associated to the regulation of molecular mechanisms involved in hepatic antioxidant defence, further analyses including measurements of antioxidant enzyme activities and intra-tissue ROS concentrations would be required to assess the actual oxidative status of the hepatic tissue. That would be of special interest since the production of ROS interacts with numerous biological processes including pathogen destruction by the innate immune system [57, 58] and the OA-induced regulation of ROS production in shellfish has been associated with a regulation of the immune system [14, 15].

Among the differentially expressed genes in the liver, 20% are directly or indirectly involved in the immune or inflammatory responses. These observations revealed an effect of OA in the fish immune mechanisms with the regulation of important pro- and anti-inflammatory mediators. The down-regulation of genes involved in pro-inflammatory response (atrn, hp, il1r1, plaa, pla2g7, serpine1) combined with up-regulation of genes involved in anti-inflammatory response (gba, nfkbiaa, sodc) suggest a fine tune regulation of the balanced pro and anti-inflammatory response. Inflammation is a characteristic feature of the hepatic wound-healing response to injury that may be beneficial in the short term for example by inducing immune responses that lead to the eradication of pathogens, but may be at long term linked to the development of fibrosis. Interestingly, our data revealed the down-regulation of mpi whose expression level is considered in vertebrates, including fish, inversely correlated with liver fibrosis [59, 60]. The liver fibrosis is characterized by the synthesis and altered deposition of extracellular matrix components. Again, the regulation of the mpi gene expression alone is not sufficient to suggest that OA regulates the fibrotic state of the liver, but together with the regulation of genes involved in cell adhesion (e.g. fbln2) that we observed in the present study, this data opens up the prospect of further functional studies to address the potential long-term impact of OA on biological processes related to liver fibrogenesis.

In line with the regulations observed on the expression of genes involved in the inflammatory response, we found several genes involved in the regulation of immune system process down-regulated in the liver of fish exposed to OA (sox4, irs2, tmem64, actb, tisb, actl6a, pck1, pla2g7, serpine1, il1r1, mfap4). Combined with the OA-induced up-regulation of other genes involved in the positive regulation of immune system such as factors involved in complement system (cfhr4), those results confirm previous data observed in the olfactory rosette of the same fish indicating molecular regulation by OA of the immune function [17]. While nrlc3l genes expression are regulated in both the liver and the rosette, we did not find many regulated hepatic genes involved in pathogen recognition receptors (except for pvrl1) and interferon-stimulated genes suggesting OA-induced tissue specific regulation of immune process. Interestingly, previous data obtained in other fish species revealed that OA stimulated the innate immunity by increasing activity of complement system and by regulating the expression of factors related to cell adhesion molecule [7, 16]. It would be interesting to determine whether this present expressional hepatic profile is correlated to difference in resistance to pathogens such as V. anguillarum which affect different organs including liver in European sea bass [61].

In conclusion, the present data revealed limited OA-induced regulation of the hepatic transcriptome following transgenerational exposure to an IPCC RCP8.5 scenario. Even limited, this plasticity is a component of molecular phenotypic acclimation of European sea bass to OA. Given the regulation occurring on the expression of genes involved in the metabolism, inflammatory and immune response observed in the present study, the question remains if these differences in phenotypes are associated with a different physiological status that may impact on the fitness of the animals. Further experiments addressing the physiological performance of the fish (sexual maturation, tolerance to pathogen, fish health status) as well as tissue-specific reaction norms towards the various climate change related drivers would greatly facilitate the ecological interpretation of this plasticity in the context of the global climate crisis.

Data Availability

Summary data containing physico-chemical parameters of the rearing seawater is available in the SEANOE repository Raw reads have been published to ENA, under the accession project number PRJEB59362. Analysis scripts used in this study are publicly available at


  1. Heuer RM, Grosell M. Physiological impacts of elevated carbon dioxide and ocean acidification on fish. Am J Physiol Regul Integr Comp Physiol. 2014;307(9):R1061–1084.

    Article  CAS  PubMed  Google Scholar 

  2. Nilsson GE, Dixson DL, Domenici P, McCormick MI, Sørensen C, Watson S-A, Munday PL. Near-future carbon dioxide levels alter fish behaviour by interfering with neurotransmitter function. Nat Clim Change. 2012;2(3):201–4.

    Article  CAS  Google Scholar 

  3. Schunter C, Ravasi T, Munday PL, Nilsson GE. Neural effects of elevated CO2 in fish may be amplified by a vicious cycle. Conserv Physiol 2019, 7(1).

  4. Esbaugh AJ. Physiological implications of ocean acidification for marine fish: emerging patterns and new insights. J Comp Physiol B Biochem systemic Environ Physiol. 2018;188(1):1–13.

    Article  CAS  Google Scholar 

  5. Stapp LS, Kreiss CM, Pörtner HO, Lannig G. Differential impacts of elevated CO2 and acidosis on the energy budget of gill and liver cells from Atlantic cod, Gadus morhua. Comp Biochem Physiol A: Mol Integr Physiol. 2015;187:160–7.

    Article  CAS  PubMed  Google Scholar 

  6. Strobel A, Leo E, Pörtner HO, Mark FC. Elevated temperature and PCO2 shift metabolic pathways in differentially oxidative tissues of Notothenia rossii. Comp Biochem Physiol B: Biochem Mol Biol. 2013;166(1):48–57.

    Article  CAS  PubMed  Google Scholar 

  7. Araújo JE, Madeira D, Vitorino R, Repolho T, Rosa R, Diniz M. Negative synergistic impacts of ocean warming and acidification on the survival and proteome of the commercial sea bream, Sparus aurata. J Sea Res. 2018;139:50–61.

    Article  Google Scholar 

  8. Ruiz-Jarabo I, Gregório SF, Alves A, Mancera JM, Fuentes J. Ocean acidification compromises energy management in Sparus aurata (Pisces: Teleostei). Comparative biochemistry and physiology part A. Mol Integr Physiol. 2021;256:110911.

    Article  CAS  Google Scholar 

  9. Pan T-CF, Applebaum SL, Manahan DT. Experimental ocean acidification alters the allocation of metabolic energy. Proc Natl Acad Sci. 2015;112(15):4696–701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Pörtner HO. Integrating climate-related stressor effects on marine organisms: unifying principles linking molecule to ecosystem-level changes. Mar Ecol Prog Ser. 2012;470:273–90.

    Article  Google Scholar 

  11. Munday PL. Transgenerational acclimation of fishes to climate change and ocean acidification. F1000prime reports 2014, 6:99.

  12. Miller GM, Watson S-A, Donelson JM, McCormick MI, Munday PL. Parental environment mediates impacts of increased carbon dioxide on a coral reef fish. Nat Clim Change. 2012;2(12):858–61.

    Article  CAS  Google Scholar 

  13. Makrinos DL, Bowden TJ. Natural environmental impacts on teleost immune function. Fish Shellfish Immunol. 2016;53:50–7.

    Article  CAS  PubMed  Google Scholar 

  14. Wang Q, Cao R, Ning X, You L, Mu C, Wang C, Wei L, Cong M, Wu H, Zhao J. Effects of ocean acidification on immune responses of the Pacific oyster Crassostrea gigas. Fish Shellfish Immunol. 2016;49:24–33.

    Article  CAS  PubMed  Google Scholar 

  15. Wu F, Xie Z, Lan Y, Dupont S, Sun M, Cui S, Huang X, Huang W, Liu L, Hu M, et al. Short-term exposure of Mytilus coruscus to decreased pH and salinity change impacts Immune Parameters of their haemocytes. Front Physiol. 2018;9:166.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Bresolin de Souza K, Asker N, Jönsson E, Förlin L, Sturve J. Increased activity of lysozyme and complement system in Atlantic halibut exposed to elevated CO2 at six different temperatures. Mar Environ Res. 2016;122:143–7.

    Article  CAS  PubMed  Google Scholar 

  17. Cohen-Rengifo M, Danion M, Gonzalez AA, Bégout ML, Cormier A, Noël C, Cabon J, Vitré T, Mark FC, Mazurais D. The extensive transgenerational transcriptomic effects of ocean acidification on the olfactory epithelium of a marine fish are associated with a better viral resistance. BMC Genomics. 2022;23(1):448.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Möller A-M, Korytář T, Köllner B, Schmidt-Posthaus H, Segner H. The teleostean liver as an immunological organ: intrahepatic immune cells (IHICs) in healthy and benzo [a] pyrene challenged rainbow trout (Oncorhynchus mykiss). Dev Comp Immunol. 2014;46(2):518–29.

    Article  PubMed  Google Scholar 

  19. Kroon F, Streten C, Harries S. A protocol for identifying suitable biomarkers to assess fish health: a systematic review. PLoS ONE. 2017;12(4):e0174762.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Passantino L, Santamaria N, Zupa R, Pousis C, Garofalo R, Cianciotta A, Jirillo E, Acone F, Corriero A. Liver melanomacrophage centres as indicators of Atlantic bluefin tuna, Thunnus thynnus L. well-being. J Fish Dis. 2014;37(3):241–50.

    Article  CAS  PubMed  Google Scholar 

  21. Hintermann E, Christen U. The Many Roles of Cell Adhesion Molecules in Hepatic Fibrosis. Cells 2019, 8(12).

  22. Jones RE, Foster DS, Hu MS, Longaker MT. Wound healing and fibrosis: current stem cell therapies. Transfusion. 2019;59(S1):884–92.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Strader ME, Wong JM, Hofmann GE. Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey. Front Zool. 2020;17:7.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Cominassi L, Moyano M, Claireaux G, Howald S, Mark FC, Zambonino-Infante JL, Peck MA. Food availability modulates the combined effects of ocean acidification and warming on fish growth. Sci Rep. 2020;10(1):2338.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Strickland JDH, Parsons TR. A practical handbook of seawater analysis. Bull Fisheries Res Board Can. 1968;167:1–311.

    Google Scholar 

  26. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Krueger f. : Babraham bioinformatics - trim galore!. https://www.bioinformaticsbabrahamacuk/projects/trim_galore/ 2012.

  28. Illumina. Best practices for read trimming for illumina stranded mrna and total rna workflows. 2020.

  29. Tine M, Kuhl H, Gagnaire PA, Louro B, Desmarais E, Martins RS, Hecht J, Knaust F, Belkhir K, Klages S, et al. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat Commun. 2014;5:5770.

    Article  CAS  PubMed  Google Scholar 

  30. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  PubMed  Google Scholar 

  31. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Irizarry RA, Love MI. rafalib: Convenience Functions for Routine Data Exploration. 2015.

  33. Kolde R. pheatmap: Pretty Heatmaps. https://CRANR-projectorg/package=pheatmap 2015.

  34. Varet H, Brillet-Guéguen L, Coppée JY, Dillies MA. SARTools: a DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data. PLoS ONE. 2016;11(6):e0157022.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Warnes G, Bolker B, Bonebakker L, Gentleman R, Huber W, Liaw A, Lumley T, Mächler M, Magnusson A, Möller S. gplots: Various R programming tools for plotting data, vol. 2; 2005.

  36. Wickham H. Elegant Graphics for Data Analysis. ggplot2, Springer New York, NY 2009.

  37. Zhu A, Ibrahim JG, Love MI. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. 2018;35(12):2084–92.

    Article  PubMed Central  Google Scholar 

  38. Klopfenstein DV, Zhang L, Pedersen BS, Ramírez F, Warwick Vesztrocy A, Naldi A, Mungall CJ, Yunes JM, Botvinnik O, Weigel M, et al. GOATOOLS: a Python library for Gene Ontology analyses. Sci Rep-Uk. 2018;8(1):10872.

    Article  CAS  Google Scholar 

  39. Murray CS, Malvezzi A, Gobler CJ, Baumann H. Offspring sensitivity to ocean acidification changes seasonally in a coastal marine fish. Mar Ecol Prog Ser. 2014;504:1–11.

    Article  Google Scholar 

  40. Schade FM, Clemmesen C, Mathias Wegner K. Within- and transgenerational effects of ocean acidification on life history of marine three-spined stickleback (Gasterosteus aculeatus). Mar Biol. 2014;161(7):1667–76.

    Article  CAS  Google Scholar 

  41. Gagnon KB, Delpire E. Physiology of SLC12 transporters: lessons from inherited human genetic mutations and genetically engineered mouse knockouts. Am J Physiol Cell Physiol. 2013;304(8):C693–714.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Porteus CS, Hubbard PC, Uren Webster TM, van Aerle R, Canário AVM, Santos EM, Wilson RW. Near-future CO2 levels impair the olfactory system of a marine fish. Nat Clim Change. 2018;8(8):737–43.

    Article  CAS  Google Scholar 

  43. Crespel A, Anttila K, Lelièvre P, Quazuguel P, Le Bayon N, Zambonino-Infante J-L, Chabot D, Claireaux G. Long-term effects of ocean acidification upon energetics and oxygen transport in the european sea bass (Dicentrarchus labrax, Linnaeus). Mar Biol. 2019;166(9):116.

    Article  Google Scholar 

  44. Dong X, Park S, Lin X, Copps K, Yi X, White MF. Irs1 and Irs2 signaling is essential for hepatic glucose homeostasis and systemic growth. J Clin Investig. 2006;116(1):101–14.

    Article  CAS  PubMed  Google Scholar 

  45. Manchenkov T, Pasillas MP, Haddad GG, Imam FB. Novel Genes Critical for Hypoxic Preconditioning in Zebrafish Are Regulators of Insulin and Glucose Metabolism. G3 (Bethesda, Md) 2015, 5(6):1107–1116.

  46. Taniguchi CM, Finger EC, Krieg AJ, Wu C, Diep AN, LaGory EL, Wei K, McGinnis LM, Yuan J, Kuo CJ, et al. Cross-talk between hypoxia and insulin signaling through Phd3 regulates hepatic glucose and lipid metabolism and ameliorates diabetes. Nat Med. 2013;19(10):1325–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Mehta MB, Shewale SV, Sequeira RN, Millar JS, Hand NJ, Rader DJ. Hepatic protein phosphatase 1 regulatory subunit 3B (Ppp1r3b) promotes hepatic glycogen synthesis and thereby regulates fasting energy homeostasis. J Biol Chem. 2017;292(25):10444–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Palma PFS, Bock C, Silva TS, Guerreiro PM, Power DM, Pörtner H-O, Canário AVM. STC1 and PTHrP modify carbohydrate and lipid metabolism in liver of a Teleost Fish. Sci Rep-Uk. 2019;9(1):723.

    Article  Google Scholar 

  49. Thorens B. GLUT2, glucose sensing and glucose homeostasis. Diabetologia. 2015;58(2):221–32.

    Article  CAS  PubMed  Google Scholar 

  50. Petochi T, Di Marco P, Priori A, Finoia MG, Mercatali I, Marino G. Coping strategy and stress response of european sea bass Dicentrarchus labrax to acute and chronic environmental hypercapnia under hyperoxic conditions. Aquaculture. 2011;315(3):312–20.

    Article  Google Scholar 

  51. Abbink W, Flik G. Parathyroid hormone-related protein in teleost fish. Gen Comp Endocrinol. 2007;152(2):243–51.

    Article  CAS  PubMed  Google Scholar 

  52. Di Santo V. Ocean acidification and warming affect skeletal mineralization in a marine fish. Proc Biol Sci. 2019;286(1894):20182187.

    PubMed  PubMed Central  Google Scholar 

  53. Servili A, Lévêque E, Mouchel O, Devergne J, Lebigre C, Roussel S, Mazurais D, Zambonino-Infante JL. Ocean acidification alters the acute stress response of a marine fish. Sci Total Environ. 2023;858(Pt 1):159804.

    Article  CAS  PubMed  Google Scholar 

  54. Nagelkerken I, Alemany T, Anquetin JM, Ferreira CM, Ludwig KE, Sasaki M, Connell SD. Ocean acidification boosts reproduction in fish via indirect effects. PLoS Biol. 2021;19(1):e3001033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Frommel AY, Hermann BT, Michael K, Lucassen M, Clemmesen C, Hanel R, Reusch TBH. Differential gene expression patterns related to lipid metabolism in response to ocean acidification in larvae and juveniles of Atlantic cod. Comp Biochem Physiol A: Mol Integr Physiol. 2020;247:110740.

    Article  CAS  PubMed  Google Scholar 

  56. Yang L, Zheng S, Kong D, Xiang S, Wu J, Wan N, Sun W, Li W. Characterization, expression, and function analysis of AKR1A1 gene from yellow catfish (Tachysurus fulvidraco). Fish Physiol Biochem. 2022;48(2):285–302.

    Article  CAS  PubMed  Google Scholar 

  57. Neumann NF, Stafford JL, Barreda D, Ainsworth AJ, Belosevic M. Antimicrobial mechanisms of fish phagocytes and their role in host defense. Dev Comp Immunol. 2001;25(8–9):807–25.

    Article  CAS  PubMed  Google Scholar 

  58. Hodgkinson JW, Grayfer L, Belosevic M. Biology of Bony Fish Macrophages. Biology. 2015;4(4):881–906.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. DeRossi C, Bambino K, Morrison J, Sakarin I, Villacorta-Martin C, Zhang C, Ellis JL, Fiel MI, Ybanez M, Lee YA, et al. Mannose phosphate isomerase and mannose regulate hepatic stellate cell activation and fibrosis in zebrafish and humans. Hepatology (Baltimore MD). 2019;70(6):2107–22.

    Article  CAS  PubMed  Google Scholar 

  60. Morrison JK, DeRossi C, Alter IL, Nayar S, Giri M, Zhang C, Cho JH, Chu J. Single-cell transcriptomics reveals conserved cell identities and fibrogenic phenotypes in zebrafish and human liver. Hepatol Commun. 2022;6(7):1711–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Crisafi F, Denaro R, Genovese M, Cappello S, Mancuso M, Genovese L. Comparison of 16SrDNA and toxR genes as targets for detection of Vibrio anguillarum in Dicentrarchus labrax kidney and liver. Res Microbiol. 2011;162(3):223–30.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank L. Madec for her technical assistance on molecular biology approach.


This study was supported by the AWI-MARUM-IFREMER AMI Partnership Programme (DEADLY TRIO project), LabexMer (ANR-10LABX-0019, OASYS project), the Ministry of Ecological and Solidarity Transition and the Foundation for Biodiversity Research (Ocean Acidification Program, PACIO project) and the Deutsche Forschungsgemeinschaft, PE 1157/8–1, MA4271/3–1 (the FITNESS project). MGX acknowledges financial support from France Génomique National infrastructure, funded as part of “Investissement d’Avenir” program managed by Agence Nationale pour la Recherche (contract ANR-10-INBS-09).

Author information

Authors and Affiliations



DM and FM conceived the scientific questions and experimental design. DM collected RNA samples and ML-F extracted the RNA and verified RNA quality and integrity. A-AG performed the RNA sequencing. PA performed the transcriptomic and Gene Ontology analysis. AS provided ideas concerning the interpretation of the data. DM led the writing of the manuscript. All authors contributed significantly to the revision of the manuscript and approved it for publication.

Corresponding author

Correspondence to David Mazurais.

Ethics declarations

Ethics approval and consent to participate

Experiments were conducted at the French Research Institute for Exploitation of the Sea (IFREMER) in Plouzané, within the facilities of the PHYTNESS Laboratory (Agreement number: B29-212-05) following the European Commission recommendation 2007/526/EC and Directive 2010/63/EU for the accommodation and care of animals used for experimental and other scientific purposes. The experimental protocol was the subject of a specific authorization issued by a French Ethics Committee for animal testing [CEEA – 074: Comité d’éthique finistérien en expérimentation animale (CEFEA): Authorization APAFIS #2018032209421223] and is re reported in accordance with ARRIVE guidelines ( for the reporting of animal experiments.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Auffret, P., Servili, A., Gonzalez, AA. et al. Transgenerational exposure to ocean acidification impacts the hepatic transcriptome of European sea bass (Dicentrarchus labrax). BMC Genomics 24, 331 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: