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.

Evidence and role of phlebitis and lipid infiltration in the onset and pathogenesis of Wooden Breast Disease in modern broiler chickens.

Wooden Breast Disease (WBD), a myopathy that frequently affects modern broiler chickens, is a disorder that has been associated with significant economic losses in the poultry industry. To examine tissue changes associated with the onset and early pathogenesis of this disorder, a time-series experiment was conducted using chickens from a high-breast-muscle-yield, purebred commercial broiler line. Birds were raised for up to seven weeks, with a subset of birds sampled weekly. Breast muscle tissues were extracted at necropsy and processed for analysis by light microscopy and transmission electron microscopy. Histologic presentation indicated localized phlebitis with lipogranulomas in Week 1, focal single-myofibril degeneration in Week 2 preceding an inflammatory response that started in Week 3. Lesions in Week 4 were characterized by multifocal to diffuse muscle fibre degeneration, necrosis, interstitial oedema accompanied by increased lipid and inflammatory cell infiltration. Lesions in Weeks 5-7 revealed diffuse muscle degeneration, necrosis, fibrosis and fatty infiltration with lipogranulomas. Ultrastructural examination showed myofibrillar splitting and degeneration, irregular, displaced and degenerated Z-lines, mitochondrial degeneration and interstitial fibrosis with dense regular collagen fibres. This study, therefore, demonstrates that WBD exhibits an earlier onset in modern broilers than when detectable by clinical examination. Further, this study shows that the disease assumes a progressive course with acute vasculitis, lipid deposition and myodegeneration occurring in the earlier stages, followed by a chronic fibrotic phase.

Detection of genomic signatures of recent selection in commercial broiler chickens.

BACKGROUND: Identification of the genomic signatures of recent selection may help uncover causal polymorphisms controlling traits relevant to recent decades of selective breeding in livestock. In this study, we aimed at detecting signatures of recent selection in commercial broiler chickens using genotype information from single nucleotide polymorphisms (SNPs). A total of 565 chickens from five commercial purebred lines, including three broiler sire (male) lines and two broiler dam (female) lines, were genotyped using the 60K SNP Illumina iSelect chicken array. To detect genomic signatures of recent selection, we applied two methods based on population comparison, cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR), and further analyzed the results to find genomic regions under recent selection in multiple purebred lines. RESULTS: A total of 321 candidate selection regions spanning approximately 1.45 % of the chicken genome in each line were detected by consensus of results of both XP-EHH and XP-CLR methods. To minimize false discovery due to genetic drift, only 42 of the candidate selection regions that were shared by 2 or more purebred lines were considered as high-confidence selection regions in the study. Of these 42 regions, 20 were 50 kb or less while 4 regions were larger than 0.5 Mb. In total, 91 genes could be found in the 42 regions, among which 19 regions contained only 1 or 2 genes, and 9 regions were located at gene deserts. CONCLUSIONS: Our results provide a genome-wide scan of recent selection signatures in five purebred lines of commercial broiler chickens. We found several candidate genes for recent selection in multiple lines, such as SOX6 (Sex Determining Region Y-Box 6) and cTR (Thyroid hormone receptor beta). These genes may have been under recent selection due to their essential roles in growth, development and reproduction in chickens. Furthermore, our results suggest that in some candidate regions, the same or opposite alleles have been under recent selection in multiple lines. Most of the candidate genes in the selection regions are novel, and as such they should be of great interest for future research into the genetic architecture of traits relevant to modern broiler breeding.

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.

Characterization of a novel chicken muscle disorder through differential gene expression and pathway analysis using RNA-sequencing.

BACKGROUND: Improvements in poultry production within the past 50 years have led to increased muscle yield and growth rate, which may be contributing to an increased rate and development of new muscle disorders in chickens. Previously reported muscle disorders and conditions are generally associated with poor meat quality traits and have a significant negative economic impact on the poultry industry. Recently, a novel myopathy phenotype has emerged which is characterized by palpably "hard" or tough breast muscle. The objective of this study is to identify the underlying biological mechanisms that contribute to this emerging muscle disorder colloquially referred to as "Wooden Breast", through the use of RNA-sequencing technology. METHODS: We constructed cDNA libraries from five affected and six unaffected breast muscle samples from a line of commercial broiler chickens. After paired-end sequencing of samples using the Illumina Hiseq platform, we used Tophat to align the resulting sequence reads to the chicken reference genome and then used Cufflinks to find significant changes in gene transcript expression between each group. By comparing our gene list to previously published histology findings on this disorder and using Ingenuity Pathways Analysis (IPA®), we aim to develop a characteristic gene expression profile for this novel disorder through analyzing genes, gene families, and predicted biological pathways. RESULTS: Over 1500 genes were differentially expressed between affected and unaffected birds. There was an average of approximately 98 million reads per sample, across all samples. Results from the IPA analysis suggested "Diseases and Disorders" such as connective tissue disorders, "Molecular and Cellular Functions" such as cellular assembly and organization, cellular function and maintenance, and cellular movement, "Physiological System Development and Function" such as tissue development, and embryonic development, and "Top Canonical Pathways" such as, coagulation system, axonal guidance signaling, and acute phase response signaling, are associated with the Wooden Breast disease. CONCLUSIONS: There is convincing evidence by RNA-seq analysis to support localized hypoxia, oxidative stress, increased intracellular calcium, as well as the possible presence of muscle fiber-type switching, as key features of Wooden Breast Disease, which are supported by reported microscopic lesions of the disease.

Linkage disequilibrium in crossbred and pure line chickens.

BACKGROUND: Both genome-wide association (GWA) studies and genomic selection depend on the level of non-random association of alleles at different loci, i.e. linkage disequilibrium (LD), across the genome. Therefore, characterizing LD is of fundamental importance to implement both approaches. In this study, using a 60K single nucleotide polymorphism (SNP) panel, we estimated LD and haplotype structure in crossbred broiler chickens and their component pure lines (one male and two female lines) and calculated the consistency of LD between these populations. RESULTS: The average level of LD (measured by r(2)) between adjacent SNPs across the chicken autosomes studied here ranged from 0.34 to 0.40 in the pure lines but was only 0.24 in the crossbred populations, with 28.4% of adjacent SNP pairs having an r(2) higher than 0.3. Compared with the pure lines, the crossbred populations consistently showed a lower level of LD, smaller haploblock sizes and lower haplotype homozygosity on macro-, intermediate and micro-chromosomes. Furthermore, correlations of LD between markers at short distances (0 to 10 kb) were high between crossbred and pure lines (0.83 to 0.94). CONCLUSIONS: Our results suggest that using crossbred populations instead of pure lines can be advantageous for high-resolution QTL (quantitative trait loci) mapping in GWA studies and to achieve good persistence of accuracy of genomic breeding values over generations in genomic selection. These results also provide useful information for the design and implementation of GWA studies and genomic selection using crossbred populations.

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.

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.

RNA-sequencing revisited data shed new light on wooden breast myopathy.

Wooden Breast (WB) abnormality represents one of the major challenges that the poultry industry has faced in the last 10 years. Despite the enormous progress in understanding the mechanisms underlying WB, the precise initial causes remain to be clarified. In this scenario, the present research is intended to characterize the gene expression profiles of broiler Pectoralis major muscles affected by WB, comparing them to the unaffected counterpart, to provide new insights into the biological mechanisms underlying this defect and potentially identifying novel genes likely involved in its occurrence. To this purpose, data obtained in a previous study through the RNA-sequencing technology have been used to identify differentially expressed genes (DEGs) between 6 affected and 5 unaffected broilers' breast muscles, by using the newest reference genome assembly for Gallus gallus (GRCg7b). Also, to deeply investigate molecular and biological pathways involved in the WB progression, pathways analyses have been performed. The results achieved through the differential gene expression analysis mainly evidenced the downregulation of glycogen metabolic processes, gluconeogenesis, and tricarboxylic acid cycle in WB muscles, thus corroborating the evidence of a dysregulated energy metabolism characterizing breasts affected by this abnormality. Also, genes related to hypertrophic muscle growth have been identified as differentially expressed (e.g., WFIKKN1). Together with that, a downregulation of genes involved in mitochondrial biogenesis and functionality has been detected. Among them, PPARGC1A and PPARGC1B chicken genes are particularly noteworthy. These genes not only have essential roles in regulating mitochondrial biogenesis but also play pivotal roles in maintaining glucose and energy homeostasis. In view of that, their downregulation in WB-affected muscle may be considered as potentially related to both the mitochondrial dysfunction and altered glucose metabolism in WB muscles, and their key involvement in the molecular alterations characterizing this muscular abnormality might be hypothesized.

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.

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.

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).

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.

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.

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.

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.

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.

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.