Transcriptome analysis of post-hatch breast muscle in legacy and modern broiler chickens reveals enrichment of several regulators of myogenic growth

Transcriptomics analysis unveils key potential genes associated with brain development and feeding behavior in the hypothalamus of L-citrulline-fed broiler chickens

High ambient temperature is a major environmental stressor affecting poultry production, especially in the tropical and subtropical regions of the world. Nutritional interventions have been adopted to combat thermal stress in poultry, including the use of amino acids. L-citrulline is a nonessential amino acid that is involved in nitric oxide generation and thermoregulation, however, the molecular mechanisms behind L-citrulline's regulation of body temperature are still unascertained. This study investigated the global gene expression in the hypothalamus of chickens fed either basal diet or L-citrulline-supplemented diets under different housing temperatures. Ross 308 broilers were fed with basal diet (CON) or 1% L-citrulline diet (LCT) from day-old, and later subjected to 2 environmental temperatures in a 2 by 2 factorial arrangement as follows; basal diet-fed chickens housed at 24°C (CON-TN); L-citrulline diet-fed chickens housed at 24°C (LCT-TN); basal diet-fed chickens housed at 35°C (CON-HS), and L-citrulline diet-fed chickens housed at 35°C (LCT-HS) from 22 to 42 d of age. At 42-days old, hypothalamic tissues were collected for mRNA analyses and RNA sequencing. A total of 1,019 million raw reads were generated and about 82.59 to 82.96% were uniquely mapped to genes. The gene ontology (GO) term between the CON-TN and LCT-TN groups revealed significant enrichments of pathways such as central nervous system development, and Wnt signaling pathway. On the other hand, GO terms between the CON-HS and LCT-HS groups revealed enrichments in the regulation of corticosteroid release, regulation of feeding behavior, and regulation of inflammatory response. Several potential candidate genes were identified to be responsible for central nervous system development (EMX2, WFIKKN2, SLC6A4 Wnt10a, and PHOX2B), and regulation of feed intake (NPY, AgRP, GAL, POMC, and NMU) in chickens. Therefore, this study unveils that L-citrulline can influence transcripts associated with brain development, feeding behavior, energy metabolism, and thermoregulation in chickens raised under different ambient temperatures.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

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.

In ovo probiotic supplementation promotes muscle growth and development in broiler embryos

In chickens, muscle development during embryonic growth is predominantly by myofiber hyperplasia. Following hatch, muscle growth primarily occurs via hypertrophy of the existing myofibers. Since myofiber number is set at hatch, production of more muscle fibers during embryonic growth would provide a greater myofiber number at hatch and potential for posthatch muscle growth by hypertrophy. Therefore, to improve performance in broilers, this study investigated the effect of in ovo spray application of probiotics on overall morphometry and muscle development in broiler embryos. For the study, fertile Ross 308 eggs were sprayed with different probiotics; Lactobacillus paracasei DUP 13076 (LP) and L. rhamnosus NRRL B 442 (LR) prior to and during incubation. The embryos were sacrificed on d 7, 10, 14, and 18 for embryo morphometry and pectoralis major muscle (PMM) sampling. Muscle sections were stained and imaged to quantify muscle fiber density (MFD), myofiber cross-sectional area (CSA), and nuclei density. Additionally, gene expression assays were performed to elucidate the effect of probiotics on myogenic genes. In ovo probiotic supplementation was found to significantly improve embryo weight, breast weight, and leg weight (P < 0.05). Further, histological analysis of PMM revealed a significant increase in MFD and nuclei number in the probiotic-treated embryos when compared to the control (P < 0.05). In 18-day-old broiler embryos, myofibers in the treatment group had a significantly smaller CSA (LP: 95.27 ± 3.28 μm², LR: 178.84 ± 15.1 μm²) when compared to the control (211.41 ± 15.67 μm²). This decrease in CSA was found to be associated with a concomitant increase in MFD (fibers/mm²) in the LP (13,647 ± 482.15) and LR (13,957 ± 463.13) group when compared to the control (7,680 ± 406.78). Additionally, this increase in myofibrillar hyperplasia in the treatment groups was associated with upregulation in the expression of key genes regulating muscle growth including MYF5, MYOD, MYOG, and IGF-1. In summary, in ovo spray application of probiotics promoted overall embryo growth and muscle development in broilers.

Chicken pituitary transcriptomic responses to acute heat stress

BACKGROUND

Poultry production is vulnerable to increasing temperatures in terms of animal welfare and in economic losses. With the predicted increase in global temperature and the number and severity of heat waves, it is important to understand how chickens raised for food respond to heat stress. This knowledge can be used to determine how to select chickens that are adapted to thermal challenge. As neuroendocrine organs, the hypothalamus and pituitary provide systemic regulation of the heat stress response.

METHODS AND RESULTS

Here we report a transcriptome analysis of the pituitary response to acute heat stress. Chickens were stressed for 2 h at 35 °C (HS) and transcriptomes compared with birds maintained in thermoneutral temperatures (25 °C).

CONCLUSIONS

The observations were evaluated in the context of ontology terms and pathways to describe the pituitary response to heat stress. The pituitaries of heat stressed birds exhibited responses to hyperthermia through altered expression of genes coding for chaperones, cell cycle regulators, cholesterol synthesis, transcription factors, along with the secreted peptide hormones, prolactin, and proopiomelanocortin.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Mining of chicken muscle growth genes and the function of important candidate gene RPL3L in muscle development

The birth weight of chickens does not significantly affect the weight at slaughter, while the different growth rate after birth was one of the important reasons for the difference in slaughter weight. Also, the increase in chickens' postnatal skeletal muscle weight is the main cause of the slaughter weight gain, but which genes are involved in this biological process is still unclear. In this study, by integrating four transcriptome datasets containing chicken muscles at different developmental times or different chicken tissues in public databases, a total of nine candidate genes that may be related to postnatal muscle development in chickens were obtained, including RPL3L, FBP2, ASB4, ASB15, CKMT2, PGAM1, YIPF7, PFKM, and LDHA. One of these candidate genes is RPL3L, whose 42 bp insertion/deletion (indel) mutation significantly correlated with multiple carcass traits in the F2 resource population from Xinghua chickens crossing with White Recessive Rock (WRR) chickens, including live weight, carcass weight, half eviscerated weight, eviscerated weight, breast meat weight, wing weight, leg muscle shear force, and breast muscle shear force. Also, there was a very significant difference between different genotypes of the RPL3L 42 bp indel mutation in these trains. Further experiments showed that RPL3L was highly expressed in chicken skeletal muscle, and its overexpression could promote the proliferation and inhibit the differentiation of chicken myoblasts by regulating ASB4 and ASB15 expression. Our findings demonstrated that the RPL3L 42 bp indel may be one of the molecular markers of chicken weight-related 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.

Thermal stress and selection for growth affect myogenic satellite cell lipid accumulation and adipogenic gene expression through mechanistic target of rapamycin pathway

Satellite cells (SCs) are multipotential stem cells having the plasticity to convert to an adipogenic lineage in response to thermal stress during the period of peak mitotic activity (the first week after hatch in poultry). The mechanistic target of rapamycin (mTOR) pathway, which regulates cellular function and fate of SCs, is greatly altered by thermal stress in turkey pectoralis major muscle SCs. The objective of the present study was to determine the effects of thermal stress, selection for growth, and the role of the mTOR pathway on SC intracellular lipid accumulation and expression of adipogenic regulatory genes. These effects were analyzed using SCs isolated from the pectoralis major muscle of 1-wk-old modern faster-growing commercial turkey line (NC) selected for increased growth and breast muscle yield as compared with SCs of a historic slower-growing Randombred Control Line 2 (RBC2) turkey. Heat stress (43 °C) of SCs during proliferation increased intracellular lipid accumulation (P < 0.001), whereas cold stress (33 °C) showed an inhibitory effect (P < 0.001) in both lines. Knockdown of mTOR reduced the intracellular lipid accumulation (P < 0.001) and suppressed the expression of several adipogenic regulatory genes: peroxisome proliferator-activated receptor-γ (PPARγ; P < 0.001), CCAAT/enhancer-binding protein-β (C/EBPβ; P < 0.001), and neuropeptide-Y (NPY; P < 0.001) during both proliferation and differentiation. The NC line SCs showed fewer reductions in lipid accumulation compared with the RBC2 line independent of temperature. Both intracellular lipid accumulation (P < 0.001) and PPARγ expression (P < 0.001) were greater at 72 h of proliferation than at 48 h of differentiation in both the RBC2 and NC lines independent of temperature. Thus, hot and cold thermal stress affected intracellular lipid accumulation in the pectoralis major muscle SCs, in part, through the mTOR pathway in wea growth-dependent manner. Altered intracellular lipid accumulation could eventually affect intramuscular fat deposition, resulting in a long-lasting effect on the structure and protein to fat ratio of the poultry pectoralis major muscle.

Thermal stress affects proliferation and differentiation of turkey satellite cells through the mTOR/S6K pathway in a growth-dependent manner

Satellite cells (SCs) are stem cells responsible for post-hatch muscle growth through hypertrophy and in birds are sensitive to thermal stress during the first week after hatch. The mechanistic target of rapamycin (mTOR) signaling pathway, which is highly responsive to thermal stress in differentiating turkey pectoralis major (p. major) muscle SCs, regulates protein synthesis and the activities of SCs through a downstream effector, S6 kinase (S6K). The objectives of this study were: 1) to determine the effect of heat (43°C) and cold (33°C) stress on activity of the mTOR/S6K pathway in SCs isolated from the p. major muscle of one-week-old faster-growing modern commercial (NC) turkeys compared to those from slower-growing Randombred Control Line 2 (RBC2) turkeys, and 2) to assess the effect of mTOR knockdown on the proliferation, differentiation, and expression of myogenic regulatory factors of the SCs. Heat stress increased phosphorylation of both mTOR and S6K in both turkey lines, with greater increases observed in the RBC2 line. With cold stress, greater reductions in mTOR and S6K phosphorylation were observed in the NC line. Early knockdown of mTOR decreased proliferation, differentiation, and expression of myoblast determination protein 1 and myogenin in both lines independent of temperature, with the RBC2 line showing greater reductions in proliferation and differentiation than the NC line at 38° and 43°C. Proliferating SCs are more dependent on mTOR/S6K-mediated regulation than differentiating SCs. Thus, thermal stress can affect breast muscle hypertrophic potential by changing satellite cell proliferation and differentiation, in part, through the mTOR/S6K pathway in a growth-dependent manner. These changes may result in irreversible effects on the development and growth of the turkey p. major muscle.

Transcriptomic Analysis of MSTN Knockout in the Early Differentiation of Chicken Fetal Myoblasts

In mammals, Myostatin (MSTN) is a known negative regulator of muscle growth and development, but its role in birds is poorly understood. To investigate the molecular mechanism of MSTN on muscle growth and development in chickens, we knocked out MSTN in chicken fetal myoblasts (CFMs) and sequenced the mRNA transcriptomes. The amplicon sequencing results show that the editing efficiency of the cells was 76%. The transcriptomic results showed that 296 differentially expressed genes were generated after down-regulation of MSTN, including angiotensin I converting enzyme (ACE), extracellular fatty acid-binding protein (EXFABP) and troponin T1, slow skeletal type (TNNT1). These genes are closely associated with myoblast differentiation, muscle growth and energy metabolism. Subsequent enrichment analysis showed that DEGs of CFMs were related to MAPK, Pl3K/Akt, and STAT3 signaling pathways. The MAPK and Pl3K/Akt signaling pathways are two of the three known signaling pathways involved in the biological effects of MSTN in mammals, and the STAT3 pathway is also significantly enriched in MSTN knock out chicken leg muscles. The results of this study will help to understand the possible molecular mechanism of MSTN regulating the early differentiation of CFMs and lay a foundation for further research on the molecular mechanism of MSTN involvement in muscle growth and development.

Domestication and Feed Restriction Programming Organ Index, Dopamine, and Hippocampal Transcriptome Profile in Chickens

The domestication process exerts different phenotypic plasticity between slow- and fast-growing breeds of chicken. Feed restriction has a critical role in production performance, physiological plasticity, and stress response. Our study aimed to explore how feed restriction programed the organ index, dopamine, and hippocampal transcriptome profile between slow- and fast-growing chickens, which were fed either ad libitum (SA and FA), or feed restricted to 70% of ad libitum (SR and FR), for 30 days. Results showed that feed restriction influenced the brain organ index (P < 0.05), but not the organ index of the heart, liver, and spleen. The slow-growing breed tested had a higher brain organ index than the fast-growing breed (P < 0.05). Under feed restriction conditions, both the slow- and fast-growing breeds had significantly elevated dopamine concentrations (P < 0.05) compared to those fed ad libitum. In the GO term, upregulated genes in the FA group were enriched in the mitochondria, respiratory chain, and energy metabolism compared to the SA group (P < 0.05). Membranes and ribosomes were enriched in the cellular component between the SR and FR groups (P < 0.05). In the KEGG functional pathways, upregulated DEGs in the FR group were enriched in the cardiovascular disease category and neurodegenerative disease category compared to the FA group (P < 0.05). Downregulated DEGs in the FA group were enriched in the oxidative phosphorylation and neurodegenerative disease categories (Parkinson's disease and Huntington's disease) compared with the SA group (P < 0.05). Upregulated DEGs in the FR group were enriched in the cardiovascular disease category, neurodegenerative disease category, and energy metabolism than the SR group (P < 0.05). In conclusion, feed restriction had profound effects on the brain organ index and plasma dopamine in the slow- and fast-growing chickens. Feed restriction may result in issues relating to cardiovascular and neurodegenerative diseases in the fast-growing breed tested, but not in the slow-growing breed.

Effective MSTN Gene Knockout by AdV-Delivered CRISPR/Cas9 in Postnatal Chick Leg Muscle

Muscle growth and development are important aspects of chicken meat production, but the underlying regulatory mechanisms remain unclear and need further exploration. CRISPR has been used for gene editing to study gene function in mice, but less has been done in chick muscles. To verify whether postnatal gene editing could be achieved in chick muscles and determine the transcriptomic changes, we knocked out Myostatin (MSTN), a potential inhibitor of muscle growth and development, in chicks and performed transcriptome analysis on knock-out (KO) muscles and wild-type (WT) muscles at two post-natal days: 3d (3-day-old) and 14d (14-day-old). Large fragment deletions of MSTN (>5 kb) were achieved in all KO muscles, and the MSTN gene expression was significantly downregulated at 14d. The transcriptomic results indicated the presence of 1339 differentially expressed genes (DEGs) between the 3d KO and 3d WT muscles, as well as 597 DEGs between 14d KO and 14d WT muscles. Many DEGs were found to be related to cell differentiation and proliferation, muscle growth and energy metabolism. This method provides a potential means of postnatal gene editing in chicks, and the results presented here could provide a basis for further investigation of the mechanisms involved in muscle growth and development.

Dynamic Transcriptomic Analysis of Breast Muscle Development From the Embryonic to Post-hatching Periods in Chickens

Skeletal muscle development and growth are closely associated with efficiency of poultry meat production and its quality. We performed whole transcriptome profiling based on RNA sequencing of breast muscle tissue obtained from Shouguang chickens at embryonic days (E) 12 and 17 to post-hatching days (D) 1, 14, 56, and 98. A total of 9,447 differentially expressed genes (DEGs) were filtered (Q < 0.01, fold change > 2). Time series expression profile clustering analysis identified five significantly different expression profiles that were divided into three clusters. DEGs from cluster I with downregulated pattern were significantly enriched in cell proliferation processes such as cell cycle, mitotic nuclear division, and DNA replication. DEGs from cluster II with upregulated pattern were significantly enriched in metabolic processes such as glycolysis/gluconeogenesis, insulin signaling pathway, calcium signaling pathway, and biosynthesis of amino acids. DEGs from cluster III, with a pattern that increased from E17 to D1 and then decreased from D1 to D14, mainly contributed to lipid metabolism. Therefore, this study may help us explain the mechanisms underlying the phenotype that myofiber hyperplasia occurs predominantly during embryogenesis and hypertrophy occurs mainly after birth at the transcriptional level. Moreover, lipid metabolism may contribute to the early muscle development and growth. These findings add to our knowledge of muscle development in chickens.

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.

Comparative transcriptomic analysis of hypothalamus-pituitary-liver axis in bighead carp (Hypophthalmichthys nobilis) with differential growth rate

Background

Growth rate is one of the most important features for aquaculture species and deciphering its regulation mechanism has great significance both in genetics and in economics. Hypothalamus-pituitary growth axis (HP growth axis) or neuro-endocrine axis plays a vital role in growth regulation in different aquaculture animals.

Results

In this study, the HP and liver transcriptomes of two female groups (H and L) with phenotypically extreme growth rate were sequenced using RNA-Seq. A total of 30,524 and 22,341 genes were found expressed in the two tissues, respectively. The average expression levels for the two tissues were almost the same, but the median differed significantly. A differential expression analysis between H and L groups identified 173 and 204 differentially expressed genes (DEGs) in HP and liver tissue, respectively. Pathway analysis revealed that DEGs in HP tissue were enriched in regulation of cell proliferation and angiogenesis while in liver tissue these genes were overrepresented in sterol biosynthesis and transportation. Genomic overlapping analyses found that 4 and 5 DEGs were within growth-related QTL in HP and liver tissue respectively. A deeper analysis of these 9 genes revealed 3 genes were functionally linked to the trait of interest. The expression of 2075 lncRNAs in HP tissue and 1490 in liver tissue were also detected, and some of lncRNAs were highly expressed in the two tissues.

Conclusions

Above all, the results of the present study greatly contributed to the knowledge of the regulation of growth and then assisted the design of new selection strategies for bighead carp with improved growth-related traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5691-4) contains supplementary material, which is available to authorized users.

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.

Transcriptome profile analysis of leg muscle tissues between slow- and fast-growing chickens

Chicken is widely favored by consumers because of some unique features. The leg muscles occupy an important position in the market. However, the specific mechanism for regulating muscle growth speed is not clear. In this experiment, we used Jinghai yellow chickens with different body weights at 300 days as research subjects. The chickens were divided into fast- and slow-growing groups, and we collected leg muscles after slaughtering for use in RNA-seq. After comparing the two groups, 87 differentially expressed genes (DEGs) were identified (fold change ≥ 2 and FDR < 0.05). The fast-growing group had 42 up-regulated genes and 45 down-regulated genes among these DEGs compared to the slow-growing group. Six items were significantly enriched in the biological process: embryo development ending in birth or egg hatching, chordate embryonic development, embryonic skeletal system development, and embryo development as well as responses to ketones and the sulfur compound biosynthetic process. Two significantly enriched pathways were found in the KEGG pathway analysis (P-value < 0.05): the insulin signaling pathway and the adipocytokine signaling pathway. This study provides a theoretical basis for the molecular mechanism of chicken growth and for improving the production of Jinghai yellow chicken.

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.

Longissimus dorsi muscle transcriptomic analysis of Yunling and Chinese simmental cattle differing in intramuscular fat content and fatty acid composition

Intramuscular fat (IMF) content and fatty acid (FA) composition vary significantly across beef cattle breeds, which play an important role in taste and nutritional value. However, the molecular mechanisms underlying these phenotypic differences remain unknown. The present study compared meat quality traits between Yunling cattle and Chinese Simmental cattle. Yunling cattle showed a lower IMF content and proportion of monounsaturated fatty acids (MUFA), as well as higher proportions of saturated fatty acids (SFA), polyunsaturated fatty acids (PUFA), and short-chain fatty acids (sc-FA) in the longissimus dorsi (LD) muscle than Chinese Simmental cattle. To further identify the candidate genes and pathways responsible for these phenotypic differences, the transcriptome of LD muscle from the two breeds were measured using RNA-seq. A total of 1347 differentially expressed genes were identified. The major metabolic pathways that were differentially modulated were lipolysis and glycometabolism. Yunling cattle showed a higher expression of lipolysis genes (ALDH9A1, ACSL5, ACADM, ACAT2, ACOT2) and a lower expression of genes related to glycometabolism (PGM1, GALM, PGM1, GPI, LDHA). This research identified candidate genes and pathways for IMF content and FA composition in the LD muscle of beef cattle, which may facilitate the design of new selection strategies to improve meat quality.

Transcriptome analysis of the pectoral muscles of local chickens and commercial broilers using Ribo-Zero ribonucleic acid sequencing

BACKGROUND

The molecular mechanisms underlying meat quality and muscle growth are not clear. The meat quality and growth rates of local chickens and commercial broilers are very different. The Ribo-Zero RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of this study was to provide a reference for studies of the differences in the meat quality and growth of different breeds of chickens.

RESULTS

Ribo-Zero RNA-Seq technology was used to analyze the pectoral muscle transcriptomes of Gushi chickens and AA broilers. Compared with AA broilers, 1649 genes with annotated information were significantly differentially expressed (736 upregulated and 913 downregulated) in Gushi chickens with Q≤0.05 (Q is the P-value corrected by multiple assumptions test) at a fold change ≥2 or ≤0.5. In addition, 2540 novel significantly differentially expressed (SDE) genes (1405 upregulated and 1135 downregulated) were discovered. The results showed that the main signal transduction pathways that differed between Gushi chickens and AA broilers were related to amino acid metabolism. Amino acids are important for protein synthesis, and they regulate key metabolic pathways to improve the growth, development and reproduction of organisms.

CONCLUSION

This study showed that differentially expressed genes in the pectoral tissues of Gushi chickens and AA broilers were related to fat metabolism, which affects meat. Additionally, a large number of novel genes were found that may be involved in fat metabolism and thus may affect the formation of meat, which requires further study. The results of this study provide a reference for further studies of the molecular mechanisms of meat formation.

Annotating long intergenic non-coding RNAs under artificial selection during chicken domestication

Background

Numerous biological functions of long intergenic non-coding RNAs (lincRNAs) have been identified. However, the contribution of lincRNAs to the domestication process has remained elusive. Following domestication from their wild ancestors, animals display substantial changes in many phenotypic traits. Therefore, it is possible that diverse molecular drivers play important roles in this process.

Results

We analyzed 821 transcriptomes in this study and annotated 4754 lincRNA genes in the chicken genome. Our population genomic analysis indicates that 419 lincRNAs potentially evolved during artificial selection related to the domestication of chicken, while a comparative transcriptomic analysis identified 68 lincRNAs that were differentially expressed under different conditions. We also found 47 lincRNAs linked to special phenotypes.

Conclusions

Our study provides a comprehensive view of the genome-wide landscape of lincRNAs in chicken. This will promote a better understanding of the roles of lincRNAs in domestication, and the genetic mechanisms associated with the artificial selection of domestic animals.

Electronic supplementary material

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

Association of growth rate with hormone levels and myogenic gene expression profile in broilers

Background

The growth rate often varies among individual broilers of the same breed under a common management condition. To investigate whether a variation in the growth rate is associated with a difference in hormone levels and myogenic gene expression profile in broilers, a feeding trial was conducted with 10,000 newly hatched Ross 308 chicks in a commercial production facility under standard management. At 38 d of age, 30 fast-, 30 medium-, and 30 slow-growing broilers were selected among 600 healthy male individuals. The levels of insulin-like growth factor-1 (IGF-1), triiodothyronine (T3), thyroxine (T4), and growth hormone in the serum or breast muscle were assayed by ELISA or RIA kits, and the expression levels of several representative pro- and anti-myogenic genes in the breast muscle were also measured by real-time PCR.

Results

Results showed that both absolute and relative weights of the breast muscle were in linear positive correlations with the body weight of broilers (P < 0.001). Fast-growing broilers had higher concentrations of IGF-1 than slow-growing broilers (P < 0.05) in both the serum and breast muscle. The serum concentration of T3 was significantly higher in fast-growing birds than in slow-growing birds (P < 0.05). However, no difference was observed in growth hormone or T4 concentration among three groups of birds. Additionally, a decreased expression of an anti-myogenic gene (myostatin) and increased expressions of pro-myogenic genes such as myogenic differentiation factor 1, myogenin, muscle regulatory factor 4, myogenic factor 5, IGF-1, and myocyte enhancer factor 2B, C, and D were observed in fast-growing broilers (P < 0.05), relative to slow-growing broilers.

Conclusions

Collectively, these findings suggested that the growth rate is linked to the hormone and myogenic gene expression levels in broiler chickens. Some of these parameters such as serum concentrations of IGF-1 and T3 could be employed to breed for enhanced growth.

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.

High concentrations of atmospheric ammonia induce alterations of gene expression in the breast muscle of broilers (Gallus gallus) based on RNA-Seq

BACKGROUND

High concentrations of atmospheric ammonia are one of the key environmental stressors affecting broiler production performance, which causes remarkable economic losses as well as potential welfare problems of the broiler industry. Previous reports demonstrated that high levels of ammonia altered body fat distribution and meat quality of broilers. However, the molecular mechanisms and metabolic pathways in breast muscle altered by high concentrations of ambient ammonia exposure on broilers are still unknown.

RESULTS

This study utilized RNA-Seq to compare the transcriptomes of breast muscles to identify differentially enriched genes in broilers exposed to high and low concentrations of atmospheric ammonia. A total of 267 promising candidate genes were identified by differential expression analysis, among which 67 genes were up-regulated and 189 genes were down-regulated. Bioinformatics analysis suggested that the up and down-regulation of these genes were involved in the following two categories of cellular pathways and metabolisms: Steroid biosynthesis (gga00100) and peroxisome proliferator-activated receptor (PPAR) signaling pathway (gga03320), which both participated in the lipid metabolism processes.

CONCLUSIONS

This study suggests that longtime exposure to high concentrations of aerial ammonia can change fat content in breast muscle, meat quality and palatability via altering expression level of genes participating in important lipid metabolism pathways. These findings have provided novel insights into our understanding of molecular mechanisms of breast muscles exposed to ammonia in broilers. This study provides new information that could be used for genetic breeding and nutritional intervention in production practice of broilers industry in the future.

WFIKKN1 and WFIKKN2: "Companion" proteins regulating TGFB activity

The WFIKKN (WAP, Follistatin/kazal, Immunoglobulin, Kunitz and Netrin domain-containing) protein family is composed of two multidomain proteins: WFIKKN1 and WFIKKN2. They were formed by domain shuffling and are likely present in deuterostoms. The WFIKKN (also called GASP) proteins are well known for their function in muscle and skeletal tissues, namely, inhibition of certain members of the transforming growth factor beta (TGFB) superfamily such as myostatin (MSTN) and growth and differentiation factor 11 (GDF11). However, the role of the WFIKKN proteins in other tissues is still poorly understood in spite of evidence suggesting possible action in the inner ear, brain and reproduction. Further, several recent studies based on next generation technologies revealed differential expression of WFIKKN1 and WFIKKN2 in various tissues suggesting that their function is not limited to MSTN and GDF11 inhibition in musculoskeletal tissue. In this review, we summarize current knowledge about the WFIKKN proteins and propose that they are "companion" proteins for various growth factors by providing localized and sustained presentation of TGFB proteins to their respective receptors, thus regulating the balance between the activation of Smad and non-Smad pathways by TGFB.