# Modulation of Energy Metabolism and Epigenetic Landscape in Rainbow Trout Fry by a Parental Low Protein/High Carbohydrate Diet

Thérèse Callet, Hongyan Li, Pascale Coste, Stéphane Glise, Cécile Heraud, Patrick Maunas, Yvan Mercier, Nicolas Turonnet, Chloé Zunzunegui, Stéphane Panserat, et al.

► To cite this version:

Thérèse Callet, Hongyan Li, Pascale Coste, Stéphane Glise, Cécile Heraud, et al. Modulation of Energy Metabolism and Epigenetic Landscape in Rainbow Trout Fry by a Parental Low Protein/High Carbohydrate Diet. Biology, 2021, 10 (7), 10.3390/biology10070585. [https://doi.org/10.32720/bio.2021.10070585](https://doi.org/10.32720/bio.2021.10070585)

HAL Id: hal-03272040

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

Submitted on 28 Jun 2021

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
Article

# Modulation of Energy Metabolism and Epigenetic Landscape in Rainbow Trout Fry by a Parental Low Protein/High Carbohydrate Diet

Thérèse Callet<fnref n="1" />,*, Hongyan Li<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Pascale Coste<fnref n="4" />, Stéphane Glise<fnref n="4" />, Cécile Heraud<fnref n="1" />,*, Patrick Maunas<fnref n="1" />, Yvan Mercier<fnref n="1" />, Nicolas Turonnet<fnref n="1" />, Chloé Zunzunegui<fnref n="1" />, Stéphane Panserat<fnref n="1" />, Valérie Bolliet<fnref n="4" /> and Lucie Marandel<fnref n="1" />,*

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

check for
updates

**Citation:** Callet, T.; Li, H.; Coste, P.; Glise, S.; Heraud, C.; Maunas, P.; Mercier, Y.; Turonnet, N.; Zunzunegui, C.; Panserat, S.; et al. Modulation of Energy Metabolism and Epigenetic Landscape in Rainbow Trout Fry by a Parental Low Protein/High Carbohydrate Diet. *Biology* **2021**, *10*, 585. https://doi.org/10.3390/biology10070585

Academic Editor: Jean-Nicolas Vollf

Received: 20 April 2021

Accepted: 20 June 2021

Published: 25 June 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. License MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** While the effects of parental diets on their progeny have been highly described in mammals, such studies are lacking in fish. To explore such a question in a high trophic level teleost fish, two-year old male and female rainbow trout were fed either a control diet (0% carbohydrate and 63.89% protein) or a high-carbohydrate diet (35% carbohydrate and 42.96% protein), for a complete reproductive cycle for females and for a period of 5 months for males. Neither the maternal nor the paternal high-carbohydrate diet alone had induced significant effects on their progeny. Nevertheless, when both parents were fed the high-carbohydrate diet, the energy metabolism and mitochondrial dynamics of their progeny were altered. Moreover, the epigenetic landscape was also highly affected. Even though, offspring growth was only slightly affected at the early stage of life; the effect of parental high-carbohydrate diet should be explored over the long term.

**Abstract:** It is now recognized that parental diets could highly affect offspring metabolism and growth. Studies in fish are, however, lacking. In particular, the effect of a parental diet high in carbohydrate (HC) and low in protein (LP) on progeny has never been examined in higher trophic level teleost fish. Thus, two-year old male and female rainbow trout (*Oncorhynchus mykiss*) were fed either a control diet (0% carbohydrate and 63.89% protein) or a diet containing 35% carbohydrate and 42.96% protein (HC/LP) for a complete reproductive cycle for females and over a 5-month period for males. Cross-fertilizations were then carried out. To evaluate the effect of the parental diet on their offspring, different phenotypic and metabolic traits were recorded for offspring before their first feeding and again three weeks later. When considering the paternal and maternal HC/LP nutrition independently, fry phenotypes and transcriptomes were only slightly affected. The combination of the maternal and paternal HC/LP diets altered the energy metabolism and mitochondrial dynamics of their progeny, demonstrating the existence of a synergistic effect. The global DNA methylation of whole fry was also highly affected by the HC/LP parental diet, indicating that it could be one of the fundamental mechanisms responsible for the effects of nutritional programming.

**Keywords:** nutritional programming; teleost fish; offspring; broodstock; mitochondria

<footnote n="1">¹ INRAE, Université de Pau et des Pays de L'Adour, E2S UPPA, NUMEA, 64310 Saint-Pée-sur-Nivelle, France; therese.callet@inrae.fr (T.C.); lihy@ihb.ac.cn (H.L.); cecile.heraud@inrae.fr (C.H.); patrick.maunas@inrae.fr (P.M.); yvan.mercier@inrae.fr (Y.M.); nicolas.turonnet@inrae.fr (N.T.); czunzunegui97@orange.fr (C.Z.); stephane.panserat@inrae.fr (S.P.)</footnote>

<footnote n="2">² State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China</footnote>

<footnote n="3">³ University of Chinese Academy of Sciences, Beijing 100190, China</footnote>

<footnote n="4">⁴ INRAE, Université de Pau et des Pays de L'Adour, E2S UPPA, ECOBIOP, 64310 Saint-Pée-sur-Nivelle, France; pascale coste-heinrich@inrae.fr (P.C.); stephane.glise@inrae.fr (S.G.); valerie.bolliet@inrae.fr (V.B.)</footnote>

<footnote>* Correspondence: lucie.marandel@inrae.fr</footnote>
# 1. Introduction

Exposure to environmental cues during phases with high plasticity (typically early life) can induce long-lasting changes in the morphology, physiology, and metabolism of an individual. This concept, referred to as programming, has been widely studied in various species and particularly for mammals with regard to the development of chronic diseases (developmental origins of health and disease (DOHaD)) [1]. Of particular interest is the nutritional programming, i.e., the consequences of nutritional stimulus that occur during the prenatal period [2], such as high carbohydrate (HC)/low protein (LP) parental diets [3–5]. Surprisingly, the numerous studies conducted on the effect of parental HC/LP diets in different species report strikingly common features [6,7]. For instance, it has been demonstrated that the offspring of parents fed HC/LP diets exhibit low birth weight [8]; metabolic disorders, such as impaired carbohydrate and lipid metabolisms [9]; and impaired cardiometabolic health [3,10].

Despite the increasing number of studies on nutritional programming, the mechanisms behind the imprinting of programming events and their subsequent effects on generations are not yet fully understood [11]. In mammals, several factors that can mediate the consequences of parental diet on the physiology of their offspring have been identified. The quality and quantity of nutrients provided via the mother’s diet are known to affect the phenotype of their offspring in the long term by directly affecting tissue development, which can lead to developmental alterations and defects in organ function [7]. Modifications of epigenetic marks are also a plausible mechanism to explain the effect of programming due to their susceptibility to environmental factors, their reversibility, and their role in the regulation of gene expression [2,12]. Dietary factors could alter different epigenetic marks, such as DNA methylation, either directly by affecting the availability of methyl donors or indirectly by affecting the availability of cofactors needed for enzymes involved in methylation processes [13]. Hence, HC/LP diets have been shown to induce either hypo- or hyper-methylation of DNA at specific sites [2]. Finally, increasing evidence also points to mitochondria as a key actor in nutritional programming. Disruption of mitochondrial activity due to nutritional events could alter the epigenome because of the role of mitochondria in producing important cofactors required for enzymes involved in epigenetic modifications [14].

In teleost fish, although only few studies have been conducted, the concept of programming has also been demonstrated in different species [15], and a particular focus has been given to aquaculture species [16]. Several studies have focused on different stimuli given during the early feeding phase and their long term consequences [17–21]. The consequences of stimuli applied earlier in development i.e., during gametogenesis, have also been tested. As such, the effects of a decrease in dietary methionine [22], complete plant-based diets [23], and selenium content [24] have been studied in rainbow trout. Nevertheless, the effects of parental HC/LP diets on their progeny have never been investigated in species with higher trophic levels. This question, however, is of importance as, for aquaculture sustainability, the diets of these species are evolving rapidly in farming practices and could trend towards an increase in the carbohydrate/protein ratio. Protein-rich fishmeal, traditionally used to formulate the diets of aquaculture species with a high protein requirement [25], could be replaced by plant-derived carbohydrates in broodstock diet [26,27]. It is thus essential to evaluate the consequences of a parental HC/LP diet on progeny in species with higher trophic levels.

To address this question, two-year old male and female trout were fed either a control diet formulated with no carbohydrate and a high protein content (NC diet, 0% carbohydrate and 63.89% protein) or a diet containing 35% carbohydrate but with a lower protein content (HC/LP, 35% carbohydrate and 42.96% protein) for an entire reproductive cycle for females and for a period of 5 months for males. Cross-fertilizations were carried out in order to obtain four groups of fry. Different phenotypic traits (mass, length, and morphometric analyses) were recorded before their first feeding. The effects on metabolism were inves-
tigated by measuring the energetic status of fry and by comparing their transcriptomes. The global DNA methylation was also assessed.

**2. Material and Methods**

**2.1. *Experimental Design***

Two-year old male and female rainbow trout (*Oncorhynchus mykiss*) were distributed into two 8 m³ tanks at 8 °C (November—year 1) and fed since the resumption of feeding in the experimental INRAE facilities of Lees-Athas (Figure 1). Broodstock were fed either a control diet containing 63.89% protein and 0% carbohydrate (NC) or a diet containing 42.96% protein and 35.30% carbohydrate (HC/LP). Fishmeal was used as the protein source, fish oil was included as the lipid sources, and carbohydrates were provided by gelatinized starch (Table S1). These two diets were prepared in our own facilities (INRAE, Donzacq, France) as extruded pellets (BC45 BisVis Clextral, France). The feeding experiment lasted until the next reproduction for females, i.e., 10 months (October-year 2). Regarding males, a *Saprolegnia* sp. infection occurred in late April (year 2) in the tank with males fed the HC/LP diet, and it triggered some deaths. To keep a sufficient amount of fish alive for reproduction, the males were re-fed the NC diet (details in Callet et al. [27]).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental design. Broodstock females and males were fasted (F) after spawn from 16 Nov (November) to 21 Dec (December). Then, they were fed for one year either the NC diet (no-carbohydrate) or the HC/LP (high carbohydrate/low protein) diet (see details in Material and Methods Section). From the 25 of October (Oct) to the 9 of November, broodstock fish were fasted again. Cross-fertilizations were then carried out to obtained four groups of offspring: NN, NH, HN, and HH. Respirometry analyses and samplings were performed after hatching and just before the first feeding. Offspring were then fed for 3 weeks, and a second sampling was performed.</figcaption></figure>

During the spawning period (November—year 2), spawns from 3 NC females and 2 HC/LP females were cross-fertilized with milts from 4 males from each experimental condition (NC and HC/LP diets). Thus, offspring from 4 different conditions were obtained: NN, fry from both males and females fed the control diet; HN, fry from only females fed the HC/LP diet and males fed the control NC diet; NH, fry from only males fed the HC/LP diet and females fed the NC diet; and HH, fry from both parents fed the HC/LP diet (Figure 1).

Eggs hatched in December (year 2), from 44 to 48 days post-fertilization (dpf) (HH: from 44 to 46 dpf; HN: from 44 to 47 dpf; NH: from 45 to 48 dpf; and NN: from 46 to 48 dpf). Fry were divided in two distinct batches. To investigate the specific effects of the parental HC/LP on metabolism, some yolk-sac fry were transferred before complete resorption to an experimental INRAE facility (Saint-Pée-sur-Nivelle, France), and they were kept at 10 °C. Experiments were run before the first feeding between 61 and 64 dpf (description below). To monitor growth, some yolk-sac fry were transferred to another experimental
INRAE facility (Donzacq, France), where they were kept at 17 °C. At complete resorption, fry were fed a commercial diet (T3-P Omega, Skretting, France). After 3 weeks (Feb-year 3), fry were weighed and individually photographed from both sides.

## 2.2. *Metabolic Rate Assays*

At 61 dpf, over a 4-day period (January—year 3), 16 yolk-sac fry were randomly selected (n = 4 per condition) to measure their metabolic rates (MR). Individuals were kept for 24 h at 10 °C (±0.1 °C) under dim light (7 lx) for acclimatization and then placed in respirometry chambers of an intermittent flow respirometer (Loligo Systems, Viborg, Denmark) as previously described by Régnier [28]. Oxygen consumption as a proxy for MR was measured over an 18 h period at 10 °C. Six hours were removed from the data, as this period is considered the acclimation period. The standard metabolic rate (SMR) was estimated by calculating the oxygen consumption during the last three hours of measurement (3 slopes). For the routine metabolism (RMR), oxygen consumption was measured over 12 h after the acclimatization period (12 slopes). SMR and RMR are expressed in mm³·O₂·g⁻¹ of wet mass. After MR measurement, individuals were anaesthetized in a benzocaine bath at 30 mg·L⁻¹, weighed (BW), and then killed in a benzocaine bath at 60 mg·L⁻¹. A photo of each yolk-sac fry was taken from both sides. The samples were then individually frozen in liquid nitrogen and stored at -80 °C for further analyses.

## 2.3. *Fry Body Morphologies*

Body and yolk-sac measurements were performed for each fry (left and right side) by two different people using the software ImageJ (https://imagej.nih.gov/ij/index.html, accessed on 1 March 2020). To characterize fry body morphologies (before complete resorption and after 3 weeks of feeding), the fork length (L), the body depth (D) the head length (HL), the pre-orbital length (POL), and the whole-body area (WBA) were measured (Figure S1). Eye morphologies were also assessed by measuring the eye diameter (ED) of each fry. The relative head length (HLr), relative pre-orbital length (POLr), and relative eye diameter (EDr) were calculated by dividing HL, POL, and ED by fry length. Fulton's condition factor, K, was calculated as follows: $K = 100 \times BW/L^3$. To describe fry yolk-sac morphologies (i.e., yolk-sac utilization), the yolk-sac length (YSL), height (YSH), and area (YSA) were measured. Then, the relative yolk-sac area (YSAr), height (YSHr), and length (YSLr) were calculated by dividing YSA, YSH, and YSL by WBA, D, and L, respectively. The yolk-sac volume (YSV) was estimated as follows: $YSV = (\pi/6) \times YSL \times YSH^2$ [29].

## 2.4. *RNA and DNA Extraction*

Whole yolk-sac fry tissues were individually homogenized in ice using an Ultra Turrax homogeniser (T25 basic IKA-WERKE) at the speed setting of 4 for 10 s until they were homogeneous. Total RNA, DNA, and protein were extracted from whole yolk-sac fry (n = 16 per condition) using the QIAGEN AllPrep DNA/RNA/Prot Preparation Kit according to the manufacturer's recommendations (Qiagen). The concentration of extracted RNA was analyzed using a spectrophotometer (Nanodrop ND1000, LabTech) by measuring absorbance at 260 nm. The quality of RNAs was checked with a Bioanalyzer (Agilent Technologies, Kista, Sweden), and 12 samples for each condition were selected for further analyses according to the RIN (RNA integrity number).

## 2.5. *Determination of Fry Sex*

The extracted DNA was used to assess the sex of the fry (n = 12 per condition), as previously described in Yano et al. [30]. Briefly, PCR was performed to amplify the master sex-determining gene (sdY gene). Then, 1 µL of DNA was mixed with 1.25 µL of each primer at 10 µm (forward: CCCAGCACTGTTTTCTTGTCTCA; reverse: CTGTTGAAGAGCATCACAGGGTC), 1 of µL dNTP mixture (Ozyme, 4 × 10 mM), and 5 µL of 5xPCR Buffer (Promega) with 0.125 µL of Taq DNA polymerase (Promega, 5 U/µL) in a total volume of
25 µL. Thermal cycling consisted of denaturation for 20 s at 94 °C followed by 35 cycles of 94 °C for 20 s, 59 °C for 20 s, and 72 °C for 20 s, with a final extension of 5 min at 72 °C. PCR products were electrophoresed on a 2% agarose gel in a 1.5× TAE buffer and stained with SYBRSafe to reveal the presence or absence of sdY.

## 2.6. Global DNA Methylation

According to manufacturer recommendations, 1.5 µg of DNA (measured by Nanodrop) was digested with 0.5 µL of RNase cocktail (Ambion, Austin, TX, USA) containing 500 U/mL of RNase A and 20,000 U/mL of RNase T1 to avoid further nucleoside contamination. Digestion was conducted for 30 min at 37 °C. Nucleosides were obtained through a single step hydrolysis process using DNA Degradase Plus (Zymo Research, Orange, CA, USA) following the manufacturer's protocol. In brief, the reaction mix consisted of 1.5 µg DNA sample, 2.5 µL 10× DNA Degradase Plus reaction buffer, 1 µL DNA Degradase Plus (5 U/µL), and H₂O ultrapure up to a total of volume of 25 µL. The reaction mixes were incubated at 37 °C for 2 h followed by heat inactivation at 70 °C for 20 min. Samples were ultracentrifuged using 0.5 mL of Amicon Ultra at 3 kDa units (Merck Millipore, Billericia, MA, USA) according to the manufacturer's instructions.

HPLC-UV analysis was based on the protocol described in Kovatsi et al. [31] with major modifications. The system consisted of an Alliance 2695 separation module (Waters, Milford, MA, USA), a 2487 dual Absorbance Detector (Waters, Milford, MA, USA), and a column oven. Chromatographic separation was performed on a Luna C8 (3 µm, 100 × 3 mm) (Phenomenex, Torrance, CA, USA) column. The mobile phase was composed of the following: solvent A: 10 mM potassium phosphate buffer, pH = 5.9 ± 0.1; solvent B: 100% methanol. Linear gradient elution was employed as follows: 0, 100% A; 0–8 min, 90% A; 8–8.5 min, 73% A; 8.5–13.5 min, 65% A. The flow was set at 0.5 mL/min. The temperature of the column oven was 25°C. The wavelength of UV detection was 277 nm.

Identification and quantification of nucleosides were conducted with external standards purchased from Berry and Associates (Dexter, MI, USA). The monitored nucleosides were 2'-deoxycytidine (dC), 5-methyl-2'-deoxycytidine (5-mdC), 5-hydroxymethyl-2'-deoxycytidine (5-hmdC), 5-formyl-2'-deoxycytidine (5-fdC), and 5-carboethoxy-2'-deoxycytidine (5-cadC). The in vivo global level of 5-mC, 5-hmC, 5-fC, and 5-caC was calculated as the percentage of each individual's molar quantity divided by the total molar quantity of all the detected cytosine forms. Taking 5-mC as an example, the percentage of 5-mC was calculated using the following equation: 5-mC% = $100 \times \frac{Q5-mC}{QC + Q5-mC + Q5-hmC + Q5-fC + Q5-caC}$, where QC, Q5-mC, Q5-hmC, Q5-fC, and Q5-caC are the molar quantities of 5mC, 5hmC, 5-fC, and 5-caC, respectively.

## 2.7. Microarrays, cDNA Labelling and Hybridization

Transcriptome profiles of yolk-sac fry were analyzed with microarray technology (n = 12 per condition). Microarray analyses were performed on an RBT-specific Agilent-based microarray platform with 8 × 60 K probes per slide. Then, 150 ng of total RNA was first amplified by a reverse transcription using a polyDT T7 primer (denaturation step: 10 min at 65 °C; reaction step: 2 h at 40 °C; inactivation step: 5 min at 70 °C). The obtained cRNA were labelled with Cy3-dye (2 h at 40 °C). Excess dye was removed using a RNeasy kit (Qiagen). The level of dye incorporation was evaluated using a spectrophotometer (Nanodrop ND1000, LabTech; yield > 0.825 µg cRNA and specific activity > 6 pmol of Cy3 per µg of cRNA). A total of 600 ng of Cy3-cRNA was then fragmented with a specific buffer (30 min at 60 °C). Cy3-cRNAs were then hybridized on a sub-array (17 h at 65 °C in a microarray hybridization oven (Agilent). Slides were washed and scanned (Agilent DNA Microarray Scanner, Agilent Technologies, Massy, France) using the standard parameters for a gene expression 8 × 60 K oligoarray (3 µm and 20 bits). Data were then collected with the Agilent Feature Extraction software (10.7.1.1) and are available in GEO (accession number: GSE171310).
## 2.8. qPCR

Expressions of some genes selected in the transcriptomic analyses and genes from different pathways of interest were assessed by qPCR. Assays were carried out according to MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) standards [32]. For each sample (n = 12 per condition), 1 µg of total RNA per condition was reverse transcribed to cDNA with SuperScript III RNase H reverse transcriptase (Invitrogen, Carlsbad, CA, USA) using random dT Primers. qPCRs were performed using the Roche Lightcycler 480 system (Roche Diagnostics, Neuilly-sur-Seine, France). Then, 2 µL of diluted cDNA (2 µL cDNA diluted in 150 µL of water) was mixed with 3 µL of LightCycler 480 SYBR Green I Master Mix and diluted to obtain a final volume of 6 µL. Forward and reverse primers were used at a final concentration of 400 nM. Thermal cycling was initiated with an incubation at 95 °C for 10 min. Forty five steps of PCR were performed, each one consisting of heating at 95 °C for 15 s for denaturing and at 60 °C for 10 s for annealing as well as a third extension step at 72 °C for 15 s. Melting curves were systematically monitored (with a gradient of 0.5 °C/10 s from 55 °C to 94 °C) to confirm the specificity of the amplification reaction. Each PCR assay was run with replicate (duplicate of reverse transcription and PCR amplification) and negative controls included. PCR efficiency was measured by the slope of a standard curve using serial dilutions of cDNA (5 dilutions of a pool of all conditions from D20 to D380 in triplicate). Primer sequences and PCR efficiency values are presented in Table S2.

## 2.9. Statistical Analyses

All the statistical analyses were performed using R Software (version 3.2.5) [33]. Linear mixed-effects models were used to analyze the different parameters measured using the packages "lme4" from R software with a significance threshold set at *p*-value = 0.05. The best model was then selected using the Akaike Information Criterion (AIC). Concerning fry phenotypes, the effects of the maternal HC/LP diet and paternal HC/LP diet as well as the interaction between these two variables were investigated via the zootechnical parameters measured and the yolk-sac measurement. The sampling day and fry identification (except for BW and K) were treated as random effects. Then, the correlations between yolk-sac measurements and body measurements were estimated using a Pearson test (*cut-off* *p* = 0.01). When a significant correlation between two traits was present, ANCOVA analyses were used to test if the relationship between those traits was affected by the parental HC/LP diet.

To investigate the potential effect of the parental HC/LP diet on fry metabolism, the maternal nutritional history, the paternal nutritional history, and the interaction between these variables were investigated using the SMR and RMR. Then, the correlations between the SMR, RMR, yolk-sac, and body measurements were also estimated (Pearson test). When a significant correlation between two traits was present, ANCOVA analyses were used to test if the relationship between them was affected by the parental HC/LP diet.

Data from the microarray analysis were transformed with a logarithmic transformation, scale normalized, and analyzed using the Limma package [34]. In order to identify the differentially expressed genes resulting from the maternal, paternal, and parental HC/LP diet, transcriptomes of HN, NH, and HH fry were successively compared with the transcriptomes of NN fry. For these three comparisons, Limma t-tests were performed with a correction for multiple tests (*cut-off* *p* = 0.05 after a Benjamini-Hochberg correction), taking into account the sex of the fry and the time of the sampling. To characterize the change induced by the parental HC/LP diet, the gene ontologies of the probes that were annotated were collected from the DAVID (Database for Annotation, Visualization and Integrated Discovery bioinformatics resource, version 6.7) [35].

The relative expressions of genes assessed by qPCR were calculated by a mathematical method based on the real-time PCR efficiencies (Table S2) using a geometric mean from three reference genes (*eef1a*, *actb* and *18S*) for normalization [36]. The relative mRNA levels obtained were then transformed with a logarithmic transformation. To investigate the
potential effect of the parental HC/LP diet on these mRNA levels, the effects of the sex of the fry, the maternal nutritional history, the paternal nutritional history, and the interaction between these last two variables were investigated using linear mixed-effects models. In situations where a significant interaction was present, a Tukey post hoc test was carried out. The day of sampling was treated as a random effect. As previously described, AIC was used to select the best fitted model.

# 3. Results

## 3.1. *Phenotypes*

Fry phenotypes, including fry body mass, body, head, and yolk-sac morphologies, were recorded before their first feeding and three weeks later. Between 61 and 64 dpf (before first feeding), the whole yolk-sac fry weighed 94.91 ± 10.49 mg (Figure 2A). Except for the fry total length and the relative eye diameter (Figure 2A), no significant differences were detected among fry body morphologies (i.e., for BW, K, D, WBA, HL, HLr, POLr, and ED, see Table S3 for statistical details). NH and HH fry had a significantly higher total length than that of NN and HN fry (+1.5%, $\chi^2 = 4.0$, df = 1, *p* = 0.046). The relative eye diameters of HN, NH, and HH fry were, respectively, 2.55, 2.32, and 3.78% lower than those for NN fry. Regarding yolk-sac morphologies (YSH, YSHr, YSL, YSLr, VYS, YSA, and YSAr), no significant differences were detected among offspring (Table S3).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Effect of the parental HC/LP on fry body mass, length, relative head length, and relative eye diameter (A) before first feeding and (B) after 3 weeks of feeding. The effect of the paternal nutritional history is depicted in blue and the effect of both the maternal and the paternal nutritional history in grey ("**" means *p*-value < 0.01 and "*" means *p*-value < 0.05).</figcaption></figure>

After 3 weeks of feeding, fry weighed 362.41 ± 106.75 mg (Figure 2B). No significant differences were detected in fry body morphologies (i.e., for BW, K, L and D). However, the head proportions of NH and HH fry were on average 2.44% lower than those of HN and NN fry ($\chi^2 = 10.7$, df = 1, *p* = 0.001). Moreover, the eye proportions of HN, NH, and HH fry were 1.62, 1.93 and 4.94% lower, respectively, than those of NN fry ($\chi^2 = 4.9$, df = 1, *p* = 0.03) (Figure 2B).

Before the first feeding, phenotypical traits related to the body (BW, K, L, D, and WBA) and head (HL, HLr, POLr, and ED) morphologies were correlated with different yolk-sac traits (Figure 3A). All but three of these relationships remained unaffected by the
parental nutrition. The relationships linking BW and yolk-sac length, volume, and area were significantly different for HH fry in comparison to those for NN fry (Figure 3B). While NN fry body mass was positively correlated with YSH, VYS and YSA, such correlations did not exist for HH fry.

## 3.2. *Metabolic Rate*

Fry metabolic rates were measured between 61 and 64 dpf (Figure 4A). The HH and HN fry had on average a 6.45% and a 5.48% lower SMR and RMR, respectively, than those of NH and NN fry ($\chi^2 = 5.4, df = 1, p = 0.02$ and $\chi^2 = 4.3, df = 1, p = 0.04$, respectively). Moreover, neither SMR nor RMR were significantly correlated with yolk-sac measurements ($p > 0.05$).

## 3.3. *Global DNA Methylation*

The five intermediates of DNA methylation were measured in whole fry between 61 and 64 dpf (Figure 5A). The C proportion was significantly decreased in HN and HH in comparison to the C proportion observed in NN fry ($\chi^2 = 12.2, df = 1, p < 0.001$). The 5-mC proportion was 8.3% and 16.8% lower in HN and HH fry, respectively, in comparison to that observed in NN fry ($\chi^2 = 9.3, df = 1, p = 0.002$). The 5-hmC and 5-fC were not detected by HPLC-UV in our samples. Finally, the 5-caC proportion significantly increased in HN, NH, and HH fry in comparison to that observed in NN fry ($\chi^2 = 11.9, df = 1, p < 0.001$).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Relationships between growth and yolk-sac measurements. (A) Matrix of correlation between traits. Correlation was investigated using Pearson tests (*cut-off* *p* = 0.01), and each value is indicated inside the matrix, while significant positive correlations are represented in red significant negative ones in blue. (B) Relationships between fry mass and yolk-sac height, volume, and area according to parental dietary history.</figcaption></figure>
<figure><img src="image_7.png" /></figure>

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

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

**Figure 4.** Effect of parental HC/LP on fry metabolism and, in particular, (A) on fry metabolic rates, (B) relative mRNA levels of genes related to energy production, and (C) relative mRNA levels of genes involved in mitochondrial dynamics. The effect of maternal nutritional history is depicted in red and the effect of paternal nutritional history in blue ("*" means *p*-value < 0.05).

3.4. *Transcriptomes*

Between HN and NN fry, no significant difference in expressions of probes was found, while 15 probes were differentially expressed between NH and HH fry, and only the ankyrin repeat domain 16 had a fold change (FC) > 1.5 (Table S4). In total, 308 probes were differentially expressed between HH and NN fry, of which 8 had a FC > 1.5 (Table S4). Among these probes, 7.1% were not annotated. Regarding the annotated probes, the most represented pathways were those involved in metabolic pathways (18.2%), regulation of transcription (17.9%), immunity (8.8%), development (including genes coding for crystallin), and the signaling pathway.

Pathways related to metabolism included genes coding for glucose transporters and genes involved in glycolysis and amino acid metabolism, energy production, and mitochondrial dynamism. Hence, these pathways were selected to be further analyzed by qPCR, thereby allowing for the validation of microarray data (Figure S2). Thus, HN and HH fry had, on average, 13.6% and 14.6% higher *hk2* and *pfkla* mRNA levels, respectively, when compared with those of NH and NN fry. On average, NH and HH fry had significantly lower *glut1ba* (−14.8%), *glut1bb* (−9.4%), *pfkmaa* (−19.5%), *pfkmab* (−26.8%), *pfkmba* (−39.3%), and *pfkmbb* (−21.9%) mRNA levels than those of HN and NN fry. HN fry had significantly higher *ldhaa* mRNA levels than those of the fry from the three other conditions.
Results regarding genes coding for glucose transporters and the 6-phosphofructokinase are in accordance with results obtained by the transcriptomes analyses. In regard to amino acid catabolism, NH and HH had, on average, significantly lower *hibadh* mRNA levels than those of NN and HN fry (−18.6%). In contrast, no differences in the levels of mRNA for genes coding for the glutamine synthetase and glutaminase were detected among all conditions, which is in contrast with results obtained from the transcriptomic analyses. Regarding energy production, mRNA levels from one of the three homologs coding for the citrate synthase were decreased in HN and HH fry in comparison with those of NN and NH (−26.1%). Moreover, mRNA levels for one of the three homologs coding for 2-oxoglutarate dehydrogenase decreased by −11.1% in NH and HH fry when compared with the levels recorded for NN and HN (Figure 4B). Together, these results are in accordance with the transcriptomic analyses. Interestingly, the comparison of transcriptomes also revealed some effects of the parental HC/LP diet on mitochondrial dynamism (*bnip3* and *mfn2*). qPCR results did not confirm the results obtained for mitofusin 2 (*mfn2*). NH fry tended to have a lower *bnip3a* mRNA level than that of NN fry. Moreover, NH and HH fry had, on average, significantly lower *opa1a* (−10.4%), *fis* (−8.4%), and *parkin* (−9.1%) mRNA levels than those of NN and HN fry (Figure 4C).

Finally, the most significant result according to the transcriptomic analyses (Table S4) concerns a gene that codes for the DNA methyltransferase 3 (log2 fold change = −2.02). Among the 8 homologs known to code for this protein in trout, *dnmt3bba2* mRNA levels were, on average, significantly decreased in HN and HH fry in comparison to those of NN and NH (−15.9%). *dnmt3bbb* mRNA levels were also −15.3% lower in HN, NH, and HH fry in comparison to those of the control NN (Figure 5B). Finally, NH, HN, and HH fry also had −20.0% lower *crybbb* mRNA levels (gene coding for crystallin) than those of NN.

<figure><img src="image_10.png" /><figcaption>B</figcaption></figure>

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

**Figure 5.** Effect of parental HC/LP on fry epigenetic landscape. (A) Proportion of intermediates of DNA methylation in their offspring and (B) relative mRNA levels of two homologs coding for Dnmt3, an enzyme involved in the maintenance of DNA methylation patterns. The effect of the maternal nutritional history is depicted in red and the effect of both the maternal and the paternal nutritional history in grey ("*" means *p*-value < 0.05). Different letters indicate significant differences between groups, which were investigated with a Tukey post hoc test, in the case of a significant interaction between the paternal and the maternal nutritional history.

**4. Discussion**

In an effort to make aquaculture more sustainable, much focus has been placed on broodstock diets in which fishmeal has been replaced by terrestrial plants products, thus increasing their carbohydrate content [26,27]. However, parental HC diets are known to
significantly impact both the phenotypes of offspring as well as metabolism in mammals [2]. The effects of such parental diets have yet to be explored in higher trophic level teleost fish. The present study sought to evaluate the impact of a maternal and paternal HC/LP diet individually in addition to the combination of the maternal and paternal HC/LP diet, namely, the parental HC/LP diet, on yolk-sac fry. During this stage, fry switch from an endogenous to exogenous feeding and thus undergo important metabolic changes, which could reveal some changes triggered by the programming led by parental nutrition.

## 4.1. Offspring Growth Was Unaffected by the Parental HC/LP Diet

First, several measurements were carried out to detect any adverse effects of the parental HC/LP diet on offspring phenotypes, especially with respect to growth.

Neither the maternal nor the paternal nutritional history affected offspring body mass. During this early stage, the yolk-sac was not yet fully resorbed, and the mean mass included not only the mass of the fry but also the mass of the remaining yolk sac. As such, some effect of the parental diet on growth could have been hidden if body mass had been considered alone. Body morphology parameters, such as total length, depth, and whole body area, needed to be used as a proxy for growth. Among the different morphological measurements that depict growth, a slight increase was recorded only for the total fry length in response to the paternal HC/LP diet. However, this effect was small and was not observed after 3 weeks of feeding. Overall, the parental HC/LP diet did not significantly affect offspring growth in rainbow trout. This result is of particular interest as it is in direct contrast with results typically found in mammals [3,37–45] and birds [46] in which both maternal and paternal LP diets are known to compromise offspring growth.

In fish, in contrast to mammals, all of the nutrients available for an individual's development and growth before the first feeding are deposited in the egg, facilitating the study of programming. It has previously been demonstrated that female rainbow trout fed a challenging diet could maintain the quality of their eggs [26,47]. Such findings have been confirmed in the present experiment. Indeed, egg size, which reflects the quantity of nutrients available and is known to be correlated with the fry mass at first feeding in rainbow trout [48], was not significantly affected by the HC/LP diet [27]. Moreover, egg macronutrient content (glycogen, free glucose, lipid, and protein) was also maintained [27]. The unimpaired quality of eggs (nature of macronutrients and quantity provided) produced by females fed the HC/LP diet accounts for the fact that offspring growth was maintained. Furthermore, we also investigated the possible effect of the parental diet on the rate of yolk-sac resorption. The negative relationship between yolk-sac measurements, fry growth, and, in particular, head length, was expected, as the more the fry grows, the more its yolk-sac diminishes [49]. These relationships were not affected by the parental HC/LP nutrition. Thus, offspring could have used their reserves at a similar rate, regardless of their parental nutritional history, explaining why no differences in fry growth were detected at this early stage. Moreover, positive relationships were observed between fry mass-fry depth and yolk-sac measurements, which probably illustrates the fact that the biggest fry had the biggest yolk-sac, regardless of their resorption rate state. Interestingly, while there are positive relationships between BW and yolk-sac height, BW and volume, and BW and area for NN fry, such relationships did not exist for the HH fry. That is, while the heaviest NN fry had the biggest yolk-sac (area and volume), this was not the case for HH fry. A large variability in HH fry size along with differences in resorption rate amongst HH fry may explain this result. However, such variability in resorption rates was not related to a greater variability in metabolic rates (Figure 4) and should be further explored.

## 4.2. Eye Development Was Affected by the Parental HC/LP

Interestingly, while fry growth did not seem to be strongly affected by the parental nutritional history, both paternal and maternal HC/LP diets significantly decreased the relative eye diameter in fry. Interestingly, this effect persisted after 3 weeks of feeding. At both time intervals, a synergistic effect of the paternal and maternal diets appeared
as the HH fry had the lowest relative eye diameter. In addition, an important number of genes differentially expressed between the HH and the control NN fry were linked to neurological development pathways. Among them, relative mRNA levels for several genes coding for crystallin, the protein responsible for the optical properties of the lens [50], were significantly decreased in HH fry (log2FC = -0.5). Interestingly, in rainbow trout, eye size is correlated with age and not fish size, as eye growth is known to be maintained during a period of restricted feeding, while somatic growth is typically decreased [51].

Despite the fact that the reduction of protein in the parental HC/LP diet did not impair offspring growth, eye development could have been affected. The effect of a maternal LP diet on hormonal signals should be explored. Hormonal signals have been identified as a potential mechanism that can explain the imprinting of maternal nutrition. Of particular interest is cortisol, which is deposited in eggs and is also known to affect fry and eye development [52,53].

### 4.3. *Energy Metabolism Was Altered by the Parental HC/LP Diet*

To assess any adverse effects of the parental HC/LP diet on yolk-sac fry, metabolic rates were measured and non a priori analyses were used to detect any pathways affected.

The paternal HC/LP diet did not impair metabolic rates in offspring, and their transcriptomes were not highly affected. These results were not awaited. It is now well recognized that sperm also carry information that could highly affect their offspring metabolism in the long term [9] via changes in sperm epigenetic status [54]. Male broodstock could have been selected by the important mortality triggered by the *Saprolegnia* sp. [27]. Following this event, a dietary shift occurred during the feeding trial, and, as a consequence, male broodstock only received the HC/LP diet for a period of 5 months. Both the potential selection of males and the restricted exposure to HP/LC diet could explain the limited effect observed for offspring phenotypes and metabolism.

In contrast, while the maternal HC/LP diet induced a shift in offspring metabolic rates, no differences were detected at the molecular level (i.e., transcripts) between the HN and NN fry.

Finally, a synergistic effect of paternal and maternal nutritional histories was observed for metabolism. HH fry had lower metabolic rates than those of the NN control. Moreover, this result is in accordance with results obtained at the molecular level. The comparison of fry transcriptomes revealed that the HH fry metabolisms were more affected than the control NN were. The combined effect of the paternal and maternal HC/LP diet led to dysregulation of energy metabolism at the transcriptomic level, as seen by the overall downregulation of the expression of genes coding for proteins involved in energy production (glycolysis and TCA) in HH fry in comparison to NN fry. Previous studies reported a similar outcome in the muscle of offspring born from mothers fed LP diets in murine species [8,55,56].

Mitochondria, the powerhouse of the cell, are dynamic organelles that can either fuse to create more elongated mitochondria or a mitochondrial network, or undergo fission, thus creating a fragmented mitochondrial network. In murine species, maternal nutrient insults, such as a high sugar diet or a LP diet, typically induce unbalanced mitochondrial fusion/fission, leading to mitochondrial dysfunction in offspring across different generations [57–61]. In addition to the previous results regarding bioenergetics, qPCR also revealed that relative mRNA levels in genes coding for Bnip3 (log2FC = -0.13 between NH and NN), Opa1 (log2FC = -0.18 between HH and NN), and Fis (log2FC = -0.14 between HH and NN), proteins that orchestrate either mitochondrial fusion or fission, were slightly affected by the parental HC/LP diet. Interestingly, even though there is a synergistic effect of paternal and maternal HC/LP diets, the effect is primarily due to the paternal HC/LP diet, suggesting that the imprinting of the programming event is more likely due to epigenetic modulation than due to the maternal transmission of defect mitochondria [61].

Gene expression can vary greatly from one tissue to another and can be regulated in the opposite direction. This is particularly true for genes involved in mitochondrial mor-
phology and dynamics [62]. Moreover, the processes of mitochondrial fusion and fission are subjected to complex regulations through transcription and both post-transcription and post-translational regulation [62–64]. In the present study, due to the size of the fry at this early stage, analyses were conducted on whole fry. Despite being relative *parkin* mRNA levels, a gene coding for a protein involved in the mitophagy pathway, which can occur after mitochondrial fission, was also significantly reduced in NH and HH fry (log2FC = −0.14). It is impossible to confirm whether fission was promoted by the paternal HC/LP diet. Despite this limit, these preliminary results illustrated that a parental HC/LP diet could affect the balance between mitochondrial fusion and fission, which need to be tightly regulated to maintain an adequate mitochondrial function and, thus, energy metabolism in offspring.

Finally and more importantly, even though mitochondrial dynamics and energy
metabolism were altered by both the maternal and the paternal HC/LP diets at the tran-
scriptomic level and metabolic rates were disturbed, fry growth were not highly affected.
Moreover, no correlation between metabolic rates and growth parameters (BW, length,
depth, and WBA) was observed. Nevertheless, metabolic rates are known to be positively
correlated with growth in fish (i.e., fish with higher SMR grew faster when fed ad libi-
tum) [65]. This striking result suggests that the changes observed in metabolism were not
strong enough to compromise fry growth. However, as eye development was affected by
the parental HC/LP diet, it could be hypothesized that the parental HC/LP modified the
allocation of fry energy.

4.4. *The Parental HC/LP Diet Induced Global DNA Hypo-Methylation*

The inheritance of epigenetic information is now recognized as the main mechanism underlying the effect of nutritional programming, and among them is the modification of DNA methylation [13]. Here, we demonstrated that the maternal HC/LP diet induced global DNA hypomethylation in offspring, with a synergistic effect of the paternal HC/LP diet. The effect of alteration of the DNA methylation landscape and the modification of relative mRNA levels of genes involved in the methylation/demethylation process by the maternal HC/LP diet have previously been demonstrated in murine species [2]. While the biological role of this intermediate of the DNA demethylation pathway still needs to be assessed [66], both the paternal and the maternal HC/LP diets induced a significant augmentation of 5-caC. This provides further evidence that the HC/LP diet could have impact on the epigenetic landscape.

Different hypotheses have been proposed to explain how a parental HC/LP diet could alter the epigenetic landscape in offspring. First, it was previously shown that maternal HC/LP diets could directly affect *dnmt1* and *dnmt3* expression by alteration of the 1-carbon metabolism [12]. These genes code for the enzymes responsible for the maintenance of DNA methylation patterns during cell replication (Dnmt1) and de novo addition of a methyl group (Dnmt3). Such outcomes have been described in the liver of the offspring in mammals [67,68]. Interestingly, in the present study, the maternal HC/LP diet also induced an alteration of the expression of two homologs coding for Dnmt3. More specifically, relative *dnmt3* mRNA levels were significantly reduced in HH and HN fry (log2FC = -0.32 and -0.38). The potential different functions of these two homologs are still currently unknown [69], but these results are in line with the reduced proportion of 5-mC in these two groups.

LP/HC diets could also alter the availability of methyl donors or the cofactors needed for enzymes involved in methylation/demethylation processes, in which mitochondria play a central role. Mitochondria, whose role is not only restricted to energy production, also produce important cofactors required for enzymes involved in epigenetic modifications, such as alpha-ketoglutarate, which mediates DNA demethylation processes catalyzed by the ten-eleven translocation enzyme. Thus, the reprogramming of the mitochondrial function in an altered nutrient environment could also strongly impact the epigenetic
landscape [14]. For this reason, additional studies are needed to further investigate and describe the effect of the parental HC/LP diet on mitochondrial dynamics.

**5. Conclusions**

In rainbow trout, a teleost fish, we showed that a parental HC/LP diet affects offspring energy metabolism in the short term. However, growth was not compromised at this early stage in contrast to results typically observed in mammals. The carbohydrate/protein ratio could be further increased in the broostock diet, and plant-derived carbohydrates thus appear as an efficient substitute to fishmeal. However, it is necessary to test the effect of the parental diet in the long term, because the epigenetic landscape (global DNA methylation) was profoundly impacted, and some effects of the nutritional programming can only be revealed later in fish life, particularly during periods of stress. Finally, at both the molecular and phenotype levels, the effects triggered by the maternal HC/LP diet seemed to accumulate with the paternal effects, despite the fact that they were only fed for a period of 5 months. Thus, particular attention should be given to the synergistic effect existing between paternal and maternal nutritional histories.

**Supplementary Materials:** The following are available at https://www.mdpi.com/article/10.3390/biology10070585/s1, Figure S1: Descriptions of the different body measurements. D: body depth; HL: head length; POL: the pre-orbital length; ED: eye diameter; YSL: yolk-sac length, and YSH: yolk-sac height. Figure S2: Log2 fold changes obtained with microarray (purple) and with qPCR (grey) analyses. Table S1: Diet composition and fatty acid profile (% total FA). Table S2: List of the primers used for qPCR analyses. Table S3: Results of the linear mixed-effects used to investigate the effect of the sex of the fry and the parental HC/LP diet on the phenotypical traits recorded, SMR/RMR, methylation pattern, and gene expression. For each trait, the chosen model (based on AIC) is indicated in bold. Table S4: Results of the transcriptomic analyses: list of probes differentially expressed between NH and HH and HH and HH fry in addition to their fold change (FC) and the corrected *p*-value obtained by Limma *t*-tests.

**Author Contributions:** Funding acquisition by L.M. Conceptualization and writing by L.M., S.P. and V.B. Fish were reared by P.M., N.T. and Y.M. Investigation by V.B. (respirometry), P.C. (respirometry), S.G. (picture acquisition), C.H. (global DNA methylation), H.L. (picture analyses), T.C. (picture analyses and molecular analyses) and C.Z. (molecular analyses). Formal analyses, visualization, and writing by T.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the French National Research Agency (ANR-18-CE20-0012-01 "SweetSex").

**Institutional Review Board Statement:** The study was conducted according to the guiding principles for the use and care of laboratory animals and in compliance with French and European regulations on animal welfare (Décret 2001-464, 29 May 2001 and Directive 2010/63/EU, respectively). This protocol and the project as a whole were approved by the French National Consultative Ethics Committee (reference numbers 201610061056842).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are available on the GEO website under the accession number: GSE171310.

**Acknowledgments:** We thank Anthony Lanuque, Frédéric Terrier, and Frank Sandres for diet preparation.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Hoffman, D.J.; Reynolds, R.M.; Hardy, D.B. Developmental origins of health and disease: Current knowledge and potential mechanisms. *Nutr. Rev.* **2017**, *75*, 951–970. [CrossRef]

2. Guo, T.; Luo, F.; Lin, Q. You are affected by what your parents eat: Diet, epigenetics, transgeneration and intergeneration. *Trends Food Sci. Technol.* **2020**, *100*, 248–261. [CrossRef]

3. Watkins, A.J.; Sinclair, K.D. Paternal low protein diet affects adult offspring cardiovascular and metabolic function in mice. *Am. J. Physiol. Heart Circ. Physiol.* **2014**, *306*, H1444–H1452. [CrossRef]
4. Jahan-Mihan, A.; Rodriguez, J.; Christie, C.; Sadeghi, M.; Zerbe, T. The role of maternal dietary proteins in development of metabolic syndrome in offspring. *Nutrients* **2015**, *7*, 9185-9217. [CrossRef]

5. Oke, S.L.; Hardy, D.B. *Effects of Protein Deficiency on Perinatal and Postnatal Health Outcomes*; Physiology and Pharmacology Publications; Western Libraries: London, ON, Canada, 2017. [CrossRef]

6. Martin-Gronert, M.; Ozanne, S.E. Experimental IUGR and later diabetes. *J. Intern. Med.* **2007**, *261*, 437-452. [CrossRef] [PubMed]

7. Langley-Evans, S.C. Nutritional programming of disease: Unravelling the mechanism. *J. Anat.* **2009**, *215*, 36-51. [CrossRef]

8. Mortensen, O.H.; Olsen, H.L.; Frandsen, L.; Nielsen, P.E.; Nielsen, F.C.; Grunnet, N.; Quistorff, B. A maternal low protein diet has pronounced effects on mitochondrial gene expression in offspring liver and skeletal muscle; protective effect of taurine. *J. Biomed. Sci.* **2010**, *17*, S38. [CrossRef]

9. Carone, B.R.; Fauquier, L.; Habib, N.; Shea, J.M.; Hart, C.E.; Li, R.; Bock, C.; Li, C.; Gu, H.; Zamore, P.D.; et al. Paternally induced transgenerational environmental reprogramming of metabolic gene expression in mammals. *Cell* **2010**, *143*, 1084-1096. [CrossRef] [PubMed]

10. Morgan, H.L.; Paganopoulou, P.; Akhtar, S.; Urquhart, N.; Philomin, R.; Dickinson, Y.; Watkins, A.J. Paternal diet impairs F1 and F2 offspring vascular function through sperm and seminal plasma specific mechanisms in mice. *J. Physiol. 2020*, *598*, 699-715. [CrossRef]

11. Safi-Stibler, S.; Gabory, A. Epigenetics and the Developmental Origins of Health and Disease: Parental environment signalling to the epigenome, critical time windows and sculpting the adult phenotype. *Semin. Cell Dev. Biol.* **2020**, *97*, 172-180. [CrossRef]

12. Burdge, G.C.; Hanson, M.A.; Slater-Jefferies, J.L.; Lillycrop, K.A. Epigenetic regulation of transcription: A mechanism for inducing variations in phenotype (fetal programming) by differences in nutrition during early life? *Br. J. Nutr.* **2007**, *97*, 1036-1046. [CrossRef] [PubMed]

13. Zhang, N. Epigenetic modulation of DNA methylation by nutrition and its mechanisms in animals. *Anim. Nutr.* **2015**, *1*, 144-151. [CrossRef] [PubMed]

14. Harvey, A.J. Mitochondria in early development: Linking the microenvironment, metabolism and the epigenome. *Reproduction* **2019**, *157*, R159-R179. [CrossRef] [PubMed]

15. Hou, Z.; Fuiman, L.A. Nutritional programming in fishes: Insights from mammalian studies. *Rev. Fish Biol. Fish.* **2020**, *30*, 67-92. [CrossRef]

16. Gotoh, T. Potential of the application of epigenetics in animal production. *Anim. Prod. Sci.* **2015**, *55*, 145-158. [CrossRef]

17. Geurden, I.; Aramendi, M.; Zambonino-Infante, J.; Panserat, S. Early feeding of carnivorous rainbow trout (*Oncorhynchus mykiss*) with a hyperglucidic diet during a short period: Effect on dietary glucose utilization in juveniles. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* **2007**, *292*, R2275-R2283. [CrossRef]

18. Geurden, I.; Borchert, P.; Balasubramanian, M.N.; Schrama, J.W.; Dupont-Nivet, M.; Quillet, E.; Kaushik, S.J.; Panserat, S.; Médale, F. The positive impact of the early-feeding of a plant-based diet on its future acceptance and utilization in rainbow trout. *PLoS ONE* **2013**, *8*, e83162. [CrossRef]

19. Geurden, I.; Mennigen, J.; Plagnes-Juan, E.; Veron, V.; Cerezo, T.; Mazurais, D.; Zambonino-Infante, J.; Gatesoupe, J.; Skiba-Cassy, S.; Panserat, S. High or low dietary carbohydrate: Protein ratios during first-feeding affect glucose metabolism and intestinal microbiota in juvenile rainbow trout. *J. Exp. Biol.* **2014**, *217*, 3396-3406. [CrossRef]

20. Balasubramanian, M.N.; Panserat, S.; Dupont-Nivet, M.; Quillet, E.; Montfort, J.; Le Cam, A.; Medale, F.; Kaushik, S.J.; Geurden, I. Molecular pathways associated with the nutritional programming of plant-based diet acceptance in rainbow trout following an early feeding exposure. *BMC Genomics* **2016**, *17*, 449. [CrossRef]

21. Séité, S.; Masagounder, K.; Heraud, C.; Véron, V.; Marandel, L.; Panserat, S.; Seiliez, I. Early feeding of rainbow trout (*Oncorhynchus mykiss*) with methionine-deficient diet over a 2 week period: Consequences for liver mitochondria in juveniles. *J. Exp. Biol.* **2019**, *222*, jeb203687. [CrossRef]

22. Seiliez, I.; Vélez, E.J.; Lutfi, E.; Dias, K.; Plagnes-Juan, E.; Marandel, L.; Panserat, S.; Geurden, I.; Skiba-Cassy, S. Eating for two: Consequences of parental methionine nutrition on offspring metabolism in rainbow trout (*Oncorhynchus mykiss*). *Aquaculture* **2017**, *471*, 80-91. [CrossRef]

23. Lazzarotto, V.; Corraze, G.; Larroquet, L.; Mazurais, D.; Médale, F. Does broodstock nutritional history affect the response of progeny to different first-feeding diets? A whole-body transcriptomic study of rainbow trout alevins. *Br. J. Nutr.* **2016**, *115*, 2079-2092. [CrossRef]

24. Wischhusen, P.; Parailloux, M.; Geraert, P.A.; Briens, M.; Bueno, M.; Mounicou, S.; Bouyssiere, B.; Prabhu, P.A.J.; Kaushik, S.J.; Fauconneau, B.; et al. Effect of dietary selenium in rainbow trout (*Oncorhynchus mykiss*) broodstock on antioxidant status, its parental transfer and oxidative status in the progeny. *Aquaculture* **2019**, *507*, 126-138. [CrossRef]

25. Jobling, M. *National Research Council (NRC): Nutrient Requirements of Fish and Shrimp*; Springer: Berlin/Heidelberg, Germany, 2012.

26. Washburn, B.S.; Frye, D.J.; Hung, S.S.; Doroshov, S.I.; Conte, F.S. Dietary effects on tissue composition, oogenesis and the reproductive performance of female rainbow trout (*Oncorhynchus mykiss*). *Aquaculture* **1990**, *90*, 179-195. [CrossRef]

27. Callet, T.; Hu, H.; Larroquet, L.; Surget, A.; Liu, J.; Plagnes-Juan, E.; Maunas, P.; Turonnet, N.; Mennigen, J.A.; Bobe, J.; et al. Exploring the impact of a low-protein high-carbohydrate diet in mature broodstock of a glucose-intolerant teleost, the rainbow trout. *Front. Physiol.* **2020**, *11*, 303. [CrossRef]

28. Régnier, T.; Bolliet, V.; Labonne, J.; Gaudin, P. Assessing maternal effects on metabolic rate dynamics along early development in brown trout (*Salmo trutta*): An individual-based approach. *J. Comp. Physiol. B* **2010**, *180*, 25-31. [CrossRef] [PubMed]
29. Kamler, E. Resource allocation in yolk-feeding fish. *Rev. Fish Biol. Fish.* **2008**, *18*, 143. [CrossRef]

30. Yano, A.; Guyomard, R.; Nicol, B.; Jouanno, E.; Quillet, E.; Klopp, C.; Cabau, C.; Bouchez, O.; Fostier, A.; Guiguen, Y. An immune-related gene evolved into the master sex-determining gene in rainbow trout, Oncorhynchus mykiss. *Curr. Biol.* **2012**, *22*, 1423–1428. [CrossRef]

31. Kovatsi, L.; Fragou, D.; Samanidou, V.; Njau, S.; Kouidou, S.; Bailey, A. Evaluation of 5-methyl-2'-deoxycytidine stability in hydrolyzed and nonhydrolyzed DNA by HPLC–UV. *Bioanalysis* **2012**, *4*, 367–372. [CrossRef] [PubMed]

32. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. *Clin. Chem.* **2009**, *55*, 611–622. [CrossRef]

33. R Core Team. *R: A Language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2013.

34. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* **2015**, *43*, e47. [CrossRef]

35. Sherman, B.T.; Lempicki, R.A.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nat. Protoc.* **2009**, *4*, 44.

36. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. *Nucleic Acids Res.* **2001**, *29*, e45. [CrossRef]

37. Dahri, S.; Snoeck, A.; Reusens-Billen, B.; Remacle, C.; Hote, J.J. Islet function in offspring of mothers on low-protein diet during gestation. *Diabetes* **1991**, *40*, 115–120. [CrossRef] [PubMed]

38. Ozanne, S.; Wang, C.; Coleman, N.; Smith, G. Altered muscle insulin sensitivity in the male offspring of protein-malnourished rats. *Am. J. Physiol. Endocrinol. Metab.* **1996**, *271*, E1128–E1134. [CrossRef]

39. Rees, W.D.; Hay, S.M.; Brown, D.S.; Antipatis, C.; Palmer, R.M. Maternal protein deficiency causes hypermethylation of DNA in the livers of rat fetuses. *J. Nutr.* **2000**, *130*, 1821–1826. [CrossRef] [PubMed]

40. Petry, C.; Ozanne, S.; Wang, C.; Hales, C. Effects of early protein restriction and adult obesity on rat pancreatic hormone content and glucose tolerance. *Horm. Metab. Res.* **2000**, *32*, 233–239. [CrossRef] [PubMed]

41. Zambrano, E.; Martinez-Samayoa, P.; Bautista, C.; Deas, M.; Guillen, L.; Rodriguez-Gonzalez, G.; Guzman, C.; Larrea, F.; Nathanielsz, P. Sex differences in transgenerational alterations of growth and metabolism in progeny (F2) of female offspring (F1) of rats fed a low protein diet during pregnancy and lactation. *J. Physiol.* **2005**, *566*, 225–236. [CrossRef]

42. Jousse, C.; Parry, L.; Lambert-Langlais, S.; Maurin, A.C.; Averous, J.; Bruhat, A.; Carraro, V.; Tost, J.; Letteron, P.; Chen, P.; et al. Perinatal undernutrition affects the methylation and expression of the leptin gene in adults: Implication for the understanding of metabolic syndrome. *FASEB J.* **2011**, *25*, 3271–3278. [CrossRef]

43. Sohi, G.; Marchand, K.; Revesz, A.; Arany, E.; Hardy, D.B. Maternal protein restriction elevates cholesterol in adult rat offspring due to repressive changes in histone modifications at the cholesterol 7 α-hydroxylase promoter. *Mol. Endocrinol.* **2011**, *25*, 785–798. [CrossRef]

44. Zheng, S.; Rollet, M.; Pan, Y.X. Protein restriction during gestation alters histone modifications at the glucose transporter 4 (GLUT4) promoter region and induces GLUT4 expression in skeletal muscle of female rat offspring. *J. Nutr. Biochem.* **2012**, *23*, 1064–1071. [CrossRef]

45. Jia, Y.; Cong, R.; Li, R.; Yang, X.; Sun, Q.; Parvizi, N.; Zhao, R. Maternal low-protein diet induces gender-dependent changes in epigenetic regulation of the glucose-6-phosphatase gene in newborn piglet liver. *J. Nutr.* **2012**, *142*, 1659–1665. [CrossRef] [PubMed]

46. Rao, K.; Xie, J.; Yang, X.; Chen, L.; Grossmann, R.; Zhao, R. Maternal low-protein diet programmes offspring growth in association with alterations in yolk leptin deposition and gene expression in yolk-sac membrane, hypothalamus and muscle of developing Langshan chicken embryos. *Br. J. Nutr.* **2009**, *102*, 848–857. [CrossRef]

47. Lazzarotto, V.; Corraze, G.; Leprevost, A.; Quillet, E.; Dupont-Nivet, M.; Médale, F. Three-year breeding cycle of rainbow trout (Oncorhynchus mykiss) fed a plant-based diet, totally free of marine resources: Consequences for reproduction, fatty acid composition and progeny survival. *PLoS ONE* **2015**, *10*, e0117609. [CrossRef] [PubMed]

48. Springate, J.; Bromage, N. Effects of egg size on early growth and survival in rainbow trout (*Salmo gairdneri* Richardson). *Aquaculture* **1985**, *47*, 163–172. [CrossRef]

49. Yang, G.; Wang, Z.; He, J.; Li, W.; Yuan, D. Development and Allometry Patterns of Fine Scale Fish Larvae at Low Temperature. *J. Phys. Conf. Ser.* **2020**, *1575*, 012202. [CrossRef]

50. Jonasova, K.; Kozmik, Z. Eye evolution: Lens and cornea as an upgrade of animal visual system. *Semin. Cell Dev. Biol.* **2008**, *19*, 71–81. [CrossRef]

51. Pankhurst, N.; Montgomery, J. Uncoupling of visual and somatic growth in the rainbow trout Oncorhynchus mykiss. *Brain Behav. Evol.* **1994**, *44*, 149–155. [CrossRef]

52. Gagliano, M.; McCormick, M.I. Hormonally mediated maternal effects shape offspring survival potential in stressful environments. *Oecologia* **2009**, *160*, 657–665. [CrossRef]

53. Nesan, D.; Vijayan, M.M. Maternal cortisol mediates hypothalamus-pituitary-interrenal axis development in zebrafish. *Sci. Rep.* **2016**, *6*, 22582. [CrossRef] [PubMed]

54. Rando, O.J. Daddy issues: Paternal effects on phenotype. *Cell* **2012**, *151*, 702–708. [CrossRef]
55. de Brito Alves, J.L.; Toscano, A.E.; da Costa-Silva, J.H.; Vidal, H.; Leandro, C.G.; Pirola, L. Transcriptional response of skeletal muscle to a low protein perinatal diet in rat offspring at different ages: The role of key enzymes of glucose-fatty acid oxidation. *J. Nutr. Biochem.* **2017**, *41*, 117-123. [CrossRef] [PubMed]

56. de Oliveira Lira, A.; de Brito Alves, J.L.; Fernandes, M.P.; Vasconcelos, D.; Santana, D.F.; da Costa-Silva, J.H.; Morio, B.; Leandro, C.G.; Pirola, L. Maternal low protein diet induces persistent expression changes in metabolic genes in male rats. *World J. Diabetes* **2020**, *11*, 182. [CrossRef]

57. Wu, L.L.; Russell, D.L.; Wong, S.L.; Chen, M.; Tsai, T.S.; St John, J.C.; Norman, R.J.; Febbraio, M.A.; Carroll, J.; Robker, R.L. Mitochondrial dysfunction in oocytes of obese mothers: Transmission to offspring and reversal by pharmacological endoplasmic reticulum stress inhibitors. *Development* **2015**, *142*, 681-691. [CrossRef] [PubMed]

58. Saben, J.L.; Boudoures, A.L.; Asghar, Z.; Thompson, A.; Drury, A.; Zhang, W.; Chi, M.; Cusumano, A.; Scheaffer, S.; Moley, K.H. Maternal metabolic syndrome programs mitochondrial dysfunction via germline changes across three generations. *Cell Rep.* **2016**, *16*, 1-8. [CrossRef]

59. Cardenas-Perez, R.E.; Fuentes-Mera, L.; de la Garza, A.L.; Torre-Villalvazo, I.; Reyes-Castro, L.A.; Rodriguez-Rocha, H.; Garcia-Garcia, A.; Corona-Castillo, J.C.; Tovar, A.R.; Zambrano, E.; et al. Maternal overnutrition by hypercaloric diets programs hypothalamic mitochondrial fusion and metabolic dysfunction in rat male offspring. *Nutr. Metab.* **2018**, *15*, 38. [CrossRef] [PubMed]

60. Santana, D.F.F.; Ferreira, D.S.; Braz, G.; Sousa, S.; Araujo, T.; Gomes, D.; Fernandes, M.P.; Andrade-Da-Costa, B.L.D.; Lagranha, C.J. Maternal Protein Restriction in Two Successive Generations Impairs Mitochondrial Electron Coupling in the Progeny's Brainstem of Wistar Rats From Both Sexes. *Front. Neurosci.* **2019**, *13*, 203. [CrossRef]

61. Ferey, J.L.; Boudoures, A.L.; Reid, M.; Drury, A.; Scheaffer, S.; Modi, Z.; Kovacs, A.; Pietka, T.; DeBosch, B.J.; Thompson, M.D.; et al. A maternal high-fat, high-sucrose diet induces transgenerational cardiac mitochondrial dysfunction independently of maternal mitochondrial inheritance. *Am. J. Physiol. Heart Circ. Physiol.* **2019**, *316*, H1202-H1210. [CrossRef]

62. Wai, T.; Langer, T. Mitochondrial dynamics and metabolic regulation. *Trends Endocrinol. Metab.* **2016**, *27*, 105-117. [CrossRef]

63. Pernas, L.; Scorrano, L. Mito-morphosis: Mitochondrial fusion, fission, and cristae remodeling as key mediators of cellular function. *Annu. Rev. Physiol.* **2016**, *78*, 505-531. [CrossRef] [PubMed]

64. Sebastián, D.; Zorzano, A. Mitochondrial dynamics and metabolic homeostasis. *Curr. Opin. Physiol.* **2018**, *3*, 34-40. [CrossRef]

65. Salin, K.; Auer, S.K.; Rey, B.; Selman, C.; Metcalfe, N.B. Variation in the link between oxygen consumption and ATP production, and its relevance for animal performance. *Proc. R. Soc. B Biol. Sci.* **2015**, *282*, 20151028. [CrossRef] [PubMed]

66. Hofer, A.; Liu, Z.J.; Balasubramanian, S. Detection, structure and function of modified DNA bases. *J. Am. Chem. Soc.* **2019**, *141*, 6420-6429. [CrossRef]

67. Lillycrop, K.A.; Slater-Jefferies, J.L.; Hanson, M.A.; Godfrey, K.M.; Jackson, A.A.; Burdge, G.C. Induction of altered epigenetic regulation of the hepatic glucocorticoid receptor in the offspring of rats fed a protein-restricted diet during pregnancy suggests that reduced DNA methyltransferase-1 expression is involved in impaired DNA methylation and changes in histone modifications. *Br. J. Nutr.* **2007**, *97*, 1064-1073. [PubMed]

68. Gong, L.; Pan, Y.X.; Chen, H. Gestational low protein diet in the rat mediates Igf2 gene expression in male offspring via altered hepatic DNA methylation. *Epigenetics* **2010**, *5*, 619-626. [CrossRef]

69. Liu, J.; Hu, H.; Panserat, S.; Marandel, L. Evolutionary history of DNA methylation related genes in chordates: New insights from multiple whole genome duplications. *Sci. Rep.* **2020**, *10*, 970. [CrossRef] [PubMed]