The rate of protein accretion and growth affect amino acid requirements in young animals. Differences in amino acid metabolism contribute to individual variations in growth rate. This study aimed at determining how amino acid needs may change with growth rates in broiler chickens. Experiment 1 consisted of testing amino acid choices in two chicken groups with extreme growth rates (the slowest –SG- or fastest –FG- growing birds in a flock). Essential (EAA) (methionine, lysine and threonine) or non-essential (NEAA) (alanine, aspartic acid and asparagine) amino acids were added to a standard control feed (13.2 MJ/kg; 21.6% crude protein). The chickens were offered simultaneous access to the control feed and a feed supplemented with one of the two amino acid mixes added at 73% above standard dietary levels. Experiment 2 consisted of the selection of the bottom 5 SG and top 5 FG chickens from a flock of 580 to study differences in amino acid metabolism using the proventriculus representing gut sensing mechanism. In this experiment, transcriptomic, proteomic, and genomic analyses were used to compare the two groups of chickens.
SG preferred NEAA, while they rejected EAA supplemented feeds (P < 0.05). However, FG rejected NEAA (P < 0.05), and they were indifferent to EAA supplemented feed (P > 0.05). Transcriptomic and proteomic analyses identified 909 differentially expressed genes and 146 differentially abundant proteins associated with differences in growth rate (P < 0.05). The integration of gene expression and protein abundance patterns showed the downregulation of sensing and transport of alanine and glucose associated with increased alanine catabolism to pyruvate in SG chickens.
Dietary preferences for NEAA in the SG group are associated with a potential cytosolic depletion of alanine following an upregulation of the catabolism into TCA cycle intermediates.
Chickens have been used as a model organism for more than a century for research into embryology, immunology, genetics, cell biology, cancer, virology and growth [1,2,3]. Broiler chickens have been selected for efficient body weight gain for decades with modern breeds reaching one of the highest growth rates across livestock and laboratory animals. However, genetic selection has not been successful in improving growth uniformity within chicken populations.
Weight gain is strongly correlated to feed intake. Many factors contribute to feed intake regulation in chickens . Imbalances in available amino acids to energy in the diet can depress or enhance feed intake . The fastest growing (FG) chickens in a flock deposit more body protein and have higher amino acid requirements than the average. These birds increase feed intake to compensate for the increased need for amino acids [6,7,8,9]. In contrast, slow growing (SG) chickens consume excessive amino acids from the diet relative to the flock average, which can supress their appetite and decrease feed intake [10, 11]. Broiler feeds are formulated primarily to meet the requirements of the limiting essential amino acids (EAA) [7, 12]. In contrast, the supply of non-essential amino acids (NEAA) is rarely considered. The NEAA are considered to be synthesised in sufficient quantity from metabolic precursors within the animal’s cell. However, mounting evidence is building in the scientific literature showing that NEAA can stimulate feed intake and growth rate of chickens, particularly when offered in low protein diets [6, 13].
Nutrient sensing is an essential mechanism enabling broilers to select diets that balance their nutrient requirements, when given a choice of feeds varying in nutrient composition [14, 15]. The molecular mechanisms that mediate choice and specific appetites for amino acids in chickens have been related to the taste receptor dimer T1R1/T1R3, amongst other receptors [15, 16]. However, a systematic study on preferences between diets varying in amino acid concentrations has not been reported in broilers. Potential differences in amino acid sensing mechanisms may explain lower feed intake and slow growth in SG compared to FG chickens.
The objective of this research was to study differences in specific amino acid appetites in broilers with extreme high or low growth rates. The hypothesis tested was that specific amino acid appetites associated with low growth rates indicate innately constrained metabolic pathways affecting feed intake and protein accretion in SG broiler chickens.
The selection of the 10% lightest and heaviest individuals in a flock resulted in 29% higher (P < 0.05) final body weight (42 days) and feed intake in FG compared to SG chickens. Double choice (DC) treatments did not significantly (P > 0.05) affect growth or total feed intake during the experimental period. A summary of the main performance parameters and statistical significance of main effects have been presented in the Table S3 in Additional file 1.
The SG birds showed a significant (P < 0.05) preference for the NEAA (Ala, Asp, Asn) supplemented feed compared to the control feed (Fig. 1A). In contrast, SG chickens rejected (P < 0.01) the EAA (Met, Lys and Thr) supplemented feed compared to the control feed. The FG chickens significantly (P < 0.01) chose against the NEAA supplemented feed while showing no preference (P > 0.05) for the EAA supplemented feed (Fig. 1A).
FG birds consumed significantly (P < 0.01) more total feed (control + test feed) than SG chickens across all treatments (Fig. 1B). The analysis of covariance showed that SG chickens consumed significantly less (P < 0.05) EAA (Met, Lys, Thr), but more NEAA (Ala, Asp, Asn) than the FG birds after adjusting for total feed intake as a covariate (Fig. 1C).
Gene and protein expression studies
RNA sequencing (RNAseq) analysis of the proventriculus found 909 DEG, and the proteomic analysis found 146 DAP when comparing SG to FG chickens (Fig. 2). The qPCR validation showed a high correlation coefficient (− 0.89; P < 0.0001) between gene expression of the reported genes measured using qPCR and RNAseq (Fig. S2 in Additional file 1). For example, both RNAseq and qPCR analyses showed that gene expression of α-gustducin was significantly down-regulated in SG compared to FG (Table S4 in Additional file 1).
The DEG significantly (P < 0.05) enriched several biological processes, cellular components and molecular functions as shown in Fig. 3 (Gene Ontology). In addition, the results showed several enriched metabolic pathways including the ‘Glyoxylate Metabolism and Glycine Degradation’ and ‘Pyruvate Metabolism’ pathways (KEGG and REACTOME; Table 1). Five DEG (DLAT: dihydrolipoamide S-acetyltransferase; LIAS: Lipoic acid synthetase; OGDH: oxoglutarate (α-ketoglutarate) dehydrogenase; PDHB: pyruvate dehydrogenase beta; PDHX: pyruvate dehydrogenase complex component) relevant to the two pathways were significantly upregulated in SG chickens (Table 2, Fig. 4). In contrast, gene expression for several amino acid sensors and transporters including Gust (Gustducin, subunit alpha), SLC38A1 (Solute carrier 38 type A1), SLC1A2 (Solute carrier 1 type A2), and the glucose transporters SLC2A1 (Solute carrier 2 type A1) and SLC2A10 (Solute carrier 2 type A10) were downregulated in SG chickens (Fig. 4). The corresponding protein names of the genes SLC38A1, SLC1A2, SLC2A1, and SLC2A10 are SNAT1, EAAT2, GLUT1, and GLUT10, respectively.
The differentially abundant proteins significantly (P < 0.05) enriched the metabolic pathways of ‘Translocation of GLUT4 to membrane’, ‘The Tricarboxylic Acid Cycle’, ‘Glycolysis’, ‘Gluconeogenesis’, ‘Biosynthesis of Amino Acid’, and ‘Arginine and Proline Metabolism’ (Table 1). Protein abundance for GLUT2 (Glucose transporter 2), YWHAE (Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein epsilon), YWHAQ (Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein theta), YWHAZ (zeta), ENO1 (Enolase 1), PKM (Pyruvate kinase), and CS (Citrate synthase) were significantly downregulated (P < 0.05) in SG chickens, whereas that of SDHA (Succinate dehydrogenase complex subunit A), SDHB (Succinate dehydrogenase complex subunit A), GOT1 (Glutamic-oxaloacetic transaminase 1), CK (Creatine Kinase), AGMAT (Agmatinase), and TPI1 (Triosephosphate isomerase 1) were upregulated (P < 0.05) (Table 2, Fig. 5, Fig. S1 in Additional file 1).
The genomic analysis identified over 900,000 variants on average in each individual (n = 5) with accuracy of 99.9% (Phred score ≥ 30). From those variants, 186 SNPs were on the genes of interest with deleterious effect on protein products based on Variant Effector Predictor analysis (SIFT score < .05). Of these, 15 SNP variants were affecting amino acid sensors/transporters and 5 variants affected glucose sensors/transporters (Table 3). One deleterious point mutation on the SLC2A4RG gene (20:9813052; A➔T) showed a significantly higher (Fisher’s Exact Test value = 0.007) allelic frequency in slow- compared to FG chickens (Table 3).
When given the choice, chickens can select diets that meet their nutrient requirements . We report the results of double choice preference tests used to assess selection of feeds with different amino acid content in FG compared to SG chickens. The control diet was a standard commercial feed formulated to meet or exceed the EAA requirements of the average bird of the flock. A single feed does not accommodate differences in nutrient requirements for chickens growing at different rates. FG chickens with high rates of protein deposition will be deficient while SG chickens with low rates of protein deposition will have excess of dietary amino acids. An excess of EAA or NEAA were added to the feed to identify potential dietary deficiencies. EAA cannot be synthesised by chickens, while NEAA have been defined to be synthesised from precursor compounds in chickens at a non-limiting rate.
The results confirmed a close relationship between feed intake and growth rate [18, 19], with the 10% fastest growing of a commercial flock consuming approximately 15% more feed than the 10% slowest growing. The potential amino acid deficiencies in FG chickens were compensated by an increased feed consumption compared to SG chickens. Previous reports [10, 20] described a loss of appetite and decreased feed intake in broiler chickens when excess EAA were added to a balanced diet. In our study, this response was observed only with the SG birds, reflecting the lower rate of protein deposition and a lower requirement for EAA for protein synthesis compared to FG chickens.
The opposite was true for NEAA. While rejecting the feed with excess EAA SG birds showed a strong preference for the feed with additional NEAA (see Fig. 1 A). This increased consumption of Ala, Asp and Asn showed a deliberate choice for NEAA by SG birds. These findings are consistent with the need for added NEAA in animal diets . Previous work in chickens identified the relevance of dietary NEAA in sustaining or improving growth performance [22, 23]. In addition, excess of EAA or NEAA resulted in depressed feed intake in chickens showing a higher response to the essential than the non-essential .
EAA are not always used efficiently as a source of nitrogen . Thus, an optimum ratio of EAA to NEAA in animal diets has been suggested for maximum growth and protein utilization [23, 25]. The optimum ratio in chicken diets has been identified as around 0.50 indicating twice the amount of EAA than NEAA [25, 26]. Our results showed that the optimum ratio of the three most limiting dietary amino acids in chickens (Met, Lys, Thr) with three non-essential (Ala, Asp, Asn) observed was 0.59 for FG and 0.57 for SG birds. The results suggest FG birds perceive a greater amino acid imbalance when diets contain higher concentrations of NEAA and SG birds when diets contain more EAA.
Recently, Hofman and co-workers  observed a decrease in growth rates in broilers fed a low crude protein diet providing adequate levels of EAA. The decrease in performance could not be compensated only by supplementing Gly (or Gly equivalent including serine) which led to the conclusion that another nutrient (possibly NEAA) was limiting the growth. In the chicken gut, neutral amino acids compete for transport. In particular, Met and Leu inhibit Ala and Gly absorption . This may explain why postprandial plasma levels of Ala decreased as dietary protein levels increased above 15% in broiler chickens .
The proventriculus is the glandular stomach in birds, which was selected as a proxy for the chicken gastrointestinal tract given that thanks to reverse peristalsis this organ plays a key role in feed digestion, nutrient sensing, and feed intake regulation. These mechanisms involve the secretion of hunger/satiety hormones such as ghrelin, gastrin, and cholecystokinin [30, 31]. For example, the proventriculus mediates gastrointestinal functions, including pancreatic secretions . The analysis of proventricular tissues comparing SG to FG birds, showed robust differences in the catabolism of NEAA, particularly for Ala which would imply a decrease on the availability of these amino acids for protein synthesis in extra-intestinal organs such as skeletal muscle . The higher catabolic rate of Ala in SG compared to FG chickens was due to a seemingly orchestrated upregulation of four metabolic pathways: the ‘glyoxylate metabolism and glycine degradation’, ‘pyruvate metabolism’, ‘gluconeogenesis’, and ‘the tricarboxylic acid cycle’ (Fig. 6).
Dietary Arg and Pro are partially metabolised in the intestine . We observed an upregulation of the ‘arginine and proline metabolism’ in SG birds resulting in the production of glyoxylate via hydroxyproline further increasing the need for Ala [35,36,37]. Shortage of Ala can lead to accumulation of glyoxylate, which converts to oxalate . Oxalate is an anti-nutritional compound present in cereal grains such as wheat and barley and in pulses such as soybean. Ala, as an amino group donor, has a key role in detoxification of oxalate and glyoxylate via Alanine-glyoxylate transaminase (AGT) [38,39,40].
Higher consumption of aspartic acid in SG chickens could be explained by the upregulation of the ‘arginine and proline metabolism’ and the ‘gluconeogenesis’ metabolic pathways (Fig. 6, Fig. S1 in Additional file 1). Aspartate-arginosuccinate shunt is the pathway that connects urea cycle to tricarboxylic acid (TCA) cycle [41, 42]. Upregulation of the key enzyme GOT1 in SG birds suggests the potential use of Asp, along with Ala as a source of energy and a key intermediate for the urea cycle.
In chickens, Ala is the primary transporter of nitrogen and hydrocarbon molecules between liver and skeletal muscles [43, 44]. A deficiency of Ala uptake in the gut lumen can reduce the availability of the amino group and hydrocarbon skeleton for amino acid and protein synthesis in extra-intestinal organs. SG compared to FG chickens had lower gene expression levels of the amino acid transporters SLC38A1 and SLC1A2. SLC38A1 belongs to the SLC38A family, which have been related to Ala and Gln transport in mammals . In addition, the alpha subunit of Gustducin, the heterotrimeric G-protein coupled to dimeric amino acid sensors T1R1/T1R3, was also downregulated in the SG birds. In chickens, the main ligand to the dimeric T1R1/T1R3 amino acid sensor known to date is Ala . These observations strongly indicate that SG chickens are less efficient in sensing and transporting Ala, making it a limiting intermediate metabolite for protein accretion and growth.
This study also highlighted significant differences in glucose transport and metabolism between SG and FG chickens. Glucose uptake in slow growers seemed to be impaired as a result of the downregulation of the glucose transporters GLUT1, GLUT2, and GLUT10 (coding genes: SLC2A1, SLC2A2, and SLC2A10, respectively). Furthermore, a deleterious mutation was identified in the SLC2A4RG gene (protein name: GLUT4RG), which is the transcription factor for the glucose transporter GLUT4, coding gene SLC2A4 [47, 48]. The insulin-stimulated translocation to the cytosolic membrane of GLUT4 may also be impaired in SG chickens. This impairment may be explained by the downregulation of the main proteins responsible for the translocation, YWHAE, YWHAQ, and YWHAZ, in SG compared to FG chickens [49,50,51]. Taken together, the data strongly suggests that slow growers have a limited capacity for glucose sensing and uptake relative to fast growers. The point mutation on SLC2A4RG warrants more investigation. It would be interesting to study if this mutation is the main cause of metabolic differences between SG and FG which has not been addressed in this study. Alternatively, for future studies, the genotyping a bigger chicken population (related to this SNP) and/or the application of gene editing technologies may help understand the function of SLC2A4RG in the chicken.
The downregulation of ‘glycolysis’, the metabolic pathway involved in the conversion of glucose to pyruvate, together with the upregulation of ‘pyruvate metabolism’ (converting pyruvate to Acetyl-CoA), ‘gluconeogenesis’ (synthesising glucose from NEAA and oxaloacetate) and ‘the tricarboxylic acid cycle’ would converge in causing a shortage of pyruvate in SG birds. An upregulated TCA cycle such as in SG individuals consumes more pyruvate, whereas the results showed that production of pyruvate is less due to downregulated glycolysis. The scarcity of pyruvate in the cell would then have implications relevant to the cellular availability and production of Ala. Ala inhibits pyruvate kinase to slow glycolysis and induce gluconeogenesis with increased production of glucose from Ala [52, 53]. First, Ala will be catabolized to pyruvate to feed the TCA cycle for energy production , and second, low pyruvate concentration will reduce the biosynthesis of Ala . Importantly, the gut uses primarily NEAA rather than glucose or EAA to produce ATP [33, 55]. Thus, a lower availability of NEAA, particularly Ala, in the intestinal cells would impair nutrient transport and absorption in the gut. Our current results provide evidence of the latter occurring in SG chickens.
A pathway illustrating the key associations between growth rate and Ala and glucose metabolism in chickens is proposed in Fig. 6. The orchestrated variations in metabolic pathways observed in SG chickens seem to lead to cytosolic glucose, pyruvate and Ala scarcity. Cereal-based diets provide a high starch content (> 50%) indicating that dietary glucose is unlikely to be a limiting nutrient in this study. The data presented is consistent with previous reports showing that feed intake and growth can be improved with dietary NEAA supplementation in chickens . In summary, the downregulation of Gustductin through action of the umami taste receptor T1R1/T1R3 and Ala sensor combined with downregulation of amino acid transporter genes SLC1A2 and SLC38A1 results in reduced sensing and reduced transport of Ala in SG birds. Decreased glucose transport into the cell and reduced glycolysis would cause low cytosolic abundance of pyruvate while reducing Ala biosynthesis. In addition, Ala catabolism was enhanced due to an increased demand for TCA cycle intermediates for gluconeogenesis. SG chickens seem to upregulate the synthesis of pyruvate from Ala to feed the TCA cycle while compensating for the lack of glucose-derived pyruvate. These metabolic changes result in a scarcity of cytosolic Ala, which is reflected in a specific appetite for Ala observed in the feeding choice tests. The results presented provide strong evidence that chickens with lower growing rates may require additional dietary NEAA, particularly Ala.
SG compared to FG chickens select diets containing less EAA (Lys, Met, Thr) and more NEAA (Ala, Asp, Asn) in double choice preference tests. The difference in preference reflects differences in the rate of protein synthesis and dietary amino acid requirements. The differential expression of genes and different abundance of proteins point to a robust impairment in the sensing and transporting of Ala and glucose in SG compared to FG chickens. In addition, Ala catabolism is upregulated to compensate the deficiency of pyruvate required for the TCA cycle in a gluconeogenic status. The subsequent depletion of cytosolic Ala would explain the preference for the NEAA in SG chickens. Dietary supplements of NEAA have the potential to decrease variability in feed intake and growth in young animals.
Body weight at 3 weeks of age was used to select FG and SG chickens from a commercial flock of Ross 308 birds (Darwalla Group, Qld, Australia). The weight of 100, randomly selected, birds was used to identify those birds weighing < 800 g as SG and those weighing > 1000 g as FG, with each category representing approximately 10% of the flock. Birds were randomly selected from the flock until there were 36 in each of the FG and SG categories. The selected birds were transferred to individual cages (45 × 35 × 35 cm) equipped with one nipple drinker and two identical feeders (16 × 10 × 7 cm) placed side-by-side. The 3-week-old chickens were kept for another 3 weeks for the experimental phase. They were offered a commercial standard grower feed for the first week of adaptation period, and a commercial standard finisher feed for the following two weeks as part of the double-choice (DC) testing period. The ingredient and nutrient composition of the two standard commercial feeds (Darwalla Group, Mt Cotton, Queensland, Australia) feeds can be found in Table S1 in Additional file 1.
The DC tests were conducted to study feed preference based on EAA or NEAA supplementation by offering ad libitum access to the two feeders one containing the control finisher feed and the other the experimental (supplemented) feed. The preparation of the experimental feeds consisted in spraying the base commercial finisher feed with water (control), or water containing a mix of three EAA (methionine, Met; lysine, Lys; and threonine, Thr) or three NEAA (alanine, Ala; aspartic acid, Asp; asparagine, Asn). The content of amino acids of interest in the control and the EAA or NEAA supplemented diets is presented in Table 4. The added EAA or NEAA were 73% above the amounts present in the control feed because this level of excess was found as the maximum inclusion of EAA not having a significant (P > 0.05) effect on feed intake or growth compared to optimum levels in broiler chickens [56, 57]. The amount of each feed consumed was measured every 24 hours. After every measurement, feeders’ position was switched to avoid errors due to side preference. Preference (%) for a feed was determined as a percentage of the experimental feed consumed divided by the total amount of both feeds consumed.
Animal and tissue sampling
The 48 lightest and 48 heaviest birds at 3 weeks of age from a batch of 580 male (feather-sexed) Ross 308 chickens (Darwalla Group, Qld, Australia) were selected as SG or FG chickens, respectively. Chickens showing any sign of pathology were not included. The birds were transferred to individual cages similar to those used for the DC experiments and offered a commercial diet until 6 weeks old. Individual cages allowed for individual feed intake and body weight recording. The five lightest and five heaviest birds were selected for the metabolic studies and euthanised using cervical dislocation. Cervical dislocation is a common method of euthanasia used in poultry research experiments where the bird’s neck is quickly pulled by hand leading to dislocation of the first cervical vertebrae with severing of the spinal cord and carotid arteries. This method of euthanasia does not involve the use of any chemical or equipment. Every bird was confirmed to be male using post-mortem and DNA sexing (MDS Australia Pty Ltd., Queanbeyan, Australia). A sample (400-500 mg) of the proventriculus taken from each bird was placed in 1 ml RNA-later solution (ThermoFisher Scientific, Waltham, USA) and stored at − 80 °C.
RNA sample preparation and sequencing
RNA was extracted from 30 to 40 mg of tissue using the Maxwell® 16 LEV Simply RNA Purification Kit manufactured by Promega Corporation, Madison, USA. The quality and quantity of the RNA samples were examined using a NanoDrop™ spectrophotometer (ThermoFisher Scientific, Waltham, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, USA). RNA samples with an Integrity Number ≥ 8, a 260/280 ratio between 2 and 2.1, and a 260/230 ratio between 2 and 2.2 were used for sequencing. An RNA sample of 500-4000 ng from each SG and FG bird was sequenced using high throughput sequencing (Queensland Brain Institute, University of Queensland, Australia). A TruSeq Stranded mRNA Library Prep Kit was used for library construction. Sequencing and coverage were based on 2 × 125 bp Paired-End Dual indexed reads and v4 illumina SBS chemistry (Illumina, San Diego, USA), with HiSeq 2000 used to produce an average of ~ 27 million sequence reads per sample.
The RNA-seq data was validated using RT-qPCR analysis. The detailed method is provided in the Additional file 1.
Protein sample preparation
Approximately 20 mg of tissue was used for protein extraction and Mass Spectrometry. The tissue was placed in ~ 300 μl of lysis buffer (6 M Guanidine Chloride, 50 mM Tris pH 8, and 10 mM DTT) and sonicated for 10 seconds before being vortexed at 30 °C for 1 hour. Acrylamide (25 mM) was added and incubated for 1 h before adding 5 mM DTT. Protein was precipitated by adding 4 volumes of 1:1 methanol acetone and stored overnight at − 20 °C before centrifuged to remove the solvent, and the pellet suspended in 0.1% SDS (Sodium dodecyl sulfate). The protein concentration of each sample was measured using NanoDrop (ThermoFisher Scientific, Waltham, USA). Samples were digested with trypsin overnight at 37 °C using an Amicon column and ammonium bicarbonate (50 mM). Samples were desalted by C-18 Zip-tip (adapted from Millipore procedure) and analyzed by liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS).
Mass spectrometry and data analysis
LC-MS/MS analysis was performed using a Prominence nanoLC system (Shimadzu, Kyoto, Japan) and TripleToF 5600 mass spectrometer with a Nanospray III interface (SCIEX, Toronto, Canada). A 70 min Liquid chromatography (LC) gradient was used to separate the peptides. The database of Protein Pilot (Uniprot, www.uniprot.org) was used to identify peptides and proteins. Sequential window acquisition of all theoretical mass spectra (SWATH) was performed on all samples. An information-dependent acquisition library (IDA) was prepared from one randomly chosen slow grower sample and one randomly chosen fast grower sample for each tissue.
Transcriptomic and genomic analysis
For transcriptomic analysis, sequenced reads were first put through a quality check (QC) using FastQC version 0.72 (Babraham Bioinformatics, Babraham Institute, Cambridge, UK) in Galaxy Australia platform with default parameters . Then, the reads were mapped to Galgal.5 reference assembly using aligners HISAT2 and Salmon with default settings [59, 60]. InteractiVenn, JMP, IGV, and MeV were used to visualize and perform cluster and Principal Component Analyses (PCA) [61, 62].
As for genomics, Galaxy Australia was used to perform the analyses . RNA sequences were mapped on the reference genome (Galgal.5) using Bowtie2, version 126.96.36.199 . Genomic variants between the 5 SG and 5 FG chickens were identified using FreeBayes, which is a Bayesian genetic variant detector designed to find small polymorphisms, specifically SNPs (single-nucleotide polymorphisms), indels (insertions and deletions), MNPs (multi-nucleotide polymorphisms), and complex events (composite insertion and substitution events). The list of candidate genes used for the study of genomic variants is presented in Table S2 in Additional file 1. The FreeBayes output was first filtered for Phred score ≥ 30 to attain accuracy of 99.9%. Phred score is an index of variant quality score indicating the confidence level in detecting a particular variant accurately . Next, the FreeBayes output was filtered for SNPs considered of interest in relations to metabolic functions of chickens growing at different rates. These SNPs were further filtered for a SIFT (Sorting Intolerant From Tolerant) score of < 0.05, which predicts the extent to which the polymorphism effects gene function, with 0.0 being deleterious and 1.0 being tolerable in Variant Effect Predictor (VEP) analysis.
Functional enrichment analysis
Differentially expressed genes (DEG) and differentially abundant proteins (DAP) between the FG and SG groups were used as input for metabolic pathway enrichment analyses. Enrichment analyses were performed in DAVID 6.8 . Enriched metabolic pathways and terms in different databases including Gene Ontology, GO , Kyoto Encyclopedia of Genes and Genomes, KEGG , and REACTOME  were used to attain insight into the function of DEGs and DAPs.
A 2 × 3 factorial design was used including two growth rates for SG or FG chickens and three preference comparisons of feed vs feed plus amino acids (Table 4). Preference for each feed in a DC test was compared to 50%, which is the neutral value, using t-student test. Main effects and interactions were analysed using the general linear model (GLM) procedure of SAS9.4 (SAS Institute, Cary, North Carolina, United States), with P < 0.05 as the level of significance in Tukey test. The statistical model to compare the preference and individual feed intakes was:
yijk = observation k in level i of group of chicken and level j of treatment; μ = the overall mean; Ai = the effect of level i of group of chicken (slow vs fast grower); Bj = the effect of level j of treatment; (AB)ij = the effect of the interaction of level i of chicken group with level j of treatment; εijk = random error; i = Number of levels of chicken group; j = Levels of treatment; k = biological replicates.
Total feed intake was used as a covariate for comparisons between treatments and controls and for amino acids consumed. The covariate allowed comparison between SG and FG chickens while adjusting the treatment effect for the variability in the feed intake.
Randomized design with a covariate was also run in GLM procedure of SAS9.4 according to the statistical model below.
yijk = observation k in treatment i group j; β0 = the intercept; β1 = the regression coefficient; xijk = a continuous independent variable of total feed intake (covariate); τi = the fixed effect of treatment (T1–3); gj = the fix effect of group (slow vs fast grower); εijk = random error.
For transcriptomics, read count and differential expression analysis were performed using Limma . Limma provides a linear model capable of multiple RNA comparisons in the context of multifactor design experiments to identify differentially expressed genes (DEG). The basic statistics used for DEG identification was t-statistics moderated across genes using a Bayesian model. This will increase the reliability of the results. Significance threshold was set as a P-value of 0.05 or less.
As for proteomics, sequential window acquisition of all theoretical mass spectra (SWATH)-MS relative quantitative data were analysed using PeakView v2.1 (SCIEX) . To identify differentially abundant proteins (DAP), statistical analyses were performed using MSstats in R as previously described [71, 72]. Values less than 0.05, adjusted for multiple testing, were considered as significant P-values.
Pathway enrichment analyses were performed on both DEG and DAP lists. Pathways that were statistically (P < 0.05) over-represented in the DEG and DEP lists were defined as significantly enriched. This was done through Fuzzy clustering algorithms and Fisher’s exact test built in DAVID bioinformatics platform .
Regarding genomic variants, allelic frequencies were compared between SG and FG chickens using Fisher’s exact test in SAS 9.2 (SAS Institute, Cary, North Carolina, United States). P-values of 0.05 or less were considered as significant for all allelic frequency comparisons.
Availability of data and materials
The raw datasets generated and analysed during the current study are available in the National Centre for Biotechnology Information (NCBI), Sequence Read Archive (SRA) under accession number PRJNA771098 (https://submit.ncbi.nlm.nih.gov/).
Conditionally essential amino acid
Differentially abundant proteins
Differentially expressed genes
Double choice test
Essential amino acid
Food and Agriculture Organization
Glucose transporter 1
Glucose transporter 2
Glucose transporter 4
Glucose transporter 10
Glucose transporter 11
Glucose transporter 13
Glucose transporter 4 Regulator
Glutamic-oxaloacetic transaminase 1
Insertions and deletions
Liquid chromatography electrospray ionization tandem mass spectrometry
Lipoic acid synthetase
Non-essential amino acid
Oxoglutarate (α-ketoglutarate) dehydrogenase
Pyruvate dehydrogenase beta
Pyruvate dehydrogenase complex component
Succinate dehydrogenase complex subunit A
Succinate dehydrogenase complex subunit B
Solute Carrier gene family
Solute carrier 1 type A2
Solute carrier 2 type A1
Solute carrier 2 type A4 Regulator
Solute carrier 2 type A10
Solute carrier 2 type A11
Solute carrier 38 type A1
Sequential window acquisition of all theoretical mass spectra
Umami taste receptor
Triosephosphate isomerase 1
Variant Effect Predictor
Tyrosine 3-monooxygenase/ tryptophan 5-monooxygenase activation protein epsilon
Tyrosine 3-monooxygenase/ tryptophan 5-monooxygenase activation protein theta
Tyrosine 3-monooxygenase/ tryptophan 5-monooxygenase activation protein zeta
Rubin CJ, Zody MC, Eriksson J, Meadows JRS, Sherwood E, et al. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010;464:587–91.
Swennen Q, Janssens GPJ, Decuypere E, Buyse J. Effects of substitution between fat and protein on feed intake and its regulatory mechanisms in broiler chickens: energy and protein metabolism and diet-induced thermogenesis. Poult Sci. 2004;83:1997–2004.
Farkhoy M, Modirsanei M, Ghavidel O, Sadegh M, Jafarnejad S. Evaluation of protein concentration and limiting amino acids including lysine and met + Cys in prestarter diet on performance of broilers. Vet Med Int. 2012;2012:394189.
Hofman P, Siegert W, Kenez A, Naranjo VD, Rodehutscord M. Very low crude protein and varying glycine concentrations in the diet affect growth performance, characteristics of nitrogen excretion, and the blood metabolome of broiler chickens. J Nutr. 2019;149:1122–32.
Oshel KM, Knight JB, Cao KT, Thai MV, Olson AL. Identification of a 30-base pair regulatory element and novel DNA binding protein that regulates the human GLUT4 promoter in transgenic mice. J Biol Chem. 2000;275:23666–73.
Reeds PJ, Burrin DG, Stoll B, van Goudoever JB. Role of the gut in the amino-acid economy of the host. In: Furst P, Young V, editors. Peptides and amino acids in enteral nutrition, vol. 3. Basel: Programme; 2000. p. 25–40.
Mack S, Bercovici D, De Groote G, Leclercq B, Lippens M, Pack M, et al. Ideal amino acid profile and dietary lysine specification for broiler chickens of 20 to 40 days of age. Brit Poult Sci. 1999;40:257–65.
Choi M, Chang CY, Clough T, Broudy D, Killeen T, MacLean B, et al. MSstats: an R package for statistical analysis of quantitative mass spectrometry-based proteomic experiments. Bioinformatics. 2014;30:2524–6.
SN did the conceptualization, data curation, formal analysis, writing manuscript, project administration; MRSF helped with formal analysis, methodology, software, revising manuscript; SC helped with experimental protocol and execution, formal analysis, revising manuscript; JLB did the editing of the manuscript, data curation; ER did the conceptualization, data interpretation, editing manuscript, funding and other resources, project management. The author(s) read and approved the final manuscript.
All the experiments and methods including raising, handling, euthanizing, and tissue sampling from the chickens were carried out according to the regulations and guidelines of 2013 edition of ‘Australian code for the care and use of animals for scientific purposes’ approved by the Australian Research Council. This study was carried out in compliance with ARRIVE Essential 10 guidelines 2.0. All the protocols and practices were approved by the Institutional Animal Care and Use Committee (ICLAS) of the University of Queensland, Australia (Approval number: CNAFS/253/14/RIRDC).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
Niknafs, S., Fortes, M.R.S., Cho, S. et al. Alanine-specific appetite in slow growing chickens is associated with impaired glucose transport and TCA cycle.
BMC Genomics23, 393 (2022). https://doi.org/10.1186/s12864-022-08625-2