Genome sequencing of drake semen micobiome with correlation with their compositions, sources and potential mechanisms affecting semen quality.

Artificial insemination (AI) technology has greatly promoted the development of the chicken industry. Recently, AI technology has also begun to be used in the duck industry, but there are some problems. Numerous researchers have shown that microbes colonizing in semen can degrade semen quality, and AI can increase the harmful microbial load in hen's reproductive tract. Different from the degraded external genitalia of roosters, drakes have well-developed external genitalia, which may cause drake semen to be more susceptible to microbial contamination. However, information on the compositions, sources, and effects of semen microbes on semen quality remains unknown in drakes. In the current study, high-throughput sequencing technology was used to detect microbial communities in drake semen, environmental swabs, cloacal swabs, and the spermaduct after quantifying the semen quality of drakes to investigate the effects of microbes in the environment, cloaca, and spermaduct on semen microbiota and the relationships between semen microbes and semen quality. Taxonomic analysis showed that the microbes in the semen, environment, cloaca, and spermaduct samples were all classified into 4 phyla and 25 genera. Firmicutes and Proteobacteria were the dominant phyla. Phyllobacterium only existed in the environment, while Marinococcus did not exist in the cloaca. Of the 24 genera present in semen: Brachybacterium, Brochothrix, Chryseobacterium, Kocuria, Marinococcus, Micrococcus, Rothia, Salinicoccus, and Staphylococcus originated from the environment; Achromobacter, Aerococcus, Corynebacterium, Desemzia, Enterococcus, Jeotgalicoccus, Pseudomonas, Psychrobacter, and Turicibacter originated from the cloaca; and Agrobacterium, Carnobacterium, Chelativorans, Devosia, Halomonas, and Oceanicaulis originated from the spermaduct. In addition, K-means clustering analysis showed that semen samples could be divided into 2 clusters based on microbial compositions, and compared with cluster 1, the counts of Chelativorans (P < 0.05), Devosia (P < 0.01), Halomonas (P < 0.05), and Oceanicaulis (P < 0.05) were higher in cluster 2, while the sperm viability (P < 0.05), total sperm number (P < 0.01), and semen quality factor (SQF) (P < 0.01) were lower in cluster 2. Furthermore, functional prediction analysis of microbes showed that the activities of starch and sucrose metabolism, phosphotransferase system, ABC transporters, microbial metabolism in diverse environments, and quorum sensing pathways between cluster 1 and cluster 2 were significantly different (P < 0.05). Overall, environmental/cloacal microbes resulted in semen contamination, and microbes from the Chelativorans, Devosia, Halomonas, and Oceanicaulis genera may have negative effects on semen quality in drakes by affecting the activities of starch and sucrose metabolism, phosphotransferase system, ABC transporters, and quorum sensing pathways that are associated with carbohydrate metabolism. These data will provide a basis for developing strategies to prevent microbial contamination of drake semen.

Comparative transcriptomics analysis identifies crucial genes and pathways during goose spleen development.

As the largest peripheral lymphoid organ in poultry, the spleen plays an essential role in regulating the body's immune capacity. However, compared with chickens and ducks, information about the age- and breed-related changes in the goose spleen remains scarce. In this study, we systematically analyzed and compared the age-dependent changes in the morphological, histological, and transcriptomic characteristics between Landes goose (LG; Anser anser) and Sichuan White goose (SWG; Anser cygnoides). The results showed a gradual increase in the splenic weights for both LG and SWG until week 10, while their splenic organ indexes reached the peak at week 6. Meanwhile, the splenic histological indexes of both goose breeds continuously increased with age, reaching the highest levels at week 30. The red pulp (RP) area was significantly higher in SWG than in LG at week 0, while the splenic corpuscle (AL) diameter was significantly larger in LG than in SWG at week 30. At the transcriptomic level, a total of 1710 and 1266 differentially expressed genes (DEGs) between week 0 and week 30 were identified in spleens of LG and SWG, respectively. Meanwhile, a total of 911 and 808 DEGs in spleens between LG and SWG were identified at weeks 0 and 30, respectively. Both GO and KEGG enrichment analysis showed that the age-related DEGs of LG or SWG were dominantly enriched in the Cell cycle, TGF-beta signaling, and Wnt signaling pathways, while most of the breed-related DEGs were enriched in the Neuroactive ligand-receptor interaction, Cytokine-cytokine receptor interaction, ECM-receptor interaction, and metabolic pathways. Furthermore, through construction of protein-protein interaction networks using significant DEGs, it was inferred that three hub genes including BUB1, BUB1B, and TTK could play crucial roles in regulating age-dependent goose spleen development while GRIA2, GRIA4, and RYR2 could be crucial for the breed-specific goose spleen development. These data provide novel insights into the splenic developmental differences between Chinese and European domestic geese, and the identified crucial pathways and genes are helpful for a better understanding of the mechanisms regulating goose immune functions.

Effects of cage vs. net-floor mixed rearing system on goose spleen histomorphology and gene expression profiles.

Due to the demands for both environmental protection and modernization of the goose industry in China, the traditional goose waterside rearing systems have been gradually transitioning to the modern intensive dryland rearing ones, such as the net-floor mixed rearing system (MRS) and cage rearing system (CRS). However, the goose immune responses to different dryland rearing systems remain poorly understood. This study aimed to investigate and compare the age-dependent effects of MRS and CRS on the splenic histomorphological characteristics and immune-related genes expression profiles among three economically important goose breeds, including Sichuan White goose (SW), Gang goose (GE), and Landes goose (LD). Morphological analysis revealed that the splenic weight and organ index of SW were higher under CRS than under MRS (p < 0.05). Histological observations showed that for SW and LD, the splenic corpuscle diameter and area as well as trabecular artery diameter were larger under MRS than under CRS at 30 or 43 weeks of age (p < 0.05), while the splenic red pulp area of GE was larger under CRS than under MRS at 43 weeks of age (p < 0.05). Besides, at 43 weeks of age, higher mRNA expression levels of NGF, SPI1, and VEGFA in spleens of SW were observed under MRS than under CRS (p < 0.05), while higher levels of HSPA2 and NGF in spleens of LD were observed under MRS than under CRS (p < 0.05). For GE, there were higher mRNA expression levels of MYD88 in spleens under CRS at 30 weeks of age (p < 0.05). Moreover, our correlation analysis showed that there appeared to be more pronounced positive associations between the splenic histological parameters and expression levels of several key immune-related genes under MRS than under CRS. Therefore, it is speculated that the geese reared under MRS might exhibit enhanced immune functions than those under CRS, particularly for SW and LD. Although these phenotypic differences are assumed to be associated with the age-dependent differential expression profiles of HSPA2, MYD88, NGF, SPI1, and VEGFA in the goose spleen, the underlying regulatory mechanisms await further investigations.

Multi-omics reveals goose fatty liver formation from metabolic reprogramming.

To comprehensively provide insight into goose fatty liver formation, we performed an integrative analysis of the liver transcriptome, lipidome, and amino acid metabolome, as well as peripheral adipose tissue transcriptome analysis using samples collected from the overfed geese and normally fed geese. Transcriptome analysis showed that liver metabolism pathways were mainly enriched in glucolipid metabolism, amino acid metabolism, inflammation response, and cell cycle; peripheral adipose tissue and the liver cooperatively regulated liver lipid accumulation during overfeeding. Liver lipidome patterns obviously changed after overfeeding, and 157 different lipids were yielded. In the liver amino acid metabolome, the level of Lys increased after overfeeding. In summary, this is the first study describing goose fatty liver formation from an integrative analysis of transcriptome, lipidome, and amino acid metabolome, which will provide a whole new dimension to understanding the mechanism of goose fatty liver formation.

Molecular mechanisms of hypothalamic-pituitary-ovarian/thyroid axis regulating age at first egg in geese.

Age at first egg (AFE) has consistently garnered interest as a crucial reproductive indicator within poultry production. Previous studies have elucidated the involvement of the hypothalamic-pituitary-ovarian (HPO) and hypothalamic-pituitary-thyroid (HPT) axes in regulating poultry sexual maturity. Concurrently, there was evidence suggesting a potential co-regulatory relationship between these 2 axes. However, as of now, no comprehensive exploration of the key pathways and genes responsible for the crosstalk between the HPO and HPT axes in the regulation of AFE has been reported. In this study, we conducted a comparative analysis of morphological differences and performed transcriptomic analysis on the hypothalamus, pituitary, thyroid, and ovarian stroma between normal laying group (NG) and abnormal laying group (AG). Morphological results showed that the thyroid index difference (D-) value (thyroid index D-value=right thyroid index-left thyroid index) was significantly (P < 0.05) lower in the NG than in the AG, while the ovarian index was significantly (P < 0.01) higher in the NG than in the AG. Furthermore, between NG and AG, we identified 99, 415, 167, and 1182 differentially expressed genes (DEGs) in the hypothalamus, pituitary, thyroid, and ovarian stroma, respectively. Gene ontology (GO) analysis highlighted that DEGs from 4 tissues were predominantly enriched in the "biological processes" category. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that 16, 14, 3, and 26 KEGG pathways were significantly enriched (P < 0.05) in the hypothalamus, pituitary, thyroid, and ovarian stroma. The MAPK signaling pathway emerged as the sole enriched pathway across all 4 tissues. Employing an integrated analysis of the protein-protein interaction (PPI) network and correlation analysis, we found GREB1 emerged as a pivotal component within the HPO axis to regulate estrogen-related signaling in the HPT axis, meanwhile, the HPT axis influenced ovarian development by regulating thyroid hormone-related signaling mainly through OPN5. Then, 10 potential candidate genes were identified, namely IGF1, JUN, ERBB4, KDR, PGF, FGFR1, GREB1, OPN5, DIO3, and THRB. These findings establish a foundation for elucidating the physiological and genetic mechanisms by which the HPO and HPT axes co-regulate goose AFE.

Comparison of the Effects of Recombinant and Native Prolactin on the Proliferation and Apoptosis of Goose Granulosa Cells.

In poultry, prolactin (PRL) plays a key role in the regulation of incubation behavior, hormone secretion, and reproductive activities. However, previous in vitro studies have focused on the actions of PRL in ovarian follicles of poultry, relying on the use of exogenous or recombinant PRL, and the true role of PRL in regulating ovarian granulosa cell (GC) functions in poultry awaits a further investigation using endogenous native PRL. Therefore, in this study, we first isolated and purified recombinant goose PRL protein (rPRL) and native goose PRL protein (nPRL) using Ni-affinity chromatography and rabbit anti-rPRL antibodies-filled immunoaffinity chromatography, respectively. Then, we analyzed and compared the effects of rPRL and nPRL at different concentrations (0, 3, 30, or 300 ng/mL) on the proliferation and apoptosis of both GCs isolated from goose ovarian pre-hierarchical follicles (phGCs) and from hierarchical follicles (hGCs). Our results show that rPRL at lower concentrations increased the viability and proliferation of both phGCs and hGCs, while it exerted anti-apoptotic effects in phGCs by upregulating the expression of Bcl-2. On the other hand, nPRL increased the apoptosis of phGCs in a concentration-dependent manner by upregulating the expressions of caspase-3 and Fas and downregulating the expressions of Bcl-2 and Becn-1. In conclusion, this study not only obtained a highly pure nPRL for the first time, but also suggested a dual role of PRL in regulating the proliferation and apoptosis of goose GCs, depending on its concentration and the stage of follicle development. The data presented here can be helpful in purifying native proteins of poultry and enabling a better understanding of the roles of PRL during the ovarian follicle development in poultry.

Effects and Mechanisms of Cage versus Floor Rearing System on Goose Growth Performance and Immune Status.

Currently, FRS and CRS are the two predominant dryland rearing systems in the goose industry. However, the effects of these two systems on goose growth performance and health, as well as the underlying mechanisms, have not been fully clarified. Thus, this study aimed to compare growth performance and immune status, as well as investigate the genome-wide transcriptomic profiles of spleen in geese, between CRS and FRS at 270 d of age. Phenotypically, the body weight and body size traits were higher in geese under FRS, while the weight and organ index of spleen were higher in geese under CRS (p < 0.05). Noticeably, the bursa of Fabricius of geese under FRS was degenerated, while that under CRS was retained. At the serum level, the immune globulin-G (IgG) and interleukin-6 (IL-6) levels were higher in geese under CRS (p < 0.05). At the transcriptomic level, we identified 251 differentially expressed genes (DEGs) in the spleen between CRS and FRS, which were mainly enriched in scavenger receptor activity, inflammatory response, immune response, neuroactive ligand-receptor interaction, phenylalanine metabolism, ECM receptor interaction, calcium signaling pathway, phenylalanine, tyrosine, and tryptophan biosynthesis, regulation of actin cytoskeleton, and MAPK signaling pathways. Furthermore, through protein-protein interaction (PPI) network analysis, ten candidate genes were identified, namely, VEGFA, FGF2, NGF, GPC1, NKX2-5, FGFR1, FGF1, MEIS1, CD36, and PAH. Further analysis demonstrated that geese in CRS could improve their immune ability through the "phenylalanine metabolism" pathway. Our results revealed that the FRS improved growth performance, whereas the CRS improved goose immune function by increasing levels of IL-6 and IgG in serum. Moreover, the phenylalanine metabolism pathway could exert positive effects on immune function of geese under CRS. These results can provide reliable references for understanding how floor and cage rearing systems affect goose growth performance and immune capacity.

Comprehensive transcriptomic and metabolomic analysis of the effect of feed restriction on duck sternal development.

Skeletal characteristics are important to the growth and development of poultry. In feeding management, constant free feeding (FF) of poultry may lead to imbalance between bone development and weight gain. Feed restriction (FR), to a certain extent, is one way to solve this problem. However, the effect of feed restriction on poultry bone development needs further elucidation at the molecular level. Therefore, in the present study, we investigated the effects of different levels of feed restriction (60% FR, 70% FR, 80% FR, and FF) on the sternum development of ducks at 7 and 8 wk old. In the seventh wk, with increasing feed restriction, the values of traits including body weight, breast muscle weight, sternal weight, keel length, and calcified keel length decreased. However, in the eighth wk, the sternum weight and keel length of ducks treated with 60% FR were unexpectedly higher than those of FF individuals, indicative of catch-up growth. Then, we conducted RNA-seq and metabolomic analysis on sterna from 7- and 8-wk-old FF and 60% FR ducks. The results identified multiple differentially expressed genes (DEGs) associated with sternum development that were influenced by feed restriction. Among them, we found that the mRNA expression levels of the chondroitin sulfate synthase 3 (CHSY3) and annexin A2 (ANXA2) which are involved in glycosaminoglycan biosynthesis and bone mineralization, had smaller changes over time under FR treatment than under FF treatment, implying that the FR treatment to a certain extent prevented the premature calcification and prolonged the development time of duck sternum. In addition, the metabolomic and integrative analyses revealed that several antiaging-related metabolites and genes were associated with sternal catch-up growth. Pyrimidine metabolism was identified as the most significant pathway in which most differential metabolites (DMs) between FF and 60% FR were enriched. The results from integrative analysis revealed that the content and expression of 4-aminobutyric acid (GABA) and its related genes showed relatively higher activity in the 60% FR group than in the FF group. The present study identifies multiple biomarkers associated with duck sternum development that are influenced by feed restriction and suggests the potential mechanism of feed restriction-associated duck sternal catch-up growth.

Comparative transcriptome analysis reveals mechanisms of restriction feeding on lipid metabolism in ducks.

Presently, excessive fat deposition is the main reason to limit the development of duck industry. In the production, the methods of restricted feeding (RF) were widely used to reduce the lipid deposition of ducks. The liver (L), abdominal adipose (AA), and subcutaneous adipose (SA) were the main tissues of lipid metabolism and deposition of ducks. However, the mechanisms of lipid metabolism and deposition of ducks under RF have not been fully clarified. In this study, in order to better understand the mechanisms of lipid metabolism and deposition in ducks under RF, a total of 120 male Nonghua ducks were randomly divided into a free feeding group (FF, n = 60) and RF group (RF, n = 60), then comparative transcriptomic analysis of L, AA, and SA between FF (n = 3) and RF (n = 3) ducks was performed at 56 d of age. Phenotypically, L, AA, and SA index of FF group was higher than that in RF group. There were 279, 390, and 557 differentially expressed genes (DEGs) in L, AA, and SA. Functional enrichment analysis revealed that ECM-receptor interaction and metabolic pathways were significantly enriched in L, AA, and SA. Lipid metabolism-related pathways including fatty acid metabolism, unsaturated fatty acid synthesis, and steroidogenesis were significantly enriched in AA and SA. Moreover, through integrated analysis weighted gene coexpression network (WGCNA) and protein-protein interaction network, 10 potential candidate genes involved in the ECM-receptor interaction and lipid metabolism pathways were identified, including 3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2), aldolase B (ALDOB), formimidoyltransferase cyclodeaminase(FTCD), phosphoenolpyruvate carboxykinase 1 (PCK1), tyrosine aminotransferase (TAT), stearoyl-CoA desaturase (SCD), squalene epoxidase (SQLE), phosphodiesterase 4B (PDE4B), choline kinase A (CHKA), and elongation of very-long-chain fatty acids-like 2 (ELOVL2), which could play a key role in lipid metabolism and deposition of ducks under RF. Our study reveals that the liver might regulate the lipid metabolism of abdominal adipose and subcutaneous adipose through ECM-receptor interaction and metabolic pathways (fatty acid metabolism, unsaturated fatty acid synthesis, and steroid synthesis), thus to reduce the lipid deposition of ducks under RF. These results provide novel insights into the avian lipid metabolism and will help better understand the underlying molecular mechanisms.

Estrogen Receptor Gene 1 (ESR1) Mediates Lipid Metabolism in Goose Hierarchical Granulosa Cells Rather than in Pre-Hierarchical Granulosa Cells.

(1) Background: The role of estrogen receptor gene 1 (ESR1) in female reproduction and lipid metabolism has been extensively investigated. However, its contribution to lipid metabolism during the development of poultry follicles remains unclear. (2) Methods: This study aimed to explore the function of ESR1 via overexpressing (ESR1ov) and interfering (ESR1si) with its expression in pre-hierarchical granulosa cells (phGCs) and hierarchical granulosa cells (poGCs). (3) Results: We successfully cloned and obtained an 1866 bp segment of the full-length CDS region of the Sichuan white goose ESR1 gene. In phGCs of the ESR1ov and ESR1si groups, there were no significant changes compared to the control group. However, in poGCs, the ESR1ov group exhibited decreased lipid deposition, triglycerides, and cholesterol compared to the control group, while the ESR1si group showed increased lipid deposition, triglycerides, and cholesterol. The expression of APOB and PPARα was significantly reduced in the ESR1ov group compared to the ESR1ov-NC group. Moreover, significant changes in the expression of ACCα, DGAT1, SCD, CPT1, and ATGL were observed between the ESR1si and ESR1si-NC group. (4) Conclusions: These findings shed light on the function and molecular mechanism of ESR1 in lipid metabolism in goose poGCs, providing a better understanding of the physiological process of goose follicular development.

Exploring right ovary degeneration in duck and goose embryos by histology and transcriptome dynamics analysis.

BACKGROUND: The development of asymmetric chick gonads involves separate developmental programs in the left and right gonads. In contrast to the left ovary developing into a fully functional reproductive organ, the right ovary undergoes gradual degeneration. However, the molecular mechanisms underlying the the degeneration of the right ovary remain incompletely understood. In the present study, we investigated the histomorphological and transcriptomic changes in the right ovary of ducks and geese during the the embryonic stage up to post-hatching day 1. RESULT: Hematoxylin-eosin stainings revealed that the right ovary developed until embryonic day 20 in ducks (DE20) or embryonic day 22 in geese (GE22), after which it started to regress. Further RNA-seq analyses revealed that both the differentially expressed genes (DEGs) in ducks and geese right ovary developmental stage were significantly enriched in cell adhesion-related pathway (ECM-receptor interaction, Focal adhesion pathway) and Cellular senescence pathway. Then during the degeneration stage, the DEGs were primarily enriched in pathways associated with inflammation, including Herpes simplex virus 1 infection, Influenza A, and Toll-like receptor signaling pathway. Moreover, duck-specific DEGs showed enrichment in Steroid hormone biosynthesis, Base excision repair, and the Wnt signaling pathway, while geese-specifically DEGs were found to be enriched in apoptosis and inflammation-related pathways, such as Ferroptosis, Necroptosis, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway. These findings suggest that the degeneration process of the right ovary in ducks occurs at a slower pace compared to that in geese. Additionally, the observation of the left ovary of the geese varying degeneration rates in the right ovary after hatching indicated that the development of the left ovary may be influenced by the degeneration of the right ovary. CONCLUSION: The data presented in this study provide valuable insights into the dynamic changes in histological structure and transcriptome during the degeneration of the right ovary in ducks and geese. In addition, through the analysis of shared characteristics in the degeneration process of the right ovary in both ducks and geese, we have uncovered the patterns of degradation and elucidated the molecular mechanisms involved in the regression of the right ovary in poultry. Furthermore, we have also made initial discoveries regarding the relationship between the degeneration of the right ovary and the development of the left ovary.

Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis.

BACKGROUND: The genetic locus responsible for duck body size has been fully explained before, but the growth trait-related genetic basis is still waiting to be explored. For example, the genetic site related to growth rate, an important economic trait affecting marketing weight and feeding cost, is still unclear. Here, we performed genome wide association study (GWAS) to identify growth rate-associated genes and mutations. RESULT: In the current study, the body weight data of 358 ducks were recorded every 10 days from hatching to 120 days of age. According to the growth curve, we evaluated the relative and absolute growth rates (RGR and AGR) of 5 stages during the early rapid growth period. GWAS results for RGRs identified 31 significant SNPs on autosomes, and these SNPs were annotated by 24 protein-coding genes. Fourteen autosomal SNPs were significantly associated with AGRs. In addition, 4 shared significant SNPs were identified as having an association with both AGR and RGR, which were Chr2: 11483045 C>T, Chr2: 13750217 G>A, Chr2: 42508231 G>A and Chr2: 43644612 C>T. Among them, Chr2: 11483045 C>T, Chr2: 42508231 G>A, and Chr2: 43644612 C>T were annotated by ASAP1, LYN and CABYR, respectively. ASAP1 and LYN have already been proven to play roles in the growth and development of other species. In addition, we genotyped every duck using the most significant SNP (Chr2: 42508231 G>A) and compared the growth rate difference among each genotype population. The results showed that the growth rates of individuals carrying the Chr2: 42508231 A allele were significantly lower than those without this allele. Moreover, the results of the Mendelian randomization (MR) analysis supported the idea that the growth rate and birth weight had a causal effect on the adult body weight, with the growth rate having a greater effect size. CONCLUSION: In this study, 41 SNPs significantly related to growth rate were identified. In addition, we considered that the ASAP1 and LYN genes are essential candidate genes affecting the duck growth rate. The growth rate also showed the potential to be used as a reliable predictor of adult weight, providing a theoretical reference for preselection.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis.

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

Comparative transcriptome analysis identified crucial genes and pathways affecting sperm motility in the reproductive tract of drakes with different libido.

Libido can affect the semen quality of male, and the sperm motility in semen quality parameters is a reliable index to evaluate the fertility of male. In drakes, the sperm motility is gradually acquired in testis, epididymis, and spermaduct. However, the relationship between libido and sperm motility in drakes has not been reported and the mechanisms of testis, epididymis, and spermaduct regulating the sperm motility of drakes are unclear. Therefore, the purpose of the present study was to compare the semen quality of drakes with libido level 4 (LL4) and libido level 5 (LL5), and tried to identify the mechanisms regulating the sperm motility in drakes by performing RNA-seq in testis, epididymis, and spermaduct. Phenotypically, the sperm motility of drakes (P < 0.01), weight of testis (P < 0.05), and organ index of epididymis (P < 0.05) in the LL5 group were significantly better than those in LL4 group. Moreover, compared with the LL4 group, the ductal square of seminiferous tubule (ST) in testis was significantly bigger in the LL5 group (P < 0.05), and the seminiferous epithelial thickness (P < 0.01) of ST in testis and lumenal diameter (P < 0.05) of ductuli conjugentes/dutus epididymidis in epididymis were significantly longer in the LL5 group. In transcriptional regulation, in addition to KEGG pathways related to metabolism and oxidative phosphorylation, lots of KEGG pathways associated with immunity, proliferation, and signaling were also significantly enriched in testis, epididymis, and spermaduct, respectively. Furthermore, through the integrated analysis of coexpression network and protein-protein interaction network, 3 genes (including COL11A1, COL14A1, and C3AR1) involved in protein digestion and absorption pathway and Staphylococcus aureus infection pathway were identified in testis, 2 genes (including BUB1B and ESPL1) involved in cell cycle pathway were identified in epididymis, and 13 genes (including DNAH1, DNAH3, DNAH7, DNAH10, DNAH12, DNAI1, DNAI2, DNALI1, NTF3, ITGA1, TLR2, RELN, and PAK1) involved in Huntington disease pathway and PI3K-Akt signaling pathway were identified in spermaduct. These genes could play crucial roles in the sperm motility of drakes with different libido, and all data the present study obtained will provide new insights into the molecular mechanisms regulating sperm motility of drakes.

miR-202-5p Inhibits Lipid Metabolism and Steroidogenesis of Goose Hierarchical Granulosa Cells by Targeting ACSL3.

miRNAs are critical for steroidogenesis in granulosa cells (GCs) during ovarian follicular development. We have previously shown that miR-202-5p displays a stage-dependent expression pattern in GCs from goose follicles of different sizes, suggesting that this miRNA could be involved in the regulation of the functions of goose GCs; therefore, in this study, the effects of miR-202-5p on lipid metabolism and steroidogenesis in goose hierarchical follicular GCs (hGCs), as well as its mechanisms of action, were evaluated. Oil Red O staining and analyses of intracellular cholesterol and triglyceride contents showed that the overexpression of miR-202-5p significantly inhibited lipid deposition in hGCs; additionally, miR-202-5p significantly inhibited progesterone secretion in hGCs. A bioinformatics analysis and luciferase reporter assay indicated that Acyl-CoA synthetase long-chain family member 3 (ACSL3), which activates long-chain fatty acids for the synthesis of cellular lipids, is a potential target of miR-202-5p. ACSL3 silencing inhibited lipid deposition and estrogen secretion in hGCs. These data suggest that miR-202-5p exerts inhibitory effects on lipid deposition and steroidogenesis in goose hGCs by targeting the ACSL3 gene.

Role of SNPs located in the exon 9 of ATAPA1 gene on goose egg production.

The meat and egg of goose is one of the main components of human food supply. The improvement of goose egg production is particularly important for the increasing human population. However, limited information is available about the effective molecular markers and mechanisms of egg production in goose. In this study, we jointly utilized the data of genome resequencing in different egg production Sichuan white goose and transcriptome at different follicle development stages to identified the molecular markers and mechanisms of egg production. The coefficient of variation of individual egg production in Sichuan white goose population is 0.42 to 0.49. Fifty individuals with the highest (laying 365 days egg number, LEN365 = 79-145) and 50 individuals with the lowest (LEN365 = 8-48) egg production were divided into high and low egg production groups. Based on whole-genome sequencing data of the selected samples, 36 SNPs (annotation novel.12.470, CELF2, ATP1A1, KCNJ6, RAB4A, UST, REV3L, DHX15, CAVN2, SLC5A9, Cldn5, MRPS23, and Tspan2) associated with the LEN365 were identified, involving multiple pathways such as metabolism and endocrinology. Notably, 5 SNPs located in the exon9 of ATP1A1 were identified by GWAS analysis. The association analysis with LEN365 showed the phenotypic variance explained of this haplotype consisting of 5 SNPs is 20.51%. Through transcriptome data analysis, we found the expression of ATP1A1 in the granular layers was increased in the stage of small yellow follicle to large yellow follicle (LYF) and LYF to F5, while decreased in F2 to F1. For the first time, we report the haplotype region formed by 5 SNPS on exon9 of ATP1A1 is associated with egg production in goose and involved in follicle selection and maturation processes.

miR-27b-3p inhibits estrogen secretion of goose granulosa cells by targeting CYP1B1 through the AMPK signaling pathway.

Although miR-27b-3p has been evidenced to regulate the proliferation, apoptosis, and differentiation of a variety of mammalian cell types, its actions and mechanisms on ovarian cell steroidogenesis remains largely unknown in both mammalian and avian species. In this study, we aimed to determine the expression profiles of miR-27b-3p in granulosa cell layers during goose ovarian follicle development and to reveal its actions on estrogen (E2) secretion of goose granulosa cells as well as the underlying regulatory mechanisms. It was observed that miR-27b-3p was ubiquitously expressed throughout follicle development but exhibited much higher levels in hierarchical- than in prehierarchical follicles. In cultured granulosa cells from the fourth through second largest preovulatory (F4-F2) follicles of goose, up- and downregulation of miR-27b-3p by using its mimic and inhibitor significantly decreased and increased E2 secretion, respectively. Meanwhile, the mRNA levels of STAR and CYP19A1 were significantly reduced while those of CYP11A1 and 3ßHSD were elevated in the mimic-transfected granulosa cells. By comparison, downregulation of miR-27b-3p enhanced the mRNA levels of STAR but had no significant effects on those of CYP19A1, CYP11A1, and 3ßHSD. Results from bioinformatic prediction and luciferase reporter assay demonstrated that CYP1B1 was a downstream target of miR-27b-3p. Although the siRNA-mediated downregulation of CYP1B1 did not significantly change E2 secretion by goose granulosa cells, it reduced the mRNA levels of STAR and CYP19A1 as well as those of LKB1 and AMPKα, which are involved in the AMPK signaling pathway. Taken together, these data suggest that miR-27b-3p plays an inhibitory role in E2 secretion by goose F4-F2 granulosa cells, at least in part, by targeting CYP1B1 through the AMPK signaling pathway.

Comparative analysis of amino acid content and protein synthesis-related genes expression levels in breast muscle among different duck breeds/strains.

Evidences have found important effects of breeds/strains on the content of amino acids (AAs) which is an important substrate for protein synthesis and contributes greatly to meat quality. Therefore, the objective of the present study was to compare the AAs content and protein synthesis-related genes expression levels in breast muscle of native breed (Jianchang duck (J)), hybrid strains (BH1, BH2, and MC♂ × (BGF2♂ × GF2♀)♀ (MC)), and commercial breed (Cherry Verry duck). Results showed that a total of 17 AAs (TAAs) was detected from breast muscle among 5 duck breeds/strains including 11 essential AAs (EAAs). Among these AAs, the contents of Proline, Threonine, Glutamine, Serine, Methionine, Phenylalanine, Histidine, and Cysteine were significant difference among 5 duck breeds/strains. The contents of EAAs, TAAs, and flavor AAs were higher in breast muscle of J and BH2 than those in other duck breeds/strains, and the ratio of EAAs/TAAs was higher in breast muscle of BH2. Furthermore, the expression levels of eukaryotic translation initiation factor 4E-binding protein 1, mammalian target of rapamycin, and proton-coupled amino acid transporter 1 were the highest in breast muscle of BH2, and that of solute carrier family 38 member 2 was the highest in breast muscle of J. Meanwhile, principal component analysis results showed that principal component 1 of BH1, principal component 3 of BH2, and principal component 2 of MC were positively corelated with EAAs/TAAs, and principal component 1 was positively correlated with flavor AAs and EAAs. In conclusion, compared to BH1, MC, and Cherry Verry duck, AA content was higher in breast muscle of BH2 and J, which might be associated with the higher expression levels of mammalian target of rapamycin, eukaryotic translation initiation factor 4E-binding protein 1, and proton-coupled amino acid transporter 1 in breast muscle of BH2 and solute carrier family 38 member 2 in breast muscle of J. The comparative analysis of AA content in breast muscle among different duck breeds/strains could provide an important basis for improving the nutritional value of duck meat in the breeding process.

Transcriptomic analysis of the thyroid and ovarian stroma reveals key pathways and potential candidate genes associated with egg production in ducks.

The importance of thyroid-related genes has been repeatedly mentioned in the transcriptome studies of poultry with different laying performance, yet there are few systematic studies to unravel the regulatory mechanisms of the thyroid-ovary axis in the poultry egg production process. In this study, we compared the transcriptome profiles in the thyroid and ovarian stroma between high egg production (GP) and low egg production (DP) ducks, and then revealed the pathways and candidate genes involved in the process. We identified 1,114 and 733 differentially expressed genes (DEGs) in the thyroid and ovarian stroma, separately. The Gene Ontology (GO) analysis showed that a total of 504 and 189 GO terms were identified in the thyroid and ovarian stroma (P < 0.05). Three common GO terms were identified from the top 5 GO terms with the highest significant level in two tissues, including extracellular space, calcium ion binding, and integral component of plasma membrane. The enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that 15 and 14 KEGG pathways were significantly (P < 0.05) enriched in the thyroid and ovarian stroma, respectively. And, there were 8 common pathways, including neuroactive ligand-receptor interaction, calcium signaling pathway, ECM-receptor interaction, PPAR signaling pathway, melanogenesis, wnt signaling pathway, vascular smooth muscle contraction, and cytokine-cytokine receptor interaction. Notably, the neuroactive ligand-receptor interaction pathway was the most significantly enriched by the DEGs both in the thyroid and ovarian stroma. The interaction among DEGs enriched in the neuroactive ligand-receptor interaction and ECM-receptor interaction suggested that the thyroid may regulate ovarian development by these genes. Through integrated analysis of the protein-protein interaction (PPI) network and KEGG pathway maps, 9 key DEGs (PTH, THBS2, THBS4, CD36, ADIPOQ, ACSL6, PRKAA2, CRH, and PCK1) were identified, which could play crucial roles in the thyroid to regulate ovarian function and then affect egg-laying performance between GP and DP. This study serves as a basis to explore the molecular mechanism of the thyroid affecting ovarian function and egg production in female ducks and may help to identify molecular markers that can be used for duck genetic selection.

Integrated Transcriptome and Metabolome Analysis Reveals the Regulatory Mechanisms of FASN in Geese Granulosa Cells.

FASN plays a critical role in lipid metabolism, which is involved in regulating ovarian follicular development. However, the molecular mechanisms of how FASN regulate the function of ovarian follicular cells still remain elusive. In this study, by overexpression or interference of FASN in pre-hierarchical follicle granulosa cells (phGCs) and hierarchical follicle granulosa cells (hGCs), we analyzed their effects on the granulosa cell transcriptome and metabolome profiles using RNA-Seq and LC-MS/MS, respectively. The results showed that overexpression of FASN promoted proinflammatory factors expression by activating TLR3/IRF7 and TLR3/NF-κB pathways in phGCs, but only by activating TLR3/IRF7 pathways in hGCs. Then, necroptosis and apoptosis were triggered through the JAK/STAT1 pathway (induced by inflammatory factors) and BAK/caspase-7 pathway, respectively. The combined analysis of the metabolome and transcriptome revealed that FASN affected the demand of GCs for 5-hydroxytryptamine (5-HT) by activating the neuroactive ligand-receptor interaction pathway in two categorized GCs and only altering the metabolic pathway of tryptophan in phGCs, and ultimately participated in regulating the physiological function of geese GCs. Taken together, this study showed that the mechanisms of FASN regulating the physiological function of geese phGCs and hGCs were similar, but they also had some different characteristics.