Next Article in Journal
Mutations Affecting HVO_1357 or HVO_2248 Cause Hypermotility in Haloferax volcanii, Suggesting Roles in Motility Regulation
Previous Article in Journal
Strawberry FaWRKY25 Transcription Factor Negatively Regulated the Resistance of Strawberry Fruits to Botrytis cinerea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 1
10.3390/genes12010057
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production

by
Xiuyu Bao
,
Yiping Song
,
Tao Li
,
Shanshan Zhang
,
Lihua Huang
,
Shuya Zhang
,
Junting Cao
,
Xiaolin Liu
and
Jianqin Zhang
*
College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
*
Author to whom correspondence should be addressed.
Genes 2021, 12(1), 57; https://doi.org/10.3390/genes12010057
Submission received: 19 November 2020 / Revised: 25 December 2020 / Accepted: 28 December 2020 / Published: 31 December 2020
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
The egg-laying rate is an important indicator for evaluating fertility of poultry. In order to better understand the laying mechanism of Muscovy ducks, gene expression profiles and pathways of ovarian tissues in high- and low-laying black (BH and BL) and white Muscovy ducks (WH and WL) during the peak production period were performed by using RNA-seq. The total number of reads produced for each ovarian sample ranged from 44,344,070 to 47,963,328. A total of 113, 619 and 87 differentially expressed genes (DEGs) were identified in BH-vs-WH, BL-vs-BH and BL-vs-WL, respectively. Among them, 54, 356 and 49 genes were up regulated and 59, 263 and 38 genes were down regulated. In addition, there were only 10 up-regulated genes in WL-vs-WH. In the comparison of DEGs in black and white Muscovy ducks, two co-expressed DEG genes were detected between BH-vs-WH and BL-vs-WL and seven DEGs were co-expressed between BL-vs-BH and WL-vs-WH. The RNA-Seq data were confirmed to be reliable by qPCR. Numerous DEGs known to be involved in ovarian development were identified, including TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Gene Ontology (GO) annotations indicated that DEGs related to ovarian development were mainly enriched in biological processes of “circadian sleep/wake cycle process,” “negative regulation of transforming growth factor-β secretion,” “positive regulation of calcium ion transport” in BH-vs-WH and “cell surface receptor signaling pathway,” “Notch signaling pathway” and “calcium ion transport” in BL-vs-BH. Besides, “steroid biosynthetic process,” “granulosa cell development” and “egg coat formation” were mainly enriched in BL-vs-WL and “reproduction,” “MAPK cascade” and “mitotic cell cycle” were mainly enriched in WL-vs-WH. KEGG pathway analysis showed that the PI3K-Akt signaling pathway and ovarian steroidogenesis were the most enriched in Muscovy duck ovary transcriptome data. This work highlights potential genes and pathways that may affect ovarian development in Muscovy duck.

1. Introduction

Reproductive performance, an important economic trait of poultry, is mainly determined by the egg production. Improving egg production ability is an important breeding objective. With traditional breeding methods, the reproductive performance of laying duck has been gradually improved but more significant improvements have been slow [1]. Ovary is a female reproductive organ, which is closely related to egg-laying characteristics and the ovarian function of poultry directly affects the egg production [2]. Ovarian follicular development is a complex process, involving a variety of endocrine, autocrine and paracrine factors, which control the proliferation and differentiation of oocytes, granulosa cells and thecal cells [3]. The number of primordial follicle pools in the ovary is fixed and the quiescence, survival and activation of follicles in the pools depends on the dynamic balance, which is regulated by many signaling molecules or pathways [4,5]. Genetic or environmental factors affect follicular development, which may lead to follicular inactivation or premature activation and even degeneration and atresia at different stages of development [6]. Therefore, the regulation of ovarian function, including the identification and analysis of reproductive-related genes, has been concerned to improve the reproductive performance of poultry. For example, transforming growth factor-β family (TGF-β), fibroblast growth factors (FGFs), insulin-like growth factor family (IGFs) [7], FSHR [8], GnRH [9], STAR [10] and PRLR [11] have been shown to be associated with poultry ovarian function and affect egg production.
Although candidate genes for avian reproduction have been identified, the specific molecular mechanisms and related signaling pathways are poorly understood. As a next-generation sequencing technology, RNA-seq is used to explore the function of genes by detecting the gene expression in animal samples [12]. In recent years, the ovarian transcriptome of livestock, such as pig [13], cattle [14], sheep [15], has been studied by using RNA-seq. Similarly, there are many studies on the transcriptome of poultry ovary. Zhang et al. (2019) found that there were 305 differentially expressed genes (DEGs) identified in the ovaries of hens with relatively greater and lesser egg production. Besides, five candidate genes related to egg production, including ZP2, WNT4, AMH, IGF1 and CYP17A1 genes, were discovered [16]. Tao et al. (2017) performed a comparative transcriptomic analysis of the ovaries of Jinding ducks with high and low egg production using RNA-seq, the results showed that 367 of 843 DEGs were down regulated in high egg production ovaries and 476 DEGs were up-regulated in low egg production ovaries [17]. Luan et al. (2014) identified the DEGs in the ovaries of Huoyan geese between the laying period and the ceased period, 344 up- and 344 down-regulated genes were classified into functional categories and 12 DEGs that were mainly involved in pathways for reproduction regulation [18]. Xu et al. (2016) compared the gene expression of pigeon ovary at different stages, there were 409 DEGs between pre- and post-ovulation. However, there are few reports on the high- and low-yield ovarian transcriptome of Muscovy duck [19].
Muscovy duck, known for its strong adaptability, fertility and high meat yield [20], is native to tropical areas of Central and South America. In China, black Muscovy and white Muscovy are mainly raised. The egg-laying peak period of Muscovy duck is between 35 and 53 weeks, covering most of the egg-laying period [21]. Although these ducks have been raised on a large scale in China, the poor egg production performance and the lack of systematic breeding have resulted in uneven production performance in Muscovy duck and seriously affects the economic benefits of many enterprises. Under the same conditions of nutrition and feeding management, the egg production performance of Muscovy duck is mainly determined by genes and related regulatory factors. Therefore, the purpose of our research was to identify the DEGs and pathways and analyze the similarities and differences of the differential expression in the ovarian tissues of Muscovy ducks with high- and low-yield by using RNA-seq. Through a comprehensive analysis of gene expression levels reflecting the ovaries of black and white Muscovy ducks, the genes and molecular mechanisms involved in this process have been found. All results provide theoretical basis for the breeding of Muscovy ducks and are useful for understanding of the molecular regulation mechanism of egg production characteristics in black and white Muscovy duck.

2. Materials and Methods

2.1. Animal and Sample Collection

The egg production of Muscovy ducks aged 216 to 280 days (65 days in total) was recorded. More than 60 eggs were considered as high-laying black Muscovy duck (BH) and white Muscovy duck (WH) and 15–20 eggs were defined as low-yield duck (BL and WL). Six black Muscovy ducks (3 high-yield and 3 low-yield) and 6 white Muscovy ducks (half each for high- and low-yield laying ducks), raised under the same environmental conditions with free access to feed and water (Table 1), were randomly selected and slaughtered quickly to collect the ovarian tissues, which comprised the whole ovary including the small and large yellow follicles. The samples were immediately frozen in liquid nitrogen and stored at −80 °C until further use. Total RNA was extracted from ovarian tissue of Muscovy duck using TRIzol reagent according to the manufacturer’s instructions (Takara, Dalian, China). The purity and concentration of RNA were detected by agarose gel electrophoresis and Nanodrop 2000 (Thermo, Waltham, MA, USA) and RNA integrity was measured using the 2100 Bioanalyzer (Agilent Technologies, San Jose, CA, USA). All samples had an RNA integrity number (RIN) ≥ 7.8 and a RNA concentration ≥ 875.6 ng/μL (Table S1). The tested Muscovy ducks were purchased from Anda Farm in Tongguan County, Shaanxi Province, China and raised according to the standard procedure. Besides this, the Muscovy ducks were known to be unrelated to each other for at least three generations. All animal procedures were performed according to guidelines provided by the China Council on Animal Care and the protocols were approved by the Experimental Animal Management Committee (EAMC) of Northwest A&F University (ethic code: #0726/2018).

2.2. Library Construction and Sequencing

Twelve RNA samples with high quality (average RIN: 8.4 ± 0.251) and concentration (average concentration: 1136.9 ± 3.828 ng/μL) were used to construct the sequencing libraries and high-throughput sequencing was performed on the Illumina novaseq 6000 platform (San Diego, CA, USA) following the manufacturer’s recommendations, generating 150 bp paired-end reads. The high-quality clean reads were obtained by filtering the raw reads and removing: (1) the sequences containing adapters; (2) the sequences with more than 10% of N bases; (3) the sequences with more than 50% base quality values less than 10. The clean reads were aligned to the Anas platyrhynchos genome (https://www.ncbi.nlm.nih.gov/genome/?term=duck) and annotated transcripts (https://www.ncbi.nlm.nih.gov/assembly/GCF_003850225.1) by way of the Tophat 2 software (http://ccb.jhu.edu/software/hisat2/index.shtml).

2.3. Differential Expression Analysis

The number of clean reads for each gene was calculated and the Fragment Per Kilobases per Million Fragment (FPKM) was used to estimate the expression abundance of transcripts from different samples [22]. Differential expression analyses of the 4 comparison groups (BL-vs-BH, WL-vs-WH, BH-vs-WH and BL-vs-WL) were performed using the DESeq R package (http://www.bioconductor.org/packages/release/bioc/html/DESeq.html) [23,24] and genes with an p-value ≤ 0.05 and an expression |log2 Fold | ≥ 1were identified as DEGs.

2.4. GO and KEGG Analysis

To annotate the function of these DEGs, Gene Ontology (GO) analysis was conducted by using the GOseq software (https://bioconductor.org/packages/release/bioc/html/goseq.html) for each of the three main categories: biological process, cellular component and molecular function. Biological pathways enriched for the identified DEGs through KEGG pathway analyses were carried out using the KOBAS software (http://kobas.cbi.pku.edu.cn/kobas3/?t=1).

2.5. Verification of RNA-seq by qPCR

To verify the reliability of the RNA-Seq results, 4 DEGs (TGFβ2, FGF18, POSTN, SMOC1) in the ovary of black Muscovy duck and 3 DEGs (POSTN, SMOC1, COL4A1) in white Muscovy duck were randomly selected, respectively. The β-actin gene of duck was used as the control (NC_040060.1), after the stability of β-actin was tested in this system. The RNA extracted from same ovaries with an RIN ≥ 7.8, A260/280 >1.8 and A260/230 > 2.0 were used for cDNA synthesis according to the manufacturer’s instructions (Takara, Dalian, China). The primers information were listed in Table 2. qPCR was run in triplicate and performed in a EcoRT48 system (OSA, London, UK), with a 10 μL reaction volume containing 5 μL SYBR Premix Ex Taq II (Takara, Dalian, China), 0.8 µL cDNA (750 ng/µL), 0.2 µL of each primer and 3.8 µL RNase-free ddH2O. The qPCR program was 95 °C for 30 s; followed by 40 cycles of 95 °C for 5 s, 60 °C for 30 s and then 95 °C for 15 s, 55 °C for 15 s, 95 °C for 15 s. Gene expression results were calculated using the 2−ΔΔCt method.

2.6. Statistical Analysis

Experiments were repeated 3 times and presented as mean ± standard deviation and conducted using one-way ANOVA with Dunnet’s t-test at p < 0.05 probability levels in SPSS18.0. “*” was considered significant difference (p < 0.05); “**” was considered extremely significant difference (p < 0.01).

3. Results

3.1. High-Throughput Sequencing and Read Mapping

In this study, a total of 12 libraries in high- and low-yield of ovarian tissue in Muscovy ducks were established by high-throughput RNA sequencing. The clean reads of each sample were more than 44 million, the mapped reads and unique mapped reads ranged from 60.28% to 67.40% and 57.81% to 64.36%, respectively. In addition, more than 16,276 genes were detected in each sample (Table 3)

3.2. Differentially Expressed Genes

According to the correlation analysis of 12 samples, the differences between biological replicates were small and the repeatability was high, which indicated that the selection of experimental samples was consistent and reliable (Figure S1). A comparison of the gene expression showed that a total of 113 genes were differentially expressed between BH and WH, including 54 up-regulated genes and 59 down-regulated genes. 619 DEGs were discovered in BL-vs-BH and 356 DEGs were up-regulated genes and 263 DEGs were down-regulated genes. In BL-vs-WL, there were 87 DEGs, of which 49 were up-regulated and 38 were down-regulated genes. In addition, compared with WH, there were only 10 up-regulated genes in WL group (Figure 1 and Figure 2).
A comparative analysis of DEGs obtained by black and white Muscovy ducks showed that two co-expressed DEG genes were found between BH-vs-WH and BL-vs-WL and 7 co-expressed DEGs between BL-vs-BH and WL-vs-WH (Figure 3).

3.3. Gene Ontology and KEGG Analysis

To further assess the functional roles of DEGs in ovarian development, a GO enrichment analysis was performed. The most significant biological processes related to ovarian development in BH-vs-WH, included “circadian sleep/wake cycle process,” “negative regulation of transforming growth factor-β secretion” and “positive regulation of calcium ion transport.” The most relevant terms for cellular component were “actin cytoskeleton” and “focal adhesion.” Significant GO terms for molecular function included “extracellular matrix structural constituent” and “calcium channel regulator activity.” In BL-vs-BH, DEGs related to ovarian development were mainly enriched in biological processes of “cell surface receptor signaling pathway,” “Notch signaling pathway” and “calcium ion transport.” “Voltage-gated calcium channel complex” and “proteasome core complex” were mainly enriched in cellular component and “calcium ion binding” and “low-density lipoprotein receptor activity” were mainly enriched in molecular function. The common DEGs related to ovarian development in BL-vs-WL were primarily enriched in biological processes of “steroid biosynthetic process,” “granulosa cell development” and “egg coat formation.” DEGs were primarily enriched in cellular component, including “voltage-gated calcium channel complex” and “cyclin-dependent protein kinase holoenzyme complex.” There were some terms mainly enriched in molecular function of “transforming growth factor β receptor binding” and “calcium ion binding.” Besides, several fundamental biological processes related to ovarian development were found to be notably enriched in WL-vs-WH, such as “reproduction,” “MAPK cascade” and “mitotic cell cycle” and several cellular component were found to be notably enriched, such as “cyclin-dependent protein kinase holoenzyme complex” and “proteasome complex.” Furthermore, DEGs were primarily enriched in molecular function, including “translation regulator activity” and “lipid kinase activity” (Figure 4).
To better understand the biological functions and interaction of genes, a KEGG pathway analysis was performed for the DEGs identified. In the first 20 pathways enriched by DEGs of the comparisons, the KEGG pathways enriched significantly (P-value < 0.05) were mainly concentrated in PI3K-Akt signaling pathway, Rap1 signaling pathway, cell adhesion molecules (CAMs), thyroid hormone signaling pathway, focal adhesion, ECM-receptor interaction, ovarian steroidogenesis, progesterone-mediated oocyte maturation, oocyte meiosis and TGF-β signaling pathway. Among them, PI3K-Akt signaling pathway and ovarian steroidogenesis were involved in the reproductive process of black and white Muscovy ducks (Figure 5).
In order to verify the expression levels of DEGs observed in our RNA-seq analysis, 4 randomly selected transcripts in BL-vs-BH (TGFβ2, FGF18, POSTN and SMOC1) and 3 randomly selected transcripts WL-vs-WH (POSTN, SMOC1 and COL4A1) shown to be differentially expressed on the basis of RPKM values were validated by qPCR. The qPCR results showed the similar regulated trend in the expression of these genes. Therefore, these results indicated the high reliability and accuracy of the RNA-seq data (Figure 6).

4. Discussion

Ovarian development and follicle formation in poultry is a dynamic and complex process, including follicle growth and selection, granulosa cell proliferation and differentiation and so forth. It is very important to understand the gene expression patterns in the ovary to improve the egg production rate in poultry [25]. In this study, the gene expression mechanisms and signal pathways in the ovaries of high- and low-yield Muscovy ducks were investigated by using RNA-seq technology. A total of 555,951,128 clean reads (46,329,260 on average) were discovered and the Mapped reads and Unique mapped reads were more than 60.28% and 57.81%, respectively. Besides, more than 16,276 genes were detected in each sample.

4.1. Analysis of DEGs

According to the DEGs in black Muscovy duck with high- and low-laying, several genes that may affect Muscovy duck egg production were identified, such as TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Moreover, some studies have shown that these regulatory transcription factors interacted with each other in regulating animal reproduction.
Transforming growth factor-β (TGF-β), one of the large family of multifunctional growth factors, regulates a wide range of biological activities related to morphogenesis, development and differentiation. There is a wide range of biological effects in TGF-β isoforms, including cell differentiation, cell proliferation, cell growth, extracellular matrix formation and immune function [26]. Some studies have found that the TGFβ2 gene may be involved in reproduction, such as Hu sheep [27]. Before maturation of porcine oocytes, the expression of TGFβ2 in large follicles was greater than that in small and medium follicles [28]. There is evidence that the up regulation of TGFβ2 may induce apoptosis during the regression of postovulatory follicle in chickens [29]. Ma et al. (2002) found that TGFβ2 was preferentially expressed in the ovaries and testes of chicken embryos, which may play a key role in regulating the development of germ cells of ovaries and testes in chicken embryos [30]. Nerve growth factor receptor (NGFR) is a multifunctional cell surface receptor, which widely exists in many cell types. The gene induces apoptosis and participates in injury, nervous system development and regeneration [31]. NGFR is expressed in female reproductive organs and plays an important role in regulating the growth, development and function of reproductive organs [32]. Studies have shown that NGFR was expressed in the follicles with different sizes in porcine ovaries during estrus, which may have an important impact on the ovarian function of sows [33]. In addition, NGFR is also expressed in the granulosa cells, thecal cells, stromal cells and luteal cells of goat ovarian [34]. The transcription factor C/EBPδ (CEBPD) is a member of the basic leucine zipper family transcription factor, which is related to cell proliferation, differentiation, growth arrest and apoptosis [35]. Becker et al. (2011) found that CEBPD had a higher proliferation rate in calf endometrium under the influence of steroid hormones [36]. All these studies have identified that the DEGs have important functions in the ovary. In this work, TGFβ2, NGFR and CEBPD were significantly expressed in high-laying Muscovy ducks, indicating that they may be very important for ovarian development in Muscovy ducks.
The cytoplasmic polyadenylation element binding protein family (CPEB), a sequence-specific RNA binding protein, regulates the physiological processes of germ cell development, cell division and cell differentiation by regulating the translation process [37,38]. The CPEB2 gene belongs to the CEPB family and has been discovered in mouse germ cells [39]. Its function is to promote meiotic maturation and embryonic development of porcine oocytes [40]. Periostin (POSTN) is a 90-kDa adhesion molecule that is secreted in the extracellular matrix. POSTN binds to cells via the Integrin/FAK/AKT signaling pathway and promotes cell migration, adhesion and survival of various cell types [41,42]. POSTN is associated with follicular granulosa cell proliferation [43,44]. POSTN affects the protein in the follicle and the insufficient protein in the follicular fluid may lead to unique characteristics of follicle development in lambs [45]. SPARC-related modular calcium binding 1 (SMOC1), one of the extracellular glycoproteins, is widely expressed in various tissues with localization to basement membranes. SMOC1 may mediate intercellular signaling and cell type-specific differentiation during gonadal and reproductive tract development of fetal [46]. In this study, the differential expression of CEPB2, POSTN and SMOC1, which was in the ovarian tissues of high- and low-laying black and white Muscovy ducks, may have an important regulatory effect on the egg-laying mechanism of Muscovy ducks.
FGF18, a member of the FGF family, participates in the development of ovarian follicles [47]. FGF18 promotes the mitosis of thecal cells and apoptosis of bovine granulosa cells [48,49] and plays an important role in the regulation of steroid production in fetal bovine ovary [50]. Besides, Zhong et al. (2006) found that FGF18 gene was also presented in mouse thecal cells and oocytes [51]. Ephrin-A5 (EFNA5) is a member of the Ephrin-Eph family, which is highly expressed in granule cells [52]. EFNA5 has pivotal roles in cellular processes such as proliferation, survival, cell migration and invasion. Besides, EFNA5 also mediates cell apoptosis, proliferation and steroid production during follicular development in female mice [53]. In our study, FGF18 and EFNA5 were mainly expressed in the ovarian tissue of high-yield black Muscovy ducks, indicating that FGF18 and EFNA5 may play a key role in ovarian granulosa cell apoptosis and follicular development.
Syndecan4 (SDC4) is a cell surface heparan sulfate proteoglycan, which participates in cell adhesion [54,55] and is widely expressed in most adult tissues [56]. Studies have shown that the expression of SDC4 is related to follicular atresia in mice [57]. In addition, compared with the dominant follicles of cattle, the upregulation of SDC4 gene in ovulatory follicles may affect the ovulation promoting or luteinizing process [58]. In this study, the expression of SDC4 in the ovarian tissue of high-laying white Muscovy ducks was significantly up-regulated, which may be related to the high-egg production of Muscovy ducks.

4.2. Analysis of GO and KEGG

To better understand the possible functions involved in ovarian development in Muscovy ducks, the DEGs from each comparison group were annotated with the GO database. the results indicated that the DEGs played an important role in regulating duck ovarian development. Functional classification and pathway assignments were based on a KEGG analysis. In this study, 2 major pathways related to ovarian development were identified, namely PI3K-Akt signaling pathway and ovarian steroidogenesis. The PI3K/Akt pathway is an intracellular signaling pathway of great importance in the cell cycle process. Its plays an important role in regulating metabolism and has been widely investigated in cell proliferation, cell transformation, paracrine function and angiogenesis [59]. In the PI3K/Akt signal transduction pathway, with the help of phosphatidylinositol dependent kinase, protein kinase B (Akt) binds to cell membrane and exert biological effects [60]. Phosphatidylinositol-4,5-diphosphate 3-kinase (PI3K) is activated by many genes and it has been shown to be closely related to follicular recruitment and ovulation, selection and development of dominant follicles in mammalian ovaries [61,62,63,64].
Ovarian cell proliferation and apoptosis are important factors affecting ovarian function. Recent studies have shown that the follicular development and granulosa cell survival and proliferation can be improved by regulating the expression of PI3K/Akt signaling pathway [65]. Mishra et al. (2020) found that DEGs associated with high rates of egg production regulate PI3K-Akt signaling pathway at hypothalamic-pituitary-ovarian axis in laying hens [66]. The ovaries with reproductive activity are composed of small pre-hierarchical and maturing preovulatory follicles showing a hierarchy according to their size [67]. In addition, ovarian function is regulated by gonadotropins, including follicle stimulating hormone (FSH), luteinizing hormone (LH) and ovarian steroids. In animals, steroid hormones are synthesized from cholesterol by LH responsive theca cells and then converted into estrogens [68]. With the maturation of follicles, the concentration of steroid hormones such as estrogen and progesterone increases to exert their reproductive functions. The pathways in this study strongly indicated that PI3K-Akt signaling pathway and ovarian steroidogenesis played an important role in follicle development and productivity of Muscovy ducks.
Finally, we validated the RNA-seq results using qPCR to measure the expression of five DEGs (TGFβ2, FGF18, POSTN, SMOC1 and COL4A1), which showed that the results were reliable. In these five DEGs, the expression of POSTN and SMOC1 were high significantly different in BL-vs-BH and WL-vs-WH, indicating that POSTN and SMOC1 played an important role in regulating ovarian development of Muscovy ducks. Besides, the expression of TGFβ2 and FGF18 in BH were significantly or extremely significantly higher than that in BL and the expression of COL4A1 in WH were extremely significantly higher than that in WL, respectively. These results suggested that TGFβ2, FGF18 and COL4A1 may have relatively more important function in the regulation of egg production in Muscovy ducks.

5. Conclusions

In this study, several genes that may be important candidate genes involved in Muscovy duck ovarian development, such as TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4, were discussed. Go term annotation and KEGG pathway enrichment analysis showed that a number of significantly enriched DEGs may be related to PI3K-Akt signaling pathway and ovarian steroidogenesis in ovaries of Muscovy ducks. These data provide comprehensive gene expression information at the transcriptional level and will certainly accelerate the study of Muscovy duck ovarian development.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/12/1/57/s1, Figure S1: The Pearson correlation coefficient between samples, Table S1: Sample quality control information

Author Contributions

X.B. and J.Z. conceived the experiments. X.B. designed the experiments, analyzed the results and wrote the manuscript. J.C. and X.L. was involved in the data interpretation and writing. T.L. and S.Z. (Shanshan Zhang) contributed reagent and materials. Y.S., S.Z. (Shuya Zhang) and L.H. performed sample analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by funds from the National Waterfowl Industrial Technology System of China (CARS-42-2) and the Research on the Effect of New Breeding of the Black Muscovy Duck of Shaanxi (2019QYPY-130).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of Experimental Animal Management Committee (EAMC) of Northwest A&F University (protocol code #0726/2018 and 26 July 2018 of approval).

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated for this study can be found in the NCBI SRA (Submission: SUB7071226). Bioproject #PRJNA641492 and Biosamples #SAMN14252184-SAMN14252193.

Acknowledgments

In this section you can acknowledge any support given which is not covered by the author contribution or funding sections. This may include administrative and technical support or donations in kind (e.g., materials used for experiments).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, Z.M.; Miao, Z.W.; Chen, H.P.; Xin, Q.W.; Li, L.; Lin, R.L.; Huang, Q.L.; Zheng, N.Z. Ovarian transcriptomic analysis of shanma ducks at peak and late stages of egg production. Asian-Australas J. Anim. Sci. 2017, 30, 1215–1224. [Google Scholar] [CrossRef] [PubMed]
  2. Yang, Y.Z.; Yao, Y.; Cao, Z.F.; Gu, T.T.; Xu, Q.; Chen, G. Histological characteristics of follicles and reproductive hormone secretion during ovarian follicle development in laying geese. Poult. Sci. 2019, 98, 6063–6070. [Google Scholar] [CrossRef] [PubMed]
  3. Woodruff, T.K.; Shea, L.D. The Role of the Extracellular Matrix in Ovarian Follicle Development. Reprod. Sci. 2007, 14, 6–10. [Google Scholar] [CrossRef] [PubMed]
  4. Reddy, P.; Zheng, W.; Liu, K. Mechanisms maintaining the dormancy and survival of mammalian primordial follicles. Trends Endocrinol. Metab. 2010, 21, 96–103. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, H.; Zhang, C.-Q. Effects of Daidzein on Messenger Ribonucleic Acid Expression of Gonadotropin Receptors in Chicken Ovarian Follicles. Poult. Sci. 2008, 87, 541–545. [Google Scholar] [CrossRef]
  6. Cao, W.N.; Dong, X.Y. Research progress of mechanism of mTOR pathway involves follicular development. J. Reprod. Med. 2016, 5, 469–472. [Google Scholar]
  7. Onagbesan, O.; Bruggeman, V.; Decuypere, E. Intra-ovarian growth factors regulating ovarian function in avian species: A review. Anim. Reprod. Sci. 2009, 111, 121–140. [Google Scholar] [CrossRef]
  8. Zhang, C.; Shimada, K.; Saito, N.; Kansaku, N. Expression of messenger ribonucleic acids of luteinizing hormone and folli-cle-stimulating hormone receptors in granulosa and theca layers of chicken preovulatory follicles. Gen. Comp. Endocrinol. 1997, 105, 402–409. [Google Scholar] [CrossRef]
  9. Wu, X.; Wan, X.P.; Lan, J.; Yan, M.; Lian, S.; Rijal, M.; Huang, Z.; Li, A. Cloning, expression and polymorphism at the 5′-flanking region of the GnRH gene and their association with laying traits in Muscovy duck (Cairina moschata). Br. Poult. Sci. 2015, 56, 531–542. [Google Scholar] [CrossRef]
  10. Johnson, J.L.; Solovieva, E.V.; Bridgham, J.T. Relationship between steroidogenic acute regulatory protein expression and proges-terone production in hen granulosa cells during follicle development. Biol. Reprod. 2002, 67, 1313–1320. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, L.; Li, D.; Liu, Y.; Wang, Y.; Zhao, X.; Zhu, Q. Genetic effect of the prolactin receptor gene on egg production traits in chickens. Genet. Mol. Res. 2012, 11, 4307–4315. [Google Scholar] [CrossRef] [PubMed]
  12. 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] [PubMed] [Green Version]
  13. Chu, Q.; Zhou, B.; Xu, F.; Chen, R.; Shen, C.; Liang, T.; Li, Y.; Schinckel, A.P. Genome-wide differential mRNA expression profiles in follicles of two breeds and at two stages of estrus cycle of gilts. Sci. Rep. 2017, 7, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hatzirodos, N.; Hummitzsch, K.; Irving-Rodgers, H.F.; Breen, J.; Perry, V.E.A.; Anderson, R.A.; Rodgers, R.J. Transcript abundance of stromal and thecal cell related genes during bovine ovarian development. PLoS ONE 2019, 14, e0213575. [Google Scholar] [CrossRef] [Green Version]
  15. Talebi, R.; Ahmadi, A.; Afraz, F.; Sarry, J.; Plisson-Petit, F.; Genet, C.; Fabre, S. Transcriptome analysis of ovine granulosa cells reveals differences between small antral follicles collected during the follicular and luteal phases. Theriogenology 2018, 108, 103–117. [Google Scholar] [CrossRef]
  16. Zhang, T.; Chen, L.; Han, K.; Zhang, X.; Zhang, G.; Dai, G.; Wang, J.; Xie, K. Transcriptome analysis of ovary in relatively greater and lesser egg producing Jinghai Yellow Chicken. Anim. Reprod. Sci. 2019, 208, 106114. [Google Scholar] [CrossRef]
  17. 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]
  18. Luan, X.; Liu, D.; Cao, Z.; Luo, L.; Liu, M.; Gao, M.; Zhang, X. Transcriptome Profiling Identifies Differentially Expressed Genes in Huoyan Goose Ovaries between the Laying Period and Ceased Period. PLoS ONE 2014, 9, e113211. [Google Scholar] [CrossRef] [Green Version]
  19. Xu, X.; Zhao, X.; Lu, L.; Duan, X.; Qin, H.; Du, X.; Li, G.; Tao, Z.; Zhong, S.; Wang, G.L. Transcriptomic analysis of different stages of pigeon ovaries by RNA-sequencing. Mol. Reprod. Dev. 2016, 83, 640–648. [Google Scholar] [CrossRef]
  20. Wu, X.; Pan, X.L.; Cao, S.M.; Xu, F.Q.; Lan, L.M.; Zhang, Y.Y.; Lian, S.Y.; Yan, M.J.; Li, A. iTRAQ-based quantitative proteomic analysis provides insights into strong broodiness in Muscovy duck (Cairina moschata) combined with metabolomics analysis. J. Proteomics 2019, 204, 103401. [Google Scholar] [CrossRef]
  21. Ye, Q.; Xu, J.; Gao, X.; Ouyang, H.; Luo, W.; Nie, Q. Associations of IGF2 and DRD2 polymorphisms with laying traits in Muscovy duck. PeerJ 2017, 5, e4083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Roberts, A.; Trapnell, C.; Donaghey, J.; Rinn, J.L.; Pachter, L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011, 12, R22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2009, 26, 136–138. [Google Scholar] [CrossRef] [PubMed]
  24. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kang, L.; Cui, X.X.; Zhang, Y.J.; Yang, C.H.; Jiang, Y.L. Identification of miRNAs associated with sexual maturity in chicken ovary by Illumina small RNA deep sequencing. BMC Genom. 2013, 14, 352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Rider, C.C.; Mulloy, B. Heparin, Heparan Sulphate and the TGF-βCytokine Superfamily. Molecules 2017, 22, 713. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, Y.; Yanli, Z.; FengZhe, L.; Yang, H.; Zhu, A.; Pang, J.; Han, L.; Zhang, T.; Yao, X.; Wang, F. Genome-wide analysis of DNA Methylation profiles on sheep ovaries associated with prolificacy using whole-genome Bisulfite sequencing. BMC Genom. 2017, 18, 1–17. [Google Scholar] [CrossRef]
  28. Jackowska, M.; Kempisty, B.; Woźna, M.; Piotrowska, H.; Antosik, P.; Zawierucha, P.; Bukowska, D.; Nowicki, M.; Jaśkowski, J.M.; Brüssow, K.-P. Differential expression of GDF9, TGFB1, TGFB2 and TGFB3 in porcine oocytes isolated from follicles of different size before and after culture in vitro. Acta Vet. Hung. 2013, 61, 99–115. [Google Scholar] [CrossRef]
  29. Sundaresan, N.R.; Saxena, V.K.; Sastry, K.V.H.; Nagarajan, K.; Jain, P.; Singh, R.; Anish, D.; Ravindra, P.V.; Saxena, M.; Ahmed, K.A. Cytokines and chemokines in postovulatory follicle regression of domestic chicken (Gallus gallus domesticus). Dev. Comp. Immunol. 2008, 32, 253–264. [Google Scholar] [CrossRef]
  30. Ma, M.H.; Furuta, H.; Hiyama, Y.; Kato, Y.; Fujihara, N. Prominent expression of transforming growth factor beta2 gene in the chicken embryonic gonad as revealed by suppressive subtraction cloning. Gen. Comp. Endocrinol. 2002, 125, 311–316. [Google Scholar]
  31. Gao, X.G.; Daugherty, R.L.; Tourtellotte, W.G. Regulation of low affinity neurotrophin receptor (NGFR) by early growth response (egr) transcriptional regulators. Mol. Cell Neurosci. 2007, 36, 501–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bu, S.Y.; Yang, J.C.; Shan, A.T.; Tian, J.J.; Zhang, Q. Progress on NGF and its receptor expression in reproductive organs of femals animals. Prog. Vet. Med. 2016, 37, 113–115. [Google Scholar]
  33. Jana, B.M.; Koszykowska, M.; Czarzasta, J. Expression of nerve growth factor and its receptors, TrkA and p75, in porcine ovaries. J. Reprod. Dev. 2011, 57, 468–474. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ren, L.; Medan, M.S.; Weng, Q.; Jin, W.; Li, C.; Watanabe, G.; Taya, K. Immunolocalization of Nerve Growth Factor (NGF) and Its Receptors (TrkA and p75LNGFR) in the Reproductive Organs of Shiba Goats. J. Reprod. Dev. 2005, 51, 399–404. [Google Scholar] [CrossRef] [Green Version]
  35. Ramji, D.P.; Foka, P. CCAAT/enhancer-binding proteins: Structure, function and regulation. Biochem. J. 2002, 365, 561–575. [Google Scholar] [CrossRef] [Green Version]
  36. Becker, C.; Riedmaier, I.; Reiter, M.; Tichopad, A.; Groot, M.J.; Stolker, L.; Pfaffl, M.W.; Nielen, M.; Meyer, H. Influence of anabolic combinations of an androgen plus an estrogen on biochemical pathways in bovine uterine endometrium and ovary. J. Steroid Biochem. Mol. Biol. 2011, 125, 192–201. [Google Scholar] [CrossRef]
  37. Groisman, I.; Ivshina, M.; Marin, V.; Kennedy, N.J.; Davis, R.J.; Richter, J.D. Control of cellular senescence by CPEB. Genes Dev. 2006, 20, 2701–2712. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, X.-P.; Cooper, N.G.F. Comparative in silico Analyses of Cpeb1–4 with Functional Predictions. Bioinform. Biol. Insights 2010, 4, 61–83. [Google Scholar] [CrossRef]
  39. Kurihara, Y.; Tokuriki, M.; Myojin, R.; Hori, T.; Kuroiwa, A.; Matsuda, Y.; Sakurai, T.; Kimura, M.; Hecht, N.B.; Uesugi, S. CPEB2, a novel putative translational regulator in mouse haploid germ cells. Biol. Reprod. 2003, 69, 261–268. [Google Scholar] [CrossRef] [Green Version]
  40. Prochazkova, B.; Komrskova, P.; Kubelka, M. CPEB2 Is Necessary for Proper Porcine Meiotic Maturation and Embryonic De-velopment. Int. J. Mol. Sci. 2018, 19, 3138. [Google Scholar] [CrossRef] [Green Version]
  41. Gillan, L.; Matei, D.; Fishman, D.; Gerbin, C.S.; Karlan, B.Y.; Chang, D.D. Periostin secreted by epithelial ovarian carcinoma is a ligand for alpha(V)beta(3) and alpha(V)beta(5) integrins and promotes cell motility. Cancer Res. 2002, 62, 5358–5364. [Google Scholar] [PubMed]
  42. Matsuzawa, M.; Arai, C.; Nomura, Y.; Murata, T.; Yamakoshi, Y.; Oida, S.; Hanada, N.; Nakamura, Y. Periostin of human periodontal ligament fibroblasts promotes migration of human mesenchymal stem cell through the alphavbeta3 integrin/FAK/PI3K/Akt pathway. J. Periodontal Res. 2015, 50, 855–863. [Google Scholar] [CrossRef] [PubMed]
  43. Kulus, M.; Sujka-Kordowska, P.; Konwerska, A.; Celichowski, P.; Kranc, W.; Kulus, J.; Piotrowska-Kempisty, H.; Antosik, P.; Bukowska, D.; Iżycki, D.; et al. New Molecular Markers Involved in Regulation of Ovarian Granulosa Cell Morphogenesis, Development and Differentiation during Short-Term Primary In Vitro Culture-Transcriptomic and Histochemical Study Based on Ovaries and Individual Separated Follicles. Int. J. Mol. Sci. 2019, 20, 3966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ożegowska, K.; Brązert, M.; Ciesiółka, S.; Nawrocki, M.J.; Kranc, W.; Celichowski, P.; Jankowski, M.; Bryja, A.; Jeseta, M.; Antosik, P.; et al. In VitroGenes Involved in the Pro-cesses of Cell Proliferation, Migration, Adhesion and Tissue Development as New Potential Markers of Porcine Granulosa Cellular Processes: A Microarray Approach. DNA Cell Biol. 2019, 38, 549–560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wu, Y.; Lin, J.; Han, B.; Wang, L.; Chen, Y.; Liu, M.; Huang, J. Proteomic profiling of follicle fluids after superstimulation in one-month-old lambs. Reprod. Domest. Anim. 2017, 53, 186–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Pazin, D.E.; Albrecht, K.H. Developmental expression of Smoc1 and Smoc2 suggests potential roles in fetal gonad and repro-ductive tract differentiation. Dev. Dyn. 2009, 238, 2877–2890. [Google Scholar] [CrossRef] [Green Version]
  47. Zhang, X.; Ibrahimi, O.A.; Olsen, S.K.; Umemori, H.; Mohammadi, M.; Ornitz, D.M. Receptor Specificity of the Fibroblast Growth Factor Family. J. Biol. Chem. 2006, 281, 15694–15700. [Google Scholar] [CrossRef] [Green Version]
  48. Portela, V.M.; Dirandeh, E.; Guerrero-Netro, H.M.; Zamberlam, G.; Barreta, M.H.; Goetten, A.; Price, C. The Role of Fibroblast Growth Factor-18 in Follicular Atresia in Cattle1. Biol. Reprod. 2015, 92, 14. [Google Scholar] [CrossRef] [Green Version]
  49. Han, P.; Guerrero-Netro, H.; Estienne, A.; Price, C. Effects of fibroblast growth factors and the transcription factor, early growth response 1, on bovine theca cells. Mol. Cell. Endocrinol. 2018, 476, 96–102. [Google Scholar] [CrossRef]
  50. Da Silva, R.B.; Yang, M.Y.; Caixeta, E.S.; Castilho, A.C.; Amorim, R.L.; Price, C.A.; Fortune, J.E.; Buratini, J. Fibroblast growth factor 18 regulates steroidogenesis in fetal bovine ovarian tissue in vitro. Mol. Reprod. Dev. 2019, 86, 166–174. [Google Scholar] [CrossRef]
  51. Zhong, W.; Wang, Q.T.; Sun, T.; Wang, F.; Liu, J.; Leach, R.; Johnson, A.; Puscheck, E.E.; Rappolee, D.A. FGF ligand family mRNA ex-pression profile for mouse preimplantation embryos, early gestation human placenta and mouse trophoblast stem cells. Mol. Reprod. Dev. 2006, 73, 540–550. [Google Scholar] [CrossRef] [PubMed]
  52. Hatzirodos, N.; Hummitzsch, K.; Irving-Rodgers, H.F.; Harland, M.L.; Morris, S.; Rodgers, R.J. Transcriptome profiling of granulosa cells from bovine ovarian follicles during atresia. BMC Genom. 2014, 15, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Worku, T.; Wang, K.; Ayers, D.; Wu, D.; Rehman, Z.U.; Zhou, H.; Yang, L. Regulatory roles of ephrinA5 and its novel signaling pathway in mouse primary granulosa cell apoptosis and proliferation. Cell Cycle 2018, 17, 892–902. [Google Scholar] [CrossRef] [PubMed]
  54. Couchman, J.R. Transmembrane signaling proteoglycans. Annu. Rev. Cell Dev. Biol. 2010, 26, 89–114. [Google Scholar] [CrossRef] [PubMed]
  55. Morgan, M.R.; Humphries, M.J.; Bass, M.D. Synergistic control of cell adhesion by integrins and syndecans. Nat. Rev. Mol. Cell Biol. 2007, 8, 957–969. [Google Scholar] [CrossRef]
  56. Bellin, R.M.; Capila, I.; Lincecum, J.; Park, P.W.; Reizes, O.; Bernfield, M.R. Unlocking the secrets of syndecans: Transgenic organisms as a potential key. Glycoconj. J. 2002, 19, 295–304. [Google Scholar] [CrossRef]
  57. Ishiguro, K.; Kojima, T.; Taguchi, O.; Saito, H.; Muramatsu, T.; Kadomatsu, K. Syndecan-4 expression is associated with follicular atresia in mouse ovary. Histochem. Cell Biol. 1999, 112, 25–33. [Google Scholar] [CrossRef]
  58. Lussier, J.G.; Diouf, M.N.; Lévesque, V.; Sirois, J.; Ndiaye, K. Gene expression profiling of upregulated mRNAs in granulosa cells of bovine ovulatory follicles following stimulation with hCG. Reprod. Biol. Endocrinol. 2017, 15, 88. [Google Scholar] [CrossRef] [Green Version]
  59. Engelman, J.A.; Luo, J.; Cantley, L.C. The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat. Rev. Genet. 2006, 7, 606–619. [Google Scholar] [CrossRef]
  60. Xie, Y.; Shi, X.; Sheng, K.; Han, G.; Li, W.; Zhao, Q.; Jiang, B.; Feng, J.; Li, J.; Gu, Y. PI3K/Akt signaling transduction pathway, erythropoiesis and glycolysis in hypoxia (Review). Mol. Med. Rep. 2018, 19, 783–791. [Google Scholar] [CrossRef] [Green Version]
  61. Datta, S.R.; Brunet, A.; Greenberg, M.E. Cellular survival: A play in three Akts. Genes Dev. 1999, 13, 2905–2927. [Google Scholar] [CrossRef] [PubMed]
  62. Ryan, K.; Casey, S.M.; Canty, M.J.; Crowe, M.A.; Martin, F.; Evans, A.C. Akt and Erk signal transduction pathways are early markers of differentiation in dominant and subordinate ovarian follicles in cattle. Reproduction 2007, 133, 617–626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Zheng, W.J.; Nagaraju, G.; Liu, Z.L.; Liu, K. Functional roles of the phosphatidylinositol 3-kinases (PI3Ks) signaling in the mamma-lian ovary. Mol. Cell Endocrinol. 2012, 356, 24–30. [Google Scholar] [CrossRef] [PubMed]
  64. Dupont, J.; Reverchon, M.; Cloix, L.; Froment, P.; Rame, C. Involvement of adipokines, AMPK, PI3K and the PPAR signaling pathways in ovarian follicle development and cancer. Int. J. Dev. Biol. 2012, 56, 959–967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Wang, S.-Q.; Lin, S.; Zhu, M.; Li, C.; Chen, S.; Pu, L.; Lin, J.; Cao, L.; Zhang, Y. Acupuncture Reduces Apoptosis of Granulosa Cells in Rats with Premature Ovarian Failure Via Restoring the PI3K/Akt Signaling Pathway. Int. J. Mol. Sci. 2019, 20, 6311. [Google Scholar] [CrossRef] [Green Version]
  66. Mishra, S.K.; Chen, B.L.; Zhu, Q.; Xu, Z.X.; Ning, C.Y.; Yin, H.D.; Wang, Y.; Zhao, X.L.; Fan, X.L.; Yang, M.Y.; et al. Transcriptome analysis reveals diferentially expressed genes associated with high rates of egg production in chicken hypothalamic-pituitary-ovarian axis. Sci. Rep. 2020, 10, 5976. [Google Scholar] [CrossRef]
  67. Etches, R.J.; Petitte, J.N. Reptilian and avian follicular hierarchies: Models for the study of ovarian development. J. Exp. Zoöl. 1990, 256, 112–122. [Google Scholar] [CrossRef]
  68. Miller, W.L.; Bose, H.S. Early steps in steroidogenesis: Intracellular cholesterol trafficking. J. Lipid Res. 2011, 52, 2111–2135. [Google Scholar] [CrossRef] [Green Version]
Genes 12 00057 g001 550
Figure 1. Statistic of differentially expressed genes (DEGs).
Figure 1. Statistic of differentially expressed genes (DEGs).
Genes 12 00057 g001
Genes 12 00057 g002 550
Figure 2. The Volcano plot of DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH.
Figure 2. The Volcano plot of DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH.
Genes 12 00057 g002
Genes 12 00057 g003 550
Figure 3. A Venn diagram of DEGs showing overlap between the four comparisons.
Figure 3. A Venn diagram of DEGs showing overlap between the four comparisons.
Genes 12 00057 g003
Genes 12 00057 g004 550
Figure 4. Gene Ontology (GO) terms for DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH.
Figure 4. Gene Ontology (GO) terms for DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH.
Genes 12 00057 g004
Genes 12 00057 g005 550
Figure 5. KEGG pathways for DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH. Note: In the figure, each circle represents a KEGG pathway, the name of which being Scheme 3. 4. Verification of Differentially Expressed Genes.
Figure 5. KEGG pathways for DEGs in (A) BH-vs-WH, (B) BL-vs-BH, (C) BL-vs-WL and (D) WL-vs-WH. Note: In the figure, each circle represents a KEGG pathway, the name of which being Scheme 3. 4. Verification of Differentially Expressed Genes.
Genes 12 00057 g005
Genes 12 00057 g006 550
Figure 6. qPCR validation of DEGs in (A) BL-vs-BH and in (B) WL-vs-WH. The results were expressed as mean ± SD, * for p < 0.05, ** for p < 0.01.
Figure 6. qPCR validation of DEGs in (A) BL-vs-BH and in (B) WL-vs-WH. The results were expressed as mean ± SD, * for p < 0.05, ** for p < 0.01.
Genes 12 00057 g006
Table
Table 1. The feed composition for black Muscovy ducks.
Table 1. The feed composition for black Muscovy ducks.
IngredientContent (%)NutrientContent (%)
Corn56.00Crude protein15.700
Soybean meal23.80Calcium0.900
Corn gluten meal10.00Total phosphorus0.680
Limestone7.00Available phosphorus0.450
CaHPO41.50Salt0.370
Premix1.00Lysine0.760
NaCl0.30Methionine0.387
Lys·HCl0.30Methionine + Cystine0.654
DL-Met0.10Isoleucine0.534
Total100.00Threonine0.579
Tryptophan0.194
Crude fiber4.100
Crude fat3.400
Crude ash5.200
Avian metabolizable energy2875 Mcal·kg−1
Table
Table 2. qPCR primers of Muscovy duck.
Table 2. qPCR primers of Muscovy duck.
PrimerSequences(5′–3′)GeneBank Accession NumberProduct Length (bp)
POSTNF:AACACGCTTGAAGTTGGC
R:TCAATGAGGTGGATAACG		XM_005018343.4107
SMOC1F:CCGCTTCAGACCTGGCAATC
R:CTCAAGAGACAGGCCCAGTTTCTAC		XM_013101462.3149
TGFβ2F:ATCTACAACAGCACCAGGGAC
R:TAGCTTGGTGGGATGGCA		XM_027453142.1155
FGF18F: ATGTTTGTTGCCGAGGAG
R:TGTTTCCCGCTTGTCCTG		XM_027468582.1125
COL4A1F:GGAGAAATGGGAGTTATGGG
R:TTGGCCTTTGAGACTAACC		XM_027444310.1213
β-actinF:TATGTCGCCCTGGATTTCG
R:CTCAAGAGACAGGCCCAGTTTCTAC		XM_013101462.3162
Table
Table 3. Summary of reads and matches.
Table 3. Summary of reads and matches.
ReadsClean ReadsMapped ReadsUnique Mapped ReadsDetected Genes
BH146,357,72063.07%60.66%17,786
BH246,947,34863.87%61.37%17,822
BH346,070,47462.75%60.33%17,699
BL147,534,37467.40%64.36%16,847
BL247,067,60664.27%61.86%17,926
BL345,252,04863.94%60.66%16,276
WH146,942,58264.33%61.76%17,371
WH244,667,72661.96%59.67%17,553
WH347,963,32863.71%61.09%17,549
WL146,411,96265.58%62.91%17,553
WL246,391,89060.28%57.81%17,417
WL344,344,07064.77%61.92%17,286
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bao, X.; Song, Y.; Li, T.; Zhang, S.; Huang, L.; Zhang, S.; Cao, J.; Liu, X.; Zhang, J. Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production. Genes 2021, 12, 57. https://doi.org/10.3390/genes12010057

AMA Style

Bao X, Song Y, Li T, Zhang S, Huang L, Zhang S, Cao J, Liu X, Zhang J. Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production. Genes. 2021; 12(1):57. https://doi.org/10.3390/genes12010057

Chicago/Turabian Style

Bao, Xiuyu, Yiping Song, Tao Li, Shanshan Zhang, Lihua Huang, Shuya Zhang, Junting Cao, Xiaolin Liu, and Jianqin Zhang. 2021. "Comparative Transcriptome Profiling of Ovary Tissue between Black Muscovy Duck and White Muscovy Duck with High- and Low-Egg Production" Genes 12, no. 1: 57. https://doi.org/10.3390/genes12010057

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