# Hepatic Glycerol Metabolism-Related Genes in Carnivorous Rainbow Trout (Oncorhynchus mykiss): Insights Into Molecular Characteristics, Ontogenesis, and Nutritional Regulation

Stéphane Panserat, Elisabeth Plagnes-Juan, Elsa Gazzola, Mariana Palma, Leonardo J Magnoni, Lucie Marandel, Ivan Viegas

► To cite this version:

Stéphane Panserat, Elisabeth Plagnes-Juan, Elsa Gazzola, Mariana Palma, Leonardo J Magnoni, et al.. Hepatic Glycerol Metabolism-Related Genes in Carnivorous Rainbow Trout (Oncorhynchus mykiss): Insights Into Molecular Characteristics, Ontogenesis, and Nutritional Regulation. Frontiers in Physiology, 2020, 11, 10.3389/fphys.2020.00882. [https://doi.org/10.3389/fphys.2020.00882](https://doi.org/10.3389/fphys.2020.00882)

HAL Id: hal-02911015

https://hal.inrae.fr/hal-02911015

Submitted on 3 Aug 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Hepatic Glycerol Metabolism-Related Genes in Carnivorous Rainbow Trout (*Oncorhynchus mykiss*): Insights Into Molecular Characteristics, Ontogenesis, and Nutritional Regulation

Stephane Panserat¹*†, Elisabeth Plagnes-Juan¹†, Elsa Gazzola¹, Mariana Palma², Leonardo J. Magnoni³, Lucie Marandel¹ and Ivan Viegas²,⁴

¹ INRAE, Université de Pau et des Pays de l'Adour, E2S UPPA, NuMéA, Saint-Pée-sur-Nivelle, France, ² Centre for Functional Ecology, Department of Life Sciences, University of Coimbra, Coimbra, Portugal, ³ Interdisciplinary Centre of Marine and Environmental Research (CIIMAR), University of Porto, Matosinhos, Portugal, ⁴ Centre for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal

OPEN ACCESS

**Edited by:**
José Luis Soengas,
University of Vigo, Spain

**Reviewed by:**
Isabel Navarro,
University of Barcelona, Spain
Jan Alexander Mennigen,
University of Ottawa, Canada

*Correspondence:
Stephane Panserat
stephane.panserat@inrae.fr;
stephane.panserat@inra.fr

† These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Aquatic Physiology,
a section of the journal
Frontiers in Physiology

**Received:** 16 May 2020
**Accepted:** 29 June 2020
**Published:** 31 July 2020

**Citation:**
Panserat S, Plagnes-Juan E, Gazzola E, Palma M, Magnoni LJ, Marandel L and Viegas I (2020) Hepatic Glycerol Metabolism-Related Genes in Carnivorous Rainbow Trout (*Oncorhynchus mykiss*): Insights Into Molecular Characteristics, Ontogenesis, and Nutritional Regulation. *Front. Physiol.* 11:882. doi: 10.3389/fphys.2020.00882

Glycerol metabolism in rainbow trout is poorly studied even though it is at the interface between lipid and glucose metabolism. Moreover, glycerol can be an important ingredient in new aquafeed formulation to decrease the catabolism of dietary amino acids. Thus, the present study aimed to characterize for the first time the different genes coding for key enzymes and proteins involved in hepatic glycerol metabolism. From the trout genomes, all the paralogous genes coding for glycerol transport (*aqp9b*), glycerol kinase (*gk2a* and *gk5*), glycerol-3-phosphate phosphatase (*pgp*), and glycerol-3-phosphate dehydrogenase (*gpd1a*, *gpd1b*, and *gpd1c*) were identified. The ontogenesis determined that the capacity to metabolize glycerol begins with the apparition of the liver during the development (stage 22) and are more expressed at the endogenous-exogenous feeding period (stage 35). The postprandial regulation of the expression of these genes in juvenile trout showed that the postprandial peak of expression is between 4 and 24 h after the last meal for many of the genes, demonstrating that glycerol metabolism could be nutritionally regulated at a molecular level. However, surprisingly, no regulation of the mRNA abundance for the glycerol metabolism-related genes by different levels of dietary glycerol (0, 2.5, and 5%) have been detected, showing that hepatic glycerol metabolism is poorly regulated at a molecular level by dietary glycerol in rainbow trout juveniles.

**Keywords:** glycerol metabolism, fish nutrition, genomics, gene expression, rainbow trout, aquaculture

**Abbreviations:** Aqp9, aquaporin 9; *aqp9b*, aquaporin 9b gene; DHAP, dihydroxyacetone phosphate; G3P, glycerol-3-phosphate; gk, glycerol kinase gene; Gk, glycerol kinase; *gpd*, glycerol-3-phosphate dehydrogenase gene; Gpd1, cytosolic glycerol-3-phosphate dehydrogenase 1; *pgp*, glycerol-3-phosphate phosphatase gene; Pgp, glycerol-3-phosphate phosphatase.
# INTRODUCTION

The production of aquafeeds based on fishmeal and fish oils has been a key element in recent aquaculture development (Klinger and Naylor, 2012; Food Agric Organ, 2018). However, the stability of fishmeal and fish oil production and the dramatic increase in their costs, as well as the challenge of improving growth while reducing these production costs, are compromising the sustainable development of aquaculture. Fish diets have evolved to include an increased proportion of plant protein sources and vegetable oils (Gatlin et al., 2007; Naylor et al., 2009; Food Agric Organ, 2018). Unfortunately, no adequate complete replacement for fishmeal and fish oil, with no negative effects on growth performance and fish quality, has yet been found (Panserat et al., 2009; Lazzarotto et al., 2018). Moreover, prices for grains and vegetable oils have increased, forcing the livestock industry to look for alternative feedstuffs. Among them, glycerol, a widely available and inexpensive source of dietary energy issue as a coproduct of biodiesel production (Ayoub and Abdullah, 2012; Quispe et al., 2013; Tan et al., 2013), has been proposed as a dietary component to replace other ingredients included in feeds.

The inclusion of glycerol in aquafeeds has been mainly investigated in the omnivorous fish Nile tilapia (*Oreochromis niloticus*), showing that growth performance, voluntary feed intake, and protein efficiency could be affected depending on the inclusion level (Gonçalves et al., 2015). However, low inclusion levels of glycerol in this species (3–6%) did not alter growth performance, protein efficiency ratio, survival rate, lipid deposition, or hepatosomatic and viscerosomatic indexes compared to fish fed a diet without glycerol (Meurer et al., 2012, 2016; Neu et al., 2013; Moesch et al., 2016). Another study in Nile tilapia has shown that this species is capable of metabolizing dietary crude glycerol into lipids, proteins, and/or carbohydrates, as well as into CO₂, to provide energy (Da Costa et al., 2017). A recent study demonstrated a positive effect of dietary glycerol for acting as a protein-sparing ingredient due to the decrease of amino acid catabolism in the carnivorous European seabass (*Dicentrarchus labrax*) (Rito et al., 2019).

Glycerol-3-phosphate (G3P) lies at the crossroad of glucose, lipid, and energy metabolism in vertebrates, actively participating in glycolysis and gluconeogenesis, lipid synthesis, and G3P electron transfer to mitochondria (Xue et al., 2017). Indeed, regarding glycerol metabolism in fish, it has been suggested that this compound is converted into glucose through the gluconeogenesis pathway, which can provide energy for cellular metabolism (Suarez and Mommsen, 1987; Brisson et al., 2001; Rito et al., 2019). Glycerol can be converted to G3P, which can enter the glycolytic pathway and the Krebs cycle to produce ATP (Pearson et al., 1990; Rito et al., 2019). In addition, glycerol is a precursor for the synthesis of triglycerides and phospholipids (Sheridan, 1988; Xue et al., 2017). Therefore, this simple sugar alcohol is pivotal for several metabolic pathways in fish as it is in other vertebrates.

The metabolic switching between feeding and fasting is central to life and involves tight regulation of metabolism including the one in the liver (Brisson et al., 2001). In the liver, aquaporin 9 (AQP9) facilitates the uptake of glycerol (Jelen et al.,

2011; Lebeck, 2014). Glycerol kinase (Gk) catalyzes the initial phosphorylation of glycerol into G3P (Sriram et al., 2008), which can be used for gluconeogenesis after being reversibly oxidized by cytoplasmic G3P dehydrogenase (Gpd1) (Xue et al., 2017) into dihydroxyacetone phosphate (DHAP), a key intermediate of the glucose metabolism. Recently, Mugabo et al. (2016) identified a new pathway of glycerol metabolism in mammals with the identification of the G3P phosphatase (Pgp) that can directly hydrolyze G3P to glycerol. Indeed, in some fish species living in cold water (such as rainbow smelt *Osmerus mordax*), the conversion of G3P to glycerol by Pgp could be highly effective for their cryoprotection through sustained glycerol production (Ditlecadet and Driezdic, 2014; Driedzic, 2015). In cold-adapted Antarctic fish species, the use of glycerol as a precursor for gluconeogenesis was of low physiological importance (Magnoni et al., 2013) underlining how little is known about the regulation of glycerol metabolism, even under such extreme environments. Many of the enzymes involved in these processes are regulated by nutritional hormones (insulin and glucagon) such as Gk (Westergaard et al., 1998), Aqp9 (Lebeck, 2014), Pgp (Mugabo et al., 2016), and cytosolic Gpd1 (Patsouris et al., 2016). On a molecular level, only one study has been published in fish species (common carp *Cyprinus carpio*), which characterized the glycerol transporter gene, i.e., the one coding for the AQP9 (Dong et al., 2016). Until now, to our knowledge, no data are available concerning the molecular regulation of glycerol metabolism by nutrition in fish species.

Rainbow trout (*Oncorhynchus mykiss*) is a suitable experimental fish model, as this carnivorous species has an extensive and solid body of literature and is considered to be one of the most exploited species for freshwater aquaculture worldwide (Food Agric Organ, 2018). The provision of dietary energy in carnivorous fish derived from compounds other than from protein, such as glycerol, may reduce the use of dietary proteins as an energy source through catabolism (Panserat et al., 2009). Rainbow trout is one of the best models for fish nutrition and metabolism (Panserat et al., 2019). A preliminary study investigating the utilization of free glycerol in rainbow trout has been performed 34 years ago (Menton et al., 1986). Despite reporting some reduction in their feeding response, Menton et al. (1986) did not observe significant alteration in the final body weight and other zootechnical parameters nor in the biochemical composition of the carcass after 12 weeks of feeding. Unlike liver glycogen and plasma protein, blood glucose in rainbow trout provided with 6 and 12% glycerol-supplemented diets was more significant than that in control fish. Curiously, while the observed hyperglycemia at 3 h post-feeding returned to basal level at 18 h post-feeding in the fish supplemented with 6% glycerol, at 12% glycerol, fish maintained significantly higher blood glucose levels. This demonstrates an impairment of their capacity to clear the excess glucose at higher levels of dietary glycerol. Recent data about the effects of glycerol supplementation have been described in rainbow trout using a metabolomics approach (Palma et al., 2019). However, as explained before, the molecular regulation of glycerol metabolism by dietary glycerol for this fish species remains to be investigated. Our hypothesis is that glycerol metabolism is regulated at a molecular level by endogenous
and/or exogenous nutrition. Thus, the objectives of the present study were (1) to characterize the glycerol metabolism-related genes using the rainbow trout genome database; (2) to analyze the ontogenesis of their gene expressions from the oocytes up to the first feeding (alevin stage); and to study the nutritional regulation of these genes (3) in the fasting-refeeding transition and the postprandial kinetic response after the last meal and (4) in continuous feeding with the inclusion of increasing levels of dietary glycerol (0, 2.5, and 5%).

# MATERIALS AND METHODS

## Fish Experimental Design and Samplings

In order to test our hypothesis about the existence of nutritional regulation of glycerol metabolism in rainbow trout, the three following trials were used to study the ontogenesis, the postprandial kinetic, and the nutritional regulation for all the glycerol metabolism-related genes.

### Fish Experimental Diet – Ontogenesis Sampling

Rainbow trout (*O. mykiss*) spawns were fertilized synchronously with neomale sperm and reared in separate tanks at 8°C at the INRA experimental facilities, Lées-Athas, France. Fish were sampled according to spawn origin before fertilization (oocyte) and then during development according to stages which are described by Vernier (1969): oocyte to stage 8, before embryonic genome activation (EGA); stage 10, EGA; stages 12 and 15, epiboly period; stage 22, primitive liver; stage 23, primitive hepatic portal vein; and stages 35–37, hatched embryo/endogenous feeding period. Endogenous feeding alevins (unfed fish) were sampled at 384°D (degree days, Vernier stage 31). Alevins at 616°D (Vernier stage 35) still exhibiting small yolk reserves were fed the first meal at 09:00 h and killed 3 h after the meal (corresponding to the postprandial period in alevins). The experimental feed, prepared in our facilities (Donzacq, France), contains crude protein (58.4% dry matter) and crude fat (17.5% dry matter), and gross energy is 21.7 kJ/g dry matter (for a detailed list of ingredients, see Marandel et al., 2016). Finally, 87-day-old alevins (696°D at 8°C), characterized by complete resorption of the yolk sac (Vernier stage 37) and maintained on a continuous feeding regime of several meals a day, were killed 3 h following the final feeding as described above. Care was taken to maintain the same daily time frame before samplings to avoid potential circadian effects, known to affect metabolism in rainbow trout (Bolliet et al., 2000). Fish were therefore consistently sampled at 12:00 h for each ontogenetic stage investigated in this study, as previously described (Mennigen et al., 2013). Embryos were directly snap-frozen, whereas alevins were killed by terminal anesthetization by bathing in benzocaine before pooling and storage in liquid nitrogen. We used *N* = 3 pools of 30 embryos, each pool coming from one female from stages O to 23 and then *N* = 3 alevins per spawn (in total, *N* = 9) from stages 31 to 37. The samples were stored at −80°C until analyzed.

### Fish Experimental Diet – Postprandial Sampling

Prior to the experiment, rainbow trout (*O. mykiss*) fish had initially been reared in our experimental facilities (INRA, Donzacq, France) at 18°C and fed a commercial diet (Skretting, France; crude protein: 49.8% dry matter, crude fat: 13.8% dry matter, gross energy: 22 kJ/g dry matter). Immediately prior to the experiment, fish were fasted for 48 h, to allow for basal metabolite plasma concentrations to be reached. In trout, these basal metabolite concentrations are typically reached more slowly compared to endothermic mammals due to slower intestinal transit and gastric emptying. Following the fast, fish were fed once *ad libitum* with the commercial diet. Six trout were sampled for each time point, starting with 48 h fasted fish at 0 h and following feeding at 2, 4, 8, 12, 16, and 24 h, as previously described (Mennigen et al., 2012). Immediately following a terminal anesthetization, the liver was dissected and frozen in liquid nitrogen prior to storage at −80°C until analyzed. Gut content of the sampled animals was checked to verify that the fish had effectively consumed the diet.

### Fish Experimental Diet – Nutritional Sampling

Juvenile rainbow trout (*O. mykiss*) were reared in tanks at the Experimental Research Station (Vila Real, Portugal) of the University of Trás-os-Montes e Alto Douro (UTAD) facilities in freshwater according to standard juvenile rearing protocols. At the beginning of the experiment, fish were fasted for 24 h and were then divided into triplicate groups of 25 fishes (225 fish in total, BW 20.22 ± 0.06 g) randomly distributed among nine fiberglass tanks in an open circulation water system (14.6 ± 0.1°C). Three experimental diets [SPAROS LDA, Olhão, Portugal; isoproteic (49% dry matter) and isoenergetic (gross energy: 21 kJ/g dry matter)] were formulated to give a contrast in the glycerol content by supplementing a basal diet (control) with this compound at 2.5 or 5% w/w; Magnoni et al., unpublished data. Glycerol used in the formulation of the diet (refined glycerine) was obtained from rapeseed (Belgosuc, 050008, Beernem, Belgium). The incorporation of cellulose at different levels in the formulation was used to create the dietary contrast for the addition of glycerol. Therefore, cellulose was used as an indigestible filler (Hilton, 1982), which does not affect the utilization of dietary energy or protein in rainbow trout (Bromley and Adkins, 1984; Wiesmann and Pfeffer, 1986). The fish were hand-fed *ad libitum* 6 days/week and twice per day (9:00 and 17:00 h). Fish survival was 99, 100, and 99% for groups assigned to control, 2.5, and 5% glycerol diets, respectively. All the fish were individually weighed at the end of the feeding trial (60 days), and no differences were detected on the final BW between groups (Magnoni et al., unpublished data). Two to three fish per tank were sampled 6 h after the morning feeding (n = 8) and 24 h after the morning feeding (n = 8). Immediately following a terminal anesthetization, the liver was dissected, weighed, and immediately frozen in liquid nitrogen and stored at −80°C until analyzed.

### In silico Analysis

Rainbow trout *aquaporin 9b* (aqp9b), *glycerol kinase (gk)*, *G3P phosphatase (pgp)*, and *G3P dehydrogenase 1* (gpd1) genes
were extracted from the *O. mykiss* genome assembly NCBI database (GCF_002163495.1, June 2017) using the zebrafish (*Danio rerio*) orthologous sequence as a query sequence in the BLAST tool (rainbow trout accession numbers are given in **Table 1**). Phylogenetic analysis (**Figure 1**) was performed using the Molecular Evolutionary Genetics Analysis (MEGA) package v6 software (Tamura et al., 2013). To do so, sequences from species other than zebrafish and trout were collected from the Ensembl Genome database (Ensembl release 99, January 2020)<fnref n="1" />

or NCBI. Phylogenetic trees, based on deduced full-length amino acid sequences, were produced using the neighbor-joining (NJ) method and confirmed by the maximum likelihood method (data not shown). Reliabilities of the interfered trees were estimated using the bootstrap method with 1,000 replications.

## mRNA Level Analysis

The relative gene mRNA level was determined by quantitative real-time RT-PCR. Samples were homogenized using Precellys® 24 (Bertin Technologies, Montigny-le-Bretonneux, France) in tubes containing TRIzol reagent (Invitrogen, Carlsbad, CA,

<table><caption>**TABLE 1** | List of the glycerol metabolism-related genes (glycerol kinase – gk; aquaporin 9b – aqp9b; glycerol-3-phosphate phosphatase – pgp; and glycerol-3-phosphate dehydrogenase 1 – gpd1) and their associated primers for quantitative RT-PCR in rainbow trout.</caption><thead><tr><th>Genes</th><th>5'/3' Forward primer</th><th>5'/3' Reverse primer</th><th>Accession numbers</th></tr></thead><tbody><tr><td><strong>Reference gene</strong></td><td></td><td></td><td></td></tr><tr><td>Luciferase control</td><td>CATTCTTCGCCAAAGCACTCTG</td><td>AGCCCATATCCTTGTCGTATCCC</td><td></td></tr><tr><td>EF1α</td><td>TCCTCTTGGTCGTTTCGCTG</td><td>ACCCGAGGGACATCCTGTG</td><td>AF498320</td></tr><tr><td><strong>Glycerol metabolism</strong></td><td></td><td></td><td></td></tr><tr><td>aqp9ba</td><td>TGCTGGATCTTGTGCTGTCT</td><td>CAATGCACAGTGGCTCCAT</td><td>XM_021586222.1</td></tr><tr><td>pgp</td><td>CTGAGGTGAAGGCTGTGGTT</td><td>GGCCACACAGTCGTACATGA</td><td>XM_021576135.1</td></tr><tr><td>gk2a</td><td>TGGGATAAAGCAGAGTTTCCC</td><td>CCGTCTCTTTGTCCCAAACCA</td><td>XM_021598077.1</td></tr><tr><td>gk5</td><td>AAGCCCAATGACTGCCTCTA</td><td>TGCATCCTCCTTCGTCCTAT</td><td>XM_021590025</td></tr><tr><td>gpd1aa</td><td>CAACTCCCAAGAAAGTCTGCA</td><td>GATGTCTGTCAGCTTGCGAC</td><td>XM_021566626.1</td></tr><tr><td>gpd1b</td><td>CAAGAACAAAGAGTGGGGCC</td><td>ATCATCTCCATCAGCCCCAG</td><td>XM_021570149.1</td></tr><tr><td>gpd1c</td><td>TTGGATCTGGAAACTGGGGC</td><td>AGTATGCTGGCCCCTTTCAC</td><td>XM_021571314.1 and XM_021608737.1</td></tr></tbody></table>

gpd1c corresponds to the global amplification of the arbitrarily named gpd1ca and gpd1cb genes.

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Continued</figcaption></figure>

<footnote n="1">¹http://www.ensembl.org</footnote>
B1

<figure><img src="image_4.png" /><figcaption>FIGURE 1 | Continued</figcaption></figure>

B2

<figure><img src="image_5.png" /></figure>

C

<figure><img src="image_6.png" /><figcaption>FIGURE 1 | Continued</figcaption></figure>
<figure><img src="image_7.png" /></figure>
United States). For ontogenesis experiment only, luciferase control RNA (Promega, Charbonnières, France), 10 pg per 1.9 mg of embryo/alevin or oocyte, was added to each sample to allow for data normalization during early development as previously described (Desvignes et al., 2011; Marandel et al., 2012). Total RNA was then extracted according to the manufacturer's instructions. The total RNA concentrations were determined using the spectrophotometer NanoDrop 2000, whereas the total RNA qualities were determined after migration on a 1% agarose gel electrophoresis. Total RNA (1 µg) was used for cDNA synthesis. For the postprandial and nutritional experiment, luciferase control RNA (Promega) was added during the reverse-transcription step (1 pg/µg of RNA). The SuperScript III RNase H-Reverse transcriptase kit (Invitrogen, Carlsbad, California, United States) was used with random primers (Promega) to synthesize cDNA.

The primer sequences used in real-time RT-PCR assays for glycerol metabolism gene analysis were designed and are listed in **Table 1**. We tried to amplify all the genes (see ahead, *In silico* Analysis, in "Results" section), but we only succeeded for genes described in **Table 1**. They were validated on a pool of cDNA, and amplified products were sequenced systematically. Sequencing the qPCR products amplified by *gpd1c*-related primers, we highlighted some polymorphism showing that the two genes *gpd1ca* (XM_021571314.1) and *gpd1cb* (XM_021608737.1) were amplified by the set of primers. These two nucleotidic sequences shared 94.8% identity, thus highly conserved. Moreover, our phylogenetic tree analysis (**Figure 1D**) strongly suggested that *gpd1ca* and *gpd1cb* were orthologous genes emerging from the recent salmonid-specific whole-genome duplication and that they potentially had not diverged. As shown in previous studies, trout ohnolog genes were mostly regulated at the transcription level in the same way (Marandel et al., 2015). We thus decided to use the set of primers that amplify both *gpd1ca* and *gpd1cb* genes.

For real-time RT-PCR assays, the Roche LightCycler 480 system was used (Roche Diagnostics, Neuilly-sur-Seine, France). The assays were performed using a reaction mix of 6 µl per sample, each of which contained 2 µl of diluted cDNA template (1:76 for the postprandial and nutritional experiment; 1:25 for the ontogenesis experiment), 0.24 µl of each primer (10 µmol/L), and 3 µl of LightCycler 480 SYBR® Green I Master mix (Roche Diagnostics). The PCR protocol was initiated at 95°C for 10 min for initial denaturation of the cDNA and hot-start enzyme activation and continued with 45 cycles of a three-step amplification program (15 s at 95°C followed by 10 s 60°C, for primer hybridization, and 15 s at 72°C to extend DNA). Melting curves were systematically monitored (temperature gradient at 0.11°C/s from 65 to 97°C; five acquisitions per degree Celsius) at the end of the last amplification cycle to confirm the specificity of the amplification reaction. Each PCR run included replicate samples (duplicate of reverse transcription and duplicate of PCR amplification) and negative controls (reverse transcriptase-free samples and RNA-free samples).

## Relative Quantification of mRNA Levels

For the analysis of mRNA levels, relative quantification of target gene expression was performed using the Roche Applied Science E-Method. For ontogenesis and postprandial experiments, the

target gene was normalized with exogenous luciferase transcript abundance (Promega). For nutritional experiments, the target gene was normalized with elongation factor 1α (EF1α) measured RNA. As the relative expression of luciferase control and EF1α did not significantly change over the sampling process (data not shown), in all cases, PCR efficiency was measured from the slope of a standard curve using serial dilutions of cDNA; PCR efficiency values ranged between 1.8 and 2.0.

## Statistical Analysis

Molecular data about ontogenesis and postprandial experiments were analyzed by one-way analysis of variance (ANOVA). When significant differences were found among the groups (*P* < 0.05), Tukey's procedure was used to rank the groups. Molecular and plasma parameters of the dietary glycerol experiment were analyzed by two-way ANOVA. In all these analyses, data sets were tested for normality (Shapiro-Wilk) and equal variance (Levene's test). Values are reported as means ± SD. Data analyses were carried out using R software (version 3.4.4).

# RESULTS

## *In silico* Analysis of Rainbow Trout *aqp9b*, *gk*, *pgp*, and *gpd1* Genes

By analyzing the assembly of the rainbow trout genome on NCBI, we identified for the first time several genes sharing high sequence homology with zebrafish glycerol pathway-related genes: two sequences were related to *aqp9b*, three to *gk*, one to *pgp*, and five to *gpd1*. Phylogenetic analyses, using full-length protein sequences, were performed to confirm the identity of trout sequences (**Figure 1**).

### Aquaporin 9b

Two sequences were identified in the trout genome sharing 70.9% (XP_021441897.1) and 71.6% (XP_021463056.1) identity with the zebrafish Aqp9b (ENSDARP00000069995.5) in terms of amino acid sequences. These two sequences were highly conserved with 97% identity in common. Our phylogenetic tree showed that the two trout sequences grouped with the vertebrate Aqp9b orthologous sequences and were close to the two Atlantic salmon (*Salmo salar*) *aqp9b* genes (**Figure 1A**). Together these results confirmed that the two sequences identified in the rainbow trout genome were *aqp9*-related ones. We arbitrarily annotated the two rainbow trout aquaporin 9 genes *aqp9ba* (NCBI ID: 110506537) and *aqp9bb* (NCBI ID: 110526415). The two sequences identified in salmonids (vs. only one in other teleost fish) suggested that this gene duplicated before or around the salmonid radiation following the salmonid whole-genome duplication (SaGD, Berthelot et al., 2014).

### Glycerol Kinase

Two zebrafish glycerol kinase-related genes were identified in the Ensembl database as *gk2* (Ensembl accession no. ENSDARG00000053456) and *gk5* (Ensembl accession no. ENSDARG00000062864). By blasting them against the NCBI rainbow trout assembly analysis, we found two sequences related to *gk2* and one related to *gk5* (NCBI ID: 110509068). Gk2-related sequences were highly conserved (86% identity) and shared 86
and 92% identity with zebrafish Gk2. Our phylogenetic tree showed that the two trout sequences grouped with vertebrate Gk2 orthologous and were close to the two Atlantic salmon (S. salar) gk2 genes (Figure 1B). As for Aqp9, this configuration was in favor of a salmonid-specific duplication occurring after SaGD. We arbitrarily called them gk2a (NCBI ID: 110520656) and gk2b (NCBI ID: 110495542).

For Gk5, the trout sequence shared 76.7% homology with the zebrafish Gk5. Again, its place in the phylogenetic tree (Figure 1B) confirmed its Gk5 identity in trout.

## G3P Phosphatase

Only one sequence was identified in the rainbow trout genome using NCBI BLAST (XP_021431810.1), which shared 79.7% identity with zebrafish Pgp (ENSDARP00000039692.6). This sequence grouped with Pgp vertebrate orthologs in our phylogenetic analysis (Figure 1C) confirming its identity as a Pgp protein in trout.

## G3P Dehydrogenase 1

Three zebrafish G3P dehydrogenase 1 genes were identified in the Ensembl database: *gpd1a* (Ensembl accession no. ENSDARG00000043701), *gpd1b* (Ensembl accession no. ENSDARG00000043180), and *gpd1c* (Ensembl accession no. ENSDARG00000036942). Blasting them against the NCBI rainbow trout assembly, we identified two sequences related to *gpd1a*, both sharing 83% identity with the zebrafish sequence; one related to *gpd1b* (NCBI ID: 110494817), sharing 85% identity with the zebrafish sequence; and two related to *gpd1c*, sharing 84.8 and 84% identity with the zebrafish one. Our phylogenetic analysis (Figure 1D) showed that the trout sequences grouped

with the corresponding orthologs in other vertebrates. Together, our results confirmed the identity of the identified sequences in trout. Finally, the two trout *gpd1a*-related sequences were arbitrarily named *gpd1aa* (NCBI ID: 110492370) and *gpd1ab* (NCBI ID: 110532036), and the two *gpd1c*-related sequences were arbitrarily named *gpd1ca* (NCBI ID: 110495844) and *gpd1cb* (NCBI ID: 110527462).

Based on these molecular findings and because we did not succeed to amplify some of the mRNAs (for *aqp9bb*, *gk2b*, and *gpd1ab* duplicated genes), transcript accession numbers corresponding to genes described above are listed in Table 1.

## Ontogenesis of Glycerol Gene Expression in Rainbow Trout

Real-time PCR was performed to determine the stage-specific expression of rainbow trout glycerol-related genes during embryogenesis and in alevins. Stages were defined based on the Vernier development table (Vernier, 1969) and chosen to target developmental periods of interest for metabolism. In Table 2A, this analysis showed that mRNA levels of *pgp*, *gk5*, *gpd1c*, and *gpd1aa* glycerol-related genes increase significantly at stage 22 (primitive liver) up to stage 30 (hatched embryo). The significant increase of *aqp9b* and *gpd1b* genes appears later at stage 30. The *gk2a* gene has a very different profile of expression during ontogenesis; indeed, the mRNA level of the *gk2a* gene decreases significantly at stage 10 (embryonic genome activation) compared to the oocyte stage. During the nutritional transition (Table 2B), we observed an increase of mRNA levels for *gk2a*, *gk5*, *gpd1b*, and *gpd1c* genes at stage 35 (stage of mixed endogenous-exogenous feeding). Compared to stage 35, there is a significant

<table><caption>TABLE 2A | Ontogenesis of glycerol metabolism-related genes in the whole body of rainbow trout before the first feeding stage.</caption><thead><tr><th rowspan="2">Target genes</th><th colspan="8">Developmental Stages</th></tr><tr><th>Oocyte</th><th>5</th><th>6</th><th>7</th><th>8</th><th>10</th><th>12</th><th>15</th></tr></thead><tbody><tr><td>aqp9ba</td><td>0.11 ± 0.00<sup>a</sup></td><td>0.01 ± 0.01<sup>a</sup></td><td>0.01 ± 0.00<sup>a</sup></td><td>0.00 ± 0.00<sup>a</sup></td><td>0.01 ± 0.01<sup>a</sup></td><td>0.00 ± 0.00<sup>a</sup></td><td>0.00 ± 0.00<sup>a</sup></td><td>0.00 ± 0.00<sup>a</sup></td></tr><tr><td>pgp</td><td>0.01 ± 0.02<sup>a</sup></td><td>0.01 ± 0.01<sup>a</sup></td><td>0.01 ± 0.01<sup>a</sup></td><td>0.00 ± 0.00<sup>a</sup></td><td>0.05 ± 0.04<sup>a</sup></td><td>0.06 ± 0.09<sup>a,b</sup></td><td>0.29 ± 0.29<sup>a,b</sup></td><td>0.39 ± 0.11<sup>b</sup></td></tr><tr><td>gk2a</td><td>3.33 ± 1.62<sup>a</sup></td><td>2.28 ± 1.48<sup>a,b</sup></td><td>1.83 ± 0.86<sup>a,b</sup></td><td>2.73 ± 1.73<sup>a,b</sup></td><td>2.68 ± 1.56<sup>a,b</sup></td><td>1.38 ± 0.50<sup>a,b</sup></td><td>0.28 ± 0.19<sup>b</sup></td><td>0.18 ± 0.02<sup>b</sup></td></tr><tr><td>gk5</td><td>0.04 ± 0.05<sup>a</sup></td><td>0.05 ± 0.05<sup>a,b</sup></td><td>0.04 ± 0.07<sup>a,b</sup></td><td>0.03 ± 0.04<sup>a</sup></td><td>0.05 ± 0.08<sup>a,b</sup></td><td>0.02 ± 0.03<sup>a</sup></td><td>0.04 ± 0.03<sup>a</sup></td><td>0.05 ± 0.02<sup>a,b</sup></td></tr><tr><td>gpd1aa</td><td>0.04 ± 0.02<sup>a</sup></td><td>0.02 ± 0.01<sup>a</sup></td><td>0.01 ± 0.00<sup>a</sup></td><td>0.02 ± 0.01<sup>a</sup></td><td>0.03 ± 0.02<sup>a</sup></td><td>0.01 ± 0.00<sup>a</sup></td><td>0.03 ± 0.02<sup>a</sup></td><td>0.04 ± 0.01<sup>a</sup></td></tr><tr><td>gpd1b</td><td>0.03 ± 0.02<sup>a</sup></td><td>0.03 ± 0.02<sup>a</sup></td><td>0.02 ± 0.01<sup>a</sup></td><td>0.03 ± 0.03<sup>a</sup></td><td>0.04 ± 0.03<sup>a</sup></td><td>0.03 ± 0.02<sup>a</sup></td><td>0.01 ± 0.00<sup>a</sup></td><td>0.01 ± 0.01<sup>a</sup></td></tr><tr><td>gpd1c</td><td>0.05 ± 0.04<sup>a</sup></td><td>0.05 ± 0.02<sup>a</sup></td><td>0.04 ± 0.04<sup>a</sup></td><td>0.06 ± 0.04<sup>a</sup></td><td>0.05 ± 0.04<sup>a</sup></td><td>0.05 ± 0.03<sup>a</sup></td><td>0.03 ± 0.01<sup>a</sup></td><td>0.07 ± 0.02<sup>a</sup></td></tr><tr><td></td><td><strong>22</strong></td><td><strong>23</strong></td><td><strong>31</strong></td><td><strong>p-value</strong></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.01 ± 0.01<sup>a</sup></td><td>0.02 ± 0.01<sup>a</sup></td><td>0.22 ± 0.03<sup>b</sup></td><td>6.30E-16</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>1.63 ± 0.17<sup>c</sup></td><td>1.75 ± 0.05<sup>c</sup></td><td>1.54 ± 0.10<sup>c</sup></td><td>&lt; 2E-16</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.68 ± 0.12<sup>a,b</sup></td><td>0.57 ± 0.08<sup>a,b</sup></td><td>0.54 ± 0.19<sup>a,b</sup></td><td>0.005</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.43 ± 0.17<sup>b,c</sup></td><td>0.47 ± 0.17<sup>c</sup></td><td>0.63 ± 0.34<sup>c</sup></td><td>1.45E-5</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>2.37 ± 0.34<sup>c</sup></td><td>3.20 ± 0.17<sup>d</sup></td><td>1.25 ± 0.19<sup>b</sup></td><td>&lt; 2E-16</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.07 ± 0.02<sup>a</sup></td><td>0.09 ± 0.03<sup>a</sup></td><td>1.09 ± 0.15<sup>b</sup></td><td>&lt; 2E-16</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>0.80 ± 0.12<sup>b</sup></td><td>0.83 ± 0.04<sup>b</sup></td><td>0.77 ± 0.20<sup>b</sup></td><td>8.44E-14</td><td></td><td></td><td></td><td></td></tr></tbody></table>

The stages are those described by Vernier (1969). Description of the stages: oocyte to stage 8, before embryonic genome activation (EGA); stage 10, EGA; stages 12 and 15, epiboly period; stage 22, primitive liver; and stage 23, primitive hepatic portal vein. *gpd1c* corresponds to the global amplification of the arbitrarily named *gpd1ca* and *gpd1cb* genes. Different letters are shown when the mRNA abundances are significantly different. Values are reported as means ± SD.
<table><caption>**TABLE 2B** | Ontogenesis of glycerol metabolism-related genes in the whole body of rainbow trout during the nutritional transition (first feeding stage).</caption><thead><tr><th rowspan="2">Target genes</th><th colspan="3">Developmental stages (first feeding stage)</th><th rowspan="2">p-value</th></tr><tr><th>31</th><th>35</th><th>37</th></tr></thead><tbody><tr><td>aqp9ba</td><td>0.22 ± 0.06<sup>a</sup></td><td>2.04 ± 1.36<sup>b</sup></td><td>1.05 ± 0.36<sup>c</sup></td><td>1.02E-7</td></tr><tr><td>pgp</td><td>1.56 ± 0.34<sup>a</sup></td><td>1.80 ± 1.38<sup>a</sup></td><td>0.79 ± 0.21<sup>b</sup></td><td>1.48E-3</td></tr><tr><td>gk2a</td><td>0.55 ± 0.22<sup>a</sup></td><td>1.48 ± 1.05<sup>b</sup></td><td>1.17 ± 0.51<sup>b</sup></td><td>6.83E-4</td></tr><tr><td>gk5</td><td>0.64 ± 0.42<sup>a</sup></td><td>1.64 ± 1.51<sup>b</sup></td><td>0.75 ± 0.36<sup>a</sup></td><td>3.99E-3</td></tr><tr><td>gpd1aa</td><td>1.26 ± 0.34<sup>a,b</sup></td><td>1.72 ± 1.14<sup>b</sup></td><td>0.70 ± 0.14<sup>a</sup></td><td>2.67E-4</td></tr><tr><td>gpd1b</td><td>1.11 ± 0.34<sup>a</sup></td><td>2.04 ± 1.31<sup>b</sup></td><td>0.78 ± 0.20<sup>a</sup></td><td>4.64E-5</td></tr><tr><td>gpd1c</td><td>0.79 ± 0.29<sup>a</sup></td><td>1.86 ± 1.64<sup>b</sup></td><td>0.96 ± 0.23<sup>a</sup></td><td>0.0033</td></tr></tbody></table>

The stages are those described by Vernier (1969). Description of the stages: stages 31–35–37, hatched embryo/endogenous feeding period. gpd1c corresponds to the global amplification of the arbitrarily named gpd1ca and gpd1cb genes. Different letters are shown when the mRNA abundances are significantly different. Values are reported as means ± SD.

decrease of expression of the glycerol genes (except the *gk2a*) when rainbow trout are fed with exogenous diet (58% protein and 17.5% lipids) (stage 37).

## Postprandial Regulation of Glycerol Metabolism-Related Gene Expression in Rainbow Trout Liver

Table 3 shows the postprandial kinetic of the glycerol metabolism-related gene expression after feeding (commercial diet: 42% of proteins and 22% of lipids). Postprandial times had a significant effect on the expressions of *pgp*, *gk2a*, *gpd1aa*, and *gpd1b* genes but no effect on *aqp9ba*, *gk5*, and *gpd1c* genes (Table 3). Indeed, hepatic expression of glycerol metabolism-related genes (*pgp*, *gk2a*, and *gpd1aa*) revealed significant postprandial increases in expression after refeeding (2 h after the last meal for *gk2a* and > 8 h after the last meal for the other genes), while *gpd1b* decreased 4 h after the last meal.

## Regulation of Glycerol Metabolism-Related Gene Expression by Dietary Glycerol in Rainbow Trout Liver

The last part of the experiments was dedicated to studying the regulation of the glycerol metabolism-related gene expression by dietary glycerol. We fed rainbow trout for 8 weeks with 0, 2.5, and 5% of glycerol. We analyzed the gene expression 6 and 24 h after the last meal. As reflected by the plasma analysis of glycerol, glucose, and triglycerides, there was an expected increase of plasma glycerol in fish fed 5% of glycerol but only 6 h after the last meal (Table 4); higher levels of nutrient absorption 6 h after the last meal are confirmed by the increase of glycemia 6 h after the last meal vs. 24 h after the last meal. Regarding the effects of dietary glycerol on the mRNA levels for glycerol metabolism-related genes, there were no significant effects on their levels (Table 5). Only an effect linked to the postprandial timing was observed for the *pgp* gene which was more highly expressed at 6 h compared to 24 h after the last meal.

## DISCUSSION

Glycerol metabolism is at the interface of lipid and glucose metabolism (Mugabo et al., 2016; Xue et al., 2017; Possick et al., 2017). Up to now, data about genes of glycerol metabolism in rainbow trout (and in other fish species) are scarce; in particular, no studies have been published concerning its molecular regulation. The recent sequencing of the rainbow trout genome should provide new tools to investigate glycerol metabolism in trout, taking into account that this genome has recently suffered a whole duplication event (Berthelot et al., 2014). The present study characterizes the rainbow trout genes coding aquaporin 9 (*aqp9b* and *aqp9b* genes), glycerol kinase (*gk2a/gk2b* and *gk5* genes), G3P phosphatase (*pgp* gene), and cytosolic G3P dehydrogenase (*gpd1* genes, i.e., *gpd1aa*, *gpd1ab*, *gpd1b*, *gpd1ca*, and *gpd1cb* genes) using the genomic sequence of the rainbow trout. Therefore, the present investigation aimed to determine, for the first time in a teleost species, the regulation of *gk2a*, *pgp*, *gk5*, *gpd1aa*, *gpd1b*, and *gpd1c* gene expressions during developmental stages and with the nutritional environment.

## mRNA Levels of Glycerol Metabolism-Related Genes Increase During Ontogenesis From the Setup of the Primitive Liver in Rainbow Trout

Following the characterization of the different genes of the glycerol pathway, the objective of the present paper was to describe the regulation of expressions of these genes during ontogenesis. We observed higher expression of these genes (*pgp*, *gk2a*, *gk5*, *gpd1aa*, *gpd1b*, and *gpd1c*) mainly from stages 22 and 23, corresponding to the setup of the primitive liver and the primitive hepatic portal vein, respectively (Vernier, 1969). These data are not surprising because glycerol metabolism-related genes are known to be mainly expressed in the liver of mammals (Possick et al., 2017). For unclear reasons (no need of glycerol transport at this stage?), the *aqp9b* and *gpd1b* gene expressions appear later at stage 30 just before the first feeding stage. On the other hand, the *gk2a* gene has a particular profile of expression during ontogenesis: indeed, higher levels of *gk2a* mRNAs were observed at the very first stages of the development (i.e., from oocyte stage up to stage 10—blastula), suggesting that this mRNA was maternally inherited. The *gk2a* gene may play an important role during the early cleavage period in rainbow trout. Indeed, an important role of *gk2a* could be to sustain embryo with G3P for its early synthesis of complex lipids (Terner et al., 1968).

## Hepatic Glycerol Metabolism-Related Gene Expressions Are Upregulated After Feeding in Trout Alevins and Juveniles

The first feeding stage is a stage of high molecular plasticity in fish metabolism (Mennigen et al., 2012; Marandel et al., 2016). Indeed, the increase of mRNA levels for all the glycerol metabolism-related genes (*aqp9ba*, *gk2a*, *gk5*, *pgp*, *gpd1aa*, *gpd1b*, and *gpd1c*) was observed during the endogenous-exogenous feeding transition (stage 35). Moreover, by contrast, when the yolk reserves were totally consumed (stage 37),
<table><caption>**TABLE 3** | Postprandial regulation of glycerol metabolism gene expression in rainbow trout liver. 0 h: fasted fish 48 h. 2, 4, 8, 12, 16, and 24 h after the last meal.</caption><thead><tr><th>Target genes</th><th>0 h</th><th>2 h</th><th>4 h</th><th>8 h</th><th>12 h</th><th>16 h</th><th>24 h</th><th>p-value</th></tr></thead><tbody><tr><td>aqp9ba</td><td>1.15 ± 0.68</td><td>0.98 ± 0.43</td><td>0.68 ± 0.35</td><td>1.04 ± 0.16</td><td>0.82 ± 0.37</td><td>0.83 ± 0.45</td><td>0.97 ± 0.32</td><td>0.543</td></tr><tr><td>pgp</td><td>0.17 ± 0.05<sup>a</sup></td><td>0.37 ± 0.06<sup>a,b</sup></td><td>0.80 ± 0.26<sup>a,c</sup></td><td>1.16 ± 0.19<sup>b,c</sup></td><td>1.24 ± 0.62<sup>b,c</sup></td><td>1.54 ± 0.77<sup>c</sup></td><td>0.93 ± 0.57<sup>a,c</sup></td><td>5.83E-05</td></tr><tr><td>gk2a</td><td>0.45 ± 0.18<sup>a</sup></td><td>1.12 ± 0.36<sup>b,c</sup></td><td>0.72 ± 0.18<sup>a,c</sup></td><td>1.34 ± 0.48<sup>c</sup></td><td>1.00 ± 0.26<sup>a,c</sup></td><td>1.10 ± 0.65<sup>a,c</sup></td><td>0.72 ± 0.24<sup>a,b</sup></td><td>0.003</td></tr><tr><td>gk5</td><td>1.29 ± 0.13</td><td>1.16 ± 2.71</td><td>0.05 ± 0.04</td><td>1.60 ± 1.25</td><td>0.94 ± 1.28</td><td>0.49 ± 1.12</td><td>1.15 ± 1.18</td><td>0.317</td></tr><tr><td>gpd1aa</td><td>0.39 ± 0.21<sup>a</sup></td><td>0.47 ± 0.17<sup>a,b</sup></td><td>0.46 ± 0.20<sup>a</sup></td><td>0.89 ± 0.19<sup>b</sup></td><td>0.75 ± 0.21<sup>a,b</sup></td><td>0.82 ± 0.28<sup>a,b</sup></td><td>0.47 ± 0.27<sup>a,b</sup></td><td>6.00E-04</td></tr><tr><td>gpd1b</td><td>0.76 ± 0.24<sup>a</sup></td><td>0.63 ± 0.18<sup>a,b</sup></td><td>0.42 ± 0.15<sup>b</sup></td><td>0.69 ± 0.17<sup>a,b</sup></td><td>0.48 ± 0.11<sup>b</sup></td><td>0.69 ± 0.25<sup>a,b</sup></td><td>0.64 ± 0.15<sup>a,b</sup></td><td>0.031</td></tr><tr><td>gpd1c</td><td>0.85 ± 0.22</td><td>0.84 ± 0.26</td><td>0.76 ± 0.31</td><td>1.13 ± 0.30</td><td>0.96 ± 0.33</td><td>1.13 ± 0.30</td><td>1.22 ± 0.41</td><td>0.101</td></tr></tbody></table>

gpd1c corresponds to the global amplification of the arbitrarily named gpd1ca and gpd1cb genes. Different letters are shown when the mRNA abundances are significantly different. Values are reported as means ± SD.

<table><caption>**TABLE 4** | Plasma metabolites in rainbow trout fed with different levels of glycerol (0, 2.5, and 5%) for 60 days.</caption><thead><tr><th rowspan="2">Metabolites</th><th colspan="3">6 h</th><th colspan="3">24 h</th><th colspan="4">p-value</th></tr><tr><th>0%</th><th>2.5%</th><th>5%</th><th>0%</th><th>2.5%</th><th>5%</th><th>Diet</th><th>Time</th><th>Diet x Time</th></tr></thead><tbody><tr><td>Glucose (g/L)</td><td>0.78 ± 0.10</td><td>0.83 ± 0.13</td><td>0.74 ± 0.18</td><td>0.73 ± 0.09</td><td>0.68 ± 0.08</td><td>0.67 ± 0.08</td><td>0.419</td><td>0.011</td><td>0.436</td></tr><tr><td>Triglyceride (g/L)</td><td>3.25 ± 0.86</td><td>3.44 ± 1.73</td><td>5.09 ± 1.68</td><td>2.97 ± 1.37</td><td>3.13 ± 1.72</td><td>3.47 ± 1.32</td><td>0.065</td><td>0.091</td><td>0.349</td></tr><tr><td>Glycerol (mM)</td><td>0.05 ± 0.00</td><td>0.06 ± 0.01</td><td>0.08 ± 0.05</td><td>0.04 ± 0.01</td><td>0.04 ± 0.01</td><td>0.04 ± 0.00</td><td>0.131</td><td>9.27E-08</td><td>1.46E-03</td></tr></tbody></table>

Values are reported as means ± SD.

<table><caption>**TABLE 5** | Regulation of glycerol metabolism gene expression in rainbow trout liver by dietary glycerol for 60 days at 6 and 24 h after the last meal.</caption><thead><tr><th rowspan="2">Target genes</th><th colspan="3">6 h</th><th colspan="3">24 h</th><th colspan="4">p-value</th></tr><tr><th>0%</th><th>2.5%</th><th>5%</th><th>0%</th><th>2.5%</th><th>5%</th><th>Diet</th><th>Time</th><th>Diet x Time</th></tr></thead><tbody><tr><td>aqp9ba</td><td>0.91 ± 0.37</td><td>0.94 ± 0.49</td><td>1.17 ± 0.72</td><td>0.94 ± 0.59</td><td>0.99 ± 0.63</td><td>1.27 ± 0.75</td><td>0.377</td><td>0.751</td><td>0.987</td></tr><tr><td>pgp</td><td>1.53 ± 0.71</td><td>1.68 ± 0.44</td><td>1.61 ± 0.43</td><td>0.47 ± 0.21</td><td>0.38 ± 0.15</td><td>0.44 ± 0.15</td><td>0.983</td><td>1.17E-11</td><td>0.657</td></tr><tr><td>gk2a</td><td>1.06 ± 0.49</td><td>0.95 ± 0.47</td><td>0.90 ± 0.39</td><td>1.21 ± 0.23</td><td>0.93 ± 0.47</td><td>1.04 ± 0.65</td><td>0.499</td><td>0.530</td><td>0.848</td></tr><tr><td>gk5</td><td>1.56 ± 1.33</td><td>0.87 ± 0.26</td><td>0.86 ± 0.68</td><td>1.04 ± 0.77</td><td>1.31 ± 0.55</td><td>0.94 ± 0.53</td><td>0.366</td><td>0.999</td><td>0.220</td></tr><tr><td>gpd1aa</td><td>0.83 ± 0.48</td><td>1.09 ± 0.37</td><td>1.2 ± 0.46</td><td>0.91 ± 0.21</td><td>0.77 ± 0.30</td><td>0.83 ± 0.29</td><td>0.587</td><td>0.076</td><td>0.206</td></tr><tr><td>gpd1b</td><td>0.89 ± 0.40</td><td>1.10 ± 0.45</td><td>1.03 ± 0.41</td><td>0.94 ± 0.30</td><td>0.84 ± 0.26</td><td>1.20 ± 0.23</td><td>0.326</td><td>0.892</td><td>0.255</td></tr><tr><td>gpd1c</td><td>0.96 ± 0.19</td><td>1.14 ± 0.33</td><td>1.23 ± 0.51</td><td>0.93 ± 0.13</td><td>0.94 ± 0.28</td><td>0.97 ± 0.14</td><td>0.372</td><td>0.082</td><td>0.574</td></tr></tbody></table>

gpd1c corresponds to the global amplification of the arbitrarily named gpd1ca and gpd1cb genes. Values are reported as means ± SD.

i.e., when the fish were fed only with inert diets, there is a decrease in the levels of many of the studied mRNAs. Just before the end of the resorption of the yolk reserves (stage 35), the majority of the remaining nutrients are lipids and especially triacylglycerol (Kamler, 2008). Yolk triglycerides would be the main endogenous stores at stage 35, with the preferential utilization by the alevins of fatty acids and the possible use of 3-carbon backbone of glycerol for glucose synthesis through gluconeogenesis (Vernier and Sire, 1976) as described in mammals and fish (Rito et al., 2019). Despite several methodological caveats, the use of glycerol with carbon isotopes, either delivered in labeled pellets (Da Costa et al., 2017) or delivered in incubation medium for hepatocytes (Lech, 1970) or IP injected (Savina and Wojtczak, 1977; Rito et al., 2019), assisted greatly in confirming that fish species are able to catabolize it as well as retain and incorporate it into carbohydrate, namely, glucose via gluconeogenesis. This use of glycerol from endogenous feeding could explain the increase of the glycerol metabolism-related gene expressions at this stage.

Later (i.e., stage 37 here), when the yolk is fully consumed and its energetic content depleted, their nutrition is from an inert commercial feed consisting of 58% of protein and 17% of lipids (Marandel et al., 2016). Contrary to stage 35, the alevins of this carnivorous species can use, and even prioritize, the catabolism of amino acids as the main energy source, being a possible explanation for the lower levels of glycerol metabolism-related gene expression.

It was also very important to know about the postprandial regulation of glycerol metabolism-related gene expression in the liver as it may reflect glycerol absorption, transport, sensing, and processing mechanisms. Indeed, it is well known that some metabolic genes can be strongly regulated in the liver after the last meal as shown for example with some hepatic genes coding key enzymes of intermediary metabolism (Mennigen et al., 2012). Here, in rainbow trout juveniles fed with a commercial diet (44% proteins and 22% lipids), the significant postprandial change of gene expressions is only observed for four hepatic genes of the glycerol metabolic pathway: pgp, gk2a, gpd1c,
and *gpd1aa*. The postprandial increase of mRNA abundances for *pgp*, *gk2a*, and *gpd1aa* is observed between 2 and 16 h after the last meal. This kinetic of postprandial regulation is similar to the one observed previously for some nutritionally dependent gene expressions (for example, for glycolytic and gluconeogenic genes) (Mennigen et al., 2012). These data strongly suggest that glycerol metabolism could be dependent on the nutritional status of fish.

## Hepatic Glycerol Metabolism-Related Gene Expressions Are Poorly Controlled by Different Levels of Dietary Glycerol

Even though the glycerol metabolism is dependent on the nutritional status (see before), when we fed the fish with different levels of dietary glycerol (0, 2.5, and 5%), we did not observe any significant effects on the hepatic glycerol metabolism-related gene expression (i.e., for the *gk2a*, *gk5*, *pgp*, *gpd1aa*, *gpd1b*, and *gpd1c* genes). Even though the significant increase 6 h after the last meal of the plasma glycerol with dietary glycerol levels proves that glycerol was absorbed at the intestine level and so, in theory, available for utilization in the liver, we must note here that the values of the plasma glycerol are low (max 0.08 mM in trout fed 5% glycerol) if we compared these data with European seabass fed the same diets (up to 6 mM of plasma glycerol; unpublished data). It is hypothesized that the levels of plasma glycerol were not sufficiently high in rainbow trout fed different levels of dietary glycerol to be able to modify the levels of hepatic gene expression from glycerol metabolism. An NMR-based metabolomic analysis using the remaining fish from the nutritional sampling here presented seemed to corroborate these observations. It revealed that the hepatic metabolome in trout was also poorly modified by these levels of dietary glycerol (Palma et al., 2019), which articulates well with the present molecular data. Despite this apparent weak effect on hepatic metabolism, it is noteworthy that the triacylglycerol turnover rate in rainbow trout was relatively high (Magnoni et al., 2008). Instead of being metabolized, glycerol can be perceived as excess carbohydrate and be diverted to support triacylglycerol synthesis. This seems to be the case in trout as in other carnivorous fish species fed with high levels of dietary carbohydrates that either reveal high triacylglycerol hepatic levels (Viegas et al., 2016; Kamalam et al., 2017) or export that excess energy to extrahepatic tissues like the muscle or perivisceral fat (Kamalam et al., 2017; Viegas et al., 2019). Such possible metabolic fate of dietary glycerol could also explain the negative effect of dietary glycerol on nitrogen efficiency detected as previously measured during the dietary trial (Magnoni et al., unpublished data).

## CONCLUSION

Overall, our data characterize, for the first time in fish species, all the genes involved in glycerol metabolism in rainbow trout and their expression during ontogenesis. Many of these genes were regulated by the nutritional status from the first feeding stage up to the juvenile stage. However, the incorporation of increased

levels of glycerol (up to 5%) did not modify the mRNA levels for the glycerol genes.

## DATA AVAILABILITY STATEMENT

The authors acknowledge that the data presented in this study must be deposited and made publicly available in an acceptable repository, prior to publication. Frontiers cannot accept a manuscript that does not adhere to our open data policies.

## ETHICS STATEMENT

The study was conducted under the supervision of accredited experts in laboratory animal science by the Portuguese Veterinary Authority (1005/92, DGV-Portugal, following FELASA category C recommendations), according to the guidelines on the protection of animals used for scientific purposes from the European directive 2010/63/UE.

## AUTHOR CONTRIBUTIONS

SP and IV contributed to the conceptualization and supervision. SP, EP-J, and LM contributed to the formal analysis. IV contributed to the funding acquisition, project administration, and resources. EP-J, EG, SP, LM, LJM, and IV contributed to the investigation. EP-J, EG, and SP contributed to the methodology. SP, LM, LJM, and IV contributed to the validation. EP-J and SP contributed to the writing of the original draft. LJM, LM, and IV contributed to the writing - review and editing. All authors contributed to the article and approved the submitted version.

## FUNDING

This work was supported by Fundação para a Ciência e Tecnologia (FCT; Portugal) through national funds with co-funding from ERDF/FEDER, within the PT2020 Partnership Agreement, and COMPETE 2020: research grant to IV (POCI-01-0145-FEDER-016828 - PTDC/CVT-NUT/2851/2014); individual grant to MP through Centro2020 (ReNATURE; Centro-01-0145-FEDER-000007); and structural funds to Center for Neuroscience and Cell Biology (UID/NEU/04539/2013) and Centre for Functional Ecology (UID/BIA/04004/2019) and Interdisciplinary Centre of Marine and Environmental Research (UID/Multi/04423/2019). The research about gene ontogenesis was internally funded by the European Commission (European project FP7-KBBE-2011-5, ARRAINA project no. 288925, Advanced Research Initiatives for Nutrition 1241 and Aquaculture).

## ACKNOWLEDGMENTS

We are grateful to Renato Ferraz (ICBAS/CIIMAR) and Sara Jorge, Sara Basto, and Paulo Rema (CECAV-UTAD) for technical
assistance in the trials and during sampling during the dietary glycerol trial. We are also grateful to Patrick Maunas and Nicolas Turonnet (INRAE fish farm in Lées-Athas, France) and Frederic

Terrier, Frank Sandres, and Anthony Lanuque (INRAE fish farm in Donzacq, France) for the ontogenesis and postprandial kinetics trials, respectively.

## REFERENCES

Ayoub, M., and Abdullah, A. Z. (2012). Critical review on the current scenario and significance of crude glycerol resulting from biodiesel industry towards more sustainable renewable energy industry. *Renew. Sustain. Energy Rev.* 16, 2671-2686. doi: 10.1016/j.rser.2012.01.054

Berthelot, C., Brunet, F., Chalopin, D., Juanchich, A., Bernard, M., Noë, I. B., et al. (2014). The rainbow trout genome provides novel insights into evolution after whole-genome duplication invertebrates. *Nat. Commun.* 5:3657.

Bolliet, V., Cheewasedtham, C., Houlihan, D., Gélineau, A., and Boujard, T. (2000). Effect of feeding time on digestibility, growth performance and protein metabolism in the rainbow trout *Oncorhynchus mykiss*: interactions with dietary fat levels. *Aquat. Living Resour.* 13, 107-113. doi: 10.1016/S0990-7440(00) 00147-9

Brisson, D., Vohl, M. C., St-Pierre, J., Hudson, T. J., and Gaudet, D. (2001). Glycerol: a neglected variable in metabolic processes? *BioEssays* 23, 534-542. doi: 10.1002/bies.1073

Bromley, P. J., and Adkins, T. C. (1984). The influence of cellulose filler on feeding, growth and utilization of protein and energy in rainbow trout, *Salmo gairdnerii* Richardson. *J. Fish Biol.* 24, 235-244. doi: 10.1111/j.1095-8649.1984.tb 04794.x

Da Costa, D. V., Dias, J., Colen, R., et al. (2017). Partition and metabolic fate of dietary glycerol in muscles and liver of juvenile tilapia. *Arch. Anim. Nutr.* 71, 165-174. doi: 10.1080/1745039x.2017.1281579

Desvignes, T., Fauvel, C., and Bobe, J. (2011). The NME gene family in zebrafish oogenesis and early development. *Naunyn Schmiedebergs Arch Pharmacol.* 384, 439-449. doi: 10.1007/s00210-011-0619-9

Ditlecadet, D. and Driedzic, W. R. (2014). Glycerol synthesis in freeze-resistant rainbow smelt: towards the characterization of a key enzyme glycerol-3-phosphatase. *Fish Physiol. Biochem.* 40, 257-266. doi: 10.1007/s10695-013-9841-3

Dong, C., Chen, L., Feng, J., Xu, J., Mahboob, S., Al-Ghanim, K., et al. (2016). Genome wide identification, phylogeny, and expression of *Aquaporin* genes in common carp (Cyprinus carpio). *PLoS One* 11:e0166160. doi: 10.1371/journal. pone.0166160

Driedzic, W. R. (2015). Rainbow smelt: the unusual case of cyprotection by sustained glycerol production in an aquatic animal. *J. Comp. Physiol. Part B* 185, 487-499. doi: 10.1007/s00360-015-0903-y

Edgar, R. C. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. *BMC Bioinformatics* 5:113. doi: 10.1186/1471-2105-5-113

Food Agric Organ (2018). *The State of the World Fisheries and Aquaculture 2018*. Rome: Food Agric Organ.

Gatlin, D. M., Barrows, F. T., Brown, P., Dabrowski, K., Gaylor, T. G., et al. (2007). Expanding the utilization of sustainable plant products in aquafeeds: a review. *Aquacult. Nutr.* 38, 551-579. doi: 10.1111/j.1365-2109.2007.01704.x

Gonçalves, L. U., Cerozi, B. D. S., Silva, T. S. C., Zanon, R. B., Cyrino, J. E. (2015). Crude glycerin as dietary energy source for Nile tilapia. *Aquaculture* 437, 230-234. doi: 10.1016/j.aquaculture.2014.12.004

Hilton, J. W. (1982). The effect of pre-fasting diet and water temperature on liver glycogen and liver weight in rainbow trout, *Salmo gairdneri* Richardson, during fasting. *J. Fish Biol.* 20, 69-78. doi: 10.1111/j.1095-8649.1982.tb03895.x

Jelen, S., Wacker, S., Aponte-Santamaría, C., Skott, M., Rojek, A., Johanson, U., et al. (2011). Aquaporin-9 protein is the primary route of hepatocyte glycerol uptake for glycerol gluconeogenesis in mice. *J. Biol. Chem.* 286, 44319-44325. doi: 10.1074/jbc.M111.297002

Kamalam, B. J., Medale, F., and Panserat, S. (2017). Utilisation of dietary carbohydrates in farmed fishes: New insights on influencing factors, biological limitations and future strategies. *Aquaculture* 467, 3-27. doi: 10.1016/j.aquaculture.2016.02.007

Kamler, E. (2008). Resource allocation in yolk-feeding fish. *Rev. Fish Biol. Fisher.* 18, 143-200. doi: 10.1007/s11160-007-9070-x

Klinger, D., and Naylor, R. (2012). Searching for solutions in aquaculture: charting a sustainable course. *Annu. Rev. Environ. Resour.* 37, 247-276. doi: 10.1146/ annurev-environ-021111-161531

Lazzarotto, V., Medale, F., Larroquet, L., and Corraze, G. (2018). Long-term dietary replacement of fishmeal and fish oil in diets for rainbow trout (*Oncorhynchus mykiss*): effects on growth, whole body fatty acids and intestinal and hepatic gene expression. *PLoS One* 13:e0190730. doi: 10.1371/journal.pone.0190730

Lebeck, J. (2014). Metabolic impact of the glycerol channels AQP7 and AQP9 in adipose tissue and liver. *J. Mol. Endocrinol.* 52, R165-R178.

Lech, J. J. (1970). Glycerol kinase and glycerol utilization in trout (Salmo gairdneri) liver. *Compar. Biochem. Physiol.* 34, 117-124. doi: 10.1016/0010-406x(70) 90058-7

Magnoni, L., Scarlato, N. A., Patricio, O. F., and Wöhler, O. C. (2013). Gluconeogenic pathway does not display metabolic cold adaptation in liver of Antarctic notothenioid fish. *Polar Biol.* 36, 661-671. doi: 10.1007/s00300-013- 1292-x

Magnoni, L., Vaillancourt, E., and Weber, J.-M. (2008). High resting triacylglycerol turnover of rainbow trout exceeds the energy requirements of endurance swimming. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* 295, R309-R315.

Marandel, L., Labbe, C., Bobe, J., and Le Bail, P.-Y. (2012). Evolutionary history of *c-myc* in teleosts and characterization of the duplicated *c-myc* genes in goldfish embryos. *Mol. Reprod. Dev.* 79, 85-96. doi: 10.1002/mrd.22004

Marandel, L., Seiliez, I., Véron, V., Skiba-Cassy, S., and Panserat, S. (2015). New insights into the nutritional regulation of gluconeogenesis in carnivorous rainbow trout (*Oncorhynchus mykiss*): a gene duplication trail. *Physiol. Genomics* 47, 253-263. doi: 10.1152/physiogenomics.00026.2015

Marandel, L., Veron, V., Surget, A., Plagnes-Juan, E., and Panserat, S. (2016). Glucose metabolism ontogenesis in rainbow trout (*Oncorhynchus mykiss*) in the light of the recently sequenced genome: new tools for intermediary metabolism programming in the light of the recently sequenced genome. *J Exp. Biol.* 219, 734-743. doi: 10.1242/jeb.134304

Mennigen, J., Skiba-Cassy, S., and Panserat, S. (2013). Ontogenetic expression of metabolic genes and microRNAs in rainbow trout alevins during the transition from the endogenous to exogenous feeding period. *J. Exp. Biol.* 216, 1597-1608. doi: 10.1242/jeb.082248

Mennigen, J. A., Panserat, S., Larquier, M., Plagnes-Juan, E., Medale, F., Seiliez, I., et al. (2012). Postprandial regulation of hepatic microRNAs predicted to target the insulin pathway in rainbow trout. *PLoS One* 7:e38604. doi: 10.1371/journal. pone.0038604

Menton, D. J., Slinger, S. J., and Hilton, J. W. (1986). Utilization of free glycerol as a source of dietary energy in rainbow trout (Salmo gairdneri). *Aquaculture* 56, 215-227. doi: 10.1016/0044-8486(86)90337-6

Meurer, F., Franzen, A., Piovesan, P., Rossato, K. A., and Santos, L. D. D. (2012). Apparent energy digestibility of glycerol from biodiesel production for Nile tilapia (*Oreochromis niloticus*, Linnaeus 1758). *Aquacult. Res.* 43, 1734-1737. doi: 10.1111/j.1365-2109.2011.02974.x

Meurer, F., Tovo, N. A., Da Silva, L. C. R., Cagol, L., Theisen, M. T., and Dos Santos, L. D. (2016). Crude glycerol in diets for Nile tilapia sex reversal (*Oreochromis niloticus*, Linnaeus 1758). *Aquacult. Res.* 47, 2682-2685. doi: 10.1111/are.1 2714

Moesch, A., Meurer, F., Zadinelo, I. V., Carneiro, W. F., Da Silva, L. C. R., and Dos Santos, L. D. (2016). Growth, body composition and hepatopancreas morphology of Nile tilapia fingerlings fed crude glycerol as a replacement for maize in diets. *Anim. Feed Sci. Technol.* 219, 122-131. doi: 10.1016/j.anifeedsci. 2016.05.009

Mugabo, Y., Zhao, S., Seifried, A., Gezzar, S., Al-Mass, A., Zhang, D., et al. (2016). Identification of a mammalian glycerol-3-phosphate phosphatase: role in metabolism and signaling in pancreatic beta-cells and hepatocytes. *Proc. Natl. Acad. Sci. U.S.A.* 113, E430-E439.

Naylor, R. L., Hardy, R. W., Bureau, D., Chiu, A., Elliott, M., Farrell, A. P., et al. (2009). Feeding aquaculture in an era of finite resources. *Proc. Natl. Acad. Sci. U.S.A.* 106, 15103-15110. doi: 10.1073/pnas.0905235106
Neu, D. H., Furuya, W. M., Boscolo, W. R., Potrich, F. R., Lui, T. A., and Feiden, A. (2013). Glycerol inclusion in the diet of Nile tilapia (*Oreochromis niloticus*) juveniles. *Aquacult. Nutr.* 19, 211–217. doi: 10.1111/j.1365-2095.2012.00698.x

Palma, M., Tavares, L. C., Rito, J., Henriques, L. F., Silva, J. G., Ozório, R., et al. (2019). Metabolic effects of dietary glycerol supplementation in muscle and liver of european seabass and rainbow trout by 1H NMR metabolomics. *Metabolites* 9:202. doi: 10.3390/metabo9100202

Panserat, S., Hortopan, G. A., Plagnes-Juan, E., Kolditz, C., Lansard, M., Skiba-Cassy, S., et al. (2009). Differential gene expression after total replacement of dietary fish meal and fish oil by plant products in rainbow trout (Oncorhynchus mykiss) liver. *Aquaculture* 294, 123–131. doi: 10.1016/j.aquaculture.2009.05.013

Panserat, S., Marandel, L., Seiliez, I., and Skiba, S. (2019). New insights on intermediary metabolism for a better understanding of nutrition in teleosts. *Ann. Rev. Anim. Biosci.* 7, 195–220. doi: 10.1146/annurev-animal-020518-115250

Patsouris, D., Mandard, S., Voshol, P. J., Escher, P., Tan, N. S., Havekes, L. M., et al. (2016). PPARα governs glycerol metabolism. *J. Clin. Investigat.* 114, 94–103.

Pearson, M. P., Spriet, L. L., and Stevens, D. (1990). Effect of sprint training on swim performance and white muscle metabolism during exercise and recovery in rainbow trout (*Salmo gairdneri*). *J. Exp. Biol.* 149, 45–60.

Possick, E., Madiraiu, M., and Prentki, M. (2017). Glycerol-3-phosphate phosphatase/PGP: role in intermediary metabolism and target for cardiometabolic diseases. *Biochimie* 143, 18–28. doi: 10.1016/j.biochi.2017.08.001

Quispe, C. A. G., Coronado, C. J. R., Carvalho, J. A. Jr. (2013). Glycerol: Production, consumption, prices, characterization and new trends in combustion. *Renew. Sustain. Energy Rev.* 27, 475–493. doi: 10.1016/j.rser.2013.06.017

Rito, J., Viegas, I., Pardal, M. Â, Meton, I., Baanante, I. V., and Jones, J. G. (2019). Utilization of glycerol for endogenous glucose and glycogen synthesis in seabass (Dicentrarchus labrax): A potential mechanism for sparing amino acid catabolism in carnivorous fish. *Aquaculture* 498, 488–495. doi: 10.1016/j.aquaculture.2018.08.066

Savina, M. V., and Wojtczak, A. B. (1977). Enzymes of gluconeogenesis and the synthesis of glycogen from glycerol in various organs of the lamprey (Lampetra fluviatilis). *Compar. Biochem. Physiol. Part B* 57, 185–190. doi: 10.1016/0305-0491(77)90141-9

Sheridan, M. A. (1988). Lipid dynamics in fish: aspects of absorption, transportation, deposition and mobilization. *Compar. Biochem. Physiol. Part B* 90, 679–690. doi: 10.1016/0305-0491(88)90322-7

Sriram, G., Rahib, L., He, J.-S., Campos, A. E., Parr, L. S., Liao, J. C., et al. (2008). Global metabolic effects of glycerol kinase overexpression in rat hepatoma cells. *Mol. Genet. Metab.* 93, 145–159. doi: 10.1016/j.ymgme.2007.09.008

Suarez, R. K., and Mommsen, T. P. (1987). Gluconeogenesis in teleost fishes. *Canad. J. Zool.* 65, 1869–1882. doi: 10.1139/z87-287

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* 30, 2725–2729. doi: 10.1093/molbev/mst197

Tan, H. W., Abdul Aziz, A. R., and Aroua, M. K. (2013). Glycerol production and its applications as a raw material: a review. *Renew. Sustain. Energy Rev.* 27, 118–127. doi: 10.1016/j.rser.2013.06.035

Terner, C., Kumar, L. A., and Choe, T. S. (1968). Studies of metabolism in embryonic development: biosynthesis of lipids in embryonated trout ova. *Compar. Biochem. Physiol.* 24, 941–950. doi: 10.1016/0010-406x(68)90806-2

Vernier, J. M. (1969). Table chronologique du développement embryonnaire de la truite arc-en-ciel. *Ann. Embryol. Morphol.* 2, 495–520.

Vernier, J. M., and Sire, M. F. (1976). Evolution of the glycogen content and of glucose-6-phosphatase activity in the liver of *Salmo gairdneri* during development. *Tissue and Cell* 8, 531–546. doi: 10.1016/0040-8166(76)90011-2

Viegas, I., Jarak, I., Rito, J., Carvalho, R. A., Metón, I., Pardal, M. A., et al. (2016). Effects of dietary carbohydrate on hepatic de novo lipogenesis in European seabass (Dicentrarchus labrax L.). *J. Lipid Res.* 57, 1264–1272. doi: 10.1194/jlr.m067850

Viegas, I., Trenkner, L. H., Rito, J., Palma, M., Tavares, L. C., Jones, J. G., et al. (2019). Impact of dietary starch on extrahepatic tissue lipid metabolism in farmed European (Dicentrarchus labrax) and Asian seabass (Lates calcarifer). *Compar. Biochem. Physiol. Part A* 231, 170–176. doi: 10.1016/j.cbpa.2019.02.025

Westergaard, N., Madsen, P., and Lundgren, K. (1998). Characterization of glycerol uptake and glycerol kinase activity in rat hepatocytes cultured under different hormonal conditions. *Biochim Biophys Acta* 1402, 261–268. doi: 10.1016/s0167-4889(98)00016-0

Wiesmann, D., and Pfeffer, E. (1986). Influence of indigestible carbohydrates on the efficiency of utilization of dietary energy and protein in growing rainbow trout (*Salmo gairdnerii*, R). *Arch. Anim. Nutr.* 36, 1145–1149. doi: 10.1080/17450398609434362

Xue, L. L., Chen, H. H., and Jiang, J. G. (2017). Implications of glycerol metabolism for lipid production. *Prog. Lipid Res.* 68, 17–25.

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Panserat, Plagnes-Juan, Gazzola, Palma, Magnoni, Marandel and Viegas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.