- Research article
- Open Access
Messenger RNA sequencing and pathway analysis provide novel insights into the biological basis of chickens’ feed efficiency
BMC Genomics volume 16, Article number: 195 (2015)
Advanced selection technologies have been developed and continually optimized to improve traits of agricultural importance; however, these methods have been primarily applied without knowledge of underlying biological changes that may be induced by selection. This study aims to characterize the biological basis of differences between chickens with low and high feed efficiency (FE) with a long-term goal of improving the ability to select for FE.
High-throughput RNA sequencing was performed on 23 breast muscle samples from commercial broiler chickens with extremely high (n = 10) and low (n = 13) FE. An average of 34 million paired-end reads (75 bp) were produced for each sample, 80% of which were properly mapped to the chicken reference genome (Ensembl Galgal4). Differential expression analysis identified 1,059 genes (FDR < 0.05) that significantly divergently expressed in breast muscle between the high- and low-FE chickens. Gene function analysis revealed that genes involved in muscle remodeling, inflammatory response and free radical scavenging were mostly up-regulated in the high-FE birds. Additionally, growth hormone and IGFs/PI3K/Akt signaling pathways were enriched in differentially expressed genes, which might contribute to the high breast muscle yield in high-FE birds and partly explain the FE advantage of high-FE chickens.
This study provides novel insights into transcriptional differences in breast muscle between high- and low-FE broiler chickens. Our results show that feed efficiency is associated with breast muscle growth in these birds; furthermore, some physiological changes, e.g., inflammatory response and oxidative stress, may occur in the breast muscle of the high-FE chickens, which may be of concern for continued selection for both of these traits together in modern broiler chickens.
Genetic selection has tremendously improved livestock and plant production over the past 50 years [1,2]. Advanced selection technologies have been developed and continually optimized to genetically improve traits of agricultural importance [1,3,4]. However, these methods have been primarily applied without knowledge of underlying biological changes that may be induced by selection [5,6]. Previous studies reported possible association between selection for improved performance and increased rate of physiological and metabolic disorders in modern breeds [7-9]. For example, chickens and turkeys selected for high growth rate have shown increased incidence of muscle disorders, heart failure syndrome and ascites [10-12]. A detailed characterization of traits of breeding interest may help to anticipate unfavorable consequences of long-term selection programs and adjust or perhaps redefine breeding objectives accordingly.
One of the most important traits in broiler chicken production is feed efficiency (FE), which defines the chicken’s ability to convert feed into body weight. As feed cost represents nearly 70% of the total cost in poultry production, improving FE has been an important goal in broiler chicken breeding programs . Selection for FE in broiler chickens can be accomplished using different measurements and procedures. A widely used measure of FE in broiler chickens is residual feed consumption (RFC), which is defined as the difference between an animal’s actual feed intake and expected feed intake on the basis of body weight and growth . Although moderate heritability, ranging from 0.42 to 0.45, for RFC has been reported in a previous study using more diverse chicken populations , to our knowledge this trait exhibits lower heritability (~0.2) in the commercial pure lines, possibly explaining the relatively slow progress in improving FE in commercial breeding programs. Insights into the biological basis of differences in chicken FE are required to develop more efficient and sustainable selection strategies.
Previous studies have revealed a link between mitochondrial function and FE in broiler chickens. Lower electron transport chain coupling and greater hydrogen peroxide (H2O2) production were observed in mitochondria of low-FE birds . A microarray gene expression analysis of breast muscle samples from high- and low-FE broiler chickens identified 782 differentially expressed genes [16,17]. Most of the genes up-regulated in high-FE chickens were related to anabolic metabolism, whereas genes up-regulated in low-FE chickens were associated with muscle fiber development, muscle function, cytoskeletal organization and stress response . With the rapid development of next-generation sequencing technologies, RNA sequencing (RNA-seq) has been replacing microarray technology for transcriptome-wide gene expression analysis. Avoiding technical issues inherent to microarray such as cross-hybridization and narrow ranges of signal detection, RNA-seq can provide more accurate and comprehensive information regarding changes in gene expression between different conditions or different phenotypes [18-21]. Therefore, a global gene expression study using RNA-seq is required for better understanding the molecular basis of FE in broiler chickens.
The objective of this study is to characterize the biological basis of differences between high- and low-FE chickens through breast muscle RNA-seq analysis. Using tissue samples from extreme high- and low-FE broiler chickens, the present study identifies genes and pathways differentially regulated in breast muscle between these two groups of chickens, providing important information toward understanding the biological basis of variation in FE in broiler chickens.
Animals and sample collection
Six groups of 400 male commercial broiler chickens from a cross between three commercial broiler pure lines were sampled at 29 days of age from the field in the Delmarva region of the United States and transferred into individual cages for feed efficiency measurement. The cages were arranged in rows at two levels, top or bottom levels, and each row had 100 cages. The birds were individually weighed at the beginning (29 days of age, BW29) and end of the FE test (46 days of age, BW46) and fed ad libitum until 47 days of age. At 47 days of age, chickens were euthanized by cervical dislocation. Breast muscle samples (~1-2 g) were obtained from high- and low-FE birds, immediately flash frozen in liquid nitrogen and kept at -80°C until further processing. Body weight post-euthanization (BW47) and breast muscle weight (BMW47) were also recorded and used for estimating the percentage of breast muscle [(BMW47/BW47)*100]. The total feed consumption of each bird was measured by subtracting the total amount of feed left at the end of the test (46 days of age) from the total amount of feed provided to each bird at the beginning of the test (29 days of age). To measure the broiler’s FE, residual feed consumption (RFC) was calculated using the following equation:
where FC represents the feed consumption of each bird; Level represents the fixed effects of row location (top or bottom level) on FC; Row (Level) represents the fixed effects of row nested within row location; BW29 is the initial (29-day) body weight; BW46 is the ending (46-day) body weight; c is the intercept; and b1 and b2 are the partial regression coefficients of FC on BW29 and BW46.
After excluding outliers and erroneous data (3.3%) and data from birds with defects (1.2%; leg and wings problem, etc.), samples from clinically healthy chickens exhibiting the highest (n = 12) and lowest (n = 13) RFC from the six groups of 400 birds were selected for cDNA library preparation. Two samples from the high-FE group did not produce enough cDNA libraries, so samples from 23 birds, 10 high- and 13 low-FE, were used for further analysis. The protocols were submitted to, and the use of the collected data and samples for research was approved by, the University of Delaware Agricultural Animal Care and Use Committee.
The frozen breast muscle samples were smashed into pieces by hammering. Pulverized tissues were stored at -80°C until RNA extraction. The total RNA was isolated from 70-100 mg of fragmented breast muscle tissues using a mirVana™ miRNA isolation kit (Ambion®; Austin, TX), according to the manufacturer’s protocol. RNA quantity and quality were assessed using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies; Wilmington, DE) and Agilent 2100 bioanalyzer (Agilent Technologies; Santa Clara, CA). The RNA integrity numbers of all the RNA samples were above 8.0.
RNA-seq library preparation and sequencing
In total, 23 cDNA libraries were constructed using an Illumina Truseq stranded RNA sample preparation kit following the manufacturer’s instruction (Illumina Inc.; San Diego, CA). Briefly, polyA containing mRNA molecules were purified by oligo (dT) magnetic beads and subsequently fragmented. The purified RNA fragments were reversely transcribed into first-strand cDNA using SuperScript II reverse transcriptase (Invitrogen™; Austin, TX). The second-strand cDNA was synthesized using dUTP instead of dTTP, as a result, the second-strand cDNA was not amplified during PCR because the polymerase can’t add nucleotide to dUTP. The double-strand cDNA was adenylated at the 3’ end and ligated to the Illumina indexing adapters. After PCR enrichment, cDNA quantity and quality were assessed using a NanoDrop ND-1000 spectrophotometer and Agilent 2100 bioanalyzer. The averaged size of synthesized cDNA fragments was approximately 260 bp. cDNA libraries were normalized to 10 nm for each sample and then pooled together and sequenced on four lanes of an Illumina Hiseq 2000 sequencer at Delaware biotechnology institute, University of Delaware. Approximate 68 million fragments per sample were sequenced by 75-bp from both ends.
Mapping reads to the chicken reference genome
Before read alignment, the quality of raw sequence reads was checked using the FastQC program, and nucleotide calls with a quality score of 28 or higher were considered very good quality . Sequencing reads from each sample were mapped to the chicken reference genome [Ensembl Galgal4 (GCA_000002315.2)] using the TopHat program . Because only the strand generated during the first-strand synthesis was sequenced, “-library-type fr-firststrand” was applied as one of the parameters in our read alignment using TopHat. Only one alignment for a given read was allowed in our analysis (i.e., -g 1), and both reads from a single sequence fragment were required to be mapped to the reference genome in a concordant manner (i.e., --no-discordant and --no-mixed). To summarize the alignments statistics, the resulting alignment files (SAM files) statistics were analyzed using SAMtools .
Differential expression analysis
Cuffdiff, a companion tool of Cufflinks (v 2.1.1), was used to quantify the gene expression levels and to perform a differential expression test . The fragment counts were normalized via a geometric method, as described previously . Genes with a false discovery rate of less than 5% (i.e., FDR < 0.05) were considered significant.
Nanostring nCounter® gene expression assay
The gene expression data was verified by NanoString nCounter® technology, as described previously . Briefly, 23 RNA samples were submitted to NanoString, Inc. (Seattle, WA USA) for gene expression assay. With 12 housekeeping genes, 192 endogenous transcripts were selected across multiple on-going RNA-seq projects in our laboratory as target sequences to be measured. Designs of specific probes for target sequences were provided by NanoString  and were screened to avoid areas of high SNP density. A total of 100 ng of each RNA sample were hybridized to the CodeSet®, which was composed of both capture and reporter probes . After 16 hours incubation, the samples were transferred to the nCounter® Prep Station and Digital Analyzer for transcript quantification. Positive control normalization factors and reference genes were used to normalize the raw data for biological analysis . Log2 ratios of gene expression levels between high- and low-FE groups were calculated to compare with the corresponding log2 ratio values from RNA-seq analysis.
Ingenuity pathway analysis
Genes differentially expressed (FDR < 0.05) between high- and low-FE birds were included in pathway and function analysis using Ingenuity pathway analysis (IPA; Ingenuity® Systems, http://www.ingenuity.com). The functional and canonical pathway analysis was used to identify the significant biological functions and pathways. Functions and pathways with P-value < 0.05 (Fischer’s exact test) were considered to be statistically significant. IPA’s upstream regulator analysis function was used to identify potential transcriptional regulators that could explain the observed changes in gene expression between high- and low-FE chickens. The activation z-score was calculated to predict activation or inhibition of transcriptional regulators based on published findings accessible through the Ingenuity knowledge base. Regulators with z-score greater than 2 or less than -2 were considered to be significantly activated or inhibited.
Results and discussion
A summary of the phenotype measurements from 23 high-FE (n = 10) and low-FE (n = 13) chickens is presented in Table 1. Although the initial bird weights (BW29) are not significantly different between these two groups (P = 0.661), the mean body weight of high-FE birds is significantly heavier than that of low-FE birds at the end of the test (P < 0.05), and the high-FE chickens consumed significantly less feed than low-FE birds during the test (P < 0.01). Consequently, the difference in mean RFC values between high- and low-FE chickens is highly significant (P < 0.001). The mean breast muscle weight and breast muscle percentage of the high-FE birds are significantly higher than those of low-FE birds (P < 0.05).
Transcriptional profile of chicken breast muscle
A total of 23 cDNA libraries were constructed using RNA samples of breast muscle tissues from 10 high- and 13 low-FE chickens and sequenced for 75 cycles from both ends on four lanes. In total, about 1.573 billion of 75-base sequence reads are obtained with an average of 393 million raw reads per lane. No significant difference in the number of reads between these four lanes is observed. The total number of reads for one sample ranges from 50 million to 88 million, with an average of 68 million reads per sample. Based on quality check reports, the averaged quality score of sequence reads is high, approximately 38, with the average GC content ranging from 49% to 51%. On average, 80% of the paired-end reads are properly mapped to the chicken reference genome (Ensembl Galgal4). The summary of alignment for all samples is shown in Additional file 1. The relative expression of a gene is normalized as fragments per kilobase of transcript per million mapped fragments (FPKM), which is proportional to the number of cDNA fragments originated from the gene transcript. The lowest limit of gene expression value is set to be 0.1 FPKM in at least one of the 23 samples. According to this limit, 14,148 genes are identified as being expressed in the breast muscle tissues. To assess the consistency of the gene expression levels between different samples, the Pearson’s correlation coefficient was calculated for each pairwise combination of samples . The averaged pairwise correlation coefficient between samples is 0.794, reflecting pretty consistent gene expression profiles.
Gene differential expression analysis
Differentially expressed genes were detected by Cuffdiff, an internal program of Cufflinks. Of 17,107 genes in the Ensemble database (Ensembl Galgal4), 1,059 were identified as significant genes with different expression levels between high- and low-FE chickens (q-value < 0.05) (Additional files 2 and 3). All of this group of 1,059 genes have a fold change greater than 1.3, and 642 genes (60.6%) have a fold change above 1.5. Among the 1,059 differentially expressed genes, 327 and 732 genes are down- and up-regulated in high-FE birds, respectively (Additional file 4). This relative imbalance in the number of down- and up-regulated genes is likely due to the increased breast muscle regeneration and inflammatory response in the high-FE chickens (discussed below). Since muscle development and inflammatory response require higher levels of activators such as growth factors, hormones and cytokines, the gene expression may be positively regulated by these activators in the breast muscle of the high-FE birds.
Confirmation of RNA-seq data
To verify the gene expression data obtained from RNA-seq analysis, we selected 192 target genes (71 significant and 121 non-significant) along with 12 housekeeping genes for the NanoString nCounter® assay. Comparison of the normalized counts from NanoString with FPKM values derived from RNA-seq shows high concordance, with pair-wise Pearson’s correlation coefficients ranging from 0.70 to 0.98. The Pearson’s correlation coefficients of fold change in gene expression levels [log2 ratio(high-FE/low-FE)] between NanoString and RNA-seq results are also high: 0.7532 for all genes and 0.8332 for the 71 significantly differentially expressed genes (Figure 1). The correlation of log2 fold-change between two analyses is notably affected by lowly expressed genes, and increases by excluding genes with low expression levels (Figure 2). This can explain why the significantly differentially expressed genes show greater correlation compared with all the selected genes, because the FPKM values of the significant target genes are equal or greater than 0.4, whereas 23 genes out of the 121 non-significant target genes have an FPKM value less than 0.4.
Overview of IPA analysis
To fully interpret the biological implications of the results from the differential expression analysis, all significant genes with their respective log2 fold-change were submitted for Ingenuity® Pathway Analysis. The top 10 up-regulated and top 10 down-regulated genes in high-FE chickens are listed in Additional file 5.
A summary of the IPA analysis, including top five biological functions and canonical pathways, are presented in Table 2. Generally, most of the differentially expressed genes are related to immune response and metabolic processes. Genes up-regulated in high-FE birds are associated with cellular function, movement, growth and proliferation, cell-to-cell signaling and interaction and cell death and survival. In contrast, genes down-regulated in high-FE birds are associated with metabolic processes including lipid, nucleic acid and carbohydrate metabolism, molecular transport and small molecule biochemistry (Table 3). Differing from the results of earlier work on chicken FE conducted using 44 K oligo microarray [16,17], genes involved in muscle fiber development, cytoskeletal organization and stress response are found to be up-regulated in the high-FE chickens (rather than in the low-FE chickens) in the current study. The discrepancy is probably from the different genetic composition of broiler chickens between two studies. Birds in previous studies are from a male broiler pure line that has been observed greater reactive oxygen species (ROS) production in the low-FE chickens than the high-FE birds , which is reverse of our findings that indicate ROS production is increased in the high-FE birds. Since ROS can act as a second messenger and mediate gene expression in the cell through signal transduction, differential gene expression is likely driven, in part, by the inherent differences that are modulated by ROS-mediated mechanisms. Further comparison of broiler chickens used in our study with pure line chickens used in these previous FE studies will be explained below.
Upstream regulator analysis through IPA predicted the cascade of upstream transcriptional regulators that can potentially explain the differences in gene expression profile between high- and low-FE chickens. A summary of the upstream regulators identified by IPA is presented in Additional file 6. A total of 27 transcriptional regulators are predicted to be activated or inhibited (24 activated and 3 inhibited) in the high-FE chickens, of which 24 regulators are considered to be significant with P-value < 0.05 (21 activated and 3 inhibited).
Increased muscle growth and remodeling in high-FE chickens
Of all differentially expressed genes, 32 are associated with muscle development (Additional file 7), supporting the increased breast muscle weight in the high-FE birds. Among them, both hepatocyte growth factor (HGF) and insulin-like growth factor 2 (IGF2) encode key growth factors that have autocrine or paracrine effects on chicken skeletal muscle development and regeneration . HGF can not only activate the proliferation of quiescent muscle satellite cells, it also can induce the migration of activated satellite cells to the injured sites . IGF2 acts as a crucial regulator in muscle regeneration by stimulating muscle cell differentiation as well as inducing muscle cell hypertrophy . Other muscle growth-related genes that are up-regulated in the high-FE chickens include myogenin (MYOG), cysteine and glycine-rich protein 3 (CSRP3), myoferlin (MYOF), glypican 1 (GPC1), protein tyrosine phosphatase, receptor type, A (PTPRA) and gap junction protein (GJA1). As a member of myogenic regulatory factors (MRFs), MYOG is essential for the fusion of myoblasts into myotubes during muscle growth and regeneration . The CSRP3 gene encodes muscle LIM protein, which is able to increase the activity of MRFs and plays a critical role in enhancing myogenesis . The MYOF-encoded protein is a fundamental modulator for myoblast fusion, highly expressed during muscle repair and regeneration . The stimulatory effects of GPC1 on muscle satellite cell differentiation and myotube formation was reported in turkeys . The protein encoded by the PTPRA gene is a signaling molecule that was found to increase myogenesis of rat muscle L6 cells . The GJA1-encoded protein is a major component in gap junctions and required for myogenesis and regeneration . Collectively, the up-regulation of genes that can positively regulate muscle growth indicates that muscle growth and development is elevated in the high-FE chickens.
In addition to genes involved in muscle development, genes associated with muscle hypertrophy, including F-box protein 32 (FBXO32; fold change = −1.879), F-box protein 40 (FBXO40; fold change = −1.879), F-box protein 9 (FBXO9; fold change = −1.347), forkhead box O3 (FOXO3; fold change = −1.540) and myotrophin (MTPN; fold change = 1.426), are found differentially expressed in the breast muscle of chickens with high versus low FE. MTPN, a well-known positive growth factor in promoting muscle growth , is up-regulated in high-FE chickens. The increased MTPN expression may indicate that myocyte growth and protein synthesis are augmented in the breast muscle of high-FE birds, accordingly, contributing to breast muscle hypertrophy in these chickens. Furthermore, the down-regulation of FOXO3 and F-box family proteins in high-FE chickens further explains muscle growth differences between high- and low-FE chickens. Protein encoded by FOXO3 is a master regulator of both autophagy-lysosome and ubiquitin-proteasomal pathways, promoting protein degradation and thereby contributing to muscle atrophy . Proteins from the F-box family mediate the interaction between substrates and ubiquitin-conjugating enzymes, which facilitate proteolysis in diverse tissues . Of them, the FBXO32-encoded protein, known as atrogin 1, is a well-recognized muscle-specific ubiquitin ligase leading to muscle atrophy in a wide range of diseases [41-43]. The FBXO40-encoded protein also has been proposed to play a role in muscle atrophy in mammals . Thus, the decreased expression of atrophy-related genes in breast muscle of high-FE birds suggests that muscle protein loss is reduced in high-FE chickens in contrast to low-FE birds. The transcription of these genes is regulated by the PI3K/Akt signaling pathway, which will be discussed later. Taken together, the up-regulation of MTPN and down-regulation of FOXO3 and FBXO family genes in the high-FE chickens suggest that birds with high FE may have elevated protein synthesis and decreased protein degradation in their breast muscle.
Genes associated with extracellular matrix (ECM) remodeling are also up-regulated in the high-FE birds. The ECM of skeletal muscle serves as a scaffold for maintaining the structure of muscle and guiding new fiber formation . Muscle regeneration is frequently accompanied by the degradation of ECM because it facilitates satellite cell migration to specific sites for proliferation and fusion into myotubes [32,46]. Therefore, the up-regulation of genes involved in ECM remodeling implies that muscle remodeling is increased in the breast muscle of high-FE chickens. Matrix metalloproteinases (MMPs) are the main endopeptidases responsible for degrading all kinds of ECM, consequently, playing an important role in mediating muscle cell migration and regeneration [47,48]. As presented in Additional file 8, six genes from the MMP family are differentially expressed in our study, all of which are up-regulated in high-FE birds. Of the proteins encoded by these genes, MMP1 and MMP13 belong to MMP collagenases that are capable of cleaving interstitial collagen types I, II and III . Through an in vitro wound-healing assay, MMP1 was able to promote myoblast migration and differentiation by increasing the expression of N-cadherin and β-catenin or pre-MMP-2 and TIMP . MMP13 also plays a role in muscle regeneration, expressing in all muscle cells during muscle regeneration, and its expression level is positively correlated with the extent of muscle damage . MMP9 is a gelatinase that also relates to muscle regeneration . Evidence showed that the expression levels of MMP9 were greatly increased in response to inflammation and the activation of satellite cells in injured muscle [52-54]. However, contrary to the positive function of MMP1 on muscle regeneration, a recent study revealed that MMP9 could lead to muscle cell necrosis, inflammation and fibrosis . Collectively, the up-regulation of MMPs in the high-FE birds suggests an augmented muscle remodeling in these birds compared with the low-FE chickens.
The IPA upstream regulator analysis provides additional support to our conclusions regarding muscle development and remodeling in the high-FE birds. According to the predictions from IPA, several transcriptional factors involved in muscle development are activated in high-FE chickens. As a regulator of postnatal muscle growth, JunB is predicted as being activated in the breast muscle of the high-FE birds (P-value = 1.2E-03, z-score = 2.200). JunB can stimulate myosin expression to elevate protein synthesis, accordingly, maintaining muscle mass and promoting muscle hypertrophy . Additionally, JunB can suppress the transcription of FOXO3 and thereby inhibit protein degradation in muscle ; therefore, activation of JunB in the high-FE birds may be one of the causes of the down-regulation of FOXO3 as well as of FBXO32 in breast muscle of these birds. In addition, through the IPA analysis, JunB is predicted to activate the transcription of MMP1, MMP9, MMP13, fibronectin 1 (FN1), heme oxygenase 1 (HMOX1), neutrophil cytosolic factor 2 (NCF2) and interleukin 1 receptor-like 1 (IL1RL1) (Figure 3A). As discussed above, FN1, MMP1 and MMP13 are all positively correlated with muscular satellite cell proliferation. Thus, JunB may also have increased myogenesis through activating transcription of these genes in the high-FE birds.
Apart from JunB, a main transcriptional factor in the formation of mature sarcomeres, myocyte-specific enhancer factor 2C (MEF2C), is predicted as being activated in breast muscle of high-FE birds [58,59]. Protein encoded by MEF2C is a member of the myocyte enhancer factor 2 (MEF2) family, which directly cooperates with MRFs and enhances skeletal muscle development . In the present study, MEF2C is predicted to be an activated upstream regulator that increases the transcription of GJA1, MMP13, MYOG, myozenin 2 (MYOZ2; fold change = 2.400) and ATPase, Ca++ transporting, cardiac muscle, slow twitch 2 (ATP2A2) (fold change = 1.390) (Figure 3B). GJA1, MMP13 and MYOG are all involved in myogenesis and exert positive effects on skeletal muscle growth and regeneration, which has been discussed above. Thus, MEF2C’s activation may augment the muscle development in high-FE birds. Moreover, MYOZ2 is also predicted as being up-regulated by MEF2C. The MYOZ2-encoded protein belongs to a family of calcineurin-interacting proteins that modulates specific skeletal muscle signaling pathways through suppressing calcineurin . It has been shown that MYOZ2 plays a role in regulating myocyte differentiation and promoting slow-oxidative fibers growth . Collectively, MYOZ2 may be more active in breast muscle of high-FE chickens and, consequently, mediates some biological pathways and leads to muscle remodeling in these birds.
Growth hormone (GH) and IGFs/PI3K/Akt signaling pathway over-represented in the differentially expressed genes
Through the IPA canonical pathway analysis, several critical pathways in the regulation of body and muscle growth are over-represented among the differentially expressed genes. One of these pathways is the GH signaling pathway, enriched by 10 genes in our dataset (P-value = 3.25E-04; ratio = 1.32E-01) (Figure 4). As a key mediator of body size, GH has an anabolic effect on skeletal muscle development . Through binding to growth hormone receptor (GHR) in muscles, biologically active GH can activate nuclear receptor STAT5 and thereby induce the synthesis and secretion of IGF-1 as well as IGF-2 (fold change = 1.657) . Furthermore, GH can stimulate signaling molecules including PI3K [phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit beta (PIK3CB; fold change = 1.663); phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit delta (PIK3CD; fold change = 1.709); phosphoinositide-3-kinase, regulatory subunit 5 (PIK3R5; fold change = 1.564)] and protein kinase C (PKC) [Protein kinase C delta type (PRKCD; fold change = 1.558)], leading to the activation of Akt/PKB signaling pathway and STAT5 . Both the PI3K/Akt/PKB pathway and IGFs are crucial contributors to muscle hypertrophy, which will be discussed later. Most of the mapped genes (8 of 10) are up-regulated in the high-FE chickens, indicating that the GH signaling pathway is more activated in the breast muscle of the high-FE birds compared with the low-FE birds. The two down-regulated genes [GHR and insulin-like growth factor binding protein (IGFBP3)] also fit this assumption. Evidence from the literature indicates that the expression of GHR is inversely correlated with the concentration of GH [66,67]. Thus, the down-regulation of GHR gene expression in the high-FE birds may be due to relatively high circulating GH levels in these birds. In spite of the stimulatory effects of GH on IGFBP3 transcription, there may be other modulators inhibiting the expression of IGFBP3 in high-FE chickens, consequently exerting an inhibitory effect on IGF-1 function .
Another important finding is that the IGFs/PI3K/Akt signaling pathway is over-represented by the differentially expressed genes. The IGFs/PI3K/Akt signaling pathway plays a key role in the regulation of muscle growth and muscle hypertrophy in a variety of organisms [69-71]. Nine differentially expressed genes are mapped to the IGFs/PI3K/Akt pathway (Figure 5). Of these, PIK3CD (fold change = 1.709), PIK3CB (fold change = 1.663) and PIK3R5 (fold change = 1.564) are up-regulated in the high-FE chickens, implying that the PI3K complex is more active in the breast muscle of these birds. The up-regulated members of the PI3K complex are predicted to increase PI3K-Akt cascade activity in the high-FE birds by IPA (Figure 5). Activated PI3K can induce the phosphorylation of phosphatidylinositol 4,5-bisphosphate (PIP2) to generate phosphatidylinositol-3,4,5-trisphosphate (PIP3). PIP3 acts as a docking site for phosphoinositide-dependent kinase 1 (PDK1) and Akt and subsequently contributes to the activation of Akt [72,73]. On the basis of the IPA prediction, the activated Akt translocates into the nucleus and then inhibits the transcription of the forkhead box O (FOXO) family, which is consistent with the down-regulation of the forkhead box O3 (FOXO3) gene (fold change = −1.540) in our results . As mentioned above, FOXO3 can promote protein degradation and muscle atrophy . Hence, considering the expression profile of the mapped genes, the protein degradation process is predicted to be reduced in the breast muscle of the high-FE chickens as a result of PI3K/Akt pathway activation. In addition, because activation of Akt can up-regulate the transcription of ATP citrate lyase (ACLY) through suppressing the activity of glycogen synthase kinase (Gsk3), the increased expression of ACLY in the high-FE birds (fold change = 1.345) lends more support to the assumption that the PI3K/Akt pathway is activated in the high-FE chickens. Apart from repressing protein degradation, the activated PI3K/Akt pathway can promote protein synthesis in muscle via inhibiting Gsk3 , which is also inferred in our analysis. Therefore, the up-regulated IGFs/PI3K/Akt pathway suggests increased protein synthesis as well as decreased protein degradation in the breast muscle of the high-FE birds, explaining in part why high-FE chickens have more breast muscle than do low-FE birds. Previous studies on chicken FE also found that gene encoding PI3K was up-regulated in the high-FE chickens and a list of differentially expressed genes associated with the Akt/mTOR pathways, which strengthened our conclusion [17,74].
Inflammatory response in the breast muscle of high-FE chickens
In the present study, a large number of the differentially expressed genes (136 genes) are involved in inflammatory response. Most of these genes (124 genes), including genes encoding for interleukin 8 (IL-8) and chemokine (C-X-C motif) ligand 14 (CXCL14), are expressed greater in the high-FE chickens. Although the cellular source of IL-8 and CXCL14 remains unknown in the current study, both not only exert direct effects on immune cell recruitment but also act as paracrine or endocrine factors in skeletal muscle. IL-8 has been recently classified as a myokine that can promote angiogenesis within the muscle [75,76]. CXCL14 is encoded by an obesity-induced gene in mice that inhibits the insulin-induced glucose uptake in cultured myocytes . In addition, the gene encoding for corticotropin-releasing hormone (CRH) is also up-regulated in the high-FE birds (fold change = 2.824). Previous studies have demonstrated that CRH is secreted from nerve terminals and epithelial cells at inflammation sites and has a local proinflammatory effect on resident immune cells . Therefore, it is likely that the elevated transcription of CRH functions to augment an immune response in the breast muscle of high-FE chickens. Apart from its immunomodulatory role, an increase in CRH may have a positive impact on thermogenesis of skeletal muscle in high-FE birds .
A series of genes encoding for cytokine receptors are also up-regulated in the high-FE chickens, further indicating an augmented immune response in the breast muscle of the high-FE chickens. These genes include chemokine (C-C motif) receptor 2 (CCR2), chemokine (C-C motif) receptor 5 (CCR5), interleukin 17 receptor A (IL17RA), interleukin 18 receptor 1(IL18R1), interleukin 1 receptor, type I (IL1R1), interleukin 1 receptor, type II (IL1R2), interleukin 1 receptor-like 1 (IL1RL1), interleukin 2 receptor, gamma (IL2RG) and interleukin 5 receptor, alpha (IL5RA). Among them, CCR2 was found to be expressed in infiltrating macrophages and playing a crucial role in muscle regeneration . This gene can recruit macrophages to injured muscle, which then produces a high level of IGF-I to promote muscle regeneration . Therefore, the up-regulation of CCR2 suggests that, compared with the low-FE chickens, macrophage infiltration and muscle regeneration are increased in the breast muscle of the high-FE birds.
The IPA canonical pathway analysis also supports our hypothesis regarding augmented immune response and active recruitment of immune cells to the breast muscle of the high-FE chickens. Several over-represented pathways involved in inflammatory response are identified in our analysis (Additional file 9). Given that nearly all genes mapped to these pathways are up-regulated in the high-FE chickens, we conclude that these immune-related pathways are activated in the breast muscle of the high-FE chickens. According to the predictions from IPA, a number of transcription factors associated with inflammatory response are also activated in the breast muscle of the high-FE chickens: v-ets erythroblastosis virus E26 oncogene homolog 1 (ETS1), spleen focus forming virus (SFFV), spi-1 proto-oncogene (SPI1), X-box binding protein 1 (XBP1) and runt-related transcription factor 1(RUNX1).
Muscle inflammation is a key step in muscle remodeling, which can clean disrupted muscle cells and promote muscle regeneration by releasing growth factors . A variety of circumstances (e.g., muscle injury, exercise and obesity) can activate transcription factors NF-kB and c-Jun/FOS in muscle cells, resulting in the expression and secretion of several factors . These factors, including cytokines and other non-protein mediators, can either directly attract circulating immune cells to the muscle or activate resident immune cells, providing chemotactic signals for recruitment . As a consequence, a number of immune cells are recruited to the muscle, phagocytizing the cellular debris and producing cytokines affecting muscle cells [83,85]. The IPA upstream regulator analysis predicts that transcription factors JUN and FOS are activated in the breast muscle of the high-FE birds (Figure 6). Protein encoded by JUN and FOS are components of activator protein 1 (AP-1), which is an important transcription factor responding to various physiological and pathological stimuli . Overall, our results suggest that, compared with the low-FE birds, the high-FE birds experienced a more intense muscle restructuring and higher inflammatory responses in the breast muscle.
Free radical scavenging enriched in the differentially expressed genes between high- and low-FE chickens
Several differentially expressed genes in our dataset are involved in the production of ROS. Genes encoding for ROS-generating enzymes, including cytochrome b-245, beta polypeptide (CYBB) [fold-change = 2.08] and NADPH oxidase organizer 1 (NOXO1) [fold-change = 2.38], are all up-regulated in the high-FE birds, suggesting that ROS production is increased in the breast muscle of these birds compared with the low-FE birds. CYBB, also known as NADPH oxidase 2 (NOX2), is a major enzyme responsible for superoxide production in the sarcoplasmic reticulum . NOXO1, a positive mediator of NOX1 and NOX3, initiates the activity of NOX1 and NOX3 for generating ROS . Moreover, the down-regulation of uncoupling protein 3 (UCP3) [fold-change = -1.67] in the high-FE birds may indicate that mitochondria from the breast muscle of the high-FE chickens have higher electron transport chain coupling compared with that from low-FE chickens. This assumption is consistent with previous findings [15,89]. Because UCP3-mediated uncoupling can attenuate the production of ROS , the down-regulation of UCP3 in the high-FE birds may also suggest a higher production of ROS from the mitochondria of the breast muscle of these birds. Collectively, our data suggest that, compared with the low-FE birds, ROS is produced at a higher level in the breast muscle of the high-FE chickens.
The IPA downstream effect analysis supports our hypothesis regarding increased ROS production in the breast muscle of high FE-chickens. Processes, including metabolism of reactive oxygen species (P-value = 5.77E-07), synthesis of reactive oxygen species (P-value = 1.75E-06), production of reactive oxygen species (P-value = 4.92E-06) and production of superoxide (P-value = 2.05E-03), are predicted to be increased in the high-FE chickens. In addition, the NRF2-mediated oxidative stress response pathway is over-represented among the differentially expressed genes, with 17 genes (P-value = 6.96E-04; ratio = 0.089) mapped to this pathway (Figure 7A). Nuclear factor (erythroid-derived 2)-like 2 (NRF2), also known as NFE2L2, is a key transcription factor in cells that responds to a range of oxidative and xenobiotic stresses . Upon exposure of cells to various stimuli such as ROS and electrophilic compounds, quiescent NRF2 in cytoplasm is activated through phosphatidylinositol 3-kinase (PI3K), RAS and PKC signaling pathways . By phosphorylation or binding to actin, activated NRF2 translocates into the nucleus and binds to the antioxidant response elements, initiating the transcription of a number of genes encoding for antioxidants and ROS detoxifying enzymes . A group of NRF2 downstream genes, including v-maf musculoaponeurotic fibrosarcoma oncogene homolog F (avian) (MAFF) [fold-change = 1.61], glutathione S-transferase A3 (GSTA3) [fold-change = 2.03], glutathione S-transferase omega 1 (GSTO1) [fold-change = 1.50], heme oxygenase (decycling) 1 (HMOX1) [fold-change = 1.36], microsomal glutathione S-transferase 1 (MGST1) [fold-change = 1.39], superoxide dismutase 3 (SOD3) [fold-change = 1.85], thioredoxin (TXN) [fold-change = 1.50], peptidylprolyl isomerase B (PPIB) [fold-change = 1.43], aldehyde oxidase 1 (AOX1) [fold-change = -1.63], DnaJ (Hsp40) homolog, subfamily A, member 1 (DNAJA1) [fold-change = -1.39] and DnaJ (Hsp40) homolog, subfamily C, member 1 (DNAJC1) [fold-change = -1.35], are differentially expressed in the current study. Genes encoding for antioxidant proteins, such as SOD3, HMOX1 and TXN, are up-regulated in the high-FE birds. Protein encoded by SOD3 is an extracellular protective enzyme against not only ROS but also inflammation, thus playing a role in tissue recovery . HMOX1 is increased in the condition of oxidative stress and has an effect on protecting cells against ROS and inflammation . The TXN-encoded protein is involved in a range of redox reactions and can decrease the quantity of ROS . The up-regulation of TXN, SOD3 and HMOX1 indicates that an NRF2-mediated antioxidant response is activated in the breast muscle of the high-FE chickens. Additionally, three members from the glutathione s-transferase (GST) group, encoded by GSTO1, GSTA1 and MGST1, are all up-regulated in the high-FE birds. GST is known for its function in detoxification of xenobiotics as well as endogenous metabolites . The increased expression of the GST superfamily also suggests that responses to oxidative stress are elevated in the breast muscle of the high-FE chickens. Although few genes in the NRF2 signaling pathway, including AOX1, DNAJA1 and DNAJC1, are down-regulated in the high-FE chickens, overall there are 14 up-regulated genes mapped to this pathway, indicating that NRF2-mediated oxidative stress response is augmented in the breast muscle of the high-FE birds. Moreover, NRF2 (NFE2L2), a transcription factor, is predicted to be activated in the high-FE chickens (P-value = 1.94E-05; z-score = 2.036) (Figure 7B). Taken together, our results suggest a higher level of ROS generated in the breast muscle of high-FE chickens.
However, in contrast to our findings, Bottje et al. (2002) reported higher amounts of ROS produced in the breast muscle of their low-FE birds . This inconsistency is likely caused by the difference in broiler chickens between two studies. Male breeders, presumably with relatively low breast muscle yield, were studied in the Bottje et al. (2002) research , whereas we study broiler chickens from a commercial cross with high breast muscle yield. The ancestors of this cross have been intensively selected for the disproportionate growth of breast muscle, and the resulting higher levels of variation in breast muscle in the broiler cross may be responsible for a significant part of the variation in FE in this cross compared to the male breeder strain in the study by Bottje et al. . In regard to broiler chickens in the current study, intensive inflammatory response is possibly a major source of increased ROS in the breast muscle of the high-FE chickens. ROS-generating enzymes, such as NOX in muscle cells, can be stimulated through extracellular inflammatory cytokines including interleukin (IL)-1, IL-6 and IL-8 in a ligand-induced pattern [99,100]. Furthermore, the implied infiltrating immune cells in the breast muscle of high-FE birds may be another cause for increased ROS. It is well known that immune cells produce a large amount of ROS to support their functions during inflammation . Hence, in our study, strong indications for elevated ROS production in the breast muscle of the high-FE chickens are likely due to augmented inflammatory response, whereas the higher level of ROS observed in the study by Bottje et al. (2002) is possibly from mitochondria of breast muscle cells. Further study of genes associated with free radical scavenging may support our assumption. Indeed, in our study, a large part of these genes (Additional file 10) are also related to inflammatory response (P-value = 1.03E-23--5.15E-06), suggesting that production of ROS in the high-FE birds is closely associated with an increased immune response in the breast muscle.
Notably, growth factors including HGF, IGF-1 and fibroblast growth factor (FGF)-2 are also found to be able to induce intracellular generation of ROS in different types of cells . As mentioned above, the breast muscle of high-FE birds have higher expression of HGF and IGF-2, which may play a role in stimulating ROS production in these birds. Moreover, such generated ROS exerts insulin-mimicking effects on the insulin/IGFs signaling pathway, which has shown to be a second messenger in insulin/IGFs signal transduction under physiological conditions . Therefore, in the breast muscle of the high-FE birds, the insulin/IGFs receptor signaling pathway may be activated, in part, because of increased ROS production.
Higher ROS production may also lead to an increase in intracellular calcium concentration. It has been found that ROS mediates the influx of extracellular Ca2+ and mobilization of intracellular Ca2+ stores [103-105]. In the present study, genes involved in calcium transport [solute carrier family 8, member B1 (SLC8B1), phospholipase C, beta 2 (PLCB2) and ATPase, Ca ++ transporting, cardiac muscle, slow twitch 2 (ATP2A2)] are all up-regulated in the high-FE birds, indicating increased calcium mobilization in the breast muscle of these birds. ATP2A2 encodes sarcoplasmic reticulum Ca2+-ATPase isoform 2 (SERCA2), which is an important pump responsible for muscle relaxation through transporting Ca2+ from the cytosol into the sarcoplasmic reticulum lumen in muscle cells . Because more SERCA2 are needed to maintain calcium homeostasis when high Ca2+ levels are present in cytosols, the up-regulation of ATP2A2 in the high-FE birds may imply a high level of cytosolic Ca2+ in the breast muscle of these chickens compared with the low-FE birds.
Transcriptional regulation of hypoxia-inducible factor-1α (HIF1α)
Hypoxia-inducible factor-1α (HIF1α) is a key transcription factor that mediates cell adaption to hypoxia through regulation of a variety of gene expression . Although HIF1α mRNA is constantly expressed in cells under both normoxic and hypoxic conditions, the HIF1α protein has a very short half-life in normoxia because of degradation through the ubiquitin-proteasome system . During hypoxia, HIF1α degradation is repressed. As a result, HIF1α translocates into the nucleus and activates downstream genes in response to low O2 tension . In our study, HIF1α mRNA, aryl-hydrocarbon receptor nuclear translocator 2 (HIF2β) mRNA as well as HIF1α inhibitor hypoxia inducible factor 1, alpha subunit inhibitor (HIF1AN) mRNA are differentially expressed in the breast muscle between high- and low-FE chickens. HIF1α and HIF2β are up-regulated in the high-FE birds (fold change = 1.341 and 1.42, respectively), whereas HIF1α inhibitor HIF1AN is down-regulated in these birds (fold change = −1.343). Although the up-regulated HIF1α and HIF2β mRNA can’t represent increased amounts of stabilized HIF1α protein in the breast muscle of the high-FE chickens, decreased expression of HIF1AN may imply that HIF1α activity is increased in the breast muscle of the high-FE chickens compared with the low-FE birds. This assumption is supported by the expression of HIF1α downstream genes. As a transcription factor, HIF1α is predicted to be activated in the high-FE birds through the IPA upstream regulator analysis (P-value = 3.85E-06; z-score = 2.332; Figure 8). Indeed, a large number of the HIF1α target genes are up-regulated in the high-FE birds, indicating the activation of HIF1α in these birds.
Moreover, the HIF1α signaling pathway is over-represented among significantly differentially expressed genes (P-value = 7.58E-04; ratio = 1.39E-01). In response to hypoxia or a variety of peptide stimulators under normoxic conditions, PI3K/Akt and MAPK signaling pathways are activated to induce the accumulation of HIF1α in human cells [108,109]. Consequently, the accumulated HIF1α is translocated to the nucleus to modulate the transcription of genes involved in angiogenesis, glucose metabolism, matrix metabolism, erythropoiesis and apoptosis . In our study, with increased expression of PIK3CB, PIK3CD, PIK3R5 and muscle RAS oncogene homolog (MRAS), both the Akt/PI3K and MAPK signaling pathways are predicted to be activated in the high-FE chickens. The activated Akt/PI3K and MAPK signaling pathways may stimulate the induction of HIF1α, as reflected by the up-regulation of HIF1α and its downstream genes [glucose transporter type 3 (GLUT3), glucose transporter-like protein 5 (GLUT5), matrix metallopeptidase 1 (MMP1), matrix metallopeptidase 7 (MMP7), matrix metallopeptidase 9 (MMP9), matrix metallopeptidase 13 (MMP13), matrix metallopeptidase 27 (MMP27), matrix metallopeptidase 28 (MMP28) and lactate dehydrogenase B (LDHB)] in the high-FE birds. Based on the gene expression profile, we conclude that, compared with the low-FE birds, the activity of HIF1α signaling pathway is increased in the breast muscle of the high-FE birds.
Although it is unclear from our results whether hypoxia and/or mediators such as IGFs induced HIF1α activation in the breast muscle of the high-FE birds, we would like to speculate here about potential mechanisms underlying this activation. It is widely accepted that inflammation and hypoxia are closely interdependent in a wide array of physiological and pathological conditions [110-114]. Inflammation is frequently accompanied with hypoxia because of the high oxygen consumption of infiltrating immune cells . Assuming an increased inflammatory response in the breast muscle of the high-FE birds, we speculate that the up-regulation of HIF1α is partly caused by an inflammation-induced hypoxia. Alternatively, the up-regulation of HIF1α may be caused by excessive muscle remodeling, which may be the result of selection for breast muscle proportion. Elevated muscle growth and rearrangements may have led to the reconstruction of vasculature, consequently reducing the blood flow and resulting in oxygen deficiency in the breast muscle of the high-FE birds . Furthermore, insulin and IGFs have shown to be modulators of HIF1α induction during both normoxia and hypoxia . Given that IGF2 is up-regulated in the breast muscle of the high-FE birds, this growth factor may also have contributed to the activation of HIF1α.
Finally, the activation of HIF1α may also be partly due to a higher production of ROS in the breast muscle of the high-FE chickens. Studies have found that ROS are essential for the stabilization of HIF1-DNA, thereby triggering HIF1α-induced transcription [116,117]. It was also proposed that cellular ROS-producing proteins could sense changes in cellular oxygen concentration . Evidence indicated that low oxygen tension inhibited mitochondrial electron transport and therefore increased ROS production. The generated ROS then acted as a second messenger that contributed to HIF1α activation . Thus, the ROS production may have been increased in the breast muscle of the high-FE chickens partly because of a relatively low oxygen concentration within this tissue, which in turn may have played a role in HIF1α activation.
The current study provides a global view of gene expression differences in the breast muscle of broiler chickens with extremely high and low FE from a population of a modern commercial high-meat-yield broiler cross. To our knowledge, this study reports for first time the RNA-seq analysis of a trait of selection and breeding importance in chickens. We identify 1,059 genes significantly differentially expressed in the breast muscle between high- and low-FE chickens based on the RNA-seq experiment. Furthermore, we achieve a large-scale validation of our RNA-seq experiment by quantifying the expression of a large number of target genes (192 transcripts + 12 house-keeping genes) using a high-sensitive non-PCR-based method, i.e. NanoString nCounter® Technology . Function and pathway analysis of the differentially expressed genes sheds light on some of the underlying mechanisms that regulate chicken FE. Birds with high FE exhibit higher expression of genes involved in muscle growth, development and remodeling, which may explain why these birds have more breast muscle than do the low-FE chickens. Pathway analysis shows that anabolic pathways, including growth hormone signaling and IGFs/PI3K/Akt signaling pathways, are more activated in the high-FE birds, which may have not only led to the increased muscle growth in the high-FE chickens but also contributed to the feed conversion advantages of these birds. Our results also suggest that transcriptional factors JunB and MEF2C play crucial roles in regulating muscle growth and remodeling in high-FE chickens.
Furthermore, most of the genes up-regulated in the high-FE birds are associated with inflammatory response and oxidative stress, suggesting augmented inflammation and oxidative stress in the breast muscle of these birds. Our results also show increased activity of HIF1α, which may be caused by a lower oxygen environment in the breast muscle of high-FE chickens. Although no clinical symptoms of sickness or muscle damage were observed in the birds used in the current study, some of the molecular changes in the high-FE chickens may be hypothesized to lead to recently reported muscle quality issues in modern broiler chickens such as white striping and wooden breast [120-122]. These disorders have been reported to be more frequent in birds with high breast muscle weight and high FE, suggesting that the susceptibility may be primarily induced by breeding for these traits. Further investigation (e.g., histological and protein analysis) would be helpful for examining inflammation and oxidative stress in the breast muscle of high-FE and high-breast-muscle-yield birds.
Availability of supporting data
The data supporting the results of this article are included within the article and its additional files. Readers may contact the corresponding author for additional information.
Hayes BJ, Lewin HA, Goddard ME. The future of livestock breeding: genomic selection for efficiency, reduced emissions intensity, and adaptation. Trends Genet. 2013;29:206–14.
Dekkers JCM, Hospital F. The use of molecular genetics in the improvement of agricultural populations. Nat Rev Genet. 2002;3:22–32.
Goddard ME, Hayes BJ. Mapping genes for complex traits in domestic animals and their use in breeding programmes. Nat Rev Genet. 2009;10:381–91.
Deeb N, Lamont SJ. Genetic architecture of growth and body composition in unique chicken populations. J Hered. 1998;93:107–18.
Emmerson DA. Commercial approaches to genetic selection for growth and feed conversion in domestic poultry. Poult Sci. 1997;76:1121–5.
Rauw W, Kanis E, Noordhuizen-Stassen E, Grommers F. Undesirable side effects of selection for high production efficiency in farm animals: a review. Livest Prod Sci. 1998;56:15–33.
Dransfield E, Sosnicki AA. Relationship between muscle growth and poultry meat quality. Poult Sci. 1999;78:743–6.
Wilson BW, Nieberg PS, Buhr RJ, Kelly BJ, Shultz FT. Turkey muscle growth and focal myopathy. Poult Sci. 1990;69:1553–62.
Petracci M, Cavani C. Muscle growth and poultry meat quality issues. Nutrients. 2012;4:1–12.
Julian RJ. Rapid growth problems: ascites and skeletal deformities in broilers. Poult Sci. 1998;77:1773–80.
Mills LJ, Mitchell MA, Mahon M. Incidence of skeletal muscle damage in selected and unselected strains of turkey. Br Poult Sci. 1998;1998(39 Suppl):S27–8.
Dunnington EA, Siegel PB. Enzyme activity and organ development in newly hatched chicks selected for high or low eight-week body weight. Poult Sci. 1995;74:761–70.
Willems OW, Miller SP, Wood BJ. Aspects of selection for feed efficiency in meat producing poultry. World’s Poult Sci J. 2013;69:77–88.
Aggrey SE, Karnuah AB, Sebastian B, Anthony NB. Genetic properties of feed efficiency parameters in meat-type chickens. Genet Sel Evol. 2010;42:25.
Bottje W, Iqbal M, Tang ZX, Cawthon D, Okimoto R, Wing T, et al. Association of mitochondrial function with feed efficiency within a single genetic line of male broilers. Poult Sci. 2002;81:546–55.
Kong B-W, Song JJ, Lee JY, Hargis BM, Wing T, Lassiter K, et al. Gene expression in breast muscle associated with feed efficiency in a single male broiler line using a chicken 44 K oligo microarray. I. Top differentially expressed genes. Poult Sci. 2011;90:2535–47.
Bottje WG, Kong B-W, Song JJ, Lee JY, Hargis BM, Lassiter K, et al. Gene expression in breast muscle associated with feed efficiency in a single male broiler line using a chicken 44K microarray. II. Differentially expressed focus genes. Poult Sci. 2012;19:2576–87.
Shendure J. The beginning of the end for microarrays? Nat Methods. 2008;5:585–7.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18:1509–17.
Wang C, Gong B, Bushel PR, Thierry-Mieg J, Thierry-Mieg D, Xu J, et al. The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance. Nat Biotechnol. 2014;32:926–32.
Babraham Bioinformatics. [http://www.bioinformatics.babraham.ac.uk/projects/fastqc/]
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.
Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008;26:317–25.
Inc. SI. JMP® 11 Multivariate Methods. Cary: NC SAS Inst Inc; 2013.
Chargé SBP, Rudnicki MA. Cellular and molecular regulation of muscle regeneration. Physiol Rev. 2004;84:209–38.
Gal-Levi R, Leshem Y, Aoki S, Nakamura T, Halevy O. Hepatocyte growth factor plays a dual role in regulating skeletal muscle satellite cell proliferation and differentiation. Biochim Biophys Acta. 1998;1402:39–51.
Kaliman P, Canicio J, Testar X, Palacín M, Zorzano A. Insulin-like growth factor-II, phosphatidylinositol 3-kinase, nuclear factor-kappaB and inducible nitric-oxide synthase define a common myogenic signaling pathway. J Biol Chem. 1999;274:17437–44.
Zanou N, Gailly P. Skeletal muscle hypertrophy and regeneration: interplay between the myogenic regulatory factors (MRFs) and insulin-like growth factors (IGFs) pathways. Cell Mol Life Sci. 2013;70:4117–30.
Kong Y, Flick MJ, Kudla AJ, Konieczny SF. Muscle LIM protein promotes myogenesis by enhancing the activity of MyoD. Mol Cell Biol. 1997;17:4750–60.
Demonbreun AR, Lapidos KA, Heretis K, Levin S, Dale R, Pytel P, et al. Myoferlin regulation by NFAT in muscle injury, regeneration and repair. J Cell Sci. 2010;123:2413–22.
Velleman SG, Song Y, Shin J, McFarland DC. Modulation of turkey myogenic satellite cell differentiation through the shedding of glypican-1. Comp Biochem Physiol - A Mol Integr Physiol. 2013;164:36–43.
Lu H, Shah P, Ennis D, Shinder G, Sap J, Le-Tien H, et al. The differentiation of skeletal muscle cells involves a protein-tyrosine phosphatase-alpha-mediated C-Src signaling pathway. J Biol Chem. 2002;277:46687–95.
Araya R, Eckardt D, Maxeiner S, Krüger O, Theis M, Willecke K, et al. Expression of connexins during differentiation and regeneration of skeletal muscle: functional relevance of connexin43. J Cell Sci. 2005;118:27–37.
Wang L, Wang Y. Molecular characterization, expression patterns and subcellular localization of Myotrophin (MTPN) gene in porcine skeletal muscle. Mol Biol Rep. 2012;39:2733–8.
Cao DJ, Jiang N, Blagg A, Johnstone JL, Gondalia R, Oh M, et al. Mechanical unloading activates FoxO3 to trigger Bnip3-dependent cardiomyocyte atrophy. J Am Heart Assoc. 2013;2:e000016.
Kipreos ET, Pagano M. Protein family review The F-box protein family. Genome Biol. 2000;1:R3002.
Sandri M, Sandri C, Gilbert A, Skurk C, Calabria E, Picard A, et al. Foxo transcription factors induce the atrophy-related ubiquitin ligase atrogin-1 and cause skeletal muscle atrophy. Cell. 2004;117:399–412.
Skurk C, Izumiya Y, Maatz H, Razeghi P, Shiojima I, Sandri M, et al. The FOXO3a transcription factor regulates cardiac myocyte size downstream of AKT signaling. J Biol Chem. 2005;280:20814–23.
Gomes MD, Lecker SH, Jagoe RT, Navon A, Goldberg AL. Atrogin-1, a muscle-specific F-box protein highly expressed during muscle atrophy. Proc Natl Acad Sci U S A. 2001;98:14440–5.
Ye J, Zhang Y, Xu J, Zhang Q, Zhu D. FBXO40, a gene encoding a novel muscle-specific F-box protein, is upregulated in denervation-related muscle atrophy. Gene. 2007;404:53–60.
Serrano AL, Muñoz-Cánoves P. Regulation and dysregulation of fibrosis in skeletal muscle. Exp Cell Res. 2010;316:3050–8.
Chen X, Li Y. Role of matrix metalloproteinases in skeletal muscle. cell Adhes Migr. 2009;3:337–41.
Carmeli E, Moas M, Reznick AZ, Coleman R. Matrix metalloproteinases and skeletal muscle: a brief review. Muscle Nerve. 2004;29:191–7.
Cornelison DDW. Context matters: in vivo and in vitro influences on muscle satellite cell activity. J Cell Biochem. 2008;105:663–9.
Brinckerhoff CE, Matrisian LM. Matrix metalloproteinases: a tail of a frog that became a prince. Nat Rev Mol Cell Biol. 2002;3:207–14.
Wang W, Pan H, Murray K, Jefferson BS, Li Y. Matrix metalloproteinase-1 promotes muscle cell migration and differentiation. Am J Pathol. 2009;174:541–9.
Wu N, Jansen ED, Davidson JM. Comparison of mouse matrix metalloproteinase 13 expression in free-electron laser and scalpel incisions during wound healing. J Invest Dermatol. 2003;121:926–32.
Kherif S, Lafuma C, Dehaupas M, Lachkar S, Fournier JG, Verdière-Sahuqué M, et al. Expression of matrix metalloproteinases 2 and 9 in regenerating skeletal muscle: a study in experimentally injured and mdx muscles. Dev Biol. 1999;205:158–70.
Rullman E, Rundqvist H, Wågsäter D, Fischer H, Eriksson P, Sundberg CJ, et al. A single bout of exercise activates matrix metalloproteinase in human skeletal muscle. J Appl Physiol. 2007;102:2346–51.
Baum O, Ganster M, Baumgartner I, Nieselt K, Djonov V. Basement membrane remodeling in skeletal muscles of patients with limb ischemia involves regulation of matrix metalloproteinases and tissue inhibitor of matrix metalloproteinases. J Vasc Res. 2007;44:202–13.
Li H, Mittal A, Makonchuk DY, Bhatnagar S, Kumar A. Matrix metalloproteinase-9 inhibition ameliorates pathogenesis and improves skeletal muscle regeneration in muscular dystrophy. Hum Mol Genet. 2009;18:2584–98.
Raffaello A, Milan G, Masiero E, Carnio S, Lee D, Lanfranchi G, et al. JunB transcription factor maintains skeletal muscle mass and promotes hypertrophy. J Cell Biol. 2010;191:101–13.
Braun T, Gautel M. Transcriptional mechanisms regulating skeletal muscle differentiation, growth and homeostasis. Nat Rev Mol Cell Biol. 2011;12:349–61.
Potthoff MJ, Arnold MA, McAnally J, Richardson JA, Bassel-Duby R, Olson EN. Regulation of skeletal muscle sarcomere integrity and postnatal muscle function by Mef2c. Mol Cell Biol. 2007;27:8143–51.
Hinits Y, Hughes SM. Mef2s are required for thick filament formation in nascent muscle fibres. Development. 2007;134:2511–9.
Molkentin JD, Black BL, Martin JF, Olson EN. Cooperative activation of muscle gene expression by MEF2 and myogenic bHLH proteins. Cell. 1995;83:1125–36.
Frey N, Barrientos T, Shelton JM, Frank D, Rütten H, Gehring D, et al. Mice lacking calsarcin-1 are sensitized to calcineurin signaling and show accelerated cardiomyopathy in response to pathological biomechanical stress. Nat Med. 2004;10:1336–43.
Wan L, Ma J, Wang N, Wang D, Xu G. Molecular cloning and characterization of different expression of MYOZ2 and MYOZ3 in Tianfu goat. PLoS One. 2013;8:e82550.
Jørgensen JOL, Jessen N, Pedersen SB, Vestergaard E, Gormsen L, Lund SA, et al. GH Receptor Signaling in Skeletal Muscle and Adipose Tissue in Human Subjects Following Exposure to an Intravenous GH Bolus. Am J Physiol Endocrinol Metab. 2006;291:E899–905.
Goodyer CG, Figueiredo RM, Krackovitch S, De Souza LL, Manalo JA, Zogopoulos G. Characterization of the growth hormone receptor in human dermal fibroblasts and liver during development. Am J Physiol Endocrinol Metab. 2001;281:E1213–20.
Carter-Su C, Schwartz J, Smit LS. Molecular mechanism of growth hormone action. Annu Rev Physiol. 1996;58:187–207.
Hull KL, Harvey S. Autoregulation of central and peripheral growth hormone receptor mRNA in domestic fowl. J Endocrinol. 1998;156:323–9.
Mao JNC, Cogburn LA, Burnside J. Growth hormone down-regulates growth hormone receptor mRNA in chickens but developmental increases in growth hormone receptor mRNA occur independently of growth hormone action. Mol Cell Endocrinol. 1997;129:135–43.
Yamada PM, Lee K-W. Perspectives in mammalian IGFBP-3 biology: local vs. systemic action. Am J Physiol Cell Physiol. 2009;296:C954–76.
Schiaffino S, Mammucari C. Regulation of skeletal muscle growth by the IGF1-Akt/PKB pathway: insights from genetic models. Skelet Muscle. 2011;1:4.
Duclos MJ. Insulin-like growth factor-I (IGF-1) mRNA levels and chicken muscle growth. J Physiol Pharmacol. 2005;56 Suppl 3:25–35.
Fuentes EN, Björnsson BT, Valdés JA, Einarsdottir IE, Lorca B, Alvarez M, et al. IGF-I/PI3K/Akt and IGF-I/MAPK/ERK pathways in vivo in skeletal muscle are regulated by nutrition and contribute to somatic growth in the fine flounder. Am J Physiol Regul Integr Comp Physiol. 2011;300:R1532–42.
Ogawa W, Matozaki T, Kasuga M. Role of binding proteins to IRS-1 in insulin signalling. Mol Cell Biochem. 1998;182:13–22.
Saltiel A, Kahn C. Insulin signalling and the regulation of glucose and lipid metabolism. Nature. 2001;414:799–806.
Bottje WG, Kong BW, Lee JY, Washington T, Baum J. Potential Roles of mTOR and Protein Degradation Pathways in the Phenotypic Expression of Feed Efficiency in Broilers. Biochem Physiol Open Access. 2014;03:1–8.
Pedersen BK, Akerström TCA, Nielsen AR, Fischer CP. Role of myokines in exercise and metabolism. J Appl Physiol. 2007;103:1093–8.
Pedersen BK. Muscles and their myokines. J Exp Biol. 2011;214:337–46.
Nara N, Nakayama Y, Okamoto S, Tamura H, Kiyono M, Muraoka M, et al. Disruption of CXC motif chemokine ligand-14 in mice ameliorates obesity-induced insulin resistance. J Biol Chem. 2007;282:30794–803.
Webster EL, Torpy DJ, Elenkov IJ, Chrousos GP. Corticotropin-releasing hormone and inflammation. Ann N Y Acad Sci. 1998;840:21–32.
Solinas G, Summermatter S, Mainieri D, Gubler M, Montani JP, Seydoux J, et al. Corticotropin-releasing hormone directly stimulates thermogenesis in skeletal muscle possibly through substrate cycling between de novo lipogenesis and lipid oxidation. Endocrinology. 2006;147:31–8.
Lu H, Huang D, Ransohoff RM, Zhou L. Acute skeletal muscle injury: CCL2 expression by both monocytes and injured muscle is required for repair. FASEB J. 2011;25:3344–55.
Lu H, Huang D, Saederup N, Charo IF, Ransohoff RM, Zhou L. Macrophages recruited via CCR2 produce insulin-like growth factor-1 to repair acute skeletal muscle injury. FASEB J. 2011;25:358–69.
Tidball JG. Inflammatory processes in muscle injury and repair. Am J Physiol Regul Integr Comp Physiol. 2005;288:R345–53.
Pillon NJ, Bilan PJ, Fink LN, Klip A. Cross-talk between skeletal muscle and immune cells: muscle-derived mediators and metabolic implications. Am J Physiol Endocrinol Metab. 2013;304:E453–65.
Philippou A, Maridaki M, Theos A, Koutsilieris M. Cytokines in Muscle Damage. Adv Clin Chem. 2012;58:49–87.
Glass DJ. Skeletal muscle hypertrophy and atrophy signaling pathways. Int J Biochem Cell Biol. 2005;37:1974–84.
Andreucci JJ, Grant D, Cox DM, Tomc LK, Prywes R, Goldhamer DJ, et al. Composition and function of AP-1 transcription complexes during muscle cell differentiation. J Biol Chem. 2002;277:16426–32.
Xia R, Webb JA, Gnall LLM, Cutler K, Abramson JJ. Skeletal muscle sarcoplasmic reticulum contains a NADH-dependent oxidase that generates superoxide. Am J Physiol Cell Physiol. 2003;285:C215–21.
Mofarrahi M, Brandes RP, Gorlach A, Hanze J, Terada LS, Quinn MT, et al. Regulation of proliferation of skeletal muscle precursor cells by NADPH oxidase. Antioxid Redox Signal. 2008;10:559–74.
Bottje WG, Carstens GE. Association of mitochondrial function and feed efficiency in poultry and livestock species. J Anim Sci. 2009;87:E48–63.
Brand MD, Affourtit C, Esteves TC, Green K, Lambert AJ, Miwa S, et al. Mitochondrial superoxide: Production, biological effects, and activation of uncoupling proteins. Free Radic Biol Med. 2004;37:755–67.
Kobayashi A, Kang M-I, Watai Y, Tong KI, Shibata T, Uchida K, et al. Oxidative and electrophilic stresses activate Nrf2 through inhibition of ubiquitination activity of Keap1. Mol Cell Biol. 2006;26:221–9.
D’Autréaux B, Toledano MB. ROS as signalling molecules: mechanisms that generate specificity in ROS homeostasis. Nat Rev Mol Cell Biol. 2007;8:813–24.
He X, Kan H, Cai L, Ma Q. Nrf2 is critical in defense against high glucose-induced oxidative damage in cardiomyocytes. J Mol Cell Cardiol. 2009;46:47–58.
Laurila JP, Castellone MD, Curcio A, Laatikainen LE, Haaparanta-Solin M, Gronroos TJ, et al. Extracellular superoxide dismutase is a growth regulatory mediator of tissue injury recovery. Mol Ther. 2009;17:448–54.
Kim S-J, Lee S-M. NLRP3 inflammasome activation in D-galactosamine and lipopolysaccharide-induced acute liver failure: role of heme oxygenase-1. Free Radic Biol Med. 2013;65:997–1004.
Wang X, Ling S, Zhao D, Sun Q, Li Q, Wu F, et al. Redox regulation of actin by thioredoxin-1 is mediated by the interaction of the proteins via cysteine 62. Antioxid Redox Signal. 2010;13:565–73.
Salinas AE, Wong MG. Glutathione S-transferases–a review. Curr Med Chem. 1999;6:279–309.
Bottje W, Pumford NR, Ojano-Dirain C, Iqbal M, Lassiter K. Feed efficiency and mitochondrial function. Poult Sci. 2006;85:8–14.
Thannickal VJ, Fanburg BL. Reactive oxygen species in cell signaling. Am J Physiol Lung Cell Mol Physiol. 2000;279:L1005–28.
Barbieri E, Sestili P. Reactive Oxygen Species in Skeletal Muscle Signaling. J Signal Transduct. 2012.
Victor VM, Rocha M, De La Fuente M. Immune cells: Free radicals and antioxidants in sepsis. Int Immunopharmacol. 2004;4:327–47.
Bashan N, Kovsan J, Kachko I, Ovadia H, Rudich A. Positive and Negative Regulation of Insulin Signaling by Reactive Oxygen and Nitrogen Species. Physiol Rev. 2009;89:27–71.
Kumasaka S, Shoji H, Okabe E. Novel mechanisms involved in superoxide anion radical-triggered Ca2+ release from cardiac sarcoplasmic reticulum linked to cyclic ADP-ribose stimulation. Antioxid Redox Signal. 1999;1:55–69.
Okabe E, Kato Y, Sasaki H, Saito G, Hess ML, Ito H. Calmodulin participation in oxygen radical-induced cardiac sarcoplasmic reticulum calcium uptake reduction. Arch Biochem Biophys. 1987;255:464–8.
Roveri A, Coassin M, Maiorino M, Zamburlini A, van Amsterdam FT, Ratti E, et al. Effect of hydrogen peroxide on calcium homeostasis in smooth muscle cells. Arch Biochem Biophys. 1992;297:265–70.
Periasamy M, Kalyanasundaram A. SERCA pump isoforms: their role in calcium transport and disease. Muscle Nerve. 2007;35:430–42.
Ke Q, Costa M. Hypoxia-inducible factor-1 (HIF-1). Mol Pharmacol. 2006;70:1469–80.
Stiehl DP. Ã WJ, Wenger RH, Hellwig-bu T, Allee R, Lu D-: Normoxic induction of the hypoxia-inducible factor 1 K by insulin and interleukin-1 L involves the phosphatidylinositol 3-kinase pathway. FEBS Lett. 2002;512:157–62.
Feldser D, Agani F, Iyer NV. Reciprocal positive regulation of hypoxia-inducible factor 1α and insulin-like growth factor 2 growth factor 2 1. Cancer Res. 1999;59:3915–8.
Eltzschig HK, Carmeliet P. Hypoxia and inflammation. N Engl J Med. 2011;364:656–65.
Nizet V, Johnson RS. Interdependence of hypoxic and innate immune responses. Nat Rev Immunol. 2009;9:609–17.
Cramer T, Yamanishi Y, Clausen BE, Förster I, Pawlinski R, Mackman N, et al. HIF-1α is essential for myeloid cell-mediated inflammation. Cell. 2003;112:645–57.
Bartels K, Grenz A, Eltzschig HK. Hypoxia and inflammation are two sides of the same coin. Proc Natl Acad Sci U S A. 2013;110:18351–2.
Frede S, Berchner-Pfannschmidt U, Fandrey J. Regulation of Hypoxia-Inducible Factors During Inflammation. Methods Enzymol. 2007;435:403–19.
Scheerer N, Dehne N, Stockmann C, Swoboda S, Baba HA, Neugebauer A, et al. Myeloid hypoxia-inducible factor-1α is essential for skeletal muscle regeneration in mice. J Immunol. 2013;191:407–14.
Chandel NS, Maltepe E, Goldwasser E, Mathieu CE, Simon MC, Schumacker PT. Mitochondrial reactive oxygen species trigger hypoxia-induced transcription. Proc Natl Acad Sci U S A. 1998;95:11715–20.
Semenza GL. Regulation of mammalian O2 homeostasis by hypoxia-inducible factor 1. Annu Rev Cell Dev Biol. 1999;15:551–78.
Dröge W. Free radicals in the physiological control of cell function. Physiol Rev. 2002;82:47–95.
Chandel NS, McClintock DS, Feliciano CE, Wood TM, Melendez JA, Rodriguez AM, et al. Reactive oxygen species generated at mitochondrial complex III stabilize hypoxia-inducible factor-1alpha during hypoxia: a mechanism of O2 sensing. J Biol Chem. 2000;275:25130–8.
Petracci M, Mudalal S, Bonfiglio A, Cavani C. Occurrence of white striping under commercial conditions and its impact on breast meat quality in broiler chickens. Poult Sci. 2013;92:1670–5.
Kuttappan VA, Shivaprasad H, Shaw DP, Valentine BA, Hargis BM, Clark FD, et al. Pathological changes associated with white striping in broiler breast muscles. Poult Sci. 2013;92:331–8.
Sihvo H-K, Immonen K, Puolanne E. Myodegeneration with fibrosis and regeneration in the pectoralis major muscle of broilers. Vet Pathol. 2014;51:619–23.
This work was partly funded by Delaware Bioscience Center for Advanced Technology and Heritage Breeders, LLC. Thanks to Brewster F. Kingham in Delaware Biotechnology Institute (DBI) for performing RNA sequencing.
The authors declare that they have no competing interests.
NZ constructed RNA-seq libraries, analyzed the RNA-seq data, performed functional and pathway analysis, and drafted the manuscript. BA and WRL designed and carried out the animal experiment. BA conceived of the study, participated in the analysis and helped to draft the manuscript. All authors read and approved the final manuscript.
The number of properly mapped, improperly mapped and unmapped reads is shown for each sample. Unmapped reads are reads not mapped onto the reference genome; Properly mapped reads are paired-end reads mapped to the reference genome and complying with the parameters (−g 1 –r 110 –no-discordant –no-mixed); Improperly mapped reads are reads mapped to the reference genome but not complying with the parameters.
Volcano plot showing differentially expressed genes between high- and low-FE chickens. Genes with q-value less than 0.05 were identified as significantly differentially expressed genes, which were shown in red, while genes with q-value larger than 0.05 were not statistically significant, shown in black.
1059 significant genes detected in the present study.
Comparison of gene expression between high- and low-FE broiler chickens. Genes up-regulated in the high- and low-FE chickens were marked with blue and red spots, respectively.
Top 10 up-regulated and down-regulated genes in high-FE chickens.
Activated or inhibited upstream transcription regulators predicted by Ingenuity® Pathway Analysis (IPA) Software.
Differentially expressed genes involved in muscle development. Ingenuity® Pathway Analysis (IPA) Software identified 31 genes in our dataset that are associated with muscle development based on Ingenuity® Knowledge Base.
Differentially expressed genes from matrix metalloproteinases (MMPs) family.
Over-represented pathways involved in immune response.
Differentially expressed genes associated with free radical scavenging are also involved in inflammatory response.
About this article
Cite this article
Zhou, N., Lee, W.R. & Abasht, B. Messenger RNA sequencing and pathway analysis provide novel insights into the biological basis of chickens’ feed efficiency. BMC Genomics 16, 195 (2015). https://doi.org/10.1186/s12864-015-1364-0
- Differential expression analysis
- Chicken feed efficiency
- Breast muscle
- Muscle remodeling
- IGFs/PI3K/Akt signaling pathway