Introduction
Low-protein diets have received significant attention in the broiler industry because of their significant role in reducing nitrogen emissions [1,2,3]. Low-protein diets were formulated to meet amino acid nutritional requirements by incorporating crystalline amino acids [4]; however, long-term feeding of low-protein diets in broilers results in an increase in fat pad rate (FR) and feed conversion ratio (FCR) [5,6,7], nevertheless. Several studies have indicated that broilers fed with varying patterns of amino acids exhibit notable enhancements in body weight (BW) and FCR, as well as significant alterations in carcass characteristics [8,9,10]. Furthermore, Macelline et al. [11] conducted a study on broilers fed diets with two different protein contents and amino acid patterns and revealed that the optimal growth of broilers necessitates distinct amino acid patterns depending on the protein content. Consequently, it proved challenging to meet bird growth requirements by formulating a low-protein diet based solely on the amino acid patterns established in a normal-protein diet.
During protein synthesis, animals require the combined participation of essential amino acids (EAAs) and non-essential amino acids (NEAAs). Numerous studies have demonstrated that animals fed low-protein diets exhibit poor growth performance, accompanied by a reduction in blood NEAA concentration [12,13,14], indicating that NEAAs or total nitrogen deficiency may limit animal growth [15]. As the total protein content decreases in a low-protein diet, the body’s supply of nitrogen used for synthesizing NEAAs also decreases, resulting in the conversion of some dietary EAAs into NEAAs through oxidative decomposition. Therefore, it is imperative to reassess EAA concentration in low-protein diets to determine the optimal EAA pattern.
In current broiler feeding experiments, the main methods for evaluating dietary amino acids are addition and deletion assays, each with its advantages and disadvantages [16]. Combining these two methods effectively allows for the determination of increasing and decreasing trends in amino acid concentrations in the diet by comparing differences in animal traits when amino acids are present or absent. Furthermore, it is important to note the interactions among amino acids of the same type; therefore, any change in the concentration of one amino acid will have a certain impact on another [17, 18]. Therefore, EAAs can be categorized into different combinations based on their structural and functional properties. The pattern of each amino acid combination in the diet was adjusted using a combination of addition and deletion methods to evaluate the changes in amino acid composition.
The indices commonly used for assessing dietary amino acid balance include the nitrogen balance assay and serum urea nitrogen or uric acid concentration [19]. These indices are based on the excretion or deposition of total nitrogen in the diet. However, the nitrogen balance assay is a cumbersome procedure with significant variation, hindering its large-scale application in production. As mentioned earlier, dietary total nitrogen is involved in NEAA synthesis, indicating that serum urea nitrogen or uric acid concentrations do not fully reflect the extent of amino acid balance. Therefore, it is necessary to identify biomarkers that accurately reflect amino acid balance. Metabolomics serves as a systematic biological analysis tool for studying amino acid nutrition [20]; it intuitively reflects metabolic differences among animals under different conditions and assesses how metabolic balance is affected by amino acid deficiency or excess [21, 22]. Subsequently, diets with optimal and suboptimal amino acid balances were selected, followed by untargeted metabolomic analysis to identify biomarkers for assessing amino acid balance.
Materials and methods
Animal ethics
This study was conducted at the Zhuozhou Experimental Base of China Agricultural University and approved by the Animal Welfare and Ethics Review Committee of China Agricultural University (AW61504202-1-2).
Experimental diets
The feeding standard for the experimental diet was based on the Feeding Standard of Chickens (2004) [23], with equal energy and electrolyte balance. Starter, grower, and finisher phases had a protein content of 23% (0–14 d), 21% (15–28 d), and 20% (29–42 d), respectively, in the positive control normal diet. Further, the negative control had a protein content of 21%, 19%, and 18% in the starter, grower, and finisher phases, respectively, and the low-protein diet, dietary composition, and nutrient concentration are shown in Table 1. Dietary EAAs were divided into five combinations according to their structure and properties: lysine and arginine, methionine, cysteine, branched-chain amino acids (BCAA), threonine, and tryptophan. According to the principles of addition and deletion methods, the five amino acid combinations were adjusted by 10% based on a low-protein diet (Table 2). The specific treatments were: the normal-protein control diet [Con(+)], the low-protein control diet [Con(-)], the lysine and arginine restriction diet [L&A(-)], the lysine and arginine supplementation diet [L&A(+)], the methionine and cysteine restriction diet [M&C(-)], the methionine and cysteine supplementation diet [M&C(+)], the BCAA restriction diet [BCAA(-)], the BCAA supplementation diet [BCAA(+)], the threonine restriction diet [Thr(-)], the threonine supplementation diet [Thr(+)], the tryptophan restriction diet [Trp(-)], and the tryptophan supplementation diet [Trp(+)], detailed dietary composition and nutrient concentration are shown in Additional file 1: Table S1–S3. All feeds were mixed evenly in a vertical spiral mixer (DLH-30, Guibao) and granulated using steam at 80 °C.
All the samples of experimental diets were polished into powder and stored at −20 °C to determine the content of crude protein, energy, and EAA concentrations.
Animals and feeding management
A total of 720 1-day-old healthy AA broilers, half male, and half female (identified by feather speed) with an initial BW of 43.5 ± 0.5 g, were randomly assigned to 12 diets according to sex with 6 pens per treatment and 10 broilers per pen. The experiment lasted for 40 d and was divided into 3 phases: the starter phase (0–14 d), grower phase (15–28 d), and finisher phase (29–40 d). The broilers were kept captive, with free access to feed and water. The illumination period of the poultry houses was gradually reduced from 24 h/d to 20 h/d. The temperature was gradually lowered from 35 to 22 °C. The immunization procedure was carried out according to the management manual, and the environment of the poultry house and health status of the broilers were checked daily.
Sample collection and analysis
All animal experiments were performed using the same procedures. The BW of each pen was measured, feed consumption was recorded at the beginning and end of each phase, and the number of deaths and the BW of dead broilers were recorded daily to correct the growth performance data. At 40 d, 4 broilers close to the average BW were selected from each pen, and 2 were used for slaughter. The slaughtering process was in accordance with the Operating Procedure of Livestock and Poultry Slaughtering—Chicken [24], and the measurements included dressing, half-eviscerated, eviscerated, fat pad, breast, and thigh rates. The remaining 2 were used for sample collection. After the broilers were euthanized by electric shock, blood was quickly collected from the jugular vein and placed in a coagulant tube. After standing at room temperature for 4 h, the blood was transferred to a 4 °C centrifuge, centrifuged at 3,000 r/min for 15 min, and the supernatant was stored at −80 °C. Liver and chest muscle tissues were collected and placed in a sterile cryopreservation tube, and after quenching with liquid nitrogen, they were transferred to −80 °C for preservation.
Serum biochemical parameter analysis
Serum concentrations of total protein (TP), albumin (ALB), globulin (GLB), triglyceride (TG), total cholesterol (TC), blood urea nitrogen (BUN), uric acid (UA), glucose (GLU), alanine transaminase (ALT), and aspartate aminotransferase (AST), were analyzed using an automatic biochemical analyzer (BK1200, Biobase, Jinan, Shandong, China). All the operating procedures were performed according to the manufacturer’s instructions.
Hepatic enzyme activities analysis
The liver tissue was mixed and ground with normal saline, and the supernatant was removed after centrifugation and transferred to 4 °C for storage. Aspartate aminotransferase (AST), alanine aminotransferase (ALT), acetyl-CoA carboxylase (ACC), fatty acid synthase (FAS), glucose-6-phosphatase (G6PC), pyruvate carboxylase (PC), HMG-CoA reductase (HMGR), and xanthine oxidase (XOD) were measured using commercial reagent kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). All procedures were performed according to the manufacturer’s instructions.
Real-time quantitative PCR
RNA extraction, reverse transcription, and qRT-PCR were performed as described by Yang et al. [25]. Total RNA was extracted from the breast muscle tissue using RNAiso Plus reagent (Takara, Beijing, China). After the purity and concentration of RNA were determined using a nucleic acid analyzer, cDNA was synthesized using a reverse transcription kit, diluted, and analyzed using qRT-PCR. The resulting mRNA expression was normalized to GAPDH, and the remaining relative mRNA expression was calculated. Primers for GAPDH, 4EBP1, S6K1, MAFbx, and MuRF1 were designed and commercially synthesized according to sequences in GenBank (Table 3) (Sangon Biotechnology, Shanghai, China).
Metabolite extraction and detection
At low temperature, 25 mg sample was weighed into an EP tube, and 500 µL extraction solution (methanol:acetonitrile:water=2:2:1; v/v) containing isotope labeled internal standard mixture was added. The mixture was vortexed and stirred for 30 s. The homogenate beads were added and homogenized in a homogenizer (35 Hz, 4 min), transferred to an ice water bath, sonicated for 5 min, and the process was repeated 3 times. The homogenate was left at −40 °C for 1 h before being transferred to a centrifuge, and centrifuged at 12,000 r/min for 15 min at 4 °C. The supernatant was collected for further examination.
Chromatographic separation of the target compounds was performed using an ultra-high performance liquid chromatograph with a chromatographic column (Waters ACQUITY UPLC BEH Amide; 2.1 mm × 50 mm, 1.7 μm). Liquid chromatography phase A was aqueous (25 mmol/L ammonium acetate and 25 mmol/L aqueous ammonia), and phase B was acetonitrile. Sample plate temperature: 4 °C. The injection volume was 2 µL. Primary and secondary mass spectrometry data were acquired using Orbitrap Exploris 120 mass spectrometer (Xcalibur, version: 4.4, Thermo).
Data preprocessing and metabolite identification
The raw metabolome data were converted into mzXML format by ProteoWizard software, and then metabolites were visualized by R language.
The identification method of metabolites is as described by Marín-García et al. [26]. Local self-built data databases and public databases were used to search databases [27]. By matching the retention time, molecular mass (molecular mass error within < 10 ppm), secondary fragmentation spectrum, collision energy and other information of metabolites in the database, the structure of metabolites in biological samples was identified, and the identification results were strictly manually checked and confirmed. The identification Level was Level 2 and above.
Statistical analysis
Each pen was considered as an experimental unit. One-way ANOVA was used for data analysis using SPSS 27 (IBM), and the results are presented as the means. P < 0.05 was considered statistically significant difference, 0.05 < P < 0.1 indicated a trend of difference. To clarify the adjustment trends of different EAA combinations, growth performance and carcass characteristics were analyzed for the same amino acid combinations. Hepatic enzyme activity and serum biochemical parameters were graphically plotted using Origin 2022 software (OriginLab).
Principal component analysis (PCA) was used to outline the population structure, and orthogonal partial least squares discriminant analysis (OPLS-DA) was used to analyze differences between treatments. The goodness of fit (R2) and goodness of prediction (Q2) were used to represent the cross-validation results, and the closer they were to 1, the more reliable the model was. The t-test and volcano plot were used for data analysis, and metabolites with VIP > 1, fold change > 2 or < 0.5 and P < 0.05 were considered differential metabolites. A Venn plot was used to show the shared differential metabolites among the four treatments, and their potential as markers was determined based on the receiver operating characteristic (ROC) curve, with an area under the curve (AUC) closer to 1 indicating greater potential. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gallus gallus (chicken) metabolome database [28] were used to enrich metabolic pathways related to differential metabolites, with higher percentages of metabolites and rich factors indicating more significant enrichment of the pathway.
Results
Growth performance
Lysine and arginine allowance
In male broilers, significant differences across the four treatment groups were observed in BW, average daily feed intake (ADFI), average daily gain (ADG), and FCR throughout the various growth phases. The lysine- and arginine-supplemented groups demonstrated the best overall performance. Specifically, this group had the highest BW during the grower and finisher phases, albeit slightly lower than that of the low-protein control group during the starter phase; ADFI was lower than that of the low-protein control group in the early phases but peaked in the finisher phase. The Lysine and Arginine supplemented groups consistently had the highest ADG, except for the starter phase. Conversely, the normal protein control group exhibited the lowest FCR except during the finishing phase.
For female broilers, the patterns of significant differences in BW, ADFI, ADG, and FCR among the treatments were similar, although the differences in the grower phase were less pronounced than those in males, yet discernible trends were present. Throughout all phases, the Lysine and Arginine supplemented groups outperformed the other treatment groups in terms of BW, ADFI, and ADG. The normal-protein control group consistently had the lowest FCR, except in the finisher phase (Table 4).
Methionine and cysteine allowance
During the starter phase, significant differences were noted in ADFI among the four treatments in male broilers. Similarly, in the grower phase, a notable divergence was observed in the FCR. Trends in the BW and ADG were discernible during the starter phase. However, no significant variations were detected in the remaining indicators. The normal protein control group had the lowest BW, ADFI, and ADG in the starter phase, and the lowest FCR in the grower phase. Female broilers showed the same performance as male broilers, but the methionine and cysteine restriction groups had the lowest FCR in the starter phase and the highest ADFI in the finisher phase (Table 5).
BCAA allowance
There were significant differences in BW, ADFI, ADG, and FCR among the four treatments at different phases in male and female broilers. However, no significant differences were observed between the BCAA supplementation and restriction groups at each phase. Broilers fed diets with adjusted BCAA concentrations exhibited poorer performance than those in the control group (Table 6).
Threonine allowance
There were significant differences in BW, ADFI, ADG, and FCR among the four treatments at different phases in male and female broilers, with the threonine supplementation group showing the worst performance. In addition, broilers that were fed diets with adjusted threonine concentrations exhibited poorer performance than those in the control groups but the threonine restriction group performed slightly better than the threonine supplementation group (Table 7).
Tryptophan allowance
There were significant differences among the four treatments in all indices at the starter and grower phases and in BW at the finisher phase in male broilers. The tryptophan-restricted group showed the worst performance during the entire lifespan. There was no significant difference between the tryptophan-supplemented and low-protein control groups; however, the tryptophan-supplemented group showed an upward trend during the finisher phase. There were differences among the four treatments in FCR at different phases, and in BW and ADFI at the starter phase in female broilers. Except for the normal protein control group, the tryptophan-supplemented group had the lowest FCR in all phases (Table 8).
In general, the growth performance of male and female broilers was similar; however, the differences among females were less significant.
Carcass characteristic
Male and female broilers had the same carcass characteristics and are therefore described uniformly in the following paragraphs.
Lysine and arginine allowance
There were significant differences in all indices except the thigh rate (TR) among the four treatments. The lysine- and arginine-restricted groups had the lowest dressing rate (DR), half-evisceration rate (HR), evisceration rate (ER), and highest fat pad rate (FR). However, the lysine and arginine supplementation group was significantly higher than the low-protein control group in breast rate (BR) (Tables 9 and 10).
Methionine and cysteine allowance
There were significant differences in the FR and BR among the four treatments. FR was higher in the methionine and cysteine restriction groups and the methionine and cysteine supplementation groups than in the normal-protein control and low-protein control groups, whereas BR was not different from the low-protein control group (Tables 9 and 10).
BCAA allowance
There were significant differences in ER, FR, and BR among the four treatments. Both the BCAA restriction group and the BCAA supplementation group were lower than the low-protein control and the normal-protein control groups in ER and higher in FR, but there was no difference in BR compared to the low-protein control group (Tables 9 and 10).
Threonine allowance
There were significant differences in FR and BR among the four treatments, with both the threonine-restriction and threonine-supplementation groups having significantly higher FR and lower BR values than the other 2 groups (Tables 9 and 10).
Tryptophan allowance
There were significant differences among the four treatments in ER, FR, and BR, and a trend in TR. Both the tryptophan restriction group and the tryptophan supplementation group had a lower ER than the normal-protein control group and the low-protein control group, and there was a significant difference in the TR between them. In addition, the tryptophan restriction group had the highest FR (Tables 9 and 10).
Hepatic enzyme activities and serum biochemical parameters
Hepatic enzyme activities and serum biochemical parameters of broilers were determined according to amino acid metabolic pathways (Fig. 1).

Amino acid metabolism pathways. AST: aspartate aminotransferase; ALT: alanine aminotransferase; TG: Triglyceride; ACC: Acetyl CoA carboxylase; FAS: Fatty acid synthase; GLU: Blood glucose; G6PC: Glucose-6-phosphatase; PC: Pyruvate carboxylase; TC: Cholesterol; HMGR: HMG CoA reductase; XOD: Xanthine oxidase; BUN: Urea nitrogen; UA: Uric acid; S6K1: The 70-kDa ribosomal protein S6 kinase; 4EBP1: Eukaryotic translation initiation factor 4E-binding protein-1
There were significant differences in hepatic enzyme activity and serum biochemical parameters among all treatments in male broilers. In the transamination pathway, the threonine restriction group had the highest, whereas the lysine and arginine supplementation groups had the lowest hepatic enzyme activities and serum metabolites (Fig. 2A). The tryptophan supplementation group had the highest, whereas the tryptophan restriction group had the lowest fat synthesis pathway (Fig. 2B). In the gluconeogenesis pathway, the methionine- and cysteine-supplemented groups had the highest levels, and the methionine- and cysteine-restricted groups had the lowest levels (Fig. 2C). Regarding the cholesterol synthesis pathway, the tryptophan supplementation group showed the highest expression, and the BCAA restriction group showed the lowest expression (Fig. 2D). In the Uric acid excretion pathway, the threonine supplementation group had the highest expression, and the normal protein control group had the lowest expression (Fig. 2E).

Effects of amino acid combinations with different patterns on hepatic enzyme activities and serum biochemical parameters of AA broilers (male). A Transamination pathway; B Fat synthesis pathway; C Gluconeogenic pathway; D Cholesterol synthesis pathway; E Uric acid excretion pathway. Con(+): The normal-protein control group; Con(-): The low-protein control group; L&A(-): The lysine and arginine restriction group; L&A(+): The lysine and arginine supplementation group; M&C(-): The methionine and cysteine restriction group; M&C(+): The methionine and cysteine supplementation group; BCAA(-): The BCAA restriction group; BCAA(+): The BCAA supplementation group; Thr(-): The threonine restriction group; Thr(+): The threonine supplementation group; Trp(-): The tryptophan restriction group; Trp(+): The tryptophan supplementation group; AST: Aspartate aminotransferase; ALT: Alanine aminotransferase; TG: Triglyceride; ACC: Acetyl CoA carboxylase; FAS: Fatty acid synthase; GLU: Blood glucose; G6PC: Glucose-6-phosphatase; PC: Pyruvate carboxylase; TC: Cholesterol; HMGR: HMG CoA reductase; XOD: Xanthine oxidase; BUN: Urea nitrogen; UA: Uric acid
In female broilers, variations in hepatic enzyme activity among the different treatment groups were fundamentally similar to those observed in males. In the transamination pathway, the threonine supplementation group had the highest, whereas the BCAA restriction group had the lowest hepatic enzyme activities and serum metabolites (Fig. 3A). The tryptophan supplementation group had the highest and the tryptophan restriction group had the lowest fat synthesis pathway (Fig. 3B). In the gluconeogenesis pathway, the methionine- and cysteine-supplemented groups had the highest values, while the lysine-and arginine-supplemented groups had the lowest values (Fig. 3C). In the Cholesterol synthesis pathway, the tryptophan supplementation group showed the highest expression, and the tryptophan restriction group showed the lowest expression (Fig. 3D). In the Uric acid excretion pathway, the threonine supplementation group had the highest expression and the lysine and arginine supplementation groups had the lowest expression (Fig. 3E).

Effects of amino acid combinations with different patterns on hepatic enzyme activities and serum biochemical parameters of AA broilers (female). A Transamination pathway; B Fat synthesis pathway; C Gluconeogenic pathway; D Cholesterol synthesis pathway; E Uric acid excretion pathway. Con(+): The normal-protein control group; Con(-): The low-protein control group; L&A(-): The lysine and arginine restriction group; L&A(+): The lysine and arginine supplementation group; M&C(-): The methionine and cysteine restriction group; M&C(+): The methionine and cysteine supplementation group; BCAA(-): The BCAA restriction group; BCAA(+): The BCAA supplementation group; Thr(-): The threonine restriction group; Thr(+): The threonine supplementation group; Trp(-): The tryptophan restriction group; Trp(+): The tryptophan supplementation group; AST: Aspartate aminotransferase; ALT: Alanine aminotransferase; TG: Triglyceride; ACC: Acetyl CoA carboxylase; FAS: Fatty acid synthase; GLU: Blood glucose; G6PC: Glucose-6-phosphatase; PC: Pyruvate carboxylase; TC: Cholesterol; HMGR: HMG CoA reductase; XOD: Xanthine oxidase; BUN: Urea nitrogen; UA: Uric acid
Expressions of mRNA in breast muscle
There were significant differences in S6K1 and MAFbx expression but no significant differences in the other indices between male (Fig. 4) and female broilers (Fig. 4 and S1). The threonine supplementation group had significantly lower levels than the other treatments, and the normalized protein control group and the lysine and arginine supplementation groups had the highest levels of S6K1. The threonine supplementation group had significantly higher MAFbx than the other treatments, and the lysine and arginine supplementation groups had the lowest MAFbx.

Effects of amino acid combinations with different patterns on alterations of mRNA expressions in the breast muscle of AA broilers. A Male-S6K1; B Male-MAFbx; C Femal-S6K1; D Female-MAFbx. Con(+): The normal-protein control group; Con(-): The low-protein control group; L&A(-): The lysine and arginine restriction group; L&A(+): The lysine and arginine supplementation group; M&C(-): The methionine and cysteine restriction group; M&C(+): The methionine and cysteine supplementation group; BCAA(-): The BCAA restriction group; BCAA(+): The BCAA supplementation group; Thr(-): The threonine restriction group; Thr(+): The threonine supplementation group; Trp(-): The tryptophan restriction group; Trp(+): The tryptophan supplementation group. S6K1: The 70-kDa ribosomal protein S6 kinase; MAFbx: Muscle atrophy f-box. a–cDissimilar letters represent significant difference among different treatments (P < 0.05)
Hepatic untargeted metabolomics analysis
Based on the screening of the data in Figs. 2, 3 and 4, as well as the sex-specific sensitivity of animals to amino acid balance. Male broiler chickens were used as models. The normal protein control group was fed a normal protein diet (NP), the low-protein control group was fed a low-protein diet (LP), the lysine and arginine supplementation group was fed a low-protein amino acid-balanced diet (LPAB), and the threonine supplementation group was fed a low-protein amino acid-imbalanced diet (LPAI). These diets were subjected to untargeted metabolomic analysis to identify amino acid balance biomarkers.
Principal component analysis
The PCA score plots revealed a significant separation between the NP and LP, with PC1 and PC2 accounting for 69.6% and 8.8% of the sample variation, respectively (Fig. 5A). Similarly, there was a notable distinction between LPAB and LPAI, as shown in the PCA score plots, where PC1 and PC2 explained 35.9% and 19.5% of sample variation, respectively (Fig. 5B).

Principal component (PC) analysis of NP vs. LP liver metabolomics data (A) and LPAB vs. LPAI liver metabolomics data (B). Volcano of NP vs. LP liver differential metabolites (C) and LPAB vs. LPAI liver differential metabolites (D). Venn plot of NP vs. LP and LPAB vs. LPAI shared liver differential metabolites (E). NP: Normal CP diet; LP: Low-CP diet; LPAB: Low-CP amino acid balance diet; LPAI: Low-CP amino acid imbalance diet
Orthogonal partial least-squares discriminant analysis
The OPLS-DA score plot demonstrated a significant separation between the NP and LP, with the abscissa and ordinate explaining intragroup differences, accounting for 40.4% and 17.2% of the sample variation, respectively (Fig. S1A). A permutation test (n = 100) was performed to validate the OPLS-DA model. The blue bands depict the original samples, where a more pronounced distribution toward the right indicates more substantial data separation between the groups. The R2Y(cum) and Q2(cum) values were 0.988 and 0.823, respectively (Fig. S1C), confirming the reliability of these results.
Similarly, the OPLS-DA score plot exhibited a significant separation between LPAB and LPAI, with inter-group differences in the abscissa and ordinate accounting for 21% and 21.6% of the sample variation, respectively (Fig. S1B). The R2Y(cum) value was calculated to be 0.994, whereas the Q2(cum) value reached 0.857 (Fig. S1D), indicating the high reliability of these results.
Identification of significant differential metabolites
Volcano plots revealed 2,528 upregulated metabolites and 3,150 downregulated metabolites between the NP and LP groups, resulting in the identification of 431 differential metabolites (Fig. 5C). Additionally, 932 upregulated and 1,872 downregulated metabolites were detected between LPAB and LPAI based on volcano plot analysis, leading to the final identification of 248 differential metabolites (Fig. 5D). Overall, 167 differential metabolites were identified in the four treatments (Fig. 5E).
Acquisition of specific potential biomarkers
To identify the relevant differential metabolites among more than 100 metabolites, the following criteria were used: P < 0.05, VIP > 1, and fold-change > 2 or < 0.5. A total of 43 metabolites met the inclusion criteria, in addition, uric acid and urea nitrogen were counted as common amino acid metabolites (Table 11). In broiler production, certain technical procedures such as handling, and blood sampling can induce stress in animals and potentially interfere with in vivo biochemical indices. Therefore, based on the identified differential metabolites, we selected those that exhibited consistent changes in direction and quantitative difference values (QD-value) > 1. The 3 metabolites that met these criteria were Oxypurinol, Pantothenic acid, and D-octopine (Table 12).
The potential of metabolites as biomarkers was evaluated using multivariate exploratory ROC curve analysis. Boxplots demonstrated significant differences in these 3 metabolites among the four treatments (Fig. 6A–C), and subsequent ROC analysis revealed that their AUCs were all equal to 1 (Fig. 6D–F), indicating their potential as biomarkers. Consequently, oxypurinol, pantothenic acid, and d-octopine are considered viable candidates for amino acid balance biomarkers.

Boxplot of oxypurinol (A), pantothenic acid (B) and D-octopine (C). Receiver operator characteristic (ROC) curve analysis of oxypurinol (D), pantothenic acid (E) and D-octopine (F). The data of 3 metabolites are expressed as relative quantification (intensity)
Metabolic pathway analysis
Pathway analysis revealed the precise impact of amino acid balance alterations on the relevant metabolic networks, and the extent of enrichment in metabolic pathways was quantified using the rich factor (RF) and percentage. Percentages were calculated by dividing the number of differential metabolites in each pathway by the total number of differential metabolites in the sample. RF represents the ratio of differential metabolites to all metabolites within the same pathway (Fig. 7). The significance level was set at P < 0.05, with RF > 0.1, and percent > 5. A total of four metabolic pathways were identified: nucleotide metabolism, pantothenate and CoA biosynthesis, fatty acid biosynthesis, and glutathione metabolism (Table 13). These pathways are associated with amino acid balance and can serve as potential targets for future studies on amino acid balance.

KEGG classification of NP vs. LP (A) and LPAB vs. LPAI (B). KEGG enrichment bubble of NP vs. LP (C) and LPAB vs. LPAI (D). NP: Normal CP diet; LP: Low-CP diet; LPAB: Low-CP amino acid balance diet; LPAI: Low-CP amino acid imbalance diet
Discussion
After adjusting for the patterns of EAA combinations in low-protein diets, significant differences were observed in various broiler performance indices, indicating the effectiveness of improving animal performance through EAA pattern adjustments. Notably, regardless of sex, broilers fed the low-protein control diet exhibited significantly poorer growth performance than those fed the normal-protein control diet. Surprisingly, a slight advantage was observed in broilers of the low-protein control group during the starter phase. This could be attributed to the presence of a yolk sac in chicks and the relatively less vigorous requirement of exogenous nutrients during the initial growth phase [29]. The high protein content in the normal-protein control group provided adequate protein nutrition with minimal feed intake; however, it also resulted in lower energy intake, which negatively affected growth performance during this phase. Once fully absorbed, a rapid improvement in growth performance was evident for broilers in the normal protein control group during the subsequent phases. In terms of carcass characteristics, broilers fed the normal-protein control diet displayed higher breast and fat pad rates, suggesting a balanced amino acid profile in which sufficient amino acids were utilized for protein synthesis rather than oxidative decomposition, leading to fat accumulation within their bodies. Although not significantly different from the normal-protein control group, broilers in the low-protein control group exhibited a reduced breast rate but comparable fat pad rates. These results indicate that although there remained a balance in amino acid ratios within the low-protein control group, insufficient amounts of amino acids were available for optimal protein synthesis, possibly because of the lower total nitrogen content [15].
The growth performance and breast rate of broilers fed the Lysine and Arginine supplemented diets were significantly improved in diets with adjusted lysine and arginine combination ratios for both male and female broilers. Notably, there was a decrease in dressing, half-evisceration, and evisceration rates after reducing the dietary concentrations of lysine and arginine, which was particularly significant in males. Lysine is the second limiting amino acid in broilers and plays an important role in protein synthesis without being involved in various complex physiological functions [30]. After supplementing low-protein diets with lysine, in addition to meeting the requirements for protein synthesis, along with other amino acids, excess lysine is oxidatively decomposed in the body to form substrates for NEAAs synthesis [31]. This leads to increased concentrations of available amino acids and enhanced protein synthesis within the body. The remaining substances that are not involved in synthesis enter the TCA cycle and eventually contribute to fat deposition, thereby explaining why the breast cancer rate increases simultaneously with fat retention when supplemented with lysine and arginine. When the dietary lysine concentration is low, priority is given to allocating amino acids to synthesize proteins essential for the functioning of vital organs and feather development [32]. Consequently, this results in decreased overall meat yield from broilers, that is, reduced dressing, half-eviscerated, and eviscerated rates. Therefore, increasing lysine and arginine levels in low-protein diets is recommended.
The reduction in methionine and cysteine concentrations in the diet resulted in an improvement in FCR in both male and female subjects and a numerical increase in BW. However, it also led to a significant increase in the fat pad rate, whereas the breast rate remained lower than that in the normal protein control group. The role of methionine in poultry has been extensively studied [33, 34]. Apart from its involvement in protein synthesis, methionine can also enhance the overall body metabolism. When the dietary methionine supply is reduced, animals increase their feed intake to meet normal growth demands, thereby improving amino acid availability within the body and subsequently reducing the FCR while increasing the breast rate. Methionine also plays a regulatory role in the lipid metabolism. Insufficient methionine affects CPT1-mediated fatty acid entry into the mitochondria, leading to excessive fat deposition [35]. Conversely, excessive methionine levels can result in hyperhomocysteinemia in broilers, impairing lipid metabolic processes [36]. This explains why a significant increase in the fat pad rate was observed regardless of whether the methionine concentration was upregulated or downregulated in broilers. Therefore, when formulating low-protein diets, large variations in methionine and cysteine concentrations should be avoided.
The growth performance and carcass characteristics of both male and female animals fed diets with adjusted ratios of BCAA combinations were significantly inferior to those of the control group. BCAAs exhibit antagonistic effects on each other [37]; this study found that adjusting the dietary concentration based on the original BCAA ratio did not prevent the occurrence of these antagonistic effects. This suggests that the interactions between the BCAAs are likely influenced by their absolute concentrations. When adjusting the ratio of BCAA combinations, it is important to consider the weighting of the individual members within the combination. Furthermore, increasing amino acid concentrations according to the current balance patterns when exploring amino acid balance patterns in low-protein diets is not recommended. Additionally, broilers fed a diet containing higher concentrations of BCAAs (the BCAA supplementation group) exhibited worse performance than those fed lower concentrations (the BCAA restriction group), indicating a more pronounced antagonistic effect with increased BCAA concentration. An increase in the grain feed proportion in low-protein diets also leads to an increase in leucine concentration [38]. Leucine is one of the most antagonistic amino acids within BCAAs [39], whereas valine is currently considered the fourth limiting amino acid in poultry [40]. Therefore, future research on BCAAs should consider the relationship between leucine and valine to minimize the negative effects caused by their antagonistic effects.
Growth performance and carcass characteristics were negatively affected by the ratio of dietary threonine after adjustment, with broilers fed a threonine supplemented diet exhibiting the poorest outcomes among all treatments, indicating an extremely unbalanced amino acid pattern. Threonine plays a crucial role in maintaining amino acid balance in the diet and can partially alleviate growth inhibition caused by excessive lysine, methionine, and tryptophan [41]. This suggests that threonine may inhibit the utilization of other amino acids. There have been previous reports on the adverse effects of excess threonine [5]. In future studies, it will be important to meet the minimum requirement of threonine without additional supplementation to maintain an optimal amino acid balance in the diet.
Despite its low dietary concentration, tryptophan cannot be disregarded for its physiological regulatory effects on broilers. Tryptophan supplementation in low-protein diets significantly improved the growth performance and carcass characteristics of both male and female broilers. Notably, the average daily feed intake decreased, which could be attributed to the presence of 5-hydroxy-tryptamine (5-HT), a tryptophan metabolite. 5-HT regulates poultry feeding behavior, and excessive levels can induce satiety and reduce feed intake [42]. Tryptophan stimulates protein synthesis by promoting insulin secretion [43], thereby increasing meat yield. This study revealed that adjusting the tryptophan ratio has a significant effect on thigh yield, particularly in males, while having a minimal effect on breast yield. Further studies are required to better understand these preferences. Moreover, there was no significant increase in the fat pad rate with high concentrations of tryptophan, indicating that supplementing low-protein diets with tryptophan improves the amino acid balance. However, maintenance of animal feed intake without reduction at high tryptophan concentrations should be prioritized in future research.
It is imperative to promptly and accurately assess the dietary amino acid balance to establish a low-protein diet amino acid balance pattern, while elucidating the adjustment trends of each EAA in low-protein diets. Therefore, we first described the entire amino acid metabolism pathway (Fig. 1) and determined different treatments at the transcriptional, protein, and metabolic levels to obtain the optimal and worst amino acid balance treatments for subsequent screening of amino acid balance biomarkers. Based on the results of hepatic enzyme activities and serum metabolites in the amino acid metabolism pathway, it was evident that the threonine supplementation group had the highest levels of both transamination and urinary nitrogen excretion pathways. When there is an imbalance of amino acids in the diet, the amino acids are oxidized and decomposed because of their inability to participate in protein synthesis in a timely manner, leading to increased enzyme activities and metabolite content in the amino acid metabolism pathways [44]. The transamination pathway represents the initiation of amino acid metabolism and decomposition, whereas the gluconeogenesis, fat synthesis, and cholesterol synthesis pathways represent the fate of decomposed amino acids within the body, and the urinary nitrogen excretion pathway indicates the degree of waste from these processes. In the threonine supplementation group, excessive threonine undergoes oxidation and decomposition within the body before being converted into glycine [45], which directly participates in uric acid synthesis [46], resulting in excess nitrogen not being retained as a functional NEAA. This might explain the high nitrogen excretion observed in this treatment. Therefore, from the perspective of protein balance and metabolism, the threonine supplementation group appeared to be the most imbalanced, while the Lysine and Arginine supplementation groups displayed lower enzyme activities and metabolite content across various pathways, indicating greater balance at these two levels. In addition, it was also observed that the tryptophan supplementation group exhibited the highest content in both the fat and cholesterol synthesis pathways. Tryptophan reduces blood glucose levels [47] and promote insulin secretion [48]. Insulin regulates liver metabolism by enhancing fat synthesis, inhibiting fatty acid oxidation, and promoting triglyceride esterification [49]. Consequently, it was anticipated that the tryptophan-supplemented group would exhibit elevated triglyceride and cholesterol levels. Based on the expression of mRNA related to protein synthesis in breast muscle, it was noted that the threonine supplementation group displayed the lowest expression for protein synthesis, while exhibiting the highest expression for degradation. Conversely, this trend was reversed in the lysine- and arginine-supplemented groups, indicating an imbalance at the transcriptional level within the threonine-supplemented group compared with the more balanced state within the lysine- and arginine-supplemented groups.
Based on the results obtained at the 3 levels, the normal-protein control group was designated as the normal control diet (NP), the low-protein control group as the low-protein control diet (LP), the lysine and arginine supplementation group as the low-protein amino acid balance diet (LPAB), and the threonine supplementation group as the low-protein amino acid imbalance diet (LPAI) for untargeted metabolomic analysis aimed at exploring the biomarkers of amino acid balance. Serum is commonly used to assess the physiological state of an animal. However, serum biochemical indicators are often influenced by animal activities, such as stress, leading to the detection of metabolites unrelated to the experiment and compromising study accuracy when exploring metabolism and identifying biomarkers. Considering that the liver serves as the primary site for amino acid metabolism, metabolites are synthesized and transferred into the bloodstream through the portal vein, and investigating changes in amino acid metabolism in the liver tissue provides a more representative approach. Therefore, liver tissue was selected for analysis in this study. Following metabolite identification and statistical analysis, 3 metabolites were identified as biomarkers of the amino acid balance: Oxypurinol, pantothenic acid, and D-octopine. Their associations with amino acid metabolism are shown in Fig. 8. Oxypurinol is an active metabolite and analog of xanthine formed through the oxidation and decomposition of its precursor, allopurinol [50]. Allopurinol, which resembles hypoxanthine, competes with xanthine oxidase (XOD) during uric acid formation. It is commonly used to treat gout [51]. The primary metabolic mechanism involves the XOD-mediated conversion of allopurinol to oxypurinol by rotation at its enzyme active site, leading to the reduction of the Mo (VI) active center to Mo (IV), thereby inhibiting uric acid production [52]. In this study, oxypurinol content significantly increased in both amino acid balance treatments (NP and LPAB), despite no exogenous supplementation of allopurinol. This can be attributed to only a small amount of surplus amino acids entering catabolic processing under amino acid balance conditions, resulting in insufficient NEAAs such as aspartic acid, glutamine, and glycine, which serve as raw materials for de novo purine synthesis [53]. Consequently, there is a decrease in hypoxanthine and xanthine content, the original substrates for XOD, which enhances the competitiveness of allopurinol, leading to an increased generation of oxypurinol. Additionally, we observed an increase in xanthine concentration (Table 12); this finding is not contradictory. The absolute utilization of XOD by allopurinol inhibits the timely excretion of xanthine, similar to that of uric acid. However, it cannot be converted back to hypoxanthine and thus accumulates in significant quantities in the body. Notably, oxypurinol exhibited the highest relative abundance among all metabolites and displayed substantial variation across different treatments. Importantly, it demonstrates a close association with the end products of amino acid metabolism. Therefore, we propose that oxypurinol is the most representative biomarker for assessing amino acid balance.

Related pathways of amino acid balance biomarkers involved in amino acid metabolism
Homeostatic regulation of pantothenic acid and vitamin B5 in the body involves their interactions with pantothenic acid and coenzyme A [54]. In the presence of pantothenate deficiency, coenzyme A is dephosphorylated by alkaline phosphatase to form dephospho-CoA, which then undergoes hydrolysis of its 5′-phosphodiester bond to yield phosphopantetheine. The subsequent step involves the synthesis of pantetheine, catalyzed by alkaline phosphatase, which is ultimately converted to pantothenic acid [55]. In the absence of coenzyme A, pantothenic acid is phosphorylated to 4′-phosphopantothenate by pantothenate-kinase kinase, which subsequently condenses with cysteine. Decarboxylase removes the carboxyl group from this compound to generate dephosphocoenzyme A, which eventually leads to coenzyme A [56]. The amino acid-balanced diets (NP and LPAB) were expected to exhibit reduced pantothenic acid content in this study. The ultimate goal of amino acid metabolism is to produce coenzyme A, which enters the tricarboxylic acid cycle and participates in other substance synthesis processes [57]. As mentioned earlier, when amino acids are balanced, fewer amino acids participate in catabolic processes, resulting in decreased levels of their decomposition products such as coenzyme A (Table 12). Under homeostatic regulation, synthesis reactions favor conversion towards increased levels of coenzyme A, leading to a decrease in the overall pantothenic acid content.
The energy metabolite d-octopine, derived from the combination of octopine dehydrogenase (OcDH) and NADH [58], undergoes a series of metabolic reactions: first, NADH binds to domain 1 of OcDH to form a complex; subsequently, arginine binds to domain 2 of OcDH, resulting in the formation of an OCDH-NADH-L-arginine complex; and finally, pyruvate binds to the last binding site within the complex and promptly condenses into D-octopine while releasing H+ and NAD+ [59, 60]. In an amino acid-balanced diet, where there is an abundance of raw materials for protein synthesis, both the quantity and rate of protein synthesis increase with higher energy requirements. Consequently, a significant increase in D-octopine content was observed in the NP and LPAB groups. During amino acid utilization, ensuring a coordinated energy supply is crucial for animal growth [61]. This is particularly true when animals are on an amino acid-balanced diet because rapid protein synthesis may lead to insufficient energy levels, causing non-participating amino acids to enter catabolic processes, thereby affecting the evaluation index values. Such circumstances can result in erroneous assessment of dietary amino acid balance patterns. Therefore, it would be meaningful to conduct experiments on dietary amino acid balance patterns while ensuring an adequate energy supply for animals.
These 3 metabolites are not commonly observed in the investigation of dietary amino acids. Therefore, we selected the common metabolites uric acid and urea from the metabolomic data for comparative analysis (Tables 11 and 12). However, these 2 metabolites exhibited relatively low evaluation scores, and their differences were insignificant when minor variations in the degree of balance of dietary amino acids were considered. In particular, the urea content was within that of the LPAB and LPAI groups. These findings suggest that the previous indicators used to assess amino acid balance lack sufficient sensitivity, especially in studies with minimal differences, which may lead to erroneous evaluations. Oxypurinol is an upstream metabolite of uric acid, whereas pantothenic acid acts upstream of tricarboxylic acid cycle products such as cholesterol and triglycerides. In contrast, D-octopine plays a central role in energy metabolism. They are highly sensitive to the amino acid balance and do not disrupt normal amino acid metabolism. Furthermore, based on the results of metabolic pathway enrichment analysis, oxypurinol corresponded to nucleotide metabolism, whereas pantothenate corresponded to pantothenate and CoA biosynthesis. Therefore, these 3 metabolites serve as reliable indicators for assessing the amino acid balance. In the subsequent amino acid balance study, the evaluation was based on the direction of the metabolite changes provided in Table 12. Specifically, higher levels of oxypurinol and d-octopine, along with lower levels of pantothenic acid, indicated a more balanced composition of dietary amino acids in this group.
Conclusion
In conclusion, adjusting the patterns of EAAs in low-protein diets is required to enhance amino acid balance and improve broiler growth performance. Untargeted metabolomic analysis identified 43 metabolites and 4 metabolic pathways associated with amino acid balance, including hepatic oxypurinol, pantothenic acid, and D-octopine, as potential biomarkers for assessing dietary amino acid balance. By detecting these biomarkers, the degree of dietary amino acid balance can be precisely evaluated. This study provides valuable insights for further exploration of amino acid balance patterns in low-protein diets of broilers.