Transcriptome Profiling Unveils Key Genes Regulating the Growth and Development of Yangzhou Goose Knob

Goose is one of the most economically valuable poultry species and has a distinct appearance due to its possession of a knob. A knob is a hallmark of sexual maturity in goose (Anser cygnoides) and plays crucial roles in artificial selection, health status, social signaling, and body temperature regulation. However, the genetic mechanisms influencing the growth and development of goose knobs remain completely unclear. In this study, histomorphological and transcriptomic analyses of goose knobs in D70, D120, and D300 Yangzhou geese revealed differential changes in tissue morphology during the growth and development of goose knobs and the key core genes that regulate goose knob traits. Observation of tissue sections revealed that as age increased, the thickness of the knob epidermis, cuticle, and spinous cells gradually decreased. Additionally, fat cells in the dermis and subcutaneous connective tissue transitioned from loose to dense. Transcriptome sequencing results, analyzed through differential expression, Weighted Gene Co-expression Network Analysis (WGCNA), and pattern expression analysis methods, showed D70-vs.-D120 (up-regulated: 192; down-regulated: 423), D70-vs.-D300 (up-regulated: 1394; down-regulated: 1893), and D120-vs.-D300 (up-regulated: 1017; down-regulated: 1324). A total of 6243 differentially expressed genes (DEGs) were identified, indicating varied expression levels across the three groups in the knob tissues of D70, D120, and D300 Yangzhou geese. These DEGs are significantly enriched in biological processes (BP) such as skin morphogenesis, the regulation of keratinocyte proliferation, and epidermal cell differentiation. Furthermore, they demonstrate enrichment in pathways related to goose knob development, including ECM-receptor interaction, NF-kappa B, and PPAR signaling. Through pattern expression analysis, three gene expression clusters related to goose knob traits were identified. The joint analysis of candidate genes associated with goose knob development and WGCNA led to the identification of key core genes influencing goose knob development. These core genes comprise WNT4, WNT10A, TCF7L2, GATA3, ADRA2A, CASP3, SFN, KDF1, ERRFI1, SPRY1, and EVPL. In summary, this study provides a reference for understanding the molecular mechanisms of goose knob growth and development and provides effective ideas and methods for the genetic improvement of goose knob traits.

Protein profile analysis of Jilin white goose testicles at different stages of the laying cycle by DIA strategy

BACKGROUND

Jilin white goose is an excellent local breed in China, with a high annual egg production and laying eggs mainly from February to July each year. The testis, as the only organ that can produce sperm, can affect the sexual maturity and fecundity of male animals. Its growth and development are affected and regulated by a variety of factors. Proteomics is generally applied to identify and quantify proteins in cells and tissues in order to understand the physiological or pathological changes that occur in tissues or cells under specific conditions. Currently, the female poultry reproductive system has been extensively studied, while few related studies focusing on the regulatory mechanism of the reproductive system of male poultry have been conducted.

RESULTS

A total of 1753 differentially expressed proteins (DEPs) were generated in which there were 594, 391 and 768 different proteins showing differential expression in three stages, Initial of Laying Cycle (ILC), Peak of Laying Cycle (PLC) and End of Laying Cycle (ELC). Furthermore, bioinformatics was used to analyze the DEPs. Gene ontology (GO) enrichment, Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein interaction (PPI) network analysis were adopted. All DEPs were found to be implicated in multiple biological processes and pathways associated with testicular development, such as renin secretion, Lysosomes, SNARE interactions in vesicle trafficking, the p53 signaling pathway and pathways related to metabolism. Additionally, the reliability of transcriptome results was verified by real-time quantitative PCR by selecting the transcript abundance of 6 selected DEPs at the three stages of the laying cycle.

CONCLUSIONS

The funding in this study will provide critical insight into the complex molecular mechanisms and breeding practices underlying the developmental characteristics of testicles in Jilin white goose.

miR-140-y targets TCF4 to regulate the Wnt signaling pathway and promote embryonic feather follicle development in Hungarian white goose

Goose down feather has become one of the most important economical products in the goose breeding industry and it provides several essential physiological roles in birds. Therefore, understanding and regulating the development of skin and feather follicles during embryogenesis is critical for avian biology and the poultry industry. MicroRNAs are known to play an important role in controlling gene expression during skin and feather follicle development. In this study, bioinformatics analysis was conducted to select miR-140-y as a potential miRNA involved in skin and feather follicle development and to predict TCF4 as its target gene. This gene was expressed at significant levels during embryonic feather follicle development, as identified by qPCR and Western blot. The targeting relationship was confirmed by a dual-luciferase assay in 293T cells. Then, the miR-140-y/TCF4 function in dermal fibroblast cells was explored. The results showed that miR-140-y could suppress the proliferation of goose embryonic dermal fibroblast cells (GEDFs) by suppressing the activity of some Wingless-types (Wnt) pathway related genes and proliferation marker genes, while miR-140-y inhibition led to the opposite effect. Similarly, the inhibition of the TCF4 gene results in blocking the proliferation of GEDFs by reducing the activity of some Wnt pathway-related genes. Finally, the co-transfection of miR-140-y inhibitor and siRNA-TCF4 results in a rescue of the TCF4 function and an increase of the Wnt signaling pathway and GEDFs proliferation. In conclusion, these results demonstrated that the miR-140-y-TCF4 axis influences the activity of the Wnt signaling pathway and works as a dynamic regulator during skin and feather follicle development.

Phenotypic and morphometric characterization of domestic geese raised in northern Benin

Documentation on the domestic geese (Anser anser) in Benin is scarce, making it objectively difficult to exploit. Its production depends on small flocks raised by livestock farmers in different areas. The aim of the study was to describe the phenotypic and morphometric characteristics of geese encountered in northern Benin. To this end, a total of 576 adult geese (353 males and 223 females) from 102 farms in 4 agro-ecological zones of northern Benin were evaluated for 11 quantitative and 5 qualitative body traits. There are 6 plumage colors with white (42.01%), white-brown magpie (24.65%) and white-grey magpie (17.19%) as dominant colors. The proportions of white-grey (75%) and multicolored (51.85%) in Far northern zone of Benin (FNZB) were higher (P < 0.05) than those in the cotton zone of northern Benin (CZNB), Food-producing zone of southern Borgou (FZSB) and West-Atacora Zone (WAZ). The red coloring of the tarsi was more dominant in the FZSB (63.09%) and the FZSB (61.79%) (P < 0.05). Orange (57.47%) and yellow (28.82%) colored beaks were dominant but did not vary (P > 0.05) from one agro-ecological zone to another. Quantitative traits such as body length, tarsus length, wing span and thoracic circumference of geese varied (P < 0.05) between 71.34 to73.22 cm, 10.08 to 10.6 cm, 131.95 to 135.42 cm and between 42.07 to 43.86 cm respectively. Males differed significantly from females (P < 0.05) for all morphometric traits. The live weight of geese in the FNZB showed higher values than those of other agro-ecological zones (P ˂ 0.05). In addition, white phenotype geese (3.76 kg) were heavier (P ˂ 0.05) than other phenotypes. All correlations between weight and body measurements of domestic geese were positive, but the correlations between live weight and wing span (r = 0.68) were the strongest. These correlations could be used to assess the live weight of the geese population studied and for selection based on live weight. This study provides a reference for morpho-biometric traits and will be complemented by molecular characterization.

The development rule of feathers and application of hair root tissue in sex identification of Yangzhou geese

Accurate gender identification is crucial for the study of bird reproduction and evolution. The current study aimed to explore and evaluate the effectiveness of a noninvasive method for gender identification in Yangzhou geese. In this experiment, 600 goose eggs were collected. Hair root tissues were used for PCR amplification, molecular sequencing, and anal inversion for early sex recognition in goslings. According to the DNA amplification results for the feather pulp tissue of 2-wk-old geese, bands appeared at 436 bp (CHD1-Z) and 330 bp (CHD1-W) upon gel electrophoresis. This method considered the base of goose feathers to accelerate the process of gender recognition. By examining the sex of anatomized poultry for verification, the accuracy rate of PCR gel electrophoresis and molecular sequencing sex identification was 100%, whereas the average accuracy rate of anal inversion was 97.41%. In the comparison of feather growth trends at 0 to 18 wk of age, the feather root weight (FRW), feather root length (FRL), feather branch length (FBL), and feather shaft diameter (FSD) of Yangzhou goose of the same age were not significantly different between males and females (P > 0.05). At 6 wk of age, the FRW, FRL, and FSD in males and FRL in females increased rapidly; their growth increased by 84.43, 67.58, 45.10, and 69.42%, respectively. At 10 wk of age, the male FRL, male FBL, and female FBL increased by 37.31, 34.81, and 21.72, respectively. The Boltzmann model was found to be the best-fitting model for the feathers of male Yangzhou geese. Early sex identification based on feather growth trends between the sexes is not feasible. This study provides a convenient and reliable technical means for early sex identification of waterfowl and serves as an ecological strategy for protecting the reproduction of poultry populations.

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.

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.

Identification, pathological, and genomic characterization of novel goose reovirus associated with liver necrosis in geese, China

Since 2021, a novel strain of goose reovirus (GRV) has emerged within the goose farming industry in Guangdong province, China. This particular viral variant is distinguished by the presence of white necrotic foci primarily localized in the liver and spleen, leading to substantial economic losses for the poultry industry. However, the etiology, prevalence and genomic characteristics of the causative agent have not been thoroughly investigated. In this study, we conducted an epidemiological inquiry employing suspected GRV samples collected from May 2021 to September 2022. The macroscopic pathological and histopathological lesions associated with GRV-infected clinical specimens were examined. Moreover, we successfully isolated the GRV strain and elucidated the complete genome sequence of the isolate GD21/88. Through phylogenetic and recombination analysis, we unveiled that the GRV strains represent a novel variant resulting from multiple reassortment events. Specifically, the μNS, λC, and σNS genes of GRV were found to have originated from chicken reovirus, while the σA gene of GRV exhibited a higher degree of similarity with a novel duck reovirus. The remaining genes of GRV were traced back to Muscovy duck reovirus. Collectively, our findings underscore the significance of GRV as a pathogenic agent impacting the goose farming industry. The insights gleaned from this study contribute to a more comprehensive understanding of the epidemiology of GRV in Southern China and shed light on the genetic reassortment events exhibited by the virus.

Integrated analysis of whole genome and transcriptome sequencing uncovers genetic differences between Zi goose and Xianghai flying goose

Zi goose is a famous indigenous breed originating from northeast China with high annual egg production. Xianghai flying goose is a composite breed and is bred by crosses of the wild swan goose and the Zi goose. Our previous study revealed significant differences in muscle fiber characteristics between the two populations. Here, we aimed to reveal the underlying genetic basis of the above phenotype differences through whole-genome and transcriptome analysis. A total of 20 blood samples (10 Zi geese and 10 Xianghai flying geese) were used for whole genome sequencing, and eight breast muscle tissue samples (four Zi geese and four Xianghai flying geese) were used for RNA sequencing. Using the FST and XP-EHH analysis, some highly differentiated genome regions annotated with egg production (RORB, WNT4, BMPR1B) and breast muscle development (WNT7B) between the two populations were detected. RNA-sequencing analysis revealed differentially expressed genes related to muscle development (IGF1, PAX7). Moreover, several genes were detected by both genome and transcriptome analysis, and some of them were reported to be associated with muscle growth (SLIT2, PREX1) and intramuscular fat (COL6A1). These findings will help researchers better understand the genetic basis related to egg production and muscle development in geese.

Transcriptome and lipidome integration unveils mechanisms of fatty liver formation in Shitou geese

Geese evolved from migratory birds, and when they consume excessive high-energy feed, glucose is converted into triglycerides. A large amount of triglyceride deposition can induce incomplete oxidation of fatty acids, leading to lipid accumulation in the liver and the subsequent formation of fatty liver. In the Chaoshan region of Guangdong, China, Shitou geese develop a unique form of fatty liver through 24 h overfeeding of brown rice. To investigate the mechanisms underlying the formation of fatty liver in Shitou geese, we collected liver samples from normally fed and overfed geese. The results showed that the liver size in the treatment group was significantly larger, weighing 3.5 times more than that in the control group. Extensive infiltration of lipid droplets was observed in the liver upon staining of tissue sections. Biochemical analysis revealed that compared to the control group, the treatment group showed significantly elevated levels of total cholesterol (T-CHO), triglycerides (TG), and glycogen in the liver. However, no significant differences were observed in the levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST), which are common indicators of liver damage. Furthermore, we performed a combined transcriptomic and lipidomic analysis of the liver samples and identified 1,510 differentially expressed genes (DEGs) and 1,559 significantly differentially abundant metabolites (SDMs). The enrichment analysis of the DEGs revealed their enrichment in metabolic pathways, cellular process-related signaling pathways, and specific lipid metabolism pathways. We also conducted KEGG enrichment analysis of the SDMs and compared them with the enriched signaling pathways obtained from the DEGs. In this study, we identified 3 key signaling pathways involved in the formation of fatty liver in Shitou geese, namely, the biosynthesis of unsaturated fatty acids, glycerol lipid metabolism, and glycerophospholipid metabolism. In these pathways, genes such as glycerol-3-phosphate acyltransferase, mitochondrial (GPAM), 1-acylglycerol-3-phosphate O-acyltransferase 2 (AGPAT2), diacylglycerol O-acyltransferase 2 (DGAT2), lipase, endothelial (LIPG), lipoprotein lipase (LPL), phospholipase D family member 4 (PLD4), and phospholipase A2 group IVF (PLA2G4F) may regulate the synthesis of metabolites, including triacylglycerol (TG), phosphatidate (PA), 1,2-diglyceride (DG), phosphatidylethanolamine (PE), and phosphatidylcholine (PC). These genes and metabolites may play a predominant role in the development of fatty liver, ultimately promoting the accumulation of TG in the liver and leading to the progression of fatty liver.

Parental betaine supplementation promotes gosling growth with epigenetic modulation of IGF gene family in the liver

Betaine is widely used as a feed additive in the chicken industry to promote laying performance and growth performance, yet it is unknown whether betaine can be used in geese to improve the laying performance of goose breeders and the growth traits of offspring goslings. In this study, laying goose breeders at 39 wk of age were fed basal (Control, CON) or betaine-supplemented diets at low (2.5 g/kg, LBT) or high (5 g/kg, HBT) levels for 7 wk, and the breeder eggs laid in the last week were collected for incubation. Offspring goslings were examined at 35 and 63 d of age. The laying rate tended to be increased (P = 0.065), and the feed efficiency of the breeders was improved by betaine supplementation, while the average daily gain of the offspring goslings was significantly increased (P < 0.05). Concentrations of insulin-like growth factor 2 (IGF-2) in serum and liver were significantly increased in the HBT group (P < 0.05), with age-dependent alterations of serum T3 levels. Concurrently, hepatic mRNA expression of the IGF gene family was significantly increased in goslings derived from betaine-treated breeders (P < 0.05). A higher ratio of proliferating cell nuclear antigen (PCNA)-immunopositive nuclei was found in the liver sections of the HBT group, which was confirmed by significantly upregulated hepatic expression of PCNA mRNA and protein (P < 0.05). Moreover, hepatic expression of thyroxine deiodinase type 1 (Dio1) and thyroid hormone receptor β (TRβ) was also significantly upregulated in goslings of the HBT group (P < 0.05). These changes were associated with significantly higher levels of global DNA 5-mC methylation, together with increased expression of methyl transfer genes (P < 0.05), including betaine-homocysteine methyltransferase (BHMT), glycine N-methyltransferase (GNMT), and DNA (cytosine-5-)-methyltransferase 1 (DNMT1). The promoter regions of IGF-2 genes, as well as the predicted TRβ binding site on the IGF-2 gene, were significantly hypomethylated (P < 0.05). These results indicate that gosling growth can be improved by dietary betaine supplementation in goose breeders via epigenetic modulation of the IGF gene family, especially IGF-2, in the liver.

Transferability of Nikita and Sukkula retrotransposons in domestic goose (Anser anser domesticus) genome

This article aimed to detect the existence of barley-specific Nikita and Sukkula retrotransposons in domestic geese samples and to evaluate the evolutionary relationships between these and other transposons belonging to the family Anatidae. Inter-retrotransposonamplified polymorphism-polymerase chain reaction (IRAP-PCR) method was performed for these retrotransposons movements in three diverse domestic goose populations (Chinese x Embden crossbred, Turkish White, and Turkish Multicolor). Polymorphism ratios were between 0 and 33% in all samples for Nikita and 0-73% in all samples for Sukkula. In addition, intrapopulation genetic polymorphism rates were also 0-15% in Chinese x Embden crossbred, 0-25% in Turkish White, 0-25% in Turkish Multicolor for Nikita; while 0-27% in Chinese x Embden, and 0-50% in Turkish Multicolor for Sukkula. There was no polymorphism for Sukkula among Turkish White samples. Moreover, the neighbour-joining method was used for phylogenetic tree construction using 38 sequences of different ducks, geese, and swans. In silico analyses supported the transitions of retrotransposons in the family Anatidae. It is concluded that transposon mobility among the phylogenetically distant species may lead to understanding evolutionary relationships. This report is one of the first studies investigating retrotransposon movements in domestic geese, revealing a new perspective on the goose genome regarding mobile genetic elements.

Effects of essential oil mixtures on expression of genes involved in lipogenesis and fatty acid oxidation in geese (Anser anser)

The use of essential oils has recently increased in the poultry sector. The aim of this study was to investigate the effects of essential oil mixture (juniper, mint, oregano and rosemary oil) on fatty acid oxidation and lipogenic gene expression in geese. Research groups were formed as C (control; no additives), EK1 (0.4 ml/l essential oil mixture supplemented) and EK2 (0.8 ml/l essential oil mixture supplemented). Relative expression levels of genes included in lipogenesis (ACCα, ChREBP, FASN, LXRα and SREBP-1) expression levels of genes included in fatty acid oxidation (ACOX1, CPT1, CPT1A, PPARα and PPARγ) were measured using RT-qPCR. Group EK1 upregulates the mRNA expression levels of genes involved in lipogenesis such as ACCα, ChREBP and SREBP-1, while it downregulates the mRNA expression in levels of all genes involved in fatty acid oxidation. Group EK2 increases the mRNA expression levels of genes involved in lipogenesis such as ACCα, FASN and SREBP-1, while it decreased mRNA expression at the levels of all genes involved in fatty acid oxidation, as in the other group. In the study, adding an essential oil mixture to drinking water is predicted to increase fatty liver because it upregulates genes related to fat synthesis (lipogenesis) and downregulates genes related to fat degradation (fatty acid oxidation).

Genomic signatures reveal selection in Chinese and European domesticated geese

The Swan goose and Greylag goose are species of geese native to East Asia and Europe, respectively, and are widely believed to be the ancestors of Chinese and European domesticated geese. The Yili goose (YL) and European domestic geese originated from the Greylag goose, but the history of domestication is unclear. In this study, we sequenced and analyzed the genome of the YL goose and the Hortobagy goose to combine with other previously sequenced goose populations for in-depth analysis. The population genetic variations in Stone geese, East Zhejiang White Geese, Taihu geese and Zi geese were also identified and compared. The results showed that admixture gene flow existed in the YL geese population, which was introgressed by Chinese geese, suggesting that gene flow events were frequent and widespread among domesticated geese. Further selected sweep analysis identified candidate genes and metabolic pathways that may be related to the differences in morphology. Several genes such as TGFBR3L, CMYA5, FOXD1, ARHGEF28 and SUCLG2 are associated with growth, reproduction and fertility traits. The results of this study will help to understand the genetic characteristics of domestic geese and the genes affecting important traits and provide a basis for the improved breed of domestic geese.

Expression profile and the G63A mutation of IGF2 gene associated with growth traits in Zhedong-White goose

Insulin-like growth factors 2 (IGF2) is an insulin-like growth factor that plays a major role in animal growth, cell proliferation and differentiation, as well as reproduction. IGF2 is well-known to be a candidate gene of growth and reproductive traits in many agricultural animals. Our previous study revealed that the G63A (Chr2: G26541617A) mutation within IGF2 exon 1 showed a significant association with egg numbers (E180d) of Sanhua goose population (p < 0.05). However, little work focus on the correlation between the IGF2 mutations and goose growth traits. In this study, qPCR indicated that the IGF2 mRNA highly expressed in leg muscle, liver, ovary and pituitary gland. Meanwhile, association analysis showed that the G63A mutation was significantly associated with the body weight of first-hatched Zhedong-White geese (BW0, p < 0.05), and strongly significantly associated with the BW2, BW4, BW6, BW8 and BW10 (p < 0.01). The GG homozygous had the lowest BW (from 4 weeks to 10 weeks) than those of AA and AG genotypes (p < 0.01), and the allele A was also positively correlated with the BW of the Zhedong-White goose population. Therefore, the G63A mutation in IGF2 may be an important genetic marker for goose breeding.

Circadian clock gene BMAL1 regulates STAR expression in goose ovarian preovulatory granulosa cells

The ovarian circadian clock plays a regulatory role in the avian ovulation-oviposition cycle. However, little is known regarding the ovarian circadian clock of geese. In this study, we investigated rhythmic changes in clock genes over a 48-h period and identified potential clock-controlled genes involved in progesterone synthesis in goose ovarian preovulatory granulosa cells. The results showed that BMAL1, CRY1, and CRY2, as well as 4 genes (LHR, STAR, CYP11A1, and HSD3B) involved in progesterone synthesis exhibited rhythmic expression patterns in goose ovarian preovulatory granulosa cells over a 48-h period. Knockdown of BMAL1 decreased the progesterone concentration and downregulated STAR mRNA and protein levels in goose ovarian preovulatory granulosa cells. Overexpression of BMAL1 increased the progesterone concentration and upregulated the STAR mRNA level in goose ovarian preovulatory granulosa cells. Moreover, we demonstrated that the BMAL1/CLOCK complex activated the transcription of goose STAR gene by binding to an E-box motif. These results suggest that the circadian clock is involved in the regulation of progesterone synthesis in goose ovarian preovulatory granulosa cells by orchestrating the transcription of steroidogenesis-related genes.

Base composition, adaptation, and evolution of goose astroviruses: codon-based investigation

Goose astroviruses (GoAstVs) are causative agents that account for fatal infection of goslings characterized by visceral urate deposition, resulting in severe economic losses in major goose-producing regions in China since 2017. In this study, we sought to unravel the intrinsic properties associated with adaptation and evolution in the host environment of GoAstVs. Consistent results from phylogenetic analysis and correspondence analysis performed on the codon usage patterns (CUPs) reveal 2 clusters of GoAstVs, namely, GoAstV-1 and GoAstV-2. However, multiple similar compositional characteristics were found, despite the high divergence between GoAstV-1 and GoAstV-2. Studies on the base composition of GoAstVs reveal an A/U bias, indicating a compositional constraint, while natural selection prevailed in determining the CUPs in the virus genome based on our neutrality plot analysis, reflecting high adaptive pressure to fit the host environment. Codon adaptation index (CAI) analysis revealed a higher degree of fitness to the CUPs of the corresponding host for GoAstVs than avian influenza virus and betacoronaviruses, which may be a favorable factor contributing to the high pathogenicity and wide distribution of GoAstVs in goslings. In addition, GoAstVs were less adapted to ducks and chickens, with significantly lower CAI values than to geese, which may be a reason for the different prevalence of GoAstVs among these species. Extensive investigations on dinucleotide distribution revealed a significant suppression of the CpG and UpA motifs in the virus genome, which may facilitate adaptation to the host's innate immune system by evading surveillance. In addition, our study reported the trends of increasing fitness to the host's microenvironment for GoAstVs through increasing adaptation to host CUPs and ongoing reduction of CpG motifs in the virus genome. The present analysis deepens our understanding of the basic biology, pathogenesis, adaptation and evolutionary pattern of GoAstVs, and contributes to the development of novel antiviral strategies.

Association between Yili goose sperm motility and expression profiles of mRNA and miRNA in testis

BACKGROUND

The study was conducted to find out the candidate microRNA (miRNA) and genes that associated with sperm motility of Yili goose through small RNA sequencing of testicular tissue of Yili goose, and provide a theoretical basis for the study of the regulation mechanism of sperm motility of Yili goose gander.

RESULTS

In this study, five male geese with high sperm motility and five male geese with low sperm motility were slaughtered to obtain their testis tissues for small RNA sequencing, and biological information methods were used for data analysis. The results showed that a total of 1575 known miRNAs and 68 novel miRNAs were identified in the testis tissue of Yili goose, and 71 differentially expressed miRNAs and 660 differentially expressed genes were screened. GO functional analysis showed that miRNAs target genes were mainly involved in the binding, kinase activity, structural constituent of cytoskeleton and intermediate filament cytoskeleton. KEGG functional analysis showed that miRNAs target genes were significantly enriched in arginine and proline metabolism, glycolysis / gluconeogenesis, fructose and mannose metabolism and beta-Alanine metabolism and other pathways. miRNAs-mRNAs interaction network suggests miR-140/miR-140-3p/miR-140-3p-NKAIN3, let-7d-BTG1 and miR-145-5p/miR -145a-5p-CLEC2E may play an important role in testis development and spermatogenesis.

CONCLUSIONS

The results of this study suggest that the sperm motility of Yili goose may be regulated by different miRNAs, and the target genes are significantly enriched in pathways related to sperm metabolism, indicating that miRNAs affect the sperm motility of Yili goose by regulating the metabolic process of sperm and the expression of related genes. This study can provide a reference for revealing the regulation mechanism of Yili goose sperm motility at the molecular level.

FGF12 regulates cell cycle gene expression and promotes follicular granulosa cell proliferation through ERK phosphorylation in geese

The granulosa cells play an important role in the fate of follicular development or atresia in poultry. Fibroblast growth factor 12 (FGF12) is downregulated in atretic follicles and may be involved in regulating granulosa cell survival in previous studies, but its molecular mechanism remains unclear. In this study, FGF12 overexpression and knockdown models of goose granulosa cells were constructed to investigate its function. The downstream expression of the cell cycle pathway was analyzed by qPCR. Granulosa cell proliferative activity and apoptosis were detected by CCK8 and TUNEL. Protein phosphorylation levels of ERK and AKT were measured using Western blotting to analyze the key pathway of FGF12 regulation of granulosa cell proliferation. ERK protein phosphorylation inhibitor was added for further verification. After overexpression of FGF12, cell proliferation activity was increased, the expressions of cell cycle pathway genes CCND1, CCNA2, MAD2, and CHK1 were upregulated, the apoptosis of granulosa cell was decreased, and Caspase 3 gene and protein expression were downregulated. After the knockdown of FGF12, cell proliferation activity decreased, the expression of downstream genes in the cell cycle pathway was downregulated, the apoptosis of granulosa cells was increased, and the Bcl-2 gene and protein were downregulated. Overexpression of FGF12 promoted the synthesis of P4 and upregulates the expression of the STAR gene. Overexpression of FGF12 promoted ERK protein phosphorylation but did not affect AKT phosphorylation. The addition of ERK phosphorylation inhibitors resulted in the elimination of the increase in cell proliferative activity caused by FGF12 overexpression. In conclusion, FGF12 could promote proliferation and inhibit apoptosis of goose granulosa cells by increasing ERK phosphorylation.

Significant association of mutations close to LCORL gene with growth performance in Zhedong white geese (Anser cygnoides)

This study aimed to investigate the role of the LCORL gene in regulating the growth performance of Zhedong white (ZDW) geese, belonging to the swan geese (Anser cygnoides), and identify possible selective signatures in diverse goose breeds. Single nucleotide polymorphisms around LCORL were genotyped, and their associations with body-size-related (BSR) traits were estimated. The results showed that the genotyped loci upstream of LCORL were significantly related to the body weight and breast width of ZDW geese aged 10 weeks (p < 0.05). A genome scan comparing expected heterozygosity among different breeds identified a ~150 kb long genomic region with extremely low heterozygosity downstream of LCORL among swan geese. Further, significant associations of variants within the low heterozygosity region among ZDW geese with BSR traits, including body weight, body length and breast width (p < 0.05) were also detected. Overall, mutations adjacent to LCORL were related to the growth performance of swan geese, and the significant effects of variants in a low-heterozygosity region on BSR traits provided valuable insights into the molecular mechanism of artificial selection reshaping body stature in swan geese.

Integrating genomics and transcriptomics to identify candidate genes for high egg production in Wulong geese (Anser cygnoides orientalis)

BACKGROUND

Wulong geese (Anser cygnoides orientalis) are known for their excellent egg-laying performance. However, they show considerable population differences in egg-laying behavior. This study combined genome-wide selection signal analysis with transcriptome analysis (RNA-seq) to identify the genes related to high egg production in Wulong geese.

RESULTS

A total of 132 selected genomic regions were screened using genome-wide selection signal analysis, and 130 genes related to high egg production were annotated in these regions. These selected genes were enriched in pathways related to egg production, including oocyte meiosis, the estrogen signaling pathway, the oxytocin signaling pathway, and progesterone-mediated oocyte maturation. Furthermore, a total of 890 differentially expressed genes (DEGs), including 340 up-regulated and 550 down-regulated genes, were identified by RNA-seq. Two genes - GCG and FAP - were common to the list of selected genes and DEGs. A non-synonymous single nucleotide polymorphism was identified in an exon of FAP.

CONCLUSIONS

Based on genome-wide selection signal analysis and transcriptome data, GCG and FAP were identified as candidate genes associated with high egg production in Wulong geese. These findings could promote the breeding of Wulong geese with high egg production abilities and provide a theoretical basis for exploring the mechanisms of reproductive regulation in poultry.

Genetic Diversity Analysis and Breeding of Geese Based on the Mitochondrial ND6 Gene

To explore the differences in body-weight traits of five goose breeds and analyze their genetic diversity and historical dynamics, we collected body-weight data statistics and used Sanger sequencing to determine the mitochondrial DNA of 100 samples of five typical goose breeds in China and abroad. The results indicated that Lion-Head, Hortobagy, and Yangzhou geese have great breeding potential for body weight. Thirteen polymorphic sites were detected in the corrected 505 bp sequence of the mitochondrial DNA (mtDNA) ND6 gene, accounting for approximately 2.57% of the total number of sites. The guanine-cytosine (GC) content (51.7%) of the whole sequence was higher than the adenine-thymine (AT) content (48.3%), showing a certain GC base preference. There were 11 haplotypes among the five breeds, including one shared haplotype. We analyzed the differences in the distribution of base mismatches among the five breeds and conducted Tajima's D and Fu's Fs neutral tests on the historical dynamics of the populations. The distribution of the mismatch difference presented an unsmooth single peak and the Tajima's D value of the neutral test was negative (D < 0) and reached a significant level, which proves that the population of the three species had expanded; the Lion-Head goose population tends to be stable. The genetic diversity of Lion-Head, Zhedong White, Yangzhou, and Taihu geese was equal to the average diversity of Chinese goose breeds. The Hortobagy goose is a foreign breed with differences in mating line breeding and hybrid advantage utilization.

Phosphoproteome Reveals Extracellular Regulated Protein Kinase Phosphorylation Mediated by Mitogen-Activated Protein Kinase Kinase-Regulating Granulosa Cell Apoptosis in Broody Geese

Geese have strong brooding abilities, which severely affect their egg-laying performance. Phosphorylation is widely involved in regulating reproductive activities, but its role in goose brooding behavior is unclear. In this study, we investigated differences in the phosphoprotein composition of ovarian tissue between laying and brooding geese. Brooding geese exhibited ovarian and follicular atrophy, as well as significant oxidative stress and granulosa cell apoptosis. We identified 578 highly phosphorylated proteins and 281 lowly phosphorylated proteins, and a KEGG pathway analysis showed that these differentially phosphorylated proteins were mainly involved in cell apoptosis, adhesion junctions, and other signaling pathways related to goose brooding behavior. The extracellular regulated protein kinase (ERK)-B-Cell Lymphoma 2(BCL2) signaling pathway was identified as playing an important role in regulating cell apoptosis. The phosphorylation levels of ERK proteins were significantly lower in brooding geese than in laying geese, and the expression of mitogen-activated protein kinase kinase (MEK) was downregulated. Overexpression of MEK led to a significant increase in ERK phosphorylation and BCL2 transcription in H2O2-induced granulosa cells (p < 0.05), partially rescuing cell death. Conversely, granulosa cells receiving MEK siRNA exhibited the opposite trend. In conclusion, geese experience significant oxidative stress and granulosa cell apoptosis during brooding, with downregulated MEK expression, decreased phosphorylation of ERK protein, and inhibited expression of BCL2.

The study of the strength and significance of four biological parameters on the body weight of goose

Alternative products such as those from high-value protien animals have increased the demand for the production of high-quality chicken meat in past few years. This study examines the impact of two distinct feeding types on goose body-weight, as well as the genetic variation of growth hormone (GH) and pituitary-specific transcription factor (Pit-1) genes in ten goose populations using single nucleotide polymorphism (SNP) markers and PCR-RFLP analysis. Both genes were seen as very important for productivity, especially in light of the COVID-19 and its effect on poultry industry at the time. The findings suggest that employing genetic indicators in these two genes in conjunction with a high-fat diet may be a feasible strategy for goose selection programme aiming to increase marketing body weight, as the high-fat diet outperformed the balanced diet. The study investigates the effect of gender, 2 types of diets, breeds and the genetic variation of the two genes, four SNPs were reported to be found: two at the GH gene exons C123T and C158T, and two at the Pit-1 gene exons G161A and T282G. Certain genotypes were found to have a substantial effect on the marketing body-weight of goose, which varied depending on the tested breeds. However, in terms of gender, males report higher and better performance levels than females. Diet, breeds and genotype interaction, and breeds, gender and genotype interaction were found to have a minor effect on goose body weight. However, diet, breeds, gender, SNP locus, diet and breeds interaction, and breeds and gender interaction were found to have a significant effect on goose body weight, as indicated by the effect size results.

Transcriptomic and Chromatin Landscape Analysis Reveals That Involvement of Pituitary Level Transcription Factors Modulate Incubation Behaviors of Magang Geese

The incubation behavior of geese seriously affects their egg production performance. Studies on incubation behavior have identified functional genes, but the regulatory architecture relationship between functional genes and chromatin accessibility remains poorly understood. Here, we present an integrated analysis of open chromatin profiles and transcriptome to identify the cis-regulatory element and their potential transcription factors involved in regulating incubation behavior in goose pituitary. Assay for transposase-accessible chromatin sequencing (ATAC-seq) revealed that open chromatin regions increased in the pituitary during the transition from incubation behavior to laying. We identified 920 significant differential accessible regions (DARs) in the pituitary. Compared to the laying stage, most DARs had higher chromatin accessibility in the brooding stage. Motif analysis of open DARs showed that the most significant transcription factor (TF) occupied sites predominantly enriched in motifs binding to the RFX family (RFX5, RFX2, and RFX1). While the majority of TF motifs enriched under sites of the nuclear receptor (NR) family (ARE, GRE, and PGR) in closed DARs at the incubation behavior stage. Footprint analysis indicated that the transcription factor RFX family exhibited higher binding on chromatin at the brooding stage. To further elucidate the effect of changes in chromatin accessibility on gene expression levels, a comparison of the transcriptome revealed 279 differentially expressed genes (DEGs). The transcriptome changes were associated with processes of steroid biosynthesis. By integrating ATAC-seq and RNA-seq, few DARs directly affect incubation behavior by regulating the transcription levels of genes. Five DAR-related DEGs were found to be closely related to maintaining the incubation behavior in geese. Footprinting analysis revealed a set of transcription factors (RFX1, RFX2, RFX3, RFX5, BHLHA15, SIX1, and DUX) which displayed the highest activity at the brooding stage. SREBF2 was predicted to be the unique differentially expressed transcription factor whose mRNA level was down-regulated and enriched in hyper-accessible regions of PRL in the broody stage. In the present study, we comprehensively profiled the transcriptome and chromatin accessibility in the pituitary related to incubation behavior. Our findings provided insight into the identification and analysis of regulatory elements in goose incubation behavior. The epigenetic alterations profiled here can help decipher the epigenetic mechanisms that contribute to the regulation of incubation behavior in birds.

Long non-coding RNA MSTRG.5970.28 regulates proliferation and apoptosis of goose follicle granulosa cells via the miR-133a-3p/ANOS1 pathway

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.

Whole Genome Resequencing Identifies Single-Nucleotide Polymorphism Markers of Growth and Reproduction Traits in Zhedong and Zi Crossbred Geese

The broodiness traits of domestic geese are a bottleneck that prevents the rapid development of the goose industry. To reduce the broodiness of the Zhedong goose and thus improve it, this study hybridized it with the Zi goose, which has almost no broody behavior. Genome resequencing was performed for the purebred Zhedong goose, as well as the F2 and F3 hybrids. The results showed that the F1 hybrids displayed significant heterosis in growth traits, and their body weight was significantly greater than those of the other groups. The F2 hybrids showed significant heterosis in egg-laying traits, and the number of eggs laid was significantly greater than those of the other groups. A total of 7,979,421 single-nucleotide polymorphisms (SNPs) were obtained, and three SNPs were screened. Molecular docking results showed that SNP11 located in the gene NUDT9 altered the structure and affinity of the binding pocket. The results suggested that SNP11 is an SNP related to goose broodiness. In the future, we will use the cage breeding method to sample the same half-sib families to accurately identify SNP markers of growth and reproductive traits.

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 (E₂) 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 E₂ 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 E₂ 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 E₂ secretion by goose F4-F2 granulosa cells, at least in part, by targeting CYP1B1 through the AMPK signaling pathway.

Characterization of circular RNA expression profiles in the age-related thymic involution of Magang goose

The thymus is a vital immune organ, but its function gradually declines with age. Circular RNAs (circRNAs) are related to the development of tissues and organs. In this study, bioinformatics analysis showed that 1329, 755, and 417 circRNAs were differentially expressed between the comparison groups of 6-month age (M6) and 20-embryo age (E20), 3-day post-hatch (P3), and 3-month age (M3) Magang geese, respectively. Among them, 167 circRNAs were differentially co-expressed between thymic development (E20, P3, and M3) and involution (M6). Functional analysis showed significant enrichment of phosphorylation and positive regulation of GTPase activity. Furthermore, pathway analysis has shown that glycerolipid metabolism and the Wnt signaling pathway are critical pathways in the thymic involution process. Finally, we constructed the competitive endogenous RNA (ceRNA) network. The results of this study suggest that circRNAs may be involved in the age-related thymic involution of the Magang goose.

Effects of Oxidative Stress on the Autophagy and Apoptosis of Granulosa Cells in Broody Geese

Broodiness is an unfavorable trait associated with the cessation of egg laying. Studies have found that excessive granulosa cell apoptosis and autophagy occur during goose broodiness. Other studies have also confirmed that oxidative stress is an important cause of apoptosis and autophagy. However, whether oxidative stress occurs during goose broodiness and whether this oxidative stress causes apoptosis and autophagy have not been fully elucidated. In this study, we investigated the effects of oxidative stress on the autophagy and apoptosis of granulosa cells in broody geese. The results showed higher mRNA expression of genes related to antioxidative stress responses (GPX, SOD-1, SOD-2, COX-2, CAT and hsp70) in pre-broody and broody geese than in laying birds. In addition, increased levels of granulosa cell apoptosis and autophagy were observed in pre-broody geese than in laying geese. Additionally, granulosa cells treated with H₂O₂ exhibited increased apoptosis and autophagy in vitro, and these effects were responsible for goose granulosa cell death. Moreover, vitamin E treatment effectively protected granulosa cells from H₂O₂-induced oxidative stress by inhibiting ROS production. Correspondingly, granulosa cell apoptosis and autophagy were greatly alleviated by vitamin E treatment. Together, our results demonstrated serious oxidative stress and granulosa cell apoptosis and autophagy in broody geese, and oxidative stress promoted apoptosis and autophagy. Vitamin E alleviated the autophagy and apoptosis of granulosa cells by inhibiting oxidative stress.

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.

Nonlinear model fitting analysis of feather growth and development curves in the embryonic stages of Jilin white geese (Anser cygnoides)

Poultry is subject to varying degrees of feather loss and feather pecking during production, which seriously affects the live appearance and carcass appearance of their commercial traits and greatly reduces the production profitability of the farming enterprise. It also has an impact on down production and quality in the case of geese. In this study, mathematical models (Logistic, Gompertz, and Von Bertalanffy) were used to assess feather growth and development during the embryonic period in Jilin white geese (Anser cygnoides) predicting the weight and length of feathers from the back, chest, and belly tracts at different embryonic ages, to determine which growth model more accurately described feather growth patterns. The result first showed that the primary feather follicles of the Jilin white goose developed at E14 and secondary feather follicles at E18; primary feather follicle density increased and then decreased, whereas secondary feather follicle density increased continuously and the primary and secondary feather follicles developed independently. Secondly, the embryonic feather growth followed a slow-fast-slow pattern, with feathers growing slowly from E12 to E18, quickly from E18 to E24, and then decreasing after E24 until just before emergence (E30). In addition, before E14, feathers were concentrated in the back tracts, and no feathers were found on the head, neck, chest, abdomen, or wings. By E22, the whole body of the embryo was covered with feathers, and the back feathers were the earliest and fastest to develop. Compared to the Gompertz and von Bertalanffy models, the logistic model fit (R2 = 0.997) was the highest, while the sum of residual squares (RSS = 25661.67), Akaike's information criterion (AIC = 77.600), Bayesian information criterion (BIC = 78.191), and mean square error (MSE = 2851.296) were the lowest. Therefore, the logistic model was more suitable for describing the changes in whole-body feather growth during the embryonic period in Jilin white geese. In conclusion, using the growth curve model to explain the relationship between feather growth and embryonic age in geese will potentially speed up the process of genetic improvement in Jilin white geese (A. cygnoides) and thus provide scientific support for molecular genetic breeding.

Glucose-induced enhanced anti-oxidant activity inhibits apoptosis in goose fatty liver

The development of mammalian nonalcoholic fatty liver disease is associated with oxidative stress, reduced mitochondrial function, and increased apoptosis in hepatocytes; however, the expressions of mitochondria-related genes are elevated in goose fatty liver, suggesting that there may be a unique protective mechanism in goose fatty liver. The aim of the study was to investigate this protective mechanism in terms of anti-oxidant capacity. Our data showed no substantial differences in the mRNA expression levels of the apoptosis-related genes including B-cell lymphoma-2 (Bcl-2), BCL2-associated X (Bax), cysteinyl aspartate-specific proteinase-3 (Caspase-3), and cysteinyl aspartate-specific proteinase-9 (Caspase-9) in the livers of the control and overfeeding Lander geese groups. The protein expression levels of Caspase-3 and cleaved Caspase-9 were not markedly different between the groups. Compared with the control group, malondialdehyde content was significantly lower (P < 0.01), glutathione peroxidase (GSH-Px) activity, glutathione (GSH) content, and mitochondrial membrane potential levels were higher (P < 0.01) in the overfeeding group. The mRNA expression levels of the anti-oxidant genes superoxide dismutase 1 (SOD1), glutathione peroxidase 1 (GPX1), and glutathione peroxidase 2 (GPX2) were increased in goose primary hepatocytes after 40 mM and 60 mM glucose treatment. Reactive oxygen species (ROS) levels were significantly reduced (P < 0.01), whereas the mitochondrial membrane potential was maintained at normal levels. The mRNA expression levels of the apoptosis-related genes Bcl-2, Bax, and Caspase-3 were not substantial. There were no significant differences in the expression levels of Caspase-3 and cleaved Caspase-9 proteins. In conclusion, glucose-induced enhanced anti-oxidant capacity may help protect the function of mitochondria and inhibit the occurrence of apoptosis in goose fatty liver.

Chromosome-level genome assembly of goose provides insight into the adaptation and growth of local goose breeds

BACKGROUND

Anatidae contains numerous waterfowl species with great economic value, but the genetic diversity basis remains insufficiently investigated. Here, we report a chromosome-level genome assembly of Lion-head goose (Anser cygnoides), a native breed in South China, through the combination of PacBio, Bionano, and Hi-C technologies.

FINDINGS

The assembly had a total genome size of 1.19 Gb, consisting of 1,859 contigs with an N50 length of 20.59 Mb, generating 40 pseudochromosomes, representing 97.27% of the assembled genome, and identifying 21,208 protein-coding genes. Comparative genomic analysis revealed that geese and ducks diverged approximately 28.42 million years ago, and geese have undergone massive gene family expansion and contraction. To identify genetic markers associated with body weight in different geese breeds, including Wuzong goose, Huangzong goose, Magang goose, and Lion-head goose, a genome-wide association study was performed, yielding an average of 1,520.6 Mb of raw data that detected 44,858 single-mucleotide polymorphisms (SNPs). Genome-wide association study showed that 6 SNPs were significantly associated with body weight and 25 were potentially associated. The significantly associated SNPs were annotated as LDLRAD4, GPR180, and OR, enriching in growth factor receptor regulation pathways.

CONCLUSIONS

We present the first chromosome-level assembly of the Lion-head goose genome, which will expand the genomic resources of the Anatidae family, providing a basis for adaptation and evolution. Candidate genes significantly associated with different goose breeds may serve to understand the underlying mechanisms of weight differences.

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.

Molecular characterization and antibacterial activities of a goose-type lysozyme gene from roughskin sculpin (Trachidermus fasciatus)

Lysozymes, acting as antimicrobial molecules, play a vital role in the host's innate immune response to pathogen infections. In the present study, a g-type lysozyme gene termed Tf-LyzG from roughskin sculpin, Trachidermus fasciatus was firstly reported. The deduced amino acid sequence of Tf-LyzG contained 188 residues and possessed conserved catalytic residues (Glu⁷¹, Asp⁸⁴, and Asp⁹⁵). Gene expression analysis revealed that Tf-LyzG was widely distributed in the tested eleven tissues with the highest expression in the gill and could be significantly induced post lipopolysaccharide (LPS) challenge. The lysozyme activity of the purified recombinant protein (rTf-LyzG) was found to be most active at pH 5.5 and 37 °C. rTf-LyzG exhibited a wide spectrum of potent bacteriolytic activity against four Gram-positive bacteria and six Gram-negative bacteria. It also displayed a high affinity to polysaccharides on bacteria surfaces including LPS, lipoteichoic acid (LTA), and peptidoglycan (PGN). rTf-LyzG was capable of binding and agglutinating all nine bacteria. Flow cytometry assay further revealed that rTf-LyzG could disrupt the membrane of Micrococcus lysodeikticus which is confirmed by scanning electron microscope (SEM) analysis that reveals blebs around the bacterial cell membrane. In summary, these data indicate that Tf-LyzG is of great importance in the fish immune response against pathogens invasion.

Single Nucleotide Polymorphism Discovery and Genetic Differentiation Analysis of Geese Bred in Poland, Using Genotyping-by-Sequencing (GBS)

Poland is the largest European producer of goose, while goose breeding has become an essential and still increasing branch of the poultry industry. The most frequently bred goose is the White Kołuda^® breed, constituting 95% of the country’s population, whereas geese of regional varieties are bred in smaller, conservation flocks. However, a goose’s genetic diversity is inaccurately explored, mainly because the advantages of the most commonly used tools are strongly limited in non-model organisms. One of the most accurate used markers for population genetics is single nucleotide polymorphisms (SNP). A highly efficient strategy for genome-wide SNP detection is genotyping-by-sequencing (GBS), which has been already widely applied in many organisms. This study attempts to use GBS in 12 conservative goose breeds and the White Kołuda^® breed maintained in Poland. The GBS method allowed for the detection of 3833 common raw SNPs. Nevertheless, after filtering for read depth and alleles characters, we obtained the final markers panel used for a differentiation analysis that comprised 791 SNPs. These variants were located within 11 different genes, and one of the most diversified variants was associated with the EDAR gene, which is especially interesting as it participates in the plumage development, which plays a crucial role in goose breeding.

Ecology and population structure of some indigenous geese breeds and the impact of four GH and Pit-1 SNPs on their body weights

This study aims to determine the genetic correlation using nine microsatellite markers to reconstruct the history of some indigenous geese populations, along with the use of four single nucleotides polymorphisms (SNPs) to investigate their correlation with the geese body weight. Microsatellite markers are mainly used to provide updated information on changes in the population structure of geese breeds. The eight goose populations reported 24% private alleles specific for each population. Expected heterozygosity (He) ranged from 0.46 to 0.70. Three breeds were reported highly polymorphic. Inbreeding coefficient (Fis) revealed that three breeds were in a minimum level of extinction danger, while one breed was in a potential endangered situation. Phylogenetic tree, principal component analysis (PCA), and self-organizing map (SOM) were constructed using MATLAB to study the population distribution and relationship among these breeds. Four SNPs were detected, two SNPs at GH gene exon (C123T and C158T), and two SNPs at Pit-1 gene exons (G161A and T282G). Four SNP loci were reported to have a significant effect on geese body weight. They were CT genotype for C123T locus, TT genotype for C158T locus, GG genotype for G161A locus, and GG genotype for T282G locus.

Transcriptome-wide study revealed m6A regulation of embryonic muscle development in Dingan goose (Anser cygnoides orientalis)

BACKGROUND

The number of myofiber is determined during the embryonic stage and does not increase during the postnatal period for birds, including goose. Thus, muscle production of adult goose is pre-determined during embryogenesis. Previous studies show N6-methyladenosine (m6A) is an important regulator for skeletal muscle development of birds and miRNAs play as a co-regulator for the skeletal muscle development in birds. Herein, we sequenced m6A and miRNA transcriptomes to investigate the profiles of m6A and their potential mechanism of regulating breast muscle development in Dingan Goose.

RESULTS

We selected embryonic 21th day (E21) and embryonic 30th day (E30) to investigate the roles of transcriptome-wide m6A modification combining with mRNAs and miRNAs in goose breast muscle development. In this study, m6A peaks were mainly enriched in coding sequence (CDS) and start codon and397 genes were identified as differentially methylated genes (DMGs). GO and KEGG analysis showed that DMGs were highly related to cellular and metabolic process and that most DMGs were enriched in muscle-related pathways including Wnt signaling pathway, mTOR signaling and FoxO signaling pathway. Interestingly, a negative correlation between m6A methylation level and mRNA abundance was found through the analysis of m6A-RNA and RNA-seq data. Besides, we found 26 muscle-related genes in 397 DMGs. We also detected 228 differentially expressed miRNAs (DEMs), and further found 329 genes shared by the target genes of DEMs and DMGs (m6A-miRNA-genes), suggesting a tightly relationship between DEMs and DMGs. Among the m6A-miRNA-genes, we found 10 genes are related to breast muscle development. We further picked out an m6A-miRNA-gene, PDK3, from the 10 genes to visualize it and the result showed differentially methylated peaks on the mRNA transcript consistent with our m6A-seq results.

CONCLUSION

GO and KEGG of DMGs between E21 and E30 showed most DMGs were muscle-related. In total, 228 DEMs were found, and the majority of DMGs were overlapped with the targets of DEGs. The differentially methylated peaks along with an m6A-miRNA-gene, PDK3, showed the similar results with m6A-seq results. Taken together, the results presented here provide a reference for further investigation of embryonic skeletal muscle development mechanism in goose.

Histomorphology and gene expression profiles during early ovarian folliculogenesis in duck and goose

In contrast to the later stages of follicle development, little is known about the characteristics and mechanisms associated with early folliculogenesis in avian species. The objectives of the present study were to examine and compare the histomorphological and molecular changes of primordial, primary, and secondary follicles from duck and goose ovaries during the first 6 post-hatching week. Morphological analysis showed that the length and width of both duck and goose ovaries increased steadily during weeks 1 to 5 but increased acutely at week 6, whereas a greater increment was observed in the ovarian length of ducks than that of geese during weeks 4 to 5. Furthermore, smaller diameters of the 3 categories of follicles were observed in ducks than those in geese at the first appearance, but they reached a similar size at week 6. More importantly, secondary follicles were found in the ovaries of ducks 1 wk earlier than in those of geese. These results indicated a more rapid growth rate for ovarian follicles in ducks than in geese during early post-hatching development. At the molecular level, it was found that the mRNAs encoding follicle stimulating hormone receptor (FSHR), anti-Müllerian hormone (AMH), B-cell leukemia/lymphoma 2, and cysteine-dependent aspartate specific protease 3 (CASPASE3) were ubiquitously expressed in all ovarian follicles of ducks and geese with different expression profiles in each follicular category during the first 6 post-hatching week. Notably, transcript levels of FSHR, AMH, and CASPASE3 changed differently between ducks and geese during weeks 5 to 6, which was postulated to be one of the mechanisms inducing more rapid growth of ovarian follicles in ducks rather than in geese. In conclusion, our results revealed, for the first time, differences in early folliculogenesis, including the rate of growth of each follicular category and the timing of transition of primary to secondary follicles, between ducks and geese, and these differences could result from different expression profiles of FSHR, AMH, and CASPASE3 during early post-hatching development.

Molecular characterization and expression analysis of the chicken-type and goose-type lysozymes from totoaba (Totoaba macdonaldi)

Lysozymes play a key role in innate immune response to bacterial pathogens, catalyzing the hydrolysis of the peptidoglycan layer of bacterial cell walls. In this study, the genes encoding the c-type (TmLyzc) and g-type (TmLyzg) lysozymes from Totoaba macdonaldi were cloned and characterized. The cDNA sequences of TmLyzg and TmLyzc were 582 and 432 bp, encoding polypeptides of 193 and 143 amino acids, respectively. Amino acid sequences of these lysozymes shared high identity (60-90%) with their counterparts of other teleosts and showed conserved functional-structural signatures of the lysozyme superfamily. Phylogenetic analysis indicated a close relationship with their vertebrate homologues but distinct evolutionary paths for each lysozyme. Expression analysis by qRT-PCR revealed that TmLyzc was expressed in stomach and pyloric caeca, while TmLyzg was highly expressed in stomach and heart. These results suggest that both lysozymes play important roles in defense of totoaba against bacterial infections or as digestive enzyme.

Pacific Biosciences assembly with Hi-C mapping generates an improved, chromosome-level goose genome

BACKGROUND

The domestic goose is an economically important and scientifically valuable waterfowl; however, a lack of high-quality genomic data has hindered research concerning its genome, genetics, and breeding. As domestic geese breeds derive from both the swan goose (Anser cygnoides) and the graylag goose (Anser anser), we selected a female Tianfu goose for genome sequencing. We generated a chromosome-level goose genome assembly by adopting a hybrid de novo assembly approach that combined Pacific Biosciences single-molecule real-time sequencing, high-throughput chromatin conformation capture mapping, and Illumina short-read sequencing.

FINDINGS

We generated a 1.11-Gb goose genome with contig and scaffold N50 values of 1.85 and 33.12 Mb, respectively. The assembly contains 39 pseudo-chromosomes (2n = 78) accounting for ∼88.36% of the goose genome. Compared with previous goose assemblies, our assembly has more continuity, completeness, and accuracy; the annotation of core eukaryotic genes and universal single-copy orthologs has also been improved. We have identified 17,568 protein-coding genes and a repeat content of 8.67% (96.57 Mb) in this genome assembly. We also explored the spatial organization of chromatin and gene expression in the goose liver tissues, in terms of inter-pseudo-chromosomal interaction patterns, compartments, topologically associating domains, and promoter-enhancer interactions.

CONCLUSIONS

We present the first chromosome-level assembly of the goose genome. This will be a valuable resource for future genetic and genomic studies on geese.

Exploring the genetic basis of fatty liver development in geese

Although geese possess an adaptive physiological capacity for lipid storage, few candidate genes contributing to this ability are characterised. By comparing the genomes of individuals with extremely high and low fatty liver weights (FLW), candidate genes were identified, including ARAP2, GABRE, and IL6. Single-nucleotide polymorphisms in or near these genes were significantly (p < 0.05) associated with carcass traits (FLW) and biochemical indexes (very-low-density lipoprotein and N-terminal procollagen III), suggesting contribution to trait variation. A common variant at the 5'-end of LCORL explained ~ 18% and ~ 26% of the phenotypic variance in body weight with/without overfeeding and had significant effects on FLW (p < 0.01). ZFF36L1, ARHGEF1 and IQCJ, involved in bile acid metabolism, blood pressure, and lipid concentration modulation, were also identified. The presence of highly divergent haplotypes within these genes suggested involvement in protection against negative effects from excessive lipids in the liver or circulatory system. Based on this and transcriptomic data, we concluded that geese hepatosteatosis results from severe imbalance between lipid accumulation and secretion, comparable to human non-alcohol fatty liver disease but involving other genes. Our results provided valuable insights into the genesis of geese fatty liver and detected potential target genes for treatment of lipid-related diseases.

A 14-bp insertion in endothelin receptor B-like (EDNRB2) is associated with white plumage in Chinese geese

BACKGROUND

Gang goose is a native species with gray plumage in Sichuan, China. As a result of overhunting, the number of gray Gang geese has decreased dramatically. To keep the species from extinction, conservation work for Gang geese was undertaken. In the process of pure breeding of gray Gang geese, approximately 2% of the offspring of each generation were white. This study aims to explain the genetic mechanism of this phenomenon and provide reliable molecular markers for goose-related plumage color breeding.

RESULTS

We used the method of pooled whole genome sequencing and Fst (fixation statistics) to identify the differentiation degree of alleles between gray Gang geese and white Gang geese from their offspring. In this way, EDNRB2, a key gene that affects the migration of melanoblasts, was identified. Then, the transcriptome was sequenced for the two geese plumage color populations, and the DEGs (differentially expressed genes) were analyzed. The results indicated that EDNRB2, as a possible candidate gene, had a significantly differential mRNA expression. In addition, a 14-bp insertion (NW_013185915.1: g. 750,748-750,735 insertion. CACAGGTGAGCTCT) in exon 3 of EDNRB2 was analyzed and found to have a significant association between gray geese and Chinese white breeds (P = 0.00), while this mutation was not found in European geese. Meanwhile, the insertion was homozygous in all the white geese we detected and heterozygous in gray geese, indicating that this mutation is recessive. Furthermore, this 14-bp insertion leads to a frameshift mutation in the EDNRB2 coding region and nonsense-mediated mRNA decay.

CONCLUSION

Our study strongly suggests that the 14-bp insertion in exon 3 of the EDNRB2 gene is associated with the white plumage phenotype in Chinese geese. This study is the first to investigate the relationship between EDNRB2 and white plumage in geese.

Molecular cloning and functional analysis of TRAF6 from Yangzhou great white goose Anser anser

TNF receptor-associated factor 6 (TRAF6) is an adaptor protein and an E3 ubiquitin ligase mediating multiple cell signaling pathway activation in a context-dependent manner. TRAF6 plays critical roles in innate immune response and regulates function of antigen-presenting cells. Here, we cloned the goose TRAF6 (goTRAF6) gene from a healthy Yangzhou great white goose (Anser anser), which had a typical TRAF structure and shared a high-sequence identity with TRAF6 of other birds. Quantitative real-time PCR revealed that goTRAF6 mRNA was broadly expressed in all the studied tissues, with highest expression in the heart and pectoral muscle. Overexpression of goTRAF6 caused NF-κB activation in a dose-dependent manner and substantially upregulated IFN-β expression in HEK293T cells. Following Toll-like receptor (TLR) ligand stimulation of goose peripheral blood mononuclear cells, goTRAF6 and downstream inflammatory cytokine mRNA levels considerably up-regulated, especially at early stages. Salmonella Enteritidis challenge caused overexpression of goTRAF6 and cytokine mRNA in all the examined organs. These findings demonstrated that goTRAF6 played a substantial role in TLR-TRAF6 signaling cascade, and further contributed to the antibacterial-responses in host.

De Novo Assembly and Comparative Transcriptome Profiling of Anser anser and Anser cygnoides Geese Species' Embryonic Skin Feather Follicles

Geese feather production and the quality of downy feathers are additional economically important traits in the geese industry. However, little information is available about the molecular mechanisms fundamental to feather formation and the quality of feathers in geese. This study conducted de novo transcriptome sequencing analysis of two related geese species using the Illumina 4000 platform to determine the genes involved in embryonic skin feather follicle development. A total of 165,564,278 for Anser anser and 144,595,262 for Anser cygnoides clean reads were generated, which were further assembled into 77,134 unigenes with an average length of 906 base pairs in Anser anser and 66,041 unigenes with an average length of 922 base pairs in Anser cygnoides. To recognize the potential regulatory roles of differentially expressed genes (DEGs) during geese embryonic skin feather follicle development, the obtained unigenes were annotated to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis. In both species, GO and KOG had shown similar distribution patterns during functional annotation except for KEGG, which showed significant variation in signaling enrichment. Anser asnser was significantly enriched in the calcium signaling pathway, whereas Anser cygnoides was significantly enriched with glycerolipid metabolism. Further analysis indicated that 14,227 gene families were conserved between the species, among which a total of 20,715 specific gene families were identified. Comparative RNA-Seq data analysis may reveal inclusive knowledge to assist in the identification of genetic regulators at a molecular level to improve feather quality production in geese and other poultry species.

Evolutional conservation of molecular structure and antiviral function of a type I interferon, IFN-kappa, in poultry

IFN-kappa (IFN-κ) is a type I IFN expressed by keratinocytes, monocytes and dendritic cells with important roles during the innate immune response period. This research was conducted to elaborate the evolution and characteristics of IFN-κ in poultry. Chicken IFN-κ is located on the sex-determining Z chromosome, which is greatly different from mammals. Poultry IFN-κ cluster together in a species-specific manner through positive selection pressure and share only 19-33% homology with mammalian IFN-κ and poultry other type I IFN. Both chicken and duck IFN-κ was constitutively expressed in spleen, skin, lung, and peripheral blood mononuclear cells (PBMC), as well as being significantly induced after treatment with virus in PBMC. Biologically, poultry IFN-κ has antiviral activity against VSV in chicken embryonic fibroblasts and duck embryonic fibroblasts (CEF and DEF) cells, and induces the expression of IFN stimulated genes (ISGs). After treatment with JAK1 inhibitor, the ISGs expression can be down-regulated. Overall, our research on poultry IFN-κ not only enriches the knowledge about IFN-κ but also facilitates further research on the role of type I IFNs in antiviral defense responses in poultry.

A Functional Mutation in KIAA1462 Promoter Decreases Glucocorticoid Receptor Affinity and Affects Egg-Laying Performance in Yangzhou Geese

The identification of genetic markers is valuable for improving the egg-laying performance in goose production. The single-nucleotide polymorphism (SNP) rs1714766362 in an intron of the goose KIAA1462 gene was found to be relevant to laying performance in our previous study. However, its function remains unclear. In this study, the full-length coding sequence of KIAA1462 gene was firstly characterized in Yangzhou geese. Q-PCR (Quantitative Real Time Polymerase Chain Reaction) results showed that KIAA1462 was highly expressed in the liver, ovary, and mature F1 follicles. For SNP rs1714766362, geese with the AA genotype showed better laying performance than the TT ones and exhibited a higher KIAA1462 expression level in the ovary. Gain- and loss-of function experiments in granulosa cells revealed that KIAA1462 affected the expression of the apoptosis marker gene caspase-3. Considering that rs1714766362 locates in an intron area, we compared the KIAA1462 promoter regions of AA and TT individuals and identified the SNP c.-413C>G (Genbank ss2137504176), which was completely linked to SNP rs1714766362. According to the transcription factor prediction results, the glucocorticoid receptor (GR) would bind to the SNP site containing the C but not the G allele. In this study, we proved this hypothesis by an electrophoretic mobility shift assay (EMSA). In summary, we identified a novel mutation in the promoter of KIAA1462 gene which can modulate GR binding affinity and affect the laying performance of geese.

Contribution of both positive selection and relaxation of selective constraints to degeneration of flyability during geese domestication

Flyability is the most discrepant trait between modern-day geese and their wild ancestors, and the degeneration of flyability is a key marker of the successful domestication of wild geese. In light of the relatively short history of domestic geese, intense artificial selection is thought to play an important role in the degeneration of flyability. However, the underlying mechanism behind this phenomenon has seldom been investigated. In this study, we applied a molecular evolutionary approach to the evaluation of partial breeds of domestic geese in order to look for genes involved in the selection pressure toward degeneration of flyability. The haplotype networks, pairwise fixation index (FST) values, and analysis of molecular variance results all clearly illustrated a population variance between Landes geese and partial Chinese domestic geese. We also detected signatures of positive artificial selection in the COX2 and COX3 genes, and related selection in the HBB gene. Our results support the independent origins of partial European domestic geese and Chinese domestic geese. In addition, both positive artificial selection and the relaxation of functional constraints appeared to play important roles in the degeneration of flyability in domestic geese.

Comparative proteomic analysis of pituitary glands from Huoyan geese between pre-laying and laying periods using an iTRAQ-based approach

In this study, we performed a comprehensive evaluation of the proteomic profile of the pituitary gland of the Huoyan goose during the laying period compared to the pre-laying period using an iTRAQ-based approach. Protein samples were prepared from pituitary gland tissues of nine pre-laying period and nine laying period geese. Then the protein samples from three randomly selected geese within each period were pooled in equal amounts to generate one biological sample pool. We identified 684 differentially expressed proteins, including 418 up-regulated and 266 down-regulated proteins. GO annotation and KEGG pathway analyses of these proteins were conducted. Some of these proteins were found to be associated with hormone and neurotransmitter secretion and transport, neuropeptide signalling and GnRH signalling pathways, among others. Subsequently, the modification of the abundance of three proteins (prolactin, chromogranin-A and ITPR3) was verified using western blotting. Our results will provide a new source for mining genes and gene products related to the egg-laying performance of Huoyan geese, and may provide important information for the conservation and utilization of local goose breeds.