Search

Abstract Breast muscle growth rate and intramuscular fat (IMF) content show apparent differences between fast-growing broilers and slow-growing indigenous chickens. However, the underlying genetic basis of these phenotypic characteristics remains elusive. In this study, we investigate the dynamic alterations of three-dimensional genome architecture and chromatin accessibility in breast muscle across four key developmental stages from embryo to starter chick in Arbor Acres (AA) broilers and Yufen (YF) indigenous chickens. The limited breed-specifically up-regulated genes (Bup-DEGs) are embedded in breed-specific A compartment, while a majority of the Bup-DEGs involving myogenesis and adipogenesis are regulated by the breed-specific TAD reprogramming. Chromatin loops allow distal accessible regions to interact with myogenic genes, and those loops share an extremely low similarity between chicken with different growth rate. Moreover, AA-specific loop interactions promote the expression of 40 Bup-DEGs, such as IGF1, which contributes to myofiber hypertrophy. YF-specific loop interactions or distal accessible regions lead to increased expression of 5 Bup-DEGs, including PIGO, PEMT, DHCR7, TMEM38B, and DHDH, which contribute to IMF deposition. These results help elucidate the regulation of breast muscle growth and IMF deposition in chickens.

Abstract Background The breeder rooster has played a pivotal role in poultry production by providing high-quality semen. Typically, fertility peaks between 30 and 40 weeks of age and then declines rapidly from 45 to 55 weeks of age. Research into improving fertility in aging roosters is essential to extend their productive life. While progress has been made, enhancing fertility in aging roosters remains a significant challenge. Methods To identify the genes related to promoting sperm remodeling in aged Houdan roosters, we combined changes in testis and semen quality with transcriptome sequencing (RNA-seq) to analyze the synchrony of semen quality and testis development. In this study, 350-day-old Houdan breeder roosters were selected for RNA-seq analysis in testis tissues from induced molting roosters (D group) and non-induced molting roosters (47DG group). All analyses of differentially expressed genes (DEGs) and functional enrichment were performed. Finally, we selected six DEGs to verify the accuracy of the sequencing by qPCR. Results Compared with the 47DG group, sperm motility (P

ABSTRACT: Valgus-varus deformity (VVD) is a common long bone deformity in broilers. Imbalance in cartilage homeostasis is the main feature of leg disease. Exosomes act as an important intercellular communication vector that regulates chondrogenesis by encapsulating specific nucleic acids and proteins. However, the exact mechanism of how plasma exosomal miRNAs regulate cartilage homeostasis in VVD broilers remains unclear. This study first demonstrated the structural disorder, growth retardation, and reduced proliferative capacity of VVD cartilage in vitro and in vivo. Subsequently, VVD and Normal broiler plasma exosomes were collected for miRNA sequencing. Cartilage-specific miR-455-5p was extraordinarily emphasized by performing bioinformatics analysis on differential miRNA target genes and further validated by tissue expression profiling. PKH67 fluorescently labeled plasma exosomes were shown to be taken up by chondrocytes, deliver miR-455-5p, inhibit chondrocyte proliferation, and disrupt their homeostasis, and these effects could be inhibited by the miR-inhibitors. Mechanistically, MiR-455-5p targets Ribosomal Protein S6 Kinase B1 (RPS6KB1) to inhibit RPS6 phosphorylation and reduce the synthesis of key proteins for cartilage proliferation, which in turn inhibits cartilage proliferation and disrupts its homeostasis. In conclusion, the present study identified abnormalities in VVD cartilage tissue and clarified the specific mechanism by which plasma exosome-derived miR-455-5p regulates cartilage homeostasis.

Conservation programs require rigorous evaluation to ensure the preservation of genetic diversity and viability of conservation populations. In this study, we conducted a comparative analysis of two indigenous Chinese chicken breeds, Gushi and Xichuan black-bone, using whole-genome SNPs to understand their genetic diversity, track changes over time and population structure. The breeds were divided into five conservation populations (GS1, 2010, ex-situ; GS2, 2019, ex-situ; GS3, 2019, in-situ; XB1, 2010, in-situ; and XB2, 2019, in-situ) based on conservation methods and generations. The genetic diversity indices of three conservation populations of Gushi chicken showed consistent trends, with the GS3 population under in-situ strategy having the highest diversity and GS2 under ex-situ strategy having the lowest. The degree of inbreeding of GS2 was higher than that of GS1 and GS3. Conserved populations of Xichuan black-bone chicken showed no obvious changes in genetic diversity between XB1 and XB2. In terms of population structure, the GS3 population were stratified relative to GS1 and GS2. According to the conservation priority, GS3 had the highest contribution to the total gene and allelic diversity in GS breed, whereas the contribution of XB1 and XB2 were similar. We also observed that the genetic diversity of GS2 was lower than GS3, which were from the same generation but under different conservation programs (in-situ and ex-situ). While XB1 and XB2 had similar levels of genetic diversity. Overall, our findings suggested that the conservation programs performed in ex-situ could slow down the occurrence of inbreeding events, but could not entirely prevent the loss of genetic diversity when the conserved population size was small, while in-situ conservation populations with large population size could maintain a relative high level of genetic diversity.

ABSTRACT: Goose astrovirus (GAstV) is a newly discovered astrovirus. GAstV causes gout and death in 4- to -16-day-old goslings. For the past few years, fatal gout, the cardinal clinical symptom of gosling infected with GAstV, has been spreading rapidly in some goose Chinese farms, which caused continuous economic losses to the goose breeding industry in China. Currently, several underlying mechanisms involved in viral replication, inflammatory reaction, virions release, and viral pathogenesis of GAstV remain to be elucidated. In this study, we explored the mechanisms of GAstV-host interactions, the transcriptome and proteome profiles of GAstV-infected LMH cells were sequenced by RNA-seq and data-independent acquisition (DIA) techniques, respectively, and followed using an integrative analysis. Compared with uninfected LMH cells, a total of 322 differentially expressed genes (DEG) (195 up-regulated, 127 down-regulated) and 36 differentially expressed proteins (DEP) (31 up-regulated, 5 down-regulated) were detected. Nine DEGs were randomly selected for further validation by quantitative real-time polymerase chain reaction (qPCR). Through GO and KEGG enrichment analysis, DEG and DEP were significantly enriched in several important cellular signaling pathways, including MAPK, PI3K-Akt, cAMP, chemokine, calcium, phospholipase D, Ras, TNF, IL-17, Rap1, NF-kappa B signaling pathways, indicating that GAstV affects cell growth and immune signaling. This study provided an overview of changes in transcriptome and proteome profiles of GAstV-infected LMH cells, therefore, providing a crucial basis to further explore the mechanisms of GAstV-host interactions.

ABSTRACT: The high-yielding Houdan chicken (GGF) is characterized by high egg production and disease resistance. This study conducted whole genome resequencing of the GGF population and compared it to data from other breeds. Genetic diversity analysis revealed higher observed heterozygosity (Ho), Polymorphism information content (PIC), number of runs of homozygosity (ROH), and inbreeding coefficient (FROH) in GGF. Linkage disequilibrium (LD) decay was slowest in GGF, indicating intensive inbreeding and strong selection. These findings suggest a need for appropriate strategies to enhance genetic diversity conservation in this breed. Population structure analysis demonstrated that GGF was genetically distinct from both the red jungle fowl (RJF) and Chinese indigenous chicken (CIC) populations, highlighting GGF as a unique genetic resource warranting intensive protection and utilization. Selective sweep analysis identified genes under selection in GGF, primarily enriched in signaling pathways related to oocyte meiosis and progesterone-mediated oocyte maturation. Key candidate genes included: CCNE1, SKP1, CDC20, CDK2, ADCY8, RPS6KA6, PPP3CB, PDE3B, HSP90AB1, and AKT3. These findings provide a theoretical foundation for their potential application in poultry breeding. Additionally, this study combined bioinformatics analysis with PCR amplification and Sanger sequencing to identify 4 SNPs that can serve as a molecular identity card (ID) for GGF: SNP1 (Chr2: 136130976), SNP3 (Chr4:11705164), SNP4 (Chr4: 63255588), and SNP5 (Chr24: 3271008). This study provides a scientific basis for effective management and conservation of GGF genetic resources, and establishes a simple, economical, and accurate set of molecular IDs to combat the proliferation of inferior breeds and protect genetic resources.

ABSTRACT: 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.

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 Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan‐genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan‐genome together with high‐depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE‐derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan‐genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock breeding.

ABSTRACT: The “Yufen 1” H line chicken (YF) has excellent characteristics including early sexual maturity and high egg production, and the conservation of its genetic diversity is the core of the breeding activity. To overcome misrepresented breeds and protect the integrity of the germplasm genetic resources, it is important to develop accurate and convenient methods to identify YF. In this study, whole genome resequencing was performed on the YF population, and bioinformatics analysis was conducted by combining the data from different breeds. Linkage disequilibrium (LD) analysis revealed that YF had the slowest LD-decay rate, suggesting strong natural and artificial selection in its history. Through selective sweep analysis, 1,126 selected regions in YF were identified, which contained 163,661 single nucleotide polymorphisms (SNPs). In particular, 5 specific SNPs (SNP1: Chr2:45509616, SNP2: Chr2:45510792, SNP3: Chr9:13788193, SNP4: Chr9:13795646, SNP5: Chr9:13798154) were found exclusively in the YF population. Subsequently, PCR amplification and Sanger sequencing confirmed the presence of these 5 SNPs in YF. Finally, 4 SNPs (SNP1, SNP2, SNP4, SNP5) were screened and verified using the Kompetitive Allele Specific PCR (KASP) typing technique. These SNPs can be used as specific molecular identity cards (IDs) for YF authentication. The present study is of great significance to ensure sustainable conservation and promotion of YF germplasm resources.

Abstract Background and Aim Primary biliary cholangitis (PBC) is a chronic autoimmune cholangiopathy, characterized by the presence of some autoantibodies in the serum. This study aimed to evaluate the clinical significance of AMA‐M2, anti‐gp210 and anti‐sp100 antibody levels detected by multiplex bead‐based flow fluorescent immunoassay (MBFFI) in PBC. Methods This study cohort included 238 PBC patients, 81 autoimmune hepatitis (AIH) patients, 62 systemic lupus erythematosus (SLE) patients, and 118 healthy controls. Serum AMA‐M2, anti‐gp210 and anti‐sp100 antibody were detected by MBFFI and immunoblotting assay (IBT). The relationship between three antibody levels and cirrhosis, liver function, cholestasis markers and therapeutic effect to ursodesoxycholic acid (UDCA) was evaluated in PBC. Results MBFFI were presented good coincidence rate (87.39%–95.38%) with IBT. The level of AMA‐M2, anti‐gp210 and anti‐sp100 antibodies in PBC patients were higher than other disease group and healthy controls (p ＜ .01). When compared with the healthy controls group, the AUC of AMA‐M2, anti‐gp210 and anti‐sp100 antibodies were 0.9245, 0.7619, and 0.6789, respectively. In addition, gp210 antibody levels have diagnostic value in patients with liver cirrhosis (AUC: 0.7567). We found that when combine detect these three antibodies, the sensitivity was higher than individually detection. High level of serum anti‐gp210 antibody could be related to worse liver function and more severe cholestasis in PBC patients. Moreover, serum antibody levels may decrease or remained flat in patients who responded well to UDCA. Conclusion The detection of AMA‐M2, anti‐gp210 and anti‐sp100 antibody levels by MBFFI showed good performance in the diagnosis of PBC. Serum anti‐gp210 antibody level is related to cirrhosis, poor liver function and severe cholestasis in PBC.

Marek's disease virus (MDV), a naturally oncogenic, highly contagious alpha herpesvirus, induces a T cell lymphoma in chickens that causes severe economic loss. Marek's disease (MD) outcome in an individual is attributed to genetic and environmental factors. Further investigation of the host-virus interaction mechanisms that impact MD resistance is needed to achieve greater MD control. This study analyzed genome-wide DNA methylation patterns in 2 highly inbred parental lines 63 and 72 and 5 recombinant congenic strains (RCS) C, L, M, N, and X strains from those parents. Lines 63 and 72, are MD resistant and susceptible, respectively, whereas the RCS have different combinations of 87.5% Line 63 and 12.5% Line 72. Our DNA methylation cluster showed a strong association with MD incidence. Differentially methylated regions (DMRs) between the parental lines and the 5 RCS were captured. MD-resistant and MD-susceptible markers of DNA methylation were identified as transgenerational epigenetic inheritable. In addition, the growth of v-src DNA tumors and antibody response against sheep red blood cells differed among the 2 parental lines and the RCS. Overall, our results provide very solid evidence that DNA methylation patterns are transgenerational epigenetic inheritance (TEI) in chickens and also play a vital role in MD tumorigenesis and other immune responses; the specific methylated regions may be important modulators of general immunity.

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.

Fasting-induced molting (FIM) is a common method used to improve the laying performance of aged laying hens. Nevertheless, this approach may impose various stresses on chickens, such as disruptions in intestinal flora and inflammation issues within the intestines. However, the impact of an imbalance in intestinal flora on intestinal health during the FIM process remains elusive. Therefore, intestinal injury, the microbiome, and the metabolome were analyzed individually and integrated to elucidate the impact of the intestinal flora on intestinal health during the FIM process. The findings indicated that fasting resulted in a notable reduction in villus height and villus/crypt ratio, coupled with elevated levels of intestinal inflammation and permeability. During the fasting period, microbiota compositions changed. The abundance of Escherichia_Shigella increased, while the abundance of Ruminococcaceae_UCG-013 and Lactobacillus decreased. Escherichia_Shigella was positively correlated with Citrinin and Sterobilin, which lead to intestinal inflammation. Ruminococcaceae_UCG-013 and Lactobacillus exhibited positive correlations with Lanthionine and reduced Glutathione, thereby reducing intestinal inflammation. This study screened the intestinal probiotics, Ruminococcaceae UCG-013 and Lactobacillus, that influence gut health during the fasting period, providing an experimental basis for improving gut microbiota and reducing intestinal inflammation during the FIM process.

ABSTRACT: The conservation of genetic resources is becoming increasingly important for the sustainable development of the poultry industry. In the present study, we systematically analyzed the population structure, conservation priority, runs of homozygosity (ROH) of chicken breeds globally, and proposed rational conservation strategies. We used a 600K Affymetrix Axiom HD genotyping SNP array dataset of 2,429 chickens from 134 populations. The chickens were divided into 5 groups based on their country of origin and sampling location: Asian chickens (AS-LOC), African chickens (AF), European local chickens (EU-LOC), Asian breeds sampled in Germany (AS-DE), and European breeds sampled in Germany (EU-DE). The results indicated that the population structure was consistent with the actual geographical distribution of the populations. AS-LOC had the highest positive contribution to the total gene (HT, 1.00%,) and allelic diversity (AT, 0.0014%), the lowest inbreeding degree and the fastest linkage disequilibrium (LD) decay rate; the lowest contribution are derived by European ex situ chicken breeds (EU-DE:HT = −0.072%, AT = −0.0014%), which showed the highest inbreeding and slowest LD decay. Breeds farmed in ex situ (AS-DE, EU-DE) conditions exhibited reduced genetic diversity and increased inbreeding due to small population size. Given limited funds, it is a better choice for government to conserve the breeds with the highest contribution to genetic diversity in each group. Therefore, we evaluated the contribution of each breed to genetic and allelic diversity in 5 groups. Among each group, KUR(AF), BANG(AS-LOC), ALxx(EU-LOC), BHwsch(AS-DE), and ARw(EU-DE) had the highest contribution to gene diversity in the order of the above grouping. Similarly, according to the allelic diversity standard (in the same order), ZIMxx, PIxx, ALxx, SHsch, and ARsch had the highest contribution. After analyzing ROH, we found a total of 144,708 fragments and 27 islands. The gene and genome regions identified by the ROH islands and QTLs indicate that chicken breeds have potential for adaptation to different production systems. Based on these findings, it is recommended to prioritize the conservation of breeds with the highest genetic diversity in each group, while paying more attention to the conservation of Asian and African breeds. Furthermore, providing a valuable reference for the conservation and utilization of chicken.

ABSTRACT: 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.

ABSTRACT: 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.

ABSTRACT: 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.

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

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.

The goal of this study was to investigate the soil particle-size distribution (PSD) and pore structure characteristics in the Kubuqi Desert in order to provide basic data for gaining insights into the soil-modifying properties of the local vegetation. Based on laser diffraction analysis, we measured the soil PSD and calculated the single and multi-fractal dimensions of the soils under typical artificial forest and natural grassland vegetation. The diameters, volumes and surface areas of the soil pores were determined via nitrogen adsorption experiments. The relationships between the fractal characteristics and pore structure parameters of the soil were investigated via redundancy analysis and stepwise regression. The clay content was less than 2.0%, and the sand content was greater than 75.0%. There was variability in the PSD and fractal dimension, which was mainly observed for the 0–20 cm, 20–80 cm, and 80–100 cm soil layers. In the 0–80 cm, the fine particle content, single fractal dimension (D), entropy dimension (D1), ratio of D1 to the capacity dimension (D1/D0), and correlation dimension (D2) were smaller for the forest soils than for the grassland soils. The fine particle (clay and silt) content of the 80–100 cm soil layer was approximately 37.8% higher for the Salix psammophila-8a than for the Salix psammophila-3a, and it was approximately 161.4% higher for the Populus popular-35a than for the Pinus sylvestris var. mongolica-8a. The silt content, D1, D1/D0, and D2 were significantly positively correlated with the specific surface area (SSA), total pore volume (TPV), and average pore diameter (APD) of the nanopores (p < 0.05, 0.01, or 0.001), and they were significantly negatively correlated with the percentage of the micropore volume (PMV) (p < 0.05 or 0.01). In the Kubuqi Desert, the fine particle content and fractal dimensions of the soil layer below the root zone of shrub and arbor vegetation increased with increasing stand age, but the trend was reversed in the shallower soil layers. The variability of the soil PSD characteristics was strongly correlated with the variability of the nanopore parameters on the microscopic scale, suggesting that the total pore volume, average pore diameter, and percentage of the micropore volume may be potential indicators of the soil structure and quality.

ABSTRACT: The reproductive performance of chicken breeders has significant economic importance in the poultry industry, and sperm motility is an indicator of reproductive performance. This study performed RNA-seq of the testes of Gushi chicken roosters with high and low sperm motility and identified differentially expressed RNAs involved in sperm motility. RNA-seq analysis showed that 73 and 67 differentially expressed mRNAs were up- and downregulated, and 47 and 56 differentially expressed long non-coding RNAs were up- and downregulated, respectively. The genes related to sperm motility and spermatogenesis included KIFC1, KCNK2, and REC8. Functional enrichment analysis revealed that the pathways related to sperm motility included oxidative phosphorylation and glycine, serine, and threonine metabolism. In addition, the MSTRG.15920.1-REC8-MSTRG.11860.2-VWC2 pathway may regulate sperm motility. This study helped elucidate the molecular genetic mechanism of sperm motility in chicken.

Abstract Introduction Previous studies of anti‐dsDNA, nucleosome (Nucl), histone (His), and C1q antibodies have revealed their clinical value in systemic lupus erythematosus (SLE). However, the correlation between four autoantibodies and SLE activity, lupus nephritis (LN) remains controversial, and data are insufficient on longitudinal monitoring. This study aimed at evaluating the value of these autoantibodies in active LN, and their performance on cross‐sectional evaluating and longitudinal monitoring of SLE disease activity. Methods Serum levels of four autoantibodies in 114 SLE patients, 219 other autoimmune disease patients (OAD), and 59 healthy controls were assayed by a quantitative immunoassay. Sera of 38 inpatients were obtained again after treatment. Results We found that serum levels of four autoantibodies were significantly higher in SLE than OAD patients (p .05). Conclusions Anti‐C1q, anti‐dsDNA, anti‐Nucl, and anti‐His perform well in diagnosing active LN and monitoring SLE disease activity. They could be indicators of active LN and SLE disease activity.

Coprophagy is an instinctive behavior in rabbit with important effects on growth and reproductive performance. The underlying mechanism of this effect in rabbit is unknown. Here, we used Elizabeth circle as a coprophagy preventing model in female rabbits and assess feed intake, growth, and reproductive performance. We found that preventing coprophagy did not affect feed intake but decreased body weight and weight of several organs and tissues and resulted in complete reproductive failure during the late pregnancy period, accompanied by reduced levels of plasma progesterone. RNA-seq analysis of rabbit ovarian tissues revealed that preventing coprophagy affected significantly 241 genes (DEGs), with the large majority being downregulated. Bioinformatic analyses revealed that those DEGs are mostly involved in apoptosis, immune response, and metabolic pathways. Among DEGs, the lysosomal cysteine protease cathepsin B (CTSB) was significantly downregulated in the coprophagy prevention group. Further studies using siRNA and adenovirus overexpression systems revealed that CTSB promotes the proliferation of rabbit granulosa cells (GCS) and prevents apoptosis. Measurement of transcripts coding for proteins related to apoptosis revealed a minor transcriptomic effect of CTSB, indicating that its effect is likely post-transcriptional. Overexpression of CTSB increased secretion of progesterone and estradiol, partly via upregulation of CYP19A1 while inhibition of CTSB decreased progesterone secretion partly via downregulation of the StAR gene. In conclusion, our study demonstrated the detrimental effect on reproduction by preventing coprophagy with a main role for this response played by CTSB on the granulosa cells of the ovary.

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.

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.

Background: Patterns of liver energy metabolism significantly differ from birth to adult in cattle undergoing change of rumen rumination. However, the genes involve in hepatic energy metabolism during bovine development and how regulate are still unclear.Methods: In this study, 0-day-old newborn calves (0W) and 9-week-old weaned calves (9W) were used to investigate differences in liver glucose metabolism at these stages of calf development. We did this primarily through the quantitation of energy metabolism indicators, then sequencing the liver transcriptome for each group of claves.Results: The transcriptome results showed 979 differentially expressed genes (DEGs), enriched in animal organ development, catabolic process, transmembrane transport. SLC16A1 involved in that and was locked to investigate. We explored the effects of SLC16A1 on glucose and lactate flux in vitro. We identified and verified its target, miR-22-3p, through bioinformatics and luciferase reporter assays. Moreover, this study found that miR-22-3p decreased cell activity by negatively regulating the SLC16A1. Importantly, our result showed the insulin-induced SLC16A1 mRNA expression decreased, regulated by promoter activity rather than miR-22-3p.Conclusions: Our study illustrates the role of SLC16A1 in the liver mediated metabolism of developing calves. These data enrich our knowledge of the regulatory mechanisms of liver mediated glucose metabolism in developing cattle.

Moult is a normal physiological phenomenon in poultry. Induced molting (IM) is the most widely used and economical molting technique. By inducing moult, the laying hens can grow new feathers during the next laying cycle and improve laying performance. However, the lack of energy supply has a huge impact on both the liver and intestines and acts on the intestines and liver through the “gut-liver axis”. More importantly, lipid metabolism in the liver is closely related to the laying performance of laying hens. Therefore, in this study, cecal metabolites and liver transcriptome data during IM of laying hens at the late stage of laying (stop feeding method) were analyzed together to reveal the regulatory mechanism of “gut-liver axis” affecting the laying performance of laying hens from the perspective of lipid metabolism. Transcriptome analysis revealed that 4,796 genes were obtained, among which 2,784 genes had significant differences (p < 0.05). Forty-nine genes were associated with lipid metabolism, and five core genes (AGPAT2, SGPL1, SPTLC1, PISD, and CYP51A1) were identified by WGCNA. Most of these differential genes are enriched in steroid biosynthesis, cholesterol metabolism, drug metabolism—cytochrome P450, synthesis and degradation of ketone bodies, PPAR signaling pathway, and bile secretion. A total of 96 differential metabolites were obtained by correlating them with metabolome data. Induced moult affects laying performance by regulating genes related to lipid metabolism, and the cecal metabolites associated with these genes are likely to regulate the expression of these genes through the “enterohepatic circulation”. This experiment enriched the theoretical basis of induced moult and provided the basis for prolonging the feeding cycle of laying hens.

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 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: Several recent studies investigated the role of the miR-29 family in muscle development. However, only a few studies focused on chicken skeletal muscle. In the present study, cell cycle, 5-ethynyl-2’-deoxyuridine (EdU), cell counting kit-8 (CCK-8), and other assays indicated that miR-29b-1-5p can inhibit the proliferation of chicken primary myoblasts (CPMs); the western blot assay and immunofluorescence detection of MYHC demonstrated that miR-29b-1-5p can promote the differentiation of myoblasts. The functional enrichment analysis revealed that the target genes of miR-29b-1-5p may be involved in muscle tissue development, muscle organ development, and striated muscle tissue development, which are biological processes related to muscle development. The correlation analysis showed that these 6 genes, that is, ankyrin repeat domain 9 (ANKRD9), lactate dehydrogenase A (LDHA), transcription factor 12 (TCF12), FAT atypical cadherin 1 (FAT1), lin-9 homolog (LIN9), and integrin beta 3 binding protein (ITGB3BP), can be used as effective candidate target genes of miR-29b-1-5p. Moreover, miR-29b-1-5p inhibits the expression of ANKRD9 by directly binding the 3′UTR of ANKRD9. Overall, these data indicate that miR-29b-1-5p inhibits the proliferation of primary chicken myoblasts, stimulates their differentiation, and is involved in the process of muscle development and that its effective target gene is ANKRD9. This study identified the molecular mechanism of miR-29b-1-5p in chicken muscle development.

To achieve sustainable development of the poultry industry, the effective conservation of genetic resources has become increasingly important. In the present study, we systematically elucidated the population structure, conservation priority, and runs of homozygosity (ROH) patterns of Chinese native chicken breeds. We used a high-density genotyping dataset of 157 native chickens from eight breeds. The population structure showed different degrees of population stratification among the breeds. Chahua chicken was the most differentiated breed from the other breeds (Nei = 0.0813), and the Wannan three-yellow chicken (WanTy) showed the lowest degree of differentiation (Nei = 0.0438). On the basis of contribution priority, Xiaoshan chicken had the highest contribution to the total gene diversity (1.41%) and the maximum gene diversity of the synthetic population (31.1%). WanTy chicken showed the highest contribution to the total allelic diversity (1.31%) and the maximum allelic diversity of the syntenic population (17.0%). A total of 5242 ROH fragments and 5 ROH island regions were detected. The longest ROH fragment was 41.51 Mb. A comparison of the overlapping genomic regions between the ROH islands and QTLs in the quantitative trait loci (QTL) database showed that the annotated candidate genes were involved in crucial economic traits such as immunity, carcass weight, drumstick and leg muscle development, egg quality and egg production, abdominal fat precipitation, body weight, and feed intake. In conclusion, our findings revealed that Chahua, Xiaoshan, and WanTy should be the priority conservation breeds, which will help optimize the conservation and breeding programs for Chinese indigenous chicken breeds.

Schizothorax o’connori (S. o’connori) is a representative tetraploid species in the subfamily Schizothoracinae and an important endemic fish in the Qinghai-Tibet Plateau. However, the domestication of S. o’connori remains challenging due to the lack of basic research. Here, we investigated the effects of artificial feeding on the oocytes and liver of S. o’connori by comparing the histological, metabolomic, and transcriptomic data. Histological results showed that the oocytes and liver of captive-reared S. o’connori had abnormal cell morphology. After comparison with the self-built database, a total of 233 metabolites were annotated. In oocytes, a total of 37 differentially accumulated metabolites (DAMs) were detected and two pathways were significantly enriched. There were obvious differences in the metabolites related to ovarian development, including pregnenolone and arachidonic acid. In liver, a total of 70 DAMs were detected and five pathways were significantly enriched. Based on the transcriptomic data, a total of 159 differentially expressed genes (DEGs) were significantly related with cell growth and death pathway in oocytes, while a total of 2841 DEGs were significantly related with 102 pathways in liver. Comparing the metabolomic and transcriptomic data showed that there were three common significant enrichment pathways in liver, including biosynthesis of unsaturated fatty acids, starch and sucrose metabolism, and fatty acid biosynthesis. These results showed that special attention should be given to the composition and intake of fatty acids during the artificial breeding of S. o’connori. In addition, many of metabolite-gene pairs were related to adenosine 5′-diphosphate, adenosine monophosphate, and pregnenolone. In summary, these data provide an overview of global metabolic and transcriptomic resources and broaden our understanding of captive-reared S. o’connori.

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

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.

Circular RNAs (circRNAs) play important roles in adipogenesis. However, studies on circRNA expression profiles associated with the development of abdominal adipose tissue are lacking in chickens. In this study, 12 cDNA libraries were constructed from the abdominal adipose tissue of Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 1,766 circRNAs were identified by Illumina HiSeq 2500 sequencing. These circRNAs were primarily distributed on chr1 through chr10 and sex chromosomes, and 84.95% of the circRNAs were from gene exons. Bioinformatic analysis showed that each circRNA has 35 miRNA binding sites on average, and 62.71% have internal ribosome entry site (IRES) elements. Meanwhile, these circRNAs were primarily concentrated in TPM < 0.1 and TPM > 60, and their numbers accounted for 18.90% and 80.51%, respectively, exhibiting specific expression patterns in chicken abdominal adipose tissue. In addition, 275 differentially expressed (DE) circRNAs were identified by comparison analysis. Functional enrichment analysis showed that the parental genes of DE circRNAs were primarily involved in biological processes and pathways related to lipid metabolism, such as regulation of fat cell differentiation, fatty acid homeostasis, and triglyceride homeostasis, as well as fatty acid biosynthesis, fatty acid metabolism, and glycerolipid metabolism. Furthermore, ceRNA regulatory networks related to abdominal adipose development were constructed. The results of this study indicated that circRNAs can regulate lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens through complex ceRNA networks between circRNAs, miRNAs, genes, and pathways. The results of this study may help to expand the number of known circRNAs in abdominal adipose tissue and provide a valuable resource for further research on the function of circRNAs in chicken abdominal adipose tissue.

Intramuscular fat (IMF), which regulated by genetics, nutrition and environment is an important factor that influencing meat quality. Up to now, the epigenetic regulation mechanism underlying poultry IMF deposition remains poorly understood. Here, we focused on the DNA methylation, which usually regulate genes in transcription level. To look into the essential role of DNA methylation on the IMF deposition, chicken intramuscular preadipocytes were isolated and cultured in vitro, and a model of intramuscular adipocyte differentiation was constructed. Combined the whole genome bisulfite sequencing (WGBS) and RNA-Seq technologies, we identified several methylated genes, which mainly affecting fatty acid metabolism and muscle development. Furthermore, we reported that DNA methylation regulate intramuscular adipogenesis by regulating the genes, such as collagen, type VI, alpha 1 (COL6A1) thus affecting IMF deposition. Overexpression of COL6A1 increases the lipid droplet and inhibits cell proliferation by regulating CHAD and CAMK2 in intramuscular adipocytes, while knockdown of COL6A1 shows the opposite effect. Taken together, our results reveal that DNA methylation plays an important role in poultry IMF deposition.

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.

Objective Polymorphisms occurring in the precursor region of microRNAs (miRNAs) affect the target gene and alter the biogenesis of miRNAs, resulting in phenotypic variation. The purpose of the study was to investigate the genetic effects of rs16681031 (C>G) mutation in the precursor region of gga-miR-1658 on the economic traits of the Gushi-Anka chicken F2 resource population. Methods To explore the effect of miR-1658 polymorphisms on chicken economic traits, the SNP was genotyped by MassArray matrix-assisted laser desorption/ionization-time of flight mass spectrometry. The association between the SNP and chicken body size, growth and carcass traits was determined by linear mixed models. Results The SNP was not only significantly associated with body weight at the age of 6, 8, 10, 12 weeks, respectively, but also with the breadth of the chicken chest, body slanting length and pelvic breadth at 4 weeks, chest depth at 8 weeks of age, and body slanting length at 12 weeks (p

The Xichuan black-bone chicken, which is a rare local chicken species in China, is an important genetic resource of black-bone chickens. Tyrosine can affect melanin production, but the molecular mechanism underlying tyrosine-induced melanin deposition in Xichuan black-bone chickens is poorly understood. Here, the blackness degree and melanin content of the breast muscle of Xichuan black-bone chickens fed a basic diet with five levels of added tyrosine (i.e., 0.2, 0.4, 0.6, 0.8, and 1.0%; these groups were denoted test groups I-V, respectively) were assessed, and the results showed that 0.8% tyrosine was the optimal level of added tyrosine. Moreover, the effects of tyrosine supplementation on the proliferation and tyrosinase content of melanocytes in Xichuan black-bone chickens were evaluated. The results revealed a dose-dependent relationship between tyrosine supplementation and melanocyte proliferation. In addition, 417 differentially expressed genes (DEGs), including 160 upregulated genes and 257 downregulated genes, were identified in a comparative analysis of the transcriptome profiles constructed using the pooled total RNA from breast muscle tissues of the control group and test group IV, respectively (fold change ≥2.0, P < 0.05). These DEGs were mainly involved in melanogenesis, the calcium signaling pathway, the Wnt signaling pathway, the mTOR signaling pathway, and vascular smooth muscle contraction. The pathway analysis of the DEGs identified some key genes associated with pigmentation, such as DCT and EDNRB2. In summary, the melanin content of breast muscle could be markedly enhanced by adding an appropriate amount of tyrosine to the diet of Xichuan black-bone chickens, and the EDNRB2-mediated molecular regulatory network could play a key role in the biological process of tyrosine-induced melanin deposition. These results have deepened the understanding of the molecular regulatory mechanism of melanin deposition in black-bone chickens and provide a basis for the regulation of nutrition and genetic breeding associated with melanin deposition in Xichuan black-bone chickens.

There is a lack of understanding surrounding the molecular mechanisms involved in the development of chicken skeletal muscle in the late postnatal stage, especially in the regulation of breast muscle development related genes, pathways, miRNAs and other factors. In this study, 12 cDNA libraries and 4 small RNA libraries were constructed from Gushi chicken breast muscle samples from 6, 14, 22, and 30 weeks. A total of 15,508 known transcripts, 25,718 novel transcripts, 388 known miRNAs and 31 novel miRNAs were identified by RNA-seq in breast muscle at the four developmental stages. Through correlation analysis of miRNA and mRNA expression profiles, it was found that 417, 370, 240, 1,418, 496, and 363 negatively correlated miRNA–mRNA pairs of W14 vs. W6, W22 vs. W6, W22 vs. W14, W30 vs. W6, W30 vs. W14, and W30 vs. W22 comparisons, respectively. Based on the annotation analysis of these miRNA–mRNA pairs, we constructed the miRNA–mRNA interaction network related to biological processes, such as muscle cell differentiation, striated muscle tissue development and skeletal muscle cell differentiation. The interaction networks for signaling pathways related to five KEGG pathways (the focal adhesion, ECM-receptor interaction, FoxO signaling, cell cycle, and p53 signaling pathways) and PPI networks were also constructed. We found that ANKRD1, EYA2, JSC, AGT, MYBPC3, MYH11, ACTC1, FHL2, RCAN1, FOS, EGR1, and FOXO3, PTEN, AKT1, GADD45, PLK1, CCNB2, CCNB3 and other genes were the key core nodes of these networks, most of which are targets of miRNAs. The FoxO signaling pathway was in the center of the five pathway-related networks. In the PPI network, there was a clear interaction among PLK1 and CDK1, CCNB2, CDK1, and GADD45B, and CDC45, ORC1 and MCM3 genes. These results increase the understanding for the molecular mechanisms of chicken breast muscle development, and also provide a basis for studying the interactions between genes and miRNAs, as well as the functions of the pathways involved in postnatal developmental regulation of chicken breast muscle.

The formation mechanism underlying the blue eggshell characteristic has been discovered in birds, and SLCO1B3 is the key gene that regulates the blue eggshell color. Insertion of an endogenous retrovirus, EAV-HP, in the SLCO1B3 5′ flanking region promotes SLCO1B3 expression in the chicken shell gland, and this expression causes bile salts to enter the shell gland, where biliverdin is secreted into the eggshell, forming a blue shell. However, at different laying stages of the same group of chickens, the color of the eggshell can vary widely, and the molecular mechanism underlying the eggshell color change remains unknown. Therefore, to reveal the molecular mechanism of the blue eggshell color variations, we analyzed the change in the eggshell color during the laying period. The results indicated that the eggshell color in Lushi chickens can be divided into three stages: 20–25 weeks for dark blue, 26–45 weeks for medium blue, and 46–60 weeks for light blue. We further investigated the expression and methylation levels of the SLCO1B3 gene at eight different weeks, finding that the relative expression of SLCO1B3 was significantly higher at 25 and 30 weeks than at other laying weeks. Furthermore, the overall methylation rate of the SLCO1B3 gene in Lushi chickens increased gradually with increasing weeks of egg production, as shown by bisulfite sequencing PCR. Pearson correlation analysis showed that methylation of the promoter region of SLCO1B3 was significantly negatively correlated with both SLCO1B3 expression in the shell gland tissue and eggshell color. In addition, we predicted that CpG5 and CpG8 may be key sites for regulating SLCO1B3 gene transcription. Our findings show that as the level of methylation increases, methylation of the CpG5 and CpG8 sites hinders the binding of transcription factors to the promoter, reducing SLCO1B3 expression during the late period and resulting in a lighter eggshell color.

Long noncoding RNAs (lncRNAs) play important roles in transcriptional and posttranscriptional regulation. However, the effects of lncRNAs on the meat quality of chicken hasn't been elucidated clearly yet. Gushi chickens are popular in China because of their superior meat quality, particularly the tender flesh, and unique flavor. Gushi chickens are popular in China because of their superior meat quality, delicate flesh, and unique flavor. We performed RNA-Seq analysis of breast muscle from Gushi chicken at two physiological stages, including juvenile (G20W) and laying (G55W). In total, 186 lncRNAs and 881 mRNAs were differentially expressed between G20W and G55W (fold change ≥ 2.0, P < 0.05). Among them, 131 lncRNAs presented upregulated and 55 were downregulated. We identified the cis and trans target genes of the differentially expressed lncRNAs, and constructed lncRNA-mRNA interaction networks. The results showed that differentially expressed mRNAs and lncRNAs were mainly involved in ECM-receptor interaction, glycerophospholipid metabolism, ubiquitin-mediated proteolysis, and the biosynthesis of amino acids. In summary, our study utilized RNA-seq analysis to predict the functions of lncRNA on chicken meat quality. Furthermore, comprehensive analysis identified lncRNAs and their target genes, which may contribute to a better understanding of the molecular mechanisms underlying in poultry meat quality and provide a theoretical basis for further research.

BACKGROUND:Coloration is one of the most recognizable characteristics in chickens, and clarifying the coloration mechanisms will help us understand feather color formation. "Yufen I" is a commercial egg-laying chicken breed in China that was developed by a three-line cross using lines H, N and D. Columbian plumage is a typical feather character of the "Yufen I" H line. To elucidate the molecular mechanism underlying the pigmentation of Columbian plumage, this study utilizes high-throughput sequencing technology to compare the transcriptome and proteome differences in the follicular tissue of different feathers, including the dorsal neck with black and white striped feather follicles (Group A) and the ventral neck with white feather follicles (Group B) in the "Yufen I" H line. RESULTS:In this study, we identified a total of 21,306 genes and 5,203 proteins in chicken feather follicles. Among these, 209 genes and 382 proteins were differentially expressed in two locations, Group A and Group B, respectively. A total of 8 differentially expressed genes (DEGs) and 9 differentially expressed proteins (DEPs) were found to be involved in the melanogenesis pathway. Additionally, a specifically expressed MED23 gene and a differentially expressed GNAQ protein were involved in melanin synthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapped 190 DEGs and 322 DEPs to 175 and 242 pathways, respectively, and there were 166 pathways correlated with both DEGs and DEPs. 49 DEPs/DEGs overlapped and were enriched for 12 pathways. Transcriptomic and proteomic analyses revealed that the following pathways were activated: melanogenesis, cardiomyocyte adrenergic, calcium and cGMP-PKG. The expression of DEGs was validated by real-time quantitative polymerase chain reaction (qRT-PCR) that produced results similar to those from RNA-seq. In addition, we found that the expression of the MED23, FZD10, WNT7B and WNT11 genes peaked at approximately 8 weeks in the "Yufen I" H line, which is consistent with the molting cycle. As both groups showed significant differences in terms of the expression of the studied genes, this work opens up avenues for research in the future to assess their exact function in determining plumage color. CONCLUSION:Common DEGs and DEPs were enriched in the melanogenesis pathway. MED23 and GNAQ were also reported to play a crucial role in melanin synthesis. In addition, this study is the first to reveal gene and protein variations in in the "Yufen I" H line during Columbian feather color development and to discover principal genes and proteins that will aid in functional genomics studies in the future. The results of the present study provide a significant conceptual basis for the future breeding schemes with the "Yufen I" H line and provide a basis for research on the mechanisms of feather pigmentation.

The excessive deposition of abdominal fat has become an important factor in restricting the production efficiency of chickens, so reducing abdominal fat deposition is important for improving growth rate. It has been proven that miRNAs play an important role in regulating many physiological processes of organisms. In this study, we constructed a model of adipogenesis by isolating preadipocytes (Ab-Pre) derived from abdominal adipose tissue and differentiated adipocytes (Ab-Ad) in vitro. Deep sequencing of miRNAs and mRNAs expressed in Ab-Pre and Ab-Ad groups was conducted to explore the effect of miRNAs and mRNAs on fat deposition. We identified 80 differentially expressed miRNAs (DEMs) candidates, 58 of which were up-regulated and 22 down-regulated. Furthermore, six miRNAs and six mRNAs were verified by qRT-PCR, and the results showed that the expression of the DEMs and differentially expressed genes (DEGs) in the two groups was consistent with our sequencing results. When target genes of miRNA were combined with mRNA transcriptome data, a total of 891 intersection genes were obtained, we predicted the signal pathways of cross genes enrichment to the MAPK signal pathway, insulin signal pathway, fatty acid metabolism, and ECM−receptor interaction. Meanwhile, we constructed miRNA and negatively correlated mRNA target networks, including 12 miRNA−mRNAs pairs, which showed a strong association with the abdominal adipocyte differentiation (miR-214−ACSBG2, NFKB2, CAMK2A, ACLY, CCND3, PLK3, ITGB2; miR-148a-5p−ROCK2; miR-10a-5p−ELOVL5; miR-146b-5p−LAMA4; miR-6615-5p−FLNB; miR-1774−COL6A1). Overall, these findings provide a background for further research on lipid metabolism. Thus, we can better understand the molecular genetic mechanism of chicken abdominal fat deposition.

In order to investigate the effects of fasting caecotrophy on hepatic lipid metabolism in rabbits, 12 weaned female New Zealand white rabbits were randomly divided into (n = 6/group) a control and fasting caecotrophy group. Rabbits in the experimental group were treated with an Elizabeth circle to prevent them from eating their own soft feces for a 60-day period. Growth and blood biochemical indices, transcriptome sequencing and histology analysis of the liver were performed. Compared with the control group, final weight, weight gain, liver weight, growth rate and feed conversion ratio, all decreased in the experimental group (p < 0.05). RNA sequencing (RNA-seq) analysis revealed a total of 301.2 million raw reads (approximately 45.06 Gb of high-quality clean data) that were mapped to the rabbit genome. After a five-step filtering process, 14,964 genes were identified, including 444 differentially expressed genes (p < 0.05, foldchange ≥ 1). A number of differently expressed genes linked to lipid metabolism were further analyzed including CYP7A1, SREBP, ABCA1, GPAM, CYP3A1, RBP4 and RDH5. The KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation of the differentially expressed genes indicated that main pathways affected were pentose and glucuronide interactions, starch and sucrose metabolism, retinol metabolism and PPAR signaling. Overall, the present study revealed that preventing caecotrophy reduced growth and altered lipid metabolism, both of which will help guide the development of new approaches for rabbits’ feeding and production. These data also provide a reference for studying the effects of soft feces in other small herbivores.

Rumen development in calves is affected by many factors, including dietary composition. MicroRNAs (miRNAs) are known to function in the development of the rumen in cattle, what is not known is how these miRNAs function in rumen development of calves fed with high and low ratios of non-fibrous carbohydrate (NFC)/neutral detergent fiber (NDF). A total of six healthy Charolais hybrids bull calves of similar weight were divided into two groups; three calves were fed a mixed diet with NFC/NDF = 1.35 (H group), and three were fed a mixed diet with NFC/NDF = 0.80 (L group). After 105 days on the diet, calves were sacrificed and rumen tissues were collected. Tissues were subjected to histological observation and miRNA expression analysis. Functional enrichment analysis was conducted on the target genes of the miRNAs. Targeting and regulatory relationships were verified by luciferase reporter assay and quantitative PCR (qPCR). We found that the length of rumen papilla in the L group was significantly greater than that in the H group, while the width of rumen papilla in H group was significantly greater than that that in L group. We identified 896 miRNAs; 540 known miRNAs, and 356 novel predicted miRNAs. After statistical testing, we identified 24 differentially expressed miRNAs (DEmiRNAs). miRNA-mRNA-cluster network analysis and literature reviews revealed that cell proliferation, differentiation, physical and nutrient stimuli processes participate in rumen development under different NFC/NDF levels. The regulatory relationships between three DEmiRNAs and five target genes were verified by examining the levels of expression. The binding sites on bta-miR-128 for the peroxisome proliferator activated receptor gamma (PPARG) and solute carrier family 16 member 1 (SLC16A1) genes were investigated using a dual luciferase assay. The results of this study provide insight into the role of miRNAs in rumen development in calves under different NFC/NDF levels.

Intramuscular fat (IMF) is the most important evaluating indicator of chicken meat quality, the content of which is positively correlated with tenderness, flavor, and succulence of the meat. Chicken IMF deposition process is regulated by many factors, including genetic, nutrition, and environment. Although large number of omics’ studies focused on the IMF deposition process, the molecular mechanism of chicken IMF deposition is still poorly understood. In order to study the role of miRNAs in chicken intramuscular adipogenesis, the intramuscular adipocyte differentiation model (IMF-preadipocytes and IMF-adipocytes) was established and subject to miRNA-Seq. A total of 117 differentially expressed miRNAs between two groups were obtained. Target genes prediction and functional enrichment analysis revealed that eight pathways involved in lipid metabolism related processes, such as fatty acid metabolism and fatty acid elongation. Meanwhile a putative miRNA, gga-miR-18b-3p, was identified be served a function in the intramuscular adipocyte differentiation. Luciferase assay suggested that the gga-miR-18b-3p targeted to the 3′UTR of ACOT13. Subsequent functional experiments demonstrated that gga-miR-18b-3p acted as an inhibitor of intramuscular adipocyte differentiation by targeting ACOT13. Our findings laid a new theoretical foundation for the study of lipid metabolism, and also provided a potential target to improve the meat quality in the poultry industry.

Lmbr1 has been reported to be associated with chicken polydactyly and carcass traits. Based on the hypothesis that chicken Lmbr1 should have multiple alternative splices as human and mouse Lmbr1/C7orf2, we successfully identified one novel chicken Lmbr1 transcript variant (designated as Lmbr1-2) by bioinformatics analysis and RT-PCR confirmation. Chicken Lmbr1-2 transcript was predicted to encode a 30 amino acid protein. Lmbr1-2 transcript was ubiquitously expressed in a range of chicken tissues. A SNP (rs14135851: G>T, also named as 288G>T reference to 603234777F1) was found to be located at the 5′ boundary (first base) of alternatively spliced exon (IN2) of chicken Lmbr1-2. Eight variants/10 haplotypes were detected from 4 breeds in a 300-bp genomic fragment surrounding 288G>T. The 288T>G and 331G>A changes were predicted resulting to the formation/loss of SRSF5 and SRSF6 motifs. The genotype/allele distribution for a G/A variant (331G>A, reference to 603234777F1) and haplotype types in polydactylous Silkie presented distinct difference with that in other four-toed breed chickens. Association analysis showed that chicken 331G>A variation had significant effect on meat quality traits including breast muscle water loss rate, leg muscle water loss rate and leg muscle fiber density, and carcass traits including evisceration weight, pancreas weight in Gushi chicken F2 resource population. These data demonstrated that chicken Lmbr1 was alternatively spliced to produce multiple splice forms as was the case in mammals, and it suggested that chicken Lmbr1-2 variation might have important effect on the carcass and meat quality traits in chicken.