Skip to main content
Download PDF

Identification of biomarkers associated with the feed efficiency by metabolomics profiling: results from the broiler lines divergent for high or low abdominal fat content

Journal of Animal Science and Biotechnology volume 13, Article number: 122 (2022) Cite this article

Abstract

Background

Improving feed efficiency (FE) is one of the main objectives in broiler breeding. It is difficult to directly measure FE traits, and breeders hence have been trying to identify biomarkers for the indirect selection and improvement of FE traits. Metabolome is the "bridge" between genome and phenome. The metabolites may potentially account for more of the phenotypic variation and can suitably serve as biomarkers for selecting FE traits. This study aimed to identify plasma metabolite markers for selecting high-FE broilers. A total of 441 birds from Northeast Agricultural University broiler lines divergently selected for abdominal fat content were used to analyze plasma metabolome and estimate the genetic parameters of differentially expressed metabolites.

Results

The results identified 124 differentially expressed plasma metabolites (P < 0.05) between the lean line (high-FE birds) and the fat line (low-FE birds). Among these differentially expressed plasma metabolites, 44 were found to have higher positive or negative genetic correlations with FE traits (|rg| ≥ 0.30). Of these 44 metabolites, 14 were found to display moderate to high heritability estimates (h2 ≥ 0.20). However, among the 14 metabolites, 4 metabolites whose physiological functions have not been reported were excluded. Ultimately, 10 metabolites were suggested to serve as the potential biomarkers for breeding the high-FE broilers. Based on the physiological functions of these metabolites, reducing inflammatory and improving immunity were proposed to improve FE and increase production efficiency.

Conclusions

According to the pipeline for the selection of the metabolite markers established in this study, it was suggested that 10 metabolites including 7-ketocholesterol, dimethyl sulfone, epsilon-(gamma-glutamyl)-lysine, gamma-glutamyltyrosine, 2-oxoadipic acid, L-homoarginine, testosterone, adenosine 5'-monophosphate, adrenic acid, and calcitriol could be used as the potential biomarkers for breeding the "food-saving broilers".

Introduction

Livestock farming serves to provide high-quality proteins to humans but requires the consumption of a lot of grain resources, intensifying the contradiction of "Competition for Food between Human and Animal" [1]. Feeding livestock with high feed efficiency (FE) can beneficially reduce food consumption, saving the expenditure on livestock production and contributing to the sustained development of the environment [2, 3]. Breeders have always focused on improving FE of broilers [4]. The direct selection of FE traits has been the most effective method for improving FE. However, direct selection is based on accurately measuring feed intake, and implementing this method in breeding farms is difficult [5]. Breeders are attempting to find easily measurable indicators for the rapid and affordable estimation of the individual FE traits.

The metabolome is often regarded as a bridge between the genome and the phenome [6]. As downstream of the gene regulatory network and the protein interaction network, metabolites may provide detailed biological terminal information on the life process [7]. In recent years, metabolomics has evolved as an increasingly popular "omics" approach for revealing the relationship between genetics and phenotypes [7]. Researchers can discover novel biomarkers by analyzing the changes in metabolite expression profiling to further understand the vital metabolic pathways associated with the traits. Emerging evidence indicated the use of metabolites as metabolic markers for selecting FE traits in mammals. A study by Novais et al. evaluated the feasibility of predicting FE traits through serum metabolomics in young Nellore cattle [8]. This study found the retinal, progesterone, and stearic acid as suitable biomarkers for indirectly selecting FE traits, and the vitamin A metabolism pathway was identified as an important pathway related to FE traits [8]. A study by Carmelo et al. found choline and pyridoxamine as the hub-metabolites related to FE traits in Duroc pigs, at the same time, cholesterol sulfate, thiamine, L-methionine, and chenodeoxycholate were identified as the hub-metabolites related to FE traits in Landrace pigs, suggesting that these metabolites served as biomarkers for selecting FE traits [9]. Therefore, metabolomics can not only help investigate the relationship between metabolites and phenotypes but also can be used for identifying biomarkers for selecting high-FE animals. However, the metabolomics study of FE traits in broilers remains poorly understood.

The Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) were selected based on the abdominal fat percentage (AFP) and plasma concentration of very low-density lipoprotein (VLDL) since 1996 [10]. A previous study identified the feed conversion ratio (FCR) and residual feed intake (RFI) to be significantly higher in the fat line than in the lean line [11]. In the present study, firstly, the differences in plasma metabolome between the lean and fat lines were investigated. Subsequently, the genetic parameters of the concentrations of differentially expressed plasma metabolites were estimated. Finally, the metabolite biomarkers that could be potentially used to select high-FE broilers were identified.

Materials and methods

Experimental population

This study used 441 birds from the 23rd generation of the NEAUHLF lines, which were selected based on AFP and plasma concentration of VLDL [10]. All the birds were housed in the similar environmental conditions and had free access to feed and water up to 7 weeks from hatching. The commercial corn-soybean-based diets, meeting all the requirements of the National Research Council (NRC, 1994) [12], were used in the study. The birds received a starter feed (3000 kcal ME/kg and 210 g/kg CP) from hatching to 3 weeks. Then, the birds were fed a grower diet (3100 kcal ME/kg and 190 g/kg CP) from 4 weeks to slaughter. Each bird was placed in an individual cage at 27 days of age, and the individual feed intake (FI) was recorded at 29—49 days of age.

Trait measurement

The total FI of each bird was recorded from 4 weeks old to 7 weeks old. The body weights were measured at 4 weeks (BW4) and 7 weeks (BW7) of age. The gain in body weights (BWG) was calculated by subtracting the BW4 from BW7. The metabolic mid-test body weight (MMBW) was calculated by taking the 0.75th power of the middle body weight during the specific period [(BW4 + BW7)/2]. The FCR and RFI were calculated using the following equations:

$$FCR=FI/BWG$$
$$RFI=FI-\left(b_0+b_1MMBW+b_2BWG\right)$$

where b0 is the intercept and b1, b2 are partial regression coefficients of FI for MMBW and BWG, respectively.

Metabolites extraction

The whole-blood samples were collected and immediately centrifuged at 3000 × g for 10 min at 25 °C. The supernatants (100 μL plasma) were placed in the EP tubes and resuspended in prechilled 80% methanol by vortexing. Then, the samples were incubated on ice for 5 min and centrifuged at 15,000 × g for 20 min at 4 °C. The supernatant was diluted to a final concentration of 53% methanol using LC–MS grade water. The samples were transferred to a fresh EP tube and centrifuged at 15,000 × g for 20 min at 4 °C. Finally, the supernatant was injected into the LC–MS/MS system for subsequent analysis.

UHPLC/MS analysis

The UHPLC/MS was analyzed using a Vanquish UHPLC system (ThermoFisher, Germany) coupled with an Orbitrap Q Exactive TMHF-X mass spectrometer (Thermo Fisher, Germany) in Novogene Co., Ltd. (Beijing, China). The samples were injected onto a Hypesil Gold column (100 mm × 2.1 mm, 1.9 μm) using a 17-min linear gradient at a flow rate of 0.2 mL/min. The eluents for the positive polarity mode were eluent A (0.1% formic acid in water) and eluent B (methanol). The eluents for the negative polarity mode were eluent A (5 mmol/L ammonium acetate, pH 9.0) and eluent B (methanol). The solvent gradient was set as follows: 2% B, 1.5 min; 2%—100% B, 12.0 min; 100% B, 14.0 min; 100%—2% B, 14.0 min; and 2% B, 17.0 min. The Q Exactive™ HF-X mass spectrometer was operated in the positive/negative polarity mode using a spray voltage of 3.2 kV, capillary temperature of 320 °C, sheath gas flow rate of 40 arb, and aux gas flow rate of 10 arb.

The raw data files generated by UHPLC/MS were analyzed using the Compound Discoverer 3.1 (CD 3.1, ThermoFisher) for peak alignment, peak picking, and quantitation for each metabolite. The main parameters were set as follows: retention time tolerance, 0.2 min; actual mass tolerance, 5 ppm; signal intensity tolerance, 30%; and signal/noise ratio, 3. Subsequently, the peak intensities were normalized to the total spectral intensity and the normalized data were used for predicting the molecular formula based on the additive ions, molecular ion peaks, and fragment ions. Then, the peaks were matched with the mzCloud, mzVault, MassList, and ChemSpider databases for obtaining accurate qualitative and relative quantitative results [13, 14]. These metabolites were annotated using the human metabolome database (HMDB) [9], ignoring the metabolites that did not correspond to the HMDB.

Statistical analysis

The differences among the groups were explored using the partial least square-discriminant analysis (PLS-DA). The quality of PLS-DA model was assessed using the goodness of fit (R2) and goodness of prediction (Q2) in sevenfold cross validation [8]. The variable importance in projection (VIP) was considered for screening out the discriminant metabolites. The statistical significance (P value) was calculated using the Student’s t-test. The metabolites with VIP > 1 and P < 0.05 were considered as differential metabolites between the lean and fat lines [15]. For heatmaps, the data were normalized using the log10 of the concentrations of differential metabolites and plotted using the Pheatmap package in the R (version 4.1.2) environment.

The genetic parameters, including the heritability of the metabolites as well as the genetic and phenotypic correlations between these metabolites and FE trait indices (FCR and RFI), were estimated using the ASReml (version 4.0) software, with line and sex treated as fixed effects. The animal model used for estimating the genetic parameters was as follows:

$$\boldsymbol Y=\boldsymbol X\boldsymbol b+\boldsymbol Z\boldsymbol a+\boldsymbol e$$

where Y represents the vector of the concentration of the plasma metabolites or the phenotypic value of the FE trait indices; b represents the vector of the fixed effect including the population mean, sex effects, and line effects; a represents the vector of the random effects including the genetic effects; e represents the vector of the random residual effects; and X and Z are the incidence matrices for b and a. The random-effects a and e were assumed to follow the normal distributions with a mean of 0. The variances of a and e were assumed to be Var(a) = Ag and Var(e) = Ir, respectively, in which A represents the numerator relationship matrix of all the animals in the pedigree file, g is the additive genetic variance, I is the identity matrix, and r is the residual variance [16]. According to the recommendation by Dong et al., the thresholds in this study for high genetic correlation and moderate to high heritability were set to be greater than 0.30 and 0.20, respectively [17].

Results

Metabolite profile differences between the lean and fat lines

In this study, a total of 941 and 882 m/z features in the positive and negative ion modes were found to match the search databases, respectively (Table 1). These features were annotated using the HMDB, and 284 and 272 metabolites were finally identified in the positive and negative ion modes, respectively, to be used for subsequent analysis. The metabolites, such as lipids and lipid-like molecules, lignans, neolignans and related compounds, nucleosides, nucleotides and analogs, alkaloids and derivatives, organic compounds, phenylpropanoids, and polyketides, were identified, while the lipids and lipid-like molecules comprised the most abundant of all the metabolites (Fig. 1). The results of PLS-DA of all the plasma metabolites indicated that the fat and lean broilers were divided into two independent groups in both the positive and negative ion modes, indicating differences in the metabolic levels between the fat and lean chickens (Fig. 2). The parameters of PLS-DA model in the positive ion mode, including R2X, R2Y, and Q2Y, were established as 0.519, 0.965, and 0.952, respectively. The R2X, R2Y, and Q2Y in the negative ion mode were established as 0.167, 0.962, and 0.947, respectively. The concentrations of the differentially expressed plasma metabolites between the lean and fat birds were shown in Fig. 3. In the positive ion mode, the concentrations of 55 plasma metabolites were found to be significantly different between the lean and fat birds. Of these, 21 metabolites were found to have significantly higher concentrations in the fat birds than in the lean birds, whereas 34 metabolites were found to have significantly lower concentrations in the fat birds than in the lean birds (P < 0.05). In the negative ion mode, the concentrations of 69 plasma metabolites were found to be significantly different between the lean and fat birds. Of these, the concentrations of 45 metabolites were significantly higher in the fat birds than in the lean birds, whereas 24 metabolites were found to have significantly lower concentrations in the fat birds than in the lean birds (P < 0.05).

Table 1 Summary of features identified by analyzing the non-targeted plasma metabolomics
Full size table
Fig. 1
figure 1

Classifications of plasma metabolites in the positive (A) and negative (B) ion modes

Full size image
Fig. 2
figure 2

PLS-DA score plots based on plasma metabolites in the positive (A) and negative (B) ion modes. Ellipse represents the 95% confidence interval of the lean lines (red) and fat lines (blue)

Full size image
Fig. 3
figure 3

Heatmap of differentially expressed plasma metabolites in the positive (A) and negative (B) ion modes. The numbers in the grid are expressed as the log10 of the plasma concentrations of metabolites. HMDB_ID is the serial number for the metabolites established by the HMDB

Full size image

Estimations of the genetic parameters of the differentially expressed plasma metabolites

A large number of differentially expressed plasma metabolites have been identified between the lean and fat lines. However, whether the expression levels of these metabolites can be related to FE traits is unknown, while heritability of differentially expressed plasma metabolites is also unknown. Thus, the genetic and phenotypic correlation coefficients were estimated between differentially expressed plasma metabolites and FE trait indices. Among these differentially expressed metabolites, the genetic and phenotypic correlation coefficients between 21 positive metabolites or 27 negative metabolites and FE trait indices could not be calculated due to no converge. The genetic parameters of the remaining 34 metabolites in the positive ion mode were shown in Table 2, while the remaining 42 metabolites in the negative ion mode were shown in Table 3.

Table 2 Genetic parameters of the differentially expressed plasma metabolites in the positive ion modea
Full size table
Table 3 Genetic parameters of the differentially expressed plasma metabolites in the negative ion modea
Full size table

In the positive ion mode (Table 2), 14 differentially expressed plasma metabolites were found to have higher positive genetic correlations with the FE trait indices (rg ≥ 0.3), whereas 5 differentially expressed plasma metabolites were found to have higher negative genetic correlations with the FE trait indices (rg ≤ –0.3). The differentially expressed plasma metabolites and the FE trait indices showed relatively lower phenotypic correlation coefficients (–0.19 ≤ rp ≤ 0.20). Among the differentially expressed plasma metabolites, the heritability estimates of biliverdin and tyrosyl-phenylalanine were high (0.57 ≤ h2 ≤ 0.70). The heritability estimates of canthaxanthin, 7-ketocholesterol, dimethyl sulfone, hypoletin 8-gentiobioside, pantothenic acid, monoolein, taurine, benzocaine, hippuric acid, epsilon-(gamma-glutamyl)-lysine, guanidine acetic acid, beta-hydroxybutyrate, and gamma-glutamyltyrosine were found to be moderate (0.22 ≤ h2 ≤ 0.38).

In the negative ion mode (Table 3), 11 differentially expressed plasma metabolites had higher positive genetic correlations with the FE trait indices (rg ≥ 0.3), whereas 14 differentially expressed plasma metabolites demonstrated higher negative genetic correlations with the FE trait indices (rg ≤ –0.3). Relatively lower phenotypic correlation coefficients were observed between the differentially expressed plasma metabolites and the FE trait indices (–0.23 ≤ rp ≤ 0.19). Among the differentially expressed plasma metabolites, the heritability estimates of asarone, 2-oxoadipic acid, 2-(formal amino) benzoic acid, L-homoarginine, (+)-alantolactone, mollicellin F, testosterone, adenosine 5′-monophosphate, eugenitin, 4-pyridoxic acid, and 3-hydroxyvaleric acid were found to be high (0.41 ≤ h2 ≤ 0.59). The heritability estimates of testate, adrenic acid, DL-homocystine, (13Z,16Z)-docosadienoic acid, 5-(3′,4′,5′-trihydroxyphenyl)-gamma-valerolactone-O-methyl-O-sulphate, aniline, mesaconic acid, bilirubin, 3b-hydroxy-5-cholenoic acid, 8-oxoguanine, xanthine, and calcitriol were found to be moderate (0.21 ≤ h2 ≤ 0.38).

Discussion

The NEAUHLF lines were successfully bred, and the AFP of broilers in the fat line was found to be 11.02 times higher than in the lean line after selecting the 23rd generation [11]. A previous study showed that FE traits changed significantly with the selection of AFP. The FCR and RFI in the lean line were 2.41 ± 0.01 and –160.08 ± 7.96 g, respectively, while those in the fat line were 2.81 ± 0.01 and 296.3 ± 11.05 g, respectively [11]. Besides, the FE trait indices showed high positive genetic correlations with the abdominal fat weight (AFW) and AFP in the NEAUHLF lines (from 0.49 to 0.58) [11]. These results indicated that the NEAUHLF lines could be used as high-FE and low-FE animal models. Based on our previous study, the present study further explored the differences in the plasma metabolic levels between the NEAUHLF lines using the same experimental chickens, and identified the potentially useful metabolite biomarkers for selecting the "food-saving broilers".

Firstly, the PLS-DA method was used for identifying the differentially expressed plasma metabolites between the lean and fat birds. PLS-DA, as a supervised and multivariate analysis approach for dimensionality reduction, was widely used for metabolomics researchs and has been established as a standard high-dimensional data analysis method [18, 19]. The value of Q2Y could explain the predictive ability of the supervised PLS-DA model and was found to be better if it was close to 1 theoretically [20]. In the present study, the Q2Y was 0.952 in the positive ion mode, and 0.947 in the negative ion mode (Fig. 2), indicating that PLS-DA model was successfully established. The results of the PLS-DA showed that the lean and fat birds were separated into two independent groups (Fig. 2), identifying 124 significantly and differentially expressed metabolites between the lean and fat birds.

The genetic and phenotypic correlations between the metabolites and the FE trait indices were estimated to identify the differentially expressed metabolites associated with the FE traits. In the positive ion mode, 19 metabolites were found to have high genetic correlations with the FE traits indices (|rg| ≥ 0.30) (Table 2). In negative ion mode, 25 metabolites were found to have high genetic correlations with the FE traits indices (|rg| ≥ 0.30) (Table 3). In addition, among the differentially expressed plasma metabolites, 15 metabolites in the positive ion mode (Table 2) and 23 metabolites in the negative ion mode (Table 3) were found to display moderate to high heritability estimates. To date, reports on estimating the genetic parameters of plasma metabolites in birds are few.

The objective of this study was to identify the metabolites that could be used as biomarkers for selecting broilers with high FE. Accordingly, a pipeline containing four criteria was established for determining whether a given metabolite could be used as a marker (Fig. 4). First, the metabolite concentration between the lean and fat birds must be significantly different. Second, the genetic correlation coefficient between the FE trait indices and the metabolite concentration should be relatively high (genetic correlation threshold: |rg| ≥ 0.30). Third, the concentration of the metabolite must be lower in the high-FE birds (lean line), if the genetic correlation coefficients with the FE trait indices were positive, and vice versa. Fourth, the heritability of the metabolite concentration should be moderate to high (heritability threshold: h2 ≥ 0.20). In the positive and negative ion modes between the lean and fat birds, 124 metabolites were found to have significantly different concentrations, meeting the first criterion. Among these 124 metabolites, 44 were found to have relatively high genetic correlations with the FE trait indices, meeting the second criterion. Of these 44 metabolites, 25 met the third criterion. Finally, of these 25 metabolites, 14 met the fourth criterion, which was shown in Table 4 and Fig. 4. In the positive ion mode, five metabolites including 7-ketocholesterol, dimethyl sulfone, hypoletin 8-gentiobioside, epsilon-(gamma-glutamyl)-lysine, and gamma-glutamyltyrosine were found to meet all the four criteria. In the negative ion mode, nine metabolites including 2-oxoadipic acid, L-homoarginine, mollicellin F, testosterone, adenosine 5′-monophosphate, adrenic acid, (13Z,16Z)-docosadienoic acid, 3b-hydroxy-5-cholenoic acid, and calcitriol were found to meet all four criteria. Therefore, these 14 metabolites were presumed to have the potential to be used as the biomarkers for selecting "food-saving broilers".

Fig. 4
figure 4

Pipeline for the selection of the metabolite biomarkers

Full size image
Table 4 Summary of 14 metabolite biomarkers for selecting the high-FE birds
Full size table

The published literature was searched, and the biological functions of 14 metabolites were summarized to rationally explain the physiological basis of the FE traits. However, among the 14 metabolites, no reports were available on the physiological functions of 4 metabolites, including hypoletin 8-gentiobioside, mollicellin F, (13Z,16Z)-docosadienoic acid, and 3b-hydroxy-5-cholenoic acid. Therefore, the effects of these 4 metabolites on FE traits of broilers could not be speculated. The remaining 10 metabolites were suggested to be biomarkers after analyzing the relationship between these metabolites and FE traits.

7-Ketocholesterol is a pro-inflammatory oxysterol that possibly activates several kinase signaling pathways via multiple transcription factors for inducing the cytokines and intracellular effectors causing cell death [21]. 7-Ketocholesterol damages the vascular endothelial cells by inducing inflammatory responses and apoptosis, elevating the risk of cardiovascular diseases [22]. Recent studies have shown the inflammatory responses to reduce FE in animals. In pigs, the inflammatory reaction could reduce FE [23]; a similar phenomenon was observed in cattle [24]. The fat line in this study had a higher concentration of 7-ketocholesterol than that in the lean line. A previous study showed that the FE in the fat line was significantly lower than that in the lean line [11]. Based on the reports from other studies and the results of this study, the fat line was suggested to have a stronger inflammatory response than the lean line and, therefore, a lower FE.

Dimethyl sulfone, also known as methylsulfonylmethane, is an organic sulfur-containing compound naturally occurring in various grains and animal tissues [25]. Dimethyl sulfone is an oxidative stress marker [26], and oxidative stress reportedly causes inflammatory responses [27, 28]. In beef cattle, positive correlations existed between the plasma dimethyl sulfone and the FE trait indices [29], which was consistent with the results of the present study (Table 4). In the present study, the fat line was found to have a higher concentration of dimethyl sulfone than the lean line. The reports of other studies and the findings of this study suggested that the fat line had stronger oxidative stress and inflammatory response compared with the lean line and, therefore, a lower FE. As a result, selecting birds with low plasma dimethyl sulfone was suggested for improving the FE.

Epsilon-(gamma-glutamyl)-lysine, belonging to the class of organic compounds known as glutamine and derivatives, was a potential biomarker in the inflammatory reaction [30, 31]. Yang et al. and Xu et al. proposed that high-FE animals have lower inflammatory responses [32, 33]. This study reported that the fat line had a higher concentration of epsilon-(gamma-glutamyl)-lysine than the lean line, demonstrating that the fat line might have a stronger inflammatory reaction than the lean line and, therefore, a lower FE. Therefore, plasma epsilon-(gamma-glutamyl)-lysine was suggested to serve as a potential biomarker for selecting FE traits.

Gamma-glutamyltyrosine is the precursor of tyrosine (Tyr) [34]. Previous studies showed that injecting gamma-glutamyltyrosine tended to increase the plasma Tyr levels in rats [35, 36]. This study reported that the gamma-glutamyltyrosine concentration in the fat line was significantly lower than that in the lean line, suggesting that the fat birds possibly had a weaker capacity for supplying Tyr than the lean birds. Tyr possessed the anti-stress biological function [37]. Therefore, the lean birds might have strong anti-stress ability and, hence, high FE. In addition, Tyr is essential for synthesizing proteins [36], suggesting that the gamma-glutamyltyrosine content can reflect the direction of the protein metabolism in animals. Our previous study found that the breast meat percentage of the fat line was significantly lower than that of the lean line (data not shown). These findings suggested the fat birds possibly had a weaker capacity for supplying Tyr and protein synthesis ability than the lean birds and, therefore, a lower FE.

2-Oxoadipic acid is a key metabolite of tryptophan and lysine [38, 39], identified as a biomarker for acute myocardial infarction and relevant pathological changes in inflammation [40]. In this study, the fat line was found to have a higher concentration of 2-oxoadipic acid than the lean line, suggesting that the fat line might have a stronger inflammatory reaction than the lean line. Therefore, the plasma 2-oxoadipic acid was suggested to serve as the potential biomarker for selecting FE traits.

L-homoarginine is a nonessential cationic amino acid synthesized from arginine and lysine in vivo of bird [41]. The genome-wide association studies identified an association between the L-homoarginine plasma concentration and SNPs related to the AGAT gene [42]. Further metabolic analysis revealed that humans and mice with AGAT deficiency were inefficient in synthesizing L-homoarginine from L-arginine and L-lysine [42, 43]. The AGAT-deficient mice were found to exhibit decreased fat deposition, attenuated gluconeogenesis, reduced cholesterol levels, and enhanced glucose tolerance [44]. Based on these advances about L-homoarginine, fat deposition was found to show a positive correlation with the concentration of L-homoarginine. A study by Ramayo-Caldas et al. reported that lipid metabolism was associated with FE in pigs [23]. In this study, the fat line was found to have a higher concentration of L-homoarginine than the lean line, which was consistent with previous findings. Therefore, plasma L-homoarginine was suggested to serve as the potential biomarker for selecting FE traits.

Testosterone is an anabolic steroid from the androstane class of steroids [45]. It is the most important androgen, stimulating muscle growth and inhibiting lipid synthesis [46, 47]. In humans, testosterone replacement therapy reduces the body fat mass, particularly in truncal adiposity, cholesterol, and triglycerides [48]. In pigs, testosterone reduces the serum cholesterol levels [49], while in broilers, testosterone promotes the proliferation of the embryonic myoblasts and development of the skeletal muscle and inhibits the deposition of the abdominal fat [50]. Zhou et al. found that high muscle yield improved FE [51]. Ramayo-Caldas et al. found that more abdominal fat deposition reduced FE [23]. The concentration of plasma testosterone was found to be significantly lower in the fat line than in the lean line in the present study. A previous study from our lab found that the breast meat percentage was significantly lower in the fat line than in the lean line (data not shown), and the fat line was found to have a higher AFP and lower FE than the lean line [11]. Based on the reports from other studies and the results of the present study, plasma testosterone was suggested to improve FE.

Adenosine 5′-triphosphate decomposition produces adenosine 5′-diphosphate, and adenosine 5′-diphosphate decomposition produces adenosine 5′-monophosphate [52]. Under normal physiological conditions, the extracellular concentrations of adenosine 5′-triphosphate and its products were found to be considerably lower than the intracellular concentrations [53]. However, the extracellular concentrations of adenosine 5′-triphosphate and its products were found to increase markedly under inflammation [54, 55]. The fat line in this study was found to have a higher plasma concentration of adenosine 5'-monophosphate than the lean line. Based on the reports from previous studies and the results of the present study, the fat line was suggested to have a stronger inflammatory reaction and, therefore, a lower FE.

Adrenic acid, an endogenously synthesized polyunsaturated free fatty acid, can induce oxidative stress accompanied by cell death [56]. Many researchers highlighted the negative effects of adrenic acid on organisms. The adrenic acid level increased in patients with nonalcoholic fatty liver disease (NAFLD) and model mice [56]. Higher hepatic and plasma levels of free adrenic acid were found in the mouse model with hepatic steatosis, inflammation, mild fibrosis, obesity, and hypercholesterolemia [57]. However, in the present study, the concentration of plasma adrenic acid was found to be lower in the fat line than in the lean line. This result was inconsistent with the results of previous studies. Considering that the fat birds had higher FE trait indices, adrenic acid was hypothesized to possibly increase the FE. Further research is needed to investigate the mechanism by which adrenic acid negatively correlates with the indices of the FE traits.

Calcitriol, produced by vitamin D metabolism, is an active molecule of vitamin D that exerts its biological activity and functions in a variety of tissues in the body [58, 59]. Calcitriol is central to calcium and phosphate homeostasis and essentially ensures the proper development and maintenance of bone [60]. Besides, calcitriol is involved in maintaining immune homeostasis [61]. Calcitriol is used in humans for the therapeutic applications in immune dysfunction [62]. Studies on humans identified lean individuals to have significantly higher serum calcitriol levels than obese ones [63]. In Angus cattle, the genetic variation in RFI was found to be associated with the immune competence traits, suggesting that high-FE animals possibly had stronger immune competence [64]. This study found that the serum calcitriol concentration was significantly lower in the fat birds than in the lean birds. Combining the results of this study with the reports of calcitriol in both humans and animals, high-FE birds were suggested to have better immune competence.

Among 10 markers, 6 plasma metabolite markers (7-ketocholesterol, dimethyl sulfone, epsilon-(gamma-glutamyl)-lysine, 2-oxoadipic acid, adenosine 5'-monophosphate, and calcitriol) were found to be involved in inflammatory reaction and immune response. The genome-wide association study of FE in the chickens reported that the inflammatory reaction and immune response might affect FE [32]. The transcriptome analysis of the chicken breast muscle [33], duodenum [65], jejunum [66], and liver [32] revealed that both inflammatory reactions and immune response might affect FE. Based on the biological functions of metabolites, the present study found that high-FE chickens had lower inflammation and higher immunity. The immune response and inflammatory studies showed that the energy requirements of animals increased dramatically during inflammation, leading to less energy available for protein deposition and, hence, a lower FE [67]. The present findings combined with the results from previous reports indicated that reducing inflammation and improving immunity could improve the broiler FE.

Conclusions

In summary, this study aimed to identify potential metabolite biomarkers for breeding "food-saving broilers" using the NEAUHLF lines, as an ideal animal model of FE traits for analyzing the relationship between FE traits and metabolites using the metabolomics strategy. According to the pipeline for the selection of the metabolite biomarkers established in this study, it was suggested that 10 metabolites, including 7-ketocholesterol, dimethyl sulfone, epsilon-(gamma-glutamyl)-lysine, gamma-glutamyltyrosine, 2-oxoadipic acid, L-homoarginine, testosterone, adenosine 5'-monophosphate, adrenic acid, and calcitriol could be used as the potential biomarkers for breeding the "food-saving broilers".

Availability of data and materials

The datasets from the current study are available from the corresponding author on reasonable request.

Abbreviations

AFP:

Abdominal fat percentage

AFW:

Abdominal fat weight

BW4:

Body weights at 4 weeks

BW7:

Body weights at 7 weeks

BWG:

Gain in body weights

FCR:

Feed conversion ratio

FE:

Feed efficiency

FI:

Feed intake

HMDB:

Human metabolome database

MMBW:

Metabolic mid-test body weight

NAFLD:

Nonalcoholic fatty liver disease

NEAUHLF:

Northeast Agricultural University broiler lines divergently selected for abdominal fat content

NRC:

National research council

PLS-DA:

Partial least square-discriminant analysis

RFI:

Residual feed intake

VIP:

Variable importance in projection

VLDL:

Very low-density lipoprotein

References

  1. Hayes BJ, Lewin HA, Goddard ME. The future of livestock breeding: genomic selection for efficiency, reduced emissions intensity, and adaptation. Trends Genet. 2013;29:206–14. https://doi.org/10.1016/j.tig.2012.11.009.

    Article  CAS  PubMed  Google Scholar 

  2. Nkrumah JD, Okine EK, Mathison GW, Schmid K, Li C, Basarab JA, et al. Relationships of feedlot feed efficiency, performance, and feeding behavior with metabolic rate, methane production, and energy partitioning in beef cattle. J Anim Sci. 2006;84:145–53. https://doi.org/10.2527/2006.841145x.

    Article  CAS  PubMed  Google Scholar 

  3. Willams OW, Miller SP, Wood BJ. Aspects of selection for feed efficiency in meat producing poultry. Worlds Poult Sci J. 2013;69:77–88. https://doi.org/10.1017/S004393391300007X.

    Article  Google Scholar 

  4. Beauclercq S, Nadal-Desbarats L, Hennequet-Antier C, Gabriel I, Tesseraud S, Calenge F, et al. Relationships between digestive efficiency and metabolomic profiles of serum and intestinal contents in chickens. Sci Rep. 2018;8:6678. https://doi.org/10.1038/s41598-018-24978-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yan W, Sun C, Wen C, Ji C, Zhang D, Yang N. Relationships between feeding behaviors and performance traits in slow-growing yellow broilers. Poult Sci. 2019;98:548–55. https://doi.org/10.3382/ps/pey424.

    Article  PubMed  Google Scholar 

  6. Hayton S, Maker GL, Mullaney I, Trengove RD. Experimental design and reporting standards for metabolomics studies of mammalian cell lines. Cell Mol Life Sci. 2017;74:4421–41. https://doi.org/10.1007/s00018-017-2582-1.

    Article  CAS  PubMed  Google Scholar 

  7. Goldansaz SA, Guo AC, Sajed T, Steele MA, Plastow GS, Wishart DS. Livestock metabolomics and the livestock metabolome: a systematic review. PLoS One. 2017;12:e177675. https://doi.org/10.1371/journal.pone.0177675.

    Article  CAS  Google Scholar 

  8. Novais FJ, Pires P, Alexandre PA, Dromms RA, Iglesias AH, Ferraz J, et al. Identification of a metabolomic signature associated with feed efficiency in beef cattle. BMC Genomics. 2019;20:8. https://doi.org/10.1186/s12864-018-5406-2.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Carmelo V, Banerjee P, Da SDW, Kadarmideen HN. Metabolomic networks and pathways associated with feed efficiency and related-traits in Duroc and Landrace pigs. Sci Rep. 2020;10:255. https://doi.org/10.1038/s41598-019-57182-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Guo L, Sun B, Shang Z, Leng L, Wang Y, Wang N, et al. Comparison of adipose tissue cellularity in chicken lines divergently selected for fatness. Poult Sci. 2011;90:2024–34. https://doi.org/10.3382/ps.2010-00863.

    Article  CAS  PubMed  Google Scholar 

  11. Chen C, Su Z, Li Y, Luan P, Wang S, Zhang H, et al. Estimation of the genetic parameters of traits relevant to feed efficiency: result from broiler lines divergent for high or low abdominal fat content. Poult Sci. 2021;100:461–6. https://doi.org/10.1016/j.psj.2020.10.028.

    Article  PubMed  Google Scholar 

  12. National Research Council (NRC). Nutrient requirements of poultry. 9th ed. Washington, DC: National Academy Press; 1994.

    Google Scholar 

  13. Brink LR, Mercer KE, Piccolo BD, Chintapalli SV, Elolimy A, Bowlin AK, et al. Neonatal diet alters fecal microbiota and metabolome profiles at different ages in infants fed breast milk or formula. Am J Clin Nutr. 2020;111:1190–202. https://doi.org/10.1093/ajcn/nqaa076.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Reisz JA, Zheng C, D’Alessandro A, Nemkov T. Untargeted and semi-targeted lipid analysis of biological samples using mass spectrometry-based metabolomics. Methods Mol Biol. 2019;1978:121–35. https://doi.org/10.1007/978-1-4939-9236-2_8.

    Article  CAS  PubMed  Google Scholar 

  15. Consolo N, Silva JD, Buarque V, Higuera-Padilla A, Barbosa L, Zawadzki A, et al. Selection for growth and precocity alters muscle metabolism in nellore cattle. Metabolites. 2020;10:58. https://doi.org/10.3390/metabo10020058.

    Article  CAS  PubMed Central  Google Scholar 

  16. Buchanan JW, Reecy JM, Garrick DJ, Duan Q, Beitz DC, Mateescu RG. Genetic parameters and genetic correlations among triacylglycerol and phospholipid fractions in Angus cattle. J Anim Sci. 2015;93:522–8. https://doi.org/10.2527/jas.2014-8418.

    Article  CAS  PubMed  Google Scholar 

  17. Dong JQ, Zhang H, Jiang XF, Wang SZ, Du ZQ, Wang ZP, et al. Comparison of serum biochemical parameters between two broiler chicken lines divergently selected for abdominal fat content. J Anim Sci. 2015;93:3278–86. https://doi.org/10.2527/jas.2015-8871.

    Article  CAS  PubMed  Google Scholar 

  18. Jasbi P, Wang D, Cheng SL, Fei Q, Cui JY, Liu L, et al. Breast cancer detection using targeted plasma metabolomics. J Chromatogr B Analyt Technol Biomed Life Sci. 2019;1105:26–37. https://doi.org/10.1016/j.jchromb.2018.11.029.

    Article  CAS  PubMed  Google Scholar 

  19. Yan Q, Liew Z, Uppal K, Cui X, Ling C, Heck JE, et al. Maternal serum metabolome and traffic-related air pollution exposure in pregnancy. Environ Int. 2019;130:104872. https://doi.org/10.1016/j.envint.2019.05.066.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Xu B, Chen Y, Chen X, Gan L, Zhang Y, Feng J, et al. Metabolomics profiling discriminates prostate cancer from benign prostatic hyperplasia within the prostate-specific antigen gray zone. Front Oncol. 2021;11:730638. https://doi.org/10.3389/fonc.2021.730638.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Pariente A, Pelaez R, Perez-Sala A, Larrayoz IM. Inflammatory and cell death mechanisms induced by 7-ketocholesterol in the retina. Implications for age-related macular degeneration. Exp Eye Res. 2019;187:0107746. https://doi.org/10.1016/j.exer.2019.107746.

    Article  CAS  Google Scholar 

  22. Chang MC, Chen YJ, Liou EJ, Tseng WY, Chan CP, Lin HJ, et al. 7-Ketocholesterol induces ATM/ATR, Chk1/Chk2, PI3K/Akt signalings, cytotoxicity and IL-8 production in endothelial cells. Oncotarget. 2016;7:74473–83. https://doi.org/10.18632/oncotarget.12578.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ramayo-Caldas Y, Ballester M, Sanchez JP, Gonzalez-Rodriguez O, Revilla M, Reyer H, et al. Integrative approach using liver and duodenum RNA-Seq data identifies candidate genes and pathways associated with feed efficiency in pigs. Sci Rep. 2018;8:558. https://doi.org/10.1038/s41598-017-19072-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Alexandre PA, Kogelman LJ, Santana MH, Passarelli D, Pulz LH, Fantinato-Neto P, et al. Liver transcriptomic networks reveal main biological processes associated with feed efficiency in beef cattle. BMC Genomics. 2015;16:1073. https://doi.org/10.1186/s12864-015-2292-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kim YH, Kim DH, Lim H, Baek DY, Shin HK, Kim JK. The anti-inflammatory effects of methylsulfonylmethane on lipopolysaccharide-induced inflammatory responses in murine macrophages. Biol Pharm Bull. 2009;32:651–6. https://doi.org/10.1248/bpb.32.651.

    Article  CAS  PubMed  Google Scholar 

  26. Ehrich J, Dubourg L, Hansson S, Pape L, Steinle T, Fruth J, et al. Serum myo-inositol, dimethyl sulfone, and valine in combination with creatinine allow accurate assessment of renal insufficiency-A proof of concept. Diagnostics (Basel). 2021;11:234. https://doi.org/10.3390/diagnostics11020234.

    Article  CAS  Google Scholar 

  27. Liu J, Wang S, Zhang Q, Li X, Xu S. Selenomethionine alleviates LPS-induced chicken myocardial inflammation by regulating the miR-128-3p-p38 MAPK axis and oxidative stress. Metallomics. 2020;12:54–64. https://doi.org/10.1039/c9mt00216b.

    Article  CAS  PubMed  Google Scholar 

  28. Wang B, Cui Y, Zhang Q, Wang S, Xu S. Selenomethionine alleviates LPS-induced JNK/NLRP3 inflammasome-dependent necroptosis by modulating miR-15a and oxidative stress in chicken lungs. Metallomics. 2021;13:048. https://doi.org/10.1093/mtomcs/mfab048.

    Article  Google Scholar 

  29. Li J, Mukiibi R, Wang Y, Plastow GS, Li C. Identification of candidate genes and enriched biological functions for feed efficiency traits by integrating plasma metabolites and imputed whole genome sequence variants in beef cattle. BMC Genomics. 2021;22:823. https://doi.org/10.1186/s12864-021-08064-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bowness JM, Tarr AH. Increase in transglutaminase and its extracellular products in response to an inflammatory stimulus by lipopolysaccharide. Mol Cell Biochem. 1997;169:157–63. https://doi.org/10.1023/a:1006846400478.

    Article  CAS  PubMed  Google Scholar 

  31. Skill NJ, Griffin M, El NA, Sanai T, Haylor JL, Fisher M, et al. Increases in renal epsilon-(gamma-glutamyl)-lysine crosslinks result from compartment-specific changes in tissue transglutaminase in early experimental diabetic nephropathy: pathologic implications. Lab Invest. 2001;81:705–16. https://doi.org/10.1038/labinvest.3780279.

    Article  CAS  PubMed  Google Scholar 

  32. Xu Z, Ji C, Zhang Y, Zhang Z, Nie Q, Xu J, et al. Combination analysis of genome-wide association and transcriptome sequencing of residual feed intake in quality chickens. BMC Genomics. 2016;17:594. https://doi.org/10.1186/s12864-016-2861-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Yang L, He T, Xiong F, Chen X, Fan X, Jin S, et al. Identification of key genes and pathways associated with feed efficiency of native chickens based on transcriptome data via bioinformatics analysis. BMC Genomics. 2020;21:292. https://doi.org/10.1186/s12864-020-6713-y.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Price ND, Magis AT, Earls JC, Glusman G, Levy R, Lausted C, et al. A wellness study of 108 individuals using personal, dense, dynamic data clouds. Nat Biotechnol. 2017;35:747–56. https://doi.org/10.1038/nbt.3870.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Berger DC, Hilton MA, Hilton FK, Duncan SD, Radmacher PG, Greene SM. Intravenous gamma-glutamyl-tyrosine elevates brain tyrosine but not catecholamine concentrations in normal rats. Metabolism. 1996;45:126–32. https://doi.org/10.1016/s0026-0495(96)90210-3.

    Article  CAS  PubMed  Google Scholar 

  36. Radmacher PG, Hilton MA, Hilton FK, Duncan SD, Adamkin DH. Use of the soluble peptide gamma-L-glutamyl-L-tyrosine to provide tyrosine in total parenteral nutrition in rats. JPEN J Parenter Enteral Nutr. 1993;17:337–44. https://doi.org/10.1177/0148607193017004337.

    Article  CAS  PubMed  Google Scholar 

  37. Kabuki Y, Mizobe Y, Yamada S, Furuse M. Dietary L-tyrosine alleviates the behavioral alterations induced by social isolation stress in mice. Brain Res Bull. 2009;80:389–96. https://doi.org/10.1016/j.brainresbull.2009.08.014.

    Article  CAS  PubMed  Google Scholar 

  38. Shibata K, Yasui M, Sano M, Fukuwatari T. Fluorometric determination of 2-oxoadipic acid, a common metabolite of tryptophan and lysine, by high-performance liquid chromatography with pre-chemical derivatization. Biosci Biotechnol Biochem. 2011;75:185–7. https://doi.org/10.1271/bbb.100723.

    Article  CAS  PubMed  Google Scholar 

  39. Hiratsuka C, Fukuwatari T, Shibata K. Fate of dietary tryptophan in young Japanese women. Int J Tryptophan Res. 2012;5:33–47. https://doi.org/10.4137/IJTR.S10497.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Jiang P, Dai W, Yan S, Chen Z, Xu R, Ding J, et al. Biomarkers in the early period of acute myocardial infarction in rat serum and protective effects of Shexiang Baoxin Pill using a metabolomic method. J Ethnopharmacol. 2011;138:530–6. https://doi.org/10.1016/j.jep.2011.09.049.

    Article  PubMed  Google Scholar 

  41. Angkanaporn K, Ravindran V, Bryden WL. De novo synthesis of homoarginine in chickens is influenced by dietary concentrations of lysine and arginine. Nutr Res. 1997;17:99–110. https://doi.org/10.1016/s0271-5317(96)00236-9.

    Article  CAS  Google Scholar 

  42. Choe CU, Atzler D, Wild PS, Carter AM, Boger RH, Ojeda F, et al. Homoarginine levels are regulated by L-arginine:glycine amidinotransferase and affect stroke outcome: results from human and murine studies. Circulation. 2013;128:1451–61. https://doi.org/10.1161/CIRCULATIONAHA.112.000580.

    Article  CAS  PubMed  Google Scholar 

  43. Davids M, Ndika JD, Salomons GS, Blom HJ, Teerlink T. Promiscuous activity of arginine:glycine amidinotransferase is responsible for the synthesis of the novel cardiovascular risk factor homoarginine. FEBS Lett. 2012;586:3653–7. https://doi.org/10.1016/j.febslet.2012.08.020.

    Article  CAS  PubMed  Google Scholar 

  44. Choe CU, Nabuurs C, Stockebrand MC, Neu A, Nunes P, Morellini F, et al. L-arginine:glycine amidinotransferase deficiency protects from metabolic syndrome. Hum Mol Genet. 2013;22:110–23. https://doi.org/10.1093/hmg/dds407.

    Article  CAS  PubMed  Google Scholar 

  45. Ye L, Su ZJ, Ge RS. Inhibitors of testosterone biosynthetic and metabolic activation enzymes. Molecules. 2011;16:9983–10001. https://doi.org/10.3390/molecules16129983.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Li D, Wang Q, Shi K, Lu Y, Yu D, Shi X, et al. Testosterone promotes the proliferation of chicken embryonic myoblasts via androgen receptor mediated PI3K/Akt signaling pathway. Int J Mol Sci. 2020;21:1152. https://doi.org/10.3390/ijms21031152.

    Article  CAS  PubMed Central  Google Scholar 

  47. Mudali S, Dobs AS. Effects of testosterone on body composition of the aging male. Mech Ageing Dev. 2004;125:297–304. https://doi.org/10.1016/j.mad.2004.01.004.

    Article  CAS  PubMed  Google Scholar 

  48. Kelly DM, Jones TH. Testosterone: a metabolic hormone in health and disease. J Endocrinol. 2013;217:25–45. https://doi.org/10.1530/JOE-12-0455.

    Article  CAS  Google Scholar 

  49. Cai Z, Xi H, Pan Y, Jiang X, Chen L, Cai Y, et al. Effect of testosterone deficiency on cholesterol metabolism in pigs fed a high-fat and high-cholesterol diet. Lipids Health Dis. 2015;14:18. https://doi.org/10.1186/s12944-015-0014-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Cui X, Cui H, Liu L, Zhao G, Liu R, Li Q, et al. Decreased testosterone levels after caponization leads to abdominal fat deposition in chickens. BMC Genomics. 2018;19:344. https://doi.org/10.1186/s12864-018-4737-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhou N, Lee WR, Abasht B. Messenger RNA sequencing and pathway analysis provide novel insights into the biological basis of chickens’ feed efficiency. BMC Genomics. 2015;16:195. https://doi.org/10.1186/s12864-015-1364-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Bours MJ, Swennen EL, Di Virgilio F, Cronstein BN, Dagnelie PC. Adenosine 5’-triphosphate and adenosine as endogenous signaling molecules in immunity and inflammation. Pharmacol Ther. 2006;112:358–404. https://doi.org/10.1016/j.pharmthera.2005.04.013.

    Article  CAS  PubMed  Google Scholar 

  53. Yoneyama Y, Suzuki S, Sawa R, Araki T. Plasma adenosine concentrations increase in women with hyperemesis gravidarum. Clin Chim Acta. 2004;342:99–103. https://doi.org/10.1016/j.cccn.2003.12.004.

    Article  CAS  PubMed  Google Scholar 

  54. Latini S, Pedata F. Adenosine in the central nervous system: release mechanisms and extracellular concentrations. J Neurochem. 2001;79:463–84. https://doi.org/10.1046/j.1471-4159.2001.00607.x.

    Article  CAS  PubMed  Google Scholar 

  55. Lazarowski ER, Boucher RC, Harden TK. Mechanisms of release of nucleotides and integration of their action as P2X- and P2Y-receptor activating molecules. Mol Pharmacol. 2003;64:785–95. https://doi.org/10.1124/mol.64.4.785.

    Article  CAS  PubMed  Google Scholar 

  56. Zhao J, Nishiumi S, Tagawa R, Yano Y, Inoue J, Hoshi N, et al. Adrenic acid induces oxidative stress in hepatocytes. Biochem Biophys Res Commun. 2020;532:620–5. https://doi.org/10.1016/j.bbrc.2020.08.102.

    Article  CAS  PubMed  Google Scholar 

  57. Horas HNS, Nishiumi S, Kawano Y, Kobayashi T, Yoshida M, Azuma T. Adrenic acid as an inflammation enhancer in non-alcoholic fatty liver disease. Arch Biochem Biophys. 2017;623–624:64–75. https://doi.org/10.1016/j.abb.2017.04.009.

    Article  CAS  Google Scholar 

  58. Hansdottir S, Monick MM, Hinde SL, Lovan N, Look DC, Hunninghake GW. Respiratory epithelial cells convert inactive vitamin D to its active form: potential effects on host defense. J Immunol. 2008;181:7090–9. https://doi.org/10.4049/jimmunol.181.10.7090.

    Article  CAS  PubMed  Google Scholar 

  59. Rigby WF, Stacy T, Fanger MW. Inhibition of T lymphocyte mitogenesis by 1,25-dihydroxyvitamin D3 (calcitriol). J Clin Invest. 1984;74:1451–5. https://doi.org/10.1172/JCI111557.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Panichi V, Migliori M, Taccola D, Consani C, Giovannini L. Effects of calcitriol on the immune system: new possibilities in the treatment of glomerulonephritis. Clin Exp Pharmacol Physiol. 2003;30:807–11. https://doi.org/10.1046/j.1440-1681.2003.03919.x.

    Article  CAS  PubMed  Google Scholar 

  61. Jones G, Strugnell SA, DeLuca HF. Current understanding of the molecular actions of vitamin D. Physiol Rev. 1998;78:1193–231. https://doi.org/10.1152/physrev.1998.78.4.1193.

    Article  CAS  PubMed  Google Scholar 

  62. Bouillon R, Garmyn M, Verstuyf A, Segaert S, Casteels K, Mathieu C. Paracrine role for calcitriol in the immune system and skin creates new therapeutic possibilities for vitamin D analogs. Eur J Endocrinol. 1995;133:7–16. https://doi.org/10.1530/eje.0.1330007.

    Article  CAS  PubMed  Google Scholar 

  63. Lu L, Yu Z, Pan A, Hu FB, Franco OH, Li H, et al. Plasma 25-hydroxyvitamin D concentration and metabolic syndrome among middle-aged and elderly Chinese individuals. Diabetes Care. 2009;32:1278–83. https://doi.org/10.2337/dc09-0209.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Herd RM, Velazco JI, Smith H, Arthur PF, Hine B, Oddy H, et al. Genetic variation in residual feed intake is associated with body composition, behavior, rumen, heat production, hematology, and immune competence traits in Angus cattle1. J Anim Sci. 2019;97:2202–19. https://doi.org/10.1093/jas/skz077.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Liu R, Liu J, Zhao G, Li W, Zheng M, Wang J, et al. Relevance of the intestinal health-related pathways to broiler residual feed intake revealed by duodenal transcriptome profiling. Poult Sci. 2019;98:1102–10. https://doi.org/10.3382/ps/pey506.

    Article  CAS  PubMed  Google Scholar 

  66. Sinpru P, Riou C, Kubota S, Poompramun C, Molee W, Molee A. Jejunal transcriptomic profiling for differences in feed conversion ratio in slow-growing chickens. Animals (Basel). 2021;11:2606. https://doi.org/10.3390/ani11092606.

    Article  Google Scholar 

  67. Fonseca LD, Eler JP, Pereira MA, Rosa AF, Alexandre PA, Moncau CT, et al. Liver proteomics unravel the metabolic pathways related to feed efficiency in beef cattle. Sci Rep. 2019;9:5364. https://doi.org/10.1038/s41598-019-41813-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors gratefully acknowledge the members of the Poultry Breeding group at Northeast Agricultural University for help in managing birds and collecting data.

Funding

This research was supported by the earmarked fund for CARS-41, National Natural Science Foundation (32102537), White Feather Broiler Breeding Joint Project of the Ministry of Agriculture and Rural Affairs in China (19190526), and “Academic Backbone” Project of Northeast Agricultural University (20XG30).

Author information

Author notes

  1. Zhiyong Su and Xue Bai contributed equally to this work.

Authors and Affiliations

  1. Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin, 150030, China

    Zhiyong Su, Xue Bai, Haoyu Wang, Shouzhi Wang, Chong Chen, Li Leng & Hui Li

  2. Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin, 150030, China

    Zhiyong Su, Xue Bai, Haoyu Wang, Shouzhi Wang, Chong Chen, Li Leng & Hui Li

  3. College of Animal Science and Technology, Northeast Agricultural University, Harbin, 150030, China

    Zhiyong Su, Xue Bai, Haoyu Wang, Shouzhi Wang, Chong Chen, Li Leng & Hui Li

  4. Fujian Sunnzer Biotechnology Development Co., Ltd, Guangze, 354100, China

    Fan Xiao, Huaishun Guo & Haihe Gao

Authors
  1. Zhiyong Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xue Bai
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Haoyu Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shouzhi Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Chong Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Fan Xiao
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Huaishun Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Haihe Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Li Leng
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Hui Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization: HL; methodology: HL, XB, LL, and ZS; software: ZS, XB, and HW; formal analysis: ZS, XB, HL, and SW; resources: HL, LL, CC, and ZS; investigation: FX, HG, and HG; funding acquisition: HL, LL, and XB; writing-original draft preparation: ZS, XB, and HL. All authors have read and approved the final manuscript. ZS and XB equally contributed to this work.

Corresponding author

Correspondence to Hui Li.

Ethics declarations

Ethics approval and consent to participate

All the animal work was conducted according to the guidelines for the Care and Use of Experimental Animals established by the Ministry of Science and Technology of the People’s Republic of China (Approval No. 2006–398) and approved by the Laboratory Animal Management Committee of Northeast Agricultural University.

Consent for publication

All authors have read and agreed to the published version of the manuscript.

Competing interests

The authors declare no competing financial interest.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Su, Z., Bai, X., Wang, H. et al. Identification of biomarkers associated with the feed efficiency by metabolomics profiling: results from the broiler lines divergent for high or low abdominal fat content. J Animal Sci Biotechnol 13, 122 (2022). https://doi.org/10.1186/s40104-022-00775-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-022-00775-3

Keywords

Journal of Animal Science and Biotechnology

ISSN: 2049-1891