Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Transcriptome Profiling of the Liver in Nellore Cattle Phenotypically Divergent for RFI in Two Genetic Groups

The identification and selection of genetically superior animals for residual feed intake (RFI) could enhance productivity and minimize environmental impacts. The aim of this study was to use RNA-seq data to identify the differentially expressed genes (DEGs), known non-coding RNAs (ncRNAs), specific biomarkers and enriched biological processes associated with RFI of the liver in Nellore cattle in two genetic groups. In genetic group 1 (G1), 24 extreme RFI animals (12 low RFI (LRFI) versus 12 high RFI (HRFI)) were selected from a population of 60 Nellore bulls. The RNA-seq of the samples from their liver tissues was performed using an Illumina HiSeq 2000. In genetic group 2 (G2), 20 samples of liver tissue of Nellore bulls divergent for RFI (LRFI, n = 10 versus HRFI, n = 10) were selected from 83 animals. The raw data of the G2 were chosen from the ENA repository. A total of 1811 DEGs were found for the G1 and 2054 for the G2 (p-value ≤ 0.05). We detected 88 common genes in both genetic groups, of which 33 were involved in the immune response and in blocking oxidative stress. In addition, seven (B2M, ADSS, SNX2, TUBA4A, ARHGAP18, MECR, and ABCF3) possible gene biomarkers were identified through a receiver operating characteristic analysis (ROC) considering an AUC > 0.70. The B2M gene was overexpressed in the LRFI group. This gene regulates the lipid metabolism protein turnover and inhibits cell death. We also found non-coding RNAs in both groups. MIR25 was up-regulated and SNORD16 was down-regulated in the LRFI for G1. For G2, up-regulated RNase_MRP and SCARNA10 were found. We highlight MIR25 as being able to act by blocking cytotoxicity and oxidative stress and RMRP as a blocker of mitochondrial damage. The biological pathways associated with RFI of the liver in Nellore cattle in the two genetic groups were for energy metabolism, protein turnover, redox homeostasis and the immune response. The common transcripts, biomarkers and metabolic pathways found in the two genetic groups make this unprecedented work even more relevant, since the results are valid for different herds raised in different ways. The results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the liver tissue transcriptome of Nellore cattle.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Reduced PGC-1β protein expression may underlie corticosterone inhibition of mitochondrial biogenesis and oxidative phosphorylation in chicken muscles

To uncover the molecular mechanism underlying glucocorticoid-induced loss of mitochondrial integrity in skeletal muscles, studies were performed to investigate whether the peroxisome proliferator-activated receptor γ coactivator 1 (PGC-1)-mediated pathway was involved in this process. In an in vivo trial, 3 groups of 30-d-old Arbor Acres male broilers were randomly subjected to one of the following treatments for 7 days: corticosterone (CORT, 30 mg/kg diet), control (blank), and pair-feeding (restricted to the same feed intake as for the CORT treatment), each with 6 replicates of 15 birds. Mitochondrial abundance, morphology, and function were determined in the pectoralis major and biceps femoris muscles. In an in vitro trial, a primary culture of embryonic chick myotubes was incubated with a serum-free medium for 24 h in the presence or absence of CORT (0, 200, and 1,000 nM). Results showed that CORT destroyed mitochondrial ultrastructure (p < 0.01), and decreased the enzymatic activity and protein expression of respiratory chain complexes (p < 0.05), leading to an inferior coupling efficiency (p < 0.05). As reflected by a decline in mitochondrial density (p < 0.01) and mitochondrial DNA copy number (p < 0.05), CORT reduced mitochondrial contents. Among all three PGC-1 family members, only PGC-1β was down-regulated by CORT at the protein level (p < 0.05). Some aspects of these responses were tissue-specific and seemed to result from the depressed feed intake. Overall, CORT may impair mitochondrial biogenesis and oxidative phosphorylation in a PGC-1β-dependent manner in chicken muscles.

Proper Immune Response Depends on Early Exposure to Gut Microbiota in Broiler Chicks

The successional changes in the early intestinal microbiota occur concomitantly with the development, expansion, and education of the mucosal immune system. Although great attention of researchers has been focused on understanding the linkage between microbiota and immune functions, many essential details of the symbiotic relationship between the intestinal pioneer microbiota and the avian immune system remain to be discovered. This study was conducted to understand the impact of different early life intestinal colonizers on innate and adaptive immune processes in chicks and further identify immune-associated proteins expressed in the intestinal tissue. To accomplish it, we performed an in ovo application of two apathogenic Enterobacteriaceae isolates and lactic acid bacteria (L) to determine their influences on the intestinal proteome profile of broilers at the day of hatch (DOH) and at 10 days old. The results indicated that there were predicted biological functions of L-treated chicks associated with the activation and balanced function of the innate and adaptive immune systems. At the same time, the Enterobacteriaceae-exposed birds presented dysregulated immunological mechanisms or downregulated processes related to immune development. Those findings suggested that a proper immune function was dependent on specific gut microbiota exposure, in which the prenatal probiotic application may have favored the fitting programming of immune functions in chicks.

An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake

The selection of cattle with enhanced feed efficiency is of importance with regard to reducing feed costs in the beef industry. Global transcriptome profiling was undertaken on liver and skeletal muscle biopsies from Simmental heifers and bulls divergent for residual feed intake (RFI), a widely acknowledged feed efficiency phenotype, in order to identify genes that may be associated with this trait. We identified 5 genes (adj. p < 0.1) to be differentially expressed in skeletal muscle between high and low RFI heifers with all transcripts involved in oxidative phosphorylation and mitochondrial homeostasis. A total of 11 genes (adj. p < 0. 1) were differentially expressed in liver tissue between high and low RFI bulls with differentially expressed genes related to amino and nucleotide metabolism as well as endoplasmic reticulum protein processing. No genes were identified as differentially expressed in either heifer liver or bull muscle analyses. Results from this study show that the molecular control of RFI in young cattle is modified according to gender, which may be attributable to differences in physiological maturity between heifers and bulls of the same age. Despite this we have highlighted a number of genes that may hold potential as molecular biomarkers for RFI cattle.

Upstream Regulator Analysis of Wooden Breast Myopathy Proteomics in Commercial Broilers and Comparison to Feed Efficiency Proteomics in Pedigree Male Broilers

In an effort to understand the apparent trade-off between the continual push for growth performance and the recent emergence of muscle pathologies, shotgun proteomics was conducted on breast muscle obtained at ~8 weeks from commercial broilers with wooden breast (WB) myopathy and compared with that in pedigree male (PedM) broilers exhibiting high feed efficiency (FE). Comparison of the two proteomic datasets was facilitated using the overlay function of Ingenuity Pathway Analysis (IPA) (Qiagen, CA, USA). We focused on upstream regulator analysis and disease-function analysis that provides predictions of activation or inhibition of molecules based on (a) expression of downstream target molecules, (b) the IPA scientific citation database. Angiopoeitin 2 (ANGPT2) exhibited the highest predicted activation Z-score of all molecules in the WB dataset, suggesting that the proteomic landscape of WB myopathy would promote vascularization. Overlaying the FE proteomics data on the WB ANGPT2 upstream regulator network presented no commonality of protein expression and no prediction of ANGPT2 activation. Peroxisome proliferator coactivator 1 alpha (PGC1α) was predicted to be inhibited, suggesting that mitochondrial biogenesis was suppressed in WB. PGC1α was predicted to be activated in high FE pedigree male broilers. Whereas RICTOR (rapamycin independent companion of mammalian target of rapamycin) was predicted to be inhibited in both WB and FE datasets, the predictions were based on different downstream molecules. Other transcription factors predicted to be activated in WB muscle included epidermal growth factor (EGFR), X box binding protein (XBP1), transforming growth factor beta 1 (TGFB1) and nuclear factor (erythroid-derived 2)-like 2 (NFE2L2). Inhibitions of aryl hydrocarbon receptor (AHR), AHR nuclear translocator (ARNT) and estrogen related receptor gamma (ESRRG) were also predicted in the WB muscle. These findings indicate that there are considerable differences in upstream regulators based on downstream protein expression observed in WB myopathy and in high FE PedM broilers that may provide additional insight into the etiology of WB myopathy.

Identification of key genes and pathways associated with feed efficiency of native chickens based on transcriptome data via bioinformatics analysis

BACKGROUND

Improving feed efficiency is one of the important breeding targets for poultry industry. The aim of current study was to investigate the breast muscle transcriptome data of native chickens divergent for feed efficiency. Residual feed intake (RFI) value was calculated for 1008 closely related chickens. The 5 most efficient (LRFI) and 5 least efficient (HRFI) birds were selected for further analysis. Transcriptomic data were generated from breast muscle collected post-slaughter.

RESULTS

The differently expressed genes (DEGs) analysis showed that 24 and 325 known genes were significantly up- and down-regulated in LRFI birds. An enrichment analysis of DEGs showed that the genes and pathways related to inflammatory response and immune response were up-regulated in HRFI chickens. Moreover, Gene Set Enrichment Analysis (GSEA) was also employed, which indicated that LRFI chickens increased expression of genes related to mitochondrial function. Furthermore, protein network interaction and function analyses revealed ND2, ND4, CYTB, RAC2, VCAM1, CTSS and TLR4 were key genes for feed efficiency. And the 'phagosome', 'cell adhesion molecules (CAMs)', 'citrate cycle (TCA cycle)' and 'oxidative phosphorylation' were key pathways contributing to the difference in feed efficiency.

CONCLUSIONS

In summary, a series of key genes and pathways were identified via bioinformatics analysis. These key genes may influence feed efficiency through deep involvement in ROS production and inflammatory response. Our results suggested that LRFI chickens may synthesize ATP more efficiently and control reactive oxygen species (ROS) production more strictly by enhancing the mitochondrial function in skeletal muscle compared with HRFI chickens. These findings provide some clues for understanding the molecular mechanism of feed efficiency in birds and will be a useful reference data for native chicken breeding.

A Proteomic View of the Cross-Talk Between Early Intestinal Microbiota and Poultry Immune System

Proteomics has been used to investigate cross-talk between the intestinal microbiome and host biological processes. In this study, an in ovo technique and a proteomics approach was used to address how early bacterial colonization in the gastrointestinal tract (GIT) could modulate inflammatory and immune responses in young broilers. Embryos at 18 embryogenic days were inoculated with saline (S), 102 CFU of Citrobacter freundii (CF), Citrobacter species (C2), or lactic acid bacteria mixture (L) into the amnion. At 10 days posthatch, ileum samples from 12 birds per treatment were selected for tandem mass spectrometry analysis. Our further findings indicated that treatment-specific influences on early GIT microbiota resulted in different immune responses in mature broilers. Predicted functional analyses revealed activation of inflammation pathways in broilers treated in ovo with L and CF. Exposure to L enhanced functional annotation related to activation, trafficking of immune cells, and skeletal growth based-network, while CF inhibited biological functions associated with immune cell migration and inflammatory response. These results highlighted that proper immune function was dependent on specific GIT microbiota profiles, in which early-life exposure to L-based probiotic may have modulated the immune functions, whereas neonatal colonization of Enterobacteriaceae strains may have led to immune dysregulation associated with chronic inflammation.

Impact of in ovo administered pioneer colonizers on intestinal proteome on day of hatch

Pioneer colonization of the gastrointestinal tract (GIT) by bacteria is thought to have major influence on neonatal tissue development. Previous studies have shown in ovo inoculation of embryos with saline (S), species of Citrobacter (C, C2), or lactic acid bacteria (L) resulted in an altered microbiome on day of the hatch (DOH). The present study investigated GIT proteomic changes at DOH in relation to different inoculations. Embryos were inoculated in ovo with S or ∼10² cfu of C, C2, or L at 18 embryonic days. On DOH, the GIT was collected, and tissue proteins were extracted for analysis via tandem mass spectrometry. A total of 493 proteins were identified for differential comparison with S at P ≤ 0.10. Different levels were noted in 107, 39, and 78 proteins in C, C2, and L groups, respectively, which were uploaded to Ingenuity Pathway Analysis to determine canonical pathways and biological functions related to these changes. Three members of the cytokine family (interleukin [IL]-1β, IL6, and Oncostatin M) were predicted to be activated in C2, indicated with Z-score ≥ 1.50, which suggested an overall proinflammatory GIT condition. This was consistent with the activation of the acute-phase response signaling pathway seen exclusively in C2 (Z-score = 2.00, P < 0.01). However, activation (Z-score = 2.00) of IL-13, upregulation of peroxiredoxin-1 and superoxide dismutase 1, in addition to activation of nitric oxide signaling in the cardiovascular system of the L treatment may predict a state of increased antioxidant capacity and decreased inflammatory status. The nuclear factor erythroid 2-related factor 2 (NRF2)-mediated oxidative stress response (Z-score = 2.00, P < 0.01) was predicted to be upregulated in C which suggested that chicks were in an inflammatory state and associated oxidative stress, but the impact of these pathways differed from that of C2. These changes in the proteome suggest that pioneer colonizing microbiota may have a strong impact on pathways associated with GIT immune and cellular development.

BOARD INVITED REVIEW: Oxidative stress and efficiency: the tightrope act of mitochondria in health and disease1,2

Oxidative stress is an unavoidable consequence of aerobic metabolism. Whereas high amounts of mitochondrial reactive oxygen species (ROS) can cause oxidation, low levels play important roles in signal transduction. In a Pedigree male (PedM) broiler model of feed efficiency (FE), the low FE phenotype was characterized by increased ROS in isolated mitochondria (muscle, liver, and duodenum) with a pervasive protein oxidation in mitochondria and tissues. Subsequent proteogenomic studies in muscle revealed evidence of enhanced mitoproteome abundance, enhanced mitochondrial phosphocreatine shuttling expression, and enhanced ribosome assembly in the high FE phenotype. Surprisingly, an enhanced infrastructure would foster greater repair of damaged proteins or organelles through the autophagy and proteosome pathways in the high FE phenotype. Although protein and organelle degradation, recycling, and reconstruction would be energetically expensive, it is possible that energy invested into maintaining optimal function of proteins and organelles contributes to cellular efficiency in the high FE phenotype. New findings in mitochondrial physiology have been reported in the last several years. Reverse electron transport (RET), once considered an artifact of in vitro conditions, now is recognized to play significant roles in inflammation, ischemia-reperfusion, muscle differentiation, and energy utilization. A topology of ROS production indicates that ROS derived from Complex I of the respiratory chain primarily causes oxidation, whereas ROS generated from Complex III are primarily involved in cell signaling. It is also apparent that there is a constant fission and fusion process that mitochondria undergo that help maintain optimal mitochondrial function and enables mitochondria to adjust to periods of nutrient limitation and nutrient excess. Understanding the balancing act that mitochondria play in health and disease will continue to be a vital biological component in health-production efficiency and disease in commercial animal agriculture.

Tissue-specific responses of antioxidant pathways to poor hygiene conditions in growing pigs divergently selected for feed efficiency

Background

Poor hygiene of housing induces a systemic inflammatory response. Because inflammation and oxidative stress are processes that can sustain each other, the ways pigs are able to activate their antioxidant defenses are critical for production performance and health during periods when the immune system is solicited. Selection for production performance can also influence reactive oxygen species (ROS) production and expression levels of genes involved in cellular response to oxidative stress in different tissues. To establish the extent by which poor hygiene and selection for feed efficiency affected redox status, pigs divergently selected for residual feed intake (RFI) were housed in poor or good hygiene during 6 weeks. At the end, blood was collected in all pigs, and half of them were killed for tissue sampling. The remaining pigs were reared in good hygiene conditions during a recovery period of 7–8 weeks.

Results

At week 6, poor hygiene was associated with a lower total antioxidant capacity assessed by plasma ferric reducing ability in all pigs, and with greater plasma levels of hydrogen peroxides in the high RFI pigs (less efficient). Adipose tissue of high RFI pigs exhibited higher activities of catalase and glutathione reductase, and greater thiobarbituric acid reactive substances (TBARS) concentrations when compared with the low RFI pigs (more efficient). Poor hygiene conditions activated the antioxidant enzymes activities (glutathione reductase, superoxide dismutase and catalase) in adipose tissue of both lines, but led to higher ROS production by mature adipocytes isolated from the high RFI pigs only. In liver and muscle, there were only minor changes in antioxidant molecules due to genetics and hygiene conditions. After the resilience period, adipose tissue of pigs previously challenged by poor hygiene maintained higher antioxidant enzyme activities, and for the high RFI line, displayed higher TBARS concentrations.

Conclusions

Pigs selected for improved feed efficiency showed a lower susceptibility to oxidative stress induced by poor hygiene conditions. This could led to a lower inflammatory response and less impaired growth when these pigs are facing sanitary challenges during the production period.

The metabolic characteristics of susceptibility to wooden breast disease in chickens with high feed efficiency

This study was conducted to characterize metabolic differences between high feed efficiency (HFE) and low feed efficiency (LFE) chickens to investigate why feed efficient chickens are more susceptible to muscle abnormalities such as wooden breast disease. Gene expression profiles were generated by RNA sequencing of pectoralis major muscle samples from 10 HFE and 13 LFE broiler chickens selected from a modern broiler population. Metabolism-associated differentially expressed genes were identified and interpreted by Ingenuity Pathway Analysis and literature mining. Our RNA-seq data indicate decreased glycolytic capacity, increased fatty acid uptake, mitochondrial oxidation of fatty acids, and several other metabolic alterations in the pectoralis major muscle of HFE chickens. We also quantified glycogen content of the pectoralis major muscle and found that the HFE chickens had a significantly (P ≤ 0.05) lower glycogen content. Collectively, this study indicates extensive metabolic differences in the pectoralis major muscle between HFE and LFE chickens and helps identify metabolic features of susceptibility to muscle disorders in modern broiler chickens.

Differential hepatic oxidative status in steers with divergent residual feed intake phenotype

Oxidative stress occurs when oxidant production exceeds the antioxidant capacity to detoxify the reactive intermediates or to repair the resulting damage. Feed efficiency has been associated with mitochondrial function due to its impact on cell energy metabolism. However, mitochondria are also recognized as a major source of oxidants. The aim of this study was to determine lipid and protein oxidative stress markers, and gene and protein expression as well as activity of antioxidant enzymes in the liver of steers of divergent residual feed intake (RFI) phenotypes. Hereford steers (n = 111) were evaluated in post-weaning 70 days standard test for RFI. Eighteen steers exhibiting the greatest (n = 9; high-RFI) and the lowest (n = 9; low-RFI) RFI values were selected for this study. After the test, steers were managed together under grazing conditions until slaughter when they reached the slaughter body weight. At slaughter, hepatic samples were obtained, were snap-frozen in liquid nitrogen and stored at -80°C until analyses. Hepatic thiobarbituric acid reactive species and protein carbonyls were greater (P = 0.05) and hepatic 4-hydroxynonenal protein adducts tended (P = 0.10) to be greater for high- than low-RFI steers. Hepatic gene expression glutathione peroxidase 4, glutamate-cysteine ligase catalytic subunit and peroxiredoxin 5 mRNA was greater (P ≤ 0.05) and glutathione peroxidase 3 mRNA tended (P = 0.10) to be greater in low- than high-RFI steers. Hepatic protein expression and enzyme activity of manganese superoxide dismutase and glutathione peroxidase enzyme activity tended (P ≤ 0.10) to be greater for low- than high-RFI steers. High-efficiency steers (low-RFI) probably had better hepatic oxidative status which was strongly associated with greater antioxidant ability near to the oxidant production site and, therefore, reduced oxidative stress of the liver. Decreased hepatic oxidative stress would reduce maintenance requirements due to a lower protein and lipid turnover and better efficiency in the use of energy.

Gene Expression Essential for Myostatin Signaling and Skeletal Muscle Development Is Associated With Divergent Feed Efficiency in Pedigree Male Broilers

Background: Feed efficiency (FE, gain to feed) is an important genetic trait as 70% of the cost of raising animals is due to feed costs. The objective of this study was to determine mRNA expression of genes involved in muscle development and hypertrophy, and the insulin receptor-signaling pathway in breast muscle associated with the phenotypic expression of FE. Methods: Breast muscle samples were obtained from Pedigree Male (PedM) broilers (8 to 10 week old) that had been individually phenotyped for FE between 6 and 7 week of age. The high FE group gained more weight but consumed the same amount of feed compared to the low FE group. Total RNA was extracted from breast muscle (n = 6 per group) and mRNA expression of target genes was determined by real-time quantitative PCR. Results: Targeted gene expression analysis in breast muscle of the high FE phenotype revealed that muscle development may be fostered in the high FE PedM phenotype by down-regulation several components of the myostatin signaling pathway genes combined with upregulation of genes that enhance muscle formation and growth. There was also evidence of genetic architecture that would foster muscle protein synthesis in the high FE phenotype. A clear indication of differences in insulin signaling between high and low FE phenotypes was not apparent in this study. Conclusion: These findings indicate that a gene expression architecture is present in breast muscle of PedM broilers exhibiting high FE that would support enhanced muscle development-differentiation as well as protein synthesis compared to PedM broilers exhibiting low FE.

Feed Restriction Reveals Distinct Serum Metabolome Profiles in Chickens Divergent in Feed Efficiency Traits

Restrictive feeding influences systemic metabolism of nutrients; however, this impact has not been evaluated in chickens of diverging feed efficiency. This study investigated the effect of ad libitum versus restrictive feeding (85% of ad libitum) on the serum metabolome and white blood cell composition in chickens of diverging residual feed intake (RFI; metric for feed efficiency). Blood samples were collected between days 33 and 37 post-hatch. While serum glucose was similar, serum uric acid and cholesterol were indicative of the nutritional status and chicken’s RFI, respectively. Feed restriction and RFI rank caused distinct serum metabolome profiles, whereby restrictive feeding also increased the blood lymphocyte proportion. Most importantly, 10 amino acids were associated with RFI rank in birds, whereas restrictive feeding affected almost all detected lysophosphatidylcholines, with 3 being higher and 6 being lower in restrictively compared to ad libitum fed chickens. As indicated by relevance networking, isoleucine, lysine, valine, histidine, and ornithine were the most discriminant for high RFI, whereas 3 biogenic amines (carnosine, putrescine, and spermidine) and 3 diacyl-glycerophospholipids (38:4, 38:5, and 40:5) positively correlated with feed intake and body weight gain, respectively. Only for taurine, feed intake mostly explained the RFI-associated variation, whereas for most metabolites, other host physiological factors played a greater role for the RFI-associated differences, and was potentially related to insulin-signaling, phospholipase A2, and arachidonic acid metabolism. Alterations in the hepatic synthesis of long-chain fatty acids and the need for precursors for gluconeogenesis due to varying energy demand may explain the marked differences in serum metabolite profiles in ad libitum and restrictively fed birds.

Enrichment of Autophagy and Proteosome Pathways in Breast Muscle of Feed Efficient Pedigree Male Broilers

Background: Feed efficiency (FE) is an important genetic trait in poultry and livestock. Autophagy (self-eating) and proteosomes are cellular processes that remove damaged cell components (e.g., proteins, organelles). As evidence of extensive protein oxidation was observed in Pedigree Male (PedM) broilers exhibiting a low FE (LFE) phenotype compared to a high FE (HFE) phenotype, the main goal of this study was to assess gene and protein expression of the autophagy and proteosome pathways in breast muscle obtained in PedM broilers exhibiting HFE and LFE phenotypes.

Methods: Feed efficiency was calculated as weight gain divided by feed intake gain in individual PedM broilers that were measured between 6 and 7 weeks of age. Targeted gene expression was conducted on breast muscle using quantitative real-time polymerase chain reaction (qPCR) to determine mRNA expression of genes associated with the autophagy pathway; AMP-activated protein kinase alpha 1 (AMPKα1), mammalian target of rapamycin (mTOR), Beclin 1, and autophagy genes (Atg) 3, Atg7, and Atg16L1. Binomial distribution analysis was conducted on transcriptomic and data obtained by RNAseq and shotgun proteomics, respectively on the same set of tissues for genes associated with autophagy, vacuole formation, and proteosome expression.

Results: Greater efficiency was attained in the HFE PedM broilers by greater weight gain on the same amount of feed consumed resulting in FEs of 0.65 ± 0.01 and 0.46 ± 0.01 in the HFE and LFE phenotypes, respectively. Targeted mRNA expression analysis revealed significant (P < 0.05) elevations in AMPKa1, mTOR, Atg16L1, and Atg7 and a marginal (P = 0.07) elevation in Beclin1. Binomial distribution analysis transcriptomic and proteomic data revealed significant skews favoring autophagy-, vacuole-, and proteosome-related genes in the HFE phenotype. These results indicate that the autophagy and proteosome expression is enhanced in the HFE compared to the LFE pedigree male broiler phenotype suggesting that protein and organelle quality control may be enhanced in high feed efficiency.

MicroRNA profiling associated with muscle growth in modern broilers compared to an unselected chicken breed

Background

Genetically selected modern broiler chickens have acquired outstanding production efficiency through rapid growth and improved feed efficiency compared to unselected chicken breeds. Recently, we analyzed the transcriptome of breast muscle tissues obtained from modern pedigree male (PeM) broilers (rapid growth and higher efficiency) and foundational Barred Plymouth Rock (BPR) chickens (slow growth and poorer efficiency). This study was designed to investigate microRNAs that play role in rapid growth of the breast muscles in modern broiler chickens.

Results

In this study, differential abundance of microRNA (miRNA) was analyzed in breast muscle of PeM and BPR chickens and the results were integrated with differentially expressed (DE) mRNA in the same tissues. A total of 994 miRNA were identified in PeM and BPR chicken lines from the initial analysis of small RNA sequencing data. After filtering and statistical analyses, the results showed miR-2131-5p, miR-221-5p, miR-126-3p, miR-146b-5p, miR-10a-5p, let-7b, miR-125b-5p, and miR-146c-5p up-regulated whereas miR-206 down-regulated in PeM compared to BPR breast muscle. Based on inhibitory regulations of miRNAs on the mRNA abundance, our computational analysis using miRDB, an online software, predicated that 118 down-regulated mRNAs may be targeted by the up-regulated miRNAs, while 35 up-regulated mRNAs appear to be due to a down-regulated miRNA (i.e., miR-206). Functional network analyses of target genes of DE miRNAs showed their involvement in calcium signaling, axonal guidance signaling, and NRF2-mediated oxidative stress response pathways suggesting their involvement in breast muscle growth in chickens.

Conclusion

From the integrated analyses of differentially expressed miRNA-mRNA data, we were able to identify breast muscle specific miRNAs and their target genes whose concerted actions can contribute to rapid growth and higher feed efficiency in modern broiler chickens. This study provides foundation data for elucidating molecular mechanisms that govern muscle growth in chickens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5061-7) contains supplementary material, which is available to authorized users.

Application of omics technologies for a deeper insight into quali-quantitative production traits in broiler chickens: A review

The poultry industry is continuously facing substantial and different challenges such as the increasing cost of feed ingredients, the European Union's ban of antibiotic as growth promoters, the antimicrobial resistance and the high incidence of muscle myopathies and breast meat abnormalities. In the last decade, there has been an extraordinary development of many genomic techniques able to describe global variation of genes, proteins and metabolites expression level. Proper application of these cutting-edge omics technologies (mainly transcriptomics, proteomics and metabolomics) paves the possibility to understand much useful information about the biological processes and pathways behind different complex traits of chickens. The current review aimed to highlight some important knowledge achieved through the application of omics technologies and proteo-genomics data in the field of feed efficiency, nutrition, meat quality and disease resistance in broiler chickens.

Identification and Differential Abundance of Mitochondrial Genome Encoding Small RNAs (mitosRNA) in Breast Muscles of Modern Broilers and Unselected Chicken Breed

Background: Although small non-coding RNAs are mostly encoded by the nuclear genome, thousands of small non-coding RNAs encoded by the mitochondrial genome, termed as mitosRNAs were recently reported in human, mouse and trout. In this study, we first identified chicken mitosRNAs in breast muscle using small RNA sequencing method and the differential abundance was analyzed between modern pedigree male (PeM) broilers (characterized by rapid growth and large muscle mass) and the foundational Barred Plymouth Rock (BPR) chickens (characterized by slow growth and small muscle mass). Methods: Small RNA sequencing was performed with total RNAs extracted from breast muscles of PeM and BPR (n = 6 per group) using the 1 × 50 bp single end read method of Illumina sequencing. Raw reads were processed by quality assessment, adapter trimming, and alignment to the chicken mitochondrial genome (GenBank Accession: X52392.1) using the NGen program. Further statistical analyses were performed using the JMP Genomics 8. Differentially expressed (DE) mitosRNAs between PeM and BPR were confirmed by quantitative PCR. Results: Totals of 183,416 unique small RNA sequences were identified as potential chicken mitosRNAs. After stringent filtering processes, 117 mitosRNAs showing >100 raw read counts were abundantly produced from all 37 mitochondrial genes (except D-loop region) and the length of mitosRNAs ranged from 22 to 46 nucleotides. Of those, abundance of 44 mitosRNAs were significantly altered in breast muscles of PeM compared to those of BPR: all mitosRNAs were higher in PeM breast except those produced from 16S-rRNA gene. Possibly, the higher mitosRNAs abundance in PeM breast may be due to a higher mitochondrial content compared to BPR. Our data demonstrate that in addition to 37 known mitochondrial genes, the mitochondrial genome also encodes abundant mitosRNAs, that may play an important regulatory role in muscle growth via mitochondrial gene expression control.

Enhanced expression of proteins involved in energy production and transfer in breast muscle of pedigree male broilers exhibiting high feed efficiency

In cells with fluctuating energy demand (e.g., skeletal muscle), a transfer system of proteins across the inner and outer mitochondrial membranes links mitochondrial oxidative phosphorylation to cytosolic phosphorylated creatine (PCr) that serves as a phosphate reservoir for rapid repletion of cytosolic adenosine triphosphate (ATP). Crucial proteins of this energy transfer system include several creatine kinase (CK) isoforms found in the cytosol and mitochondria. In a recent proteomic study (Kong et al., 2016), several components of this system were up-regulated in high feed efficiency (FE) compared to low FE breast muscle; notably adenine nucleotide translocase (ANT), voltage dependent activated channel (VDAC), the brain isoform of creatine kinase (CK-B), and several proteins of the electron transport chain. Reexamination of the original proteomic dataset revealed that the expression of two mitochondrial CK isoforms (CKMT1A and CKMT2) had been detected but were not recognized by the bioinformatics program used by Kong et al. (2016a). The CKMT1A isoform was up-regulated (7.8-fold, P = 0.05) in the high FE phenotype but there was no difference in CKMT2 expression (1.1-fold, P = 0.59). From these findings, we hypothesize that enhanced expression of the energy production and transfer system in breast muscle of the high FE pedigree broiler male could be fundamentally important in the phenotypic expression of feed efficiency.

Proteogenomics Reveals Enriched Ribosome Assembly and Protein Translation in Pectoralis major of High Feed Efficiency Pedigree Broiler Males

Background: In production animal agriculture, the cost of feed represents 60–70% of the total cost of raising an animal to market weight. Thus, development of viable biomarkers for feed efficiency (FE, g gain/g feed) to assist in genetic selection of breeding stock remains an important goal in commercial breeding programs.

Methods: Global gene (cDNA microarray, RNAseq) and protein expression (shotgun proteomics) analyses have been conducted on breast muscle samples obtained from pedigree broiler males (PedM) exhibiting high and low FE phenotypes. Using the entire datasets (i.e., no cutoffs for significance or fold difference in expression) the number of genes or proteins that were expressed numerically higher or lower in the high FE compared to the low FE phenotype for key terms or functions, e.g., ribosomal, mitochondrial ribosomal, tRNA, RNA binding motif, RNA polymerase, small nuclear ribonucleoprotein, and protein tyrosine phosphatase, were determined. Bionomial distribution analysis (exact) was then conducted on these datasets to determine significance between numerically up or down expression.

Results: Processes associated with mitochondrial proteome expression (e.g., mitochondrial ribosomal proteins, mitochondrial transcription, mitochondrial tRNA, and translation) were enriched in breast muscle from the high FE compared to the low FE pedigree male broiler phenotype. Furthermore, the high FE phenotype exhibited enrichment of ribosome assembly (e.g., RNA polymerase, mitochondrial and cytosolic ribosomes, small, and heterogeneous nuclear ribonucleoproteins), as well as nuclear transport and protein translation processes compared to the low FE phenotype. Quality control processes (proteosomes and autophagy) were also enriched in the high FE phenotype. In contrast, the low FE phenotype exhibited enrichment of cytoskeletal proteins, protein tyrosine phosphatases, and tyrosine kinases compared to the high FE phenotype. These results suggest that processes of mitochondrial and cytosolic ribosomal construction, activity, and protein translation would be enhanced in high FE breast muscle, and that phosphorylation of tyrosine moieties of proteins could be prolonged in the high compared to low FE phenotype. The results indicate the presence of a proteogenomic architecture that could enhance ribosome construction, protein translation, and quality control processes and contribute to the phenotypic expression of feed efficiency in this PedM broiler model.

A transcriptome multi-tissue analysis identifies biological pathways and genes associated with variations in feed efficiency of growing pigs

BACKGROUND

Animal's efficiency in converting feed into lean gain is a critical issue for the profitability of meat industries. This study aimed to describe shared and specific molecular responses in different tissues of pigs divergently selected over eight generations for residual feed intake (RFI).

RESULTS

Pigs from the low RFI line had an improved gain-to-feed ratio during the test period and displayed higher leanness but similar adiposity when compared with pigs from the high RFI line at 132 days of age. Transcriptomics data were generated from longissimus muscle, liver and two adipose tissues using a porcine microarray and analyzed for the line effect (n = 24 pigs per line). The most apparent effect of the line was seen in muscle, whereas subcutaneous adipose tissue was the less affected tissue. Molecular data were analyzed by bioinformatics and subjected to multidimensional statistics to identify common biological processes across tissues and key genes participating to differences in the genetics of feed efficiency. Immune response, response to oxidative stress and protein metabolism were the main biological pathways shared by the four tissues that distinguished pigs from the low or high RFI lines. Many immune genes were under-expressed in the four tissues of the most efficient pigs. The main genes contributing to difference between pigs from the low vs high RFI lines were CD40, CTSC and NTN1. Different genes associated with energy use were modulated in a tissue-specific manner between the two lines. The gene expression program related to glycogen utilization was specifically up-regulated in muscle of pigs from the low RFI line (more efficient). Genes involved in fatty acid oxidation were down-regulated in muscle but were promoted in adipose tissues of the same pigs when compared with pigs from the high RFI line (less efficient). This underlined opposite line-associated strategies for energy use in skeletal muscle and adipose tissue. Genes related to cholesterol synthesis and efflux in liver and perirenal fat were also differentially regulated in pigs from the low vs high RFI lines.

CONCLUSIONS

Non-productive functions such as immunity, defense against pathogens and oxidative stress contribute likely to inter-individual variations in feed efficiency.

Progesterone signalling in broiler skeletal muscle is associated with divergent feed efficiency

Background

We contrast the pectoralis muscle transcriptomes of broilers selected from within a single genetic line expressing divergent feed efficiency (FE) in an effort to improve our understanding of the mechanistic basis of FE.

Results

Application of a virtual muscle model to gene expression data pointed to a coordinated reduction in slow twitch muscle isoforms of the contractile apparatus (MYH15, TPM3, MYOZ2, TNNI1, MYL2, MYOM3, CSRP3, TNNT2), consistent with diminishment in associated slow machinery (myoglobin and phospholamban) in the high FE animals. These data are in line with the repeated transition from red slow to white fast muscle fibres observed in agricultural species selected on mass and FE. Surprisingly, we found that the expression of 699 genes encoding the broiler mitoproteome is modestly–but significantly–biased towards the high FE group, suggesting a slightly elevated mitochondrial content. This is contrary to expectation based on the slow muscle isoform data and theoretical physiological capacity arguments. Reassuringly, the extreme 40 most DE genes can successfully cluster the 12 individuals into the appropriate FE treatment group. Functional groups contained in this DE gene list include metabolic proteins (including opposing patterns of CA3 and CA4), mitochondrial proteins (CKMT1A), oxidative status (SEPP1, HIG2A) and cholesterol homeostasis (APOA1, INSIG1). We applied a differential network method (Regulatory Impact Factors) whose aim is to use patterns of differential co-expression to detect regulatory molecules transcriptionally rewired between the groups. This analysis clearly points to alterations in progesterone signalling (via the receptor PGR) as the major driver. We show the progesterone receptor localises to the mitochondria in a quail muscle cell line.

Conclusions

Progesterone is sometimes used in the cattle industry in exogenous hormone mixes that lead to a ~20% increase in FE. Because the progesterone receptor can localise to avian mitochondria, our data continue to point to muscle mitochondrial metabolism as an important component of the phenotypic expression of variation in broiler FE.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0396-2) contains supplementary material, which is available to authorized users.

RNA sequencing for global gene expression associated with muscle growth in a single male modern broiler line compared to a foundational Barred Plymouth Rock chicken line

Background

Modern broiler chickens exhibit very rapid growth and high feed efficiency compared to unselected chicken breeds. The improved production efficiency in modern broiler chickens was achieved by the intensive genetic selection for meat production. This study was designed to investigate the genetic alterations accumulated in modern broiler breeder lines during selective breeding conducted over several decades.

Methods

To identify genes important in determining muscle growth and feed efficiency in broilers, RNA sequencing (RNAseq) was conducted with breast muscle in modern pedigree male (PeM) broilers (n = 6 per group), and with an unselected foundation broiler line (Barred Plymouth Rock; BPR). The RNAseq analysis was carried out using Ilumina Hiseq (2 x 100 bp paired end read) and raw reads were assembled with the galgal4 reference chicken genome. With normalized RPM values, genes showing >10 average read counts were chosen and genes showing <0.05 p-value and >1.3 fold change were considered as differentially expressed (DE) between PeM and BPR. DE genes were subjected to Ingenuity Pathway Analysis (IPA) for bioinformatic functional interpretation.

Results

The results indicate that 2,464 DE genes were identified in the comparison between PeM and BPR. Interestingly, the expression of genes encoding mitochondrial proteins in chicken are significantly biased towards the BPR group, suggesting a lowered mitochondrial content in PeM chicken muscles compared to BPR chicken. This result is inconsistent with more slow muscle fibers bearing a lower mitochondrial content in the PeM. The molecular, cellular and physiological functions of DE genes in the comparison between PeM and BPR include organismal injury, carbohydrate metabolism, cell growth/proliferation, and skeletal muscle system development, indicating that cellular mechanisms in modern broiler lines are tightly associated with rapid growth and differential muscle fiber contents compared to the unselected BPR line. Particularly, PDGF (platelet derived growth factor) signaling and NFE2L2 (nuclear factor, erythroid 2-like 2; also known as NRF2) mediated oxidative stress response pathways appear to be activated in modern broiler compared to the foundational BPR line. Upstream and network analyses revealed that the MSTN (myostatin) –FST (follistatin) interactions and inhibition of AR (androgen receptor) were predicted to be effective regulatory factors for DE genes in modern broiler line. PRKAG3 (protein kinase, AMP-activated, gamma 3 non-catalytic subunit) and LIPE (lipase E) are predicted as core regulatory factors for myogenic development, nutrient and lipid metabolism.

Conclusion

The highly upregulated genes in PeM may represent phenotypes of subclinical myopathy commonly observed in the commercial broiler breast tissue, that can lead to muscle hardening, named as woody breast. By investigating global gene expression in a highly selected pedigree broiler line and a foundational breed (Barred Plymouth Rock), the results provide insight into cellular mechanisms that regulate muscle growth, fiber composition and feed efficiency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3471-y) contains supplementary material, which is available to authorized users.

Mitochondrial metabolism: a driver of energy utilisation and product quality?

High feed efficiency is a very desirable production trait as it positively influences resource utilisation, profitability and environmental considerations, albeit at the possible expense of product quality. The modern broiler is arguably the most illustrative model species as it has been transformed over the past half century into an elite feed converter. Some producers are currently reporting that 42-day-old birds gain 1 kg of wet weight for every 1.35 kg of dry weight consumed. Its large breast muscle is exclusively composed of large, low mitochondrial-content Type IIB fibres, which may contribute to low maintenance costs and high efficiency. In an effort to gain a better understanding of individual variation in chicken feed efficiency, our group has been exploring the biology of the mitochondrion at multiple levels of organisation. The mitochondrion is the organelle where much biochemical energy transformation occurs in the cell. Using Cobb-Vantress industrial birds as our primary experimental resource, we have explored the tissue content, structure and function of the mitochondrion and its relationship to growth, development, efficiency and genetic background. While much remains to be understood, recent highlights include (1) variation in muscle mitochondrial content that is associated with performance phenotypes, (2) altered muscle mitochondrial gene and protein expression in birds differing in feed efficiency, (3) variation in isolated mitochondrial function in birds differing in feed efficiency and (4) evidence for an unexpected role for the mitochondrially localised progesterone receptor in altering bird muscle metabolism. Mitochondrial function is largely conserved across the vertebrates, so the same metabolic principles appear to apply to the major production species, whether monogastric or ruminant. A speculative role for the mitochondria in aspects of meat quality and in influencing postmortem anaerobic metabolism will conclude the manuscript.