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.

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.

Transcriptome identification and characterization of long non-coding RNAs in the ovary of hens at four stages.

Long non-coding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation. LncRNAs, which are defined as non-coding RNAs more than 200 bp in length, are involved in key biological processes, such as cell proliferation and differentiation, epigenetic regulation, and gene transcriptional translation. Recent studies have shown that lncRNAs also play major regulatory roles in the reproduction of mammals. However, knowledge of the roles of lncRNAs in the chicken ovary lacking. In this study, we performed RNA-seq analyses of ovarian tissue from Hy-Line brown laying hens at four physiological stages [15, 20, 30, and 68 weeks of age (W)]. We identified 657 lncRNA transcripts that were differentially expressed during ovarian development, the number of down-regulated lncRNAs was higher than the number of up-regulated lncRNAs during development. We predicted the cis and trans target genes of the DE lncRNAs and constructed a lncRNA-mRNA interaction network, which indicated that the DE genes (DEGs) and the target genes of the DE lncRNAs are mainly involved in signaling pathways associated with ovarian development, including oocyte meiosis, calcium signaling pathways, ECM-receptor interactions, and ribosome and focal adhesion. Overall, we found that twelve lncRNAs were strongly involved in ovarian development: LNC_013443, LNC_001029, LNC_005713, LNC_016762, ENSGALT00000101857, LNC_003913, LNC_013692, LNC_012219, LNC_004140, ENSGALT00000096941, LNC_009356, and ENSGALT00000098716. In summary, our study utilized RNA-seq analysis of hen ovaries to explore key lncRNAs involved in ovarian development and function. Furthermore, the comprehensive analysis identified the target genes of these lncRNAs providing a better understanding of the mechanisms underlying ovarian development in hens and a theoretical basis for further research.

Genome-wide association study identifies SNPs for growth performance and serum indicators in Valgus-varus deformity broilers (Gallus gallus) using ddGBS sequencing.

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.

Weighted gene co-expression network indicates that the DYNLL2 is an important regulator of chicken breast muscle development and is regulated by miR-148a-3p.

BACKGROUND: The characteristics of muscle fibers determine the growth and meat quality of poultry. In this study, we performed a weighted gene co-expression network analysis (WGCNA) on the muscle fiber characteristics and transcriptome profile of the breast muscle tissue of Gushi chicken at 6, 14, 22, and 30 weeks. RESULTS: A total of 27 coexpressed biological functional modules were identified, of which the midnight blue module had the strongest correlation with muscle fiber and diameter. In addition, 7 hub genes were found from the midnight blue module, including LC8 dynein light chain 2 (DYNLL2). Combined with miRNA transcriptome data, miR-148a-3p was found to be a potential target miRNA of DYNLL2. Experiments on chicken primary myoblasts (CPMs) demonstrated that miR-148a-3p promotes the expression of myosin heavy chain (MYHC) protein by targeting DYNLL2, proving that it can promote differentiation of myoblasts. CONCLUSIONS: This study proved that the hub gene DYNLL2 and its target miR-148-3p are important regulators in chicken myogenesis. These results provide novel insights for understanding the molecular regulation mechanisms related to the development of chicken breast muscle.

Genetic effect of an InDel in the promoter region of the NUDT15 and its effect on myoblast proliferation in chickens.

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.

The Chicken Pan-Genome Reveals Gene Content Variation and a Promoter Region Deletion in IGF2BP1 Affecting Body Size.

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional approximately 66.5-Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression levels are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based genome-wide association studies identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size quantitative trait locus located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Genome-wide identification evolution and expression of vestigial-like gene family in chicken.

Vestigial-like (Vgll) genes are widespread in vertebrates and play an important role in muscle development. In this study, we used bioinformatics methods to systematically identify the chicken VGLL family in the whole genome and investigated its evolutionary history and gene structure features. Tissue expression spectra combined with real-time PCR data were used to analyze the organizational expression pattern of the genes. Based on the maximum likelihood method, a phylogenetic tree of the VGLL family was constructed, and 94 VGLL genes were identified in 24 breeds, among which four VGLL family genes were identified in the chicken genome. Ten motifs were detected in the VGLL genes, and the analysis of introns combined with gene structure revealed that the family was conserved during evolution. Tissue expression analysis suggested that the expression profiles of the VGLL family genes in 16 tissues differed between LU Shi and AA broilers. In addition, a single gene (VGLL2) showed increased expression in chickens at embryonic days 10-16 and was involved in the growth and development of skeletal muscle in chickens in the embryonic stage. In summary, VGLL genes are involved in chicken muscle growth and development, which provides useful information for subsequent functional studies of VGLL genes.

LncRNAs and their regulatory networks in breast muscle tissue of Chinese Gushi chickens during late postnatal development.

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.

Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population.

Chicken growth traits are economically important, but the relevant genetic mechanisms have not yet been elucidated. Herein, we performed a genome-wide association study to identify the variants associated with growth traits. In total, 860 chickens from a Gushi-Anka F2 resource population were phenotyped for 68 growth and carcass traits, and 768 samples were genotyped based on the genotyping-by-sequencing (GBS) method. Finally, 734 chickens and 321,314 SNPs remained after quality control and removal of the sex chromosomes, and these data were used to carry out a GWAS analysis. A total of 470 significant single-nucleotide polymorphisms (SNPs) for 43 of the 68 traits were detected and mapped on chromosomes (Chr) 1-6, -9, -10, -16, -18, -23, and -27. Of these, the significant SNPs in Chr1, -4, and -27 were found to be associated with more than 10 traits. Multiple traits shared significant SNPs, indicating that the same mutation in the region might have a large effect on multiple growth or carcass traits. Haplotype analysis revealed that SNPs within the candidate region of Chr1 presented a mosaic pattern. The significant SNPs and pathway enrichment analysis revealed that the MLNR, MED4, CAB39L, LDB2, and IGF2BP1 genes could be putative candidate genes for growth and carcass traits. The findings of this study improve our understanding of the genetic mechanisms regulating chicken growth and carcass traits and provide a theoretical basis for chicken breeding programs.

A novel model to predict mental distress among medical graduate students in China.

BACKGROUND: Poor mental health was reported among medical graduate students in some studies. Identification of risk factors for predicting the mental health is capable of reducing psychological distress among medical graduate students. Therefore, the aim of the study was to identify potential risk factors relating to mental health and further create a novel prediction model to calculate the risk of mental distress among medical graduate students. METHODS: This study collected and analyzed 1079 medical graduate students via an online questionnaire. Included participants were randomly classified into a training group and a validation group. A model was developed in the training group and validation of the model was performed in the validation group. The predictive performance of the model was assessed using the discrimination and calibration. RESULTS: One thousand and fifteen participants were enrolled and then randomly divided into the training group (n = 508) and the validation group (n = 507). The prevalence of severe mental distress was 14.96% in the training group, and 16.77% in the validation group. The model was developed using the six variables, including the year of study, type of student, daily research time, monthly income, scientific learning style, and feeling of time stress. The area under the receiver operating characteristic curve (AUROC) and calibration slope for the model were 0.70 and 0.90 (95% CI: 0.65 ~ 1.15) in the training group, respectively, and 0.66 and 0.80 (95% CI, 0.51 ~ 1.09) in the validation group, respectively. CONCLUSIONS: The study identified six risk factors for predicting anxiety and depression and successfully created a prediction model. The model may be a useful tool that can identify the mental status among medical graduate students. TRIAL REGISTRATION: No. ChiCTR2000039574 , prospectively registered on 1 November 2020.

Breeding history and candidate genes responsible for black skin of Xichuan black-bone chicken.

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.

Detection of CNV in the SH3RF2 gene and its effects on growth and carcass traits in chickens.

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.

α7 Nicotinic Acetylcholine Receptor Mediates the Neuroprotection of Remote Ischemic Postconditioning in a Rat Model of Asphyxial Cardiac Arrest.

BACKGROUND: Remote ischemic postconditioning (RIPost) has been shown to reduce the ischemia-reperfusion injury of the heart and brain. However, the protection mechanisms have not yet been fully elucidated. We have observed that RIPost could alleviate the brain injury after cardiac arrest (CA). The aim of this study was to explore whether α7 nicotinic acetylcholine receptor (α7nAChR) mediates the neuroprotection of RIPost in a rat model of asphyxial CA. MATERIALS AND METHODS: Asphyxial CA model was induced by occlusion of the tracheal tube for 8 min and resuscitated later. RIPost produced by three cycles of 15-min occlusion and 15-min release of the right hind limb by a tourniquet was performed respectively at the moment and the third hour after restoration of spontaneous circulation. The α7nAChR agonist PHA-543613 and the antagonist methyllycaconitine (MLA) were used to investigate the role of α7nAChR in mediating neuroprotective effects. RESULTS: Results showed that α7nAChR was decreased in hippocampus and cortex after resuscitation, whereas RIPost could attenuate the reduction. The use of PHA-543613 provided neuroprotective effects against cerebral injury after CA. Furthermore, RIPost decreased the levels of neuron-specific enolase, inflammatory mediators, the number of apoptotic cells, and phosphorylation of nuclear factor-κB while increased the phosphorylation of signal transducer and activator of transcription-3. However, the above effects of RIPost were attenuated by α7nAChR antagonist methyllycaconitine. CONCLUSIONS: Neuroprotection of RIPost was related with the activation of α7nAChR, which could suppress nuclear factor-κB and activate signal transducer and activator of transcription-3 in a rat asphyxial CA model.

Gut microbiota profiles of commercial laying hens infected with tumorigenic viruses.

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.

Target gene identification and functional characterization of miR-1704 in chicken.

MiRNAs are small non-coding RNAs that negatively regulate gene expression at the post-transcriptional level. SNPs in miRNA genes may lead to phenotypic variation by altering miRNA expression and their targets. In this study, miR-1704 expression profiles in nine tissues at 1 d, 6 weeks and 16 weeks old Gushi chickens were detected. MiR-1704 was widely expressed in the detection of tissues. The expression in 1 d and 6 weeks old was low abundance, while its expression at 16 weeks was very high. An rs14668705 (C > G) SNP was detected within the pre-miR-1704 in an F2 resource population of Gushi chicken crossed with Anka broiler. Bioinformatic analysis indicated that the C > G mutation could introduce a base-pair mismatch and cause the change of free energy. Experiments further revealed that the rs14668705 in precursor miR-1704 could significantly affect mature miR-1704 biogenesis and was significantly associated with body weight at the age of 0, 6, 8, 10, and 12 weeks, shank circumference at 4, 8, and 12 weeks, carcass weight, and semi-evisceration weight (p < 0.05). Insulin receptor 2 (IRS2) gene, one of the potential targets of miR-1704 was identified and further confirmed. These data suggested that miR-1704 targeted IRS2 and have an effect on body weight in chicken.

Estrogen Abolishes the Repression Role of gga-miR-221-5p Targeting ELOVL6 and SQLE to Promote Lipid Synthesis in Chicken Liver.

Few studies have been conducted regarding the biological function and regulation role of gga-miR-221-5p in the liver. We compared the conservation of miR-221-5p among species and investigated the expression pattern of gga-miR-221-5p, validating the direct target genes of gga-miR-221-5p by dual luciferase reporter assay, the biological function of gga-miR-221-5p in the liver was studied by gga-miR-221-5p overexpression and inhibition. Furthermore, we explored the regulation of gga-miR-221-5p and its target genes by treatment with estrogen and estrogen antagonists in vivo and in vitro. The results showed that miR-221-5p was highly conserved among species, expressed in all tested tissues and significantly downregulated in peak-laying hen liver compared to pre-laying hen liver. Gga-miR-221-5p could directly target the expression of elongase of very long chain fatty acids 6 (ELOVL6) and squalene epoxidase (SQLE) genes to affect triglyceride and total cholesterol content in the liver. 17ß-estradiol could significantly inhibit the expression of gga-miR-221-5p but promote the expression of ELOVL6 and SQLE genes. In conclusion, the highly conservative gga-miR-221-5p could directly target ELOVL6 and SQLE mRNAs to affect the level of intracellular triglyceride and total cholesterol. Meanwhile, 17ß-estradiol could repress the expression of gga-miR-221-5p but increase the expression of ELOVL6 and SQLE, therefore promoting the synthesis of intracellular triglyceride and cholesterol levels in the liver of egg-laying chicken.

MicroRNAs and their regulatory networks in Chinese Gushi chicken abdominal adipose tissue during postnatal late development.

BACKGROUND: Abdominal fat is the major adipose tissue in chickens. The growth status of abdominal fat during postnatal late development ultimately affects meat yield and quality in chickens. MicroRNAs (miRNAs) are endogenous small noncoding RNAs that regulate gene expression at the post-transcriptional level. Studies have shown that miRNAs play an important role in the biological processes involved in adipose tissue development. However, few studies have investigated miRNA expression profiles and their interaction networks associated with the postnatal late development of abdominal adipose tissue in chickens. RESULTS: We constructed four small RNA libraries from abdominal adipose tissue obtained from Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 507 known miRNAs and 53 novel miRNAs were identified based on the four small RNA libraries. Fifty-one significant differentially expressed (SDE) miRNAs were identified from six combinations by comparative analysis, and the expression patterns of these SDE miRNAs were divided into six subclusters by cluster analysis. Gene ontology enrichment analysis showed that the SDE miRNAs were primarily involved in the regulation of fat cell differentiation, regulation of lipid metabolism, regulation of fatty acid metabolism, and unsaturated fatty acid metabolism in the lipid metabolism- or deposition-related biological process categories. In addition, we constructed differentially expressed miRNA-mRNA interaction networks related to abdominal adipose development. The results showed that miRNA families, such as mir-30, mir-34, mir-199, mir-8, and mir-146, may have key roles in lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens. CONCLUSIONS: This study determined the dynamic miRNA transcriptome and characterized the miRNA-mRNA interaction networks in Gushi chicken abdominal adipose tissue for the first time. The results expanded the number of known miRNAs in abdominal adipose tissue and provide novel insights and a valuable resource to elucidate post-transcriptional regulation mechanisms during postnatal late development of abdominal adipose tissue in chicken.

Identification of differentially expressed genes and pathways between intramuscular and abdominal fat-derived preadipocyte differentiation of chickens in vitro.

BACKGROUND: The distribution and deposition of fat tissue in different parts of the body are the key factors affecting the carcass quality and meat flavour of chickens. Intramuscular fat (IMF) content is an important factor associated with meat quality, while abdominal fat (AbF) is regarded as one of the main factors affecting poultry slaughter efficiency. To investigate the differentially expressed genes (DEGs) and molecular regulatory mechanisms related to adipogenic differentiation between IMF- and AbF-derived preadipocytes, we analysed the mRNA expression profiles in preadipocytes (0d, Pre-) and adipocytes (10d, Ad-) from IMF and AbF of Gushi chickens. RESULTS: AbF-derived preadipocytes exhibited a higher adipogenic differentiation ability (96.4% + 0.6) than IMF-derived preadipocytes (86.0% + 0.4) (p < 0.01). By Ribo-Zero RNA sequencing, we obtained 4403 (2055 upregulated and 2348 downregulated) and 4693 (2797 upregulated and 1896 downregulated) DEGs between preadipocytes and adipocytes in the IMF and Ad groups, respectively. For IMF-derived preadipocyte differentiation, pathways related to the PPAR signalling pathway, ECM-receptor interaction and focal adhesion pathway were significantly enriched. For AbF-derived preadipocyte differentiation, the steroid biosynthesis pathways, calcium signaling pathway and ECM-receptor interaction pathway were significantly enriched. A large number of DEGs related to lipid metabolism, fatty acid metabolism and preadipocyte differentiation, such as PPARG, ACSBG2, FABP4, FASN, APOA1 and INSIG1, were identified in our study. CONCLUSION: This study revealed large transcriptomic differences between IMF- and AbF-derived preadipocyte differentiation. A large number of DEGs and transcription factors that were closely related to fatty acid metabolism, lipid metabolism and preadipocyte differentiation were identified in the present study. Additionally, the microenvironment of IMF- and AbF-derived preadipocyte may play a significant role in adipogenic differentiation. This study provides valuable evidence to understand the molecular mechanisms underlying adipogenesis and fat deposition in chickens.

Weighted gene coexpression network analysis identifies specific transcriptional modules and hub genes related to intramuscular fat traits in chicken breast muscle.

Intramuscular fat (IMF) traits are important factors that influence meat quality. However, the molecular regulatory mechanisms that underlie this trait in chickens are still poorly understood at the gene coexpression level. Here, we performed a weighted gene coexpression network analysis between IMF traits and transcriptome profile in breast muscle in the Chinese domestic Gushi chicken breed at 6, 14, 22, and 30 weeks. A total of 26 coexpressed gene modules were identified. Six modules, which included the dark gray, purple, cyan, pink, light cyan, and blue modules, showed a significant positive correlation (P < 0.05) with IMF traits. The strongest correlation was observed between the dark gray module and IMF content (r = 0.85; P = 4e-04) and between the blue module and different fatty acid content (r = 0.87~0.91; P = 5e-05~2e-04). Enrichment analysis showed that the enrichment of biological processes, such as fatty acid metabolic process, fat cell differentiation, acylglycerol metabolic process, and glycerolipid metabolism were significantly different in the six modules. In addition, the 32, 24, 4, 7, 6, and 25 hub genes were identified from the blue, pink, light cyan, cyan, dark gray, and purple modules, respectively. These hub genes are involved in multiple links to fatty acid metabolism, phospholipid metabolism, cholesterol metabolism, diverse cellular behaviors, and cell events. These results provide novel insights into the molecular regulatory mechanisms for IMF-related traits in chicken and may also help to uncover the formation mechanism for excellent meat quality traits in local breeds of Chinese chicken.