- Open Access
Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq
- Qing Quan†1, 3,
- Qi Zheng†1, 2,
- Yinghui Ling1, 2Email author,
- Fugui Fang1, 2,
- Mingxing Chu4,
- Xiaorong Zhang1, 2,
- Yong Liu5 and
- Wenyong Li5
© The Author(s) 2019
- Received: 10 December 2017
- Accepted: 22 April 2019
- Published: 6 May 2019
A multitude of genes tightly regulate ovarian follicular development and hormone secretion. These complex and coordinated biological processes are altered during pregnancy. In order to further understand the regulatory role of these genes during pregnancy, it is important to screen the differentially expressed genes (DEGs) in the ovaries of pregnant and nonpregnant mammals. To detect the genes associated with the development of pregnancy in goats, DEGs from the ovaries from pregnant and nonpregnant Anhui white goats (pAWGs and nAWGs, respectively) were analyzed using RNA sequencing technology (RNA-Seq).
In this study, 13,676,394 and 13,549,560 clean reads were generated from pAWGs and nAWGs, respectively, and 1724 DEGs were identified between the two libraries. Compared with nAWGs, 1033 genes were upregulated and 691 genes were downregulated in pAWGs, including PGR, PRLR, STAR and CYP19A1, which play important roles in goat reproduction. Gene Ontology analysis showed that the DEGs were enriched for 49 functional GO terms. Kyoto Encyclopedia of Genes and Genomes analysis revealed that 397 DEGs were significantly enriched in 13 pathways, including “cell cycle”, “cytokine–cytokine receptor interaction” and “steroid biosynthesis”, suggesting that the genes may be associated with cell cycle regulation, follicular development and hormone secretion to regulate the reproduction process. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data.
The data obtained in this work enrich the genetic resources of goat and provide a further understanding of the complex molecular regulatory mechanisms occurring during the development of pregnancy and reproduction in goats. The DEGs screened in this study may play an important role in follicular development and hormone secretion and they would provide scientific basis for related research in the future.
- Differentially expressed genes
- Pregnant and nonpregnant
The goat (Capra hircus) can provide high-quality wool, meat, and other products and it is thus an important domestic and commercial animal. However, compared with other commercial animals, goat fecundity is relatively low, which has hindered the development of the goat industry. Anhui white goats (AWGs) have been used as a model to study goat fecundity in recent years  because of its precocious puberty, higher fertility, and higher leather quality than other types of goats .
Pregnancy is one of the most complex reproductive activities of mammals and is tightly regulated by the external environment, various endocrine factors, and the expression of a large number of genes [3, 4]. The ovary plays a vital role during pregnancy. There are significant differences in the activity and endocrine characteristics of the ovary during pregnancy and nonpregnancy . In the nonpregnant phase, ovulation is normal and estrogen secretion dominates to maintain normal estrous function; in contrast, with the formation of the corpus luteum, ovulation is suspended during pregnancy and progesterone secretion is gradually increased to maintain pregnancy . These differences between pregnancy and nonpregnancy may be associated with the distinct regulation of gene expression, which is implicated in various reproductive activities, including follicular development [7, 8], hormone secretion , ovulation , luteinization , and pregnancy maintenance .
Many studies have indicated that the differential expression of ovarian mRNA hormone receptor genes (follicle-stimulating hormone receptor [FSHR] and luteinizing hormone/choriogonadotropin receptor [LHR]) and growth factor genes (insulin-like growth factor 1 [IGF-I] and bone morphogenetic protein 6 [BMP6]) during follicular development significantly impact the development of preantral follicles [13–15]. Pramod et al. showed that the differential expression of bone morphogenetic protein 15 (BMP-15) and growth/differentiation factor (GDF-9) might be involved in higher prolificacy in goats . Suppression subtractive hybridization has been used to screen differentially expressed genes (DEGs) in ovarian tissues between polytocous and monotocous goats, obtaining 29 differentially expressed sequence tags . Microarray analysis has been used to analyze gene expression during early folliculogenesis in goats, and 2466 genes and three cell death-related pathways were explored . In 2014, the different expression of genes in the ovaries from uniparous and multiparous AWGs were analyzed by high-throughput sequencing. In total, 2201 DEGs were identified, and 12 genes that may be associated with AWG high prolificacy were explored . However, the information on the functions of specific genes in regulating goat pregnancy and reproduction is insufficient.
RNA sequencing technology (RNA-Seq) is mainly used for quantitative gene expression analysis of biological processes in a particular tissue or cells in specific species. With the development of high-throughput technologies, RNA-Seq has been extensively used to study specific genes in numerous species, including mice, pig, cattle, human, sheep, and goat [19–24]. It was reported that 338 genes were upregulated in the Jining grey goats and 404 were upregulated in Laiwu blank goats, and these genes might be involved in the regulation of goat fecundity and prolificacy . Ling et al. used RNA-Seq to study the DEGs between the ovaries of uniparous and multiparous goats, and identified candidate genes for goat prolificacy .
In the present study, we analyzed the gene expression levels in ovaries from pregnant and nonpregnant AWGs (pAWGs and nAWGs, respectively) using an Ion Proton™ RNA-Seq platform. The results shed further light on the role of specific genes in reproductive biological processes, including hormone secretion, ovulation, luteinization formation, and pregnancy maintenance, which will help to identify genes that could potentially be used to regulate hircine reproduction and breeding practices in the future.
In the present study, AWGs raised under identical feeding conditions at the College of Animal Science and Technology (Anhui Agricultural University, Hefei, China) were used as animal samples. Six ewes, which were 2 years old and had a similar appearance, were selected from two groups for this study. Three of them had been pregnant for 3 months, and the rest were not pregnant and were undergoing anestrus, which was determined by hormone levels. After slaughter, both whole ovaries from each ewe were immediately collected and frozen in liquid nitrogen, followed by storage at − 80 °C. Half of the ovaries were used for RNA-Seq, and the other half were used for real-time PCR. All experimental procedures involving AWGs had been approved by the ethics committee of Anhui Agricultural University, Anhui, China (permit no. AHAU20101025).
Total RNA extraction and analysis for RNA-Seq
To minimize differences among individuals, one ovary from each pair of ovaries from the three nAWGs and three pAWGs were grouped for RNA extraction using RNAiso Plus (Takara, Dalian, China) following the manufacturer’s protocol. Briefly, the correct amount of frozen tissue sample was dissolved in 1 ml RNAiso Plus. After incubation at room temperature for 5 min, the homogenate was centrifuged at 12,000g for 5 min at 4 °C. The supernatant was then transferred to a new centrifuge tube. Chloroform (1/5 volume of RNAiso Plus) was then added to the supernatant. The solution was fully emulsified and incubated for 5 min at room temperature before centrifugation at 12,000g for 15 min at 4 °C. The supernatant was then carefully transferred to another centrifuge tube. An equal amount of isopropanol was added to the supernatant. The solution was fully mixed and incubated for 10 min at room temperature before an additional centrifugation at 12,000g for 10 min at 4 °C. Then, 1 ml of 75% ethanol was slowly added to the white precipitate, which was centrifuged at 12,000g for 5 min at 4 °C. Finally, the supernatant was discarded, and the precipitate was dried and dissolved in RNase-free ddH2O. The quantity and quality of the total RNA were measured using the Agilent 2100 Bioanalyzer system (Santa Clara, CA), and the samples were stored at − 80 °C until analysis.
Library preparation and sequencing
Two groups of total RNA were used for library preparation and sequencing by pooling equal quantities (10 μg) of total RNA isolated from 3 individual pregnant or 3 individual nonpregnant goat ovaries. Briefly, the total RNA samples were treated with RNase-free DNase I (Ambion Inc., Austin, TX) to degrade any possible genomic DNA contamination. The digested products were then purified by magnetic beads. Next, the mRNA was enriched using oligo (dT) magnetic beads (for eukaryotes) and fragmented into short fragments (approximately 200 bp) using fragmentation buffer. First-strand cDNA was then synthesized with random hexamer primers using the mRNA fragments as templates. Buffer, dNTPs, RNase H, and DNA polymerase I were added to synthesize the second strand cDNA. The double-stranded cDNA was purified with the QIAQuick PCR extraction kit (Qiagen, Germany) and eluted with elution buffer for end repair and poly (A) addition. Sequencing adapters were then ligated to the 5′ and 3′ ends of the fragments. Next, the ligation products were size selected and purified by agarose gel electrophoresis. Finally, the fragments were enriched by PCR amplification (95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 45 s), purified by magnetic beads, and dissolved in the appropriate amount of elution buffer. During the quality control steps, an Agilent 2100 Bioanalyzer (Agilent, USA) was used for qualification and quantification of the sample library. Finally, the cDNA libraries were sequenced using an Ion Proton platform (Life Technologies, USA) at the Beijing Genomics Institute (BGI) in Shenzhen, China.
Sequencing data analysis
According to the requirements of this experiment and the principles of standard bioinformatics analysis, some contaminant reads, such as low-quality reads, adaptor sequences and reads with a length less than the threshold (30 bp), should be removed from the raw reads . After filtering, the valid sequences, also known as the clean reads, were retained for further analysis. To identify the gene expression patterns in each genotype from the two libraries, the clean reads were mapped to goat reference gene sequences (http://goat.kiz.ac.cn/GGD/download.htm) and goat reference genome sequences (http://goat.kiz.ac.cn/GGD/download.htm) using SOAPaligner/SOAP2 . No more than two mismatches were allowed in the alignment.
Identification of DEGs
Real-time quantitative PCR
The primers used for real-time quantitative PCR analysis
GO function enrichment and KEGG pathway analysis of the DEGs
Using nr annotation, the gene annotation of the DEGs was conducted using Blast2GO . After obtaining the GO annotation for the DEGs, Web Gene Ontology Annotation Plot (WEGO) software was used to perform GO functional classification of the DEGs and understand the distribution of the gene functions from the species at the macro level .
The Kyoto Encyclopedia of Genes and Genomes (KEGG) is the major public pathway-related database . Pathway enrichment analysis identified significantly enriched metabolic pathways or signal transduction pathways for the DEGs compared with the whole genome background. The formula for calculating the pathways is the same as for the GO analysis, with N as the number of all genes with KEGG annotations, n as the number of DEGs in N, M as the number of genes annotated to specific pathways, and m as the number of DEGs in M. Pathways with Q values ≤ 0.05 are significantly enriched for the DEGs.
Overview of sequencing data
According to the mapping analysis, approximately 98% of the reads (98.01% for nAWGs and 97.95% for pAWGs) could be mapped to the goat genome sequence. When the unique matches were regarded, slightly more than 92% (92.50% for nAWGs and 92.72% for pAWGs) of the reads corresponded to the genome, and approximately 45% of the reads (45.03% for nAWGs and 43.47% for pAWGs) could be perfectly matched to the genome.
When the reference genes were used to analyze the data, only 75.91% of the reads in nAWGs and 75.64% of the reads in pAWGs could be mapped to reference genes. Regarding the unique matches, 54.14% of the reads in nAWGs and the 54.54% in pAWGs corresponded to the reference genes. Additionally, 37.83% of the reads in nAWGs and 36.43% of the reads in pAWGs could be perfectly matched to the reference genes. The major characteristics of the two libraries are described in Additional file 1.
Identification and analysis of the DEGs in the ovaries from pregnant and nonpregnant goats
In total, 20,651 genes were detected in the study, among which 17,181 genes were co-expressed in the two libraries and 2002 and 1468 genes were specifically expressed in the ovary libraries from pAWGs and nAWGs, respectively (Additional file 2).
Analysis of the gene coverage distribution indicated that 16% of genes had 70–90% coverage, while 26% of genes in the nAWG and 27% of genes in the pAWG libraries had 90–100% coverage (Additional file 3: Fig. S1), which means the read distributions were similar between the two libraries.
Details of the differentially expressed genes (DEGs) associated with reproduction
log2 ratio (pregnant ovary/non-pregnant ovary)
gi|339460411|gb|AEJ76924.1|/0/prolactin receptor long form
gi|57163979|ref|NP_001009243.1|/5.24994e-157/steroidogenic acute regulatory protein, mitochondrial precursor
gi|440906469|gb|ELR56725.1|/0/3 beta-hydroxysteroid dehydrogenase/Delta 5–>4-isomerase, partial
gi|426223863|ref|XP_004006093.1|/1.84864e-105/PREDICTED:3-oxo-5-alpha-steroid 4-dehydrogenase 2
gi|426226402|ref|XP_004007333.1|/0/PREDICTED:lutropin-choriogonadotropic hormone receptor
gi|193876570|gb|ACF24869.1|/0/bone morphogenetic protein 1B
gi|1911705|gb|AAB50810.1|/0/11 beta-hydroxysteroid dehydrogenase type 2
gi|426247178|ref|XP_004017363.1|/0/PREDICTED: scavenger receptor class B member 1
gi|426228997|ref|XP_004008580.1|/0/PREDICTED: low-density lipoprotein receptor
gi|426222344|ref|XP_004005354.1|/0/PREDICTED: very low-density lipoprotein receptor
gi|426217043|ref|XP_004002763.1|/0/PREDICTED: 3-keto-steroid reductase
gi|296473962|tpg|DAA16077.1|/1.29128e-162/TPA: bone morphogenetic protein 6-like
gi|254892264|dbj|BAH86639.1|/2.61573e-94/vascular endothelial growth factor 188
gi|426245339|ref|XP_004016470.1|/4.26298e-173/PREDICTED:estradiol 17-beta-dehydrogenase 12
gi|338807808|gb|AEJ07667.1|/0/inhibin beta A subunit
gi|172088165|ref|NP_001116472.1|/0/cytochrome P450 19A1
gi|440900738|gb|ELR51808.1|/1.70145e-143/Inhibin beta B chain, partial
gi|440903234|gb|ELR53921.1|/1.18142e-66/Estradiol 17-beta-dehydrogenase 1
gi|193876578|gb|ACF24872.1|/0/follicle stimulating hormone receptor
gi|193876584|gb|ACF24875.1|/0/estrogen receptor beta
GO functional enrichment and KEGG pathway analysis of the DEGs
The significantly enriched pathways among the DEGs between the pregnant and nonpregnant
DEGs with pathway annotation (1413)
All genes with pathway annotation (22,885)
Complement and coagulation cascades
Staphylococcus aureus infection
Hematopoietic cell lineage
Cytokine–cytokine receptor interaction
Systemic lupus erythematosus
Vitamin digestion and absorption
We successfully identified the gene expression profiles of ovarian tissues from pregnant and nonpregnant AWGs using RNA-Seq (Quantification) and analyzed the gene expression differences between the two libraries. Furthermore, the GO enrichment and KEGG pathways for the DEGs were analyzed.
In our present study, 1724 genes were found to be significantly differentially expressed between the two libraries. Among them, a large number of DEGs were associated with prolificacy processes, follicular development, ovulation control, hormone secretion, and luteinization or pregnancy maintenance. For example, genes previously associated with reproduction, including PRLR, LHR, STAR, SPP1, BMP6, 3β-hydroxysteroid dehydrogenase estradiol (3BHSD), 17-beta-dehydrogenase 12 (HSD17B12), and BMPR1B, were found to be upregulated in pAWGs compared to nAWGs. The upregulated expression of PRLR mRNA during pregnancy may be associated with corpus luteum function for maintaining pregnancy . LHR is a specific membrane receptor on the granulosa and theca cells that binds to luteinizing hormone (LH), resulting in androgen and progesterone production. LHR may impact follicular development from the primary follicle stage onwards . 3BHSD plays a crucial role in the biosynthesis of all classes of hormonal steroids . HSD17B12 plays an important role in female fertility . STAR, which is the key rate-limiting enzyme in the production of steroid hormones, is involved in the regulation of steroid hormone biosynthesis and follicular development in the mammalian ovary . SPP1, which is also known as OPN (osteopontin), is localized to luteal cells during the luteal phase of the estrous cycle and is involved in the development and regression of the corpus luteum . Previous research has indicated that BMP6 may regulate ovary granulosa cell steroidogenesis, ovarian hormone secretion, follicular growth, and luteinization [38, 39]. A BMPRIB mutation may decrease the ability of BMP to inhibit differentiation, increase INHBA (inhibin beta A) mRNA expression in granular cells, and enhance the responsiveness of granule cells to gonadotropin . PGR, FSHR, CYP19A1, ESR2, INHA, NHBA, estradiol 17-beta-dehydrogenase 1 (HSD17B1), and inhibin beta B (INHBB) were found to be downregulated in pAWG. Recent studies showed that PGR, ESR2, FSHR, CYP19A1, HSD17B1, INHA, INHBA and INHBB played a very important role in hormone secretion, follicle growth, cell proliferation and apoptosis [41–43].
Our data showed some uniquely transcribed genes in nAWGs and pAWGs. Among these specifically expressed genes, two genes (DEFB112 and IFNE) from nAWGs and two genes (SST and SSTR1) from pAWGs may be involved in related events during pregnancy. DEFB112 is a type of B-defensin that is preferentially expressed in the male reproductive system and actively involved in sperm maturation and capacitation but is poorly documented. Further study is required to determine whether DEFB112 is involved in the reproductive process in female animals . IFNE is a novel type I IFN that has been detected in the lung, reproductive tissue, intestine and so on. IFNE may protect the female reproductive tract from viral and bacterial infection. Its expression was found to be at the lowest levels during diestrus and the highest levels at estrus . IFNE can activate the JAK-STAT signaling pathway, which plays a very important role in reproduction . SST is a 14-amino acid polypeptide with a short half-life that may inhibit the secretion of various hormones including GH and LH. SSTR1 is somatostatin receptor subtype 1. SST and its receptors are present in the ovary in various species. It was reported that SSTR2 and 5 activation modulate ovarian steroidogenesis by upregulating endogenous BMP activity in growing follicles . Previous studies have shown that SSTR1 inhibits cAMP accumulation by coupling to pertussis toxin-sensitive G proteins . However, the specific functions of SSTR1 are unclear and require further study.
Regarding GO enrichment analysis, we found that some of the correlated genes were enriched into reproduction-related GO terms. For example, CYP19A1, ESR2, PGR INHA, INHBA and STAR were enriched into GO terms for reproduction, reproduction process, molecular transducer activity, and transporter activity. The results indicated that these genes may play a specific role in goat reproduction. For the uniquely transcribed genes, the two genes (ENC1 and RGS17) from nAWGs were enriched for signaling, enzyme regulator activity, and molecular transducer activity GO terms, while the six genes (DCDC2, TBC1D15, PTHRP, SLC1A2, RDH16 and SST) from pAWGs were enriched for biological regulation, enzyme regulator activity, growth, metabolic process, and response to stimulus GO terms. Our data indicated that these genes may be related to cell proliferation and hormone secretion.
Three pathways, including “cytokine–cytokine receptor interaction”, “steroid biosynthesis” and “cell cycle”, were selected from the 13 significantly enriched pathways and were found to be intimately associated with mammalian reproduction. Most of the DEGs participating in cytokine–cytokine receptor interaction and steroid biosynthesis pathways were upregulated, including EPOR, KIT, IL1B, PRLR and BMPR1B, indicating that these pathways are activated in the ovaries of pregnant goats. It was shown that EPOR and IL1B might play a role in establishing pregnancy [49, 50]. KIT plays a very important role in the maintenance of the primordial follicle reserve and in primary to secondary follicle transition . Most DEGs were downregulated in the cell cycle pathway, such as CCNA2, CCNB1 and so on, and only one gene (cyclin D1, CCND1) was unregulated. CCNA2 regulates cell cycle progression by interacting with CDK kinases and is involved in the G2/M transition. CCNB1 is a regulatory protein involved in mitosis. Research indicates that the downregulated expression of CCNA2 may promote primordial follicle activation, while CCNB1 plays a significant role in promoting the maturation of large-follicle oocytes [52, 53]. CCND1 is a protein required for progression through G1 phase of the cell cycle and its main function is to promote cell proliferation . It can be inferred that these genes may be involved in cell cycle regulation, granulosa cell proliferation and follicular development.
Furthermore, some other signaling pathways associated with reproduction were also identified in our current analysis. These pathways included the ECM-receptor interaction, focal adhesion, and Jak-STAT, Wnt, insulin, and p53 signaling pathways, suggesting that some molecules in these pathways might also be involved in goat reproduction or maintaining the physiological activity of the ovary. However, the exact relationship between these signaling networks is not fully understood and further research is needed. Generally, the current KEGG pathway analyses have provided a solid basis for future work.
In conclusion, this study presents the transcriptomic profiles of ovarian tissue from pregnant and nonpregnant AWGs using RNA-Seq technology. 1724 genes were detected to be significantly differentially expressed in pAWGs compared with nAWGs. Among them, 1033 genes were upregulated and 691 were downregulated in pAWGs, suggesting that DEGs may play an important role in the regulation of goat ovarian function. GO annotation and KEGG pathway analyses were conducted on the DEGs between the two libraries, and 13 pathways were pointed out for significantly enriched, in which “cytokine–cytokine receptor interaction”, “cell cycle”, and “steroid biosynthesis” were found to be associated with reproduction. The DEGs and pathways identified could facilitate further to elucidate the functions of DEGs in goat follicular development and hormone secretion, which will help us to understand the relationship between genes and mammalian reproduction and provide a solid foundation for future studies.
These studies were designed by QQ, QZ and YHL. QQ and QZ carried out all the experimental analyses and prepared all figures and tables. YHL and QQ analyzed the data and drafted the manuscript. FFG, CMX, ZXR LY, and LWY contributed to revisions of the manuscript. YHL assisted in explaining the results and revised the final version of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
All authors have consented the manuscript been published.
Ethics approval and consent to participate
All procedures involving animals were approved by the Ethics Committee of the Institute of Animal Science, Chinese Academy of Agricultural Sciences (IAHCAAS201535). We comply with all the points claimed by the committee.
This research was supported by the National Natural Science Foundation of China (31772566), the Natural Science Foundation of Anhui Province (1708085MC61), the Key Research Projects of Natural Science in Anhui Colleges and Universities (KJ2017A334) and the Agricultural Science and Technology Innovation Program of China (ASTIP-IAS13).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ling YH, Ren CH, Guo XF, Xu LN, Huang YF, Luo JC, et al. Identification and characterization of microRNAs in the ovaries of multiple and uniparous goats (Capra hircus) during follicular phase. BMC Genomics. 2014;15:339.View ArticleGoogle Scholar
- Chen S, Cheng GL, Zhu DJ, Jiang XC, Zhao HL. Determination on the body properties and meat performance of Anhui white goat. Anim Husb Feed Sci. 2009;30:150–2.Google Scholar
- Gram A, Boos A, Kowalewski MP. Uterine and placental expression of canine oxytocin receptor during pregnancy and normal and induced parturition. Reprod Domest Anim. 2014;49:41–9.View ArticleGoogle Scholar
- Agaoglu OK, Agaoglu AR, Guzeloglu A, Kurar E, Kayis SA, Ozmen O, et al. Expression of hypoxia-inducible factors and vascular endothelial growth factor during pregnancy in the feline uterus. Theriogenology. 2015;84:24–33.View ArticleGoogle Scholar
- Stephens SM, Moley KH. Follicular origins of modern reproductive endocrinology. Am J Physiol Endocrinol Metab. 2009;297:E1235–6.View ArticleGoogle Scholar
- Ling YH, Guo XF, Chen T, Ding JP, Ma YH, Chu MX, et al. Characterization and analysis of differentially expressed microRNAs in hircine ovaries during the follicular and luteal phases. Anim Reprod Sci. 2016;166:47–57.View ArticleGoogle Scholar
- Orisaka M, Jiang JY, Orisaka S, Kotsuji F, Tsang BK. Growth differentiation factor 9 promotes rat preantral follicle growth by up-regulating follicular androgen biosynthesis. Endocrinology. 2009;150:2740–8.View ArticleGoogle Scholar
- Sargent KM, Lu N, Clopton DT, Pohlmeier WE, Brauer VM, Ferrara N, et al. Loss of vascular endothelial growth factor A (VEGFA) isoforms in granulosa cells using pDmrt-1-Cre or Amhr2-Cre reduces fertility by arresting follicular development and by reducing litter size in female mice. PLoS ONE. 2015;10:e0116332.View ArticleGoogle Scholar
- Yamashita H, Murayama C, Takasugi R, Miyamoto A, Shimizu T. BMP-4 suppresses progesterone production by inhibiting histone H3 acetylation of StAR in bovine granulosa cells in vitro. Mol Cell Biochem. 2011;348:183–90.View ArticleGoogle Scholar
- Galloway SM, McNatty KP, Cambridge LM, Laitinen MP, Juengel JL, McLaren RJ, et al. Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner. Nat Genet. 2000;25:279–83.View ArticleGoogle Scholar
- Poole DH, Ndiaye K, Pate JL. Expression and regulation of secreted phosphoprotein 1 in the bovine corpus luteum and effects on T cell lymphocyte chemotaxis. Reproduction. 2013;146:527–37.View ArticleGoogle Scholar
- Piccinni MP, Maggi E, Romagnani S. Role of hormone-controlled T-cell cytokines in the maintenance of pregnancy. Biochem Soc Trans. 2000;28:212–5.View ArticleGoogle Scholar
- Saraiva MV, Celestino JJ, Araújo VR, Chaves RN, Almeida AP, Lima-Verde IB, et al. Expression of follicle-stimulating hormone receptor (FSHR) in goat ovarian follicles and the impact of sequential culture medium on in vitro development of caprine preantral follicles. Zygote. 2011;19:205–14.View ArticleGoogle Scholar
- Magalhães-Padilha DM, Duarte ABG, Araújo VR, Saraiva MVA, Almeida AP, Rodrigues GQ, et al. Steady-state level of insulin-like growth factor-I (IGF-I) receptor mRNA and the effect of IGF-I on the in vitro culture of caprine preantral follicles. Theriogenology. 2012;77:206–13.View ArticleGoogle Scholar
- Frota IM, Leitão CC, Costa JJ, van den Hurk R, Saraiva MV, Figueiredo JR, et al. Levels of BMP-6 mRNA in goat ovarian follicles and in vitro effects of BMP-6 on secondary follicle development. Zygote. 2013;21:270–8.View ArticleGoogle Scholar
- Pramod RK, Sharma SK, Singhi A, Pan S, Mitra A. Differential ovarian morphometry and follicular expression of BMP15, GDF9 and BMPR1B influence the prolificacy in goat. Reprod Domest Anim. 2013;48:803–9.View ArticleGoogle Scholar
- An XP, Hou JX, Li G, Peng JY, Liu XQ, Liu HY, et al. Analysis of differentially expressed genes in ovaries of polytocous versus monotocous dairy goats using suppressive subtractive hybridization. Reprod Domest Anim. 2012;47:498–503.View ArticleGoogle Scholar
- Magalhães-Padilha DM, Geisler-Lee J, Wischral A, Gastal MO, Fonseca GR, Eloy YR, et al. Gene expression during early folliculogenesis in goats using microarray analysis. Biol Reprod. 2013;89:19.View ArticleGoogle Scholar
- Ling YH, Xiang H, Li YS, Liu Y, Zhang YH, Zhang ZJ, et al. Exploring differentially expressed genes in the ovaries of uniparous and multiparous goats using the RNA-Seq (Quantification) method. Gene. 2014;550:148–53.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticleGoogle Scholar
- Gunawan A, Sahadevan S, Neuhoff C, Große-Brinkhaus C, Gad A, Frieden L, et al. RNA deep sequencing reveals novel candidate genes and polymorphisms in boar testis and liver tissues with divergent androstenone levels. PLoS ONE. 2013;8:e63259.View ArticleGoogle Scholar
- Huang W, Khatib H. Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq. BMC Genomics. 2010;11:711.View ArticleGoogle Scholar
- Hackett NR, Butler MW, Shaykhiev R, Salit J, Omberg L, Rodriguez-Flores JL, et al. RNA-Seq quantification of the human small airway epithelium transcriptome. BMC Genomics. 2012;13:82.View ArticleGoogle Scholar
- Chen HY, Shen H, Jia B, Zhang YS, Wang XH, Zeng XC. Differential gene expression in ovaries of Qira Black sheep and Hetian sheep using RNA-Seq technique. PLoS ONE. 2015;10:e0120170.View ArticleGoogle Scholar
- Miao X, Luo Q, Qin X. Genome-wide transcriptome analysis in the ovaries of two goats identifies differentially expressed genes related to fecundity. Gene. 2016;582:69–76.View ArticleGoogle Scholar
- Kerpedjiev P, Frellsen J, Lindgreen S, Krogh A. Adaptable probabilistic mapping of short reads using position specific scoring matrices. BMC Bioinform. 2014;15:100.View ArticleGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25:1966–7.View ArticleGoogle Scholar
- Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–95.View ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-Seq data. Bioinformatics. 2010;26:136–8.View ArticleGoogle Scholar
- Ishida-Takagishi M, Enomoto A, Asai N, Ushida K, Watanabe T, Hashimoto T, et al. The Dishevelled-associating protein Daple controls thenon-canonical Wnt/Rac pathway and cell motility. Nat Commun. 2012;3:859.View ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.View ArticleGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:W293–7.View ArticleGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.View ArticleGoogle Scholar
- Zi XD, Chen DW, Wang HM. Molecular characterization, mRNA expression of prolactin receptor (PRLR) gene during pregnancy, nonpregnancy in the yak (Bos grunniens). Gen Comp Endocrinol. 2012;175:384–8.View ArticleGoogle Scholar
- Phoophitphong D, Srisuwatanasagul S, Tummaruk P. Immunohistochemical localization of luteinizing hormone receptor in the cyclic gilt ovary. Anat Histol Embryol. 2017;46:94–100.View ArticleGoogle Scholar
- Lachance Y, Luu-The V, Labrie C, Simard J, Dumont M, de Launoit Y, et al. Characterization of human 3 beta-hydroxysteroid dehydrogenase/delta 5-delta 4-isomerase gene and its expression in mammalian cell. J Biol Chem. 1990;265:20469–75.PubMedGoogle Scholar
- Kemiläinen H, Adam M, Mäki-Jouppila J, Damdimopoulou P, Damdimopoulos AE, Kere J, et al. The hydroxysteroid (17β) dehydrogenase family Gene HSD17B12 is involved in the prostaglandin synthesis pathway, the ovarian function, and regulation of fertility. Endocrinology. 2016;157:3719–30.View ArticleGoogle Scholar
- Campbell BK, Kendall NR, Baird DT. Effect of direct ovarian infusion of bone morphogenetic protein 6 (BMP6) on ovarian function in sheep. Biol Reprod. 2009;81:1016–23.View ArticleGoogle Scholar
- Khalaf M, Morera J, Bourret A, Reznik Y, Denoual C, Herlicoviez M, et al. BMP system expression in GCs from polycystic ovary syndrome women and the in vitro effects of BMP4, BMP6, and BMP7 on GC steroidogenesis. Eur J Endocrinol. 2013;168:437–44.View ArticleGoogle Scholar
- Fabre S, Pierre A, Pisselet C, Mulsant P, Lecerf F, Pohl J, et al. The Booroola mutation in sheep is associated with an alteration of the bone morphogenetic protein receptor-IB functionality. J Endocrinol. 2003;177:435–44.View ArticleGoogle Scholar
- Komatsu K, Masubuchi S. The concentration-dependent effect of progesterone on follicle growth in the mouse ovary. J Reprod Dev. 2017;63:271–7.View ArticleGoogle Scholar
- Chakravarthi VP, Khristi V, Ghosh S, Yerrathota S, Dai E, Roby KF, et al. ESR2 Is essential for gonadotropin-induced Kiss1 expression in granulosa cells. Endocrinology. 2018;159:3860–73.View ArticleGoogle Scholar
- Cadoret V, Frapsauce C, Jarrier P, Maillard V, Bonnet A, Locatelli Y, et al. Molecular evidence that follicle development is accelerated in vitro compared to in vivo. Reproduction. 2017;153:493–508.View ArticleGoogle Scholar
- Patil AA, Cai Y, Sang Y, Blecha F, Zhang G. Cross-species analysis of the mammalian beta-defensin gene family: presence of syntenic gene clusters and preferential expression in the male reproductive tract. Physiol Genomics. 2005;23:5–17.View ArticleGoogle Scholar
- Fung KY, Mangan NE, Cumming H, Horvat JC, Mayall JR, Stifter SA, et al. Interferon-ε protects the female reproductive tract from viral and bacterial infection. Science. 2013;339:1088–92.View ArticleGoogle Scholar
- Yang L, Xu L, Li Y, Li J, Bi Y, Liu W. Molecular and functional characterization of canine interferon-epsilon. J Interferon Cytokine Res. 2013;33:760–8.View ArticleGoogle Scholar
- Nakamura E, Otsuka F, Inagaki K, Tsukamoto N, Ogura-Ochi K, Miyoshi T, et al. Involvement of bone morphogenetic protein activity in somatostatin actions on ovarian steroidogenesis. J Steroid Biochem Mol Biol. 2013;134:67–74.View ArticleGoogle Scholar
- Hadcock JR, Strnad J, Eppler CM. Rat somatostatin receptor type 1 couples to G proteins and inhibition of cyclic AMP accumulation. Mol Pharmacol. 1994;45:410–6.PubMedGoogle Scholar
- Geisert R, Fazleabas A, Lucy M, Mathew D. Interaction of the conceptus and endometrium to establish pregnancy in mammals: role of interleukin 1β. Cell Tissue Res. 2012;349:825–38.View ArticleGoogle Scholar
- Ji YQ, Zhang YQ, Li MQ, Du MR, Wei DD, Li DJ. EPO improves the proliferation and inhibits apoptosis of trophoblast and decidual stromal cells through activating STAT-5 and inactivating p38 signal in human early pregnancy. Int J Clin Exp Pathol. 2011;4:765–74.PubMedPubMed CentralGoogle Scholar
- John GB, Shidler MJ, Besmer P, Castrillon DH. Kit signaling via PI3K promotes ovarian follicle maturation but is dispensable for primordial follicle activation. Dev Biol. 2009;331:292–9.View ArticleGoogle Scholar
- Tong Y, Li F, Lu Y, Cao Y, Gao J, Liu J. Rapamycin-sensitive mTORC1 signaling is involved in physiological primordial follicle activation in mouse ovary. Mol Reprod Dev. 2013;80:1018–34.View ArticleGoogle Scholar
- Liu HL, Gao Y, Zhai B, Jiang H, Ding Y, Zhang L, et al. The Effects of polyadenylation status on MPFs during in vitro porcine oocyte maturation. Cell Physiol Biochem. 2016;39:1735–45.View ArticleGoogle Scholar
- Fu M, Wang C, Li Z, Sakamaki T, Pestell RG. Minireview: cyclin D1: normal and abnormal functions. Endocrinology. 2004;145:5439–44.View ArticleGoogle Scholar