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

OMICs approaches and technologies for understanding low-high feed efficiency traits in chicken: implication to breeding

In poultry production, there has been a trend of continuous increase in cost of feed ingredients which represents the major proportion of the production costs. Feed costs can be reduced by improving feed efficiency traits which increase the possibility of using various indigestible feed sources and decrease the environmental impact of the enhanced poultry production. Therefore, feed efficiency has been used as one of the most important economic traits of selection in the breeding program of chickens. Recently, many OMICs experimental studies have been designed to characterize biological differences between the high and low feed efficiency chicken phenotypes. Biological complexity cannot be fully captured by main individual OMICs such as genomics, transcriptomics, proteomics and metabolomics. Therefore, researchers have combined multiple assays from the same set of samples to create multi-OMICs datasets. OMICs findings are crucial in improving existing approaches to poultry breeding. The current review aimed to highlight the components of feed efficiency and general OMICs approaches and technologies. Besides, individual and multi-OMICs based understanding of chicken feed efficiency traits and the application of the acquired knowledge in the chicken breeding program were addressed.

Multi-breed host rumen epithelium transcriptome and microbiome associations and their relationship with beef cattle feed efficiency

Understanding host-microbial interactions in the rumen and its influence on desirable production traits may lead to potential microbiota manipulation or genetic selection for improved cattle feed efficiency. This study investigated the host transcriptome and its correlation with the rumen archaea and bacteria differential abundance of two pure beef cattle breeds (Angus and Charolais) and one composite beef hybrid (Kinsella) divergent for residual feed intake (RFI; low-RFI vs. high-RFI). Using RNA-Sequencing of rumen tissue and 16S rRNA gene amplicon sequencing, differentially expressed genes (FDR ≤ 0.05, |log2(Fold-change) >|2) and differentially abundant (p-value < 0.05) archaea and bacteria amplicon sequence variants (ASV) were determined. Significant correlations between gene expression and ASVs (p-value < 0.05) were determine using Spearman correlation. Interesting associations with muscle contraction and the modulation of the immune system were observed for the genes correlated with bacterial ASVs. Potential functional candidate genes for feed efficiency status were identified for Angus (CCL17, CCR3, and CXCL10), Charolais (KCNK9, GGT1 and IL6), and Kinsella breed (ESR2). The results obtained here provide more insights regarding the applicability of target host and rumen microbial traits for the selection and breeding of more feed efficient beef cattle.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

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.

Genome-Wide Association Study Identifies 12 Loci Associated with Body Weight at Age 8 Weeks in Korean Native Chickens

Meat from Korean native chickens (KNCs) has high consumer demand; however, slow growth performance and high variation in body weight (BW) of KNCs remain an issue. Genome-wide association study (GWAS) is a powerful method to identify quantitative trait-associated genomic loci. A GWAS, based on a large-scale KNC population, is needed to identify underlying genetic mechanisms related to its growth traits. To identify BW-associated genomic regions, we performed a GWAS using the chicken 60K single nucleotide polymorphism (SNP) panel for 1328 KNCs. BW was measured at 8 weeks of age, from 2018 to 2020. Twelve SNPs were associated with BW at the suggestive significance level (p < 2.95 × 10⁻⁵) and located near or within 11 candidate genes, including WDR37, KCNIP4, SLIT2, PPARGC1A, MYOCD and ADGRA3. Gene set enrichment analysis based on the GWAS results at p < 0.05 (1680 SNPs) showed that 32 Gene Ontology terms and two Kyoto Encyclopedia of Genes and Genomes pathways, including regulation of transcription, motor activity, the mitogen-activated protein kinase signaling pathway, and tight junction, were significantly enriched (p < 0.05) for BW-associated genes. These pathways are involved in cell growth and development, related to BW gain. The identified SNPs are potential biomarkers in KNC breeding.

A 31-bp indel in the 5' UTR region of GNB1L is significantly associated with chicken body weight and carcass traits

BACKGROUND

G-protein subunit beta 1 like (GNB1L) encodes a G-protein beta-subunit-like polypeptide. Chicken GNB1L is upregulated in the breast muscle of high feed efficiency chickens, and its expression is 1.52-fold that in low feed efficiency chickens. However, no report has described the effects of GNB1L indels on the chicken carcass and growth traits.

RESULTS

This study identified a 31-bp indel in the 5' untranslated region (UTR) of GNB1L and elucidated the effect of this gene mutation on the carcass and growth traits in chickens. The 31-bp indel showed a highly significant association with the body weight at 8 different stages and was significantly correlated with daily gains at 0 to 4 weeks and 4 to 8 weeks. Similarly, the mutation was significantly associated with small intestine length, breast width, breast depth and breast muscle weight. Moreover, DD and ID were superior genotypes for chicken growth and carcass traits.

CONCLUSIONS

These results show that the 31-bp indel of GNB1L significantly affects chicken body weight and carcass traits and can serve as a candidate molecular marker for chicken genetics and breeding programs.

Genes associated with body weight gain and feed intake identified by meta-analysis of the mesenteric fat from crossbred beef steers

Mesenteric fat is a visceral fat depot that increases with cattle maturity and can be influenced by diet. There may be a relationship between the accumulation of mesenteric fat and feed efficiency in beef cattle. The purpose of this study was to identify genes that may be differentially expressed in steers with high and low BW gain and feed intake. RNA-Seq was used to evaluate the transcript abundance of genes in the mesenteric fat from a total of 78 steers collected over 5 different cohorts. A meta-analysis was used to identify genes involved with gain, feed intake or the interaction of both phenotypes. The interaction analysis identified 11 genes as differentially expressed. For the main effect of gain, a total of 87 differentially expressed genes (DEG) were identified (P_ADJ<0.05), and 24 were identified in the analysis for feed intake. Genes identified for gain were involved in functions and pathways including lipid metabolism, stress response/protein folding, cell proliferation/growth, axon guidance and inflammation. The genes for feed intake did not cluster into pathways, but some of the DEG for intake had functions related to inflammation, immunity, and/or signal transduction (JCHAIN, RIPK1, LY86, SPP1, LYZ, CD5, CD53, SRPX, and NF2). At P_ADJ<0.1, only 4 genes (OLFML3, LOC100300716, MRPL15, and PUS10) were identified as differentially expressed in two or more cohorts, highlighting the importance of evaluating the transcriptome of more than one group of animals and incorporating a meta-analysis. This meta-analysis has produced many mesenteric fat DEG that may be contributing to gain and feed intake in cattle.

Distinct cell proliferation, myogenic differentiation, and gene expression in skeletal muscle myoblasts of layer and broiler chickens

Myoblasts play a central role during skeletal muscle formation and growth. Precise understanding of myoblast properties is thus indispensable for meat production. Herein, we report the cellular characteristics and gene expression profiles of primary-cultured myoblasts of layer and broiler chickens. Broiler myoblasts actively proliferated and promptly differentiated into myotubes compared to layer myoblasts, which corresponds well with the muscle phenotype of broilers. Transcriptomes of layer and broiler myoblasts during differentiation were quantified by RNA sequencing. Ontology analyses of the differentially expressed genes (DEGs) provided a series of extracellular proteins as putative markers for characterization of chicken myogenic cells. Another ontology analyses demonstrated that broiler myogenic cells are rich in cell cycle factors and muscle components. Independent of these semantic studies, principal component analysis (PCA) statistically defined two gene sets: one governing myogenic differentiation and the other segregating layers and broilers. Thirteen candidate genes were identified with a combined study of the DEGs and PCA that potentially contribute to proliferation or differentiation of chicken myoblasts. We experimentally proved that one of the candidates, enkephalin, an opioid peptide, suppresses myoblast growth. Our results present a new perspective that the opioids present in feeds may influence muscle development of domestic animals.

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.

Comparative transcriptome analysis reveals regulators mediating breast muscle growth and development in three chicken breeds

Objective: The goal of this study was to investigate the mechanisms of muscle growth and development of three chicken breeds. Participants: Eighteen chickens, including three different breeds with different growth speeds (White Broiler, Daheng, and Commercial Layers of Roman), were used. Methods: Total RNA from breast muscle of these chickens was subjected to a gene expression microarray. Differentially expressed genes (DEGs) were screened and functional enrichment analysis was performed using DAVID. Seven DEGs were confirmed by quantitative reverse transcription PCR. Results: Overall, 8,398 DEGs were found among the different lines. The DEGs between each two lines that were unique for a developmental stage were greater than those that were common during all stages. Functional analysis revealed that DEGs across the entire developmental process were primarily involved in positive cell proliferation, growth, cell differentiation, and developmental processes. Genes involved in muscle regulation, muscle construction, and muscle cell differentiation were upregulated in the faster-growing breed compared to the slower-growing breed. DEGs including myosin heavy chain 15 (MYH15), myozenin 2 (MYOZ2), myosin-binding protein C (MYBPC3), insulin-like growth factor 2 (IGF2), apoptosis regulator (BCL-2), AP-1 transcription factor subunit (JUN), and AP-1 transcription factor subunit (FOS) directly regulated muscle growth or were in the center of the protein-protein interaction network. Pathways, including the extracellular matrix (ECM)-receptor interaction, mitogen-activated protein kinase (MAPK) signaling pathway, and focal adhesion, were the most enriched DEGs between lines or within lines under different developmental stages. Conclusions: Genes involved in muscle construction and cell differentiation were differentially expressed among the three breeds.

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.

Transcriptional shifts account for divergent resource allocation in feed efficient broiler chickens

Considerable variation in feed efficiency (FE) has been observed in indigenous and selected meat-type chicken populations. Although this variation could be partially linked to extrinsic factors like diet, housing environment and microbiota, it further illustrates the existence of strong molecular mechanisms enabling the differential allocation of resources for various physiological processes. To further deepen the molecular basis of individual allocation capacity in male and female broilers, an RNA-seq experiment was conducted which based on a phenotyped chicken population divergent in FE. Transcriptional differences linked to FE were pronounced in intestinal and muscular tissue sites of male animals. Specifically, signalling pathways of farnesoid X receptor (FXR) and retinoid X receptor (RXR) might contribute to mediate individual FE. The transcriptional profiles suggested ACSBG2 (muscular lipid utilisation), ASBT (intestinal bile salt transport), CLEC2B (natural killer cell activation), HMGCS2 (jejunal, duodenal and muscular ketogenesis), and SCARB1 (jejunal lipid uptake) as potential mediators driving FE. Results indicate that improvements in FE exploit shifts in resource allocation which might occur at the expense of general immune responsiveness in high efficient male chickens. Consequently, to further improve FE traits and to explore causative molecular patterns, effects originating from sex-dimorphism in chickens need to be taken into consideration.

Pathway analysis identifies altered mitochondrial metabolism, neurotransmission, structural pathways and complement cascade in retina/RPE/ choroid in chick model of form-deprivation myopia

Purpose

RNA sequencing analysis has demonstrated bidirectional changes in metabolism, structural and immune pathways during early induction of defocus induced myopia. Thus, the aim of this study was to investigate whether similar gene pathways are also related to the more excessive axial growth, ultrastructural and elemental microanalytic changes seen during the induction and recovery from form-deprivation myopia (FDM) in chicks and predicted by the RIDE model of myopia.

Methods

Archived genomic transcriptome data from the first three days of induction of monocularly occluded form deprived myopia (FDMI) in chicks was obtained from the GEO database (accession # GSE6543) while data from chicks monocularly occluded for 10 days and then given up to 24 h of normal visual recovery (FDMR) were collected. Gene set enrichment analysis (GSEA) software was used to determine enriched pathways during the induction (FDMI) and recovery (FDMR) from FD. Curated gene-sets were obtained from open access sources.

Results

Clusters of significant changes in mitochondrial energy metabolism, neurotransmission, ion channel transport, G protein coupled receptor signalling, complement cascades and neuron structure and growth were identified during the 10 days of induction of profound myopia and were found to correlate well with change in axial dimensions. Bile acid and bile salt metabolism pathways (cholesterol/lipid metabolism and sodium channel activation) were significantly upregulated during the first 24 h of recovery from 10 days of FDM.

Conclusions

The gene pathways altered during induction of FDM are similar to those reported in defocus induced myopia and are established indicators of oxidative stress, osmoregulatory and associated structural changes. These findings are also consistent with the choroidal thinning, axial elongation and hyperosmotic ion distribution patterns across the retina and choroid previously reported in FDM and predicted by RIDE.

Integrative approach using liver and duodenum RNA-Seq data identifies candidate genes and pathways associated with feed efficiency in pigs

This study aims identifying candidate genes and pathways associated with feed efficiency (FE) in pigs. Liver and duodenum transcriptomes of 37 gilts showing high and low residual feed intake (RFI) were analysed by RNA-Seq. Gene expression data was explored through differential expression (DE) and weighted gene co-expression network analyses. DE analysis revealed 55 and 112 differentially regulated genes in liver and duodenum tissues, respectively. Clustering genes according to their connectivity resulted in 23 (liver) and 25 (duodenum) modules of genes with a co-expression pattern. Four modules, one in liver (with 444 co-expressed genes) and three in duodenum (gathering 37, 126 and 41 co-expressed genes), were significantly associated with FE indicators. Intra-module analyses revealed tissue-specific candidate genes; 12 of these genes were also identified as DE between individuals with high and low RFI. Pathways enriched by the list of genes showing DE and/or belonging to FE co-expressed modules included response to oxidative stress, inflammation, immune response, lipid metabolism and thermoregulation. Low overlapping between genes identified in duodenum and liver tissues was observed but heat shock proteins were associated to FE in both tissues. Our results suggest tissue-specific rather than common transcriptome regulatory processes associated with FE in pigs.

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.

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

As feed represents 60 to 70% of the cost of raising an animal to market weight, feed efficiency (the amount of dry weight intake to amount of wet weight gain) remains an important genetic trait in animal agriculture. To gain greater understanding of cellular mechanisms of feed efficiency (FE), shotgun proteomics was conducted using in-gel trypsin digestion and tandem mass spectrometry on breast muscle samples obtained from pedigree male (PedM) broilers exhibiting high feed efficiency (FE) or low FE phenotypes (n = 4 per group). The high FE group had greater body weight gain (P = 0.004) but consumed the same amount of feed (P = 0.30) from 6 to 7 wk resulting in higher FE (P < 0.001). Over 1800 proteins were identified, of which 152 were different (P < 0.05) by at least 1.3 fold and ≤ 15 fold between the high and low FE phenotypes. Data were analyzed for a modified differential expression (DE) metric (Phenotypic Impact Factors or PIF) and interpretation of protein expression data facilitated using the Ingenuity Pathway Analysis (IPA) program. In the entire data set, 228 mitochondrial proteins were identified whose collective expression indicates a higher mitochondrial expression in the high FE phenotype (binomial probability P < 0.00001). Within the top up and down 5% PIF molecules in the dataset, there were 15 mitoproteome proteins up-regulated and only 5 down-regulated in the high FE phenotype. Pathway enrichment analysis also identified mitochondrial dysfunction and oxidative phosphorylation as the number 1 and 5 differentially expressed canonical pathways (up-regulated in high FE) in the proteomic dataset. Upstream analysis (based on DE of downstream molecules) predicted that insulin receptor, insulin like growth receptor 1, nuclear factor, erythroid 2-like 2, AMP activated protein kinase (α subunit), progesterone and triiodothyronine would be activated in the high FE phenotype whereas rapamycin independent companion of target of rapamycin, mitogen activated protein kinase 4, and serum response factor would be inhibited in the high FE phenotype. The results provide additional insight into the fundamental molecular landscape of feed efficiency in breast muscle of broilers as well as further support for a role of mitochondria in the phenotypic expression of FE.

Funding provided by USDA-NIFA (#2013–01953), Arkansas Biosciences Institute (Little Rock, AR), McMaster Fellowship (AUS to WB) and the Agricultural Experiment Station (Univ. of Arkansas, Fayetteville).

RNA-Seq Analysis of Abdominal Fat Reveals Differences between Modern Commercial Broiler Chickens with High and Low Feed Efficiencies

For economic and environmental reasons, chickens with superior feed efficiency (FE) are preferred in the broiler chicken industry. High FE (HFE) chickens typically have reduced abdominal fat, the major adipose tissue in chickens. In addition to its function of energy storage, adipose tissue is a metabolically active organ that also possesses endocrine and immune regulatory functions. It plays a central role in maintaining energy homeostasis. Comprehensive understanding of the gene expression in the adipose tissue and the biological basis of FE are of significance to optimize selection and breeding strategies. Through gene expression profiling of abdominal fat from high and low FE (LFE) commercial broiler chickens, the present study aimed to characterize the differences of gene expression between HFE and LFE chickens. mRNA-seq analysis was carried out on the total RNA of abdominal fat from 10 HFE and 12 LFE commercial broiler chickens, and 1.48 billion of 75-base sequence reads were generated in total. On average, 11,565 genes were expressed (>5 reads/gene/sample) in the abdominal fat tissue, of which 286 genes were differentially expressed (DE) at q (False Discover Rate) < 0.05 and fold change > 1.3 between HFE and LFE chickens. Expression levels from RNA-seq were confirmed with the NanoString nCounter analysis system. Functional analysis showed that the DE genes were significantly (p < 0.01) enriched in lipid metabolism, coagulation, and immune regulation pathways. Specifically, the LFE chickens had higher expression of lipid synthesis genes and lower expression of triglyceride hydrolysis and cholesterol transport genes. In conclusion, our study reveals the overall differences of gene expression in the abdominal fat from HFE and LFE chickens, and the results suggest that the divergent expression of lipid metabolism genes represents the major differences.

Messenger RNA sequencing and pathway analysis provide novel insights into the biological basis of chickens' feed efficiency

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

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

A replication study for genome-wide gene expression levels in two layer lines elucidates differentially expressed genes of pathways involved in bone remodeling and immune responsiveness

The current replication study confirmed significant differences in gene expression profiles of the cerebrum among the two commercial layer lines Lohmann Selected Leghorn (LSL) and Lohmann Brown (LB). Microarray analyses were performed for 30 LSL and another 30 LB laying hens kept in the small group housing system Eurovent German. A total of 14,103 microarray probe sets using customized Affymetrix ChiGene-1_0-st Arrays with 20,399 probe sets were differentially expressed among the two layer lines LSL and LB (FDR adjusted P-value <0.05). An at least 2-fold change in expression levels could be observed for 388 of these probe sets. In LSL, 214 of the 388 probe sets were down- and 174 were up-regulated and vice versa for the LB layer line. Among the 174 up-regulated probe sets in LSL, we identified 51 significantly enriched Gene ontology (GO) terms of the biological process category. A total of 63 enriched GO-terms could be identified for the 214 down-regulated probe sets of the layer line LSL. We identified nine genes significantly differentially expressed between the two layer lines in both microarray experiments. These genes play a crucial role in protection of neuronal cells from oxidative stress, bone mineral density and immune response among the two layer lines LSL and LB. Thus, the different regulation of these genes may significantly contribute to phenotypic trait differences among these layer lines. In conclusion, these novel findings provide a basis for further research to improve animal welfare in laying hens and these layer lines may be of general interest as an animal model.

Microarray resources for genetic and genomic studies in chicken: a review

Advent of microarray technologies revolutionized the nature and scope of genetic and genomic research in human and other species by allowing massively parallel analysis of thousands of genomic sites. They have been used for diverse purposes such as for transcriptome analysis, CNV detection, SNP and CNV genotyping, studying DNA-protein interaction, and detection of genome methylation. Microarrays have also made invaluable contributions to research in chicken which is an important model organism for studying embryology, immunology, oncology, virology, evolution, genetics, and genomics and also for other avian species. Despite their huge contributions in life science research, the future of microarrays is now being questioned with the advent of massively parallel next generation sequencing (NGS) technologies, which promise to overcome some of the limitations of microarray platforms. In this article we review the various microarray resources developed for chicken and their past and potential future applications. We also discuss about the future of microarrays in the NGS era particularly in the context of livestock genetics. We argue that even though NGS promises some major advantages-in particular, offers the opportunity to discover novel elements in the genome-microarrays will continue to be major tools for research and practice in the field of livestock genetics/genomics due to their affordability, high throughput nature, mature established technologies and ease of application. Moreover, with advent of new microarray technologies like capture arrays, the NGS and microarrays are expected to complement each other in future research in life science.

Cell Biology Symposium: feed efficiency: mitochondrial function to global gene expression

Understanding the cellular basis of feed efficiency (FE) is instrumental to helping poultry and livestock industries continue to provide high-quality protein for an increasingly crowded world. To understand relationships of FE and gene expression, global RNA transcription was investigated in breast muscle obtained from a male broiler line fed the same diet and individually phenotyped for FE. In these studies, RNA samples obtained from broilers that exhibited either high FE (0.65 ± 0.01) or low FE (0.46 ± 0.01) were analyzed with an Agilent 44K chicken oligoarray. A 1.3-fold cutoff in expression (30% difference between groups) resulted in 782 genes that were differentially expressed (P < 0.05) in muscle between the high- and low-FE phenotypes. Ingenuity Pathway Analysis, an online software program, was used to identify genes, gene networks, and pathways associated with the phenotypic expression of FE. The results indicate that the high-FE phenotype exhibited increased expression of genes associated with 1) signal transduction pathways, 2) anabolic activities, and 3) energy-sensing and energy coordination activities, all of which would likely be favorable to cell growth and development. In contrast, the low-FE broiler phenotype exhibited upregulation of genes 1) associated with actin-myosin filaments, cytoskeletal architecture, and muscle fibers and 2) stress-related or stress-responsive genes. Because the low-FE broiler phenotype exhibits greater oxidative stress, it would appear that the low-FE phenotype is the product of inherent gene expression that is modulated by oxidative stress. The results of these studies begin to provide a comprehensive picture of gene expression in muscle, a major organ of energy demand in an animal, associated with phenotypic expression of FE.