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

HAL
open science

# Evolution of oxidative stress markers in livers of ducks during force-feeding

Herve Remignon, Pierre Burgues

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

## To cite this version:

Herve Remignon, Pierre Burgues. Evolution of oxidative stress markers in livers of ducks during force-feeding. Scientific Reports, 2023, 13 (1), pp.1046. 10.1038/s41598-022-27271-y. hal-03956709

HAL Id: hal-03956709

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

Submitted on 25 Jan 2023

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_3.png" /></figure>

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

OPEN Evolution of oxidative stress markers in livers of ducks during force-feeding

Herve Remignon<fnref n="1" />,²✉ & Pierre Burgues<fnref n="1" />

Mule ducks have been force-fed to develop a hepatic steatosis, also called “foie gras”, which is similar to the non-alcoholic fatty liver disease (NAFLD) described in humans and mammals. However, in hepatic steatosis resulting from force-feeding of ducks, very little is known about the fine biochemical events that occur due to the enormous and very rapid increase in total lipids that mainly accumulate in hepatocytes. To begin to reduce this lack of knowledge associated with the development of this specific hepatic steatosis, liver samples were taken at different times to follow the overall biochemical transformation of the liver as well as different markers of oxidative stress, hypoxia and apoptosis. The results indicate that the lipid content increases rapidly in the liver throughout the force-feeding period while the protein content decreases. The amount of hydroxyproline remains constant indicating that no liver fibrosis develops during the force-feeding period. On the contrary, all the tested biomarkers of cellular oxidative stress increase rapidly but without any visible disorder in the coordination of paired activities. At the same time, hypoxia-inducible factors also increase indicating that a hypoxia situation is gradually occurring in hepatocytes. This leads, in addition to the lipotoxicity induced by the accumulation of lipids, to an increased number of liver cells to enter into apoptosis. A relative variability in the level of these cellular responses was also observed indicating that, probably, certain animals support the development of this steatosis differently. This leads us to imagine that the physiological status of these birds may differ widely for reasons that remain to be clarified.

Abbreviations

<table><tr><td>ANOVA</td><td>ANalysis Of VAriance</td></tr><tr><td>B</td><td>Group of ducks at baseline, before the beginning of force-feeding</td></tr><tr><td>BFF</td><td>Group of BFF (Beginning of Force-Feeding) ducks</td></tr><tr><td>Cat</td><td>Catalase</td></tr><tr><td>CV</td><td>Coefficient of variation</td></tr><tr><td>DM</td><td>Dry matter</td></tr><tr><td>EFF</td><td>Group of EFF (End of Force-Feeding) ducks</td></tr><tr><td>ELISA</td><td>Enzyme-linked ImmunoSorbent assay</td></tr><tr><td>FFA</td><td>Free fatty acids</td></tr><tr><td>GSH</td><td>Reduced glutathione</td></tr><tr><td>GSSG</td><td>Oxidized glutathione</td></tr><tr><td>GSH/GSSG</td><td>Ratio GSH/GSSG</td></tr><tr><td>GPX1</td><td>Glutathione peroxidase 1</td></tr><tr><td>HIF1α</td><td>Hypoxia inducible factor 1 alpha</td></tr><tr><td>HIF2α</td><td>Hypoxia inducible factor 2 alpha</td></tr><tr><td>NAFLD</td><td>Non-alcoholic fatty liver disease</td></tr><tr><td>NASH</td><td>Non-alcoholic SteatoHepatitis</td></tr><tr><td>NQO1</td><td>Nicotinamide adenine dinucleotide (phosphate) hydrogen dehydrogenase quinone-1</td></tr><tr><td>OH-Pro</td><td>Hydroxy-PROline</td></tr><tr><td>Proc GLM</td><td>General linear model procedure</td></tr><tr><td>RMSE</td><td>Root mean square error</td></tr><tr><td>ROS</td><td>Reactive oxygen species</td></tr><tr><td>SD</td><td>Standard deviation</td></tr></table>

<footnote n="1">¹Toxalim (Research Centre in Food Toxicology), Université de Toulouse, INRAE, ENVT, UPS, 31300 Toulouse, France. ²INP-ENSAT, Université de Toulouse, 31320 Castanet-Tolosan, France.✉email: herve.remignon@toulouse-inp.fr</footnote>
<table><tr><td>SOD</td><td>Super oxide dismutase</td></tr><tr><td>TBARS</td><td>Thio barbituric acid reactive substances</td></tr></table>

The production of “foie gras” is traditional in France and is typical of French gastronomy which has been included on UNESCO’s list of Intangible Heritage of Humanity since 2010. To obtain this food product, waterfowl (either geese or ducks) are force-fed until they develop a liver steatosis. The metabolic mechanisms associated with the development of this hepatic steatosis are very similar to those described in human as non-alcoholic fatty liver disease (NAFLD) and/or non-alcoholic steatohepatitis (NASH). In short, it is the daily consumption of large quantities of carbohydrates that places the liver of force-fed ducks in a permanent state of synthesis of lipids which are only partially exported into the peripheral tissues, mainly the abdominal and subcutaneous adipose tissues. Numerous studies¹⁻⁶ have been conducted to study the development of this hepatic steatosis in force-fed ducks and mainly indicate that the de novo lipogenesis is dramatically increased while other hepatic metabolisms are also modified (proteolysis and proteogenesis, lipids exportation, carbohydrates metabolism, and programmed cell deaths, etc.). However, none have focused on the oxidative stress that could be at the origin of those large metabolic dysfunctions of the hepatic cells.

Even in humans, NAFLD remains a complex and multifactorial disease associated with a myriad of genetic, epigenetic, and environmental factors, and its pathogenesis is still not entirely clear⁷. Hepatic steatosis increases liver lipogenesis but hinders the breakdown of free fatty acids (FFA). Therefore, lipid accumulation sensitizes the liver to induce inflammation and cell death by pathogenic aggression that promotes oxidative stress, culminating in NASH and fibrosis. Oxidative stress is thus considered to be the main contributor to liver damage and disease progression in NAFLD⁸. According to⁹, this first hypothesis has been subsequently revised in a “multiple parallel-hit” model in which, in the presence of a significant accumulation of fat in hepatocytes and systemic and hepatic insulin resistance, multiple simultaneous alterations lead to an imbalance between the antilipotoxic protection system of the liver (mitochondrial β-oxidation) and the free radical production in gut and adipose tissue, resulting in endoplasmic stress, oxidative stress, and hepatocyte apoptosis.

For the present study, we induced the development of hepatic steatosis in ducks by force-feeding them for 10 days. Then, we measured throughout the development of steatosis several hepatic biomarkers to follow its biochemical transformation possibly associated with a restricted supply of oxygen, the development of cellular oxidative stress and finally its possible transition to an apoptotic state. The objective of the study was therefore to draw up an initial overview of the oxidative stress that can occur in the hepatocytes of force-fed ducks and to possibly associate it with the high accumulation of lipids associated with a restricted supply of oxygen and resulting in an increase in apoptotic events.

## Material and methods

**Animals and livers sampling.** A flock of about 1000 male mule ducks (Caïrina moschata x Anas platyrhynchos) was reared for 12 weeks, from hatching, according to standard commercial rules. At 12 weeks of age, 12 birds (live weights similar to the whole flock mean live weight) were randomly selected and slaughtered to constitute the group of ducks at baseline (B) before the beginning of force-feeding. The remaining birds entered in the force-feeding program for 20 meals (twice a day during 10 days as farmers usually do) according to a standard force-feeding program based on moistened corn flour (97.5% of corn). At the beginning of the force-feeding period, ducks received a quantity of 225 g of dry flour/meal, which was progressively increased to a final value of 510 g for the last meal. As a whole, during the force-feeding period, ducks ingested an average quantity of 7.0 kg of feed. After 5 meals of force-feeding, 12 animals were randomly (similarly than described upper for B ducks) selected and slaughtered to constitute the group of BFF (Beginning of Force-Feeding) ducks. The remaining birds continued the force-feeding program for 15 supplementary meals. At that time, all the birds were slaughtered and 24 of them were selected to constitute the group of EFF (End of Force-Feeding) ducks. At this final step, because that data was available for the whole flock, the selection of the EFF birds was done to represent the observed variability in fatty livers weights (mean ± SD = 571 ± 114 g, CV% = 20). At each step, (B, BFF, EFF) ducks were slaughtered approximately 11 h hours after the last meal in a commercial slaughter-house according to its standardized operations of slaughtering (electronarcosis, bleeding, scalding, plucking). At the end of the slaughtering line, 20 min post mortem, livers were harvested and weighed. Then, 50 g of tissue were collected from the median lobe and directly frozen in liquid nitrogen before storage at -80 °C.

All the biochemical measures were performed in duplicate after the grinding of the tissue in liquid nitrogen.

**Gross biochemical composition of livers.** The dry matter (DM) content was determined by desiccation of crushed livers in an oven at 105 °C for 24 h¹⁰. The total lipids content was measured according to¹¹ after extraction with chloroform:methanol (2:1, v/v). The total proteins content was determined according to the procedure described by the manufacturer (Pierce™ BCA protein assay kit) after an extraction with a phosphate buffered saline solution. The hydroxy-proline (OH-Pro) content was determined according to¹² on the delipidated and dry residue obtained after the total lipids extraction.

**Oxidative status.** GSH/GSSG analysis: The reduced to oxidized glutathione ratio (GSH/GSSG) was determined according to the protocol described by the manufacturer (catalog #: 239709, Abcam, Cambridge, UK).

The activities of the superoxide dismutase (SOD, catalog #: 19160, Sigma, St Louis, MO, USA) and catalase (Cat, catalog #: KB03012, BioQuoChem, Llanera-Asturias, Spain) enzymes were determined according to the procedures described by the manufacturer.

Results are expressed in U/mg of proteins.
<table><caption>**Table 1.** Chemical composition of the analyzed livers harvested before (B), at the beginning (BFF) or at the end (EFF) of the force-feeding period. RMSE root-mean-square error, OH-Pro hydroxyproline. Values are means ± SD. Means with the same superscripts (a, or b, or c) are not different (p < 0.05). *** = p < 0.001, * = p < 0.05.</caption><thead><tr><th></th><th>B</th><th>BFF</th><th>EFF</th><th>RMSE</th><th>p</th></tr></thead><tbody><tr><td>Liver weight (g)</td><td>150<sup>c</sup> ± 32</td><td>276<sup>b</sup> ± 43</td><td>632<sup>a</sup> ± 120</td><td>89.70</td><td>***</td></tr><tr><td>Dry matter (% raw tissue)</td><td>34.08<sup>c</sup> ± 1.23</td><td>45.29<sup>b</sup> ± 1.46</td><td>69.34<sup>a</sup> ± 3.59</td><td>4.77</td><td>***</td></tr><tr><td>Lipids (% of raw tissue)</td><td>4.92<sup>c</sup> ± 1.44</td><td>23.49<sup>b</sup> ± 1.85</td><td>52.91<sup>a</sup> ± 6.48</td><td>2.74</td><td>***</td></tr><tr><td>Proteins (% raw tissue)</td><td>12.7<sup>a</sup> ± 0.94</td><td>6.23<sup>b</sup> ± 1.21</td><td>3.09<sup>c</sup> ± 0.58</td><td>0.87</td><td>***</td></tr><tr><td>OH-Pro (mg / g delipidated and dried tissue)</td><td>0.16<sup>a</sup> ± 0.05</td><td>0.11<sup>b</sup> ± 0.03</td><td>0.15<sup>a</sup> ± 0.06</td><td>0.05</td><td>*</td></tr></tbody></table>

**ELISA tests.** HIF1α (Hypoxia Inducible Factor 1 Alpha), HIF2α (Hypoxia Inducible Factor 2 Alpha), GPX1 (Glutathione Peroxidase 1) and NQO1 (Nicotinamide Adenine Dinucleotide (Phosphate) Hydrogen Dehydrogenase (Quinone) 1) contents were determined with ELISA tests by using assay kits from MyBioSource (San Diego, CA, USA) according to the manufacturer's protocols.

Results are expressed in pg/mg of proteins (GPX1, HIF1α, HIF2α) or in ng/mg of proteins (NQO1).

**Caspases activities assays.** Activities of caspase-3, caspase-8 and caspase-9 enzymes were colorimetrically estimated via the CaspACE testing framework acquired from Biovision (Mountain View, CA, USA) in respect to the producer's guidelines. For each enzyme, the total activity of a pool of liver extracts composed with all the available samples was measured concomitantly with the individual samples. Then, the measured activity of this pool was arbitrary set to a value of 100 and the activity of each individual sample was expressed as a percentage of this value to finally indicate the enzyme activities in arbitrary units.

**Statistics.** Statistical analyses were performed with the SAS software, version 9.4 of the SAS System for Windows. Analysis of variance (ANOVA) were performed with the General Linear Model (Proc GLM) completed with the Student–Newman–Keuls' post-hoc test to compare the means of each groups. Values are expressed as mean ± standard deviation (SD). We set the significant level at p < 0.05.

**Ethics approval.** Ethical review and approval were not required for this animal study because we collected samples from a commercial slaughter house after the regular operations of anesthesia, bleeding, scalding, plucking and evisceration. We do not manipulate live animals and no specific experimental manipulations were performed during the time of rearing. All the samples we harvested were commercial ones, ready to be purchased by consumers.

## Results and discussion

As a whole, the rearing and the force-feeding programs took place without any noticeable events (mortality < 1.5%). As reported in the Table 1, after the 5 first meals of the force-feeding program (period B to BFF), the liver weights increased by 1.8 (average meal gain = 25 g) while it increased by 2.3 during the rest of the force-feeding time (period BFF to EFF, average meal gain = 24 g). This huge increase in liver weight was expected because the feeding program delivered very large quantities of corn which contain mainly carbohydrates (62% of starch). Those sugars are easily absorbed and brought to the liver where they are mainly metabolized in fatty acids by the de novo lipogenesis pathway and stored as triglycerides. This increase of the liver weight during force-feeding has already been largely documented in previous studies²,³,¹³.

During the force-feeding period, the dry matter and lipid contents dramatically increased while the protein one decreased. Those variations in the gross biochemical composition of the livers are conform to what were expected and simply trade the huge accumulation of lipids, and the relative dilution of proteins, during the development of hepatic steatosis¹,⁴,⁶,¹⁴ induced by the force-feeding. Considering the duration of the gavage period and the age of the animals between stages B and EFF, we assume that all the changes observed in the chemical composition of the liver are mainly due to the gavage and not to the increase in 1 age of the ducks (+10d for ducks 12w old). This is consistent with what has already been reported by²,³,⁶,¹⁴.

Even if the OH-Pro content of the livers was different from BFF to B and EFF groups (p < 0.05, Table 1), no big variation in the hepatic OH-Pro content was observed. This indicates that, despite its high intensity, the duration of the force-feeding is not sufficient to induce a significant increase in the development of the hepatic connective tissue. On the contrary, in mouse, during the shift from non-alcoholic fatty liver disease (NAFLD) to non-alcoholic steatohepatitis (NASH), an increase in OH-Pro content of the liver was reported¹⁵,¹⁶ and considered as a sign of the development of a fibrosis characteristic of the NASH state. To our knowledge, it is the first time that hydroxyproline content of the liver was determined during the development of the hepatic steatosis in force-fed ducks. Our results indicate that no NASH occurs during that period even if a large quantity of lipids accumulates in hepatocytes. This is probably due to the short time (10 d) of development of the hepatic steatosis in force-fed ducks which is not enough to generate a fibrosis per se.

The development of the hepatic steatosis is due to the imbalance between the capacities of synthesis and export of lipids by hepatocytes. This leads to a high accumulation of newly synthesized lipids despite the increase in the export of these same lipids by lipoproteins¹⁷,¹⁸. This is facilitated by the fact that corn contains few lipotropic
factors such as choline or methionine and a lot of vitamins B8 (biotin) that promotes the synthesis of fatty acids. According to¹⁹, one of the major issues with this imbalance in lipids metabolism is a subsequent metabolic stress due to the accumulation of free fatty acids (FFAs), which serve as a substrate to generate lipotoxic species that cause endoplasmic reticulum stress and produce reactive oxygen species (ROS). Those ROS are highly toxic for the cell and it tries to eliminate them with several mechanisms. Among them the coupling activity of superoxide dismutase (SOD) and catalase (Cat) is efficient to neutralise the superoxide ion $O_2^-$, one of the most powerful ROS. In the present study, SOD and Cat activities rapidly increase at the beginning (BFF) and all along the force-feeding period (Fig. 1). This indicates that force-feeding induced rapidly an oxidative stress in liver cells but the constancy of the SOD/Cat ratio also indicates that no major dysregulation of these coupled activities is present. In different pathological²⁰,²¹ or feeding²² or environmental conditions²³, the activities of SOD and Cat enzymes also increase to reduce the oxidative stress and are therefore considered to be one of the predominant antioxidant mechanisms in duck hepatocytes. Other enzymatic antioxidant systems (GPX1 and NQO1) activities are presented in Fig. 2. NQO1 removes quinone from biological systems through a detoxification reaction. This reaction oxidises the substrate without the formation of damaging semiquinone and oxygen free radical species²⁴. GPX1 acts in combination with reduced glutathione (GSH) to catalyze the reduction of $H_2O_2$ in $H_2O$ and oxidized glutathione (GSSG). In this study, both GPX1 and NQO1 increase significantly along the force-feeding period. This supports the idea that from the beginning until the end of the development of the steatosis, hepatocytes have to face an important oxidative stress and react by enhancing their antioxidant activities.

Glutathione is a small tripeptide which is present in all cells. Reduced glutathione (GSH) plays critical roles in protecting cells from oxidative damage and maintaining redox homeostasis²⁵. GSH is considered the most important intracellular antioxidant²⁶. It allows the reduction of ROS thanks to the action of glutathione peroxidases (GPx) which allows the transition from the reduced form (GSH) to the oxidized form (GSSG). The GSSG/GSH ratio is therefore a good indicator of good cellular health and a decrease in its value is a sign of major cellular dysfunction²⁷,²⁸. The cell keeps a permanent control on the GSSG/GSH ratio and in a state of oxidative stress, the consumption of GSH by GPx increasing, the cell activates several pathways to renew the quantity of GSH. The increased glutathione redox status is a compensatory mechanism in response to the increased oxidative burden due to the development of steatosis. In this study (Fig. 3), the level of the GSH content remained nearly constant while GSSG and GPX1 contents varied in opposite directions from the beginning to the end of the force-feeding period. During the development of the steatosis, the accumulation of lipids in hepatocytes induces an increase in the activity of the $\beta$-oxidation pathway and consequently a higher leak of electrons which link to the oxygen and generate ROS. If the level of response for producing more antioxidant molecules is not sufficient, excessive ROS induce strong dysregulations of mitochondrial functions²⁹. According to³⁰, in humans with severe NAFLD or NASH syndromes, the levels of hepatic SOD, catalase and GSH are low because the level of ROS is too high and exceeded the capacity of regulation of the hepatocytes. In the present study, these markers increase throughout the force-feeding period but ultimately do not fall which may indicate that hepatocytes are still able to cope with the rapid development of steatosis despite the observed increase in liver oxidative stress. The effective protective role of cellular antioxidants has also been suggested by⁴ who measured no increase in TBARS (thiobarbituric acid reactive substances) values in livers of different weights from force-fed ducks.

HIF1α and HIF2α are protein transcription factors involved in the response to initial or acute and delayed or chronic hypoxia, respectively³¹. In this study, HIF1α increased rapidly from start and throughout the force-feeding period while HIF2α only increased from BFF to EFF (Fig. 2). These increases indicate that during the period of force-feeding, the livers encountered hypoxic conditions that could be related to physical disability due to liver enlargement and circulatory problems linked to the elevated lipemia typically associated with the development of the hepatic steatosis. According to³², an increase in ROS can result from a reduction of an oxygen supply efficiency. According to³³, in mammals, impaired hepatic metabolism and tissue remodeling seen in fatty liver disease disrupt hepatic oxygen homeostasis, resulting in severe hepatic hypoxia. Reference³⁴ also reported that HIF-2 is overexpressed in humans or mice developing liver damage steatosis of various origins, which leads to a hypoxic microenvironment and mitochondrial dysfunctions. Then, it must be concluded that throughout the force-feeding period, livers are submitted to increased hypoxic situations that will favor the development of the observed cellular oxidative stress.

Caspase-3 is considered as a key effector enzyme in inducing cell apoptosis. Caspase-8 and caspase-9 are initiator caspases of extrinsic and intrinsic apoptosis, respectively³⁵. As reported in Fig. 4, caspase-8, -9 and -3 have a low level of activity in livers from B ducks. On the contrary, after 5 meals (BFF stage) or 20 meals (EFF stage) of force-feeding, activities of all the studied caspases significantly increase. This indicates that the development of the observed hepatic steatosis induced by the feeding program is associated with the development of apoptosis in the liver. This increase was previously reported in the liver of force-fed ducks⁴,⁵,³⁶. According to³⁷,³⁸, the production of reactive oxygen species, and the development of oxidative stress and endoplasmic reticulum stress are involved in the progression of apoptosis associated with the activation of NAFLD/NASH. Thus, in the present study, the hypoxic and oxidative environments of hepatocytes in force-fed ducks might also be involved in the observed onset of apoptosis development, while evidence for the transition from NAFLD to NASH is not still visible. Moreover, given that both activities of caspase-8 and -9 increased during the force-feeding period, none of the extrinsic or intrinsic activation pathways can be favored to explain the development of apoptosis.

## Conclusion

During the development of force-feeding induced hepatic steatosis in ducks, hepatocytes accumulate huge amounts of lipids. Reduced oxygen supply, which also occurs during the development of this hepatic steatosis, as well as lipotoxicity due to lipid accumulation, both contribute to the development of cellular oxidative stress
SOD (U /mg proteins)

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

Catalase (U / mg proteins)

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

SOD / Catalase

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

**Figure 1.** Activities of superoxide dismutase (SOD) and catalase enzymes in livers collected before (B) or at the beginning (BFF) or at the end (EFF) of the force-feeding period. Values are means ± SD. Bars with different letters differ (p < 0.05).
<figure><img src="image_8.png" /></figure>

GPX1 (pg / mg proteins)

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

NQO1 (ng / mg proteins)

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

HIF1α (pg / mg proteins)

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

HIF2α (pg / mg proteins)

**Figure 2.** Biochemical contents of glutathione peroxidase 1 (GPX1), nicotinamide adenine dinucleotide (phosphate) hydrogen dehydrogenase (quinone) 1 (NQO1), hypoxia inducible factor 1 alpha (HIF1α) and hypoxia inducible factor2 alpha (HIF2α) in livers collected before (B) or at the beginning (BFF) or at the end (EFF) of the force-feeding period. Values are means ± SD. Bars with different letters differ (p < 0.05).

conditions. Consequently, at the end of this period of force-feeding, many markers of oxidative stress are radi-cally modified but without any notable deregulation of their functionalities. This indicates that the hepatic cells are still actively coping with the metabolic challenge induced by the imposed feeding program. Nevertheless, the increased levels of activities of various caspases also indicate that these stressful conditions pushed an increased number of cells into the apoptotic cell death pathway.

In this experiment we demonstrate, for the first time, that oxidative stress sets in rapidly in the hepatocytes of ducks subjected to force-feeding. This could affect the normal functions of liver cells and it would be interest-ing to additionally test whether certain nutrients, such as natural antioxidants (vitamins C or E, beta-carotene, polyphenols, etc.) added in the diet before or during the force-feeding period, or synthetic products (Coenzyme Q10, glutathione, taurine, etc.) could attenuate this phenomenon (Supplementary Information).
GSSG (ng / mg of raw material)

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

GSH (ng / mg of raw material)

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

GSSG / GSH

<figure><img src="image_14.png" /><figcaption>**Figure 3.** Values of oxidized glutathione (GSSG) and reduced glutathione (GSH) contents in livers collected before (B) or at the beginning (BFF) or at the end (EFF) of the force-feeding period. Values are means ± SD. Bars with different letters differ (p < 0.05).</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>Caspase 9 (relative units)</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Caspase 3 (relative units)</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 4.** Caspase-3, caspase-8 and caspase-9 relative activities in livers collected before (B) or at the beginning (BFF) or at the end (EFF) of the force-feeding period. Values are means ± SD. Bars with different letters differ (*p* < 0.05).</figcaption></figure>
## Data availability

Details on the materials used are given in the M&M section. Raw data are available in the file data.xls for reviewers and on request to the corresponding author after publication.

Received: 4 November 2022; Accepted: 29 December 2022

Published online: 19 January 2023

## References

1. Herault, F. et al. Liver gene expression in relation to hepatic steatosis and lipid secretion in two duck species. *Anim. Genet.* **41**(1), 12–20. https://doi.org/10.1111/j.1365-2052.2009.01959.x (2010).

2. Bax, M. L. et al. Proteomic profile evolution during steatosis development in ducks. *Poult. Sci.* **91**, 112–120. https://doi.org/10.3382/ps.2011-01663 (2012).

3. Bonnefont, C. et al. Evolution of liver fattening and foie gras technological yield during the overfeeding period in mule duck. *Poult. Sci.* **98**(11), 5724–5733. https://doi.org/10.3382/ps/pez359 (2019).

4. Lo, B. et al. Cellular markers of mule duck livers after force-feeding. *Poult. Sci.* **99**(7), 3567–3573. https://doi.org/10.1016/j.psj.2020.03.048 (2020).

5. Lo, B. et al. Proteomic analysis of two weight classes of mule duck "foie gras" at the end of an overfeeding period. *Front. Physiol.* **11**, 569329. https://doi.org/10.3389/fphys.2020.569329 (2020).

6. Pioche, T. et al. Kinetic study of the expression of genes related to hepatic steatosis, glucose and lipid metabolism, and cellular stress during overfeeding in mule ducks. *Am. J. Physiol.-Regul. Integr. Comp. Physiol.* **318**(2), R453–R467. https://doi.org/10.1152/ajpregu.00198.2019 (2020).

7. Chen, X. et al. The effects of dietary Se on productive and reproductive performance, tibial quality, and antioxidant capacity in laying duck breeders. *Poult. Sci.* **99**(8), 3971–3978. https://doi.org/10.1016/j.psj.2020.04.006 (2020).

8. Friedman, S. L., Neuschwander-Tetri, B. A., Rinella, M. & Sanyal, A. J. Mechanisms of NAFLD development and therapeutic strategies. *Nat. Med.* **24**, 908–922. https://doi.org/10.1038/s41591-018-0104-9 (2018).

9. Masarone, M. et al. Role of oxidative stress in pathophysiology of nonalcoholic fatty liver disease. *Oxid. Med. Cell. Longev.* https://doi.org/10.1155/2018/9547613 (2018).

10. J.O.C.E. (1971). Dosage de l'humidité. *Journal Officiel des Communautés Européennes*. L 279/8.

11. Folch, J., Lee, M. & Sloane Stanley, G. H. A simple method for the isolation and purification of total lipides from animal tissues. *J. Biol. Chem.* **226**(1), 497–509 (1957).

12. Woessner, J. F. The determination of hydroxyproline in tissue and protein samples containing small proportions of this amino acid. *Arch. Biochem. Biophys.* **93**, 440–447. https://doi.org/10.1016/0003-9861(61)90291-0 (1961).

13. Gabarrou, J. F., Salichon, M. R., Guy, G. & Blum, J. C. Hybrid ducks overfed with boiled corn develop an acute hepatic steatosis with decreased choline and polyunsaturated fatty acid level in phospholipids. *Reprod. Nutr. Dev.* **36**, 473–484. https://doi.org/10.1051/rnd:19960503 (1996).

14. Tavernier, A. et al. Inter genotype differences in expression of genes involved in glucose metabolism in the establishment of hepatic steatosis in Muscovy, Pekin and mule ducks. *Mol. Biol. Rep.* **47**(2), 1527–1533. https://doi.org/10.1007/s11033-019-05182-y (2020).

15. Matsumoto, M. et al. An improved mouse model that rapidly develops fibrosis in non-alcoholic steatohepatitis. *Int. J. Exp. Pathol.* **94**(2), 93–103. https://doi.org/10.1111/iep.12008 (2013).

16. Hartimath, S. V. et al. Imaging fibrogenesis in a diet-induced model of nonalcoholic steatohepatitis (NASH). *Contrast Media Mol. Imaging* https://doi.org/10.1155/2019/6298128 (2019).

17. Hermier, D. et al. Differential channelling of liver lipids in relation to susceptibility to hepatic steatosis in two species of ducks. *Comp. Biochem. Physiol. B* **135**, 663–675. https://doi.org/10.1016/s1096-4959(03)00146-5 (2003).

18. Tavernier, A., Ricaud, K., Bernardet, M.-D., Gontier, K. & Davail, S. Pre- and post-prandial expression of genes involved in lipid metabolism at the end of the overfeeding period of mule ducks. *Mol. Cell. Biochem.* **438**(1–2), 111–121. https://doi.org/10.1007/s11010-017-3118-6 (2018).

19. Chen, Z., Yu, Y., Cai, J. & Li, H. Emerging molecular targets for treatment of nonalcoholic fatty liver disease. *Trends Endocrinol. Metab.* **30**(12), 903–914. https://doi.org/10.1016/j tem.2019.08.006 (2019).

20. Du, H. et al. A flavone-polysaccharide based prescription attenuates the mitochondrial dysfunction induced by duck hepatitis A virus type 1. *PLoS ONE* **12**(4), e0175495. https://doi.org/10.1371/journal.pone.0175495 (2017).

21. Zhang, L. et al. Ameliorating effects of *Bacillus subtilis* ANSB060 on growth performance, antioxidant functions, and aflatoxin residues in ducks fed diets contaminated with aflatoxins. *Toxins* **9**, 1. https://doi.org/10.3390/toxins9010001 (2017).

22. Liu, Y. et al. Effects of tryptophan supplementation on growth performance, antioxidative activity, and meat quality of ducks under high stocking density. *Poult. Sci.* **94**(8), 1894–1901. https://doi.org/10.3382/ps/pev155 (2015).

23. Zeng, T. et al. Effects of heat stress on antioxidant defense system, inflammatory injury, and heat shock proteins of Muscovy and Pekin ducks: Evidence for differential thermal sensitivities. *Cell Stress Chaperones* **19**(6), 895–901. https://doi.org/10.1007/s12192-014-0514-7 (2014).

24. Ross, D. & Siegel, D. Functions of NQO1 in cellular protection and CoQ10 metabolism and its potential role as a redox sensitive molecular switch. *Front. Physiol.* https://doi.org/10.3389/fphys.2017.00595 (2017).

25. Forman, H. J., Zhang, H. & Rinna, A. Glutathione: Overview of its protective roles, measurement and biosynthesis. *Mol. Aspects Med.* **30**(1–2), 1–12. https://doi.org/10.1016/j.mam.2008.08.006 (2009).

26. Sahoo, S., Awasthi, J. P., Sunkar, R. & Panda, S. K. Determining glutathione levels in plants. In *Plant Stress Tolerance* (ed. Sunkar, R.) 273–277 (Springer, New York, 2017). https://doi.org/10.1007/978-1-4939-7136-7_16.

27. Chen, Y. et al. Glutathione defense mechanism in liver injury: Insights from animal models. *Food Chem. Toxicol.* **60**, 38–44. https://doi.org/10.1016/j.fct.2013.07.008 (2013).

28. Alkazemi, D., Rahman, A. & Habra, B. Alterations in glutathione redox homeostasis among adolescents with obesity and anemia. *Sci. Rep.* **11**, 3034. https://doi.org/10.1038/s41598-021-82579-5 (2021).

29. Guo, C., Sun, L., Chen, X. & Zhang, D. Oxidative stress, mitochondrial damage and neurodegenerative diseases. *Neural Regen. Res.* **8**(21), 2003–2014. https://doi.org/10.3969/j.issn.1673-5374.2013.21.009 (2013).

30. Sumida, Y., Niki, E., Naito, Y. & Yoshikawa, T. Involvement of free radicals and oxidative stress in NAFLD/NASH. *Free Radical Res.* **47**, 869–880. https://doi.org/10.3109/10715762.2013.837577 (2013).

31. Downes, N. L., Laham-Karam, N., Kaikkonen, M. U. & Ylä-Herttuala, S. Differential but complementary HIF1α and HIF2α transcriptional regulation. *Mol. Ther.* **26**(7), 1735–1745. https://doi.org/10.1016/j.ymthe.2018.05.004 (2018).

32. McGarry, T., Biniecka, M., Veale, D. J. & Fearon, U. Hypoxia, oxidative stress and inflammation. *Free Radical Biol. Med.* **125**, 15–24. https://doi.org/10.1016/j.freeradbiomed.2018.03.042 (2018).

33. Suzuki, T., Shinjo, S., Arai, T., Kanai, M. & Goda, N. Hypoxia and fatty liver. *World J. Gastroenterol.* **20**, 15087–15097. https://doi.org/10.3748/wjg.v20.i41.15087 (2014).

34. Morello, E. et al. Hypoxia-inducible factor 2α drives nonalcoholic fatty liver progression by triggering hepatocyte release of histidine-rich glycoprotein. *Hepatology* **67**, 2196–2214. https://doi.org/10.1002/hep.29754 (2018).
35. McIlwain, D. R., Berger, T. & Mak, T. W. Caspase functions in cell death and disease. *Cold Spring Harbor Perspect. Biol.* 5(4), a008656. https://doi.org/10.1101/cshperspect.a008656 (2013).

36. Remignon, H., Yahia, R., Marty-Gasset, N. & Wilkesman, J. Apoptosis during the development of the hepatic steatosis in force-fed ducks and cooking yield implications. *Poult. Sci.* **97**(6), 2211–2217. https://doi.org/10.3382/ps/pey054 (2018).

37. Majmundar, A. J., Wong, W. J. & Simon, M. C. Hypoxia-inducible factors and the response to hypoxic stress. *Mol. Cell* **40**(2), 294–309. https://doi.org/10.1016/j.molcel.2010.09.022 (2010).

38. Kanda, T. et al. Apoptosis and non-alcoholic fatty liver diseases. *World J. Gastroenterol.* **24**(25), 2661–2672. https://doi.org/10.3748/wjg.v24.i25.2661 (2018).

**Author contributions**

H.R. contributed to the conception, the design of the work and the acquisition, analysis and interpretation of data for the work. P.B. contributed to the acquisition and analysis of data for the work. Both authors approved that version of manuscript submitted for publication.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-27271-y.

**Correspondence** and requests for materials should be addressed to H.R.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

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

© The Author(s) 2023