# Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?

Stéphane Beauclercq, Antoine Lefèvre, Lydie Nadal-Desbarats, Karine Germain, Christophe Praud, Patrick Emond, Elisabeth Le Bihan-Duval, Sandrine Mignon-Grasteau

► To cite this version:

Stéphane Beauclercq, Antoine Lefèvre, Lydie Nadal-Desbarats, Karine Germain, Christophe Praud, et al.. Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?. Poultry Science, 2019, 98 (3), pp.1425-1431. 10.3382/ps/pey483. hal-02619205

HAL Id: hal-02619205

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

Submitted on 18 Jan 2024

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

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

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

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

HAL
open science

# Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?

Stéphane Beauclercq, Antoine Lefèvre, Lydie Nadal-Desbarats, Karine

Germain, Christophe Praud, Patrick Emond, Elisabeth Le Bihan-Duval,

Sandrine Mignon-Grasteau

► To cite this version:

Stéphane Beauclercq, Antoine Lefèvre, Lydie Nadal-Desbarats, Karine Germain, Christophe Praud, et al.. Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?. Poultry Science, 2019, 98 (3), pp.1425-1431. 10.3382/ps/pey483. hal-02619205

HAL Id: hal-02619205

https://hal.inrae.fr/hal-02619205v1

Submitted on 18 Jan 2024

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

Distributed under a Creative Commons Attribution 4.0 International License
# Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?

Stéphane Beauclercq ID,* Antoine Lefèvre,† Lydie Nadal-Desbarats,†,‡ Karine Germain,§ Christophe Praud,* Patrick Emond,†,‡,# Elisabeth Le Bihan-Duval,* and Sandrine Mignon-Grasteau*,**

*BOA, INRA, Université de Tours, 37380 Nouzilly, France; †Université de Tours, PST Analyse des systèmes biologiques, 37044 Tours Cedex 9, Tours, France; ‡UMR 1253, iBrain, Université de Tours, Inserm, 37032 Tours cedex 1, Tours, France; §INRA, UE 1206 Elevage Alternatif et Santé des Monogastriques, Le Magneraud, Saint-Pierre-d'Amilly, BP 52, 17700 Surgères, France; and #CHRU de Tours, Service de Médecine Nucléaire In Vitro, 37032 Tours cedex, Tours, France

**ABSTRACT** The increasing cost of conventional feedstuffs used in poultry diets has bolstered interest in genetic selection for digestive efficiency (DE) to improve the adaptation of the birds to various alternative feedstuffs. However, DE measurement through AMEn is time-consuming and constraining. To simplify selection for DE, the potential of serum composition to predict AMEn was evaluated based on 40 birds from two broiler lines (D+ and D−) divergently selected on the fecal AMEn of a difficult-to-digest wheat-based diet. Differences in serum coloration were suspected between the two lines, and thus a spectrophotometric analysis was carried out, revealing a signif-

icant difference in absorption between 430 nm and 516 nm, corresponding to the signature of orange-red lipophilic pigments such as xanthophylls. To go further, the liposoluble fraction of the serum was explored for its lipidome by mass spectrometry. Discriminant analysis revealed that a pattern of 10 metabolites, including zeaxanthin/lutein, can explain 82% of the lipidomic differences between the two lines. Corimetry combined with lipidomics studies confirmed the relationship between digestive efficiency and serum composition, which opens up new possibilities for using it as a quick and easy proxy of digestive efficiency.

**Key words:** digestive efficiency, lipidomics, serum, biomarker, carotenoid

2019 Poultry Science 98:1425–1431

http://dx.doi.org/10.3382/ps/pey483

## INTRODUCTION

The intense selection performed during the past decades for chicken production-related traits has led to huge genetic progress when considering feed conversion ratio and growth rate, with broiler chickens being slaughtered 19% earlier, and requiring 24% less feed to reach the same weight in 2010 than in 1970 (Riffard et al., 2011). This is a success if considering solely the purpose of increasing animal protein production at a lower cost, which has been the goal of genetic selection in livestock until now. Nevertheless, this advance relies on the use of optimal feedstuffs in poultry diets, which are becoming more and more expensive and are

increasingly in competition with human food. Therefore, digestive efficiency (DE), a component of feed conversion ratio, must also be improved in order to include less digestible feedstuffs in the diets. DE of energy (measured by AMEn), proteins, lipids, and starch have a genetic origin and have been successfully selected in a divergent selection experiment in chickens (Mignon-Grasteau et al., 2004). After eight generations of selection, the high DE (D+) and the low DE (D−) lines showed a difference of 40% in AMEn with no impact on the body weight (de Verdal et al., 2011). However, the only available way to assess DE is through balance trials, which are time-consuming and detrimental to animal welfare as animals are kept in individual cages, which consequently limits its integration in genetic selection schemes. Using indirect, non-invasive indicators of DE such as metabolomic biomarkers would thus be useful. However, currently none are available and no extensive metabolomics profiling of DE has been performed.

In a prior study, we showed that the metabolomics profile of the hydrosoluble fraction of serum was linked to DE (Beauclercq et al., 2018). In this study, we observed that the color of serum seemed to differ

© The Author(s) 2018. Published by Oxford University Press on behalf of Poultry Science Association. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals permissions@oup.com.

Received February 21, 2018.

Accepted September 19, 2018.

<footnote n="1">¹Corresponding author: E-mail: sandrine.grasteau@inra.fr</footnote>
# Does lipidomic serum analysis support the assessment of digestive efficiency in chickens?

Stéphane Beauclercq ID,* Antoine Lefèvre,† Lydie Nadal-Desbarats,†,‡ Karine Germain,§ Christophe Praud,* Patrick Emond,†,‡,# Elisabeth Le Bihan-Duval,* and Sandrine Mignon-Grasteau*,¹

*BOA, INRA, Université de Tours, 37380 Nouzilly, France; †Université de Tours, PST Analyse des systèmes biologiques, 37044 Tours Cedex 9, Tours, France; ‡UMR 1253, iBrain, Université de Tours, Inserm, 37032 Tours cedex 1, Tours, France; §INRA, UE 1206 Elevage Alternatif et Santé des Monogastriques, Le Magneraud, Saint-Pierre-d'Amilly, BP 52, 17700 Surgères, France; and #CHRU de Tours, Service de Médecine Nucléaire In Vitro, 37032 Tours cedex, Tours, France

**ABSTRACT** The increasing cost of conventional feedstuffs used in poultry diets has bolstered interest in genetic selection for digestive efficiency (DE) to improve the adaptation of the birds to various alternative feedstuffs. However, DE measurement through AMEn is time-consuming and constraining. To simplify selection for DE, the potential of serum composition to predict AMEn was evaluated based on 40 birds from two broiler lines (D+ and D−) divergently selected on the fecal AMEn of a difficult-to-digest wheat-based diet. Differences in serum coloration were suspected between the two lines, and thus a spectrophotometric analysis was carried out, revealing a signif-

icant difference in absorption between 430 nm and 516 nm, corresponding to the signature of orange-red lipophilic pigments such as xanthophylls. To go further, the liposoluble fraction of the serum was explored for its lipidome by mass spectrometry. Discriminant analysis revealed that a pattern of 10 metabolites, including zeaxanthin/lutein, can explain 82% of the lipidomic differences between the two lines. Corimetry combined with lipidomics studies confirmed the relationship between digestive efficiency and serum composition, which opens up new possibilities for using it as a quick and easy proxy of digestive efficiency.

**Key words:** digestive efficiency, lipidomics, serum, biomarker, carotenoid

2019 Poultry Science 98:1425–1431

http://dx.doi.org/10.3382/ps/pey483

## INTRODUCTION

The intense selection performed during the past decades for chicken production-related traits has led to huge genetic progress when considering feed conversion ratio and growth rate, with broiler chickens being slaughtered 19% earlier, and requiring 24% less feed to reach the same weight in 2010 than in 1970 (Riffard et al., 2011). This is a success if considering solely the purpose of increasing animal protein production at a lower cost, which has been the goal of genetic selection in livestock until now. Nevertheless, this advance relies on the use of optimal feedstuffs in poultry diets, which are becoming more and more expensive and are

increasingly in competition with human food. Therefore, digestive efficiency (DE), a component of feed conversion ratio, must also be improved in order to include less digestible feedstuffs in the diets. DE of energy (measured by AMEn), proteins, lipids, and starch have a genetic origin and have been successfully selected in a divergent selection experiment in chickens (Mignon-Grasteau et al., 2004). After eight generations of selection, the high DE (D+) and the low DE (D−) lines showed a difference of 40% in AMEn with no impact on the body weight (de Verdal et al., 2011). However, the only available way to assess DE is through balance trials, which are time-consuming and detrimental to animal welfare as animals are kept in individual cages, which consequently limits its integration in genetic selection schemes. Using indirect, non-invasive indicators of DE such as metabolomic biomarkers would thus be useful. However, currently none are available and no extensive metabolomics profiling of DE has been performed.

In a prior study, we showed that the metabolomics profile of the hydrosoluble fraction of serum was linked to DE (Beauclercq et al., 2018). In this study, we observed that the color of serum seemed to differ

© The Author(s) 2018. Published by Oxford University Press on behalf of Poultry Science Association. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals permissions@oup.com.

Received February 21, 2018.

Accepted September 19, 2018.

¹Corresponding author: E-mail: sandrine.grasteau@inra.fr
between efficient and non-efficient birds, thus suggesting that pigments also differed between the two lines. As

pigments are located in the lipophilic fraction of serum, we proposed to screen the liposoluble fraction of D+ and D− broiler serum to search for potential biomarkers of DE within this fraction. Liquid chromatography-high-resolution mass spectrometry (LC-HRMS) was used to characterize the lipidomic composition of the lipophilic part of the serum.

## MATERIALS AND METHODS

All animal care and experimental procedures needed for this study were approved by the Ethics Committee for Animal Experimentation of Val de Loire (Authorization No. 01047.02). This ethics committee is registered by the National Committee under the number C2EA 19.

All chemicals were bought from Sigma Aldrich (Saint-Quentin Fallavier, France) unless otherwise specified.

### Birds and Sampling

This study was conducted on birds originating from the 20th generation of two genetic lines divergently selected for low (D−) or high (D+) DE through faecal AMEn of a difficult to digest wheat-based diet at 3 wk (Mignon-Grasteau et al., 2004; Carré et al., 2008). The animals were reared at UE1206 EASM (Unité Expérimentale Elevage Alternatif et Santé des Mono-gastriques, INRA, Le Magneraud, France), on the floor during the first 2 wk and then in individual cages to measure AMEn. During the whole experiment, the animals were fed ad libitum with the same diet as that used during the initial selection experiment, i.e., a diet adapted to nutritional requirements (2,943 kcal.kg⁻¹ dry matter, 21% crude protein; 6% lipids), but of low digestibility due to its high content (55%) of Rialto wheat, a variety with high viscosity (Carré et al., 2002).

Birds were weighed at 17, 23, 25, and 26 d and their feed intake was recorded individually between 17 and 23 d, 23 and 25 d, and 25 and 26 d. Their feed conversion ratio between 17 and 26 d was calculated as the ratio of feed intake to growth during this period.

The balance trial was done between 23 and 25 d to get fecal AMEn and coefficients of digestive use of starch (CDU~S~), nitrogen (CDU~N~) and lipids (CDU~L~) from 20 birds of each line (males and females) using a method based on total excreta collection, as described by Bourdillon et al. (1990). Gross energy, lipid, starch, and nitrogen content of individual freeze-dried excreta were measured for all birds using the Near Infrared Spectroscopy procedure (Foss, Hilleroed, Denmark) described by Bastianelli et al. (2010).

After the end of the balance trial, blood was drawn at the occipital sinus for serum composition analyses. Serum was prepared by keeping the blood at room tem-

perature for 15 min until coagulation and centrifugation (3,000 g for 10 min). Sera were aliquoted and stored at -20°C until further analysis.

The equality of medians between the D− and D+ (n = 20) samples was tested for the DE traits with the non-parametric Mann-Whitney's test.

### Serum Color

Samples of 200 µL from the serum of the 40 birds were transferred to a transparent 96-well plate (Greiner Bio-One, Kremsmünster, Austria) and their absorption spectra were acquired between 400 and 600 nm (2 nm steps) using an Infinite M200 spectrophotometer (Tecan, Männedorf, Switzerland). The equality of means between the D− (n = 20) and D+ (n = 20) serum absorption was tested for each 2 nm step at the 95% confidence level with Welch's t-test.

### LC-HRMS Analysis

The liposoluble fraction of the serum of each bird was extracted based on the Bligh and Dyer method (Bligh and Dyer, 1959). Briefly, 375 µL of chloroform/methanol (1:2) were added to 100 µL of serum, the mix was then thoroughly vortexed before addition of 125 µL chloroform. A volume of 125 µL of water was subsequently added followed by mixing and centrifugation (15,000 g, 10 min, 4°C). The lower phase (200 µL) under the protein disk corresponding to the nonpolar fraction was recovered and put in glass tubes for further solvent evaporation in a SpeedVac (Thermo Fisher Scientific, Waltham, MA) at room temperature. The residue was then reconstituted with 100 µL of a 6:3:1 mix of acetonitrile (ACN)/water/isopropanol followed by centrifugation (15,000 g, 10 min, 4°C) before mass spectrometry analysis.

LC-HRMS analysis was performed on a UHPLC Ultimate 3000 system (Dionex, Sunnyvale, CA), coupled to a Q-Exactive mass spectrometer (Thermo Fisher Scientific) and operated in positive ionization mode (ESI+). Chromatography was carried out with a 1.7 µm XB-C18 (150 mm × 2.10 mm, 100 Å) UHPLC column (Kinetex, Phenomenex, Torrance, CA) heated at 55°C. The solvent system comprised mobile phase A [isopropanol/ACN (9:1) + 0.1% (vol/vol) formic acid + 10 mM ammonium formate], and mobile phase B [ACN/water (6:1) + 0.1% (vol/vol) formic acid + 10 mM ammonium formate]; the gradient operated at a flow rate of 0.26 mL/min over a run time of 24 min. The multisteps gradient was programmed as follows: 0 to 1.5 min, 32% to 45% A; 1.5 to 5 min, 45% to 52% A; 5 to 8 min, 52% to 58% A; 8 to 11 min, 58% to 66% A; 11 to 14 min, 66% to 70% A; 14 to 18 min, 70% to 75% A; 18 to 21 min, 75% to 97% A; 21 to 24 min, 97% A. The autosampler temperature (Ultimate WPS-3000 UHPLC system, Dionex) was set at 4°C, and the injection volume for each sample was 5 µL. The heated
between efficient and non-efficient birds, thus suggesting that pigments also differed between the two lines. As

pigments are located in the lipophilic fraction of serum, we proposed to screen the liposoluble fraction of D+ and D− broiler serum to search for potential biomarkers of DE within this fraction. Liquid chromatography-high-resolution mass spectrometry (LC-HRMS) was used to characterize the lipidomic composition of the lipophilic part of the serum.

## MATERIALS AND METHODS

All animal care and experimental procedures needed for this study were approved by the Ethics Committee for Animal Experimentation of Val de Loire (Authorization No. 01047.02). This ethics committee is registered by the National Committee under the number C2EA 19.

All chemicals were bought from Sigma Aldrich (Saint-Quentin Fallavier, France) unless otherwise specified.

### Birds and Sampling

This study was conducted on birds originating from the 20th generation of two genetic lines divergently selected for low (D−) or high (D+) DE through faecal AMEn of a difficult to digest wheat-based diet at 3 wk (Mignon-Grasteau et al., 2004; Carré et al., 2008). The animals were reared at UE1206 EASM (Unité Expérimentale Elevage Alternatif et Santé des Mono-gastriques, INRA, Le Magneraud, France), on the floor during the first 2 wk and then in individual cages to measure AMEn. During the whole experiment, the animals were fed ad libitum with the same diet as that used during the initial selection experiment, i.e., a diet adapted to nutritional requirements (2,943 kcal.kg⁻¹ dry matter, 21% crude protein; 6% lipids), but of low digestibility due to its high content (55%) of Rialto wheat, a variety with high viscosity (Carré et al., 2002).

Birds were weighed at 17, 23, 25, and 26 d and their feed intake was recorded individually between 17 and 23 d, 23 and 25 d, and 25 and 26 d. Their feed conversion ratio between 17 and 26 d was calculated as the ratio of feed intake to growth during this period.

The balance trial was done between 23 and 25 d to get fecal AMEn and coefficients of digestive use of starch (CDU~S~), nitrogen (CDU~N~) and lipids (CDU~L~) from 20 birds of each line (males and females) using a method based on total excreta collection, as described by Bourdillon et al. (1990). Gross energy, lipid, starch, and nitrogen content of individual freeze-dried excreta were measured for all birds using the Near Infrared Spectroscopy procedure (Foss, Hilleroed, Denmark) described by Bastianelli et al. (2010).

After the end of the balance trial, blood was drawn at the occipital sinus for serum composition analyses. Serum was prepared by keeping the blood at room tem-

perature for 15 min until coagulation and centrifugation (3,000 g for 10 min). Sera were aliquoted and stored at −20°C until further analysis.

The equality of medians between the D− and D+ (n = 20) samples was tested for the DE traits with the non-parametric Mann-Whitney's test.

### Serum Color

Samples of 200 µL from the serum of the 40 birds were transferred to a transparent 96-well plate (Greiner Bio-One, Kremsmünster, Austria) and their absorption spectra were acquired between 400 and 600 nm (2 nm steps) using an Infinite M200 spectrophotometer (Tecan, Männedorf, Switzerland). The equality of means between the D− (n = 20) and D+ (n = 20) serum absorption was tested for each 2 nm step at the 95% confidence level with Welch's t-test.

### LC-HRMS Analysis

The liposoluble fraction of the serum of each bird was extracted based on the Bligh and Dyer method (Bligh and Dyer, 1959). Briefly, 375 µL of chloroform/methanol (1:2) were added to 100 µL of serum, the mix was then thoroughly vortexed before addition of 125 µL chloroform. A volume of 125 µL of water was subsequently added followed by mixing and centrifugation (15,000 g, 10 min, 4°C). The lower phase (200 µL) under the protein disk corresponding to the nonpolar fraction was recovered and put in glass tubes for further solvent evaporation in a SpeedVac (Thermo Fisher Scientific, Waltham, MA) at room temperature. The residue was then reconstituted with 100 µL of a 6:3:1 mix of acetonitrile (ACN)/water/isopropanol followed by centrifugation (15,000 g, 10 min, 4°C) before mass spectrometry analysis.

LC-HRMS analysis was performed on a UHPLC Ultimate 3000 system (Dionex, Sunnyvale, CA), coupled to a Q-Exactive mass spectrometer (Thermo Fisher Scientific) and operated in positive ionization mode (ESI+). Chromatography was carried out with a 1.7 µm XB—C18 (150 mm × 2.10 mm, 100 Å) UHPLC column (Kinetex, Phenomenex, Torrance, CA) heated at 55°C. The solvent system comprised mobile phase A [isopropanol/ACN (9:1) + 0.1% (vol/vol) formic acid + 10 mM ammonium formate], and mobile phase B [ACN/water (6:1) + 0.1% (vol/vol) formic acid + 10 mM ammonium formate]; the gradient operated at a flow rate of 0.26 mL/min over a run time of 24 min. The multisteps gradient was programmed as follows: 0 to 1.5 min, 32% to 45% A; 1.5 to 5 min, 45% to 52% A; 5 to 8 min, 52% to 58% A; 8 to 11 min, 58% to 66% A; 11 to 14 min, 66% to 70% A; 14 to 18 min, 70% to 75% A; 18 to 21 min, 75% to 97% A; 21 to 24 min, 97% A. The autosampler temperature (Ultimate WPS-3000 UHPLC system, Dionex) was set at 4°C, and the injection volume for each sample was 5 µL. The heated
ESI source parameters were a spray voltage of 3.5 kV, capillary temperature of 350°C, heater temperature of 250°C, sheath gas flow of 35 arbitrary units (AU), auxiliary gas flow of 10 AU, spare gas flow of 1 AU, and tube lens voltage of 60 V for C18. During the full-scan acquisition, which ranged from 250 to 1600 m/z, the instrument operated at 70,000 resolution, with an automatic gain control target of 1 × 10⁶ charges and a maximum injection time of 250 ms.

The instrumental stability was evaluated by multiple injections (n = 5) of a quality control (QC) sample obtained from a pool of 10 µL of all samples analyzed. This QC sample was injected once at the beginning of the analysis, between every 10 sample injections, and at the end of the run.

## LC-HRMS Data Processing

XCMS online (https://xcmsonline.scripps.edu) was used to process the raw data for feature detection, retention time correction, and alignment using the default parameter for UHPLC Q-Extractive (Tautenhahn et al., 2012). Features with variance intensities greater than 30% in QC samples or with QC variance superior to sample variance were removed as well as those identified as background noise or poorly integrated after visual inspection (Want et al., 2012). The peak intensities were normalized by the median area of the spectra for each sample, transformed into a logarithm using the automatic transformation tool of SIMCA 13 (Umetrics, Umeå, Sweden), and centered to maximize the separation between the D+ and D− lines. The putative assignment of the peaks was based on a METLIN search of the accurate mass (Smith et al., 2005). Lipids were annotated according the LIPID MAPS Lipid Classification System (Fahy et al., 2005, 2009).

## Chemometric Analyses

Principal Component Analysis (PCA) was performed on the data set as exploratory unsupervised analysis. Individuals outside of the 95% confidence interval of the Hotelling's T-square were considered to be outliers and excluded from the subsequent analyses.

Partial least squares discriminant analysis (PLS-DA), a method of supervised classification, was performed to predict the chicken group (i.e., D+ or D−) by explanatory quantitative variables (defined here by the peaks), using the SIMCA 13 Software (Barker and Rayens, 2003). The selection of metabolomics variables was carried out by iteratively excluding the variables with low regression coefficients and wide confidence intervals derived from jackknifing combined with low variable importance in the projection (VIP) until maximal improvement of the quality of the models. The model quality was evaluated after 7-fold cross validation by R²Y (goodness of fit), Q² (goodness of prediction), and CV-ANOVA (cross validation-analysis of

variance). CV-ANOVA is a diagnostic tool for assessing the reliability of PLS-DA models; the returned *P* value is indicative of the statistical significance of the fitted model (Eriksson et al., 2008). The contribution of each predictor in the model was evaluated through the variable contribution plot.

## RESULTS AND DISCUSSION

### Bird Digestive Efficiency Traits

The composition of the liposoluble fraction of serum was investigated in 20 individuals from the D+ line and 20 from the D- line. Their average growth, feed conversion ratio, and DE values are described in Table 1. The individuals from the two lines have a similar body weight at 17 d. However, the feed conversion ratio of D− was 16% higher due to a lower DE (−16% for the CDU<sub>L</sub>, −7% for CDU<sub>N</sub>, and −12% in AMEn).

### Serum Pigmentation and Digestive Efficiency

The coloring of the sera was noticeably different between the D+ and D− lines, which was confirmed by spectrophotometry. Indeed, the absorption spectra between 430 and 516 nm, corresponding to a yellow-orange color, were different (*P* < 0.05) between the two lines. The difference was maximal at 492 nm with a serum absorption 31% lower in D− than in D+ birds. The significance of differences between the lines for the absorption profile of the pigment molecules are presented in Figure 1a. This profile was characterized by two peaks at around 466 nm and 492 nm, which are close to the xanthophyll and carotenoid absorption spectra (Bernstein et al., 2010). It is worth noticing that the absorption spectra of those molecules are influenced by the solvent (Zang et al., 1997), as well as by the proportion of the different xanthophylls and carotenoids in the mix, which could explain the slight difference between our spectrum and the Bernstein's et al. spectra acquired in olive oil.

These differences in pigmentation may be an ideal diagnostic tool to discriminate between chickens with high or low DE because of its ease of implementation. It has previously been used in chickens to determine their status regarding coccidiosis infection, a parasitic disease of the intestine (Hamzic et al., 2015). In cases of infection, the intestinal absorption of pigments is disturbed and pigments are partially lost by leakage through the damaged wall of the cecum, leading to a reduced level of xanthophylls and carotenoids in the plasma from infected broilers (Ruff et al., 1974; Ruff and Fuller, 1975). The mechanism leading to the difference in serum pigmentation between D+ and D− could be similar as absorption capacity is probably lower in D− than in D+, including for pigments. Indeed, de Verdal et al. (2010) showed that the structure of the intestine was more
ESI source parameters were a spray voltage of 3.5 kV, capillary temperature of 350°C, heater temperature of 250°C, sheath gas flow of 35 arbitrary units (AU), auxiliary gas flow of 10 AU, spare gas flow of 1 AU, and tube lens voltage of 60 V for C18. During the full-scan acquisition, which ranged from 250 to 1600 m/z, the instrument operated at 70,000 resolution, with an automatic gain control target of 1 × 10⁶ charges and a maximum injection time of 250 ms.

The instrumental stability was evaluated by multiple injections (n = 5) of a quality control (QC) sample obtained from a pool of 10 µL of all samples analyzed. This QC sample was injected once at the beginning of the analysis, between every 10 sample injections, and at the end of the run.

## LC-HRMS Data Processing

XCMS online (https://xcmsonline.scripps.edu) was used to process the raw data for feature detection, retention time correction, and alignment using the default parameter for UHPLC Q-Extractive (Tautenhahn et al., 2012). Features with variance intensities greater than 30% in QC samples or with QC variance superior to sample variance were removed as well as those identified as background noise or poorly integrated after visual inspection (Want et al., 2012). The peak intensities were normalized by the median area of the spectra for each sample, transformed into a logarithm using the automatic transformation tool of SIMCA 13 (Umetrics, Umeå, Sweden), and centered to maximize the separation between the D+ and D− lines. The putative assignment of the peaks was based on a METLIN search of the accurate mass (Smith et al., 2005). Lipids were annotated according the LIPID MAPS Lipid Classification System (Fahy et al., 2005, 2009).

## Chemometric Analyses

Principal Component Analysis (PCA) was performed on the data set as exploratory unsupervised analysis. Individuals outside of the 95% confidence interval of the Hotelling's T-square were considered to be outliers and excluded from the subsequent analyses.

Partial least squares discriminant analysis (PLS-DA), a method of supervised classification, was performed to predict the chicken group (i.e., D+ or D−) by explanatory quantitative variables (defined here by the peaks), using the SIMCA 13 Software (Barker and Rayens, 2003). The selection of metabolomics variables was carried out by iteratively excluding the variables with low regression coefficients and wide confidence intervals derived from jackknifing combined with low variable importance in the projection (VIP) until maximal improvement of the quality of the models. The model quality was evaluated after 7-fold cross validation by R²Y (goodness of fit), Q² (goodness of prediction), and CV-ANOVA (cross validation-analysis of

variance). CV-ANOVA is a diagnostic tool for assessing the reliability of PLS-DA models; the returned *P* value is indicative of the statistical significance of the fitted model (Eriksson et al., 2008). The contribution of each predictor in the model was evaluated through the variable contribution plot.

## RESULTS AND DISCUSSION

### Bird Digestive Efficiency Traits

The composition of the liposoluble fraction of serum was investigated in 20 individuals from the D+ line and 20 from the D- line. Their average growth, feed conversion ratio, and DE values are described in Table 1. The individuals from the two lines have a similar body weight at 17 d. However, the feed conversion ratio of D− was 16% higher due to a lower DE (−16% for the CDU<sub>L</sub>, −7% for CDU<sub>N</sub>, and −12% in AMEn).

### Serum Pigmentation and Digestive Efficiency

The coloring of the sera was noticeably different between the D+ and D− lines, which was confirmed by spectrophotometry. Indeed, the absorption spectra between 430 and 516 nm, corresponding to a yellow-orange color, were different (*P* < 0.05) between the two lines. The difference was maximal at 492 nm with a serum absorption 31% lower in D− than in D+ birds. The significance of differences between the lines for the absorption profile of the pigment molecules are presented in Figure 1a. This profile was characterized by two peaks at around 466 nm and 492 nm, which are close to the xanthophyll and carotenoid absorption spectra (Bernstein et al., 2010). It is worth noticing that the absorption spectra of those molecules are influenced by the solvent (Zang et al., 1997), as well as by the proportion of the different xanthophylls and carotenoids in the mix, which could explain the slight difference between our spectrum and the Bernstein's et al. spectra acquired in olive oil.

These differences in pigmentation may be an ideal diagnostic tool to discriminate between chickens with high or low DE because of its ease of implementation. It has previously been used in chickens to determine their status regarding coccidiosis infection, a parasitic disease of the intestine (Hamzic et al., 2015). In cases of infection, the intestinal absorption of pigments is disturbed and pigments are partially lost by leakage through the damaged wall of the cecum, leading to a reduced level of xanthophylls and carotenoids in the plasma from infected broilers (Ruff et al., 1974; Ruff and Fuller, 1975). The mechanism leading to the difference in serum pigmentation between D+ and D− could be similar as absorption capacity is probably lower in D− than in D+, including for pigments. Indeed, de Verdal et al. (2010) showed that the structure of the intestine was more
<table><caption>**Table 1.** Growth, feed conversion ratio, and DE traits of the D+ and D− line samples.</caption><thead><tr><th>Traits<sup>1</sup></th><th>D+ (n = 20)<sup>2</sup></th><th>D− (n = 20)</th><th>P value<sup>3</sup></th></tr></thead><tbody><tr><td>Body Weight at 17d (g)</td><td>329 ± 30</td><td>350 ± 31</td><td>0.06</td></tr><tr><td>FCR 17/26 d (g.g<sup>-1</sup>)</td><td>1.85 ± 0.12</td><td>2.14 ± 0.21</td><td>7.57 × 10<sup>-6</sup></td></tr><tr><td>CDU<sub>S</sub> (%)</td><td>98.32 ± 0.65</td><td>89.30 ± 9.41</td><td>2.76 × 10<sup>-10</sup></td></tr><tr><td>CDU<sub>L</sub> (%)</td><td>90.48 ± 3.71</td><td>74.17 ± 15.38</td><td>8.34 × 10<sup>-8</sup></td></tr><tr><td>CDU<sub>N</sub> (%)</td><td>81.53 ± 2.54</td><td>74.75 ± 4.28</td><td>1.69 × 10<sup>-7</sup></td></tr><tr><td>AMEn (kcal.kg<sup>-1</sup> DM)</td><td>3477.58 ± 92.27</td><td>3053.65 ± 345.76</td><td>4.85 × 10<sup>-7</sup></td></tr></tbody></table>

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

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Significance of the difference between the two lines for the absorption profiles of lipophilic part of serum (1a) and reference profile of the xanthophylls reprinted from Bernstein et al. (Copyright 2010, with permission from Elsevier; 1b).</figcaption></figure>

fragile in D−, which could favor losses by leakage in this line. However, it was unlikely that these differences were the result of a difference in sensitivity to coccidiosis since the diet included 0.05% of coccidiostats.

Spectrophotometry cannot identify with certainty the molecules responsible for these differences between the two lines without further biochemical characterization. For example, bilirubin, a product of heme catabolism present in serum, also absorbs at 460 nm and cannot be distinguished from xanthophylls and carotenoids by this method. However, xanthophyll and carotenoid, but not bilirubin, transfer to nonpolar solvents such as hexane (Brodersen, 1979). To obtain a better understanding, the LC-HRMS spectra of the

lipophilic fraction of the sera from the same animals were thus acquired.

## Serum Lipidome and Digestive Efficiency

The LC-HRMS (ESI+) spectra were processed using XCMS and 13,420 spectral features were identified (not including isotopic peaks). A total of 2,260 peaks were conserved after removal of those with poor qualities. The PCA analysis, when using all the data samples, showed that the five QC samples were clustered, which validates the stability of the analysis quality during the LC-HRMS acquisition. The PCA plots revealed that three samples could be considered to be outliers (2 D− and 1 D+) and thus were removed from the subsequent PLS-DA analysis (data not shown).

The PLS-DA model resulted in two predictive components with a cross-validated predictive ability (Q<fnref n="2" />) of 0.80 and an overall proportion of the variation in DE explained by the model (R<fnref n="2" />_Y) of 0.82. The reliability of this model assessed by CV-ANOVA was good, with a *P* value of 9.31 × 10⁻¹¹. The scatter plot resulting from the PLS-DA is shown in Figure 2a. The PLS-DA showed a clear separation between the D+ and the D− birds with high quality of prediction (Q<fnref n="2" />), suggesting lipidome differences between the two lines. The model included 10 features each characterized by an exact mass and a retention time, which were listed in Table 2 as well as their importance in the model (VIP) and their loadings on the two components. The annotation of those peaks was performed based on the exact mass. That might have affected the precision of the annotation and could not distinguish some isomers or isobars. However, the difference between the exact mass and theoretical mass of the metabolite (expressed in ppm in Table 2) was low, which bring more confidence into this annotation of the metabolites. The putative molecules included in the model were lipids [ceramide (Cer), diglyceride (DG), phosphocholine (PC), phosphoethanolamide (PE), 3-O-(2-O-(2E-decenoyl)-alpha-L-rhamnopyranosyl- (1-2)-alpha-L-rhamnopyranosyl)-3-hydroxydecanoic acid, kamlolenic acid] as well as a terpenoid (24-norolea-3,12-diene), a steroid (22:5 cholesteryl ester), and xanthophylls (lutein/zeaxanthin), the latter being the most

<footnote n="1">¹FCR 17/26: feed conversion ratio between 17 and 26 d; CDU<sub>S, L, N</sub>: coefficients of digestive use of starch, lipids and nitrogen; AMEn: metabolizable energy corrected to zero nitrogen balance.</footnote>

<footnote n="2">²Values are expressed as mean ± SD.</footnote>

<footnote n="3">³Two-tailed Mann-Whitney's median value equality t-test.</footnote>
<table><caption>**Table 1.** Growth, feed conversion ratio, and DE traits of the D+ and D− line samples.</caption><thead><tr><th>Traits<sup>1</sup></th><th>D+ (n = 20)<sup>2</sup></th><th>D− (n = 20)</th><th>P value<sup>3</sup></th></tr></thead><tbody><tr><td>Body Weight at 17d (g)</td><td>329 ± 30</td><td>350 ± 31</td><td>0.06</td></tr><tr><td>FCR 17/26 d (g.g<sup>-1</sup>)</td><td>1.85 ± 0.12</td><td>2.14 ± 0.21</td><td>7.57 × 10<sup>-6</sup></td></tr><tr><td>CDU<sub>S</sub> (%)</td><td>98.32 ± 0.65</td><td>89.30 ± 9.41</td><td>2.76 × 10<sup>-10</sup></td></tr><tr><td>CDU<sub>L</sub> (%)</td><td>90.48 ± 3.71</td><td>74.17 ± 15.38</td><td>8.34 × 10<sup>-8</sup></td></tr><tr><td>CDU<sub>N</sub> (%)</td><td>81.53 ± 2.54</td><td>74.75 ± 4.28</td><td>1.69 × 10<sup>-7</sup></td></tr><tr><td>AMEn (kcal.kg<sup>-1</sup> DM)</td><td>3477.58 ± 92.27</td><td>3053.65 ± 345.76</td><td>4.85 × 10<sup>-7</sup></td></tr></tbody></table>

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

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Significance of the difference between the two lines for the absorption profiles of lipophilic part of serum (1a) and reference profile of the xanthophylls reprinted from Bernstein et al. (Copyright 2010, with permission from Elsevier; 1b).</figcaption></figure>

fragile in D−, which could favor losses by leakage in this line. However, it was unlikely that these differences were the result of a difference in sensitivity to coccidiosis since the diet included 0.05% of coccidiostats.

Spectrophotometry cannot identify with certainty the molecules responsible for these differences between the two lines without further biochemical characterization. For example, bilirubin, a product of heme catabolism present in serum, also absorbs at 460 nm and cannot be distinguished from xanthophylls and carotenoids by this method. However, xanthophyll and carotenoid, but not bilirubin, transfer to nonpolar solvents such as hexane (Brodersen, 1979). To obtain a better understanding, the LC-HRMS spectra of the

lipophilic fraction of the sera from the same animals were thus acquired.

## Serum Lipidome and Digestive Efficiency

The LC-HRMS (ESI+) spectra were processed using XCMS and 13,420 spectral features were identified (not including isotopic peaks). A total of 2,260 peaks were conserved after removal of those with poor qualities. The PCA analysis, when using all the data samples, showed that the five QC samples were clustered, which validates the stability of the analysis quality during the LC-HRMS acquisition. The PCA plots revealed that three samples could be considered to be outliers (2 D− and 1 D+) and thus were removed from the subsequent PLS-DA analysis (data not shown).

The PLS-DA model resulted in two predictive components with a cross-validated predictive ability (Q<fnref n="2" />) of 0.80 and an overall proportion of the variation in DE explained by the model (R<fnref n="2" />_Y) of 0.82. The reliability of this model assessed by CV-ANOVA was good, with a *P* value of 9.31 × 10⁻¹¹. The scatter plot resulting from the PLS-DA is shown in Figure 2a. The PLS-DA showed a clear separation between the D+ and the D− birds with high quality of prediction (Q<fnref n="2" />), suggesting lipidome differences between the two lines. The model included 10 features each characterized by an exact mass and a retention time, which were listed in Table 2 as well as their importance in the model (VIP) and their loadings on the two components. The annotation of those peaks was performed based on the exact mass. That might have affected the precision of the annotation and could not distinguish some isomers or isobars. However, the difference between the exact mass and theoretical mass of the metabolite (expressed in ppm in Table 2) was low, which bring more confidence into this annotation of the metabolites. The putative molecules included in the model were lipids [ceramide (Cer), diglyceride (DG), phosphocholine (PC), phosphoethanolamide (PE), 3-O-(2-O-(2E-decenoyl)-alpha-L-rhamnopyranosyl- (1-2)-alpha-L-rhamnopyranosyl)-3-hydroxydecanoic acid, kamlolenic acid] as well as a terpenoid (24-norolea-3,12-diene), a steroid (22:5 cholesteryl ester), and xanthophylls (lutein/zeaxanthin), the latter being the most

<footnote n="1">¹FCR 17/26: feed conversion ratio between 17 and 26 d; CDU<sub>S, L, N</sub>: coefficients of digestive use of starch, lipids and nitrogen; AMEn: metabolizable energy corrected to zero nitrogen balance.</footnote>

<footnote n="2">²Values are expressed as mean ± SD.</footnote>

<footnote n="3">³Two-tailed Mann-Whitney's median value equality t-test.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Metabolome score plot (a) and contribution plot indicating the contribution of the metabolites identified in PLS-DA (b). The individuals from the D− and D+ lines are represented by blue diamonds and red triangles, respectively, in the score plot. The features outside of 3 standard deviations in the raw data distribution were coloured in black in the contribution plot.</figcaption></figure>

differential between lines with a fold change of 4.63. The contribution plots generated by SIMCA software (Figure 2b) were used to identify the metabolites allowing chicken-line discrimination. The M395T10 feature (24-Noroleana-3,12-diene) was the only one to highly and positively contribute (outside of 3 standard deviations in the raw data distribution) to the D− line, while

M551T5 (lutein/zeaxanthin) was the highest contributor to the D+ line.

## Biomarkers of Digestive Efficiency

The spectrophotometric approach as well as the PLS-DA model fitted on a selection of LC-HRMS features were two sources of biomarkers related to DE with their advantages and inconveniences. Spectrophotometric measurement of the absorption of the sera at 492 nm is easy to implement, but it is not as robust as the measurement of the metabolites included in the PLS-DA, which on the other hand requires a quite expensive analytical procedure. However, the number of features screened could be reduced according to the precision of the diagnosis needed or technical constraints for practical use. For instance, reducing the set of markers to the three features with the highest VIP (M551T5, M395T10, M295T18) did not hugely impact the explicative ability of the model ($R^2_Y = 0.73$) or the predictive ability ($Q^2 = 0.63$). Interestingly, M551T5, the variable with the highest VIP in the PLS-DA, was assigned to lutein/zeaxanthin, which has absorption spectra similar to the differential absorption profile between D+ and D−. Lutein and zeaxanthin, two isomers of $C_{40}H_{56}O_2$, represent 61% to 73% of yellow-orange pigment in corn kernels (O'Hare et al., 2015). They are also present at a lower concentration in wheat and soy kernels (Pardio et al., 2005; Roosevelt et al., 2009).

The LC-HRMS approach applied to lipidic fractions led to a PLS-DA model with similar explicative abilities $R^2_Y$ (0.82 vs 0.77) and higher predictive abilities $Q^2$ (0.80 vs 0.48) than the previous PLS model fitted to metabolites from serum polar extracts identified by NMR to predict the quantitative value of AMEn (Beaucercq et al., 2018). This could be explained by a higher sensitivity of LC-HRMS (nM to pM) in comparison to NMR ($\mu$M to mM) (Goldansaz et al.,

<table><caption>**Table 2.** Serum lipophilic fraction metabolites identified by PLS-DA discriminating between the D+ and D− lines.</caption><thead><tr><th>ID</th><th>Retention time (min)</th><th>Exact mass (m/z)</th><th>VIP<sup>1</sup></th><th>Loading comp.1</th><th>Loading comp.2</th><th>Putative metabolite<sup>2</sup></th></tr></thead><tbody><tr><td>M551T5</td><td>5.31</td><td>551.4250</td><td>1.95</td><td>-0.6269</td><td>0.3988</td><td>Lutein (1 ppm) Zeaxanthin (1 ppm)</td></tr><tr><td>M395T10</td><td>9.85</td><td>395.3669</td><td>1.20</td><td>0.3843</td><td>0.3529</td><td>24-Noroleana-3,12-diene (1 ppm)</td></tr><tr><td>M295T18</td><td>17.91</td><td>295.2265</td><td>0.87</td><td>-0.2848</td><td>-0.1813</td><td>Kamlolenic acid (1 ppm)</td></tr><tr><td>M804T15</td><td>14.99</td><td>803.5135</td><td>0.87</td><td>0.2836</td><td>0.1625</td><td>3-O-(2-O-(2E-decenoyl)-alpha-L-rhamnopyranosyl-(1-2)-alpha-L-rhamnopyranosyl)-3-hydroxydecanoic acid (2 ppm)</td></tr><tr><td>M995T14</td><td>13.57</td><td>994.6450</td><td>0.81</td><td>-0.2673</td><td>-0.0139</td><td>Unknown</td></tr><tr><td>M657T15-3</td><td>14.75</td><td>656.6186</td><td>0.80</td><td>-0.2646</td><td>-0.0462</td><td>Ceramide (18:0/22:0) + isomers (0 ppm) Diglyceride (17:0/20:0) + isomers (0 ppm)</td></tr><tr><td>M738T19-2</td><td>18.63</td><td>737.5611</td><td>0.74</td><td>-0.2387</td><td>-0.2213</td><td>22:5 Cholesteryl ester (3 ppm)</td></tr><tr><td>M845T9-1</td><td>9.35</td><td>844.5229</td><td>0.70</td><td>-0.2286</td><td>-0.1024</td><td>Phosphocholine (16:0/22:6) + isomers (3 ppm)</td></tr><tr><td>M779T12-1</td><td>12.43</td><td>778.5738</td><td>0.68</td><td>-0.2152</td><td>-0.2445</td><td>Phosphoethanolamide(22:4/18:1) + isomers (1 ppm)</td></tr><tr><td>M579T12-2</td><td>12.00</td><td>578.5141</td><td>0.67</td><td>0.0273</td><td>0.7467</td><td>Ceramide (14:2/22:1) + isomers (0 ppm)</td></tr></tbody></table>

<footnote n="1">¹Variable importance in the projection.</footnote>

<footnote n="2">²Putative annotation of the peaks based on their exact masses. The differences between the mass determined by spectrometry and the theoretical mass of the metabolites are indicated in ppm.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Metabolome score plot (a) and contribution plot indicating the contribution of the metabolites identified in PLS-DA (b). The individuals from the D− and D+ lines are represented by blue diamonds and red triangles, respectively, in the score plot. The features outside of 3 standard deviations in the raw data distribution were coloured in black in the contribution plot.</figcaption></figure>

differential between lines with a fold change of 4.63. The contribution plots generated by SIMCA software (Figure 2b) were used to identify the metabolites allowing chicken-line discrimination. The M395T10 feature (24-Noroleana-3,12-diene) was the only one to highly and positively contribute (outside of 3 standard deviations in the raw data distribution) to the D− line, while

M551T5 (lutein/zeaxanthin) was the highest contributor to the D+ line.

## Biomarkers of Digestive Efficiency

The spectrophotometric approach as well as the PLS-DA model fitted on a selection of LC-HRMS features were two sources of biomarkers related to DE with their advantages and inconveniences. Spectrophotometric measurement of the absorption of the sera at 492 nm is easy to implement, but it is not as robust as the measurement of the metabolites included in the PLS-DA, which on the other hand requires a quite expensive analytical procedure. However, the number of features screened could be reduced according to the precision of the diagnosis needed or technical constraints for practical use. For instance, reducing the set of markers to the three features with the highest VIP (M551T5, M395T10, M295T18) did not hugely impact the explicative ability of the model ($R^2_Y = 0.73$) or the predictive ability ($Q^2 = 0.63$). Interestingly, M551T5, the variable with the highest VIP in the PLS-DA, was assigned to lutein/zeaxanthin, which has absorption spectra similar to the differential absorption profile between D+ and D−. Lutein and zeaxanthin, two isomers of $C_{40}H_{56}O_2$, represent 61% to 73% of yellow-orange pigment in corn kernels (O'Hare et al., 2015). They are also present at a lower concentration in wheat and soy kernels (Pardio et al., 2005; Roosevelt et al., 2009).

The LC-HRMS approach applied to lipidic fractions led to a PLS-DA model with similar explicative abilities $R^2_Y$ (0.82 vs 0.77) and higher predictive abilities $Q^2$ (0.80 vs 0.48) than the previous PLS model fitted to metabolites from serum polar extracts identified by NMR to predict the quantitative value of AMEn (Beaucercq et al., 2018). This could be explained by a higher sensitivity of LC-HRMS (nM to pM) in comparison to NMR ($\mu$M to mM) (Goldansaz et al.,

<table><caption>**Table 2.** Serum lipophilic fraction metabolites identified by PLS-DA discriminating between the D+ and D− lines.</caption><thead><tr><th>ID</th><th>Retention time (min)</th><th>Exact mass (m/z)</th><th>VIP<sup>1</sup></th><th>Loading comp.1</th><th>Loading comp.2</th><th>Putative metabolite<sup>2</sup></th></tr></thead><tbody><tr><td>M551T5</td><td>5.31</td><td>551.4250</td><td>1.95</td><td>-0.6269</td><td>0.3988</td><td>Lutein (1 ppm) Zeaxanthin (1 ppm)</td></tr><tr><td>M395T10</td><td>9.85</td><td>395.3669</td><td>1.20</td><td>0.3843</td><td>0.3529</td><td>24-Noroleana-3,12-diene (1 ppm)</td></tr><tr><td>M295T18</td><td>17.91</td><td>295.2265</td><td>0.87</td><td>-0.2848</td><td>-0.1813</td><td>Kamlolenic acid (1 ppm)</td></tr><tr><td>M804T15</td><td>14.99</td><td>803.5135</td><td>0.87</td><td>0.2836</td><td>0.1625</td><td>3-O-(2-O-(2E-decenoyl)-alpha-L-rhamnopyranosyl-(1-2)-alpha-L-rhamnopyranosyl)-3-hydroxydecanoic acid (2 ppm)</td></tr><tr><td>M995T14</td><td>13.57</td><td>994.6450</td><td>0.81</td><td>-0.2673</td><td>-0.0139</td><td>Unknown</td></tr><tr><td>M657T15-3</td><td>14.75</td><td>656.6186</td><td>0.80</td><td>-0.2646</td><td>-0.0462</td><td>Ceramide (18:0/22:0) + isomers (0 ppm) Diglyceride (17:0/20:0) + isomers (0 ppm)</td></tr><tr><td>M738T19-2</td><td>18.63</td><td>737.5611</td><td>0.74</td><td>-0.2387</td><td>-0.2213</td><td>22:5 Cholesteryl ester (3 ppm)</td></tr><tr><td>M845T9-1</td><td>9.35</td><td>844.5229</td><td>0.70</td><td>-0.2286</td><td>-0.1024</td><td>Phosphocholine (16:0/22:6) + isomers (3 ppm)</td></tr><tr><td>M779T12-1</td><td>12.43</td><td>778.5738</td><td>0.68</td><td>-0.2152</td><td>-0.2445</td><td>Phosphoethanolamide(22:4/18:1) + isomers (1 ppm)</td></tr><tr><td>M579T12-2</td><td>12.00</td><td>578.5141</td><td>0.67</td><td>0.0273</td><td>0.7467</td><td>Ceramide (14:2/22:1) + isomers (0 ppm)</td></tr></tbody></table>

<footnote n="1">¹Variable importance in the projection.</footnote>

<footnote n="2">²Putative annotation of the peaks based on their exact masses. The differences between the mass determined by spectrometry and the theoretical mass of the metabolites are indicated in ppm.</footnote>
2017) or a lipophilic fraction more related to DE vari-
ations.

The colorimetric and lipidomic approaches opens up
new possibilities to obtain a quick and easy proxy of DE
that could be generalized to other chicken populations
as metabolome and metabolic responses could be less
dependent on the genetics than single nucleotide poly-
morphism or even the transcriptome in a given physio-
logical state. Indeed, metabolomics and lipidomics are
powerful means to non-invasively predict the quality of
the production (Straadt et al., 2014; Beauclercq et al.,
2016), feed efficiency and residual feed intake, ascer-
tain disease propensity, evaluate dietary responses to
different feeds, fertility, and other important economic
or breeding traits associated with livestock (Goldansaz
et al., 2017). This preliminary study is a first step in
the discovery of DE biomarkers in chickens. The set
of marker (exact mass, retention time, or visible wave-
length) presented in this study might be used as a proxy
for DE. However, their chemical natures should be con-
firmed at the molecular level to get a better understand-
ing of the metabolic pathways influenced by DE; more-
over, the biomarker set needs to be validated on a larger
and independent population and diets with other com-
positions. This metabolic pattern will certainly assist
in the development of more sustainable poultry produc-
tion schemes by allowing a better combined nutritional
and genetic approach to production within the context
of a more variable and complex composition of broiler
diets.

ACKNOWLEDGMENTS

We thank the staff of the Poultry Breeding Facil-
ities (INRA, Unité Expérimentale Pôle Expérimental
Avicole de Tours, INRA, Nouzilly, France) and of the
Rearing Facilities (INRA, Unité Expérimentale Elevage
Alternatif et Santé des Monogastriques, INRA, Le Mag-
neraud, France). This study was supported by Feed-a-
Gene, a project that has received funding from the Eu-
ropean Union's Horizon 2020 research and innovation
program under grant agreement No. 633531.

REFERENCES

Barker, M., and W. Rayens 2003. Partial least squares for discrimi-
nation. J. Chemom. 17:166-173.

Bastianelli, D., L. Bonnal, H. Juin, S. Mignon-Grasteau, F.
Davrieux, and B. Carré 2010. Prediction of the chemical com-
position of poultry excreta by near infrared spectroscopy. J. Near
Infrared Spectrosc. 18:69-77.

Beauclercq, S., L. Nadal-Desbarats, C. Hennequet-Antier, A. Collin, S. Tesseraud, M. Bourin, E. Le Bihan-Duval, and C. Berri 2016. Serum and Muscle Metabolomics for the Prediction of Ultimate pH, a Key Factor for Chicken-Meat Quality. J. Proteome Res. 15:1168-1178.

Beauclercq, S., L. Nadal-Desbarats, C. Hennequet-Antier, I. Gabriel, S. Tesseraud, F. Calenge, E. Le Bihan-Duval, and S. Mignon-Grasteau 2018. Relationships between digestive efficiency and metabolomic profiles of serum and intestinal contents in chickens. Sci. Rep. 8:6678.

Bernstein, P. S., F. C. Delori, S. Richer, F. J. M. van Kuijk, and A.
J. Wenzel 2010. The value of measurement of macular carotenoid

pigment optical densities and distributions in age-related macular
degeneration and other retinal disorders. Vision Res. 50:716-728.

Bligh, E. G., and W. J. Dyer 1959. A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 37:911-917.

Bourdillon, A., B. Carré, L. Conan, J. Duperray, G. Huyghebaert, B. Leclercq, M. Lessire, J. McNab, and J. Wiseman 1990. European reference method for the *in vivo* determination of metabolisable energy with adult cockerels: Reproducibility, effect of food intake and comparison with individual laboratory methods. Br. Poult. Sci. 31:557-565.

Brodersen, R. 1979. Bilirubin. Solubility and interaction with albu-
min and phospholipid. J. Biol. Chem. 254:2364-2369.

Carré, B., A. Idi, S. Maisonnier, J.-P. Melcion, F.-X. Oury, J. Gomez, and P. Pluchard 2002. Relationships between digestibilities of food components and characteristics of wheats (Triticum aestivum) introduced as the only cereal source in a broiler chicken diet. Br. Poult. Sci. 43:404-415.

Carré, B., S. Mignon-Grasteau, J. Besnard, N. Rougière, H. Juin, and D. Bastianelli 2008. The D+ and D-"Digestion" chicken lines selected for divergent digestion efficiency on a wheat-based died. Page 210 in World Poultry Congress. Brisbane, Australia.

de Verdal, H., S. Mignon-Grasteau, C. Jeulin, E. Le Bihan-Duval, M. Leconte, S. Mallet, C. Martin, A. Narcy, H. J. M., and C. B. 2010. Digestive tract measurements and histological adaptation in broiler lines divergently selected for digestive efficiency. Poult. Sci. 89:1955-1961.

de Verdal, H., A. Narcy, D. Bastianelli, H. Chapuis, N. Même, S. Urvoix, E. Le Bihan-Duval, and S. Mignon-Grasteau 2011. Improving the efficiency of feed utilization in poultry by selection. 1. Genetic parameters of anatomy of the gastro-intestinal tract and digestive efficiency. BMC Genet. 12:59.

Eriksson, L., J. Trygg, and S. Wold 2008. CV-ANOVA for signifi-
cance testing of PLS and OPLS models. J. Chemom. 22:594-600.

Fahy, E., S. Subramaniam, H. A. Brown, C. K. Glass, A. H. Merrill, R. C. Murphy, C. R. H. Raetz, D. W. Russell, Y. Seyama, W. Shaw, T. Shimizu, F. Spener, G. van Meer, M. S. VanNieuwenhze, S. H. White, J. L. Witztum, and E. A. Dennis 2005. A comprehensive classification system for lipids. J. Lipid Res. 46:839-861.

Fahy, E., S. Subramaniam, R. C. Murphy, M. Nishijima, C. R. H. Raetz, T. Shimizu, F. Spener, G. van Meer, M. J. O. Wakelam, and E. A. Dennis 2009. Update of the LIPID MAPS comprehensive classification system for lipids. J. Lipid Res. 50:9-14.

Goldansaz, S. A., A. C. Guo, T. Sajed, M. A. Steele, G. S. Plastow, and D. S. Wishart 2017. Livestock metabolomics and the livestock metabolome: A systematic review (MFW te Pas, Ed.). PLoS One 12:e0177675.

Hamzic, E., B. Bed'Hom, H. Juin, R. Hawken, M. S. Abrahamsen, J. M. Elsen, B. Servin, M. H. Pinard-van der Laan, and O. Demeure 2015. Large-scale investigation of the parameters in response to challenge in broilers. J. Anim. Sci. 93:1830.

Mignon-Grasteau, S., N. Muley, D. Bastianelli, J. Gomez, A. Péron, N. Sellier, N. Millet, J. Besnard, J. M. Hallouis, and B. Carré 2004. Heritability of digestibilities and divergent selection for digestion ability in growing chicks fed a wheat diet. Poult. Sci. 83:860-867.

O'Hare, T. J., K. J. Fanning, and I. F. Martin 2015. Zeaxanthin biofortification of sweet-corn and factors affecting zeaxanthin accumulation and colour change. Arch. Biochem. Biophys. 572:184-187.

Pardio, V. T., L. A. Landin, K. N. Waliszewski, F. Perez-Gil, L. Diaz, and B. Hernandez 2005. The effect of soybean soapstock on the quality parameters and fatty acid composition of the hen egg yolk. Poult. Sci. 84:148-157.

Riffard, C., S. Gallot, and P. Magdelaine 2011. Performances techniques et coûts de production en volailles de chair, poulettes et pondeuses. Résultats 2010. Tech. Rep., ITAVI 57.

Roose, M., J. Kahl, and A. Ploeger 2009. Influence of the farming system on the xanthophyll content of soft and hard wheat. J. Agric. Food Chem. 57:182-188.

Ruff, M. D., and H. L. Fuller 1975. Some mechanisms of reduction of carotenoid levels in chickens infected with Eimeria acervulina or E. tenella. J. Nutr. 105:1447-1456.
2017) or a lipophilic fraction more related to DE vari-
ations.

The colorimetric and lipidomic approaches opens up
new possibilities to obtain a quick and easy proxy of DE
that could be generalized to other chicken populations
as metabolome and metabolic responses could be less
dependent on the genetics than single nucleotide poly-
morphism or even the transcriptome in a given physio-
logical state. Indeed, metabolomics and lipidomics are
powerful means to non-invasively predict the quality of
the production (Straadt et al., 2014; Beauclercq et al.,
2016), feed efficiency and residual feed intake, ascer-
tain disease propensity, evaluate dietary responses to
different feeds, fertility, and other important economic
or breeding traits associated with livestock (Goldansaz
et al., 2017). This preliminary study is a first step in
the discovery of DE biomarkers in chickens. The set
of marker (exact mass, retention time, or visible wave-
length) presented in this study might be used as a proxy
for DE. However, their chemical natures should be con-
firmed at the molecular level to get a better understand-
ing of the metabolic pathways influenced by DE; more-
over, the biomarker set needs to be validated on a larger
and independent population and diets with other com-
positions. This metabolic pattern will certainly assist
in the development of more sustainable poultry produc-
tion schemes by allowing a better combined nutritional
and genetic approach to production within the context
of a more variable and complex composition of broiler
diets.

ACKNOWLEDGMENTS

We thank the staff of the Poultry Breeding Facil-
ities (INRA, Unité Expérimentale Pôle Expérimental
Avicole de Tours, INRA, Nouzilly, France) and of the
Rearing Facilities (INRA, Unité Expérimentale Elevage
Alternatif et Santé des Monogastriques, INRA, Le Mag-
neraud, France). This study was supported by Feed-a-
Gene, a project that has received funding from the Eu-
ropean Union's Horizon 2020 research and innovation
program under grant agreement No. 633531.

REFERENCES

Barker, M., and W. Rayens 2003. Partial least squares for discrimi-
nation. J. Chemom. 17:166-173.

Bastianelli, D., L. Bonnal, H. Juin, S. Mignon-Grasteau, F.
Davrieux, and B. Carré 2010. Prediction of the chemical com-
position of poultry excreta by near infrared spectroscopy. J. Near
Infrared Spectrosc. 18:69-77.

Beauclercq, S., L. Nadal-Desbarats, C. Hennequet-Antier, A. Collin, S. Tesseraud, M. Bourin, E. Le Bihan-Duval, and C. Berri 2016. Serum and Muscle Metabolomics for the Prediction of Ultimate pH, a Key Factor for Chicken-Meat Quality. J. Proteome Res. 15:1168-1178.

Beauclercq, S., L. Nadal-Desbarats, C. Hennequet-Antier, I. Gabriel, S. Tesseraud, F. Calenge, E. Le Bihan-Duval, and S. Mignon-Grasteau 2018. Relationships between digestive efficiency and metabolomic profiles of serum and intestinal contents in chickens. Sci. Rep. 8:6678.

Bernstein, P. S., F. C. Delori, S. Richer, F. J. M. van Kuijk, and A. J. Wenzel 2010. The value of measurement of macular carotenoid

pigment optical densities and distributions in age-related macular
degeneration and other retinal disorders. Vision Res. 50:716–728.

Bligh, E. G., and W. J. Dyer 1959. A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 37:911-917.

Bourdillon, A., B. Carré, L. Conan, J. Duperray, G. Huyghebaert, B. Leclercq, M. Lessire, J. McNab, and J. Wiseman 1990. European reference method for the *in vivo* determination of metabolisable energy with adult cockerels: Reproducibility, effect of food intake and comparison with individual laboratory methods. Br. Poult. Sci. 31:557-565.

Brodersen, R. 1979. Bilirubin. Solubility and interaction with albu-
min and phospholipid. J. Biol. Chem. 254:2364-2369.

Carré, B., A. Idi, S. Maisonnier, J.-P. Melcion, F.-X. Oury, J. Gomez, and P. Pluchard 2002. Relationships between digestibilities of food components and characteristics of wheats (Triticum aestivum) introduced as the only cereal source in a broiler chicken diet. Br. Poult. Sci. 43:404–415.

Carré, B., S. Mignon-Grasteau, J. Besnard, N. Rougière, H. Juin, and D. Bastianelli 2008. The D+ and D-“Digestion” chicken lines selected for divergent digestion efficiency on a wheat-based died. Page 210 in World Poultry Congress. Brisbane, Australia.

de Verdal, H., S. Mignon-Grasteau, C. Jeulin, E. Le Bihan-Duval, M. Leconte, S. Mallet, C. Martin, A. Narcy, H. J. M., and C. B. 2010. Digestive tract measurements and histological adaptation in broiler lines divergently selected for digestive efficiency. Poult. Sci. 89:1955-1961.

de Verdal, H., A. Narcy, D. Bastianelli, H. Chapuis, N. Même, S. Urvoix, E. Le Bihan-Duval, and S. Mignon-Grasteau 2011. Improving the efficiency of feed utilization in poultry by selection. 1. Genetic parameters of anatomy of the gastro-intestinal tract and digestive efficiency. BMC Genet. 12:59.

Eriksson, L., J. Trygg, and S. Wold 2008. CV-ANOVA for signifi-
cance testing of PLS and OPLS models. J. Chemom. 22:594–600.

Fahy, E., S. Subramaniam, H. A. Brown, C. K. Glass, A. H. Merrill, R. C. Murphy, C. R. H. Raetz, D. W. Russell, Y. Seyama, W. Shaw, T. Shimizu, F. Spener, G. van Meer, M. S. VanNieuwenhze, S. H. White, J. L. Witztum, and E. A. Dennis 2005. A comprehensive classification system for lipids. J. Lipid Res. 46:839–861.

Fahy, E., S. Subramaniam, R. C. Murphy, M. Nishijima, C. R. H. Raetz, T. Shimizu, F. Spener, G. van Meer, M. J. O. Wakelam, and E. A. Dennis 2009. Update of the LIPID MAPS comprehensive classification system for lipids. J. Lipid Res. 50:9–14.

Goldansaz, S. A., A. C. Guo, T. Sajed, M. A. Steele, G. S. Plastow, and D. S. Wishart 2017. Livestock metabolomics and the livestock metabolome: A systematic review (MFW te Pas, Ed.). PLoS One 12:e0177675.

Hamzic, E., B. Bed'Hom, H. Juin, R. Hawken, M. S. Abrahamsen, J. M. Elsen, B. Servin, M. H. Pinard-van der Laan, and O. Demeure 2015. Large-scale investigation of the parameters in response to challenge in broilers. J. Anim. Sci. 93:1830.

Mignon-Grasteau, S., N. Muley, D. Bastianelli, J. Gomez, A. Péron, N. Sellier, N. Millet, J. Besnard, J. M. Hallouis, and B. Carré 2004. Heritability of digestibilities and divergent selection for digestion ability in growing chicks fed a wheat diet. Poult. Sci. 83:860-867.

O'Hare, T. J., K. J. Fanning, and I. F. Martin 2015. Zeaxanthin biofortification of sweet-corn and factors affecting zeaxanthin accumulation and colour change. Arch. Biochem. Biophys. 572:184-187.

Pardio, V. T., L. A. Landin, K. N. Waliszewski, F. Perez-Gil, L. Diaz, and B. Hernandez 2005. The effect of soybean soapstock on the quality parameters and fatty acid composition of the hen egg yolk. Poult. Sci. 84:148-157.

Riffard, C., S. Gallot, and P. Magdelaine 2011. Performances techniques et coûts de production en volailles de chair, poulettes et pondeuses. Résultats 2010. Tech. Rep., ITAVI 57.

Roose, M., J. Kahl, and A. Ploeger 2009. Influence of the farming system on the xanthophyll content of soft and hard wheat. J. Agric. Food Chem. 57:182-188.

Ruff, M. D., and H. L. Fuller 1975. Some mechanisms of reduction of carotenoid levels in chickens infected with Eimeria acervulina or E. tenella. J. Nutr. 105:1447-1456.
Ruff, M. D., W. M. Reid, and J. K. Johnson 1974. Lowered blood carotenoid levels in chickens infected with coccidia. Poult. Sci. 53:1801-1809.

Smith, C. A., G. O'Maille, E. J. Want, C. Qin, S. A. Trauger, T. R. Brandon, D. E. Custodio, R. Abagyan, and G. Siuzdak 2005. METLIN: a metabolite mass spectral database. Ther. Drug Monit. 27:747-751.

Straadt, I. K., M. D. Aaslyng, and H. C. Bertram 2014. An NMR-based metabolomics study of pork from different crossbreeds and relation to sensory perception. Meat Sci. 96:719-728.

Tautenhahn, R., G. J. Patti, D. Rinehart, and G. Siuzdak 2012. XCMS online: a web-based platform to process untargeted metabolomic data. Anal. Chem. 84:5035-5039.

Want, E. J., P. Masson, F. Michopoulos, I. D. Wilson, G. Theodoridis, R. S. Plumb, J. Shockcor, N. Loftus, E. Holmes, and J. K. Nicholson 2012. Global metabolic profiling of animal and human tissues via UPLC-MS. Nat. Protoc. 8:17-32.

Zang, L. Y., O. Sommerburg, and F. J. van Kuijk 1997. Absorbance changes of carotenoids in different solvents. Free Radic. Biol. Med. 23:1086-1089.
Ruff, M. D., W. M. Reid, and J. K. Johnson 1974. Lowered blood carotenoid levels in chickens infected with coccidia. Poult. Sci. 53:1801-1809.

Smith, C. A., G. O'Maille, E. J. Want, C. Qin, S. A. Trauger, T. R. Brandon, D. E. Custodio, R. Abagyan, and G. Siuzdak 2005. METLIN: a metabolite mass spectral database. Ther. Drug Monit. 27:747-751.

Straadt, I. K., M. D. Aaslyng, and H. C. Bertram 2014. An NMR-based metabolomics study of pork from different crossbreeds and relation to sensory perception. Meat Sci. 96:719-728.

Tautenhahn, R., G. J. Patti, D. Rinehart, and G. Siuzdak 2012. XCMS online: a web-based platform to process untargeted metabolomic data. Anal. Chem. 84:5035-5039.

Want, E. J., P. Masson, F. Michopoulos, I. D. Wilson, G. Theodoridis, R. S. Plumb, J. Shockcor, N. Loftus, E. Holmes, and J. K. Nicholson 2012. Global metabolic profiling of animal and human tissues via UPLC-MS. Nat. Protoc. 8:17-32.

Zang, L. Y., O. Sommerburg, and F. J. van Kuijk 1997. Absorbance changes of carotenoids in different solvents. Free Radic. Biol. Med. 23:1086-1089.