Research Note: Chicken breast muscle satellite cell function: effect of expression of CNN1 and PHRF1.

The Wooden Breast myopathy results in the necrosis and fibrosis of breast muscle fibers in fast-growing heavy weight meat-type broiler chickens. Myogenic satellite cells are required to repair and regenerate the damaged muscle fibers. Using Genome Wide Association, candidate genes affected with Wooden Breast have been previously reported. The effect of these genes on satellite cell proliferation, differentiation, and the synthesis of lipids by satellite cells is unknown. Satellite cells isolated from the pectoralis major muscle from commercial Ross 708 broilers and a Randombred chicken (RBch) line were used. Expression of calponin 1 (CNN1) and PHD and ring fingers domains 1 (PHRF1) were knocked down by silent interfering RNA to determine their effect on satellite cell-mediated proliferation, differentiation, and lipid accumulation. CNN1 and PHRF1 affected satellite cell activity and lipid accumulation in both lines. Proliferation was reduced in the Ross 708 and RBch lines by knocking down the expression of both genes, and differentiation was affected with a line and treatment interaction when gene expression was reduced at the beginning of proliferation. During differentiation lipid accumulation was decreased with knocking down the expression of CNN1 and PHRF1. Both CNN1 and PHRF1 have not been reported previously in skeletal muscle and further research is required to determine their effect on satellite cell-mediated growth and regeneration of the pectoralis major (breast) muscle.

Spatial transcriptomics reveals alterations in perivascular macrophage lipid metabolism in the onset of Wooden Breast myopathy in broiler chickens.

This study aims to use spatial transcriptomics to characterize the cell-type-specific expression profile associated with the microscopic features observed in Wooden Breast myopathy. 1 cm3 muscle sample was dissected from the cranial part of the right pectoralis major muscle from three randomly sampled broiler chickens at 23 days post-hatch and processed with Visium Spatial Gene Expression kits (10X Genomics), followed by high-resolution imaging and sequencing on the Illumina Nextseq 2000 system. WB classification was based on histopathologic features identified. Sequence reads were aligned to the chicken reference genome (Galgal6) and mapped to histological images. Unsupervised K-means clustering and Seurat integrative analysis differentiated histologic features and their specific gene expression pattern, including lipid laden macrophages (LLM), unaffected myofibers, myositis and vasculature. In particular, LLM exhibited reprogramming of lipid metabolism with up-regulated lipid transporters and genes in peroxisome proliferator-activated receptors pathway, possibly through P. Moreover, overexpression of fatty acid binding protein 5 could enhance fatty acid uptake in adjacent veins. In myositis regions, increased expression of cathepsins may play a role in muscle homeostasis and repair by mediating lysosomal activity and apoptosis. A better knowledge of different cell-type interactions at early stages of WB is essential in developing a comprehensive understanding.

Proteomic insight into human directed selection of the domesticated chicken Gallus gallus.

Chicken domestication began at least 3,500 years ago for purposes of divination, cockfighting, and food. Prior to industrial scale chicken production, domestication selected larger birds with increased egg production. In the mid-20th century companies began intensive selection with the broiler (meat) industry focusing on improved feed conversion, rapid growth, and breast muscle yield. Here we present proteomic analysis comparing the modern broiler line, Ross 708, with the UIUC legacy line which is not selected for growth traits. Breast muscle proteome analysis identifies cellular processes that have responded to human directed artificial selection. Mass spectrometry was used to identify protein level differences in the breast muscle of 6-day old chicks from Modern and Legacy lines. Our results indicate elevated levels of stress proteins, ribosomal proteins and proteins that participate in the innate immune pathway in the Modern chickens. Furthermore, the comparative analyses indicated expression differences for proteins involved in multiple biochemical pathways. In particular, the Modern line had elevated levels of proteins affecting the pentose phosphate pathway, TCA cycle and fatty acid oxidation while proteins involved in the first phase of glycolysis were reduced compared to the Legacy line. These analyses provide hypotheses linking the morphometric changes driven by human directed selection to biochemical pathways. These results also have implications for the poultry industry, specifically Wooden Breast disease which is linked to rapid breast muscle growth.

Integrative transcriptomic and metabolomic analysis reveals alterations in energy metabolism and mitochondrial functionality in broiler chickens with wooden breast.

This integrative study of transcriptomics and metabolomics aimed to improve our understanding of Wooden Breast myopathy (WB). Breast muscle samples from 8 WB affected and 8 unaffected male broiler chickens of 47 days of age were harvested for metabolite profiling. Among these 16 samples, 5 affected and 6 unaffected also underwent gene expression profiling. The Joint Pathway Analysis was applied on 119 metabolites and 3444 genes exhibiting differential abundance or expression between WB affected and unaffected chickens. Mitochondrial dysfunctions in WB was suggested by higher levels of monoacylglycerols and down-regulated genes involved in lipid production, fatty acid beta oxidation, and oxidative phosphorylation. Lower levels of carnosine and anserine, along with down-regulated carnosine synthase 1 suggested decreased carnosine synthesis and hence impaired antioxidant capacity in WB. Additionally, Weighted Gene Co-expression Network Analysis results indicated that abundance of inosine monophosphate, significantly lower in WB muscle, was correlated with mRNA expression levels of numerous genes related to focal adhesion, extracellular matrix and intercellular signaling, implying its function in connecting and possibly regulating multiple key biological pathways. Overall, this study showed not only the consistency between transcript and metabolite profiles, but also the potential in gaining further insights from analyzing multi-omics data.

Identification of circulating metabolites associated with wooden breast and white striping.

Current diagnostic methods for wooden breast and white striping, common breast muscle myopathies of modern commercial broiler chickens, rely on subjective examinations of the pectoralis major muscle, time-consuming microscopy, or expensive imaging technologies. Further research on these disorders would benefit from more quantitative and objective measures of disease severity that can be used in live birds. To this end, we utilized untargeted metabolomics alongside two statistical approaches to evaluate plasma metabolites associated with wooden breast and white striping in 250 male commercial broiler chickens. First, mixed linear modeling was employed to identify metabolites with a significant association with these muscle disorders and found 98 metabolites associated with wooden breast and 44 metabolites associated with white striping (q-value < 0.05). Second, a support vector machine was constructed using stepwise feature selection to determine the smallest subset of metabolites with the highest categorization accuracy for wooden breast. The final support vector machine achieved 94% accuracy using only 6 metabolites. The metabolite 3-methylhistidine, which is often used as an index of myofibrillar breakdown in skeletal muscle, was the top metabolite for both wooden breast and white striping in our mixed linear model and was also the metabolite with highest marginal prediction accuracy (82%) for wooden breast in our support vector machine. Overall, this study identified a candidate set of metabolites for an objective measure of wooden breast or white striping severity in live birds and expanded our understanding of these muscle disorders.

3'UTR-Seq analysis of chicken abdominal adipose tissue reveals widespread intron retention in 3'UTR and provides insight into molecular basis of feed efficiency.

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3' UTR-sequencing (3' UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3' UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3' UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2).

Effect of expression of PPARG, DNM2L, RRAD, and LINGO1 on broiler chicken breast muscle satellite cell function.

Disorders affecting the breast muscle of modern commercial broiler chickens have increased in recent years. Wooden Breast (WB) myopathy is characterized by a palpably hard breast muscle with increased fat deposition. WB is a metabolic disorder with lipid accumulation considered to be a primary causal factor. The adult myoblasts, satellite cells, are a partially differentiated stem cell population and primarily function in muscle growth and regeneration. The satellite cells also express adipogenic genes. The objective of this study was to determine the expression of the adipogenic genes PPARG, DNM2L, RRAD, and LINGO1 in commercial Ross 708 (708) and Randombred (RBch) satellite cells. RBch satellite cells are from commercial 1995 broilers before WB and 708 broilers are a modern commercial line. In general, expression of these genes was different between the 708 and RBch satellite cells during proliferation and differentiation. Expression of PPARG and RRAD were both significantly increased during both proliferation and differentiation in the 708 cells (P ≤ 0.05). Knocking down the expression of these genes with small interfering RNAs did not greatly affect either proliferation or differentiation. Lipid accumulation was affected by the knockdown of these genes with significant line effects from 48 h of proliferation through 72 h of differentiation. In general, 708 satellite cells had higher lipid levels. Knockdown treatment effect was specific to each gene. The results demonstrate that lipid biosynthesis has been affected in breast muscle satellite cells which may contribute to the increased lipid deposition in modern day commercial broiler chickens.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species.

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

Genetic basis and identification of candidate genes for wooden breast and white striping in commercial broiler chickens.

Wooden breast (WB) and white striping (WS) are highly prevalent and economically damaging muscle disorders of modern commercial broiler chickens characterized respectively by palpable firmness and fatty white striations running parallel to the muscle fiber. High feed efficiency and rapid growth, especially of the breast muscle, are believed to contribute to development of such muscle defects; however, their etiology remains poorly understood. To gain insight into the genetic basis of these myopathies, a genome-wide association study was conducted using a commercial crossbred broiler population (n = 1193). Heritability was estimated at 0.5 for WB and WS with high genetic correlation between them (0.88). GWAS revealed 28 quantitative trait loci (QTL) on five chromosomes for WB and 6 QTL on one chromosome for WS, with the majority of QTL for both myopathies located in a ~ 8 Mb region of chromosome 5. This region has highly conserved synteny with a portion of human chromosome 11 containing a cluster of imprinted genes associated with growth and metabolic disorders such as type 2 diabetes and Beckwith-Wiedemann syndrome. Candidate genes include potassium voltage-gated channel subfamily Q member 1 (KCNQ1), involved in insulin secretion and cardiac electrical activity, lymphocyte-specific protein 1 (LSP1), involved in inflammation and immune response.

Investigation of allele specific expression in various tissues of broiler chickens using the detection tool VADT.

Differential abundance of allelic transcripts in a diploid organism, commonly referred to as allele specific expression (ASE), is a biologically significant phenomenon and can be examined using single nucleotide polymorphisms (SNPs) from RNA-seq. Quantifying ASE aids in our ability to identify and understand cis-regulatory mechanisms that influence gene expression, and thereby assist in identifying causal mutations. This study examines ASE in breast muscle, abdominal fat, and liver of commercial broiler chickens using variants called from a large sub-set of the samples (n = 68). ASE analysis was performed using a custom software called VCF ASE Detection Tool (VADT), which detects ASE of biallelic SNPs using a binomial test. On average ~ 174,000 SNPs in each tissue passed our filtering criteria and were considered informative, of which ~ 24,000 (~ 14%) showed ASE. Of all ASE SNPs, only 3.7% exhibited ASE in all three tissues, with ~ 83% showing ASE specific to a single tissue. When ASE genes (genes containing ASE SNPs) were compared between tissues, the overlap among all three tissues increased to 20.1%. Our results indicate that ASE genes show tissue-specific enrichment patterns, but all three tissues showed enrichment for pathways involved in translation.

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.

Evidence of vascular endothelial dysfunction in Wooden Breast disorder in chickens: Insights through gene expression analysis, ultra-structural evaluation and supervised machine learning methods.

Several gene expression studies have been previously conducted to characterize molecular basis of Wooden Breast myopathy in commercial broiler chickens. These studies have generally used a limited sample size and relied on a binary disease outcome (unaffected or affected by Wooden Breast), which are appropriate for an initial investigation. However, to identify biomarkers of disease severity and development, it is necessary to use a large number of samples with a varying degree of disease severity. Therefore, in this study, we assayed a relatively large number of samples (n = 96) harvested from the pectoralis major muscle of unaffected (U), partially affected (P) and markedly affected (A) chickens. Gene expression analysis was conducted using the nCounter MAX Analysis System and data were analyzed using four different supervised machine-learning methods, including support vector machines (SVM), random forests (RF), elastic net logistic regression (ENET) and Lasso logistic regression (LASSO). The SVM method achieved the highest prediction accuracy for both three-class (U, P and A) and two-class (U and P+A) classifications with 94% prediction accuracy for two-class classification and 85% for three-class classification. The results also identified biomarkers of Wooden Breast severity and development. Additionally, gene expression analysis and ultrastructural evaluations provided evidence of vascular endothelial cell dysfunction in the early pathogenesis of Wooden Breast.

Blood Gas Disturbances and Disproportionate Body Weight Distribution in Broilers With Wooden Breast.

Wooden breast syndrome is a widespread and economically important myopathy and vasculopathy of fast growing, commercial broiler chickens, primarily affecting birds with high feed efficiency and large breast muscle yield. To investigate potential systemic physiological differences between birds affected and unaffected by wooden breast, a total of 103 market-age Cobb 500 broilers were sampled for 13 blood parameters and the relative weights of the pectoralis major muscle, pectoralis minor muscle, external oblique muscle, wing, heart, lungs, liver, and spleen. Blood analysis was performed on samples taken from the brachial vein of live birds and revealed significant differences in venous blood gases between affected and unaffected chickens. Chickens with wooden breast exhibited significantly higher potassium (K+) and lower partial pressure of oxygen (pO2), oxygen saturation (sO2), and pH. Additionally, affected males had significantly higher partial pressure of carbon dioxide (pCO2) and total carbon dioxide (TCO2) than unaffected males. Wooden breast affected broilers also possessed a significantly heavier pectoralis major muscle and whole feathered wing compared to unaffected broilers. Blood gas disturbances characterized by high pCO2 and low pH are indicative of insufficient respiratory gas exchange, suggesting that wooden breast affected broilers have an elevated metabolic rate that may also be inadequately compensated due to cardiovascular deficiencies such as poor venous return or respiratory insufficiency. Lung tissues from 12 birds with extreme sO2 values were subsequently examined to assess whether lung pathology contributed to the observed blood gas disturbance. Comparison of lung morphology between affected and unaffected birds revealed no apparent differences that could contribute to decreased parabronchial gas exchange. However, an interesting finding was the detection of pulmonary phlebitis in one of the wooden breast-affected samples consistent with vascular changes observed in pectoralis major muscle exhibiting the wooden breast phenotype. Our results suggest that the effects of wooden breast are not limited to the pectoralis major muscle and further indicate the importance of research into metabolic changes associated with the myopathy.

Glucolipotoxicity: A Proposed Etiology for Wooden Breast and Related Myopathies in Commercial Broiler Chickens.

Wooden breast is one of several myopathies of fast-growing commercial broilers that has emerged as a consequence of intensive selection practices in the poultry breeding industry. Despite the substantial economic burden presented to broiler producers worldwide by wooden breast and related muscle disorders such as white striping, the genetic and etiological underpinnings of these diseases are still poorly understood. Here we propose a new hypothesis on the primary causes of wooden breast that implicates dysregulation of lipid and glucose metabolism. Our hypothesis addresses recent findings that have suggested etiologic similarities between wooden breast and type 2 diabetes despite their phenotypic disparities. Unlike in mammals, dysregulation of lipid and glucose metabolism is not accompanied by an increase in plasma glucose levels but generates a unique skeletal muscle phenotype, i.e., wooden breast, in chickens. We hypothesize that these phenotypic disparities result from a major difference in skeletal muscle glucose transport between birds and mammals, and that the wooden breast phenotype most closely resembles complications of diabetes in smooth and cardiac muscle of mammals. Additional basic research on wooden breast and related muscle disorders in commercial broiler chickens is necessary and can be informative for poultry breeding and production as well as for human health and disease. To inform future studies, this paper reviews the current biological knowledge of wooden breast, outlines the major steps in its proposed pathogenesis, and examines how selection for production traits may have contributed to its prevalence.

Dysregulation of lipid metabolism and appearance of slow myofiber-specific isoforms accompany the development of Wooden Breast myopathy in modern broiler chickens.

Previous transcriptomic studies have hypothesized the occurrence of slow myofiber-phenotype, and dysregulation of lipid metabolism as being associated with the development of Wooden Breast (WB), a meat quality defect in commercial broiler chickens. To gain a deep understanding of the manifestation and implication of these two biological processes in health and disease states in chickens, cellular and global expression of specific genes related to the respective processes were examined in pectoralis major muscles of modern fast-growing and unselected slow-growing chickens. Using RNA in situ hybridization, lipoprotein lipase (LPL) was found to be expressed in endothelial cells of capillaries and small-caliber veins in chickens. RNA-seq analysis revealed upregulation of lipid-related genes in WB-affected chickens at week 3 and downregulation at week 7 of age. On the other hand, cellular localization of slow myofiber-type genes revealed their increased expression in mature myofibers of WB-affected chickens. Similarly, global expression of slow myofiber-type genes showed upregulation in affected chickens at both timepoints. To our knowledge, this is the first study to show the expression of LPL from the vascular endothelium in chickens. This study also confirms the existence of slow myofiber-phenotype and provides mechanistic insights into increased lipid uptake and metabolism in WB disease process.

Variant analysis pipeline for accurate detection of genomic variants from transcriptome sequencing data.

The wealth of information deliverable from transcriptome sequencing (RNA-seq) is significant, however current applications for variant detection still remain a challenge due to the complexity of the transcriptome. Given the ability of RNA-seq to reveal active regions of the genome, detection of RNA-seq SNPs can prove valuable in understanding the phenotypic diversity between populations. Thus, we present a novel computational workflow named VAP (Variant Analysis Pipeline) that takes advantage of multiple RNA-seq splice aware aligners to call SNPs in non-human models using RNA-seq data only. We applied VAP to RNA-seq from a highly inbred chicken line and achieved high accuracy when compared with the matching whole genome sequencing (WGS) data. Over 65% of WGS coding variants were identified from RNA-seq. Further, our results discovered SNPs resulting from post transcriptional modifications, such as RNA editing, which may reveal potentially functional variation that would have otherwise been missed in genomic data. Even with the limitation in detecting variants in expressed regions only, our method proves to be a reliable alternative for SNP identification using RNA-seq data. The source code and user manuals are available at https://modupeore.github.io/VAP/.

Increased Expression of Lipid Metabolism Genes in Early Stages of Wooden Breast Links Myopathy of Broilers to Metabolic Syndrome in Humans.

Wooden breast is a muscle disorder affecting modern commercial broiler chickens that causes a palpably firm pectoralis major muscle and severe reduction in meat quality. Most studies have focused on advanced stages of wooden breast apparent at market age, resulting in limited insights into the etiology and early pathogenesis of the myopathy. Therefore, the objective of this study was to identify early molecular signals in the wooden breast transcriptional cascade by performing gene expression analysis on the pectoralis major muscle of two-week-old birds that may later exhibit the wooden breast phenotype by market age at 7 weeks. Biopsy samples of the left pectoralis major muscle were collected from 101 birds at 14 days of age. Birds were subsequently raised to 7 weeks of age to allow sample selection based on the wooden breast phenotype at market age. RNA-sequencing was performed on 5 unaffected and 8 affected female chicken samples, selected based on wooden breast scores (0 to 4) assigned at necropsy where affected birds had scores of 2 or 3 (mildly or moderately affected) while unaffected birds had scores of 0 (no apparent gross lesions). Differential expression analysis identified 60 genes found to be significant at an FDR-adjusted p-value of 0.05. Of these, 26 were previously demonstrated to exhibit altered expression or genetic polymorphisms related to glucose tolerance or diabetes mellitus in mammals. Additionally, 9 genes have functions directly related to lipid metabolism and 11 genes are associated with adiposity traits such as intramuscular fat and body mass index. This study suggests that wooden breast disease is first and foremost a metabolic disorder characterized primarily by ectopic lipid accumulation in the pectoralis major.

RNA-Seq Analysis Reveals Spatial and Sex Differences in Pectoralis Major Muscle of Broiler Chickens Contributing to Difference in Susceptibility to Wooden Breast Disease.

Wooden Breast Disease (WBD) is a novel myopathy affecting the pectoralis major muscle of modern broiler chickens. The etiology of WBD is not currently known, but has been linked to increased feed efficiency, growth rate, and muscle yield in broiler chickens. Differential effect of WBD has been detected between regions of the P. major and between sexes of broilers-male birds and the cranial aspect of the muscle tend to be more severely affected by the disease than females and the caudal aspect. This study aimed to characterize biological differences in the P. major between regions of the muscle and sexes of birds. Samples were taken from the cranial and caudal aspects of P. major muscles of 3-week-old, unaffected male and female birds for RNA sequencing. RNA was extracted and used for preparation of cDNA libraries, which were sequenced by the Delaware Biotechnology Institute (DBI) using HiSeq2500. Sequence reads were aligned to the chicken reference genome with HISAT, and genes were analyzed for differential expression between regions of the breast muscle and sexes of birds using CuffDiff. Functional analysis was performed on differentially expressed genes (DEGs) between sex groups using DAVID and Ingenuity Pathway Analysis (IPA). There were 12 DEGs between cranial and caudal samples, and 260 between male and female birds. Out of the 260 genes differentially expressed between sexes, 189 were upregulated in males. Of this subset, 103 genes (55%) were located on the Z-chromosome. There was increased expression of genes involved in fat metabolism and oxidative stress responses in the cranial region of the P. major muscle, as well as increased expression of fat metabolism, oxidative stress response, antiangiogenesis, and connective tissue proliferation genes in male broilers. These results support the hypothesis that there are biological characteristics in male broilers and the cranial region of the breast muscle that may make them more susceptible to WBD, as well as raising the possibility of a metabolic switch in modern broiler chickens that may be more prominent in males.

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.

RNA-Seq Analyses Identify Additivity as the Predominant Gene Expression Pattern in F1 Chicken Embryonic Brain and Liver.

The superior performance of hybrids to parents, termed heterosis, has been widely utilized in animal and plant breeding programs, but the molecular mechanism underlying heterosis remains an enigma. RNA-Seq provides a novel way to investigate heterosis at the transcriptome-wide level, because gene expression functions as an intermediate phenotype that contributes to observable traits. Here we compared embryonic gene expression between chicken hybrids and their inbred parental lines to identify inheritance patterns of gene expression. Inbred Fayoumi and Leghorn were crossed reciprocally to obtain F1 fertile eggs. RNA-Seq was carried out using 24 brain and liver samples taken from day 12 embryos, and the differentially expressed (DE) genes were identified by pairwise comparison among the hybrids, parental lines, and mid-parent expression values. Our results indicated the expression levels of the majority of the genes in the F1 cross are not significantly different from the mid-parental values, suggesting additivity as the predominant gene expression pattern in the F1. The second and third prevalent gene expression patterns are dominance and over-dominance. Additionally, we found only 7⁻20% of the DE genes exhibit allele-specific expression in the F1, suggesting that trans regulation is the main driver for differential gene expression and thus contributes to heterosis effect in the F1 crosses.