Comparative transcriptome analysis identifies crucial candidate genes and pathways in the hypothalamic-pituitary-gonadal axis during external genitalia development of male geese.

BACKGROUND: All birds reproduce via internal fertilization, but only ~3% of male birds possess the external genitalia that allows for intromission. Waterfowl (e.g., duck and goose) are representatives of them, and the external genitalia development of male geese is directly related to mating ability. Notably, some male geese show abnormal external genitalia development during ontogenesis. However, until now little is known about the molecular mechanisms of the external genitalia development in goose. In the present study, comparative transcriptomic analyses were performed on the hypothalamus, pituitary gland, testis, and external genitalia isolated from the 245-day-old male Tianfu meat geese showing normal (NEGG, n = 3) and abnormal (AEGG, n = 3) external genitals in order to provide a better understanding of the mechanisms controlling the development of the external genitalia in aquatic bird species. RESULTS: There were 107, 284, 2192, and 1005 differentially expressed genes (DEGs) identified in the hypothalamus, pituitary gland, testis and external genitalia between NEGG and AEGG. Functional enrichment analysis indicated that the DEGs identified in the hypothalamus were mainly enriched in the ECM-receptor interaction pathway. The ECM-receptor interaction, focal adhesion, and neuroactive ligand-receptor interaction pathways were significantly enriched by the DEGs in the pituitary gland. In the testis, the DEGs were enriched in the neuroactive ligand-receptor interaction, cell cycle, oocyte meiosis, and purine metabolism. In the external genitalia, the DEGs were enriched in the metabolic, neuroactive ligand-receptor interaction, and WNT signaling pathways. Furthermore, through integrated analysis of protein-protein interaction (PPI) network and co-expression network, fifteen genes involved in the neuroactive ligand-receptor interaction and WNT signaling pathways were identified, including KNG1, LPAR2, LPAR3, NPY, PLCB1, AVPR1B, GHSR, GRM3, HTR5A, FSHB, FSHR, WNT11, WNT5A, WIF1, and WNT7B, which could play crucial roles in the development of goose external genitalia. CONCLUSIONS: This study is the first systematically comparing the hypothalamus, pituitary gland, testis, and external genitalia transcriptomes of male geese exhibiting normal and abnormal external genitals. Both bioinformatic analysis and validation experiments indicated that the neuroactive ligand-receptor interaction pathway could regulate the WNT signaling pathway through PLCB1 to control male goose external genitalia development.

Oestrogen promotes lipids transportation through oestrogen receptor α in hepatic steatosis of geese in vitro.

Evidence has shown that oestrogen suppresses lipids deposition in the liver of mammals. However, the molecular mechanism of oestrogen action in hepatic steatosis of geese liver has yet to be determined. This study aimed to investigate the effect of oestrogen on lipid homeostasis at different states of geese hepatocytes in vitro. The results showed that an in vitro model of hepatic steatosis was induced by 1.5 mM sodium oleate via detecting the viability of hepatocytes and content of lipids. When the normal hepatocytes were administrated with different concentrations of oestrogen (E2 ), the expression levels of diacylglycerol acyltransferase 2 (DGAT2), microsomal triglyceride transfer protein (MTTP) and oestrogen receptors (ERs, alpha and beta) were up-regulated only at high concentrations of E2 , whereas the lipid content was not a significant difference. In goose hepatocytes of hepatic steatosis, however, the expression levels of MTTP, apolipoprotein B (apoB) and ERα/ß significantly increased at 10-7 or 10-6  M E2 . Meanwhile, the lipids content significantly increased at 10-9 and 10-8  M E2 and decreased at 80 µM E2 . Further heatmap analysis showed that ERα was clustered with apoB and MTTP in either normal hepatocytes or that of hepatic steatosis. Taken together, E2  might bind to ERα to up-regulate the expression levels of apoB and MTTP, promoting the transportation of lipids and alleviating lipids overload in hepatic steatosis of geese in vitro.

Whole transcriptome analysis revealed the regulatory network and related pathways of non-coding RNA regulating ovarian atrophy in broody hens.

Poultry broodiness can cause ovarian atresia, which has a detrimental impact on egg production. Non-coding RNAs (ncRNAs) have become one of the most talked-about topics in life sciences because of the increasing evidence of their novel biological roles in regulatory systems. However, the molecular mechanisms of ncRNAs functions and processes in chicken ovarian development remain largely unknown. Whole-transcriptome RNA sequencing of the ovaries of broodiness and laying chickens was thus performed to identify the ncRNA regulatory mechanisms associated with ovarian atresia in chickens. Subsequent analysis revealed that the ovaries of laying chickens and those with broodiness had 40 differentially expressed MicroRNA (miRNAs) (15 up-regulated and 25 down-regulated), 379 differentially expressed Long Noncoding RNA (lncRNAs) (213 up-regulated and 166 down-regulated), and 129 differentially expressed circular RNA (circRNAs) (63 up-regulated and 66 down-regulated). The competing endogenous RNAs (ceRNA) network analysis further revealed the involvement of ECM-receptor interaction, AGE-RAGE signaling pathway, focal adhesion, cytokine-cytokine receptor interaction, inflammatory mediator regulation of TRP channels, renin secretion, gap junction, insulin secretion, serotonergic synapse, and IL-17 signaling pathways in broodiness. Upon further analysis, it became evident that THBS1 and MYLK are significant candidate genes implicated in the regulation of broodiness. The expression of these genes is linked to miR-155-x, miR-211-z, miR-1682-z, gga-miR-155, and gga-miR-1682, as well as to the competitive binding of novel_circ_014674 and MSTRG.3306.4. The findings of this study reveal the existence of a regulatory link between non-coding RNAs and their competing mRNAs, which provide a better comprehension of the ncRNA function and processes in chicken ovarian development.

Integrative metabolomics and transcriptomics analysis revealed specific genes and metabolites affecting meat quality of chickens under different rearing systems.

Different rearing systems have varying effect on animal welfare and meat quality of poultry. Currently, there are no established standards for the rearing systems of Chinese indigenous chickens. Our study aimed to investigate the effects of different rearing systems on the meat quality, gene profiles, and metabolites of Chinese indigenous chickens (Nanchuan chicken). 10-wk-old Nanchuan chickens (n=360) were randomly divided into 3 groups (cage, net, and free-range groups), with 6 replicates per group (20 chickens per replicate). The experiment lasted for 12 wk. At 154-days-old, 36 healthy chickens (6 males and 6 females per group) were randomly selected, euthanized, and their breast muscles were collected to assess the meat quality parameters and histomorphological characteristics. Additionally, breast muscles from 18 random hens (3 males and 3 females per group) were used for metabolomics and RNA-seq analysis. The results showed that rearing systems significantly affected the meat quality and myofiber characteristics. The meat quality of breast muscles from free-range chickens was superior to that of caged chickens, characterized by more tender meat and smaller myofiber cross-sectional areas. Integrative metabolomics and transcriptomics analysis revealed that the differentially expressed genes of chicken breast muscles were primarily involved in the myofiber differentiation. Mechanically, the improved meat quality of breast muscle in free-range chickens were mainly associated with enhanced skeletal muscle differentiation facilitated by fibromodulin, increased levels of up-regulated Acetyl-L-carnitine and Propionylcarnitine level, and decreased levels of Nonanoic acid and Elaidic acid abundance (Graphical abstract). This provides a comprehensive understanding of the most effective and sustainable breeding, production, and rearing systems for Chinese indigenous chickens. It also contributes to the current knowledge of the molecular mechanisms underlying the effects of rearing systems on growth performance and meat quality of chickens.

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.

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.

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.

Research Note: Integrated transcriptomic and metabolomic analysis reveals potential candidate genes and regulatory pathways associated with egg weight in ducks.

Egg weight is an important indicator of egg phenotypic traits, which directly affects the economic benefits of the poultry industry. In the present research, laying ducks were classified into high egg weight (HEW) and light egg weight (LEW) groups. To reveal the underlying mechanism that may be responsible for the egg weight difference, the integrated analysis of transcriptomes and serum metabolomics was performed between the two groups. The results showed extremely significant differences (P < 0.01) in the total egg weight at 300 d, and average egg weight between the HEW and LEW groups. 733, 591, 82, and 74 differentially expressed genes (DEGs) were identified in the liver, magnum, F1, and F5 (hierarchical follicles) follicle membrane, respectively. The candidate genes were screened further from the perspective of forming an egg. In terms of egg yolk formation, the functional analysis revealed fatty acid metabolism-related pathways account for 36% of the liver's top pathways, including fatty acid biosynthesis, folate biosynthesis, fatty acid metabolism, and glycerol lipid metabolism pathways. FASN gene was identified as the key candidate gene by comprehensive analysis of gene expression and protein-protein interaction (PPI) network. In the follicle membrane, the DEGs were mainly enriched in protein processing in the endoplasmic reticulum, and MAPK signaling pathway, and HSPA2, HSPA8, BAG3 genes were identified as crucial candidate genes. In terms of egg white formation, the functional analysis revealed protein metabolism-related pathways account for 40% of the magnum's top pathways, which includes protein processing in the endoplasmic reticulum pathway. HSP90AA1 and HSPA8 genes were identified as key candidate genes. In addition, the integrated transcriptomic and metabolomic analysis showed that arginine and proline metabolism pathways could contribute to differences in egg weight. Thus, we speculated that the potential candidate genes, regulatory pathways, and metabolic biomarkers mentioned above might be responsible for the egg weight difference. These findings might provide a theoretical basis for improving the egg weight of ducks.

Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese.

The lower egg production of geese (20~60 eggs per year) compared with chicken and duck limits the development of the industry, while the yolk weight and fatty liver susceptibility of geese was higher than that of other poultry. Therefore, the relationship between lipid metabolism and the laying performance of geese remains to be explored. Phenotypically, we observed that the liver fat content of the high-, low-, and no-egg production groups decreased in turn, while the abdominal fat weight increased in turn. For transcriptional regulation, the KEGG pathways related to lipid metabolism were enriched in all pairwise comparisons of abdominal fat and liver through functional analysis. However, some KEGG pathways related to inflammation and the circadian rhythm pathway were enriched by DEGs only in abdominal fat and the liver, respectively. The DEGs in ovarian stroma among different groups enriched some KEGG pathways related to ovarian steroidogenesis and cell adhesion. Our research reveals that lipid metabolism regulated by the circadian rhythm of the liver may directly or indirectly affect ovarian function through the inflammation and hormone secretion of abdominal fat. These results offer new insights into the regulation mechanisms of goose reproductive traits.

Effects of cage versus floor rearing system on goose intestinal histomorphology and cecal microbial composition.

Due to the demand for modern goose production and the pressure of environmental protection, the rearing systems of geese are changing from traditional waterside rearing to intensive rearing systems such as floor rearing (FR) and cage rearing (CR) systems. However, little is known about the effects of different rearing systems on goose intestinal functions and cecal microbial composition. Therefore, this study aimed to compare intestinal histomorphology and cecal microbial composition differences in geese reared under CR and FR at 270 d of age. Histomorphological analysis showed that the ileal villus height (VH) to crypt depth (CD) ratio was significantly greater in CR than in FR (P < 0.001). Taxonomic analysis showed that the dominant bacteria of cecal microorganisms in both rearing systems were roughly similar, with Bacteroidota, Firmicutes, Fusobacteriota, and Proteobacteria being the dominant phyla while Bacteroides, Fusobacterium, and uncultured_bacterium_o_Bacteroidales being the dominant genera. Differentially abundant taxa between CR and FR were also identified using Linear Discriminant Analysis Effect Size (LEfSe) analysis (P < 0.05, LDA score > 3.5). Megamonas and Anaerobiospirillum were significantly enriched in the CR group at the genus level, while uncultured_bacterium_f_Rikenellaceae and Sutterella were significantly enriched in the FR group. Notably, we found that the relative abundance of uncultured_bacterium_f_Rikenellaceae was significantly negatively correlated with the ileal VH and VH/CD (P < 0.05). The relative abundance of Megamonas and Anaerobiospirillum were significantly negatively correlated with abdominal fat weight and relative abdominal fat weight (P < 0.01), whereas that of Sutterella was significantly positively correlated with abdominal fat weight and relative abdominal fat weight (P < 0.01). Furthermore, PICRUSt2 analysis indicated that the lipid metabolism pathways of cecal microorganisms were lower enriched in CR than in FR. In conclusion, compared with FR, the CR significantly changed goose ileal histomorphological characteristics and cecal microbial composition, thereby affecting goose physiological functions and production performance.

Comparative Transcriptome Analysis Provided a New Insight into the Molecular Mechanisms of Epididymis Regulating Semen Volume in Drakes.

Semen volume is an important factor in artificial insemination (AI) of ducks. In drakes, seminal plasma that is produced by the epididymis determines the semen volume. However, the mechanism of epididymis regulating semen volume of drakes remains unclear. Therefore, the aim of the present study was to preliminarily reveal the mechanism regulating the semen volume through comparing the epididymal histomorphology and mRNA expression profiles between drakes with high-volume semen (HVS) and low-volume semen (LVS). Phenotypically, drakes in the HVS group produced more sperm than drakes in the LVS group. In addition, compared with the HVS group, the ductal square of ductuli conjugentes (DC) and dutus epididymidis (DE) in epididymis was significantly smaller in the LVS group, and the lumenal diameter and epithelial thickness of DC/DE were significantly shorter in the LVS group. In transcriptional regulation, 72 different expression genes (DEGs) were identified from the epididymis between HVS and LVS groups. Gene Ontology (GO) analysis indicated that the DEGs were mainly related to hormone secretion, neurotransmitter synthesis/transport, transmembrane signal transduction, transmembrane transporter activity, and nervous system development (p < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis showed that the DEGs were significantly enriched in pathways associated with hormone and neurotransmitter transmission (p < 0.05). In addition, further analysis of the top five pathways enriched by KEGG, nine key candidate genes (including SLC18A2, SNAP25, CACNA1B, GABRG2, DRD3, CAMK2A, NR5A1, and STAR) were identified, which could play a crucial role in the formation of semen. These data provide new insights into the molecular mechanism regulating semen volume of drakes and make feasible the breeding of drakes by semen volume.

A comparative proteomic study of high and low semen quality seminal plasma in drakes.

Semen quality is the most important indicator in evaluating drake fecundity. At present, the low semen quality has become a major factor restricting the development of artificial insemination (AI) technology in ducks. Numerous studies have indicated that seminal plasma proteins play a crucial role in semen quality, but the mechanism of seminal plasma proteins regulating semen quality of drakes remains unclear. Thus, the objective of this study was to identify seminal plasma proteins associated with semen quality by comparing the seminal plasma proteomic profile of drakes with high-quality semen (HQS) and low-quality semen (LQS). Using a label-free MS-based method, a total of 745 seminal plasma proteins were identified. Of these, 55 differentially expressed proteins (DEPs) were identified (40 up-regulated and 15 down-regulated). Gene Ontology (GO) analysis showed that the DEPs were mainly enriched in transmembrane transport, extracellular matrix structural constituent, transferase activity, transferring acyl groups other than amino-acyl groups, transmembrane transporter activity, and integral component of membrane (P < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis indicated that the DEPs were significantly enriched in apoptosis, tyrosine metabolism, glycerophospholipid metabolism, and sulfur metabolism pathways (P < 0.05). Moreover, through protein-protein interaction (PPI) network analysis, eight potential candidate proteins were identified, including P19140 (Alpha-enolase), R0KUV7 (Calreticulin), R0K3X3 (Solute carrier family 2, facilitated glucose transporter member 5), R0L6V0 (Proteasome subunit beta), R0JKW0 (Cytochrome c), R0JMC5 (Tubulin alpha chain), R0LCK1 (Cathepsin C), and R0JUP6 (Cathepsin D), which could play crucial roles in semen quality. Notably, further analysis demonstrated that key protein P19140 (Alpha-enolase) might can control the semen quality of drakes by regulating the expression of proteins related to apoptosis pathway. This study is the first systematically comparing the seminal plasma proteome of drakes exhibiting high and low semen quality. These results provide novel insights into the mechanisms regulating semen quality of drakes.

Identification of polymorphic loci in the deiodinase 2 gene and their associations with head dimensions in geese.

OBJECTIVE: This study was conducted to clone and compare the molecular characteristics of the deiodinase 2 (DIO2) gene between Sichuan White geese and Landes geese, and to analyze the association between polymorphisms of the DIO2 gene and head dimensions in Tianfu meat geese. METHODS: The coding sequence of the DIO2 gene was cloned by polymerase chain reaction and vector ligation and aligned by DNAMAN software. A total of 350 Tianfu meat geese were used to genotype the polymorphisms of the DIO2 gene and measure the head dimensions. Association analysis between the polymorphisms of the DIO2 gene and head dimensions was carried out. RESULTS: An 840-bp coding sequence of the DIO2 gene was obtained and comparison analysis identified four polymorphic loci between Sichuan White geese and Landes geese. Further analysis showed that the dominant alleles for the four polymorphic loci were G, G, A, and T and the frequency of the heterozygous genotype was higher than that of the homozygous genotype in Tianfu meat geese. Compared to that in the population of non-knob geese of Tianfu meat geese, the head dimensions in the population of knob geese were significantly higher except for nostril height. However, in the non-knob geese, beak width 1, beak width 2, nostril length, cranial width 1, and maxillary length had significant differences among different genotypes or haplotypes/diplotypes. CONCLUSION: These results suggested that polymorphisms of the DIO2 gene could be considered molecular markers to select larger heads of geese in the population of non-knob geese.

Role of stearyl-coenzyme A desaturase 1 in mediating the effects of palmitic acid on endoplasmic reticulum stress, inflammation, and apoptosis in goose primary hepatocytes.

OBJECTIVE: Unlike mammals, goose fatty liver shows a strong tolerance to fatty acids without obvious injury. Stearyl-coenzyme A desaturase 1 (SCD1) serves crucial role in desaturation of saturated fatty acids (SAFs), but its role in the SAFs tolerance of goose hepatocytes has not been reported. This study was conducted to explore the role of SCD1 in regulating palmitic acid (PA) tolerance of goose primary hepatocytes. METHODS: 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide was examined to reflect the effect of PA on hepatocytes viability, and quantitative polymerase chain reaction was used to detect the mRNA levels of several genes related to endoplasmic reticulum (ER) stress, inflammation, and apoptosis, and the role of SCD1 in PA tolerance of goose hepatocytes was explored using RNA interfere. RESULTS: Our results indicated that goose hepatocytes exhibited a higher tolerant capacity to PA than human hepatic cell line (LO2 cells). In goose primary hepatocytes, the mRNA levels of fatty acid desaturation-related genes (SCD1 and fatty acid desaturase 2) and fatty acid elongate enzyme-related gene (elongase of very long chain fatty acids 6) were significantly upregulated with 0.6 mM PA treatment. However, in LO2 cells, expression of ER stressrelated genes (x box-binding protein, binding immunoglobulin protein, and activating transcription factor 6), inflammatory response-related genes (interleukin-6 [IL-6], interleukin- 1ß [IL-1ß], and interferon-γ) and apoptosis-related genes (bcl-2-associated X protein, b-cell lymphoma 2, Caspase-3, and Caspase-9) was significantly enhanced with 0.6 mM PA treatment. Additionally, small interfering RNA (siRNA) mediated downregulation of SCD1 significantly reduced the PA tolerance of goose primary hepatocytes under the treatment of 0.6 mM PA; meanwhile, the mRNA levels of inflammatory-related genes (IL-6 and IL-1ß) and several key genes involved in the phosphoinositide 3-kinase (PI3K)/ protein kinase B (AKT), forkhead box O1 (FoxO1), mammalian target of rapamycin and AMPK pathways (AKT1, AKT2, FoxO1, and sirtuin 1), as well as the protein expression of cytochrome C and the apoptosis rate were upregulated. CONCLUSION: In conclusion, our data suggested that SCD1 was involved in enhancing the PA tolerance of goose primary hepatocytes by regulating inflammation- and apoptosisrelated genes expression.