Next Article in Journal
CFD Study of a Tunnel-Ventilated Compost-Bedded Pack Barn Integrating an Evaporative Pad Cooling System
Next Article in Special Issue
Presence, Tissue Localization, and Gene Expression of the Adiponectin Receptor 1 in Testis and Accessory Glands of Male Rams during the Non-Breeding Season
Previous Article in Journal
A Prospective Cohort Study on the Periparturient Muscle Tissue Mobilisation in High Producing Dairy Cows
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
*
Author to whom correspondence should be addressed.
These authors contribute equally to this work.
Animals 2022, 12(14), 1775; https://doi.org/10.3390/ani12141775
Submission received: 16 June 2022 / Revised: 5 July 2022 / Accepted: 7 July 2022 / Published: 11 July 2022
(This article belongs to the Special Issue Adipokines: Hormones Integrating Energy Homeostasis and Reproduction)

Abstract

:

Simple Summary

The importance of lipid metabolism in the egg production of poultry has been widely reported. Meanwhile, geese have lower egg production and unique lipid metabolism patterns compared with chicken and duck. It is of great significance to further improve egg laying performance to explore the differences of fat metabolism and the molecular mechanisms in geese with different egg laying performance. This study compared the phenotypic differences of liver and abdominal fat, as well as the transcriptome level differences of liver, abdominal fat, and ovarian stroma among high-, low-, and no-egg production groups. The results reveal 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.

Abstract

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.

1. Introduction

Compared with chicken and duck, the lower egg production of geese (20~60 eggs per year) limits the development of the industry. With the development of cage breeding technology of geese, we found that different individuals have different abilities to lay eggs, and there are even individuals that never lay eggs. Previous studies showed that lipids gradually accumulate from small white, small yellow to large yellow during the follicle development stage of geese [1,2]. Meanwhile, the weight of goose egg yolk is much greater than that of chicken and duck, and it may take more time to deposit enough yolk. Therefore, it is speculated that lipid metabolism (synthesis, transport, deposition) may be one of the important factors affecting the egg production of geese.
It has been proven that the ovaries of geese have the ability of de novo lipid synthesis [3] Meanwhile, most of the lipids in follicles (more than 90%) come from the transport of liver [4]. Compared with chicken, geese have a stronger capacity for lipid metabolism in the liver [5]. Geese are commonly used to produce fatty liver, and some scholars have pointed out that the difference in lipid metabolism caused by overfeeding is more obvious in the liver than in adipose tissue. The liver involves pathways such as lipid metabolism and immune response, while the adipose tissue involves pathways such as protein binding and gluconeogenesis after overfeeding. Part of the lipids synthesized by the liver are transported to the ovary to participate in the formation of egg yolk, and the others are stored in the form of sebum or abdominal fat. For obese individuals with a large amount of fat deposited in the abdominal cavity, fat will compress the ovary and fallopian tube, reducing the function of the reproductive system and affecting the formation of eggs and ovulation rate, resulting in a decline in egg production [6]. Too much or too little fat deposition has a negative effect on the initiation of egg production [7].
Studies have shown that the increase in egg production of broilers under feed restriction is due to the different amounts of lipid transport from the liver to the ovaries or sebum/abdominal fat [8]. It has been reported that lipid metabolism differences in ovarian stroma are involved in regulating egg production in geese [9]. Therefore, this study aims to explore the pathways of different lipid metabolisms affecting the reproductive performance of geese by constructing the mRNA expression profiles of liver, abdominal fat, and ovarian stroma of geese among the no-, low-, and high-egg production geese. These data will provide new insights into the regulation mechanism of goose reproductive traits.

2. Materials and Methods

2.1. Sample Collection

A total of 232 female Sichuan white geese hatched in the same batch were reared in single cages after the end of the brood period. The laying performance of 232 geese was recorded until 365 days of age. The geese that no eggs laid were divided into the no egg production group (n = 40), the geese with the top 30% egg production were divided into the high egg production group (n = 54), and the geese with the bottom 30% of egg production were divided into the low egg production group (n = 54). Meanwhile, we weighed 232 geese. Then, five individuals with similar body weight and consistent physiological states were selected from the three groups for slaughter. Abdominal fat and liver weight were recorded. Livers were dissected and placed onto a tissue mold, filled with an optimal cutting temperature compound (Fisher Scientific, Pittsburgh, PA, USA), and frozen at −20 °C. Next, abdominal fat, the partial liver, and ovarian stroma were collected and placed in liquid nitrogen.

2.2. Histological Observation

Each liver tissue was taken to a size of 24 × 24 × 2 mm, frozen and placed on a tissue-supporting device dripping with OTC embedding agent. Tissue blocks were frozen and leveled. The tissue was cut into 10 m thick slices using a CM1520 slicer (Leica, Weztlar, Germany) and attached to the anti-slip slide. Distilled water was washed for 2 min, 60% isopropanol aqueous solution (Keshi, Shanghai, China) for 2 s, oil red O working staining solution (Sigma, Shanghai, China) for 15 min, 60% isopropanol aqueous solution for 2 s, running water for 2 min, and hematoxylin (Thermo Fisher Scientific, Waltham, MA, USA) for 30 s. Then, neutral resin adhesive was used to seal the film. The staining of tissue samples was observed under a microscope (OLMPUS, Tokyo, Japan). The microscopic imaging system was used to take photos and record the staining results. Image Pro Plus 6.0 software (Media Cybernetics, Inc., Rockville, MD, USA) was used to analyze the average optical density of the positive results of oil red staining.

2.3. RNA-Seq and Bioinformatics Analysis

Liver, abdominal fat, and ovarian stroma were selected from 3 individuals in each group. The Trizol kit (Invitrogen, Santa Clara, CA, USA) was used to extract the total RNA of liver, abdominal fat, and ovarian stroma according to the manufacturer’s instructions. The RNA integrity was determined by an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). The RNA samples were used for library construction. The mRNA libraries were sequenced by Novogene Co., Ltd. (Beijing, China) using Nova-PE150 (Illumina, San Diego, CA, USA). The clean reads were obtained after the filtration of low-quality reads using standard quality control by FastaQC software. Clean reads were mapped to the Anser cygnoides domestication reference genome (data being published) using the HISAT2 (version 2.2.1) software [10]. The output SAM (sequencing alignment/mapping) file was converted to a BAM (binary alignment/mapping) file and sorted using SAMtools (version 1.10) [11]. Subsequently, the expression of each transcript was calculated by featureCounts (version 1.6.0) [12]. DEseq2 was used to identify the different expression genes (DEGs) among different groups, the screening criteria were |log2Foldchange| > 1, p-adjust < 0.05. Functional analysis used KOBAS 3.0 online (http://kobas.cbi.pku.edu.cn/kobas3/?t=1, accessed on 26 August 2021) [13].

2.4. Quantitative Real-Time PCR Validation

Total RNA extracted from the tissues was reverse transcribed into cDNA using a Taq Pro Universal SYBR qPCR Master Mix (Vazyme Co., Ltd., Nanjing, China). Primer 5.0 was used to design the primers (Table 1). A BLAST search against the reference genome was then carried out to confirm that primers were specific for the intended target genes. Taq Pro Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) and a Bio-Rad CFX96 real-time PCR detection system (Bio-Rad, Hercules, CA, USA) were used for RT-PCR, and each sample was assayed three times. β-actin and GAPDH were used as housekeeping genes. The 2−∆∆CT method was used for normalization of the qPCR results, after which the normalized data were used for statistical analysis, and p < 0.05 was considered significantly different.

2.5. Statistical Analysis

The abdominal fat weight/index, liver weight/index, fat count of liver, and expression levels of DEGs were expressed as the mean ± standard deviation (SD). One-way ANOVA was used for statistical analysis in this study. p < 0.05 was considered statistically significant. All statistical analyses were carried out using the SPSS 27.0 software.

3. Results

3.1. The Lipid Deposition Patterns of Geese with Different Laying Performance

There was no significant difference in body weight among the three groups (Figure 1A). However, the abdominal fat weight and abdominal fat index of geese in the low- and no-egg production groups were significantly (p < 0.05) higher than those in the high-egg production group (Figure 1B). In contrast to the weight of abdominal fat, the content of lipids in the liver of the high-egg production group was significantly (p < 0.05) higher than that of the low- and no-egg production groups.

3.2. Overview of the mRNA Transcriptome with Different Egg Production Performance Geese

A total of 672,265,178 raw reads were obtained from 27 samples through mRNA sequencing, and 90.65% of the clean reads were aligned to the goose reference genome (Supplementary Table S1). As shown in Figure 2A, the number of DEGs in ovarian stroma with different egg performances was lower than that of liver and abdominal fat. Meanwhile, the number of unique DEGs between the no- and high-egg production group was the largest in the three tissues (Figure 2B–D).

3.3. Functional Analysis of DEGs among Different Egg Production Performances in Abdominal Fat, Liver, and Ovarian Stroma

The DEGs in almost groups were significantly enriched in several KEGG pathways (Supplementary Table S2). Of note, the KEGG pathways related to lipid metabolism (e.g., cholesterol metabolism, fat digestion and absorption, bile secretion) were enriched in all pairwise comparisons of abdominal fat and liver (Figure 3). However, some KEGG pathways related to inflammation (IL-17 signaling pathway, NF-kappa B signaling pathway, inflammatory mediator regulation of TRP channels, leukocyte transendothelial migration, Toll-like receptor signaling pathway, inflammatory bowel disease) and the circadian rhythm pathway were enriched by DEGs only in abdominal fat and liver, respectively. Meanwhile, the hormone-related KEGG pathways (thyroid hormone synthesis, GnRH signaling pathway, insulin signaling pathway, and estrogen signaling pathway) were enriched by DEGs in abdominal fat. The DEGs in ovarian stroma enriched some KEGG pathways related to ovarian steroidogenesis (ovarian steroidogenesis and oxytocin signaling pathway) and cell adhesion (focal adhesion, ECM-receptor interaction, and adherens junction).

3.4. Network Construction of Liver, Abdominal Fat, and Ovarian Stroma Regulating Laying Performance in Geese

We chose key DEGs in liver (enriched in lipid metabolism and circadian rhythm), abdominal fat (enriched in inflammatory response, lipid metabolism, and hormone), and ovarian stroma (enriched in steroid hormone and cell adhesion) to construct a protein–protein interaction network (PPI) (Figure 4A). It was found that the DEGs in liver, abdominal fat, and ovarian stroma have a mutual regulatory relationship through the PPI. DEGs with more than 10 nodes were used for qPCR verification in this network (Figure 4B). As shown in Figure 4B, expression of almost all these selected DEGs displayed changes in the same direction with those observed using RNA-seq, indicating the true reliability of our sequencing and analysis methods.

4. Discussion

Phenotypically, we found that liver lipid content and abdominal fat weight had opposite trends in the high-, low-, and no-egg production groups. Meanwhile, a previous study of broilers showed a positive relationship between liver lipid content and abdominal fat weight [14,15] This difference between chicken and geese may be due to the unique lipid metabolism pattern of geese (domesticated by migratory birds and high lipid storage capacity of liver) [16]. The abdominal fat weight of the no-egg production group was higher than that of the other groups. Consistent with that, the previous results showed that excessive abdominal fat delays the start of production in poultry [17], and that excessive abdominal fat reduces the egg production performance of female poultry [18]. Our results also confirm this point.
Through the transcriptome results, we found that the number of DEGs in liver and abdominal fat was greater than that in ovarian stroma. These results showed that energy metabolism organs play a vital role in female poultry reproduction activities [19,20]. The liver and fat are the most important tissues of lipid metabolism [21,22]. Meanwhile, we found that the KEGG pathways related to lipid metabolism were enriched in all pairwise comparisons of abdominal fat and liver by the functional analysis of DEGs, while the circadian rhythm pathway was enriched by DEGs only in the liver. Circadian rhythms are physiologic and behavioral cycles that control a variety of biological processes, including feeding, the sleep–wake cycle, and the female reproductive cycle [23]. Many studies have shown that circadian rhythms in the ovary is important for reproduction [24,25]. However, our study first showed that the circadian rhythms in liver is related in poultry reproduction. The dynamic balance of liver energy depends on enzymes, and these rate-limiting enzymes are expressed in a circadian rhythm [26,27]. In view of the circadian rhythm of the liver in mammals, Reinke et al. reviewed its functions in carbohydrate metabolism, lipid metabolism, amino acid metabolism, detoxification, synthesis of plasma proteins, and bile acid metabolism [28]. The important role of NPAS2 [29,30], RORG [31], and PER3 [32] in the circadian regulation of lipid metabolism has been repeatedly demonstrated. Studies in mammals have shown that NPAS2 not only activates the expression of PER but is also the target gene of ROR [33]. In our study, NPAS2, RORG, and PER3 were differentially expressed in the liver of geese among different laying groups. Hence, we speculate that the difference of circadian rhythm in the liver may cause the differences of lipid metabolism.
Meanwhile, we found that some KEGG pathways related to inflammation, including the IL-17 signaling pathway, NF-kappa B signaling pathway, inflammatory mediator regulation of TRP channels, leukocyte transendothelial migration, Toll-like receptor signaling pathway, and inflammatory bowel disease KEGG pathways, were enriched by DEGs in abdominal fat. The positive relationship between obesity and inflammation has been widely established [34,35]. In our results, multiple inflammation-related genes (IL8 and CCL4) showed a trend of low to high expression in the abdominal fat of the high-, low-, and no-egg production groups. The negative effects of inflammation are multifaceted, including for reproduction [36]. Adipose tissue has been considered as an endocrine organ capable of secreting hormones that travel through the bloodstream to reach their target tissues [37]. In addition, the DEGs in abdominal fat are also enriched in the hormone-related KEGG pathways. Therefore, the change in hormone-related gene expression in abdominal fat may directly affect the laying performance of geese.
The ovarian stroma is directly related to follicle development and laying, and DEGs in ovarian stroma are enriched in the KEGG pathways associated with steroid hormones and cell adhesion. Lipids are not only important precursors of steroid hormones but are also important components of cell membranes. Steroid hormones play a crucial role in both the initiation of egg production and the level of egg production [38,39]. Meanwhile, cell adhesion plays an important role in the selection and development of poultry follicles [40]. Our results support that steroid hormones and cell adhesion determines the laying performance of geese.

5. Conclusions

In conclusion, we found that the liver fat content of the high-, low-, and no-egg production group decreased in turn, while the abdominal fat weight increased in turn by phenotypic measurement. For transcriptional regulation, the number of DEGs in liver and abdominal fat with different laying performance was greater than that in ovarian stroma. Through functional analysis, the KEGG pathways related to lipid metabolism were enriched in all pairwise comparisons of abdominal fat and liver. 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 enriched some KEGG pathways related to ovarian steroidogenesis and cell adhesion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12141775/s1, Table S1: Basic information of mRNA sequencing data of all samples in this study; Table S2: KEGG pathways enriched by DEGs of different groups.

Author Contributions

Q.O., S.H. and J.W. conceptualized and designed this study. Q.O. and S.H. performed the main experiments and analyzed the data. B.T., B.H., J.H., H.H. and L.L. participated in experimental animal management, tissue sampling, and data analysis. Q.O. and S.H. drafted this manuscript. J.W. reviewed this manuscript. J.W. acquired the funding and supervised this study. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the China Agriculture Research System of MOF and MARA (CARS-42-4) and the Key Technology Support Program of Sichuan Province (2021YFYZ0014).

Institutional Review Board Statement

All geese were obtained from the Waterfowl Breeding Experimental Farm of Sichuan Agricultural University. All experimental procedures involved in animal manipulation were approved by the Institutional Animal Care and Use Committee (IACUC) of Sichuan Agricultural University (Chengdu Campus, Sichuan, China) under approval no. 20180034.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original sequencing data for this study can be found in the Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra; accessed on 22 January 2022) at NCBI with the BioProject ID: PRJNA747742.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Gao, S.; Gan, X.; He, H.; Hu, S.; Deng, Y.; Chen, X.; Li, L.; Hu, J.; Li, L.; Wang, J. Dynamic characteristics of lipid metabolism in cultured granulosa cells from geese follicles at different developmental stages. Biosci. Rep. 2019, 39, BSR20192188. [Google Scholar] [CrossRef] [PubMed]
  2. Li, Q.; Hu, S.; Wang, Y.; Deng, Y.; Yang, S.; Hu, J.; Li, L.; Wang, J. mRNA and miRNA Transcriptome Profiling of Granulosa and Theca Layers From Geese Ovarian Follicles Reveals the Crucial Pathways and Interaction Networks for Regulation of Follicle Selection. Front. Genet. 2019, 10, 988. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wen, R.; Gan, X.; Hu, S.; Gao, S.; Deng, Y.; Qiu, J.; Sun, W.; Li, L.; Han, C.; Hu, J.; et al. Evidence for the existence of de novo lipogenesis in goose granulosa cells. Poult. Sci. 2019, 98, 1023–1030. [Google Scholar] [CrossRef] [PubMed]
  4. Hermier, D. Lipoprotein metabolism and fattening in poultry. J. Nutr. 1997, 127, 805s–808s. [Google Scholar] [CrossRef]
  5. Wang, G.; Jin, L.; Li, Y.; Tang, Q.; Hu, S.; Xu, H.; Gill, C.A.; Li, M.; Wang, J. Transcriptomic analysis between Normal and high-intake feeding geese provides insight into adipose deposition and susceptibility to fatty liver in migratory birds. BMC Genom. 2019, 20, 372. [Google Scholar] [CrossRef] [Green Version]
  6. Milisits, G.; Szentirmai, E.; Donko, T.; Budai, Z.; Sütő, Z. Effect of starting body fat content and genotype of laying hens on the changes in their live weight, body fat content, egg production and egg composition during the first egg-laying period. Br. Poult. Sci. 2015, 56, 666–672. [Google Scholar] [CrossRef]
  7. Bornstein, S.; Plavnik, I.; Lev, Y. Body weight and/or fatness as potential determinants of the onset of egg production in broiler breeder hens. Br. Poult. Sci. 1984, 25, 323–341. [Google Scholar] [CrossRef]
  8. Wei, Z.; Li, P.; Huang, S.; Lkhagvagarav, P.; Zhu, M.; Liang, C.; Jia, C. Identification of key genes and molecular mechanisms associated with low egg production of broiler breeder hens in ad libitum. BMC Genom. 2019, 20, 408. [Google Scholar] [CrossRef]
  9. Ouyang, Q.; Hu, S.; Wang, G.; Hu, J.; Zhang, J.; Li, L.; Hu, B.; He, H.; Liu, H.; Xia, L.; et al. Comparative Transcriptome Analysis Suggests Key Roles for 5-Hydroxytryptamlne Receptors in Control of Goose Egg Production. Genes 2020, 11, 455. [Google Scholar] [CrossRef] [Green Version]
  10. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [Green Version]
  11. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  12. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  13. Xie, C.; Mao, X.; Huang, J.; Ding, Y.; Wu, J.; Dong, S.; Kong, L.; Gao, G.; Li, C.Y.; Wei, L. KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39, W316–W322. [Google Scholar] [CrossRef] [Green Version]
  14. Liang, M.J.; Wang, Z.P.; Xu, L.; Leng, L.; Li, H. Estimating the genetic parameters for liver fat traits in broiler lines divergently selected for abdominal fat. Genet. Mol. Res. 2015, 14, 9646–9654. [Google Scholar] [CrossRef]
  15. Leng, L.; Zhang, H.; Dong, J.Q.; Wang, Z.P.; Zhang, X.Y.; Wang, S.Z.; Cao, Z.P.; Li, Y.M.; Li, H. Selection against abdominal fat percentage may increase intramuscular fat content in broilers. J. Anim. Breed. Genet. 2016, 135, 422–428. [Google Scholar] [CrossRef]
  16. Lu, L.; Chen, Y.; Wang, Z.; Li, X.; Chen, W.; Tao, Z.; Shen, J.; Tian, Y.; Wang, D.; Li, G. The goose genome sequence leads to insights into the evolution of waterfowl and susceptibility to fatty liver. Genome Biol. 2015, 16, 89. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, X.Y.; Wu, M.Q.; Wang, S.Z.; Zhang, H.; Du, Z.Q.; Li, Y.M.; Cao, Z.P.; Luan, P.; Leng, L.; Li, H. Genetic selection on abdominal fat content alters the reproductive performance of broilers. Anim. Int. J. Anim. Biosci. 2017, 12, 1232–1241. [Google Scholar] [CrossRef]
  18. Pan, Y.E.; Liu, Z.C.; Chang, C.J.; Huang, Y.F.; Lai, C.Y.; Walzem, R.L.; Chen, S.E. Feed restriction ameliorates metabolic dysregulation and improves reproductive performance of meat-type country chickens. Anim. Reprod. Sci. 2014, 151, 229–236. [Google Scholar] [CrossRef]
  19. Sauerwein, H.; Bendixen, E.; Restelli, L.; Ceciliani, F. The adipose tissue in farm animals: A proteomic approach. Curr. Protein Pept. Sci. 2014, 15, 146–155. [Google Scholar] [CrossRef]
  20. Gloux, A.; Duclos, M.J.; Brionne, A.; Bourin, M.; Nys, Y.; Réhault-Godbert, S. Integrative analysis of transcriptomic data related to the liver of laying hens: From physiological basics to newly identified functions. BMC Genom. 2019, 20, 821. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, G.; Kim, W.K.; Cline, M.A.; Gilbert, E.R. Factors affecting adipose tissue development in chickens: A review. Poult. Sci. 2017, 96, 3687–3699. [Google Scholar] [CrossRef]
  22. Zaefarian, F.; Abdollahi, M.R.; Cowieson, A.; Ravindran, V. Avian Liver: The Forgotten Organ. Animals 2019, 9, 63. [Google Scholar] [CrossRef] [Green Version]
  23. Loh, D.H.; Kuljis, D.A.; Azuma, L.; Wu, Y.; Truong, D.; Wang, H.B.; Colwell, C.S. Disrupted reproduction, estrous cycle, and circadian rhythms in female mice deficient in vasoactive intestinal peptide. J. Biol. Rhythm. 2014, 29, 355–369. [Google Scholar] [CrossRef] [Green Version]
  24. Di Cara, F.; King-Jones, K. The Circadian Clock Is a Key Driver of Steroid Hormone Production in Drosophila. Curr. Biol. 2016, 26, 2469–2477. [Google Scholar] [CrossRef] [Green Version]
  25. Tao, Z.; Song, W.; Zhu, C.; Xu, W.; Liu, H.; Zhang, S.; Huifang, L. Comparative transcriptomic analysis of high and low egg-producing duck ovaries. Poult. Sci. 2017, 96, 4378–4388. [Google Scholar] [CrossRef]
  26. Panda, S.; Antoch, M.P.; Miller, B.H.; Su, A.I.; Schook, A.B.; Straume, M.; Schultz, P.G.; Kay, S.A.; Takahashi, J.S.; Hogenesch, J.B. Coordinated transcription of key pathways in the mouse by the circadian clock. Cell 2002, 109, 307–320. [Google Scholar] [CrossRef] [Green Version]
  27. Storch, K.F.; Lipan, O.; Leykin, I.; Viswanathan, N.; Davis, F.C.; Wong, W.H.; Weitz, C.J. Extensive and divergent circadian gene expression in liver and heart. Nature 2002, 417, 78–83. [Google Scholar] [CrossRef]
  28. Reinke, H.; Asher, G. Circadian Clock Control of Liver Metabolic Functions. Gastroenterology 2016, 150, 574–580. [Google Scholar] [CrossRef] [Green Version]
  29. Lee, S.M.; Zhang, Y.; Tsuchiya, H.; Smalling, R.; Jetten, A.M.; Wang, L. Small heterodimer partner/neuronal PAS domain protein 2 axis regulates the oscillation of liver lipid metabolism. Hepatology 2015, 61, 497–505. [Google Scholar] [CrossRef] [Green Version]
  30. Dudley, C.A.; Erbel-Sieler, C.; Estill, S.J.; Reick, M.; Franken, P.; Pitts, S.; McKnight, S.L. Altered patterns of sleep and behavioral adaptability in NPAS2-deficient mice. Science 2003, 301, 379–383. [Google Scholar] [CrossRef]
  31. Hojo, H.; Enya, S.; Arai, M.; Suzuki, Y.; Nojiri, T.; Kangawa, K.; Koyama, S.; Kawaoka, S. Remote reprogramming of hepatic circadian transcriptome by breast cancer. Oncotarget 2017, 8, 34128–34140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. van der Veen, D.R.; Laing, E.E.; Bae, S.E.; Johnston, J.D.; Dijk, D.J.; Archer, S.N. A Topological Cluster of Differentially Regulated Genes in Mice Lacking PER3. Front. Mol. Neurosci. 2020, 13, 15. [Google Scholar] [CrossRef] [PubMed]
  33. Crumbley, C.; Wang, Y.; Kojetin, D.J.; Burris, T.P. Characterization of the core mammalian clock component, NPAS2, as a REV-ERBalpha/RORalpha target gene. J. Biol. Chem. 2010, 285, 35386–35392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Izaola, O.; de Luis, D.; Sajoux, I.; Domingo, J.C.; Vidal, M. Inflammation and obesity (lipoinflammation). Nutr. Hosp. 2015, 31, 2352–2358. [Google Scholar] [CrossRef] [Green Version]
  35. Karczewski, J.; Śledzińska, E.; Baturo, A.; Jończyk, I.; Maleszko, A.; Samborski, P.; Begier-Krasińska, B.; Dobrowolska, A. Obesity and inflammation. Eur. Cytokine Netw. 2018, 29, 83–94. [Google Scholar] [CrossRef]
  36. Mattina, G.F.; van Lieshout, R.J.; Steiner, M. Inflammation, depression and cardiovascular disease in women: The role of the immune system across critical reproductive events. Ther. Adv. Cardiovasc. Dis. 2019, 13, 1753944719851950. [Google Scholar] [CrossRef]
  37. Frigolet, M.E.; Gutiérrez-Aguilar, R. The colors of adipose tissue. Gac. Med. Mex. 2020, 156, 142–149. [Google Scholar] [CrossRef] [Green Version]
  38. Yu, J.; Lou, Y.; Zhao, A. Transcriptome analysis of follicles reveals the importance of autophagy and hormones in regulating broodiness of Zhedong white goose. Sci. Rep. 2016, 6, 36877. [Google Scholar] [CrossRef] [Green Version]
  39. Brady, K.; Liu, H.C.; Hicks, J.A.; Long, J.A.; Porter, T.E. Transcriptome Analysis During Follicle Development in Turkey Hens with Low and High Egg Production. Front. Genet. 2021, 12, 619196. [Google Scholar] [CrossRef]
  40. Zhu, G.; Mao, Y.; Zhou, W.; Jiang, Y. Dynamic Changes in the Follicular Transcriptome and Promoter DNA Methylation Pattern of Steroidogenic Genes in Chicken Follicles throughout the Ovulation Cycle. PLoS ONE 2015, 10, e0146028. [Google Scholar] [CrossRef]
Figure 1. The lipid deposition patterns of geese with different laying performance. The body weight (A) and abdominal fat weight (index) (B) of different laying performance geese. (C) From left to right, the liver morphology and oil red staining sections of the high-egg production group, low-egg production group and no-egg group were in turn. Meanwhile, the OD values of oil red staining sections were counted. Data are shown as the mean ± standard deviation (SD). * means p < 0.05, and ** means p < 0.01.
Figure 1. The lipid deposition patterns of geese with different laying performance. The body weight (A) and abdominal fat weight (index) (B) of different laying performance geese. (C) From left to right, the liver morphology and oil red staining sections of the high-egg production group, low-egg production group and no-egg group were in turn. Meanwhile, the OD values of oil red staining sections were counted. Data are shown as the mean ± standard deviation (SD). * means p < 0.05, and ** means p < 0.01.
Animals 12 01775 g001
Figure 2. The number of DEGs in different groups. (A) Histogram of the number of DEGs in different groups. Venn diagram of the common DEGs between three pairwise comparisons in abdominal fat (B), liver (C), and ovarian stroma (D). “AF” means abdominal fat, “L” means liver, “O” means ovarian stroma.
Figure 2. The number of DEGs in different groups. (A) Histogram of the number of DEGs in different groups. Venn diagram of the common DEGs between three pairwise comparisons in abdominal fat (B), liver (C), and ovarian stroma (D). “AF” means abdominal fat, “L” means liver, “O” means ovarian stroma.
Animals 12 01775 g002
Figure 3. Functional annotation of the DEGs between different performances of geese in abdominal fat (A), liver (B), and ovarian stroma (C). “AF” means abdominal fat, “L” means liver, “O” means ovarian stroma.
Figure 3. Functional annotation of the DEGs between different performances of geese in abdominal fat (A), liver (B), and ovarian stroma (C). “AF” means abdominal fat, “L” means liver, “O” means ovarian stroma.
Animals 12 01775 g003
Figure 4. The network of liver, abdominal fat, and ovarian stroma regulating laying performance in geese and the expression validation of hub genes. The network of L–AF–O (liver, abdominal fat, ovarian stroma) regulating laying performance in geese (A). qRT-PCR validation of expression of the main DEGs involving the network of egg production performance regulation by liver, abdominal fat, and ovarian stroma (B). The results of qRT-PCR and RNA-Seq were expressed as the mean ± SD (n = 3, three biological replications at each group). “AF” means abdominal fat, “L” means liver.
Figure 4. The network of liver, abdominal fat, and ovarian stroma regulating laying performance in geese and the expression validation of hub genes. The network of L–AF–O (liver, abdominal fat, ovarian stroma) regulating laying performance in geese (A). qRT-PCR validation of expression of the main DEGs involving the network of egg production performance regulation by liver, abdominal fat, and ovarian stroma (B). The results of qRT-PCR and RNA-Seq were expressed as the mean ± SD (n = 3, three biological replications at each group). “AF” means abdominal fat, “L” means liver.
Animals 12 01775 g004
Table 1. Primers used for qRT-PCR in this study.
Table 1. Primers used for qRT-PCR in this study.
Primer NameSequence (5′–3′)Product Length (bp)
GAPDH-FGCTGATGCTCCCATGTTCGTGAT86
GAPDH-RGTGGTGCAAGAGGCATTGCTGAC
β-ACTIN-FCAACGAGCGGTTCAGGTGT92
β-ACTIN-RTGGAGTTGAAGGTGGTCTCGT
RORG-FTGTGCCAGAACGACCAGAT102
RORG-RAGAGGACGGTCCGGTTGT
PER3-FGAGCAGTGCCTTTGTTGGGT276
PER3-RTCAGAGGGCTTGTTCGGACT
NPAS2-FTCACAGAGCACCACCGATTA148
NPAS2-RATAGCAACACGACTTCCCCT
NR1H4-FGCCTCAGATTTCATCGCCAC228
NR1H4-RGCTTTGTCACCACAGACCACG
LCAT-FCAGCGTGTCTTCCTCATTGC187
LCAT-RACATAAGTGGGATGCCCTGAT
DGAT1-FGCCTACCCCGACAACCTCAC180
DGAT1-RCACCATCCACTGCTGGATCA
IL8-FCCTGGTAAGGATGGGAAACG168
IL8-RGGGTCCAAGCACACCTCTCT
CCL4-FATGAAGGTCTCTGTGGCTGC119
CCL4-RTCCCGTTGGATGTAGGTGAA
DGAT2-FACCCACAATCTGCTGACCAC239
DGAT2-RGATAAGATGTAGTCTATGCTGTCGC
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ouyang, Q.; Hu, S.; Tang, B.; Hu, B.; Hu, J.; He, H.; Li, L.; Wang, J. Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese. Animals 2022, 12, 1775. https://doi.org/10.3390/ani12141775

AMA Style

Ouyang Q, Hu S, Tang B, Hu B, Hu J, He H, Li L, Wang J. Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese. Animals. 2022; 12(14):1775. https://doi.org/10.3390/ani12141775

Chicago/Turabian Style

Ouyang, Qingyuan, Shenqiang Hu, Bincheng Tang, Bo Hu, Jiwei Hu, Hua He, Liang Li, and Jiwen Wang. 2022. "Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese" Animals 12, no. 14: 1775. https://doi.org/10.3390/ani12141775

APA Style

Ouyang, Q., Hu, S., Tang, B., Hu, B., Hu, J., He, H., Li, L., & Wang, J. (2022). Comparative Transcriptome Analysis Provides Novel Insights into the Effect of Lipid Metabolism on Laying of Geese. Animals, 12(14), 1775. https://doi.org/10.3390/ani12141775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop