Background: Since the domestication of chicken, various breeds have been developed for food production, entertainment, and so on. Compared to indigenous chicken breeds which generally do not show elite production performance, commercial breeds or lines are selected intensely for meat or egg production. In the present study, in order to understand the molecular mechanisms underlying the dramatic differences of egg number between commercial egg-type chickens and indigenous chickens, we performed a genome-wide association study (GWAS) in a mixed linear model. Results: We obtained 148 single nucleotide polymorphisms (SNPs) associated with egg number traits (57 significantly, 91 suggestively). Among them, 4 SNPs overlapped with previously reported quantitative trait loci (QTL), including 2 for egg production and 2 for reproductive traits. Furthermore, we identified 32 candidate genes based on the function of the screened genes. These genes were found to be mainly involved in regulating hormones, playing a role in the formation, growth, and development of follicles, and in the development of the reproductive system. Some genes such as NELL2 (neural EGFL like 2), KITLG (KIT ligand), GHRHR (Growth hormone releasing hormone receptor), NCOA1 (Nuclear receptor coactivator 1), ITPR1 (inositol 1, 4, 5-trisphosphate receptor type 1), GAMT (guanidinoacetate N-methyltransferase), and CAMK4 (calcium/calmodulin-dependent protein kinase IV) deserve our attention and further study since they have been reported to be closely related to egg production, egg number and reproductive traits. In addition, the most significant genomic region obtained in this study was located at 48.61–48.84 Mb on GGA5. In this region, we have repeatedly identified four genes, in which YY1 (YY1 transcription factor) and WDR25 (WD repeat domain 25) have been shown to be related to oocytes and reproductive tissues, respectively, which implies that this region may be a candidate region underlying egg number traits. Conclusion: Our study utilized the genomic information from various chicken breeds or populations differed in the average annual egg number to understand the molecular genetic mechanisms involved in egg number traits. We identified a series of SNPs, candidate genes, or genomic regions that associated with egg number, which could help us in developing the egg production trait in chickens.

The study was conducted to investigate the transcriptomic differences of the hypothalamic-pituitary-gonadal axis between Xinjiang Yili geese with high and low egg production and to find candidate genes regulating the egg production of Xinjiang Yili geese. The 8 selected Xinjiang Yili Geese with high or low egg production (4 for each group) were 3 years old, with good health, and under the same feeding condition. High-throughput sequencing technology was used to sequence cDNA libraries of the hypothalami, pituitary glands, and ovaries. The sequencing data were compared and analyzed, and the transcripts with significant differences were identified and analyzed with bioinformatics. The study showed that the transcriptome sequencing data of the 24 samples contained a total of 1,176,496,146 valid reads and 176.47 gigabase data. Differential expression analyses identified 135, 56, and 331 genes in the hypothalami, pituitary glands, and ovaries of Xinjiang Yili geese with high and low egg production. Further annotation of these differentially expressed genes in the non-redundant protein sequence database (Nr) revealed that 98, 52, and 309 genes were annotated, respectively. Through the annotations of GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases, 30 candidate genes related to the egg production of Xinjiang Yili geese were preliminarily selected. The gap junction, focal adhesion, and ECM-receptor interaction signaling pathways were enriched with the hypothalamic, pituitary, and ovarian differentially expressed genes, and the calcium signaling pathway was enriched with the pituitary and ovarian differentially expressed genes. Thus, these pathways in the hypothalamic-pituitary-gonadal axis may play an important role in regulating egg production of Xinjiang Yili geese. The results provided the transcriptomic information of the hypothalamic-pituitary-gonadal axis of Xinjiang Yili geese and laid the theoretical basis for revealing the molecular mechanisms regulating the egg-laying traits of Xinjiang Yili geese.

The development of follicles in the ovaries is a critical determinant of poultry egg production. There are existing studies on the follicular development patterns in poultry, but the specific regulatory mechanisms still need further study. In a previous study, we identified long non-coding RNA (lncRNA) MSTRG.5970.28, anosmin 1 (ANOS1), and its predicted target miR-133a-3p that may be associated with goose ovary development. However, the function of MSTRG.5970.28 in goose granulosa cells and its regulatory mechanisms affecting granulosa cell proliferation and apoptosis have not been reported. In the present study, MSTRG.5970.28 and miR-133a-3p overexpression and interference vectors were constructed. Combined with reverse-transcription real-time quantitative PCR (RT-qPCR), a dual luciferase activity assay, Cell Counting Kit-8 (CCK-8), and flow cytometric analysis, we investigated the role of the MSTRG.5970.28-miR-133a-3p-ANOS1 axis in goose follicular granulosa cells and the associated regulatory mechanisms. MSTRG.5970.28 was found to be localized in the cytoplasm and its expression was influenced by reproductive hormones. The targeting relationship among MSTRG.5970.28, ANOS1, and miR-133a-3p were verified by a dual luciferase activity assay. CCK-8 and apoptosis assays showed that MSTRG.5970.28 inhibited the proliferation and promoted apoptosis of goose granulosa cells. The regulatory role of miR-133a-3p on granulosa cell proliferation and apoptosis was opposite to MSTRG.5970.28. We found that the proliferative and apoptotic effects of granulosa cells caused by MSTRG.5970.28 overexpression were attenuated by miR-133a-3p. MSTRG.5970.28 functions as a competitive endogenous RNA that regulates ANOS1 expression by sponging miR-133a-3p and thus exerts regulatory functions in granulosa cells. In sum, the present study identified lncRNA MSTRG.5970.28 as associated with goose ovary development, which affects the expression of ANOS1 by targeting miR-133a-3p, thereby influencing the proliferation and apoptosis of goose granulosa cells.

Molting in birds provides us with an ideal genetic model for understanding aging and rejuvenation since birds present younger characteristics for reproduction and appearance after molting. Forced molting (FM) by fasting in chickens causes aging of their reproductive system and then promotes cell redevelopment by providing water and feed again. To reveal the genetic mechanism of rejuvenation, we detected blood hormone indexes and gene expression levels in the hypothalamus and ovary of hens from five different periods during FM. Three hormones were identified as participating in FM. Furthermore, the variation trends of gene expression levels in the hypothalamus and ovary at five different stages were found to be basically similar using transcriptome analysis. Among them, 45 genes were found to regulate cell aging during fasting stress and 12 genes were found to promote cell development during the recovery period in the hypothalamus. In addition, five hub genes (INO80D, HELZ, AGO4, ROCK2, and RFX7) were identified by WGCNA. FM can restart the reproductive function of aged hens by regulating expression levels of genes associated with aging and development. Our study not only enriches the theoretical basis of FM but also provides insights for the study of antiaging in humans and the conception mechanism in elderly women.

Background: The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development. Results: In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145–5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed. Conclusions: Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

Plumage color is an artificially and naturally selected trait in domestic ducks. Black, white, and spotty are the main feather colors in domestic ducks. Previous studies have shown that black plumage color is caused by MC1R, and white plumage color is caused by MITF. We performed a genome-wide association study (GWAS) to identify candidate genes associated with white, black, and spotty plumage in ducks. Two non-synonymous SNPs in MC1R (c.52G>A and c.376G>A) were significantly related to duck black plumage, and three SNPs in MITF (chr13:15411658A>G, chr13:15412570T>C and chr13:15412592C>G) were associated with white plumage. Additionally, we also identified the epistatic interactions between causing loci. Some ducks with white plumage carry the c.52G>A and c.376G>A in MC1R, which also compensated for black and spotty plumage color phenotypes, suggesting that MC1R and MITF have an epistatic effect. The MITF locus was supposed to be an upstream gene to MC1R underlying the white, black, and spotty colors. Although the specific mechanism remains to be further clarified, these findings support the importance of epistasis in plumage color variation in ducks.