RESEARCH ARTICLE

Genome-wide association study for the free amino acid and nucleotide components of breast meat in an F2 crossbred chicken population

Minjun Kim1https://orcid.org/0000-0002-8173-8431, Eunjin Cho2https://orcid.org/0000-0003-4800-1603, Jean Pierre Munyaneza1https://orcid.org/0000-0003-2521-3339, Thisarani Kalhari Ediriweera2https://orcid.org/0000-0001-9649-6925, Jihye Cha3https://orcid.org/0000-0002-9705-2979, Daehyeok Jin4https://orcid.org/0000-0001-5091-4271, Sunghyun Cho5https://orcid.org/0000-0002-2657-9159, Jun Heon Lee1,2,*https://orcid.org/0000-0003-3996-9209
Author Information & Copyright
1Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea
2Department of Bio-AI Convergence, Chungnam National University, Daejeon 34134, Korea
3Animal Genome & Bioinformatics, National Institute of Animal Science, Rural Development Administration, Wanju 55365, Korea
4Animal Genetic Resources Research Center, National Institute of Animal Science, Rural Development Administration, Hamyang 50000, Korea
5Research and Development Center, Insilicogen Inc., Yongin 19654, Korea
*Corresponding author: Jun Heon Lee, Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea. Tel: +82-42-821-7031, E-mail: junheon@cnu.ac.kr

© Copyright 2022 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Sep 15, 2022; Revised: Oct 26, 2022; Accepted: Oct 31, 2022

Published Online: Jan 31, 2023

Abstract

Flavor is an important sensory trait of chicken meat. The free amino acid (FAA) and nucleotide (NT) components of meat are major factors affecting meat flavor during the cooking process. As a genetic approach to improve meat flavor, we performed a genome-wide association study (GWAS) to identify the potential candidate genes related to the FAA and NT components of chicken breast meat. Measurements of FAA and NT components were recorded at the age of 10 weeks from 764 and 767 birds, respectively, using a White leghorn and Yeonsan ogye crossbred F2 chicken population. For genotyping, we used 60K Illumina single-nucleotide polymorphism (SNP) chips. We found a total of nine significant SNPs for five FAA traits (arginine, glycine, lysine, threonine content, and the essential FAAs and one NT trait (inosine content), and six significant genomic regions were identified, including three regions shared among the essential FAAs, arginine, and inosine content traits. A list of potential candidate genes in significant genomic regions was detected, including the KCNRG, KCNIP4, HOXA3, THSD7B, and MMUT genes. The essential FAAs had significant gene regions the same as arginine. The genes related to arginine content were involved in nitric oxide metabolism, while the inosine content was possibly affected by insulin activity. Moreover, the threonine content could be related to methylmalonyl-CoA mutase. The genes and SNPs identified in this study might be useful markers in chicken selection and breeding for chicken meat flavor.

Keywords: Chicken; Breast meat; Free amino acids; Nucleotides; Genome-wide association study

INTRODUCTION

Yeonsan ogye is a chicken breed native to Korea. The nutritional quality, taste, overall acceptance, and texture of the meat of the Korean native chicken breed, including the content of taste-related compounds, are superior to those of commercial broilers [1,2]; thus, consumers have a higher preference for it. Therefore, it would be a good strategy to crossbreed native chickens to produce broilers with tastier meat by taking advantage of the unique meat quality characteristics of native chickens. A study of the characteristics of the Yeonsan ogye chicken meat reported higher crude fat content and shearing power in Yeonsan ogye compared to other Korean native chickens, and the color of the meat was expressed as black [3]. In addition, it has been reported that the breast meat extract of Korean native chickens can alleviate inflammation [4].

To improve the palatability of chicken products, many studies have been conducted to determine the factors that influence meat quality. Studies have identified several parameters that influence meat quality, including tenderness, water-holding capacity, and meat color, in relation to the physical composition, while pH, dry matter, moisture, and intramuscular fat content are related to the chemical composition of the meat [58]. Taste-active precursors in meat, such as free amino acid (FAA) and nucleotide (NT), undergo chemical reactions during cooking, which determines many aspects of flavor, one of the main sensory features of meat [911].

As reported in [9], glycine, alanine, serine, threonine, and proline in the form of FAAs provide the sweet taste to the meat after cooking. Regarding NTs, inosine 5’-monophosphate, alone or in combination with other flavor compounds, provides an umami taste, greatly affecting the flavor of meat [12].

Many studies have been conducted to identify the genes related to the physicochemical properties of chickens. In [13], a study of the association of the fatty acid composition in chicken was conducted, and in [14], a gene related to the lipid profile of chicken muscle was reported. In [15,16], association studies were conducted on body composition. However, there have been no attempts to identify the genes related to FAA and NT content, which determine the flavor of chicken, in the whole genome. Furthermore, the traits present in the chicken Quantitative Trait Loci Database do not include the FAA and NT contents of breast meat [17]. Therefore, the objective of this study was to identify the potential candidate genes related to different FAA and NT contents in chicken breast through a whole-genome association analysis.

MATERIALS AND METHODS

Experimental animals

For the construction of the experimental chicken population, a reciprocal crossbred F2 generation was produced using the Yeonsan ogye breed and the White Leghorn breed owned by the National Institute of Animal Science (NIAS, Korea). A total of 810 birds from the F2 generation were made up of two lines according to the sex combination of F0 birds and were raised until 10 weeks of age with the mixing of males and females. All F2 birds were raised on a farm at the NIAS Center for Livestock Genetic Resources under the same environmental and nutritional conditions. Overall experimental outline and the construction of the F2 generation is illustrated in Fig. 1.

jast-65-1-57-g1
Fig. 1. Outline of the experimental workflow and the construction of the F2 generation chicken population. (A) Diagram of overall experimental contents. (B) The construction of the reciprocal crossbred F2 generation chickens using Yeonsan ogye and White Leghorn. The male bird of Yeonsan ogye and the female bird of White Leghorn produced F1 generation birds. And F1 birds mated between F1 siblings so that a part of F2 birds was generated. Another part of F2 birds was made up of the F2 progeny of a different sex combination of Yeonsan ogye and White Leghorn chickens as the F0 generation. FAA, free amino acid content; NT, nucleotide content; SNP, single nucleotide polymorphism; GWAS, genome-wide association study; GRM, genetic relationship matrix; GCTA, genome-wide complex trait analysis.
Download Original Figure
Phenotype measurement and quality control

At 10 weeks of age, all birds were slaughtered on the same day at a commercial abattoir. Viscera were removed from the chicken carcass immediately after slaughter, rapidly frozen at −35°C, and then stored at −20°C at NIAS (Wanju, Korea) on the same day. Frozen samples were stored for 2–6 weeks, and then transported to Chungnam National University (Deajeon, Korea) and stored at −80°C until deboning. Next, samples were thawed at 4°C for 20 h, the breast meat was separated from the carcass, and then the breast meat was divided into samples for the FAA and NT component analyses. The FAA analysis was performed using a liquid chromatography system (UltiMate 3000RS, Thermo Fisher Scientific, Waltham, MA, USA) to quantify the content of 21 FAAs from a sample (Table 1). From the quantified FAA results, the essential FAAs and non-essential FAAs contents of chickens were added and each phenotype was used in an association analysis. For the NT content analysis, high-performance liquid chromatography (1200 Series, Agilent, Santa Clara, CA, USA) was applied. The resulting phenotypes were used in the association test as traits for the content of each substance and are summarized in Table 1.

Table 1. Descriptive statistics of the measured phenotypes in the F2 chicken population
Trait Type unit Mean SD CV (%) Min Max Trans
Alanine Nonessential µmol/g 1.5772 0.7439 0.4716 0.032 5.139 Sqrt
Arginine Essential µmol/g 9.0896 5.9595 0.6556 0.059 42.62 Sqrt
Asparagine Nonessential µmol/g 0.2946 0.158 0.5363 0 0.834 Sqrt
Aspartic acid Nonessential µmol/g 0.4188 0.2583 0.6169 0.012 1.93 Sqrt
Cysteine Nonessential µmol/g 3.4547 1.9151 0.5544 0.09 12.61 Sqrt
GABA Nonessential µmol/g 0.0564 0.0632 1.1207 0 0.68 Sqrt
Glutamic acid Nonessential µmol/g 1.2466 0.661 0.5302 0.031 9.389 Sqrt
Glutamine Nonessential µmol/g 1.1161 0.6518 0.584 0.02 5.257 Sqrt
Glycine Nonessential µmol/g 0.8388 0.5757 0.6863 0.01 4.91 Sqrt
Histidine Essential µmol/g 0.34 0.1819 0.5349 0.01 1.07 Sqrt
Isoleucine Essential µmol/g 0.3973 0.1789 0.4503 0.008 1.12 Min-Max
Leucine Essential µmol/g 0.8782 0.4067 0.4631 0.019 2.38 No trans
Lysine Essential µmol/g 0.0299 0.1128 3.7769 0 1.44 Sqrt
Methionine Essential µmol/g 0.6394 0.5242 0.8198 0.013 4.169 Log
Phenylalanine Essential µmol/g 0.4195 0.1998 0.4762 0.006 1.58 Min-Max
Proline Nonessential µmol/g 0.3564 0.2792 0.7834 0.004 6.29 Sqrt
Serine Nonessential µmol/g 1.342 0.5624 0.419 0.039 3.26 Min-Max
Threonine Essential µmol/g 0.8727 0.4369 0.5006 0.014 3.007 Sqrt
Tryptophan Essential µmol/g 0.0343 0.0376 1.0964 0 0.89 Sqrt
Tyrosine Nonessential µmol/g 0.6248 0.5362 0.8582 0.019 2.97 Log
Valine Essential µmol/g 0.3802 0.2194 0.5771 0.002 1.707 Sqrt
Essential FAAs µmol/g 15.2358 7.4925 0.4918 0.32 52.59 Sqrt
Nonessential FAAs µmol/g 9.108 3.7094 0.4073 0.74 25.98 Sqrt
AMP NT mg/100g 10.2654 4.4131 0.4299 0.03 30.61 Sqrt
IMP NT mg/100g 297.7067 86.4317 0.2903 24.4 572.61 No trans
Inosine NT mg/100g 49.8131 20.7431 0.4164 4.6749 146.47 Sqrt
Hypoxanthine NT mg/100g 15.7305 5.7408 0.3649 4.74 42.52 Log

CV, coefficient of variation; Min, minimum; Max, maximum; Trans, transformation; Sqrt, square root scaling; GABA, γ-aminobutyric acid; Min-Max, minimum-maximum scaling; Log, Log scaling; FAA, free amino acid content; AMP, adenosine monophosphate; NT, nucleotide content; IMP, inosine 5’-monophosphate.

Download Excel Table

To normalize the phenotypic data, a transformation was carried out to accurately estimate the effect of single-nucleotide polymorphisms (SNPs) in the association test. For each phenotype in the dataset, a Shapiro–Wilk normality test was performed using R software with the “shapiro.test()” function, and a transformation step was applied to the phenotype with a p-value of < 0.05. To find the optimal transformation method, the dataset with the highest p-value was adopted after performing normality tests on the raw data and data transformed by log, square root, square, and minimum/maximum (min/max) methods (Table 1). Then, to remove outliers, values of three standard deviations or more from the mean of each data point were treated as missing data, and all subjects whose analysis results were identified as missing data were excluded from the analysis. The FAA and NT component datasets were filtered in the quality control process, leaving 764 and 767 birds for the later analysis of FAA and NT components, respectively.

Genotype data generation and quality control

For genomic DNA (gDNA) extraction, blood collection was performed from the brachial veins of birds at 8 weeks of age. gDNA was extracted from the collected blood using the Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA), and the quality and concentration of the DNA were checked using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific), with samples stored at −20°C. A total of 57,636 SNP genotypes were produced from the extracted gDNA using the Illumina 60K Chicken SNP BeadChip. The SNPs were removed if any of the following conditions were not met: call rate > 0.9, minor allele frequency > 0.01, Hardy–Weinberg equilibrium test p-value of > 0.01 × 10−6. After filtering, 29,175 SNPs were available for the subsequent analysis. PLINK 1.9 software was used for genotype quality control procedure [18].

The genome-wide association study (GWAS) and identification of the candidate genes

A genetic relationship matrix (GRM) was calculated using the genome-wide complex trait analysis (GCTA) tool[19]. An association analysis was performed using GCTA’s mixed linear model leave-one-chromosome-out statistical method, and the following model was used:

y = a + b x + g + e

where, y is a vector of the phenotype value, including FAAs and NTs; a is the mean value of phenotypic record; b is the additive effect of the SNP to be tested for association; x is the genotype of the SNP to be tested; g is the effect of all SNPs captured by the GRM, except those on the chromosome where SNP testing is located; and e is the residual effect vector. To accurately estimate the marker effect, sex, experimental batch, half-sibling family, and body weight were used as covariates. After the association analysis, the significance of the SNP effect was evaluated by a Bonferroni corrected p-value of 5%: p < 0.05/number of SNPs use for the association test. The GWAS result was visualized with a Manhattan plot using R software.

Based on the chicken reference genome (GRCg6a), we assumed that the genomic region around 1 Mbp was the significant region centered on the position of the most significant SNP for each chromosome, and known genes included in that region were reported as candidate genes. Location-based gene information was obtained by searching the Ensembl (https://asia.ensembl.org) and NCBI (https://www.ncbi.nlm.nih.gov) databases.

RESULTS

Basic statistics of phenotypic values

The descriptive statistics of the phenotype measurements of FAAs and NTs from the F2 chicken population are presented in Table 1. For these phenotypic values, basic statistical analyses were performed after transformation and outlier trimming. Histograms and quantile–quantile plots are shown in Fig. 2. Lysine, for which most values were measured as zero, was excluded from the test because its distribution was judged to be inappropriate for the association test. It was determined that the cause of this inappropriate data was an error in phenotype measurement.

jast-65-1-57-g2
Fig. 2. Phenotype distribution of free amino acids and nucleotide contents used for the association test. Histogram and quantile–quantile plot contents, respectively: (A) essential free amino acids, (B) arginine, (C) glycine, (D) threonine, (E) inosine, (F) lysine. Note that most values of lysine are zero in (F).
Download Original Figure
GWAS results and related genomic regions

Table 2 shows the significant SNPs for each trait as a result of the association test, with results presented as summary statistics of the GWAS. The SNP markers were located at chicken (Gallus gallus) chromosome 1 (GGA1), GGA2, GGA3, GGA4, and GGA7. Of the three SNPs associated with the essential FAAs, Gga_rs15943371 and Gga_rs14162297 were located at bp 33,096,495 and 33,432,575 on GGA2, respectively, and Gga_rs16609168 was located at bp 31,241,647 on GGA7. In arginine, Gga_rs14162297 did not reach the significance level, but it was the second most significant SNP among the markers located on GGA2, and Gga_rs16609168 was also the SNP with the highest significance in GGA7. Another SNP related to arginine was Gga_rs13971906, which was located at bp 170,582,107 on GGA1, and was significant for inosine. As other SNPs associated with inosine, GGaluGA265714, GGaluGA265712, and GGaluGA265969 were located on GGA4 at bp 74,314,272, 74,285,145, and 75,827,200, respectively. The SNP associated with glycine was GGaluGA028030 and was located at bp 81,831,988 on GGA1. The SNP associated with threonine was Gga_rs14412837 and was located at bp 108,537,571 on GGA3.

Table 2. Genome-wide significant single-nucleotide polymorphisms for the tested traits
Traits SNP ID Chr Position (bp) Allele 1 Allele 2 MAF SNP effect p-value
Essential FAAs Gga_rs15943371 2 33,096,495 T C 0.38 0.25 6.08 × 10−7
Essential FAAs Gga_rs14162297 2 33,432,575 T C 0.42 0.24 7.12 × 10−7
Essential FAAs Gga_rs16609168 7 31,241,647 C T 0.26 0.29 8.09 × 10−7
Arginine Gga_rs16609168 7 31,241,647 C T 0.26 0.32 2.16 × 10−8
Arginine Gga_rs13971906 1 170,582,107 A G 0.36 0.24 1.54 × 10−6
Inosine GGaluGA265714 4 74,314,272 G A 0.39 −0.32 1.63 × 10−8
Inosine GGaluGA265712 4 74,285,145 A G 0.39 −0.32 2.08 × 10−8
Inosine GGaluGA265969 4 75,827,200 T C 0.30 −0.34 4.08 × 10−8
Inosine Gga_rs13971906 1 170,582,107 A G 0.36 0.30 1.27 × 10−6
Glycine GGaluGA028030 1 81,831,988 T C 0.06 0.16 9.37 × 10−8
Threonine Gga_rs14412837 3 108,537,571 G A 0.50 −0.06 3.47 × 10−7

Chr, chromosome; bp, base pair; Allele 1, minor allele; Allele2, major allele; MAF, minor allele frequency; SNP, single nucleotide polymorphism; FAA, free amino acids.

Download Excel Table

The Manhattan plot in Fig. 3 presents only the traits with significant SNPs. The most significant SNP on GGA1 was rs13971906 for the essential FAAs, arginine, and inosine. Gga_rs14162297, which had the second highest significance for the essential FAAs, also had the second highest significance for arginine. Table 3 shows the significantly related genomic regions for each trait and the genes included therein.

jast-65-1-57-g3
Fig. 3. Genome-wide Manhattan plot of the −Log10 (p-value) for significant single-nucleotide polymorphisms. (A) the essential free amino acids, (B) arginine, (C) inosine, (D) glycine, and (E) threonine components. The x-axis shows the chromosome number and the y-axis indicates the −Log10 (p-value). The red horizontal line indicates a Bonferroni-corrected significance level of 5%. Circles with the same color indicate the same genomic region.
Download Original Figure
Table 3. Genes located in the genomic regions for the significant single-nucleotide polymorphisms
Trait Chr Start position (bp) End position (bp) Genes
Arginine, Inosine 1 170,065,610 170,700,171 RB1, LPAR6, RCBTB2, CYSLTR2, FNDC3A, MLNR, CDADC1, CAB39L, SETDB2, PHF11, RCBTB1, ARL11, KPNA3, SPRYD7, TRIM13, KCNRG
Arginine, Essential FAAs 2 32,588,887 33,940,447 HOXA3, HOXA4, HOXA5, HOXA6, HOXA7, HOXA9, HOXA10, HOXA11, HOXA13, EVX1, HIBADH, TAX1BP1, JAZF1, CREB5, TRIL, CHN2, WIPF3, SCRN1, FKBP14, PLEKHA8
Arginine, Essential FAAs 7 30,761,616 31,457,806 CXCR4, THSD7B, HNMT, SPOPL, NXPH2
Inosine 4 74,767,885 75,920,712 KCNIP4, ADGRA3, PPARGC1A, LDB2, QDPR, CLRN2, LAP3, MED28, FAM184B, LCORL, NCAPG
Glycine 1 81,317,906 82,292,366 CHD1L, FAAHL, TMEM39A, B4GALT4, B4GALTL, B3GAT1L, UPK1B, IGSF11, LSAMP
Threonine 3 108,221,189 108,846,159 TFAP2B, TFAP2D, CRISP3, CRISP2, RHAG, CYP2AC1, CYP2AC2, CENPQ, MMUT, OPN5L2, CDC5L

Chr, chromosome; bp, base pair; FAA, free amino acid.

Download Excel Table

DISCUSSION

The shared significant genomic region between the arginine and the essential free amino acids content

In this study, the candidate genes related to FAA and NT contents of chicken breast meat were identified through the GWAS analysis. In the GWAS result, the trend of significance for the essential FAAs was similar to that of arginine (Fig. 3). The reason why the two traits shared a significant genomic region in GGA1, GGA2, and GGA7 was that arginine content accounted for the largest portion of the essential FAAs in the phenotype. The arginine content in the breast meat of the F2 crossbred chicken was greater than the content of all other FAAs when the mean values were compared.

Candidate genes in the shared significant region of the arginine and inosine contents

From the GWAS results of inosine and arginine, it was observed that the genomic region on GGA1 was commonly significantly related (Table 2). One of the genes included in this region, the KCNRG gene, encoded a voltage-gated potassium channel regulator protein. The voltage-gated potassium channel can regulate the nitric oxide synthesis pathway, which is a process in which L-arginine is converted to nitric oxide. This might have an effect on the arginine content in the body [20,21]. In addition, the voltage-gated potassium channel can modulate the insulin sensitivity of human skeletal muscle, which modulates the nitric oxide pathway and affects L-arginine transport and metabolism. Therefore, the potassium channel may be related to the amount of L-arginine in muscle [2224]. Interestingly, if the content of adenosine triphosphate (ATP), a glucose metabolite in the muscle, changes as insulin manages skeletal muscle glucose uptake, the content of inosine, which is a repeatedly degraded form of ATP, may also be affected [25,26]. For the same reason, it could be concluded that the KCNIP4 gene on the significant region in GGA4 encodes a potassium channel interacting protein.

Diet and nutritional level in birds are factors affecting the chemical composition of muscles. Among the candidate genes related to dry matter intake in cattle reported in Javaid et al. [27] and Brunes et al. [28], the RB1, LPAR6, RCBTB2, CYSLTR2, FNDC3A, MLNR, PHF11, and KPNA3 genes were identical to the genes in the genomic region commonly related to inosine and arginine in this study. Therefore, it is possible that the difference in nutritional condition according to the feed intake of chickens raised under ad libitum conditions in this study affected the chemical composition of breast meat [27].

Candidate genes related to arginine and threonine

It has been reported that the arginine content can be increased by reducing nitric oxide synthase when the expression of the HOXA3 gene among the GGA2 genes related to arginine content in Table 3 is regulated [29]. It has also been confirmed that mutations in the THSD7B gene of GGA7 affect nitric oxide metabolism [30]. In the GWAS result of the threonine content, the MMUT gene encoding the methylmalonyl-CoA mutase located on GGA3 was included as a candidate gene. In animals, including chickens, threonine is catabolized to alpha-ketobutyrate by serine/threonine dehydratase, then converted to propionyl-CoA and carboxylated to methylmalonyl-CoA [3133]. The catalyzed product is finally converted to succinyl-CoA by methylmalonyl-CoA mutase and used in the tricarboxylic acid cycle [34]. Therefore, threonine content may be affected by methylmalonyl-CoA mutase involved pathway which use threonine as a substrate.

Further improvements from the present study

The molecules with significantly related SNP markers in this study were arginine, inosine, threonine, and glycine. As taste-active precursors that determine meat flavor during cooking, arginine and inosine give a bitter taste, and threonine and glycine give a sweet taste [9]. This study had limited success in finding a significant gene region related to the contents of FAA and NT molecules that give umami and sour tastes in chicken breast. The reason that significant SNPs were not found through GWAS in most traits is thought to be due to the large impact of environmental effects on the FAA and NT contents. The FAA content increases as protease hydrolyzes meat protein, and this degradation proceeds even in frozen meat [35]. A significant change in the amino acid profile of meat according to the storage duration at freezing temperature has also been reported [36]. In addition, in the process of thawing frozen meat, several NT-metabolizing enzymes perform a chain degradation to change the NT content [3739]. Based on the results of this study, it could be concluded that the period of storage under freezing conditions before analysis after the slaughter of chickens was long enough to result in significant changes in the FAA and NT contents. Future studies may capture more genetic effects when association tests are conducted by performing phenotyping on meat with a shorter storage time than in the present study.

CONCLUSION

Changes in the concentrations of FAAs and NTs in chicken are the main factors that change the flavor. To our knowledge, this was the first GWAS to explore genetic markers related to the content of FAAs and NTs in chicken. Through the GWAS, nine significant SNPs related to essential FAAs, arginine, glycine, threonine, and inosine components were found, and six significant genomic regions and candidate genes included therein were reported. These genomic regions were mapped on chromosome 1, 4, and 7 of chicken, and the potential candidate genes were KCNRG for arginine and inosine, HOXA3 and THSD7B for arginine, KCNIP4 for inosine, and MMUT for threonine. The reported candidate genes might be important genetic markers in chicken selection and breeding to help meet consumer expectations for healthier, tastier, and flavored meat.

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ0157852022)” Rural Development Administration, Korea and also supported by Chungnam National University.

Acknowledgements

Not applicable.

Availability of data and material

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Authors’ contributions

Conceptualization: Cha J, Jin D, Lee JH.

Data curation: Kim M, Cho E, Cha J, Jin D, Cho S.

Formal analysis: Kim M, Cho E.

Methodology: Kim M, Cho E, Cho S.

Software: Kim M, Cho E.

Validation: Cho E, Cho S, Lee JH.

Investigation: Kim M, Cho E, Munyaneza JP, Ediriweera TK.

Writing - original draft: Kim M.

Writing - review & editing: Kim M, Cho E, Munyaneza JP, Ediriweera TK, Cha J, Jin D, Cho S, Lee JH.

Ethics approval and consent to participate

Experimental use of animal was approved by Animal Ethics Committee of Chungnam National University (202103A-CNU-061).

REFERENCES

1.

Choe JH, Nam K, Jung S, Kim BN, Yun HJ, Jo CR. Differences in the quality characteristics between commercial Korean native chickens and broilers. Korean J Food Sci Anim Resour. 2010; 30:13-9

2.

Jung YK, Jeon HJ, Jung S, Choe JH, Lee JH, Heo KN, et al. Comparison of quality traits of thigh meat from Korean native chickens and broilers. Korean J Food Sci Anim Resour. 2011; 31:684-92

3.

Lee SG, Utama DT, Baek KH, Park YH, Han JY, Lee SK. Comparison of physicochemical characteristics of the meat in four lines of Korean native chickens. Korean J Poult Sci. 2015; 42:335-45

4.

Lee KC, Leem KH, Kim MG, Kim HK. Comparison of chemical composition and immune-enhancing activity of the four lines of Korean native chickens. Korean J Poult Sci. 2016; 43:135-42

5.

Lee Y, Saha A, Xiong R, Owens CM, Meullenet JF. Changes in broiler breast fillet tenderness, water-holding capacity, and color attributes during long-term frozen storage. J Food Sci. 2008; 73:E162-8

6.

Kirkpinar F, Ünlü HB, Serdaroğlu M, Turp GY. Effects of dietary oregano and garlic essential oils on carcass characteristics, meat composition, colour, pH and sensory quality of broiler meat. Br Poult Sci. 2014; 55:157-66

7.

Park MN, Hong EC, Kang BS, Kim HK, Kim JH, Na SH, et al. Chemical composition and meat quality of crossbred Korean native chickens (KNC). Korean J Poult Sci. 2010; 37:415-21

8.

Tougan PU, Dahouda M, Salifou CFA, Ahounou SGA, Kpodekon MT, Mensah GA, et al. Conversion of chicken muscle to meat and factors affecting chicken meat quality: a review. Int J Agron Agric Res. 2013; 3:1-20

9.

Dashdorj D, Amna T, Hwang I. Influence of specific taste-active components on meat flavor as affected by intrinsic and extrinsic factors: an overview. Eur Food Res Technol. 2015; 241:157-71

10.

Mottram D. Some aspects of the chemistry of meat flavour.In In: Shahidi F, editor.editor Flavor of meat and meat products. Boston, MA: Springer. 1994; p p. 210-30

11.

Liu Y, Xu X, Zhou G. Changes in taste compounds of duck during processing. Food Chem. 2007; 102:22-6

12.

Jung S, Bae YS, Kim HJ, Jayasena DD, Lee JH, Park HB, et al. Carnosine, anserine, creatine, and inosine 5′-monophosphate contents in breast and thigh meats from 5 lines of Korean native chicken. Poult Sci. 2013; 92:3275-82

13.

Jin S, Lee SH, Lee DH, Manjula P, Lee SH, Lee JH. Genetic association of DEGS1, ELOVL6, FABP3, FABP4, FASN and SCD genes with fatty acid composition in breast and thigh muscles of Korean native chicken. Anim Genet. 2020; 51:344-5

14.

d’André Hirwa C, Yan W, Wallace P, Nie Q, Luo C, Li H, et al. Effects of the thyroid hormone responsive spot 14α gene on chicken growth and fat traits. Poult Sci. 2010; 89:1981-91

15.

Liu R, Sun Y, Zhao G, Wang F, Wu D, Zheng M, et al. Genome-wide association study identifies loci and candidate genes for body composition and meat quality traits in Beijing-You chickens. PLOS ONE. 2013; 8e61172

16.

Nie Q, Fang M, Xie L, Peng X, Xu H, Luo C, et al. Molecular characterization of the ghrelin and ghrelin receptor genes and effects on fat deposition in chicken and duck. J Biomed Biotechnol. 2009; 2009:567120

17.

Hu ZL, Park CA, Reecy JM. Bringing the animal QTLdb and CorrDB into the future: meeting new challenges and providing updated services. Nucleic Acids Res. 2022; 50:D956-61

18.

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559-75

19.

Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011; 88:76-82

20.

Pérez-García MT, Cidad P, López-López JR. The secret life of ion channels: Kv1.3 potassium channels and proliferation. Am J Physiol Cell Physiol. 2018; 314:C27-42

21.

Senbel AM, Abd Elmoneim HM, Sharabi FM, Mohy El-Din MM. Neuronal voltage gated potassium channels may modulate nitric oxide synthesis in corpus cavernosum. Front Pharmacol. 2017; 8:297

22.

Xu J, Wang P, Li Y, Li G, Kaczmarek LK, Wu Y, et al. The voltage-gated potassium channel Kv1.3 regulates peripheral insulin sensitivity. Proc Natl Acad Sci. 2004; 101:3112-7

23.

Petrie JR, Ueda S, Webb DJ, Elliott HL, Connell JMC. Endothelial nitric oxide production and insulin sensitivity: a physiological link with implications for pathogenesis of cardiovascular disease. Circulation. 1996; 93:1331-3

24.

Rajapakse NW, Chong AL, Zhang WZ, Kaye DM. Insulin-mediated activation of the L-arginine nitric oxide pathway in man, and its impairment in diabetes. PLOS ONE. 2013; 8e61840

25.

Uemoto Y, Ohtake T, Sasago N, Takeda M, Abe T, Sakuma H, et al. Effect of two non- synonymous ecto-5′-nucleotidase variants on the genetic architecture of inosine 5′-monophosphate (IMP) and its degradation products in Japanese black beef. BMC Genom. 2017; 18:874

26.

Liu W, Zhao J. Insights into the molecular mechanism of glucose metabolism regulation under stress in chicken skeletal muscle tissues. Saudi J Biol Sci. 2014; 21:197-203

27.

Javaid S, Javid A, Farooq U, Kiran U, Akmal T. Variations in meat chemical composition of some captive avian species. J Rural Dev Agric. 2017; 2:57-65

28.

Brunes LC, Baldi F, Lopes FB, Lôbo RB, Espigolan R, Costa MFO, et al. Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle. J Anim Breed Genet. 2021; 138:23-44

29.

Al Sadoun H, Burgess M, Hentges KE, Mace KA. Enforced expression of Hoxa3 inhibits classical and promotes alternative activation of macrophages in vitro and in vivo. J Immunol. 2016; 197:872-84

30.

Yao Z, Lin A, Yi Y, Shen W, Zhang J, Luo P. THSD7B Mutation induces platinum resistance in small cell lung cancer patients. Drug Des Devel Ther. 2022; 16:1679-95

31.

Kidd MT, Kerr BJ. L-threonine for poultry: a review. J Appl Poult Res. 1996; 5:358-67

32.

Macelline SP, Chrystal PV, Liu SY, Selle PH. Implications of elevated threonine plasma concentrations in the development of reduced-crude protein diets for broiler chickens. Anim Prod Sci. 2021; 61:1442-8

33.

Tang Q, Tan P, Ma N, Ma X. Physiological functions of threonine in animals: beyond nutrition metabolism. Nutrients. 2021; 13:2592

34.

Shane B. Folate and vitamin B12 function.In In: Lennarz WJ, Lane MD, editors.editors Encyclopedia of biological chemistry. 2nd ed London: Academic Press. 2013; p p. 324-8

35.

Baranenko D, Kolodyaznaya V, Broyko Y. Effect of cold treatment on the amino acid composition of veal. Agron Res. 2014; 12:705-16

36.

Al-Sabagh ESH, Taha NM, Saleh EA, El-Far AH, Sadek KM. Effect of freezing and frozen storage on amino acid profile and fatty acid pattern in imported and local meat. Alex J Vet Sci. 2016; 49:113-21

37.

Qi J, Xu Y, Zhang W, Xie X, Xiong G, Xu X. Short-term frozen storage of raw chicken meat improves its flavor traits upon stewing. LWT. 2021; 142:111029

38.

Zhang T, Lu H, Wang L, Yin M, Yang L. Specific expression pattern of IMP metabolism related-genes in chicken muscle between cage and free range conditions. PLOS ONE. 2018; 13e0201736

39.

Hong H, Regenstein JM, Luo Y. The importance of ATP-related compounds for the freshness and flavor of post-mortem fish and shellfish muscle: a review. Crit Rev Food Sci Nutr. 2017; 57:1787-98