Search

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

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

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

Abstract Background In molecular biology studies, the selection of optimal reference genes is of vital importance for accurately quantifying gene expression. The purpose of the present study was to screen the most stable reference genes in different muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits. Methods and results Results indicated that the most stable reference genes in the muscle tissues of New Zealand white rabbits were HPRT1, ACTB and PPIC, while HPRT1, PPIC, and RPL13A were the most stable reference genes in muscle tissues of Yufeng yellow rabbits. However, in the longissimus dorsi muscle and the abdominal wall muscle of both varieties, the most stable reference genes were HPRT1, RPL13A, and SDHA. In the quadriceps femoris muscle, the most stable reference genes were ACTB, HPRT1, and SDHA. Furthermore, the relative abundance of MYOG, MYH3 and MSTN was used to confirm the suitability and reliability of the selected most stable reference genes and the most unstable reference gene. Results revealed the same expression patterns of these myogenic genes when normalized according to the most stable genes, while normalization against the unstable reference gene altered the observed expression patterns. Conclusions Taken together, our results demonstrated that the most stable reference genes varied among different muscle tissues and different breeds of rabbits. However, HPRT1, PPIC and SDHA presented high stability among all examined reference genes; thus, the combined analysis of HPRT1/ PPIC/ SDHA gene provides the best reference for RT-qPCR in muscle tissues of New Zealand white rabbits and Yufeng yellow rabbits, while HPRT1 is a better choice than other reference genes when using a single reference gene to assess target gene expression. Our results provide basic data for better measuring target gene expression profiles in muscle tissues of rabbits.

Abstract Background Elongases of very long chain fatty acids (ELOVLs), a family of first rate-limiting enzymes in the synthesis of long-chain fatty acids, play an essential role in the biosynthesis of complex lipids. Disrupting any of ELOVLs affects normal growth and development in mammals. Genetic variations in ELOVLs are associated with backfat or intramuscular fatty acid composition in livestock. However, the effects of ELOVL gene family on breeding selection and lipid deposition in different tissues are still unknown in chickens. Results Genetic variation patterns and genetic associations analysis showed that the genetic variations of ELOVL genes were contributed to breeding selection of commercial varieties in chicken, and 14 SNPs in ELOVL2-6 were associated with body weight, carcass or fat deposition traits. Especially, one SNP rs17631638T > C in the promoter of ELOVL3 was associated with intramuscular fat content (IMF), and its allele frequency was significantly higher in native and layer breeds compared to that in commercial broiler breeds. Quantitative real-time PCR (qRT-PCR) determined that the ELOVL3 expressions in pectoralis were affected by the genotypes of rs17631638T > C. In addition, the transcription levels of ELOVL genes except ELOVL5 were regulated by estrogen in chicken liver and hypothalamus with different regulatory pathways. The expression levels of ELOVL1-6 in hypothalamus, liver, abdominal fat and pectoralis were correlated with abdominal fat weight, abdominal fat percentage, liver lipid content and IMF. Noteworthily, expression of ELOVL3 in pectoralis was highly positively correlated with IMF and glycerophospholipid molecules, including phosphatidyl choline, phosphatidyl ethanolamine, phosphatidyl glycerol and phospholipids inositol, rich in ω-3 and ω-6 long-chain unsaturated fatty acids, suggesting ELOVL3 could contribute to intramuscular fat deposition by increasing the proportion of long-chain unsaturated glycerophospholipid molecules in pectoralis. Conclusions In summary, we demonstrated the genetic contribution of ELOVL gene family to breeding selection for specialized varieties, and revealed the expression regulation of ELOVL genes and their potential roles in regulating lipid deposition in different tissues. This study provides new insights into understanding the functions of ELOVL family on avian growth and lipid deposition in different tissues and the genetic variation in ELOVL3 may aid the marker-assisted selection of meat quality in chicken.

Abstract Background Molecular breeding accelerates the speed of animal breeding. Screening molecular markers that can affect economic traits through genome-wide association studies (GWAS) can provide a theoretical basis for molecular breeding. At present, a large number of molecular markers have been screened in poultry research, but few reports on how molecular markers affect economic traits exist. It is particularly important to reveal the action mechanisms of molecular markers, which can provide more accurate information for molecular breeding. Results The aim of this study was to investigate the relationships between two indels (NUDT15-indel-2777 and NUDT15-indel-1673) in the promoter region of NUDT15 and growth and carcass traits in chickens and to explore the regulatory mechanism of NUDT15. Significant differences were found in genotype and allele frequencies among commercial broilers, commercial laying hens and dual-purpose chickens. The results of association analyses showed that these two indel loci could significantly affect growth traits, such as body weight, and carcass traits. Tissue expression profiling at E12 showed that the expression of NUDT15 was significantly higher in skeletal muscle, and time-expression profiling of leg muscle showed that the expression of NUDT15 in myoblasts was significantly higher in the E10 and E12 proliferation stages than in other stages. Promoter activity analysis showed that pro-1673-I and pro-1673-D significantly inhibited promoter activity, and the promoter activity of pro-1673-D was significantly lower than that of pro-1673-I. In addition, when NUDT15 was overexpressed or underwent interference in chicken primary myoblasts (CPMs), NUDT15 could inhibit the proliferation of CPMs. Conclusion The results suggest that the studied indels in the promoter region of NUDT15 may regulate the proliferation of CPMs by affecting NUDT15 expression, ultimately affecting the growth and carcass traits of chickens. These indel polymorphisms may be used together as molecular markers for improving economic traits in chickens.

Abstract Background Valgus-varus deformity (VVD) is a lateral or middle deviation of the tibiotarsus or tarsometatarsus, which is associated with compromised growth, worse bone quality and abnormal changes in serum indicators in broilers. To investigate the genetic basis of VVD, a genome wide association study (GWAS) was performed to identify candidate genes and pathways that are responsible for VVD leg disease, serum indicators and growth performance in broilers. Results In total, VVD phenotype, seven serum indicators and three growth traits were measured for 126 VVD broilers (case group) and 122 sound broilers (control group) based on a high throughput genome wide genotyping-by-sequencing (GBS) method. After quality control 233 samples (113 sound broilers and 120 VVD birds) and 256,599 single nucleotide polymorphisms (SNPs) markers were used for further analysis. As a result, a total of 5 SNPs were detected suggestively significantly associated with VVD and 70 candidate genes were identified that included or adjacent to these significant SNPs. In addition, 43 SNPs located on Chr24 (0.22 Mb - 1.79 Mb) were genome-wide significantly associated with serum alkaline phosphatase (ALP) and 38 candidate genes were identified. Functional enrichment analysis showed that these genes are involved in two Gene Ontology (GO) terms related to bone development (cartilage development and cartilage condensation) and two pathways related to skeletal development (Toll−like receptor signaling pathway and p53 signaling pathway). BARX2 (BARX homeobox 2) and Panx3 (Pannexin 3) related to skeleton diseases and bone quality were obtained according to functional analysis. According to the integration of GWAS with transcriptome analysis, HYLS1 (HYLS1 centriolar and ciliogenesis associated) was an important susceptibility gene. Conclusions The results provide some reference for understanding the relationship between metabolic mechanism of ALP and pathogenesis of VVD, which will provide a theoretical basis for disease-resistant breeding of chicken leg soundness.

Abstract Background Estrogen plays an essential role in female development and reproductive function. In chickens, estrogen is critical for lipid metabolism in the liver. The regulatory molecular network of estrogen in chicken liver is poorly understood. To identify estrogen-responsive genes and estrogen functional sites on a genome-wide scale, we determined expression profiles of mRNAs, lncRNAs, and miRNAs in estrogen-treated ((17β-estradiol)) and control chicken livers using RNA-Sequencing (RNA-Seq) and studied the estrogen receptor α binding sites by ChIP-Sequencing (ChIP-Seq). Results We identified a total of 990 estrogen-responsive genes, including 962 protein-coding genes, 11 miRNAs, and 17 lncRNAs. Functional enrichment analyses showed that the estrogen-responsive genes were highly enriched in lipid metabolism and biological processes. Integrated analysis of the data of RNA-Seq and ChIP-Seq, identified 191 genes directly targeted by estrogen, including 185 protein-coding genes, 4 miRNAs, and 2 lncRNAs. In vivo and in vitro experiments showed that estrogen decreased the mRNA expression of PPARGC1B, which had been reported to be linked with lipid metabolism, by directly increasing the expression of miR-144-3p. Conclusions These results increase our understanding of the functional network of estrogen in chicken liver and also reveal aspects of the molecular mechanism of estrogen-related lipid metabolism.

Abstract Background Chicken skeletal muscle is an important economic product. The late stages of chicken development constitute the main period that affects meat production. LncRNAs play important roles in controlling the epigenetic process of growth and development. However, studies on the role of lncRNAs in the late stages of chicken breast muscle development are still lacking. In this study, to investigate the expression characteristics of lncRNAs during chicken muscle development, 12 cDNA libraries were constructed from Gushi chicken breast muscle samples from 6-, 14-, 22-, and 30-week-old chickens. Results A total of 1252 new lncRNAs and 1376 annotated lncRNAs were identified. Furthermore, 53, 61, 50, 153, 117, and 78 DE-lncRNAs were found in the W14 vs. W6, W22 vs. W14, W22 vs. W6, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparison groups, respectively. After GO enrichment analysis of the DE-lncRNAs, several muscle development-related GO terms were found in the W22 vs. W14 comparison group. Moreover, it was found that the MAPK signaling pathway was one of the most frequently enriched pathways in the different comparison groups. In addition, 12 common target DE-miRNAs of DE-lncRNAs were found in different comparison groups, some of which were muscle-specific miRNAs, such as gga-miR-206, gga-miR-1a-3p, and miR-133a-3p. Interestingly, the precursors of four newly identified miRNAs were found to be homologous to lncRNAs. Additionally, we found some ceRNA networks associated with muscle development-related GO terms. For example, the ceRNA networks contained the DYNLL2 gene with 12 lncRNAs that targeted 2 miRNAs. We also constructed PPI networks, such as IGF-I-EGF and FZD6-WNT11. Conclusions This study revealed, for the first time, the dynamic changes in lncRNA expression in Gushi chicken breast muscle at different periods and revealed that the MAPK signaling pathway plays a vital role in muscle development. Furthermore, MEF2C and its target lncRNA may be involved in muscle regulation through the MAPK signaling pathway. This research provided valuable resources for elucidating posttranscriptional regulatory mechanisms to promote the development of chicken breast muscles after hatching.

Abstract Background Domesticated chickens have a wide variety of phenotypes, in contrast with their wild progenitors. Unlike other chicken breeds, Xichuan black-bone chickens have blue-shelled eggs, and black meat, beaks, skin, bones, and legs. The breeding history and the economically important traits of this breed have not yet been explored at the genomic level. We therefore used whole genome resequencing to analyze the breeding history of the Xichuan black-bone chickens and to identify genes responsible for its unique phenotype. Results Principal component and population structure analysis showed that Xichuan black-bone chicken is in a distinct clade apart from eight other breeds. Linkage disequilibrium analysis showed that the selection intensity of Xichuan black-bone chickens is higher than for other chicken breeds. The estimated time of divergence between the Xichuan black-bone chickens and other breeds is 2.89 ka years ago. Fst analysis identified a selective sweep that contains genes related to melanogenesis. This region is probably associated with the black skin of the Xichuan black-bone chickens and may be the product of long-term artificial selection. A combined analysis of genomic and transcriptomic data suggests that the candidate gene related to the black-bone trait, EDN3, might interact with the upstream ncRNA LOC101747896 to generate black skin color during melanogenesis. Conclusions These findings help explain the unique genetic and phenotypic characteristics of Xichuan black-bone chickens, and provide basic research data for studying melanin deposition in animals.

Abstract Background Studies have shown that some viral infections cause structural changes in the intestinal microflora, but little is known about the effects of tumorigenic viral infection on the intestinal microflora of chickens. Results A 29-week commercial layer flock positive for avian leukosis virus-J (ALV-J), Marek’s disease virus (MDV) and avian reticuloendotheliosis virus (REV) was selected, and fresh fecal samples were collected and examined for the composition of the gut microflora by Illumina sequencing of the V3-V4 region of the 16S rRNA gene. The operational taxonomic units (OTUs) of the fecal microbiota differentiated the chickens infected with only ALV-J and those coinfected with ALV-J and MDV or REV from infection-negative chickens. The enrichment and diversity of cloacal microflora in chickens infected with ALV-J alone were slightly different from those in the infection-negative chickens. However, the diversity of cloacal microflora was significantly increased in chickens coinfected with both ALV-J and MDV or REV. Conclusions The intestinal microbiota was more strongly disturbed in chickens after coinfection with ALV-J and MDV or REV than after infection with ALV-J alone, and there may be underlying mechanisms by which the capacity for the stabilization of the intestinal flora was impaired due to viral infection and tumorigenesis.

Abstract Background The SH3RF2 gene is a protein-coding gene located in a quantitative trait locus associated with body weight, and its deletion has been shown to be positively associated with body weight in chickens. Results In the present study, CNV in the SH3RF2 gene was detected in 4079 individuals from 17 populations, including the “Gushi ×Anka” F2 resource population and populations of Chinese native chickens, commercial layers, and commercial broilers. The F2 resource population was then used to investigate the genetic effects of the chicken SH3RF2 gene. The results showed that the local chickens and commercial layers were all homozygous for the wild-type allele. Deletion mutation individuals were detected in all of the commercial broiler breeds except Hubbard broiler. A total of, 798 individuals in the F2 resource group were used to analyze the effects of genotype (DD/ID/II) on chicken production traits. The results showed that CNV was associated with 2-, 6-, 10-, and 12-week body weight (P = 0.026, 0.042, 0.021 and 0.039 respectively) and significantly associated with 8-week breast bone length (P = 0.045). The mutation was significantly associated with 8-week body weight (P = 0.007) and 4-week breast bone length (P = 0.010). CNV was significantly associated with evisceration weight, leg muscle weight, carcass weight, breast muscle weight and gizzard weight (P = 0.032, 0.033, 0.045, 0.004 and 0.000, respectively). Conclusions CNV of the SH3RF2 gene contributed to variation in the growth and weight gain of chickens.

Abstract Background The molecular mechanisms underlying stress-influenced immune function of chicken (Gallus Gallus) are not clear. The stress models can be established effectively by feeding chickens corticosterone (CORT) hormone. The bursa of Fabricius is a unique central immune organ of birds. RNA-Seq technology was used to investigate differences in the expression profiles of immune-related genes and associated pathways in the bursa of Fabricius to clarify molecular mechanisms. The aim of this study was to broaden the understanding of the stress-influenced immune function in chickens. Results Differentially expressed genes (DEGs) in the bursa of Fabricius between experimental group (basal diet with added CORT 30 mg/kg; C_B group) and control group (basal diet; B_B group) were identified by using RNA-seq technology. In total, we found 1434 significant DEGs (SDEGs), which included 199 upregulated and 1235 downregulated genes in the C_B group compared with the B_B group. The immune system process GO term was the top significantly GO term, including MYD88, TLR4, IL15, VEGFA gene and so on. The cytokine-cytokine receptor interaction pathway and the Toll-like receptor signaling pathway were the key pathways affected by stress. The protein-protein interaction (PPI) analysis of the SDEGs showed that VEGFA, MyD88 and IL15 were hub genes and module analysis showed that MYD88, TLR4 and VEGFA play important roles in response to stress. Conclusion This study showed that the VEGFA and ILs (such as IL15) via the cytokine-cytokine receptor interaction pathway, MYD88 and TLR4 via the Toll-like receptor signaling pathway may play important roles in the regulation of immune function under stress condition with CORT administration. The results of this study provide a reference for further studies of the molecular mechanisms of stress-influenced immune function.

Abstract Background The intracerebroventricular injection of visfatin increases feed intake. However, little is known about the molecular mechanism in chicks. This study was conducted to assess the effect of visfatin on the feeding behavior of chicks and the associated molecular mechanism. Results In response to the intraventricular injection of 40 ng and 400 ng visfatin, feed intake in chicks was significantly increased, and the concentrations of glucose, insulin, TG, HDL and LDL were significantly altered. Using RNA-seq, we identified DEGs in the chick hypothalamus at 60 min after injection with various doses of visfatin. In total, 325, 85 and 519 DEGs were identified in the treated chick hypothalamus in the LT vs C, HT vs C and LT vs HT comparisons, respectively. The changes in the expression profiles of DEGs, GO functional categories, KEGG pathways, and PPI networks by visfatin-mediated regulation of feed intake were analyzed. The DEGs were grouped into 8 clusters based on their expression patterns via K-mean clustering; there were 14 appetite-related DEGs enriched in the hormone activity GO term. The neuroactive ligand-receptor interaction pathway was the key pathway affected by visfatin. The PPI analysis of DEGs showed that POMC was a hub gene that interacted with the maximum number of nodes and ingestion-related pathways, including POMC, CRH, AgRP, NPY, TRH, VIP, NPYL, CGA and TSHB. Conclusion These common DEGs were enriched in the hormone activity GO term and the neuroactive ligand-receptor interaction pathway. Therefore, visfatin causes hyperphagia via the POMC/CRH and NPY/AgRP signaling pathways. These results provide valuable information about the molecular mechanisms of the regulation of food intake by visfatin.