Next Article in Journal
Effect of a Low-Methane Diet on Performance and Microbiome in Lactating Dairy Cows Accounting for Individual Pre-Trial Methane Emissions
Next Article in Special Issue
Comparison of the Automatic and Manual Broiler Pre-Slaughter Chain Based on Trailer Microclimate during Transportation and Its Effect on m. pectoralis major
Previous Article in Journal
Neurohormonal Profiles of Assistance Dogs Compared to Pet Dogs: What Is the Impact of Different Lifestyles?
Previous Article in Special Issue
Comparison of Quality and Sensory Characteristics of Spent Hen and Broiler in South Korea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 9
10.3390/ani11092596
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine

by
Satoshi Kubota
*,
Kasarat Promkhun
,
Panpradub Sinpru
,
Chanadda Suwanvichanee
,
Wittawat Molee
and
Amonrat Molee
School of Animal Technology and Innovation, Institute of Agricultural Technology, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
*
Author to whom correspondence should be addressed.
Animals 2021, 11(9), 2596; https://doi.org/10.3390/ani11092596
Submission received: 29 July 2021 / Revised: 28 August 2021 / Accepted: 31 August 2021 / Published: 3 September 2021
(This article belongs to the Special Issue Poultry Production: Factors Affecting Meat Quality)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Carnosine is a bioactive food component with several potential health benefits for humans due to its physiological functions. Dietary supplementation with β-alanine or L-histidine can increase the carnosine content of skeletal muscles in chickens. Dietary supplementation with β-alanine or L-histidine has produced a slow-growing chicken variety with high carnosine content in the breast meat; however, the supplementation with L-histidine alone softens the meat toughness, which may affect consumers’ willingness to buy the meat. Gene expression is a key factor that influences meat quality. Understanding the molecular mechanisms that affect carnosine content and meat toughness would allow the production of more value-added slow-growing chickens. We compared global gene expression in chicken breast muscles with differing carnosine contents and meat toughness produced through dietary supplementation with β-alanine or L-histidine. We identified differentially expressed genes involved in regulating myosin, collagen, intramuscular fat, and calpain—factors that may affect meat tenderness. Pathway enrichment analysis indicated that the insulin-related and adipocytokine signaling pathways were altered by dietary supplementation with β-alanine or L-histidine. These data will be useful for future studies on carnosine content and meat toughness in slow-growing chickens.

Abstract

Korat chicken (KRC) is a slow-growing chicken bred in Thailand, whose meat exhibits a unique toughness. A previous study produced KRC breast meat containing high carnosine content through dietary supplementation with β-alanine or L-histidine; however, the KRC that were fed an L-histidine-supplemented diet produced meat that was significantly more tender. Herein, we performed RNA-Seq to identify candidate genes involved in the regulation of carnosine content and meat toughness. Total RNA was isolated from five female KRC breast muscles in each treatment group that KRC fed diets without supplementation, supplemented with β-alanine or L-histidine. Compared to the non-supplemented group, we identified 118 and 198 differentially expressed genes (DEGs) in the β-alanine or L-histidine supplementation groups, respectively. Genes potentially related to meat tenderness—i.e., those regulating myosin, collagen, intramuscular fat, and calpain—were upregulated (LOC107051274, ACSBG1, and CAPNS2) and downregulated (MYO7B, MYBPH, SERPINH1, and PGAM1). However, carnosine synthase gene was not identified. Functional enrichment analysis identified pathways affected by dietary supplementation, including the insulin signaling pathway (β-alanine supplementation) and the insulin resistance and adipocytokine signaling pathways (L-histidine supplementation). The FoxO signaling pathway was identified as a regulatory network for both supplementation groups. The identified genes can be used as molecular markers of meat tenderness in slow-growing chickens.

1. Introduction

Korat chicken (KRC) is produced by crossing the males of Thai indigenous chickens with the female Suranaree University of Technology (SUT) breeder line, and has a slow growth rate. KRC meat contains lower levels of fat and higher quantities of collagen compared to the meat from a commercial broiler. It also has a firm and chewy texture, which attracts domestic consumers to a greater extent than broiler meat in Thailand [1]. However, female KRCs are not as preferred by farmers because they grow at a slower rate than males, but have the same selling price. Thus, adding more value to female KRCs is desirable. Intarapichet and Maikhunthod [2] reported that meat from Thai indigenous female chickens contains higher levels of carnosine than meat from males. Carnosine is well-known for its physiological functions, and intake of carnosine has been shown to contribute to human health in many age-associated diseases and conditions [3,4]. Carnosine concentration in the skeletal muscles of broiler chickens can be enhanced by dietary supplementation with the constituent amino acids of carnosine (i.e., β-alanine and L-histidine) [5,6]. Suwanvichanee et al. [7] reported that compared to KRCs fed a diet without supplementation, female KRCs fed a diet supplemented with β-alanine or L-histidine have significantly higher carnosine content in breast muscles (by 26.42 and 32.76%, respectively). Growth and slaughtering performances were not significantly different between dietary treatments; however, meat toughness (as measured by shear force value) was significantly reduced only by supplementation with L-histidine [7]. These findings suggest that female KRCs with increased carnosine content produced through supplemented diets can increase the consumers’ demand as a functional food, which motivates farmers’ willingness to raise them. However, dietary supplementation with L-histidine also reduces the toughness of the meat, which is a unique property of KRC meat.
Among the traits that strongly affect the quality of meat and consumers’ willingness to purchase poultry products, tenderness is the most important attribute [8]. Many factors contribute to the wide variation in tenderness of poultry meat, including bird strain, age, sex, environment, antemortem and/or postmortem handling conditions, and cooking method [9]. Several studies have identified the genetic basis of tenderness—including genes encoding the proteins of the calpain system (CAPN1)—in broiler chickens [10]. Using a whole transcriptomic analysis, Piórkowska et al. [11] identified biological processes and genes that are potentially associated with factors affecting chicken breast meat tenderness, such as degradation of filamins (ASB2), lipogenesis (PLIN1 and THRSP), and collagen synthesis (COL16A1, COL20A1, LEPREL4, PCOLCE2, P4HA3, and VWA1). These findings suggest that such candidate genes can be used as markers to determine meat tenderness; however, it is unclear how the actions of these genes are affected in chickens fed a supplemented diet to increase carnosine content.
Transcriptome sequencing technology (RNA-Seq) can help elucidate the genetic and physiological processes that regulate the phenotype of particular quantitative traits [12]. Dietary supplementation has been shown to alter the transcriptomic profile of breast muscles in chickens, and RNA-Seq can be used to identify key genes and molecular regulatory mechanisms related to factors affecting meat quality, such as color and taste [13], and melanin deposition [14]. In this study, we performed RNA-Seq in female KRC breast muscle to identify differentially expressed genes and underlying molecular mechanisms affected by the diets supplemented with β-alanine or L-histidine. To the best of our knowledge, this is the first study on global gene expression in chicken breast muscle with differing carnosine contents and meat toughness produced through dietary supplementation. The identified genes can be used as a molecular marker associated with carnosine content and meat tenderness in slow-growing chickens.

2. Materials and Methods

2.1. Animals and Sample Collection

All animal procedures applied in this research were approved by the Ethics Committee on Animal Use of the Suranaree University of Technology (user application ID: U1-08780-2563).
In this study, we used the 15 breast muscles of female KRCs collected in a previous study by Suwanvichanee et al. [7]. Briefly, a total of 300 one-day-old female KRCs were generated by crossbreeding male Thai indigenous chickens and the female SUT breeder line. All chicks were reared in open-sided housing from 0 to 3 weeks of age. KRC feeding management is designed as three different phases based on the age of the chicken: 0 to 3, 4 to 6, and 7 to 10 weeks. The basal diet at the initial, second, and final phases contains 21, 19, and 17% of protein, respectively. Birds were raised on a basal diet at the initial phase before starting the experiments. Chickens were fed an experimental diet containing β-alanine (Sigma-Aldrich, St. Louis, MO, USA, 146064) or L-histidine (AppliChem GmbH, Darmstadt, Germany, A3738) from 4 to 10 weeks of age to increase carnosine content. Group A was the control without supplements. The other groups were fed by diets supplemented with 1% (1 g/100 g diet) β-alanine (B) or 0.5% (0.5 g/100 g diet) L-histidine (C). Each treatment group contained five replicates (cages) with 20 chickens per replicate. All chickens were given feeds and water ad libitum during the experiment. The composition of the experimental diet is presented in Table S1. At the end of the feeding trial, seven chickens were removed from each cage. Of these seven chickens, five were fasted and euthanized for carcass and meat quality measurements and the other two chickens were sacrificed for carnosine content measurements. For transcriptome analysis, a small amount of the left breast muscle tissue was extracted from similar parts of the two birds used for the carnosine content measurement. The muscle tissue was collected in tubes, frozen in liquid nitrogen immediately, and kept at −80 °C until the usage for transcriptome analysis.

2.2. RNA Extraction

The breast muscle from one of the two KRCs was subjected to RNA extraction. Total RNA was extracted from the 15 chickens—i.e., one chicken from each replicate of each treatment group—using TRIzol Reagent (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s protocol. The purity of the extracted RNA was analyzed using 1% agarose gel electrophoresis and measuring the absorbance using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). The RNA integrity number (RIN) was verified using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and 1 µg of total RNA with RIN ≥ 6.7 was used for subsequent transcriptome analysis.

2.3. Library Preparation and Sequencing

Library preparation, sequencing, and analysis were performed by Vishuo Biomedical (Thailand) Ltd. according to the manufacturer’s instructions. Poly (A) mRNA was isolated using the Poly (A) mRNA Magnetic Isolation Module (NEB, County Road Ipswich, MA, USA) or the VAHTS mRNA-seq V3 Library Prep Kit for Illumina (Vazyme, Nanjing, China). A first-strand synthesis reaction buffer and random primers were applied for mRNA fragmentation and priming. First-strand and second-strand cDNAs were synthesized respectively using ProtoScript II Reverse Transcriptase (NEB) or Second Strand Synthesis Enzyme Mix (NEB). The double-stranded cDNA was then purified using beads. Following DNA repair and dA-tailing using the End Prep Enzyme Mix (NEB), an adaptor was added to both ends by T-A ligation. Adaptor-ligated DNA was then selected using beads, and fragments of ~400 bp with an insert size of approximately 300 bp were recovered. The size-selected DNA was amplified by bridge PCR using the P5 and P7 primers that contain an index for multiplexing. The PCR products were purified using beads, validated, and quantified using a Qsep100 system (BiOptic, New Taipei City, Taiwan) or Qubit 3.0 fluorometer (Invitrogen, Carlsbad, CA, USA), respectively. Libraries were analyzed using an Illumina Novaseq 6000 instrument (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Sequencing was performed using a 2 × 150 bp paired-end configuration. Base calling was performed using the NovaSeq software.

2.4. Sequencing Quality Assessment and Differential Expression Analysis

RNA-Seq reads were obtained from the breast muscles of KRC fed the experimental diets. Low-quality and technical sequences such as adapters and primers were removed from the raw sequencing reads using Cutadapt [15] version 1.9.1, and high-quality data (clean reads) were generated for further analysis. Hisat2 [16] version 2.0.1 was used to align the clean reads against the chicken reference genome, GRCg6a (GenBank Assembly ID: GCA_000002315.5). Gene expression levels were estimated based upon fragments per kilobase of transcript per million fragments mapped (FPKM) using HTSeq [17] version 0.6.1. The DESeq2 Bioconductor package [18] was used for differential expression analysis between the control group and each group fed with dietary supplements. The Benjamini–Hochberg method was applied to control the false discovery rate. Transcripts with adjusted values of p < 0.05 and an absolute fold change of ≥1.5 were considered to be differentially expressed. The raw data conferred in this publication have been deposited in NCBI’s Gene Expression Omnibus which can be access through GEO Series accession number GSE179643 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE179643, accessed on 7 July 2021).

2.5. Gene Ontology and Pathway Enrichment Analysis

To identify the biological functions associated with the differentially expressed genes (DEGs), gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were conducted using the database for annotation, visualization, and integrated discovery (DAVID) [19] version 6.8 with default parameters (p < 0.10). DEGs found in the β-alanine or L-histidine supplementation groups were independently subjected to GO enrichment and KEGG pathway analysis. The GO database comprises three main branches: biological process (BP), cellular component (CC), and molecular function (MF).

2.6. Verification of RNA-Seq Results by Quantitative PCR (qPCR)

Differential expression results from RNA-Seq were confirmed by qPCR. The extracted RNA was reverse transcribed into cDNA using a High-Capacity cDNA reverse transcription kit (Thermo Fisher Scientific) following the manufacturer’s protocols. A total of five transcripts (i.e., FOXO1, IGF2BP3, PDE10A, GRIP2, and PPM1K) were randomly selected from DEGs that showed differences in both comparison groups, and primers were designed by applying NCBI Primer-BLAST on the website of NCBI (http://www.ncbi.nlm.nih.gov/tools/primer-blast/, accessed on 14 May 2021). β-Actin was used as an internal reference gene to normalize the expression data. The qPCR reactions were carried out on a LightCycler 480 System (Roche Diagnostics GmbH, Mannheim, Germany). The qPCR reactions had a total volume of 20 µL with 10 µL of SYBR Green I Master Mix (Roche Diagnostics GmbH), 1 µL of each primer (10 µM), 6 µL of nuclease-free water, and 2 µL of cDNA template. The thermocycling program consisted of initial denaturation at 95 °C for 5 min followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 61 °C for 30 s, and extension at 70 °C for 1 min for FOXO1, PDE10A, and PPM1K, or 40 cycles of denaturation at 95 °C for 30 s, annealing at 62 °C for 30 s, and extension at 72 °C for 1 min for GRIP2 and IGF2BP3. All qPCR reactions were run with three replicates for each sample. Melting curve analysis was performed following the qPCR reactions. Relative gene expression values were calculated using the 2−∆∆CT method [20]. The mean 2−∆∆CT value for each gene was transformed into a fold change value to compare the qPCR results with those of RNA-Seq. The primer sequences used are listed in Table S2.

3. Results

3.1. Quality of RNA-Seq Reads

An average of 46,769,860 raw sequencing reads were generated. After excluding low-quality reads, 15 high-quality samples were further mapped onto the chicken reference genome. The quality report for RNA-Seq showed that at least 94.31% of reads had a sequencing quality score exceeding Q30 (percentage of bases with a Phred value ≥ 30) in 15 libraries. The average ratio of high-quality reads to the reference genome was 82.43%. The results of quality assessment of sequencing data are presented in Table 1.

3.2. Significantly Differentially Expressed Transcripts

Followin g RNA-Seq of the breast muscle of female KRCs in the non-supplemented group (A), β-alanine or L-histidine supplementation groups (B and C), we compared global gene expression in the B and C groups to that in group A so as to identify the underlying molecular mechanisms affected by diet. A total of 118 and 198 DEGs were identified in the groups B and C, respectively, compared to the group A (Figure 1). Compared with group A, group B had 55 upregulated genes and 63 downregulated genes, and group C had 89 upregulated genes and 109 downregulated genes among these DEGs. All information regarding the DEGs is provided in Supplementary Materials (Tables S3 and S4).

3.3. GO Annotation and KEGG Pathways Analyses of DEGs

The GO functional analysis identified 14 and 16 GO terms in the β-alanine or L-histidine supplementation groups, respectively (p < 0.10; Table S5), of which six and eight GO terms were significantly enriched (p < 0.05) (Table 2 and Table 3). No GO terms categorized into MF and CC branches were detected in the β-alanine or L-histidine supplementation groups. The identified KEGG pathways of the DEGs are listed in Table 4. Compared to group A, the insulin signaling pathway (gga04910) was identified as being enriched in the β-alanine supplementation group, whereas the insulin resistance (gga04910) and adipocytokine signaling (gga04920) pathways were enriched in the L-histidine supplementation group. The forkhead box O (FoxO) signaling pathway (gga04068) was commonly identified as a regulatory network in both dietary supplement groups.

3.4. Validation of DEGs by qPCR

Of the detected DEGs, 77 and 157 DEGs were only identified in the β-alanine or L-histidine supplementation groups, respectively, compared to the non-supplemented group (Figure 2). Forty-one genes were found to be common in the two comparison groups. Five DEGs, including three upregulated genes (FOXO1, IGF2BP3, and PDE10A) and two downregulated genes (GRIP2 and PPM1K), were randomly selected from the 41 DEGs and their expression levels were verified by qPCR. The qPCR results showed that the expression trend of the DEGs was consistent in the RNA-Seq results in both comparison groups (Figure 3), suggesting that the RNA-Seq results in this study were accurate and reproducible.

4. Discussion

A previous study showed the breast meat from female KRCs fed diets without supplementation (A), supplemented with β-alanine (B) or L-histidine (C) contained 2756, 3484, and 3659 µg/g of carnosine, respectively, with shear force values of 3.31, 3.50, and 2.94 kg/g [7]. The carnosine content in group A and shear force value in group C were significantly lower than those in the other groups, respectively (p < 0.05) [7]. In this study, we performed RNA-Seq to identify differentially expressed genes and underlying molecular mechanisms affected by the diets supplemented with β-alanine or L-histidine. The quality report for RNA-Seq data obtained in this study (Table 1) indicated that the data were replicable and could be used in future analyses. In the differential expression analysis of transcripts (Figure 1, Tables S3 and S4), we identified some genes associated with myosin, collagen, intramuscular fat (IMF), and calpain, all of which may be related to meat quality.
Carnosine can be synthesized from β-alanine and L-histidine under in a reaction catalyzed by carnosine synthase. Increased skeletal carnosine concentrations in chickens by supplementation of dietary carnosine, β-alanine, or L-histidine have been reported elsewhere [5,6,21,22], and the upregulated expression of carnosine synthase corresponds with high carnosine content [6,21]. Qi et al. [6] reported that dietary supplementation with 250–2000 mg/kg β-alanine quadratically increased carnosine content, promoted mRNA expression of carnosine synthase 1 (CARNS1) and solute carrier family 6 member 6 (SLC6A6), and reduced the shear force of the breast muscle in broilers (p < 0.05). However, CARNS1 and SLC6A6 were not detected in this study. Carnosine content in chickens is strongly influenced by their breed [2]. Tian et al. [23] found that the carnosine content observed in the breast meat of Black-Bone Silky Fowl—a unique breed of chicken native to South China—was approximately 1.8-fold higher than that observed in the meat of White Plymouth Rock reared under the same conditions (p < 0.01). A recent study has reported significantly higher carnosine content in the breast muscles of Thai indigenous slow-growing chickens than in broiler chickens [24]. KRC is a crossbred chicken with Thai indigenous chickens and has a slow growth rate. The carnosine synthase gene was not detected as DEGs in this study, likely because slow-growing chickens have constantly high expression levels of this gene and are not affected by the supplemented diet. Downregulation of acidic amino acid decarboxylase GADL1 isoform X1 (GADL1), which is involved with β-alanine metabolism in chicken, was detected in the L-histidine supplementation group. GADL1 is considered to be an enzyme with multiple biological activities and is responsible for converting aspartic acid to β-alanine in carnosine biosynthesis [25]. In mice, GADL1 knockout mice were deficient in β-alanine and carnosine in the skeletal muscle [25]. However, the association of the GADL1 and carnosine contents has not been reported in chickens. It is still unclear how the downregulation of the GADL1 gene by L-histidine supplementation contributes to an increase in the carnosine contents of female KRC breast muscle.
Muscle structure and composition (i.e., muscle fibers, intramuscular connective tissue, and intramuscular fat) in livestock are related to meat sensory quality, including tenderness [26]. The characterization of muscle fibers by contractile and metabolic properties suggested that the myosin heavy chain isoforms (MyHCs) existing in thick filaments has a significant effect on the muscle’s contractile properties. MyHCs appear to be the most appropriate markers for muscle fiber-type delineation [27]. Indeed, upregulation of myosin heavy chains 1 (MYH1) and 15 (MYH15) have been detected in tough meat in previous studies on beef cattle and chicken, respectively [11,28]. Although myosin-related genes were found to be downregulated in this study—i.e., VIIB (MYO7B) in the β-alanine supplementation group and myosin binding protein H like (MYBPH) in the L-histidine supplementation group—no differentially-expressed MyHC gene was identified in this study. In chickens, myosin binding protein H has been speculated to be involved in the regulation of muscle contraction [29]; however, its function remains unclear. The composition and structure of connective tissue also affect meat tenderness [30], and collagen is considered the primary determinant of the shear force, particularly in cattle. Piórkowska et al. [11] found that a gene encoding collagen protein or involved in their synthesis (i.e., COL16A1, COL20A1, LEPREL4, PCOLCE2, P4HA3, and VWA1) was highly expressed in chickens with lower levels of tenderness. This suggests that components of extracellular matrices, such as collagen, glycoproteins, and proteoglycans, are involved in determining the tenderness of broiler chicken meat. In our study, two genes related to collagen—serpin family H member 1 (SERPINH1), and collagen alpha-2(IV) chain-like (LOC107051274)—were downregulated and upregulated, respectively, in the β-alanine supplementation group. SERPINH1 encodes Hsp47, a collagen-specific molecular chaperone. Hsp47 is essential for the molecular maturation of collagen, and Hsp47 expression strongly correlates with the expression of collagens in several types of cells and tissues [31]. Zhang et al. [32] reported upregulation of the SERPINH1 gene in male Qinchuan cattle that had a significantly higher shear force than females. Downregulation of the SERPINH1 was also detected in the present study; however, the shear force values did not differ significantly between groups supplemented and not supplemented with β-alanine. Although upregulation of the LOC107051274 gene was detected after β-alanine supplementation, the function of this gene has not been reported yet. Future studies are required to clarify how altering the expression of these genes affects collagen synthesis in chickens. A certain amount of IMF can enhance meat quality traits in chickens [33]. For instance, Wen et al. [34] reported significant differences in tenderness between pectoral muscles with < 0.5% IMF and those with >0.5% IMF in partridges (Alectoris chukar) (p < 0.01). Several studies have reported DEGs related to IMF content in chickens [33,35]. Of these genes, upregulation of the acyl-CoA synthetase bubblegum family member 1 (ACSBG1) gene was detected in the L-histidine supplementation group. Although this gene mediates the activation of long-chain fatty acids to synthesize cellular lipids [35], a study on IMF deposition is required to clarify the involvement of this gene in KRCs.
Shear force depends not only on the structural and metabolic properties of a muscle at the time of slaughter, but also on muscle metabolism during rigor mortis and aging [34,36]. The kinetics of postmortem muscle acidification are faster in poultry than in cattle and strongly affect tenderness of the meat [26]. Calpain is an intracellular Ca2+-dependent cysteine protease that is closely related to the tenderness of meat in livestock [37]. Okumura et al. [38] identified four calpain genes (i.e., CAPN1, CAPN2, CAPN3, and CAPN1.5) that were ubiquitously expressed in chickens. Among these, a single nucleotide polymorphism (SNP) in the CAPN1 gene was significantly related to meat tenderness in chickens and Japanese quail [10,39]. Moreover, Piórkowska et al. [40] found that higher expression of the CAPN1 gene leads to higher tenderness in breast muscles in fast-growing broiler chickens. The CAPNS1 and CAPNS2 genes encode the small subunits one and two of calpain, respectively [41], and upregulation of the CAPNS2 gene was detected in the β-alanine supplementation group in this study. Although the CAPNS1 protein is known to be important for maintaining the stability of both conventional catalytic subunits of calpain, its physiological role is yet unknown [41]. Recently, a significant association has been reported between SNPs in CAPNS2 and CAPN1 genes and meat tenderness in Iberian pigs, and the G allele in CAPN1_rs8135866G > A SNP was associated with tender meat [42]. These results suggest that the calpain enzyme is related to meat tenderness in chickens and the upregulation of the CAPNS2 gene by β-alanine supplementation may associate with meat toughness.
The GO analysis showed that compared to the non-supplemented group, the DEGs were significantly enriched in six and eight functions in the β-alanine or L-histidine supplementation groups, respectively (p < 0.05, Table 2 and Table 3). These functions did not overlap between the groups. Among the DEGs that may be related to meat tenderness (i.e., MYO7B, MYBPH, LOC107051274, SERPINH1, ACSBG1, and CAPNS2), SERPINH1 was identified in the “extracellular space” term of the CC category, along with another 11 DEGs (Table 2). Although this gene is associated with collagen, no BP terms involving collagen were detected in either experimental group. Listrat et al. [26] reported that collagen has a limited effect on the sensory quality of chicken meat because the animals are slaughtered at an early physiological stage when intramuscular collagen is not cross-linked significantly. El-Senousey et al. [43] reported that dietary alpha-lipoic acid significantly decreased collagen content and mRNA expression of the COL3A1 gene in broiler chickens, although it did not change meat tenderness. Katemala et al. [44] reported that the breast meat of 10-week-old KRC exhibited higher shear force values than that of commercial broiler chickens (p < 0.05), although there was no difference in the collagen content. These findings suggest that the higher tenderness detected in the L-histidine supplementation group was not caused by changes in the extracellular matrix (e.g., due to changes in collagen synthesis). Of the BP terms listed in Table 3, wound healing and gluconeogenesis were also identified among the top 10 enriched GO terms in a transcriptome study on the breast muscles of two broiler chicken lines that showed significant differences in meat quality (ultimate pH, muscle glycolytic potential, and toughness) [45]. In this study, we found that the L-histidine supplementation group showed downregulation of the 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 1 and phosphoglycerate mutase 1 (PGAM1) genes, upregulation of the cryptochrome circadian regulator 1 gene, and enrichment in the gluconeogenesis function. A high correlation has been reported between three proteins—pyruvate kinase, PGAM1, and triosephosphate isomerase 1—and the tenderness of breast meat in Thai indigenous chickens [46]. Moreover, significantly higher levels of these three proteins have also been observed in the tough meat. PGAM catalyzes the interconversion of 2-phosphoglycerate and 3-phosphoglycerate in the glycolytic cycle [47]. These findings suggest that the DEGs involved in gluconeogenesis are associated with meat characteristics and that the PGAM1 gene can be a functional marker for meat tenderness in KRC.
The KEGG pathway analysis comparing the non-supplemented and supplemented groups identified the insulin signaling pathway in the β-alanine-supplemented group and insulin resistance and adipocytokine signaling pathway in the L-histidine-supplemented group (p < 0.05, Table 4). In animals, adipogenesis occurs earliest in visceral fat deposits, followed by adipogenesis in subcutaneous, intermuscular, and intramuscular fat deposits [48]. Insulin affects a series of biological reactions in muscle and adipose tissues, including lipid metabolism [49]. By comparing DEGs between muscle and fat tissues in broiler chickens, Zhang et al. [50] identified lipid metabolism signaling pathways involved in the regulation of fat deposition (i.e., adipocytokine signaling pathway, insulin signaling pathway, and pyruvate metabolism). These findings suggest that the insulin-related and adipocytokine signaling pathways affected by the supplemented diet likely regulate fat deposition in intermuscular adipose tissue in the breast muscle of KRC. In this study, the FoxO signaling pathway was identified as a common regulatory network affected by the supplemented diet. FoxO acts as a regulator to maintain tissue homeostasis over time and affects many cellular processes involved in apoptosis, cell cycle regulation, DNA damage repair, metabolism, and oxidative stress [51,52]. The FoxO family consists of four members—i.e., FOXO1, FOXO3, FOXO4, and FOXO6—and upregulation of the FOXO1 gene was observed in this study. FoxOs have been previously identified as potential candidate genes regulating growth in chickens [53,54]. Willson et al. [55] investigated gene expression in the liver among meat birds, the layer strain, and their F1 crosses, and found that the FoxO signaling pathway—through its effects on cell regulation and altered metabolism—is associated with growth variations in chickens. Furthermore, Li et al. [56] studied miRNA and mRNA expression in slow-growing chickens at 6, 14, 22, and 30 weeks of age and found that FoxO signaling is a key pathway underlying the postnatal development of breast muscles in slow-growing chickens. Postnatal growth of skeletal muscle is achieved primarily by increasing the hypertrophy of existing muscle fibers [56]. This suggests that the FoxO signaling pathway is related to muscle fiber hypertrophy. FoxO transcription factors act as negative regulators of type I fiber-related gene expression in muscle atrophy [57]. The thigh muscles of broiler chickens fed a high-carnosine diet show a higher density of type I fiber and higher carnosine content along with a decrease in type I fiber diameter and meat toughness (p < 0.05) [21]. These findings suggest that the expression levels of FoxO transcription factors may correlate with muscle fiber size after diet supplementation with β-alanine or L-histidine. In broiler, increased skeletal muscle fiber diameter and area could increase shear force and reduce meat tenderness [58]. A significant association has been reported between the genotype of the CAPN1 gene and breast muscle fiber diameter in KRC, and KRCs carrying the BB genotype showed a significantly smaller muscle fiber diameter at 10 weeks of age [59]. Genetic selection of female KRCs for large muscle fiber size combined with feeding them diets supplemented with β-alanine or L-histidine may produce meat containing high carnosine concentrations without losing the toughness unique to this meat.

5. Conclusions

This study analyzed differences in global gene expression in the breast muscle of female KRCs with differing carnosine contents and meat toughness that were produced by dietary supplementation with β-alanine or L-histidine. Compared to the non-supplemented group, 118 and 198 DEGs were identified in the β-alanine or L-histidine supplementation groups, respectively. Although dietary supplementation with β-alanine or L-histidine produced breast meat with significantly higher carnosine content, genes previously reported to be associated with increasing carnosine content, such as CARNS1, were not activated. Functional enrichment analysis of the DEGs identified different underlying genetic mechanisms that were altered by each supplemented diet and likely regulated meat tenderness. The identified genes, including MYO7B, MYBPH, LOC107051274, SERPINH1, ACSBG1, CAPNS2, and PGAM1, may be potential functional markers for meat tenderness in slow-growing chickens. However, we did not examine the structural and metabolic properties of a muscle in this study, future studies are required to clarify how altering the expression of these genes affects meat tenderness.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11092596/s1, Table S1: Composition of the experimental diet, Table S2: List of primer sequences for quantitative PCR, Table S3: Differentially expressed genes in the β-alanine supplementation group compared to the non-supplemented group, Table S4: Differentially expressed genes in the L-histidine supplementation group compared to the non-supplemented group, Table S5: Gene ontology terms for the differentially expressed genes in the β-alanine or L-histidine supplementation groups compared to the non-supplemented group.

Author Contributions

Conceptualization, S.K., W.M. and A.M.; methodology and investigation, S.K., K.P., P.S. and C.S.; writing, S.K.; supervision, W.M. and A.M.; project administration, S.K.; funding acquisition, S.K. and A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Suranaree University of Technology, Thailand Science Research and Innovation (TSRI) grant number MRG6280120, and the Center of Excellence on Technology and Innovation for Korat chicken Business Development.

Institutional Review Board Statement

All animal procedures were approved by the Ethics Committee on Animal Use of the Suranaree University of Technology (user application ID: U1-08780-2563).

Data Availability Statement

Raw sequencing data are accessible through GEO Series accession number GSE179643 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE179643). All data generated or analyzed during this study available from the corresponding authors on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maliwan, P.; Khempaka, S.; Molee, W. Evaluation of various feeding programmes on growth performance, carcass and meat qualities of Thai indigenous crossbred chickens. S. Afr. J. Anim. Sci. 2017, 47, 16–25. [Google Scholar]
  2. Intarapichet, K.O.; Maikhunthod, B. Genotype and gender differences in carnosine extracts and antioxidant activities of chicken breast and thigh meats. Meat Sci. 2005, 71, 634–642. [Google Scholar]
  3. Boldyrev, A.A.; Aldini, G.; Derave, W. Physiology and pathophysiology of carnosine. Physiol. Rev. 2013, 93, 1803–1845. [Google Scholar]
  4. Hipkiss, A.R. Carnosine and its possible roles in nutrition and health. Adv. Food Nutr. Res. 2009, 57, 87–154. [Google Scholar]
  5. Kai, S.; Watanabe, G.; Kubota, M.; Kadowaki, M.; Fujimura, S. Effect of dietary histidine on contents of carnosine and anserine in muscles of broilers. Anim. Sci. J. 2015, 86, 541–546. [Google Scholar]
  6. Qi, B.; Wang, J.; Ma, Y.; Wu, S.; Qi, G.; Zhang, H. Effect of dietary β-alanine supplementation on growth performance, meat quality, carnosine content, and gene expression of carnosine-related enzymes in broilers. Poult. Sci. 2018, 97, 1220–1228. [Google Scholar]
  7. Suwanvichanee, C.; Kubota, S.; Molee, W.; Yongsawatdigul, J.; Molee, A. Effect of β-Alanine and L-Histidine supplementation on growth performance and meat quality in female Korat chicken. In Proceedings of the 2nd International Conference on Tropical Animal Science and Production and the 2nd International Conference on Native Chicken, Nakhon Ratchasima, Thailand, 9–12 July 2019; p. 163. [Google Scholar]
  8. Maltin, C.; Balcerzak, D.; Tilley, R.; Delday, M. Determinants of meat quality: Tenderness. Proc. Nutr. Soc. 2003, 62, 337–347. [Google Scholar]
  9. Ismail, I.; Joo, S.T. Poultry meat quality in relation to muscle growth and muscle fiber characteristics. Korean J. Food Sci. Anim. Resour. 2017, 37, 873–883. [Google Scholar]
  10. Shu, J.T.; Zhang, M.; Shan, Y.J.; Xu, W.J.; Chen, K.W.; Li, H.F. Analysis of the genetic effects of CAPN1 gene polymorphisms on chicken meat tenderness. Genet. Mol. Res. 2015, 14, 1393–1403. [Google Scholar]
  11. Piórkowska, K.; Żukowski, K.; Nowak, J.; Połtowicz, K.; Ropka-Molik, K.; Gurgul, A. Genome-wide RNA-Seq analysis of breast muscles of two broiler chicken groups differing in shear force. Anim. Genet. 2016, 47, 68–80. [Google Scholar]
  12. Malone, J.H.; Oliver, B. Microarrays, deep sequencing and the true measure of the transcriptome. BMC Biol. 2011, 9, 1–9. [Google Scholar]
  13. Li, W.; Zhang, X.; Du, J.; Li, Y.; Chen, Y.; Cao, Y. RNA-seq-based quantitative transcriptome analysis of meat color and taste from chickens administered by eucalyptus leaf polyphenols extract. J. Food Sci. 2020, 85, 1319–1327. [Google Scholar]
  14. Li, D.; Wang, X.; Fu, Y.; Zhang, C.; Cao, Y.; Wang, J.; Zhang, Y.; Li, Y.; Chen, Y.; Li, Z.; et al. Transcriptome analysis of the breast muscle of Xichuan black-bone chickens under tyrosine supplementation revealed the mechanism of tyrosine-induced melanin deposition. Front. Genet. 2019, 10, 457. [Google Scholar]
  15. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar]
  16. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar]
  17. Anders, S.; Pyl, P.T.; Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar]
  18. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar]
  19. Dennis, G.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4, 1–11. [Google Scholar]
  20. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)). Methods 2001, 25, 402–408. [Google Scholar]
  21. Cong, J.; Zhang, L.; Li, J.; Wang, S.; Gao, F.; Zhou, G. Effects of dietary supplementation with carnosine on growth performance, meat quality, antioxidant capacity and muscle fiber characteristics in broiler chickens. J. Sci. Food Agric. 2017, 97, 3733–3741. [Google Scholar]
  22. Kralik, G.; Sak-Bosnar, M.; Kralik, Z.; Galović, O. Effect of β-alanine dietary supplementation on concentration of carnosine and quality of broiler muscle tissue. J. Poult. Sci. 2014, 51, 151–156. [Google Scholar]
  23. Tian, Y.; Xie, M.; Wang, W.; Wu, H.; Fu, Z.; Lin, L. Determination of carnosine in Black-Bone Silky Fowl (Gallus gallus domesticus Brisson) and common chicken by HPLC. Eur. Food Res. Technol. 2007, 226, 311–314. [Google Scholar]
  24. Lengkidworraphiphat, P.; Wongpoomchai, R.; Taya, S.; Jaturasitha, S. Effect of genotypes on macronutrients and antioxidant capacity of chicken breast meat. Asian-Australas. J. Anim. Sci. 2020, 33, 1817–1823. [Google Scholar]
  25. Mahootchi, E.; Homaei, S.C.; Kleppe, R.; Winge, I.; Hegvik, T.A.; Megias-Perez, R.; Totland, C.; Mogavero, F.; Baumann, A.; Glennon, J.C.; et al. GADL1 is a multifunctional decarboxylase with tissue-specific roles in β-alanine and carnosine production. Sci. Adv. 2020, 6, eabb3713. [Google Scholar]
  26. Listrat, A.; Lebret, B.; Louveau, I.; Astruc, T.; Bonnet, M.; Lefaucheur, L.; Picard, B.; Bugeon, J. How muscle structure and composition influence meat and flesh quality. Sci. World J. 2016, 2016, 1–14. [Google Scholar]
  27. Pette, D.; Staron, R.S. Myosin isoforms, muscle fiber types, and transitions. Microsc. Res. Tech. 2000, 50, 500–509. [Google Scholar]
  28. Muniz, M.M.M.; Fonseca, L.F.S.; dos Santos Silva, D.B.; de Oliveira, H.R.; Baldi, F.; Chardulo, A.L.; Ferro, J.A.; Cánovas, A.; de Albuquerque, L.G. Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle. Meat Sci. 2021, 173, 108378. [Google Scholar]
  29. Vaughan, K.T.; Weber, F.E.; Einheber, S.; Fischman, D.A. Molecular cloning of chicken myosin-binding protein (MyBP) H (86-kDa protein) reveals extensive homology with MyBP-C (C-protein) with conserved immunoglobulin C2 and fibronectin type III motifs. J. Biol. Chem. 1993, 268, 3670–3676. [Google Scholar]
  30. Liu, A.; Nishimura, T.; Takahasi, K. Relationship between structural properties of intramuscular connective tissue and toughness of various chicken skeletal muscles. Meat Sci. 1996, 43, 43–49. [Google Scholar]
  31. Ito, S.; Nagata, K. Biology of Hsp47 (Serpin H1), a collagen-specific molecular chaperone. Semin. Cell Dev. Biol. 2017, 62, 142–151. [Google Scholar]
  32. Zhang, Y.Y.; Zan, L.S.; Wang, H.B. Screening candidate genes related to tenderness trait in Qinchuan cattle by genome array. Mol. Biol. Rep. 2011, 38, 2007–2014. [Google Scholar]
  33. Cui, H.X.; Liu, R.R.; Zhao, G.P.; Zheng, M.Q.; Chen, J.L.; Wen, J. Identification of differentially expressed genes and pathways for intramuscular fat deposition in pectoralis major tissues of fast- and slow-growing chickens. BMC Genom. 2012, 13, 213. [Google Scholar]
  34. Wen, Y.; Liu, H.; Liu, K.; Cao, H.; Mao, H.; Dong, X.; Yin, Z. Analysis of the physical meat quality in partridge (Alectoris chukar) and its relationship with intramuscular fat. Poult. Sci. 2020, 99, 1225–1231. [Google Scholar]
  35. Yang, S.H.; He, C.S.; Li, C.H.; Liu, G.Q. RNA-Seq reveals differentially expressed genes and pathways affecting intramuscular fat metabolism in Huangshan black chicken population. J. Agric. Sci. 2020, 12, 117–127. [Google Scholar]
  36. Girolami, A.; Marsico, I.; D'Andrea, G.; Braghieri, A.; Napolitano, F.; Cifuni, G.F. Fatty acid profile, cholesterol content and tenderness of ostrich meat as influenced by age at slaughter and muscle type. Meat Sci. 2003, 64, 309–315. [Google Scholar]
  37. Lian, T.; Wang, L.; Liu, Y. A new insight into the role of calpains in post-mortem meat tenderization in domestic animals: A review. Asian-Australas. J. Anim. Sci. 2013, 26, 443–454. [Google Scholar]
  38. Okumura, F.; Shimogiri, T.; Shinbo, Y.; Yoshizawa, K.; Kawabe, K.; Mannen, H.; Okamoto, S.; Cheng, H.H.; Maeda, Y. Linkage mapping of four chicken calpain genes. Anim. Sci. J. 2005, 76, 121–127. [Google Scholar]
  39. Rasouli, Z.; Zerehdaran, S.; Azari, M.A.; Shargh, M.S. Genetic polymorphism of the CAPN1 gene is associated with meat quality traits in Japanese quail. Br. Poult. Sci. 2013, 54, 171–175. [Google Scholar]
  40. Piórkowska, K.; Nowak, J.; Połtowicz, K. The normalisation of CAPN gene expression in M. pectoralis superficialis in broiler lines differing in growth rate and their relationship to breast muscle tenderness. Br. Poult. Sci. 2015, 56, 452–458. [Google Scholar]
  41. Sorimachi, H.; Hata, S.; Ono, Y. Impact of genetic insights into calpain biology. J. Biochem. 2011, 150, 23–37. [Google Scholar]
  42. Fernández-Barroso, M.Á.; Silió, L.; Rodríguez, C.; Palma-Granados, P.; López, A.; Caraballo, C.; Sánchez-Esquiliche, F.; Gómez-Carballar, F.; García-Casco, J.M.; Muñoz, M. Genetic parameter estimation and gene association analyses for meat quality traits in open-air free-range Iberian pigs. J. Anim. Breed Genet. 2020, 137, 581–598. [Google Scholar]
  43. El-Senousey, H.K.; Fouad, A.M.; Yao, J.H.; Zhang, Z.G.; Shen, Q.W. Dietary alpha lipoic acid improves body composition, meat quality and decreases collagen content in muscle of broiler chickens. Asian-Australas. J. Anim. Sci. 2013, 26, 394–400. [Google Scholar]
  44. Katemala, S.; Molee, A.; Thumanu, K.; Yongsawatdigul, J. Meat quality and Raman spectroscopic characterization of Korat hybrid chicken obtained from various rearing periods. Poult. Sci. 2021, 100, 1248–1261. [Google Scholar]
  45. Beauclercq, S.; Hennequet-Antier, C.; Praud, C.; Godet, E.; Collin, A.; Tesseraud, S.; Metayer-Coustard, S.; Bourin, M.; Moroldo, M.; Martins, F.; et al. Muscle transcriptome analysis reveals molecular pathways and biomarkers involved in extreme ultimate pH and meat defect occurrence in chicken. Sci. Rep. 2017, 7, 6447. [Google Scholar]
  46. Teltathum, T.; Mekchay, S. Relationships between Pectoralis muscle proteomes and shear force in Thai indigenous chicken meat. Agric. Nat. Resour. 2010, 44, 53–60. [Google Scholar]
  47. Doherty, M.K.; McLean, L.; Hayter, J.R.; Pratt, J.M.; Robertson, D.H.L.; EI-Shafei, A.; Gaskell, S.J.; Beynon, R.J. The proteome of chicken skeletal muscle: Changes in soluble protein expression during a growth in layer strain. Proteomics 2004, 4, 2082–2093. [Google Scholar]
  48. Hausman, G.J.; Dodson, M.V.; Ajuwon, K.; Azain, M.; Barnes, K.M.; Guan, L.L.; Jiang, Z.; Poulos, S.P.; Sainz, R.D.; Smith, S.; et al. Board-invited review: The biology and regulation of preadipocyte and adipocytes in meat animals. J. Anim. Sci. 2009, 87, 1218–1246. [Google Scholar]
  49. Dimitriadis, G.; Mitrou, P.; Lambadiari, V.; Maratou, E.; Raptis, S.A. Insulin effects in muscle and adipose tissue. Diabetes Res. Clin. Pract. 2011, 93, S52–S59. [Google Scholar]
  50. Zhang, M.; Zheng, D.; Peng, Z.; Zhu, Y.; Li, R.; Wu, Q.; Li, Y.; Li, H.; Xu, W.; Zhang, M.; et al. Identification of differentially expressed genes and lipid metabolism signaling pathways between muscle and fat tissues in broiler chickens. J. Poult. Sci. 2021, 58, 131–137. [Google Scholar]
  51. Huang, H.; Tindall, D.J. Dynamic FoxO transcription factors. J. Cell Sci. 2007, 120, 2479–2487. [Google Scholar]
  52. Eijkelenboom, A.; Burgering, B.M.T. FOXOs: Signalling integrators for homeostasis maintenance. Nat. Rev. Mol. Cell Biol. 2013, 14, 83–97. [Google Scholar]
  53. Xie, L.; Luo, C.; Zhang, C.; Zhang, R.; Tang, J.; Nie, Q.; Ma, L.; Hu, X.; Li, N.; Da, Y.; et al. Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits. PLoS ONE 2012, 7, e30910. [Google Scholar]
  54. Chen, B.; Xu, J.; He, X.; Xu, H.; Li, G.; Du, H.; Nie, Q.; Zhang, X. A genome-wide mRNA screen and functional analysis reveal FOXO3 as a candidate gene for chicken growth. PLoS ONE 2015, 10, e0137087. [Google Scholar]
  55. Willson, N.L.; Forder, R.E.A.; Tearle, R.; Williams, J.L.; Hughes, R.J.; Nattrass, G.S.; Hynd, P.I. Transcriptional analysis of liver from chickens with fast (meat bird), moderate (F1 layer x meat bird cross) and low (layer bird) growth potential. BMC Genom. 2018, 19, 309. [Google Scholar]
  56. Li, Y.; Chen, Y.; Jin, W.; Fu, S.; Li, D.; Zhang, Y.; Sun, G.; Jiang, R.; Han, R.; Li, Z.; et al. Analyses of microRNA and mRNA expression profiles reveal the crucial interaction networks and pathways for regulation of chicken breast muscle development. Front. Genet. 2019, 10, 197. [Google Scholar]
  57. Sanchez, A.M.J.; Candau, R.B.; Bernardi, H. FoxO transcription factors: Their roles in the maintenance of skeletal muscle homeostasis. Cell Mol. Life Sci. 2014, 71, 1657–1671. [Google Scholar]
  58. Felício, A.M.; Gaya, L.G.; Ferraz, J.B.S.; Moncau, C.T.; Mattos, E.C.; Santos, N.P.; Filho, T.M.; Balieiro, J.C.C.; Eler, J.P. Heritability and genetic correlation estimates for performance, meat quality and quantitative skeletal muscle fiber traits in broiler. Livest. Sci. 2013, 157, 81–87. [Google Scholar]
  59. Kubota, S.; Vandee, A.; Keawnakient, P.; Molee, W.; Yongsawatdigul, J.; Molee, A. Effects of the MC4R, CAPN1, and ADSL genes on body weight and purine content in slow-growing chickens. Poult. Sci. 2019, 98, 4327–4337. [Google Scholar]
Animals 11 02596 g001 550
Figure 1. Scatter plot of differentially expressed genes: (a) volcano plot of the β-alanine supplementation; (b) volcano plot of the L-histidine supplementation.
Figure 1. Scatter plot of differentially expressed genes: (a) volcano plot of the β-alanine supplementation; (b) volcano plot of the L-histidine supplementation.
Animals 11 02596 g001
Animals 11 02596 g002 550
Figure 2. Number of differentially expressed genes in the β-alanine or L-histidine supplementation group compared to the non-supplemented group. A, B, and C represent the non-supplemented group, the β-alanine supplementation group, and the L-histidine supplementation group, respectively.
Figure 2. Number of differentially expressed genes in the β-alanine or L-histidine supplementation group compared to the non-supplemented group. A, B, and C represent the non-supplemented group, the β-alanine supplementation group, and the L-histidine supplementation group, respectively.
Animals 11 02596 g002
Animals 11 02596 g003 550
Figure 3. Validation of five DEGs detected from RNA-Seq analysis by qPCR. The x-axis represents the gene name and the y-axis represents the mRNA expression levels expressed in a fold change value. The expression levels obtained by RNA-Seq and qPCR are represented by solid fill and pattern fill columns, respectively. A, B, and C represent the non-supplemented group, the β-alanine supplementation group, and the L-histidine supplementation group, respectively. Gene abbreviation represents forkhead box O1 (FOXO1), insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3), phosphodiesterase 10A (PDE10A), glutamate receptor interacting protein 2 (GRIP2), and protein phosphatase Mg2+/Mn2+ dependent 1K (PPM1K), respectively.
Figure 3. Validation of five DEGs detected from RNA-Seq analysis by qPCR. The x-axis represents the gene name and the y-axis represents the mRNA expression levels expressed in a fold change value. The expression levels obtained by RNA-Seq and qPCR are represented by solid fill and pattern fill columns, respectively. A, B, and C represent the non-supplemented group, the β-alanine supplementation group, and the L-histidine supplementation group, respectively. Gene abbreviation represents forkhead box O1 (FOXO1), insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3), phosphodiesterase 10A (PDE10A), glutamate receptor interacting protein 2 (GRIP2), and protein phosphatase Mg2+/Mn2+ dependent 1K (PPM1K), respectively.
Animals 11 02596 g003
Table
Table 1. RNA-Seq reads and mapping rate of breast muscle samples from 15 chickens fed an experimental diet.
Table 1. RNA-Seq reads and mapping rate of breast muscle samples from 15 chickens fed an experimental diet.
Sample ID 1Raw ReadsClean ReadsQ30 (%) 2GC
Content (%)		Total
Mapped Read		Uniquely Mapped ReadsMultiple Mapped ReadsMapping Rage (%)
A146,462,62046,288,55494.9455.0137,652,12733,451,2124,200,91581.34
A247,166,51847,013,35695.0454.0139,871,96835,898,9643,973,00484.81
A347,881,99447,713,85094.3156.6136,339,76632,824,8933,514,87376.16
A443,994,62443,820,89294.4956.2533,585,94228,618,8634,967,07976.64
A548,893,70448,746,18095.0254.8741,053,31935,969,2635,084,05684.21
B146,163,11646,008,44694.7654.5538,295,42334,456,4083,839,01583.24
B248,282,90648,121,86094.8155.0340,616,99535,710,6644,906,33184.4
B347,914,04447,776,80094.6355.3639,450,43636,438,5543,011,88282.57
B447,830,61447,678,70494.6654.139,195,22935,730,6023,464,62782.21
B546,867,32246,705,71094.6154.538,611,38935,737,4112,873,97882.67
C146,346,77846,172,55894.853.9638,610,45535,584,4663,025,98983.62
C246,252,97246,087,26494.8953.7839,318,31436,417,3622,900,95285.31
C347,149,78047,007,51094.3654.7738,372,18534,911,8273,460,35881.63
C448,618,94048,452,69095.054.4640,454,30436,297,7074,156,59783.49
C541,721,98041,582,68494.6753.5235,046,74132,462,3852,584,35684.28
Average46,769,860.846,611,803.894.7354.7138,431,639.534,700,705.43,730,934.182.43
1 Sample ID represents treatment group and biological replicate number. Treatment groups are A (non-supplemented), B (β-alanine supplementation), and C (L-histidine supplementation), respectively. A1, A2, A3, A4, and A5 indicate the five biological replicate breast muscle tissues in treatment group A, the rest of the groups (B and C) share the same name rules. 2 Q30 indicates the percentage of bases with a Phred value ≥ 30.
Table
Table 2. Significantly enriched gene ontology terms in the β-alanine supplementation group compared to the non-supplemented group.
Table 2. Significantly enriched gene ontology terms in the β-alanine supplementation group compared to the non-supplemented group.
CategoryTerm IDTermCountp-ValueGenes 1
Biological processGO:0003151Outflow tract morphogenesis46.75 × 10−4BMP4↑, HEYL↓, DHRS3↑, SFRP2↓
GO:0007160Cell-matrix adhesion30.031OTOA↑, TSC1↑, ITGB1BP1↑
GO:0002043Blood vessel endothelial cell proliferation involved in sprouting angiogenesis20.034BMP4↑, ITGB1BP1↑
Cellular componentGO:0009986Cell surface87.13 × 10−4ADAMTS15↓, RTN4RL1↓, SDC1↓, CFTR↑, VASN↓, PLA2R1↑, OTOA↑, THBD↓
GO:0005615Extracellular space120.001PXDNL↓, ADAMTS15↓, TST↓, BMP4↑, PLA2G15↑, VASN↓, TNFSF10↑, CPN1↓, SERPINH1↓, SFRP2↓, PPFIBP2↑, CTSV↑
GO:0016324Apical plasma membrane40.029CNTFR↓, CFTR↑, OTOA↑, AMOTL1↑
1 Genes with up and down arrows represent those upregulated or downregulated in the supplementation group, respectively.
Table
Table 3. Significantly enriched gene ontology terms in the L-histidine supplementation group compared to the non-supplemented group.
Table 3. Significantly enriched gene ontology terms in the L-histidine supplementation group compared to the non-supplemented group.
CategoryTerm IDTermCountp-ValueGenes 1
Biological processGO:0035556Intracellular signal Transduction105.68 × 10−4ASB8↓, ASB2↑, SGK1↑, PRKAA2↓, RGS9↑, ASB5↑, NUAK1↑, DAPK1↑, SPSB1↑, ARHGEF7↑
GO:0042752Regulation of circadian rhythm30.025PRKAA2↓, NOCT↓, PPARA↑
GO:0006094Gluconeogenesis30.032PFKFB4↓, PGAM1↓, CRY1↑
GO:0010629Negative regulation of gene expression40.038CDKN1A↓, CTGF↑, TIPARP↓, NOCT↓
GO:0009267Cellular response to Starvation30.039PRKAA2↓, MYOD1↑, PIK3C2B↓
GO:0042060Wound healing30.039PECAM1↑, PPARA↑, SLC11A1↓
GO:0032922Circadian regulation of gene expression30.045NOCT↓, PPARA↑, CRY1↑
Molecular functionGO:0004722Protein serine/threonine phosphatase activity30.029PPM1K↓, PPM1J↓, PDP1↓
1 Genes with up and down arrows represent those upregulated or downregulated in the supplementation group, respectively.
Table
Table 4. Significantly enriched KEGG pathway in the β-alanine or L-histidine supplementation groups compared to the non-supplemented group.
Table 4. Significantly enriched KEGG pathway in the β-alanine or L-histidine supplementation groups compared to the non-supplemented group.
TermCountp-ValueGenes 1
β-alanine supplementation
gga04068: FoxO signaling pathway50.008FOXO1↑, CCNG2↑, TNFSF10↑, G6PC3↓, IRS2↑
gga04910: Insulin signaling pathway40.044FOXO1↑, G6PC3↓, IRS2↑, TSC1↑
L-histidine supplementation
gga04068: FoxO signaling pathway60.008CDKN1A↓, SGK1↑, PRKAA2↓, FOXO1↑, TNFSF10↑, IRS2↑
gga04931: Insulin resistance50.018PRKAA2↓, FOXO1↑, PPARA↑, CREB5↑, IRS2↑
gga04920: Adipocytokine signaling pathway40.029ACSBG1↑, PRKAA2↓, PPARA↑, IRS2↑
1 Genes with up and down arrows represent that upregulated or downregulated in the supplementation group, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kubota, S.; Promkhun, K.; Sinpru, P.; Suwanvichanee, C.; Molee, W.; Molee, A. RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine. Animals 2021, 11, 2596. https://doi.org/10.3390/ani11092596

AMA Style

Kubota S, Promkhun K, Sinpru P, Suwanvichanee C, Molee W, Molee A. RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine. Animals. 2021; 11(9):2596. https://doi.org/10.3390/ani11092596

Chicago/Turabian Style

Kubota, Satoshi, Kasarat Promkhun, Panpradub Sinpru, Chanadda Suwanvichanee, Wittawat Molee, and Amonrat Molee. 2021. "RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine" Animals 11, no. 9: 2596. https://doi.org/10.3390/ani11092596

APA Style

Kubota, S., Promkhun, K., Sinpru, P., Suwanvichanee, C., Molee, W., & Molee, A. (2021). RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine. Animals, 11(9), 2596. https://doi.org/10.3390/ani11092596

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop