Whole transcriptome profiling reveals a lncMDP1 that regulates myogenesis by adsorbing miR-301a-5p targeting CHAC1.

Myoblast proliferation and differentiation are essential for skeletal muscle development. In this study, we generated the expression profiles of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in different developmental stages of chicken primary myoblasts (CPMs) using RNA sequencing (RNA-seq) technology. The dual luciferase reporter system was performed using chicken embryonic fibroblast cells (DF-1), and functional studies quantitative real-time polymerase chain reaction (qPCR), cell counting kit-8 (CCK-8), 5-Ethynyl-2'-deoxyuridine (EdU), flow cytometry cycle, RNA fluorescence in situ hybridization (RNA-FISH), immunofluorescence, and western blotting assay. Our research demonstrated that miR-301a-5p had a targeted binding ability to lncMDP1 and ChaC glutathione-specific gamma-glutamylcyclotransferase 1 (CHAC1). The results revealed that lncMDP1 regulated the proliferation and differentiation of myoblasts via regulating the miR-301a-5p/CHAC1 axis, and CHAC1 promotes muscle regeneration. This study fulfilled the molecular regulatory network of skeletal muscle development and providing an important theoretical reference for the future improvement of chicken meat performance and meat quality.

LncRNA lncMGR regulates skeletal muscle development and regeneration by recruiting CDK9 and sponging miRNAs.

Long non-coding RNAs (lncRNAs) play an essential role in vertebrate myogenesis and muscle diseases. However, the dynamic expression patterns, biological functions, and mechanisms of lncRNAs in skeletal muscle development and regeneration remain largely unknown. In this study, a novel lncRNA (named lncMGR) was differentially expressed during breast muscle development in fast- and slow-growing chickens. Functionally, lncMGR promoted myoblast differentiation, inhibited myoblast proliferation in vitro, and promoted myofiber hypertrophy and injury repair in vivo. Mechanistically, lncMGR increased the mRNA and protein expression of skeletal muscle myosin heavy chain 1 A (MYH1A) via both transcriptional and post-transcriptional regulation. Nuclear lncMGR recruited cyclin-dependent kinase 9 (CDK9) to the core transcriptional activation region of the MYH1A gene to activate MYH1A transcription. Cytoplasmic lncMGR served as a competitive endogenous RNA (ceRNA) to competitively absorb miR-2131-5p away from MYH1A and subsequently protected the MYH1A from miR-2131-5p-mediated degradation. Besides miR-2131-5p, cytoplasmic lncMGR could also sponge miR-143-3p to reconcile the antagonist between the miR-2131-5p/MYH1A-mediated inhibition effects and miR-143-3p-mediated promotion effects on myoblast proliferation, thereby inhibiting myoblast proliferation. Collectively, lncMGR could recruit CDK9 and sponge multiple miRNAs to regulate skeletal muscle development and regeneration, and could be a therapeutic target for muscle diseases.

Whole-transcriptome sequencing reveals a melanin-related ceRNA regulatory network in the breast muscle of Xichuan black-bone chicken.

The economic losses incurred due to reduced muscle pigmentation highlight the crucial role of melanin-based coloration in the meat of black-bone chickens. Melanogenesis in the breast muscle of black-bone chickens is currently poorly understood in terms of molecular mechanisms. This study employed whole-transcriptome sequencing to analyze black and white breast muscle samples from black-bone chickens, leading to the identification of 367 differentially expressed (DE) mRNAs, 48 DElncRNAs, 104 DEcircRNAs, and 112 DEmiRNAs involved in melanin deposition. Based on these findings, a competitive endogenous RNA (ceRNA) network was developed to better understand the complex mechanisms of melanin deposition. Furthermore, our analysis revealed key DEmRNAs (TYR, DCT, EDNRB, MLPH and OCA2) regulated by DEmiRNAs (gga-miR-140-5p, gga-miR-1682, gga-miR-3529, gga-miR-499-3p, novel-m0012-3p, gga-miR-200b-5p, gga-miR-203a, gga-miR-6651-5p, gga-miR-7455-3p, gga-miR-31-5p, miR-140-x, miR-455-x, novel-m0065-3p, gga-miR-29b-1-5p, miR-455-y, novel-m0085-3p, and gga-miR-196-1-3p). These DEmiRNAs competitively interacted with DElncRNAs including MSTRG.2609.2, MSTRG.4185.1, LOC112530666, LOC112533366, LOC771030, LOC107054724, LOC121107411, LOC100859072, LOC101750037, LOC121108550, LOC121109224, LOC121110876, and LOC101749016, as well as DEcircRNAs, such as novel_circ_000158, novel_circ_000623, novel_001518, and novel_circ_003596. The findings from this study provide insight into the mechanisms that regulate lncRNA, circRNA, miRNA, and mRNA expression in chicken melanin deposition.

RRM2 promotes the proliferation of chicken myoblasts, inhibits their differentiation and muscle regeneration.

During myogenesis and regeneration, the proliferation and differentiation of myoblasts play key regulatory roles and may be regulated by many genes. In this study, we analyzed the transcriptomic data of chicken primary myoblasts at different periods of proliferation and differentiation with protein‒protein interaction network, and the results indicated that there was an interaction between cyclin-dependent kinase 1 (CDK1) and ribonucleotide reductase regulatory subunit M2 (RRM2). Previous studies in mammals have a role for RRM2 in skeletal muscle development as well as cell growth, but the role of RRM2 in chicken is unclear. In this study, we investigated the effects of RRM2 on skeletal muscle development and regeneration in chickens in vitro and in vivo. The interaction between RRM2 and CDK1 was initially identified by co-immunoprecipitation and mass spectrometry. Through a dual luciferase reporter assay and quantitative real-time PCR, we identified the core promoter region of RRM2, which is regulated by the SP1 transcription factor. In this study, through cell counting kit-8 assays, 5-ethynyl-2'-deoxyuridine incorporation assays, flow cytometry, immunofluorescence staining, and Western blot analysis, we demonstrated that RRM2 promoted the proliferation and inhibited the differentiation of myoblasts. In vivo studies showed that RRM2 reduced the diameter of muscle fibers and slowed skeletal muscle regeneration. In conclusion, these data provide preliminary insights into the biological functions of RRM2 in chicken muscle development and skeletal muscle regeneration.

miR-19b-3p regulated by estrogen controls lipid synthesis through targeting MSMO1 and ELOVL5 in LMH cells.

miR-19b-3p is reported to undertake various biological role, while its function and action mechanism in chicken hepatic lipid metabolism is unclear. Conservation analysis and tissue expression pattern of miR-19b-3p and its target gene were evaluated, respectively. Dual luciferase reporter system and Western blot technologies were adopted to validate miR-19b-3p target gene. Overexpression and knockdown assays were done to explore the biological functions of miR-19b-3p and target gene in Leghorn Male Hepatoma cell line (LMH). Regulatory approaches of estrogen on miR-19b-3p and target gene expressions are analyzed through site-directed mutation combined with estrogen receptors antagonist treatment assays. The results showed that chicken miR-19b-3p mature sequences are highly conserved among Capra hircus, Columba livia, Rattus norvegicus, Mus musculus, Cricetulus griseus, Danio rerio, Danio novaehollandiae, Orycodylus porosus, Crocodylus porosus, Gadus morhua, and widely expressed in lung, ovary, spleen, duodenum, kidney, heart, liver, leg muscle, and pectoral muscle tissues. miR-19b-3p could significantly increase intracellular triglyceride (TG) content and decrease intracellular cholesterol (TC) content via targeting methylsterol monooxygenase 1 (MSMO1) and elongase of very long chain fatty acids 5 (ELOVL5), which are highly conserved among species, in both mRNA and protein levels. Estrogen could inhibit miR-19b-3p expression, but directly promoted MSMO1 transcription via estrogen receptor α (ERα) and indirectly regulated ELOVL5 expression at the transcription level. Meanwhile, estrogen could also upregulate MSMO1 and ELOVL5 expression through inhibiting miR-19b-3p expression at the post-transcription level. Taken together, these results highlight the role and regulatory mechanism of miR-19b-3p in hepatic lipid metabolism in chicken, and might produce useful comparative information for human obesity studies and biomedical research.

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.

Multi-Stage Transcriptome Analysis Revealed the Growth Mechanism of Feathers and Hair Follicles during Induction Molting by Fasting in the Late Stage of Egg Laying.

Induced molting is a common method to obtain a new life in laying hens, in which periodic changes in feathers are the prominent feature. Nevertheless, its precise molecular mechanism remains unclear. In this study, feather and hair follicle samples were collected during fasting-induced physiological remodeling for hematoxylin-eosin staining, hormone changes and follicle traits, and transcriptome sequencing. Feather shedding was observed in F13 to R25, while newborns were observed in R3 to R32. Triiodothyronine and tetraiodothyronine were significantly elevated during feather shedding. The calcium content was significantly higher, and the ash content was significantly lower after the changeover. The determination of hair follicle traits revealed an increasing trend in pore density and a decrease in pore diameter after the resumption of feeding. According to RNA-seq results, several core genes were identified, including DSP, CDH1, PKP1, and PPCKB, which may have an impact on hair follicle growth. The focus was to discover that starvation may trigger changes in thyroid hormones, which in turn regulate feather molting through thyroid hormone synthesis, calcium signaling, and thyroid hormone signaling pathways. These data provide a valuable resource for the analysis of the molecular mechanisms underlying the cyclical growth of hair follicles in the skin during induced molting.

LncHLEF promotes hepatic lipid synthesis through miR-2188-3p/GATA6 axis and encoding peptides and enhances intramuscular fat deposition via exosome.

Long noncoding RNAs (lncRNAs) have emergingly been implicated in mammalian lipid metabolism. However, their biological functions and regulatory mechanisms underlying adipogenesis remain largely elusive in chicken. Here, we systematically characterized the genome-wide full-length lncRNAs in the livers of pre- and peak-laying hens, and identified a novel intergenic lncRNA, lncHLEF, an RNA macromolecule with a calculated molecular weight of 433 kDa. lncHLEF was primarily distributed in cytoplasm of chicken hepatocyte and significantly up-regulated in livers of peak-laying hens. Functionally, lncHLEF could promote hepatocyte lipid droplet formation, triglycerides and total cholesterol contents. Mechanistically, lncHLEF could not only serve as a competitive endogenous RNA to modulate miR-2188-3p/GATA6 axis, but also encode three small functional polypeptides that directly interact with ACLY protein to enable its stabilization. Importantly, adeno-associated virus-mediated liver-specific lncHLEF overexpression resulted in increased hepatic lipid synthesis and intramuscular fat (IMF) deposition, but did not alter abdominal fat (AbF) deposition. Furthermore, hepatocyte lncHLEF could be delivered into intramuscular and abdominal preadipocytes via hepatocyte-secreted exosome to enhance intramuscular preadipocytes differentiation without altering abdominal preadipocytes differentiation. In conclusion, this study revealed that the lncHLEF could promote hepatic lipid synthesis through two independent regulatory mechanisms, and could enhance IMF deposition via hepatocyte-adipocyte communications mediated by exosome.

miR-128-3p inhibits intramuscular adipocytes differentiation in chickens by downregulating FDPS.

BACKGROUND: Intramuscular fat (IMF) content is the major indicator for evaluating chicken meat quality due to its positive correlation with tenderness, juiciness, and flavor. An increasing number of studies are focusing on the functions of microRNAs (miRNAs) in intramuscular adipocyte differentiation. However, little is known about the association of miR-128-3p with intramuscular adipocyte differentiation. Our previous RNA-seq results indicated that miR-128-3p was differentially expressed at different periods in chicken intramuscular adipocytes, revealing a possible association with intramuscular adipogenesis. The purpose of this research was to investigate the biological functions and regulatory mechanism of miR-128-3p in chicken intramuscular adipogenesis. RESULTS: The results of a series of assays confirmed that miR-128-3p could promote the proliferation and inhibit the differentiation of intramuscular adipocytes. A total of 223 and 1,050 differentially expressed genes (DEGs) were identified in the mimic treatment group and inhibitor treatment group, respectively, compared with the control group. Functional enrichment analysis revealed that the DEGs were involved in lipid metabolism-related pathways, such as the MAPK and TGF-ß signaling pathways. Furthermore, target gene prediction analysis showed that miR-128-3p can target many of the DEGs, such as FDPS, GGT5, TMEM37, and ASL2. The luciferase assay results showed that miR-128-3p targeted the 3' UTR of FDPS. The results of subsequent functional assays demonstrated that miR-128-3p acted as an inhibitor of intramuscular adipocyte differentiation by targeting FDPS. CONCLUSION: miR-128-3p inhibits chicken intramuscular adipocyte differentiation by downregulating FDPS. Our findings provide a theoretical basis for the study of lipid metabolism and reveal a potential target for molecular breeding to improve meat quality.

Integrated LC/MS-based lipidomics and transcriptomics analyses revealed lipid composition heterogeneity between pectoralis intramuscular fat and abdominal fat and its regulatory mechanism in chicken.

Intramuscular fat (IMF) content is conducive to multiple meat quality properties, while abdominal fat (AF) is treated as waste product in chicken industry. However, the heterogeneity and distinct regulatory mechanisms of lipid composition between the IMF and AF are still unclear. In this study, we carried out non-targeted lipidomics analyses of pectoralis IMF and AF, and detected a total of 423 differential lipid molecules (DLMs) between chicken IMF and AF, including 307 up-regulated and 116 down-regulated DLMs in pectoral IMF. These DLMs exhibited the definite alteration of lipid composition. The up-reglated DLMs in IMF were mainly glycerophospholipids (GPs), including the bulk of phosphatidylcholines (PC, PC (P) and PC (O)), phosphatidylethanolamines (PE, PE (P) and PE (O)), phosphatidylglycerols (PG) and phosphatidylinositol (PI), while the up-reglated DLMs in AF were mainly glycerolipids (GLs), including most of triacylglycerols (TG) and diacylglycerols (DG). We further identified 28 main DLMs contributing to the heterogeneous deposition of IMF and AF, including 11 TGs common to IMF and AF, 12 PCs/PC (P)s specific to IMF and 5 DGs specific to AF. Further integration of transcriptome with the main DLMs by weighted gene co-expression network analysis (WGCNA), we found five key gene sets that included 386 unique genes promoting IMF deposition in pectoralis, 213 unique genes promoting AF deposition, 6 unique genes detrimental to AF deposition, 7 common genes that promote IMF deposition in pectoralis while adversely affect AF deposition, and 28 genes that only promoted IMF deposition in pectoralis but had no effect on AF deposition. In addition, we also observed the expression characteristics of key genes in vivo and in vitro, and found that transmembrane protein family gene TMEM164 might be mainly involved in the positive regulation of intramuscular fat deposition in pectoralis and zinc finger protein family gene ZNF488 had a potential unique positive regulatory function on abdominal fat deposition. These findings provide new perspectives for understanding IMF and AF heterodeposition and will serve as a valuable information resource for improving meat quality via breeding selection in chicken.

Effects of SLC45A2 and GPNMB on Melanin Deposition Based on Transcriptome Sequencing in Chicken Feather Follicles.

As an essential genetic and economic trait, chicken feather color has long been an important research topic. To further understand the mechanism of melanin deposition associated with coloration in chicken feathers, we selected feather follicle tissues from the neck and wings of chickens with differently colored feathers (yellow, sub-Columbian, and silver) for transcriptome analysis. We focused on genes that were expressed in both the wings and neck and were expressed with the same trends in breeds with two different plumage colors, specifically, SLC45A2, GPNMB, MLPH, TYR, KIT, WNT11, and FZD1. GO and KEGG enrichment analyses showed the DEGs were enriched in melanin-related pathways, such as tyrosine metabolic pathway and melanogenesis, and PPI analysis highlighted the genes SLC45A2 and GPNMB as associated with melanin deposition. Verification experiments in chicken melanocytes demonstrated that these two genes promote melanocyte melanin deposition. These data enrich our knowledge of the mechanisms that regulate chicken feather color.

Transcriptomic analysis of mechanism underlying the effect of induced molting on semen quality and reproductive performance in aged Houdan roosters.

The reproductive performance of breeder roosters has significant economic importance in the poultry industry. Breeder roosters have severely reduced semen quality with age and will be at risk of culling in the following years. In order to extend the use of breeder roosters, we drew on the induced molting model of hens and selected 35 Houdan roosters aged 50 wk for induced molting. By comparing the body weight, testicular weight, semen quality, and reproductive performance before and after induced molting, we found that induced molting could restore the body weight and testicular weight to the levels before molting (P > 0.05). At the same time, it significantly improved sperm motility (P < 0.05) and also improved reproductive performance such as fertilization rate and hatching rate. To further reveal the mechanism underlying the effects of induced molting on semen quality and reproductive performance in aged Houdan roosters, we collected testes from 3 periods: 1 d before fasting (F0), 15 d after fasting (F15), and 32 d after recovery feeding (R32) for transcriptome sequencing analysis. A total of 5,671 genes were detected in F0, F15, and R32, and trend analysis of the 5,671 differential genes showed 2 significant trends (profile 5 and profile 2). KEGG enrichment analysis of the genes in the 2 profiles, revealed significantly enriched pathway regulation of actin cytoskeleton. In the regulation of actin cytoskeleton pathway, we found a protein kinase gene (SRC) and a senescence gene (ROCK2). SRC was highly expressed at F15, leading to the phosphorylation of key substrates, which in turn disrupted the Sertoli cell spermatid connection and the spermiogenesis process, resulting in no mature spermatozoa produced from F15, SRC expression was inhibited at R32, the expression level was reduced, and mature spermatozoa reappeared. The senescence gene ROCK2 was highly expressed at F15 compared to F0 and R32, which may have been responsible for inducing senescence atrophy in the testes.

Transcriptome-Based Identification of the Muscle Tissue-Specific Expression Gene CKM and Its Regulation of Proliferation, Apoptosis and Differentiation in Chicken Primary Myoblasts.

Skeletal muscle is an essential tissue in meat-producing animals, and meat-producing traits have been a hot topic in chicken genetic breeding research. Current research shows that creatine kinase M-type-like (CKM) is one of the most abundant proteins in skeletal muscle and plays an important role in the growth and development of skeletal muscle, but its role in the development of chicken skeletal muscle is still unclear. Via RNA sequencing (RNA-seq), we found that CKM was highly expressed in chicken breast muscle tissue. In this study, the expression profile of CKM was examined by quantitative real-time PCR (qPCR), and overexpression and RNA interference techniques were used to explore the functions of CKM in the proliferation, apoptosis and differentiation of chicken primary myoblasts (CPMs). It was shown that CKM was specifically highly expressed in breast muscle and leg muscle and was highly expressed in stage 16 embryonic muscle, while CKM inhibited proliferation, promoted the apoptosis and differentiation of CPMs and was involved in regulating chicken myogenesis. Transcriptome sequencing was used to identify genes that were differentially expressed in CPMs after CKM disruption, and bioinformatics analysis showed that CKM was involved in regulating chicken myogenesis. In summary, CKM plays an important role in skeletal muscle development during chicken growth and development.

Effect of induced molting on ovarian function remodeling in laying hens.

Induced molting (IM) can restore the laying rate of aged laying hens to the peak level of laying and rejuvenate ovarian function for the second laying cycle. To explore the mechanism of ovarian function remodeling during IM in laying hens, in this study, ninety 71-wk-old laying lady hens with 60% laying rate and uniform weight were selected for molting induction by fasting. Samples (serum and fresh ovarian tissue) were collected on the day before fasting (F0), the 3rd and 16th days of fasting (F3, F16), and the 6th, 15th, 32nd days of refeeding (R6, R15, and R32), and the number of follicles in each period was counted. Then, the reproductive hormone levels in serum and antioxidant levels in ovarian tissues were detected at different stages, and the gene expression of the KIT-PI3K-PTEN-AKT pathway and GDF-9 in ovaries was measured by qRT-PCR. The results showed that the laying rate increased rapidly after refeeding and returned to the prefasting level by R32. At F16 and R6, the number of mature follicles significantly decreased; the number of primary and secondary follicles significantly increased; the contents of follicle stimulating hormone (FSH), luteinizing hormone (LH), estradiol (E2), and progesterone (P4) in serum decreased; the relative expression of KIT, PI3K, AKT, and GDF-9 significantly increased; and that of PTEN significantly decreased. At R15 and R32, except for GDF-9, which maintained a high expression state, other indicators showed opposing trends to those observed at F16 and R6. In conclusion, IM activated the KIT-PI3K-PTEN-AKT signaling pathway and promoted the activation of primordial follicles during the fasting period and early resumption of feeding; gonadotropin secretion increased gradually, which promoted the rapid development of primary and secondary follicles to mature follicles and ovulation. This study explained the mechanism of ovarian function remodeling in the process of IM and provided a theoretical basis for improving the ovarian function of laying hens and optimizing the IM program.

Research Note: Combined analysis of BSA-seq based mapping and RNA-seq reveals candidate genes associated with sub-Columbian plumage in H line chickens.

Columbian coloration patterns in plumage are widespread phenomena in several standard breeds of poultry, such as the Columbian Plymouth Rock chicken. H line chicken plumage is generally a pure white except in the hackle, wing, and tail plumage, which coloration is very similar to the Columbian plumage pattern, but with the barring substituting for the black vertical striping. Thus, we refer to this plumage coloration as "sub-Columbian" pattern. However, the genetic basis of this phenotype remains unknown. Here, a F3 cross population between yellow plumage roosters and sub-Columbian plumage hens was constructed, for verifying sub-Columbian plumage was sex-linked dominant inheritance. To identify the candidate regions, F2 generation sub-Columbian plumage hens and yellow plumage hens with their parental lines were used for BSA-seq, and sub-Columbian plumage genes were mapped to a 10.46 Mb interval on chromosome Z. Remarkably, by transcriptome analysis of the neck and wing tip follicle tissues of the 2 plumage colors, we demonstrated that within the interval, only 1 gene, SLC45A2 expressed significant differently (P < 0.05). Through KASP, we identified L347M and A10331272T in solute carrier family 45 member 2 (SLC45A2), and B2 haplotype of cyclin-dependent kinase inhibitor 2A (CDKN2A), showed near complete association with the phenotype. Eventually, we designed a hybridization experiment for verifying the locus of sub-Columbian plumage, which is inherited through Z-linked dominant inheritance and is controlled by SLC45A2 and CDKN2A.

Carcass composition, meat quality, leg muscle status, and its mRNA expression profile in broilers affected by valgus-varus deformity.

Valgus-varus deformity (VVD) is a common leg disease in commercial broilers, which seriously affects animal welfare and causes economic losses. Up to now, most of the studies on VVD have been on skeleton, whereas there are fewer studies on VVD muscle. In this study, carcass composition and meat quality of 35-day-old normal and VVD Cobb broilers assess the effect of VVD on broiler growth. Molecular biology, morphology, and RNA sequencing (RNA-seq) were used to study the difference between normal and VVD gastrocnemius muscle. In comparison with the normal broilers, the breast muscle and leg muscle of the VVD broilers had lower shear force, notably lower crude protein, lower water content, cooking loss, and deeper meat color (P < 0.05). The morphological results showed that the weight of skeletal muscle was significantly higher in the normal broilers than that in the VVD broilers (P < 0.01), the diameter and area of myofibrils in the affected VVD were smaller than in the normal broilers (P < 0.01). Quantitative real-time PCR (qPCR) of gastrocnemius muscle revealed that the expression of myasthenic marker genes, fast myofiber marker genes, and apoptosis-related factors were significantly higher in the VVD broilers than in the normal broilers (P < 0.01). In total, 736 differentially expressed genes (DEGs) were identified firstly in the normal and VVD leg muscle by RNA-seq. Gene ontology (GO) enrichment indicated that these DEGs were mainly involved in the multicellular organismal process and anatomical structure development. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that DEGs are significantly enriched in proteasome. Protein interaction analysis obtained that DEGs with high interaction were proteasome-related coding genes and ubiquitin-related genes, these DEGs were closely associated with muscle atrophy. These show that VVD has an adverse effect on growth characteristics, slaughter characteristics, and meat quality in broilers, which may cause leg muscle atrophy. This study provides some reference values and basis for studying the pathogenesis of VVD in broilers.

Genome-wide association study of 17 serum biochemical indicators in a chicken F2 resource population.

BACKGROUND: Serum biochemical indicators are often regarded as direct reflections of animal metabolism and health. The molecular mechanisms underlying serum biochemical indicators metabolism of chicken (Gallus Gallus) have not been elucidated. Herein, we performed a genome-wide association study (GWAS) to identify the variation associated with serum biochemical indicators. The aim of this research was to broaden the understanding of the serum biochemical indicators in chickens. RESULTS: A GWAS of serum biochemical indicators was carried out on 734 samples from an F2 Gushi× Anka chicken population. All chickens were genotyped by sequencing, 734 chickens and 321,314 variants were obtained after quality control. Based on these variants, a total of 236 single-nucleotide polymorphisms (SNPs) on 9 chicken chromosomes (GGAs) were identified to be significantly (-log10(P) > 5.72) associated with eight of seventeen serum biochemical indicators. Ten novel quantitative trait locis (QTLs) were identified for the 8 serum biochemical indicator traits of the F2 population. Literature mining revealed that the ALPL, BCHE, GGT2/GGT5 genes at loci GGA24, GGA9 and GGA15 might affect the alkaline phosphatase (AKP), cholinesterase (CHE) and γ-glutamyl transpeptidase (GGT) traits, respectively. CONCLUSION: The findings of the present study may contribute to a better understanding of the molecular mechanisms of chicken serum biochemical indicator regulation and provide a theoretical basis for chicken breeding programs.

Population Structure and Genetic Diversity Analysis of "Yufen 1" H Line Chickens Using Whole-Genome Resequencing.

The effective protection and utilization of poultry resources depend on an accurate understanding of the genetic diversity and population structure. The breeding of the specialized poultry lineage "Yufen 1", with its defined characteristics, was approved by the China Poultry Genetic Resource Committee in 2015. Thus, to investigate the relationship between the progenitor H line and other poultry breeds, the genetic diversity and population structure of "Yufen 1" H line (YF) were investigated and compared with those of 2 commercial chicken breeds, the ancestor breed Red Jungle Fowls, and 11 Chinese indigenous chicken breeds based on a whole-genome resequencing approach using 8,112,424 SNPs. The genetic diversity of YF was low, and the rate of linkage disequilibrium decay was significantly slower than that of the other Chinese indigenous breeds. In addition, it was shown that the YF population was strongly selected during intensive breeding and that genetic resources have been seriously threatened, which highlights the need to establish a systematic conservation strategy as well as utilization techniques to maintain genetic diversity within YF. Moreover, a principal component analysis, a neighbor-joining tree analysis, a structure analysis, and genetic differentiation indices indicated that YF harbors a distinctive genetic resource with a unique genetic structure separate from that of Chinese indigenous breeds at the genome level. The findings provide a valuable resource and the theoretical basis for the further conservation and utilization of YF.

Integration Analysis of circRNA-miRNA-mRNA and Identification of Critical Networks in Valgus-Varus Deformity (Gallus gallus).

Valgus-valgus deformity (VVD) is a common leg deformity in broilers with inward or outward deviation of the tibiotarsus and tarsometatarsus. The competing endogenous RNA (ceRNA) network plays an essential role in the study of leg disease. However, its role in the etiology and pathogenesis of VVD remains unclear. Here, based on case (VVD) and control (normal) group design, we performed analyses of differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs) and differentially expressed mRNAs (DEmRNAs). Transcriptome data derived 86 DEcircRNAs, 13 DEmiRNAs and 410 DEmRNAs. Functional analysis showed that DEmRNAs were significantly enriched in cell cycle, apoptosis, ECM-receptor interaction, FoxO signaling pathway and protein processing synthesis. DEcirc/miRNA-associated DEmRNAs were associated with skeletal and muscle growth and development pathways, including mTOR, Wnt, and VEGF signaling pathways. Subsequently, a circRNA-miRNA-mRNA regulatory network was constructed based on the ceRNA hypothesis, including 8 circRNAs, 6 miRNAs, and 31 mRNAs, which were significantly enriched in the skeletal developmental pathway. Finally, two key mRNAs (CDC20 and CTNNB1) and their regulatory axes were screened by the PPI network and cytohubba. The expression levels of CDC20 and CTNNB1 in cartilage and seven other tissues were also quantified by qPCR. In conclusion, we analyzed the functions of DEmRNA, DEcircRNA and DEmiRNA and constructed the hub ceRNA regulatory axis, and obtained two hub genes, CDC20 and CTNNB1. The study more deeply explored the etiology and pathogenesis of VVD and lays the foundation for further study of the role of the ceRNA network on skeletal development.

Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds.

There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.