Abstract
Ghrelin, an orexigenic gut-derived peptide, is gaining increasing attention due to its multifaceted role in a number of physiological functions, including reproduction. Ghrelin exists in circulation primarily as des-acylated and acylated ghrelin. Des-acyl ghrelin, until recently considered to be an inactive form of ghrelin, is now known to have independent physiological functionality. However, the relative contribution of acyl and des-acyl ghrelin to reproductive development and function is currently unknown. Here we used ghrelin-O-acyltransferase (GOAT) knockout (KO) mice that have no measurable levels of endogenous acyl ghrelin and chronically high levels of des-acyl ghrelin, to characterize how the developmental and life-long absence of acyl ghrelin affects ovarian development and reproductive capacity. We combined the assessment of markers of reproductive maturity and the capacity to breed with measures of ovarian morphometry, as well as with ovarian RNA sequencing analysis. Our data show that while GOAT KO mice retain the capacity to breed in young adulthood, there is a diminished number of ovarian follicles (per mm3) in the juvenile and adult ovaries, due to a significant reduction in the number of small follicles, particularly the primordial follicles. We also show pronounced specific changes in the ovarian transcriptome in the juvenile GOAT KO ovary, indicative of a potential for premature ovarian development. Collectively, these findings indicate that an absence of acyl ghrelin does not prevent reproductive success but that appropriate levels of acyl and des-acyl ghrelin may be necessary for optimal ovarian maturation.
Introduction
Since its initial discovery and characterization (1), ghrelin has received increasing attention in metabolic, cardiovascular, stress, motivation, and memory research as having multiple important roles in these fields (2–10). Ghrelin has also been implicated in supporting reproductive function, acting at all levels of the hypothalamic-pituitary-gonadal (HPG) axis (11–14), including both directly and indirectly on the mammalian ovary (11, 15–17). In doing so it mediates changes in the metabolic state and in the levels of stress and reward on puberty, fertility, and fecundity (18).
Ghrelin is a 28 amino acid peptide, produced in the gastrointestinal tract. It is modified by n-octanoylation of serine at the third position (Ser-3) by the enzyme ghrelin-O-acyltransferase (GOAT), to form acylated ghrelin (19, 20). Ghrelin and GOAT are highly conserved in vertebrates (21). In circulation, ghrelin exists in at least two major bioactive forms: acylated and des-acylated ghrelin (22). It is the acylated form that acts at the growth hormone secretagogue receptor (GHSR). The receptor for des-acyl ghrelin has not yet been identified, nevertheless, des-acyl ghrelin has been shown to suppress the effects of acyl ghrelin (23, 24) and to exert independent biological effects on metabolism (25–27), cardiovascular function (9, 28), stress (10), and reproduction (29–31) and is thus an additional important target for investigation.
Both high and low levels of ghrelin appear to be detrimental for fertility, suggesting that a certain balance between circulating acyl and des-acyl ghrelin is important for reproductive potential (18). As such, acute administration of acyl ghrelin in rats impairs folliculogenesis, induces morphometric changes in the ovary, and reduces ovarian volume (32). Chronic administration of acyl or des-acyl ghrelin, or the combination of both, delays follicle maturation and reduces ovarian weight, suggesting the inhibitory effects of ghrelin on the ovary may not be solely dependent on the GHSR-mediated signaling pathway (31). In mice, both administration of a high dose of acyl ghrelin and GHSR antagonism during peri-implantation and early gestation impair fertilization, implantation, and embryo development (33). Human data show that while acyl ghrelin inhibits ovarian steroidogenesis (16), endometrial expression of the ghrelin gene and GHSR1a are decreased in infertile women (34), supporting the hypothesis that an adequate balance within the ghrelin system is required to maintain healthy reproductive function.
The focus on the physiological role of des-acyl ghrelin has only recently begun to gain attention and a large number of studies report the levels and the effects of either total or acyl ghrelin, with only a limited number of studies assessing des-acyl ghrelin (18, 35). Given that both high and low concentrations of ghrelin exert negative effects on fertility, and that some of these effects may be mediated through GHSR-independent pathways, it is imperative to further our understanding of the role des-acyl ghrelin plays in reproduction.
Genetic deletion of GOAT removes the capacity for ghrelin acylation and results in undetectable concentrations of acyl and chronically high levels of des-acyl ghrelin during development and throughout life (36, 37). These GOAT KO mice have been used to investigate the role of the GOAT-ghrelin system in metabolism and stress (36–39), and we used it here to test the hypothesis that acyl ghrelin plays a significant role in the development and function of the reproductive system. GOAT KO mice do not display developmental or overt anatomical differences from WT (40). They do, however, respond differently to conditions of altered metabolic state and stress, as evident from findings in males (10, 37). Thus, we hypothesized that the absence of GOAT (and so the absence of acyl ghrelin and high levels of des-acyl ghrelin) would alter the development of the ovary, leading to detrimental changes in the ovarian transcriptomic profile and impaired development and function of the ovary. We also tested whether these changes were maintained into adulthood and were reflected in differences in the number of ovarian follicles per mm3 of ovary or in major reproductive endpoints, including puberty onset and fecundity.
Materials and Methods
Animals
In these experiments we used female mice. GOAT KO mice on a C57/Bl6 background were obtained from Regeneron Pharmaceuticals (Tarrytown, NY) and bred (het × het; 1 male × 2 females) in the Monash Animal Services to generate WT and KO littermates, as previously described (10). Mice were group-housed under standard laboratory conditions with ad libitum access to food and water at 23°C in a 12 h light/dark cycle. All procedures described here were in accordance with the National Health and Medical Research Council Australia Code of Practice for the Care of Experimental Animals and the Monash University Animal Ethics Committee guidelines.
Ovarian Tissue Collection
To assess the effects of GOAT deletion on ovarian morphology and transcriptome, we collected ovaries from GOAT KO and WT juvenile (3 weeks old) and adult (10 weeks old) mice. Mice were deeply anesthetized by isoflurane inhalation and ovaries were excised. One ovary from each animal was snap frozen in liquid nitrogen and stored at −80°C for gene analysis, and one ovary was fixed in Bouin's solution (Sigma-Aldrich, St Louis, MO, USA) overnight, rinsed four times in 70% ethanol and stored in ethanol until processing.
Characterization of Ovarian Morphometry
Exogenous acyl ghrelin suppresses follicle maturation and reduces ovarian volume in the prepubertal ovary (31, 32). It also disrupts granulosa cell steroidogenesis (16). We therefore investigated if the deletion of GOAT and thus a change in acyl and des-acyl ghrelin concentrations would induce morphometric changes in the GOAT KO ovary. We thus assessed the number of ovarian follicles in juvenile and adult mice. As previously described (41, 42), fixed ovaries were dehydrated, embedded in paraffin and sectioned at 4 μm. For morphometric analysis, 20 sections on 10 slides, 36 μm apart, were stained with haematoxylin-eosin (H&E). Two sections per slide were assessed on the basis of an 8 μm distance between the sections, allowing a complete assessment of primordial follicle counts at this location, as per our previous publications (42, 43). Follicles were classified as: (a) primordial: an oocyte surrounded by a single layer of flattened pregranulosa cells; (b) early primary: an oocyte surrounded by a single layer of flattened pregranulosa cells with at least two cuboidal granulosa cells; (c) primary: an oocyte surrounded by cuboidal granulosa cells; (d) preantral: follicles with no antral cavity and two or more layers of cuboidal granulosa cells; (e) antral: an antral cavity visible, with at least two layers of cuboidal granulosa cells. We scanned whole ovarian sections using an Olympus VS120 slide scanner (Olympus, Tokyo, Japan). Only follicles with visible nuclei and nucleoli were counted to prevent counting the same follicle more than once. Area measurements were obtained using ImageJ (National Institutes of Health, MD, USA). Follicle counts were adjusted per total section volume, calculated as area multiplied by section thickness, according to Bernal et al. (44), Aiken et al. (45, 46), Tsoulis et al. (47), Asadi-Azarbaijani et al. (48), and Chan et al. (49), and presented as counts per mm3 (n = 4–6 animals per group). Follicles were classified as atretic if they presented with one or more of the following criteria: oocyte degeneration; granulosa cell degeneration, disorganization and retraction from the oocyte; appearance of pyknosis in more than 10% of granulosa cells (49–51).
Immunohistochemistry
We used proliferating cell nuclear antigen (PCNA) immunolabeling to visualize follicle activation and growth, as previously described (42, 52). We de-waxed paraffin-embedded sections (4 μm) in histolene and rehydrated them in ethanol washes. Antigen retrieval was carried out by microwaving sections in sodium citrate buffer for 15 min (10 mM sodium citrate, pH = 6). Slides were then cooled down to room temperature (RT) and blocked in 3% bovine serum albumin (BSA)/0.03% Triton X-100/ phosphate-buffered saline (PBS) for 1 h at RT. Sections were then incubated overnight at 4°C with mouse monoclonal anti-PCNA (1:200; #ab29, lot #GR201287, Abcam, Cambridge, UK). We then washed the slides in PBS/0.1% Triton X-100 and incubated them with Alexa Fluor 594 donkey anti-mouse IgG fluorescent-conjugated secondary antibody (1:200; A21203 Thermo Scientific, Rockford, IL, USA). Sections were then counterstained with 4′,6-diamidino-2-phenylindole (DAPI) using Fluoroshield with DAPI mounting medium (Sigma-Aldrich, St Louis, MO, USA) and viewed under an Olympus BX61 fluorescent microscope fitted with a Nikon DS-Ri2 camera. A minimum of two randomly selected slides were evaluated for each animal. We assessed fluorescence intensity in PCNA-positive follicles with a visible nucleus using ImageJ (National Institutes of Health, MD, USA). Mean fluorescence intensity was calculated using a corrected total cell fluorescence (CTCF) formula (CTCF = integrated density—(area of selected section × mean fluorescence of background readings), as previously described (53). Data were normalized to the mean fluorescence intensity of the control group and expressed in arbitrary units (AU) (43, 54, 55).
For detection of follicular apoptosis, we de-paraffinised and rehydrated paraffin-embedded sections as described above. We pre-treated the sections with 20 μg/mL proteinase K for 15 min at RT and performed analysis for terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) using ApopTag Fluorescein In situ Apoptosis Detection Kit (Merck Millipore, Burlington, MA, USA) according to the manufacturer's instructions. We then counterstained the sections using Fluoroshield with DAPI mounting medium (Sigma-Aldrich). Positive controls consisted of sections treated with DNAse I (Qiagen, Carlsbad, CA, USA) to induce non-specific DNA fragmentation, and negative control staining was performed without active Terminal deoxynucleotidyl Transferase (TdT) but including proteinase K digestion, to control for non-specific incorporation of nucleotides. Slides were viewed under a fluorescent microscope and follicles were classified as apoptotic if they contained a TUNEL-positive oocyte and/or ≥4 TUNEL-positive granulosa cells for primary, secondary and antral follicles, or >1 TUNEL-positive granulosa cells for primordial follicles (42, 56, 57).
RNA Isolation
We isolated total RNA using QIAzol reagent and RNeasy Mini Kits (QIAGEN, Valencia, CA, USA). RNA concentrations were determined using a spectrophotometer, NanoDrop 2000/2000c (Thermo Fisher Scientific, Wilmington, DE USA) and 1 μg RNA was transcribed to cDNA using an iScript cDNA synthesis kit (Bio-Rad Laboratories, Hercules, CA, USA), according to the manufacturer's instructions.
Next-Generation Sequencing
Total RNA (1 μg) samples isolated from juvenile GOAT KO and WT mouse ovaries (n = 4 per group) were submitted to the Australian Gene Research Facility (AGRF; Melbourne, VIC, Australia) for processing and bioinformatics analysis. The quality of the total RNA was assessed by a Bioanalyser (RNA Integrity Number = 10 for all samples). The samples were then sequenced on an Illumina HiSeq 2000 platform (Illumina, San Diego, CA) generating 50 bp single-reads per lane. The primary bioinformatics analysis involved de-multiplexing and quality control. The per base sequence quality indicated that the Phred quality score was above 30 for >96% of bases across all samples (58, 59). The reads were also screened for the presence of any Illumina adaptor/overrepresented sequences and cross-species contamination. The cleaned sequence reads were then aligned against the Mus musculus genome (Build version mm10). The Tophat aligner (v2.0.14) was used to map reads to the genomic sequences. The counts of reads mapping to each known gene were summarized. The transcripts were assembled with the Stringtie tool v1.0.4 (http://ccb.jhu.edu/software/stringtie/) utilizing the reads alignment with reference annotation based assembly option (RABT). The GENCODE annotation containing both coding and non-coding annotations for mouse genome version GRcm38 (Ensemble release 81) was used as a guide.
To estimate differences in gene counts between groups, differential expression analysis was undertaken using specialized R libraries from Bioconductor version 3.2 (http://www.bioconductor.org) (60). A multidimensional scaling plot revealed that two samples (one from each of the GOAT KO and WT groups) did not cluster with the rest of the samples from that group and were thus considered as outliers for further analysis. The data filter was set to 0.5 < logFC < −0.5 difference and p < 0.05. A test for over-representation of Gene Ontology (GO) terms was performed using the GOANA method (https://www.bioconductor.org/packages/devel/bioc/manuals/limma/man/limma.pdf). The clusterProfiler software package was used to analyse and visualize functional profiles (GO and Kyoto Encyclopedia of Genes and Genomes, KEGG) of gene and gene clusters (http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html). The clusterProfiler supports enrichment analysis of Reactome and KEGG with either hypergeometric test or Gene Set Enrichment Analysis (GSEA) (61).
In addition to these analyses, we also used the Ingenuity Pathway Analysis (IPA; Qiagen Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis) platform to explore further downstream and upstream effects of GOAT deletion in our dataset. The recommended set size for IPA core analysis from gene expression data is 200–3,000 molecules, and different fold change cut-offs are routinely used to allow inclusion of more differentially expressed genes for meaningful pathway analyses [e.g., (62, 63)]. We therefore performed the core analysis at 0.3 < logFC < −0.3 difference and p < 0.05, resulting in 656 analysis-ready molecules. These cut-off criteria allowed us to predict directionality of change in downstream functions and upstream regulators, accounting for a potential dilution of information as a result of whole ovary sequencing. Statistical significance was calculated using the right-tailed Fisher's Exact Test. The activation z-score, a statistical measure that assesses the match between observed and predicted upstream or downstream regulation patterns based on previous literature was also used to evaluate significance of effects on diseases and biological functions, as well as the activation and inhibition states of predicted upstream regulators [see (64) for further information on the IPA core analysis]. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE106339.
Real-Time Quantitative PCR Array
We also used a Custom RT2 PCR array (CLAM26350; Qiagen) designed to specifically examine, in juvenile and adult ovaries, the changes in the top ten over- and under-expressed genes (Table 1) that were identified by RNA sequencing in the ovaries of juvenile GOAT KO mice. Pseudogenes were excluded from the analysis. We also used this ovary to confirm the absence of GOAT (Mboat4) in the ovaries of GOAT KO mice. Total RNA, 400 ng, extracted as detailed above, was transcribed to cDNA using the RT2 First Strand Kit (Qiagen), according to the manufacturer's instructions. Samples were then diluted as per manufacturer's instructions in RT2 SYBR Green Mastermix, loaded onto 384-well PCR array plates and amplified on the QuantStudio™ 7 Flex Real-Time PCR System (Applied Biosystems, Carlsbad, CA, USA), including an initial activation step at 95°C for 10 min followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. Actb and Gapdh were used as endogenous controls. The C(t) values for these genes were averaged and used for the comparative threshold cycle (ΔΔ C(t)) calculations, where C(t) is ≤ 40. Fold changes were then calculated using the 2−ΔΔC(t) Equation (98).
Table 1
| Gene symbol | Gene name NCBI | Log2FC | Summary of function |
|---|---|---|---|
| Gm10036 | Predicted gene 10036 | 2.75 | Ribosomal protein L11 pseudogene (65). |
| Grik3 | Glutamate receptor, ionotropic, kainate 3 | 1.74 | Also known as Glur7. Plays a role in neuroactive ligand-receptor interaction, underlying glutamate-mediated excitatory synaptic transmission, and is expressed in the ovary (66, 67) |
| Spocd1 | SPOC domain containing 1 | 1.72 | Involved in negative regulation of phosphatase activity. SPOC domain (Spen paralogue and ortholog C terminal) plays a role in developmental signaling (68, 69). |
| Grem1 | Gremlin 1, DAN family BMP antagonist | 1.60 | Expressed in granulosa cells. Regulates folliculogenesis and primordial to primary follicle transition (70, 71) |
| Cyp19a1 | Cytochrome P450, family 19, subfamily a, polypeptide 1 | 1.40 | Encodes aromatase, the key enzyme in estrogen biosynthesis. Significantly increased during preovulatory follicle development (72–74). |
| Inhba | Inhibin, beta A | 1.37 | Encodes activin β A subunit that negatively regulates pituitary follicle-stimulating hormone (FSH) synthesis. Prominently expressed in granulosa cells of preantral and antral follicles. Deletion causes neonatal lethality, significant craniofacial defects and abnormal folliculogenesis (75–77). |
| Hspb7 | Heat shock protein family B (small) member 7 | 1.29 | Potent suppressor of protein aggregation, assists in the clearance of stress-induced misfolded proteins (78, 79). |
| Gm15421 | Predicted gene 15421 | 1.17 | Ribosomal protein L22 like 1 pseudogene (80). |
| Sohlh1 | Spermatogenesis and oogenesis specific basic helix-loop-helix 1 | 1.08 | Expression is confined to primordial oocytes and is required for their differentiation. In adult ovary transcript expression is decreased as the oocytes are recruited to form primary and secondary follicles (81, 82). |
| Drd4 | Dopamine receptor D4 | 1.05 | Androgen-dependent gene (83–85). Strongly expressed in mouse adult ovary, with no known function (67). |
| Leprotl1 | Leptin receptor overlapping transcript-like 1 | −0.98 | Negatively regulates growth hormone (GH) receptor expression and is overexpressed during fasting (86). Transgenic mice overexpressing leprotl1 show growth retardation (87). Expressed in granulosa cells throughout follicular development from small to ovulatory follicles and significantly increased in corpus luteum, compared to small follicles (88). |
| Dcdc2b | Doublecortin domain containing 2b | −1.02 | Encodes a member of the doublecortin family. The doublecortin domain binds tubulin and increases microtubule polymerisation (89). |
| Gm5620 | Predicted gene 5620 | −1.06 | tubulin, alpha 1B pseudogene. Non-protein coding. |
| Cited4 | Cbp/p300-interacting transactivator, with Glu/Asp-rich carboxy-terminal domain, 4 | −1.07 | Luteinising hormone (LH) target gene during ovulation. The pre-ovulatory LH surge induces Cited4 expression in cumulus and granulosa cells and this expression is required for cumulus-oocyte complex expansion and ovulation. Regulates histone acetylation in the promoters of the LH-induced target genes as a histone acetyltransferase in mouse granulosa cells undergoing luteinization after the ovulatory LH surge (90). |
| Gm13152 | Zinc finger protein 982 | −1.09 | Expressed in mouse fetal ovary in meiotic prophase, with increasing expression between E12.5 to E16.5 (91). No specific function in adult ovary has been attributed. |
| Gm13103 | Predicted gene 13103 | −1.11 | Highly expressed in adult mouse ovary. Specifically expressed in fully grown oocytes (92) |
| Rab3b | RAB3B, member RAS oncogene family | −1.15 | A marker for regulated secretion, expressed in cells with a high activity of regulated exocytosis. In the pituitary, Rab3b is essential for GnRH-regulated exocytosis in gonadotrophs (93). In the sheep ovary, Rab3b has been co-localized with oxytocin to the same luteal staining granules of the corpus luteum during the luteal phase of the estrous cycle (94) |
| Rab6b | RAB6B, member RAS oncogene family | −1.53 | Controls retrograde transport from the Golgi body to the endoplasmic reticulum and is predominantly expressed in neuronal cells (95, 96). |
| Gm6166 | Predicted gene 6166 | −1.81 | Fatty acid-binding protein, epidermal-like. Non-protein coding. |
| Myh6 | Myosin, heavy polypeptide 6, cardiac muscle, alpha | −2.05 | Involved in protein dimerization activity. Overexpressed in the ovaries of 5α-dihydrotestosterone treated rats, mimicking the hyperandrogenic state in women with polycystic ovarian syndrome (97). |
Top 10 over- and under-expressed genes in GOAT KO mouse ovary (p < 0.05) and supporting literature.
Onset of Puberty and Breeding Capacity
Changes in the availability of acyl ghrelin have been implicated in the timing of puberty onset (99, 100). We therefore examined if the deletion of GOAT and thus an absence of acyl ghrelin affected the onset of puberty in our study. Mice (n = 7–8 per group) were inspected daily, beginning at postnatal day (P)25, to determine the day of vaginal opening (a physical marker of puberty onset). When vaginal opening occurred, mice were weighed and left undisturbed until adulthood. It has been previously reported that GOAT KO mice are capable of breeding normally (37), but we wanted to test if this capacity was affected in our housing facility. We therefore also evaluated historic breeding records from Monash Animal Research Platform of 15 WT and 16 GOAT KO females ranging between 2.5 and 6 months of age at first mating, that were continuously mated for 3–5 months with WT and GOAT KO male studs, respectively. Breeding success was indicated by the mean number of pups per litter per dam, as well as by the mean number of pups in the first litter, since C57/Bl6 dams have been shown to have a higher pup mortality rate in their first than in subsequent litters (101).
Statistical Analysis
In addition to the analysis of RNA sequencing data described above, we used Student's unpaired t-tests for the assessment of ovarian morphometry, RT2 PCR array, puberty onset, and breeding data between GOAT KO and WT mice. We also used Pearson's correlation analysis to assess the relationship between RT2 PCR array and RNA sequencing fold changes. Data are presented as the mean ± SEM. Statistical significance was assumed when p ≤ 0.05.
Results
Juvenile GOAT KO Mice Have a Reduced Number of Ovarian Follicles (Per mm3)
Body weights in the juvenile phase were not affected by the absence of GOAT (data not shown). However, the juvenile GOAT KO mice had a significant reduction in their small follicle numbers per mm3 of ovarian tissue, such that GOAT KO mice had a reduction of more than 50% in the number of primordial [t(7) = 2.52, p = 0.039, Figures 1A,D] and early primary follicles [t(7) = 2.72, p = 0.029, Figures 1A,D; expressed per mm3], compared to age-matched WT controls. There were no differences in the numbers of large healthy follicles or atretic follicles (Figures 1B,C; expressed per mm3).
Figure 1
Consistent with a lack of differences in large healthy and atretic follicles, there were no significant differences in the number of TUNEL-positive follicles between WT and GOAT KO mice, in the juvenile phase (Figures 1E,F). TUNEL-positive granulosa cells were primarily expressed in large follicles (i.e., secondary and antral), consistent with previous studies showing that the current commonly used apoptotic markers are unable to detect primordial and primary follicle atresia and therefore follicle counts provide the most accurate assessment of primordial and primary follicle loss (102, 103). PCNA, a marker of follicle growth, was also highly expressed in the granulosa cells and oocytes of large follicles and some primary follicles, as previously described (104), with no differences in PCNA expression between juvenile WT and GOAT KO mice ovaries (Figures 1G,H).
Juvenile GOAT KOs Have Differences in the Ovarian Transcriptome Compared to WT Mouse Ovaries
Since reliable markers of apoptosis and proliferation in the small follicle population remain to be developed, we assessed if the reduction in the follicle numbers was associated with changes in ovarian genes and pathways related to reproduction. We thus performed RNA sequencing, characterizing the ovarian transcriptome of juvenile GOAT KO compared to WT mice. The counts of reads mapped to each known gene are summarized (Figure 2A). Differential expression analysis for estimating differences in transcripts across groups identified 14,573 genes, with a significant difference in the expression of 252 genes of at least 1-fold change (−0.5 < logFC < 0.5), p < 0.05. A summary of the RNA sequencing results and the distribution of differentially expressed genes (DEGs) are presented in Figures 2B,C. The top ten over- and under-expressed genes in the ovaries of GOAT KO mice are presented in Table 1 along with a summary of their known functions. These DEGs included several genes associated with major biological processes and functions regulating reproductive development. For instance, Grem1, Cyp19a, Inhba, and Sohlh1 play a critical role in folliculogenesis. Grem1, regulates primordial to primary follicle transition. Sohlh1, is required for oogenesis and is essential for primordial follicle activation. Inhba expression is associated with follicular growth, regulating cell proliferation, and follicle stimulating hormone (FSH) action in the ovary. Cyp19a1, encodes aromatase cytochrome P450, catalyzing a critical step in ovarian estrogen biosynthesis (71, 72, 81, 105).
Figure 2
DEGs were annotated by association with three GO term categories: biological process, molecular function and cellular component. The top 22 GO terms for each category are presented in Figure 3. These included biological processes regulating reproduction, particularly its positive regulation; immune response (e.g., defense response, regulation of immune system process, positive regulation of antigen processing and presentation, myeloid leukocyte migration, innate immune response); cell signaling and transport, and others (Figure 3A).
Figure 3
In the ontology of cellular component, GO categories of extracellular space, plasma membrane, trans-Golgi network, and major histocompatibility complex (MHC) class II protein complex, were among the most significant overrepresented GO terms in the DEGs (Figure 3B). Molecular functions, such as retinoid binding, activin receptor binding, endopeptidase/peptidase regulator and inhibitor activity and growth factor activity were among the top 22 most significant GO terms (Figure 3C). Using enrichment analysis for Reactome and KEGG pathways we identified several pathways involved in the immune response and cell signaling. Pathway enrichment results are summarized in Table 2.
Table 2
| Pathway | ID | Statistics | Annotated genes |
|---|---|---|---|
| (A) ENRICHED REACTOME PATHWAYS | |||
| Neuronal System | 5604671 | GeneRatio = 6/55; BgRatio = 237/6,598; p-value = 0.013; p.adjust = 0.121 | Kcnab3, Kcnn3, Camkk1, Gls2, Gabrb1, Kcnmb2 |
| Immune System | 5604803 | GeneRatio = 14/55; BgRatio = 1007/6,598; p-value = 0.033; p.adjust = 0.150 | Prkar2b, C3, Sh3gl2, Fbxo44, Kif5a, Ctss/H2-Ab1, Peli3, C8g, Trem2, Il6ra, Tlr7, Itgb2, Cd74 |
| Adaptive Immune System | 5604808 | GeneRatio = 8/55; BgRatio = 552/6,598; p-value = 0.085; p.adjust = 0.193 | C3, Sh3gl2, Fbxo44, Kif5a, Ctss, H2-Ab1, Itgb2, Cd74 |
| Hemostasis | 5605036 | GeneRatio = 7/55; BgRatio = 460/6,598; p-value = 0.086; p.adjust = 0.193 | Prkar2b, Igf2, Kif5a, H3f3a, Serpinc1, Itgb2, Kcnmb2 |
| Innate immune system | 5604802 | GeneRatio = 7/55; BgRatio = 487/6,598; p-value = 0.108; p.adjust = 0.194 | Prkar2b, C3, Ctss, C8g, Trem2, Tlr7, Itgb2 |
| (B) ENRICHED KEGG PATHWAYS | |||
| Drug metabolism—cytochrome P450 | mmu00982 | GeneRatio = 4/58; BgRatio = 29/3,997; p-value < 0.001; p.adjust < 0.05 | Gsta2, Gsta4, Fmo1, Adh1 |
| Retinol metabolism | mmu00830 | GeneRatio = 3/58; BgRatio = 15/3,997; p-value < 0.01; p.adjust < 0.05 | Adh1, Lrat, Cyp26b1 |
| Complement and coagulation cascades | mmu04610 | GeneRatio = 4/58; BgRatio = 35/3,997; p-value < 0.01; p.adjust < 0.05 | Serpina5, C3, C8g, Serpinc1 |
| Systemic lupus erythematosus | mmu05322 | GeneRatio = 4/58; BgRatio = 37/3,997; p-value < 0.01; p.adjust < 0.05 | C3, H2-Ab1, H3f3a, C8g |
| Staphylococcus aureus infection | mmu05150 | GeneRatio = 3/58; BgRatio = 21/3,997; p-value < 0.01; p.adjust < 0.05 | C3, H2-Ab1, Itgb2 |
| Metabolism of xenobiotics by cytochrome P450 | mmu00980 | GeneRatio = 3/58; BgRatio = 25/3,997; p-value < 0.01; p.adjust < 0.05 | Gsta2, Gsta4, Adh1 |
Pathway enrichment analysis according to Reactome and Kyoto Encyclopedia of genes and Genomes (KEGG).
GeneRatio: ratio between the number of DEGs in the pathway and the number of DEGs. BgRatio: ratio between the number of genes in the pathway and the total examined background of genes. p.adjust: p-value for hypergeometric test adjusted for Benjamini–Hochberg correction.
Ingenuity Pathway Analysis of DEGs From Juvenile GOAT KO and WT Mice
Using IPA, we found several canonical pathways were affected by GOAT deletion. The most common of these pathways were those involved in the immune response [Complement System, interleukin (IL)-6 Signaling, Dendritic Cell Maturation, Acute Phase Response Signaling; Figure 4, Supplementary Table 1], similar to the results of the pathway enrichment analysis for Reactome and KEGG pathways, as described above.
Figure 4
DEGs between the GOAT KO and WT juveniles were found to be mostly related to diseases and disorders associated with inflammation (Inflammatory Disease, Fisher's Exact Test p-value range = 1.28E-03 to 1.81E-10, 112 molecules; Inflammatory Response, p-value range 1.20E-03 to 1.37E-09, 133 molecules) and organismal injury (Organismal Injury and Abnormalities, p-value range 1.32E-03 to 1.81E-10, 456 molecules; Connective Tissue Disorders, p-value range 1.00E-03 to 1.81E-10, 111 molecules), as well as to contribute to Organismal Development (p-value range 1.30E-03 to 2.90E-07, 166 molecules). This latter biological process included 26 functions associated with Reproductive System Development and Function. The functions in this category with an absolute z-score of >1 are presented in Table 3.
Table 3
| Diseases or function annotation | p-value | Activation z-score | Genes and their direction of change | # Genes |
|---|---|---|---|---|
| Development of genital organ | 0.001 | 2.226 | Bmp7, C14orf39, Cyp11a1, Cyp19a1, Cyp26b1, Dnah9, Dnajc19, Dnd1, Fads2, Fkbp6, Fshr, Fst, H3f3a, Inhba, Inhbb, Mael, Mcmdc2, Mov10l1, Ngf, Nppc, Piwil2, Serpina5, Smc1b, Sohlh1, Sox9, Stra6, Tarbp2, Tbpl2, Tk1, Ubb, Wt1, Zmynd15 | 32 |
| Development of female reproductive tract | 0.001 | 2.113 | Cyp11a1, Cyp19A1, Fads2, Fshr, Fst, Inhba, Inhbb, Ngf, Nppc, Sohlh1, Stra6, Tbpl2, Ubb | 13 |
| Quantity of antral follicle | < 0.001 | 1.342 | Bmp7, Cyp19A1, Inhba, Ngf, Ngfr | 5 |
| Fertility | < 0.001 | 1.280 | Ctfr, Chdh, Cyp19a1, Dio3, Dppa3, Fads2, Fshr, H3f3a, H2-Ab1, Hsd17b1, Inhba, Inhbb, Khdc3l, Kmt2c, Pip5k1b, Rarg, Scara5, Serpina5, Spink13, Srd5a1, Stc1, Tk1, Wt1 | 23 |
| Quantity of gonad | < 0.001 | 1.173 | Bmp7, Cdkn2d, Cftr, Cyp19a1, Cyp26b1, Fads2, Fshr, Hsd17b1, Igf2, Il6r, Inhba, Inhbb, Itpa, Ngf, Wt1 | 16 |
| Quantity of ovarian follicle | < 0.001 | 1.122 | Bmp7, Cftr, Cyp19a1, Fads2, Fshr, Hsd17b1, Inhba, Inhbb, Ngf, Ngfr | 10 |
| Estrous cycle | 0.001 | 1 | C3, Cftr, Cyp19a1, Fads2, Fshr, Inhba, Ngf, Ngfr | 8 |
| Quantity of corpus luteum | < 0.001 | 0.896 | Cftr, Cyp19a1, Fshr, Hsd17b1, Inhba, Inhbb | 6 |
| Quantity of primary ovarian follicle | < 0.001 | −1 | Bmp7, Cyp19a1, Fshr, Ngf | 4 |
Top functions associated with Reproductive System Development and Function, identified by Ingenuity Pathway Analysis platform.
The activation z-score makes predictions by using information about the direction of gene regulation in the dataset. Positive z-score predicts an increase in the biological process or disease, while a negative z-score predicts a decrease (inhibition). −2 ≥ z-score ≥ 2 indicates a significant change. A red up arrow indicates an increase in the dataset; A green down arrow indicates a decrease in the dataset. p < 0.05 indicate a statistically significant, non-random association between a set of genes in the dataset and a related function.
To gain further insight into the biological impact of DEGs in the dataset, we performed a Regulator Effects analysis. The Regulator Effects algorithm connects potential upstream regulators with DEGs in and downstream functions that are affected in the dataset. This algorithm thus aims to provide a hypothesis that may explain how an upstream regulator affects the downstream gene expression and the impact of this activation or inhibition on biological functions and diseases (64). Upstream regulators were limited to genes, RNAs and proteins, while the diseases and functions category was limited to include the previously identified diseases and disorders associated with inflammation and organismal injury and development (Table 3). A cut-off setting of p < 0.05 and an absolute z-score of > 2 were applied. The analysis identified $10 potential regulators. FSH was identified as the main potential regulator of several genes that are likely to be involved in the development of female reproductive tract and development of genital organ, the top biological functions that are likely to be affected by GOAT deletion (See Figure 4B). This networks summary is presented in Table 4.
Table 4
| Regulators | Consistency score | Target genes in the dataset | Diseases and functions | Predicted relationship | Known regulator/disease relationship |
|---|---|---|---|---|---|
| FSH | 3 | Cyp11a1, Cyp19a1, Fshr, Fst, Inhba, Inhbb, Nppc, Tk1, Wt1 | Development of female reproductive tract; Development of genital organ | Activation | 2/2 |
| ELF4 | 1.155 | Cdkn1a, Cxcl2, Spp1 | Cell movement of phagocytes | Inhibition | 0/1 |
| STAT1 | −4.536 | Angpt2, Apoe, C3, Cdkn1a, Cxcl2, Edn1, Serpina3 | Cell movement of phagocytes | Inhibition | 1/1 |
| TNF | −6.791 | Apoe, Bmper, Cdkn11a, Col1a1, Cyp26b1, Edn1, Frzb, Fst, Mndal, Igf2, Naglu, Ngfr, Sox9, Th, Thbs2, Tnfrsf1b, Vegfc | Development of sensory organ | Activation | 1/1 |
| IL1A | −7.506 | Bcl3, Il18, Tac1 | Inflammation of the limb | Activation | 0/1 |
| CSF2 | −8.89 | C3, Cd74, Cdkn1a, Cxcl2, Edn1, Inhba, Spp1, Tnfrsf1b | Cell movement of phagocytes | Inhibition | 1/1 |
| EPO | −11.023 | Angpt2, Cdkn1a, Edn1, Il18, Spp1, Vegfc | Cell movement of phagocytes | Inhibition | 1/1 |
| LIF | −11.431 | Cdkn1a, Cxcl2, Cyp19a1, Serpina3, Spp1, Tac1 | Cell movement of phagocytes | Inhibition | 0/1 |
| IFNα | −11.839 | C3, Cdkn11a, Cxcl2, Il17ra, Il6r, Mmp28 | Cell movement of phagocytes | Inhibition | 1/1 |
| IFNγ | −14.500 | Bcl3, Il18, Itgb2, Tac1 | Inflammation of the limb | Activation | 1/1 |
Regulator Effects networks, identified by Ingenuity Pathway Analysis platform.
The networks are scaled by a consistency score, a measure of how causally consistent and densely connected a network is. A red up arrow indicates an increase in the dataset; A green down arrow indicates a decrease in the dataset.
qRT-PCR Analysis of Key DEGs in Juvenile GOAT KO Mouse Ovaries
We next performed a qRT-PCR assessment of the top DEGs in juvenile WT and GOAT KO mice, to more specifically identify individual genes that might be influenced by the absence of GOAT. While expression of certain ovarian genes can be influenced by estrous cycle stage, cyclicity was not controlled in the study to avoid additional handling and stress associated with vaginal smearing, and the clear role for ghrelin in regulation of the stress response (106). Gene expression analysis confirmed the absence of Mboat4 transcript in the juvenile GOAT KO mouse ovaries, in which no amplification was observed (data not shown). In the juvenile ovary, there was no significant correlation between the RNA sequencing and RT2 PCR array in the fold changes of the six top under-expressed genes assessed. This absence of correlation was specifically due to a different direction of change in the expression of Myh6, which was under-expressed in the RNA sequencing and over-expressed in the RT2 PCR array (Table 5). Leprotl1, a gene closely involved in the negative regulation of growth hormone (GH) receptor expression and intracellular protein trafficking (86), was under-expressed in both the RNA sequencing and the RT2 PCR array.
Table 5
| Gene symbol | Gene name NCBI | Juvenile ovary fold change (p-value) | Adult ovary fold change (p-value) |
|---|---|---|---|
| Grik3 | Glutamate receptor, ionotropic, kainate 3 | 2.7 (*0.01) | 1.30 (0.14) |
| Spocd1 | SPOC domain containing 1 | 2.6 (*0.01) | 0.74 (0.36) |
| Grem1 | Gremlin 1, DAN family BMP antagonist | 1 (0.97) | 1.67 (0.53) |
| Cyp19a1 | Cytochrome P450, family 19, subfamily a, polypeptide 1 | 0.6 (0.07) | 0.81 (0.76) |
| Inhba | Inhibin, beta A | 0.8 (0.24) | 0.9 (0.81) |
| Hspb7 | Heat shock protein family B (small) member 7 | 1.1 (0.63) | 2.85 (0.33) |
| Sohlh1 | Spermatogenesis and oogenesis specific basic helix-loop-helix 1 | 1.7 (*0.05) | 0.91 (0.66) |
| Drd4 | Dopamine receptor D4 | 0.6 (0.14) | 1.11 (0.73) |
| Leprotl1 | Leptin receptor overlapping transcript-like 1 | 0.8 (*0.02) | 0.73 (*0.04) |
| Dcdc2b | Doublecortin domain containing 2b | 1 (0.90) | 1.08 (0.68) |
| Cited4 | Cbp/p300-interacting transactivator, with Glu/Asp-rich carboxy-terminal domain, 4 | 0.8 (0.45) | 0.97 (0.93) |
| Rab3b | RAB3B, member RAS oncogene family | 1 (0.92) | 0.35 (*0.03) |
| Rab6b | RAB6B, member RAS oncogene family | 0.9 (0.63) | 2.02 (0.09) |
| Myh6 | Myosin, heavy polypeptide 6, cardiac muscle, alpha | 3.5 (* 0.01) | 0.35 (0.14) |
RT2-PCR array gene expression in juvenile and adult GOAT KO mice.
Gene expression changes in GOAT KO mice ovaries relative to WT controls.
p < 0.05.
There was a significant positive correlation between RNA sequencing and RT2 PCR array fold changes of the top over-expressed genes (r = 0.66, p = 0.04). This was driven principally by significant increases in Grik3, a glutamate ionotropic receptor encoding gene and an excitatory target in the ovary (66, 107), and Spocd1, a gene involved in negative regulation of phosphatase activity (68, 69).
Effects of GOAT on the Mature Ovary and the Capacity to Breed
Since GOAT KO mice had a pronounced reduction in the number of ovarian follicles (per mm3 of ovarian tissue) as juveniles, coupled with changes in genes and gene pathways closely involved in the regulation of reproductive development and function, we next examined if these effects were likely to be carried through into adulthood and if the breeding capacity of these mice was likely to be altered.
The reduction in the number of follicles seen in juvenile GOAT KO mice persisted into adulthood, with significantly reduced numbers of primordial [t(9) = 3.46, p = 0.007, Figure 5A; expressed per mm3] and primary follicles [t(9) = 2.39, p = 0.04, Figure 5A; expressed per mm3]. There were, again, no differences in the populations of large healthy and atretic follicles (Figures 5B,C; expressed per mm3), and no effect of GOAT deletion on apoptosis (TUNEL) or proliferation (PCNA) markers able to detect changes primarily in large follicles (Figures 5D,E). Of all the top 10 over- and under-expressed DEGs identified in GOAT KO juvenile ovaries by RNA sequencing, only two were also altered in adult GOAT KO ovaries: Rab3b, involved in exocytosis (93, 94), and Leprotl1, which was also under-expressed in the juvenile phase and is involved in the regulation of GH action (86, 87).
Figure 5
Despite these remaining subtle effects of GOAT deletion on the mature ovary, the age at puberty onset was not affected in these mice (Figure 5F). Encouragingly, there were also no differences between GOAT KO and WT mice in the number of pups born in the dam's first litter (Figure 5G) or in all litters (data not shown), suggesting pregnancy was not compromised. There were also no differences in the numbers of litters produced by continuously mated dams within a 3–5 month period (Figure 5H) or in the time between litters (data not shown). These data suggest that despite acyl ghrelin's regulatory role in reproductive development and pregnancy (31, 33, 34), it is not necessary for successful reproduction, at least not under institutional breeding facility conditions during the peak of the reproductive lifespan.
Discussion
Acylated and des-acylated ghrelin peptides regulate multiple physiological functions, including reproduction [reviewed in (18)]. While the function of des-acyl ghrelin is not yet fully elucidated, and its receptor is currently unknown, there is now substantial evidence to support its independent role in a number of physiological conditions (9, 10, 24, 108). Here, we show that genetic deletion of GOAT, an enzyme responsible for the acylation of ghrelin that thus leads to an absence of acyl and a chronic increase in the levels of des-acyl ghrelin, resulted in long-term changes in ovarian morphology, as well as changes in gene pathways associated with reproductive development and function. These changes were not reflected in the reproductive maturation timeline or breeding capacity, suggesting that while GOAT KO mice do not have an overt reproductive phenotype, some of their underlying biological functionality is notably different from that in WT mice. These findings therefore have important implications for future studies employing this global knockout model, as well as for the greater understanding of ghrelin's role in reproductive development. These data may also indicate a degree of functional redundancy within the ghrelin system to ensure reproductive success, similar to a substantial functional redundancy that exists within the gonadotropin-releasing hormone (GnRH) neuronal population, where the presence of only 12% of GnRH neurons is sufficient for pulsatile gonadotropin release and puberty onset, and 12% to 34% are sufficient for the control of estrous cyclicity in the mouse (109). It may therefore be possible that while the presence of GOAT and acyl ghrelin at certain levels is essential for optimal development of the ovary, these may not be essential in maintaining fertility. It is also important to consider that in GOAT KO mice a certain compensation may occur in the context of life-long absence of acyl ghrelin. These mice have significantly elevated levels of des-acyl ghrelin, compared to WT controls (37). While the receptor for des-acyl ghrelin remains to be discovered, it is now acknowledged to have an independent bioactivity, that alternatively counteracts or mimics the actions of acyl ghrelin. As such, des-acyl ghrelin has been shown to reduce the levels of acyl ghrelin (24) and to normalize acyl-ghrelin induced changes in insulin and glucose levels (29), while it does not affect acyl ghrelin-induced GH, prolactin or adrenocorticotropic hormone production (29, 30). It does, however, mimic the inhibitory effects of acyl ghrelin on luteinising hormone (LH) release (30). It is plausible that at least to some degree the elevated des-acyl ghrelin levels in GOAT KO mice exert compensatory effects driven by the absence of acyl ghrelin, through GHSR-independent pathways.
Our characterization of the ovarian transcriptome in juvenile mice revealed that although the number of DEGs between the two genotypes represented a relatively small subset of genes, these genes were associated with several biological processes and functions regulating reproductive development. As such, Grem1, Cyp19a1, Inhba, and Sohlh1, that were among the top ten upregulated transcripts, play a critical role in folliculogenesis. Grem1, expressed in the granulosa cells of developing follicles (110), regulates primordial to primary follicle transition, by antagonizing the members of the bone morphogenetic protein (BMP) family (111), such as the anti-Müllerian hormone (AMH) that controls the activation of primordial follicles into the growing follicle pool (71). Sohlh1, that is required for oogenesis, is also expressed in postnatal ovary where it is confined to primordial oocytes (81, 112). Sohlh1 expression, together with other oocyte-specific transcription factors, is essential for primordial follicle activation (113), with its absence leading to follicular arrest (81). Inhba expression is associated with follicular growth, regulating cell proliferation and FSH action in the ovary (76, 105, 114, 115). Together with an aromatase-encoding gene, Cyp19a1, that was also significantly upregulated in the ovaries of GOAT KOs, these transcripts modulate endocrine signaling (116).
Overall, increased expression of the above transcripts suggests the GOAT KO juvenile ovary may exhibit advanced follicle maturation, growth, and recruitment of primordial follicles into the growing pool. When we assessed the numbers of ovarian follicles (per mm3 of ovarian tissue), we saw a significant reduction in the presence of small follicles (primordial and early primary in juveniles, and primordial and primary follicles in adults). Secondary and antral were not affected by GOAT deletion at any age and this was further confirmed by the absence of follicular atresia, apoptosis, and proliferation in these follicle populations. These data suggest that by three weeks of age the number of primordial follicles (at least when expressed per mm3 of ovarian tissue) are already significantly reduced in GOAT KO mice, but these primordial follicles are not excessively recruited to grow, at least not at this age. It remains to be established whether the reduction in the number of primordial follicles in these mice is driven by a reduction in the number of embryonic germ cells; by excessive apoptosis during the mitotic-meiotic transition [embryonic days (E) 13.5–15.5]; or during the nest breakdown and primordial follicle pool formation (E17.5-P5), typically associated with a significant wave of germ cell loss and oocyte death [reviewed in (117)]; or whether this occurs later during postnatal development. It would be also of interest to examine, in future studies, if any alterations in gonadal development are also evident in male GOAT KO mice. Nevertheless, the decline in the small follicle populations in GOAT KO mice (per mm3 of ovarian tissue), particularly primordial follicles, is possibly indicative of an accelerated exhaustion of the ovarian reserve and a shortened reproductive lifespan (118). While this also remains to be explored in future studies, the reduction in the number of the primordial follicles was not associated with changes in reproductive development and function. We found no differences in the onset of puberty in GOAT KOs, as well as no changes to the reproductive capacity of these mice, as also noted in initial studies using this global knockout model (37, 40). It is important to note, however, that our assessments were conducted under standard non-stressed laboratory housing conditions and that the mice were not assessed into the period of expected reproductive senescence. We have previously shown that GOAT KO mice are more anxious than WT animals, under stressed conditions (10). GOAT is also essential for survival in a calorie-restricted environment (37). It therefore remains to be established whether the reproductive capacity of GOAT KOs is affected in a suboptimal environment, and how the increased depletion of the primordial follicles affects the timing of cessation of the reproductive lifespan.
In addition to differences in ovary specific genes and processes in juvenile GOAT KO mice, we observed significant changes in genes contributing to cell signaling and immune pathways, as identified by the pathway enrichment analysis, using the Reactome and KEGG databases, as well as by the IPA platform. For instance, two top enriched canonical pathways included EIF2 Signaling and Complement System. eIF2 (eukaryotic initiation factor-2) initiates protein translation and synthesis in ribosomes. Phosphorylation of eIF2 is among the first steps in response to cellular stress and apoptosis (119), and the EIF2 Signaling pathway is significantly enriched in human primordial oocytes during the transition from primordial to primary stage (120). The complement system integrates the interaction between the innate and adaptive immune responses, and its major role is the clearance of immune complexes and apoptotic cells (121). An upregulation of the Complement System pathway in xenobiotic-treated neonatal mouse ovaries has been suggested to underlie xenobiotic-induced ovotoxicity and primordial follicle apoptosis (122), and may be associated with the reduction in the numbers of small follicles in the GOAT KOs in our study.
Our use of the IPA platform, in addition to the Reactome and KEGG databases of pathway enrichment analysis, allowed us to make predictions for what potential upstream regulators may be modulating the DEGs in the juvenile GOAT KO ovaries, and what downstream biological functions they affect. In this analysis we focused on downstream functions associated with inflammatory diseases and disorders, as well as functions associated with organismal development, as we identified these to be most reflective of the changes in our dataset. As the result of this analysis, FSH was predicted as the main upstream regulator, affecting the expression of several genes in the dataset (Wt1, Tk1, Cyp11a1, Fst, Cyp19a1, Fshr, Nppc, Inhba, Inhbb), subsequently driving the development of genital organ and the development of female reproductive tract, the top biological functions associated with Reproductive System Development and Function, as identified by the IPA platform in our dataset. The overall predicted activation state of these downstream biological functions is once again suggestive of premature ovarian development in the GOAT KOs, which may be the cause of the significant reduction in the number of primordial follicles (per mm3) in these mice. The number of primordial follicles is a major predictor of the female reproductive lifespan. In the mammalian ovary, the vast majority of germ cells are lost before the primordial follicle formation. In the mouse ovary, the establishment of primordial follicle pool is completed by P7 (123). By P19, approximately half of these follicles are already depleted and in the post-pubertal ovary at P45, only a third of the initial population of primordial follicles are left (124). After this phenomenal loss, only a small proportion of primordial follicles will be recruited into the growing follicle pool and reach ovulation, while the remainder of the primordial follicle pool continues to gradually decline during the period of sexual maturity, until only ~4% of the primordial follicle population is left at 12 months of age (124). The period of pubertal development therefore represents an important milestone of primordial follicle depletion. This extensive depletion appears to be gonadotropin driven, and while the exact mechanisms are unknown, pubertal increases in the levels of FSH and LH that drive folliculogenesis are also likely to indirectly drive primordial follicle depletion, since GnRH antagonism during the peri-pubertal period prevents the significant primordial follicle loss that typically occurs during this time (103, 125, 126). Our sequencing data from the juvenile ovary did not reveal changes in the expression of BCL-2 modifying factor (BMF), which has recently been identified as a critical promoter of fetal oocyte and prepubertal primordial follicle loss (103, 127). However, differential expression of genes driving reproductive development and their potential regulation by FSH, as indicated by our upstream regulator analysis, warrant investigation of the quantity and quality of the follicle pool in the fetal and early postnatal GOAT KO ovary in future studies. Importantly, mutations in several genes that had significant contribution to the pathway analyses in our study, such as inhibin genes, including INHBA, as well as SOHLH1 and FSHR are associated with premature ovarian failure in humans, and these genes are among the potential candidate genes responsible for this condition (128–130).
In summary, here, for the first time, we have characterized the ovarian gene and follicle profiles, as well as the reproductive potential of female GOAT KO mice, a model that through a genetic deletion of GOAT results in an absence of circulating acyl and high levels of des-acyl ghrelin (37). Our findings indicate that while the ovarian transcriptome and follicles in these animals are affected by the global deletion of GOAT, their reproductive capacity is unchanged. Although global gene knockout may induce widespread developmental effects, these data suggest that while a presence of acyl ghrelin supports ovarian development, as is the case in WT mice, its absence is not detrimental for successful reproduction. Our data also suggest that substantial reduction in ovarian follicle numbers (per mm3 of ovarian tissue) can be sustained without overt detrimental effects on the ability to reproduce at least not during the peak of the reproductive capacity.
Statements
Author contributions
LS, ZA, and SS conceived of and designed the work. LS, JG, ZA, and SS (i.e., all authors) made substantial contributions to the acquisition, analysis and interpretation of data. LS and SS wrote the manuscript. All authors critically revised it for important intellectual content.
Funding
LS is a recipient of an RMIT Vice-Chancellor's Postdoctoral Fellowship. JG was a recipient of a funding from Conselho Nacional de Desenvolvimento Cientifico e Tecnologico-CNPq Brazil. ZA is supported by a Career Development Fellowship II from the National Health and Medical Research Council of Australia to ZA (APP1084344). SS is a recipient of a National Health and Medical Research Council Career Development Fellowship II (APP1128646) and was also supported by a Brain Foundation Research Gift.
Conflict 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2018.00815/full#supplementary-material
References
1.
KojimaMHosodaHDateYNakazatoMMatsuoHKangawaK. Ghrelin is a growth-hormone-releasing acylated peptide from stomach. Nature (1999) 402:656–60. 10.1038/45230
2.
TschopMSmileyDLHeimanML. Ghrelin induces adiposity in rodents. Nature (2000) 407:908–13. 10.1038/35038090
3.
TschopMWeyerCTataranniPADevanarayanVRavussinEHeimanML. Circulating ghrelin levels are decreased in human obesity. Diabetes (2001) 50:707–9. 10.2337/diabetes.50.4.707
4.
ChangLRenYLiuXLiWGYangJGengBet al. Protective effects of ghrelin on ischemia/reperfusion injury in the isolated rat heart. J Cardiovasc Pharmacol. (2004) 43:165–70. 10.1097/00005344-200402000-00001
5.
DianoSFarrSABenoitSCMcNayECDa SilvaIHorvathBet al. Ghrelin controls hippocampal spine synapse density and memory performance. Nat Neurosci. (2006) 9:381–8. 10.1038/nn1656
6.
SpencerSJXuLClarkeMALemusMReichenbachAGeenenBet al. Ghrelin regulates the hypothalamic-pituitary-adrenal axis and restricts anxiety after acute stress. Biol Psychiatry (2012) 72:457–65. 10.1016/j.biopsych.2012.03.010
7.
HeppnerKMPiechowskiCLMullerAOttawayNSisleySSmileyDLet al. Both acyl and des-acyl ghrelin regulate adiposity and glucose metabolism via central nervous system ghrelin receptors. Diabetes (2014) 63:122–31. 10.2337/db13-0414
8.
KuJMAndrewsZBBarsbyTReichenbachALemusMBDrummondGRet al. Ghrelin-related peptides exert protective effects in the cerebral circulation of male mice through a nonclassical ghrelin receptor(s). Endocrinology (2015) 156:280–90. 10.1210/en.2014-1415
9.
KuJMTaherMChinKYBarsbyTAustinVWongCHet al. Protective actions of des-acylated ghrelin on brain injury and blood-brain barrier disruption after stroke in mice. Clin Sci. (2016) 130:1545–58. 10.1042/CS20160077
10.
StarkRSantosVVGeenenBCabralADinanTBaylissJAet al. Des-acyl ghrelin and ghrelin O-acyltransferase regulate hypothalamic-pituitary-adrenal axis activation and anxiety in response to acute stress. Endocrinology (2016) 157:3946–57. 10.1210/en.2016-1306
11.
CaminosJETena-SempereMGaytanFSanchez-CriadoJEBarreiroMLNogueirasRet al. Expression of ghrelin in the cyclic and pregnant rat ovary. Endocrinology (2003) 144:1594–602. 10.1210/en.2002-221058
12.
GaytanFBarreiroMLCaminosJEChopinLKHeringtonACMoralesCet al. Expression of ghrelin and its functional receptor, the type 1a growth hormone secretagogue receptor, in normal human testis and testicular tumors. J Clin Endocrinol Metab. (2004) 89:400–9. 10.1210/jc.2003-031375
13.
Fernandez-FernandezRTena-SempereMNavarroVMBarreiroMLCastellanoJMAguilarEet al. Effects of ghrelin upon gonadotropin-releasing hormone and gonadotropin secretion in adult female rats: in vivo and in vitro studies. Neuroendocrinology (2005) 82:245–55. 10.1159/000092753
14.
Fernandez-FernandezRTena-SempereMRoaJCastellanoJMNavarroVMAguilarEet al. Direct stimulatory effect of ghrelin on pituitary release of LH through a nitric oxide-dependent mechanism that is modulated by estrogen. Reproduction (2007) 133:1223–32. 10.1530/REP-06-0227
15.
GaytanFBarreiroMLChopinLKHeringtonACMoralesCPinillaLet al. Immunolocalization of ghrelin and its functional receptor, the type 1a growth hormone secretagogue receptor, in the cyclic human ovary. J Clin Endocrinol Metab. (2003) 88:879–87. 10.1210/jc.2002-021196
16.
VianiIVotteroATassiFCremoniniGSartoriCBernasconiSet al. Ghrelin inhibits steroid biosynthesis by cultured granulosa-lutein cells. J Clin Endocrinol Metab. (2008) 93:1476–81. 10.1210/jc.2007-2063
17.
RakASzczepankiewiczDGregoraszczukEL. Expression of ghrelin receptor, GHSR-1a, and its functional role in the porcine ovarian follicles. Growth Horm IGF Res. (2009) 19:68–76. 10.1016/j.ghir.2008.08.006
18.
SominskyLHodgsonDMMcLaughlinEASmithRWallHMSpencerSJ. Linking stress and infertility: a novel role for ghrelin. Endocr Rev. (2017) 38:432–67. 10.1210/er.2016-1133
19.
GutierrezJASolenbergPJPerkinsDRWillencyJAKniermanMDJinZet al. Ghrelin octanoylation mediated by an orphan lipid transferase. Proc Natl Acad Sci USA. (2008) 105:6320–5. 10.1073/pnas.0800708105
20.
YangJBrownMSLiangGGrishinNVGoldsteinJL. Identification of the acyltransferase that octanoylates ghrelin, an appetite-stimulating peptide hormone. Cell (2008) 132:387–96. 10.1016/j.cell.2008.01.017
21.
MullerTDNogueirasRAndermannMLAndrewsZBAnkerSDArgenteJet al. Ghrelin. Mol Metab. (2015) 4:437–60. 10.1016/j.molmet.2015.03.005
22.
HosodaHKojimaMMatsuoHKangawaK. Ghrelin and des-acyl ghrelin: two major forms of rat ghrelin peptide in gastrointestinal tissue. Biochem Biophys Res Commun. (2000) 279:909–13. 10.1006/bbrc.2000.4039
23.
InhoffTMonnikesHNoetzelSStengelAGoebelMDinhQTet al. Desacyl ghrelin inhibits the orexigenic effect of peripherally injected ghrelin in rats. Peptides (2008) 29:2159–68. 10.1016/j.peptides.2008.09.014
24.
OzcanBNeggersSJMillerARYangHCLucaitesVAbribatTet al. Does des-acyl ghrelin improve glycemic control in obese diabetic subjects by decreasing acylated ghrelin levels?Eur J Endocrinol. (2014) 170:799–807. 10.1530/EJE-13-0347
25.
ToshinaiKYamaguchiHSunYSmithRGYamanakaASakuraiTet al. Des-acyl ghrelin induces food intake by a mechanism independent of the growth hormone secretagogue receptor. Endocrinology (2006) 147:2306–14. 10.1210/en.2005-1357
26.
DelhantyPJHuismanMBaldeon-RojasLYVan Den BergeIGrefhorstAAbribatTet al. Des-acyl ghrelin analogs prevent high-fat-diet-induced dysregulation of glucose homeostasis. FASEB J. (2013) 27:1690–700. 10.1096/fj.12-221143
27.
SominskyLZikoINguyenTXAndrewsZBSpencerSJ. Early life disruption to the ghrelin system with over-eating is resolved in adulthood in male rats. Neuropharmacology (2016) 113:21–30. 10.1016/j.neuropharm.2016.09.023
28.
BaldanziGFilighedduNCutrupiSCatapanoFBonissoniSFubiniAet al. Ghrelin and des-acyl ghrelin inhibit cell death in cardiomyocytes and endothelial cells through ERK1/2 and PI 3-kinase/AKT. J Cell Biol. (2002) 159:1029–37. 10.1083/jcb.200207165
29.
BroglioFGotteroCProdamFGaunaCMuccioliGPapottiMet al. Non-acylated ghrelin counteracts the metabolic but not the neuroendocrine response to acylated ghrelin in humans. J Clin Endocrinol Metab. (2004) 89:3062–5. 10.1210/jc.2003-031964
30.
MartiniACFernandez-FernandezRTovarSNavarroVMVigoEVazquezMJet al. Comparative analysis of the effects of ghrelin and unacylated ghrelin on luteinizing hormone secretion in male rats. Endocrinology (2006) 147:2374–82. 10.1210/en.2005-1422
31.
RoaJGarcia-GalianoDCastellanoJMGaytanFPinillaLTena-SempereM. Metabolic control of puberty onset: new players, new mechanisms. Mol Cell Endocrinol. (2010) 324:87–94. 10.1016/j.mce.2009.12.018
32.
KheradmandARoshangarLTaatiMSirotkinAV. Morphometrical and intracellular changes in rat ovaries following chronic administration of ghrelin. Tissue Cell (2009) 41:311–7. 10.1016/j.tice.2009.01.002
33.
LuqueEMTorresPJDe LoredoNVincentiLMStutzGSantillanMEet al. Role of ghrelin in fertilization, early embryo development, and implantation periods. Reproduction (2014) 148:159–67. 10.1530/REP-14-0129
34.
AghajanovaLRummanAAltmaeSWanggrenKStavreus-EversA. Diminished endometrial expression of ghrelin and ghrelin receptor contributes to infertility. Reprod Sci. (2010) 17:823–32. 10.1177/1933719110371683
35.
StarkRReichenbachALockieSHPrachtCWuQTupsAet al. Acyl ghrelin acts in the brain to control liver function and peripheral glucose homeostasis in male mice. Endocrinology (2014) 156:858–68. 10.1210/en.2014-1733
36.
KirchnerHGutierrezJASolenbergPJPflugerPTCzyzykTAWillencyJAet al. GOAT links dietary lipids with the endocrine control of energy balance. Nat Med. (2009) 15:741–5. 10.1038/nm.1997
37.
ZhaoTJLiangGLiRLXieXSleemanMWMurphyAJet al. Ghrelin O-acyltransferase (GOAT) is essential for growth hormone-mediated survival of calorie-restricted mice. Proc Natl Acad Sci USA. (2010) 107:7467–72. 10.1073/pnas.1002271107
38.
XieTYNgoSTVeldhuisJDJefferyPLChopinLKTschopMet al. Effect of deletion of ghrelin-O-acyltransferase on the pulsatile release of growth hormone in mice. J Neuroendocrinol. (2015) 27:872–86. 10.1111/jne.12327
39.
KounoTAkiyamaNItoTOkudaTNanchiINotoyaMet al. Ghrelin O-acyltransferase knockout mice show resistance to obesity when fed high-sucrose diet. J Endocrinol. (2016) 228:115–25. 10.1530/JOE-15-0330
40.
KangKZmudaESleemanMW. Physiological role of ghrelin as revealed by the ghrelin and GOAT knockout mice. Peptides (2011) 32:2236–41. 10.1016/j.peptides.2011.04.028
41.
SominskyLMeehanCLWalkerAKBobrovskayaLMcLaughlinEAHodgsonDM. Neonatal immune challenge alters reproductive development in the female rat. Horm Behav. (2012) 62:345–55. 10.1016/j.yhbeh.2012.02.005
42.
SominskyLZikoISochASmithJTSpencerSJ. Neonatal overfeeding induces early decline of the ovarian reserve: implications for the role of leptin. Mol Cell Endocrinol. (2016) 431:24–35. 10.1016/j.mce.2016.05.001
43.
FullerEASominskyLSutherlandJMRedgroveKAHarmsLMcLaughlinEAet al. Neonatal immune activation depletes the ovarian follicle reserve and alters ovarian acute inflammatory mediators in neonatal rats. Biol Reprod. (2017) 97:719–30. 10.1093/biolre/iox123
44.
BernalABVickersMHHamptonMBPoyntonRASlobodaDM. Maternal undernutrition significantly impacts ovarian follicle number and increases ovarian oxidative stress in adult rat offspring. PLoS ONE (2010) 5:e15558. 10.1371/journal.pone.0015558
45.
AikenCETarry-AdkinsJLOzanneSE. Transgenerational developmental programming of ovarian reserve. Sci. Rep. (2015) 5:16175. 10.1038/srep16175
46.
AikenCETarry-AdkinsJLPenfoldNCDeardenLOzanneSE. Decreased ovarian reserve, dysregulation of mitochondrial biogenesis, and increased lipid peroxidation in female mouse offspring exposed to an obesogenic maternal diet. FASEB J. (2016) 30:1548–56. 10.1096/fj.15-280800
47.
TsoulisMWChangPEMooreCJChanKAGohirWPetrikJJet al. Maternal high-fat diet-induced loss of fetal oocytes is associated with compromised follicle growth in adult rat offspring. Biol Reprod. (2016) 94:94. 10.1095/biolreprod.115.135004
48.
Asadi-AzarbaijaniBSantosRRJahnukainenKBraberSVan DuursenMBMToppariJet al. Developmental effects of imatinib mesylate on follicle assembly and early activation of primordial follicle pool in postnatal rat ovary. Reprod Biol. (2017) 17:25–33. 10.1016/j.repbio.2016.11.003
49.
ChanKAJazwiecPAGohirWPetrikJJSlobodaDM. Maternal nutrient restriction impairs young adult offspring ovarian signaling resulting in reproductive dysfunction and follicle loss. Biol Reprod. (2018) 98:664–82. 10.1093/biolre/ioy008
50.
PangasSARademakerAWFishmanDAWoodruffTK. Localization of the activin signal transduction components in normal human ovarian follicles: implications for autocrine and paracrine signaling in the ovary. J Clin Endocrinol Metab. (2002) 87:2644–57. 10.1210/jcem.87.6.8519
51.
TomicDMillerKPKennyHAWoodruffTKHoyerPFlawsJA. Ovarian follicle development requires Smad3. Mol Endocrinol. (2004) 18:2224–40. 10.1210/me.2003-0414
52.
SominskyLSobinoffAPJoblingMSPyeVMcLaughlinEAHodgsonDM. Immune regulation of ovarian development: programming by neonatal immune challenge. Front Neurosci. (2013) 7:100. 10.3389/fnins.2013.00100
53.
McCloyRARogersSCaldonCELorcaTCastroABurgessA. Partial inhibition of Cdk1 in G 2 phase overrides the SAC and decouples mitotic events. Cell Cycle (2014) 13:1400–12. 10.4161/cc.28401
54.
CamlinNJSobinoffAPSutherlandJMBeckettELJarnickiAGVandersRLet al. Maternal smoke exposure impairs the long-term fertility of female offspring in a murine model. Biol Reprod. (2016) 94:39. 10.1095/biolreprod.115.135848
55.
MihalasBPDe IuliisGNRedgroveKAMcLaughlinEANixonB. The lipid peroxidation product 4-hydroxynonenal contributes to oxidative stress-mediated deterioration of the ageing oocyte. Sci Rep (2017) 7:6247. 10.1038/s41598-017-06372-z
56.
AlmogBGoldRTajimaKDantesASalimKRubinsteinMet al. Leptin attenuates follicular apoptosis and accelerates the onset of puberty in immature rats. Mol Cell Endocrinol. (2001) 183:179–91. 10.1016/S0303-7207(01)00543-3
57.
Roti RotiECLeismanSKAbbottDHSalihSM. Acute doxorubicin insult in the mouse ovary is cell- and follicle-type dependent. PLoS ONE (2012) 7:e42293. 10.1371/journal.pone.0042293
58.
EwingBGreenP. Base-calling of automated sequencer traces using phred. II Error probabilities. Genome Res. (1998) 8:186–94. 10.1101/gr.8.3.186
59.
EwingBHillierLWendlMCGreenP. Base-calling of automated sequencer traces using phred. I Accuracy assessment. Genome Res. (1998) 8:175–85. 10.1101/gr.8.3.175
60.
McCarthyDJChenYSmythGK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. (2012) 40:4288–97. 10.1093/nar/gks042
61.
YuGWangLGHanYHeQY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS (2012) 16:284–7. 10.1089/omi.2011.0118
62.
DalmanMRDeeterANimishakaviGDuanZH. Fold change and p-value cutoffs significantly alter microarray interpretations. BMC Bioinformatics (2012) (Suppl. 2) 13:S11. 10.1186/1471-2105-13-S2-S11
63.
WillemsEGuerrero-BosagnaCDecuypereEJanssensSBuyseJBuysNet al. Differential expression of genes and DNA methylation associated with prenatal protein undernutrition by albumen removal in an avian model. Sci Rep. (2016) 6:20837. 10.1038/srep20837
64.
KramerAGreenJPollardJJrTugendreichS. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics (2014) 30:523–30. 10.1093/bioinformatics/btt703
65.
ZambrowiczBPAbuinARamirez-SolisRRichterLJPiggottJBeltrandelrioHet al. Wnk1 kinase deficiency lowers blood pressure in mice: a gene-trap screen to identify potential targets for therapeutic intervention. Proc Natl Acad Sci USA. (2003) 100:14109–14. 10.1073/pnas.2336103100
66.
GillSBarkerMPulidoO. Neuroexcitatory targets in the female reproductive system of the nonhuman primate (Macaca fascicularis). Toxicol Pathol. (2008) 36:478–84. 10.1177/0192623308315663
67.
YueFChengYBreschiAVierstraJWuWRybaTet al. A comparative encyclopedia of DNA elements in the mouse genome. Nature (2014) 515:355–64. 10.1038/nature13992
68.
HendrickxABeullensMCeulemansHDen AbtTVan EyndeANicolaescuEet al. Docking motif-guided mapping of the interactome of protein phosphatase-1. Chem Biol. (2009) 16:365–71. 10.1016/j.chembiol.2009.02.012
69.
FardilhaMEstevesSLKorrodi-GregorioLVintemAPDominguesSCRebeloSet al. Identification of the human testis protein phosphatase 1 interactome. Biochem Pharmacol. (2011) 82:1403–15. 10.1016/j.bcp.2011.02.018
70.
MyersMTripuraniSKMiddlebrookBEconomidesANCanalisEPangasSA. Loss of gremlin delays primordial follicle assembly but does not affect female fertility in mice. Biol Reprod. (2011) 85:1175–82. 10.1095/biolreprod.111.091728
71.
NilssonEELarsenGSkinnerMK. Roles of Gremlin 1 and Gremlin 2 in regulating ovarian primordial to primary follicle transition. Reproduction (2014) 147:865–74. 10.1530/REP-14-0005
72.
FisherCRGravesKHParlowAFSimpsonER. Characterization of mice deficient in aromatase (ArKO) because of targeted disruption of the cyp19 gene. Proc Natl Acad Sci USA. (1998) 95:6965–70. 10.1073/pnas.95.12.6965
73.
RichardsJS. Perspective: the ovarian follicle-a perspective in 2001. Endocrinology (2001) 142:2184–93. 10.1210/endo.142.6.8223
74.
FanHYLiuZShimadaMSterneckEJohnsonPFHedrickSMet al. MAPK3/1 (ERK1/2) in ovarian granulosa cells are essential for female fertility. Science (2009) 324:938–41. 10.1126/science.1171396
75.
BrownCWHouston-HawkinsDEWoodruffTKMatzukMM. Insertion of Inhbb into the Inhba locus rescues the Inhba-null phenotype and reveals new activin functions. Nat Genet. (2000) 25:453–7. 10.1038/78161
76.
PangasSAJorgezCJTranMAgnoJLiXBrownCWet al. Intraovarian activins are required for female fertility. Mol Endocrinol. (2007) 21:2458–71. 10.1210/me.2007-0146
77.
RichardsJSPangasSA. The ovary: basic biology and clinical implications. J Clin Invest. (2010) 120:963–72. 10.1172/JCI41350
78.
VosMJZijlstraMPKanonBVanWaarde-Verhagen MABruntEROosterveld-HutHMet al. HSPB7 is the most potent polyQ aggregation suppressor within the HSPB family of molecular chaperones. Hum Mol Genet. (2010) 19:4677–93. 10.1093/hmg/ddq398
79.
KeLMeijeringRAHoogstra-BerendsFMackovicovaKVosMJVan GelderICet al. HSPB1, HSPB6, HSPB7 and HSPB8 protect against RhoA GTPase-induced remodeling in tachypaced atrial myocytes. PLoS ONE (2011) 6:e20395. 10.1371/journal.pone.0020395
80.
VanburenVPiaoYDudekulaDBQianYCarterMGMartinPRet al. Assembly, verification, and initial annotation of the NIA mouse 7.4K cDNA clone set. Genome Res. (2002) 12:1999–2003. 10.1101/gr.633802
81.
PangasSAChoiYBallowDJZhaoYWestphalHMatzukMMet al. Oogenesis requires germ cell-specific transcriptional regulators Sohlh1 and Lhx8. Proc Natl Acad Sci USA. (2006) 103:8090–5. 10.1073/pnas.0601083103
82.
ToyodaSYoshimuraTMizutaJMiyazakiJ. Auto-regulation of the Sohlh1 gene by the SOHLH2/SOHLH1/SP1 complex: implications for early spermatogenesis and oogenesis. PLoS ONE (2014) 9:e101681. 10.1371/journal.pone.0101681
83.
O'shaughnessyPJAbelMCharltonHMHuBJohnstonHBakerPJ. Altered expression of genes involved in regulation of vitamin A metabolism, solute transportation, and cytoskeletal function in the androgen-insensitive tfm mouse testis. Endocrinology (2007) 148:2914–24. 10.1210/en.2006-1412
84.
SchauwaersKDe GendtKSaundersPTAtanassovaNHaelensACallewaertLet al. Loss of androgen receptor binding to selective androgen response elements causes a reproductive phenotype in a knockin mouse model. Proc Natl Acad Sci USA. (2007) 104:4961–6. 10.1073/pnas.0610814104
85.
AbelMHBabanDLeeSCharltonHMO'shaughnessyPJ. Effects of FSH on testicular mRNA transcript levels in the hypogonadal mouse. J Mol Endocrinol. (2009) 42:291–303. 10.1677/JME-08-0107
86.
WuSGrunwaldTKharitonenkovADamJJockersRDe LucaF. Increased expression of fibroblast growth factor 21 (FGF21) during chronic undernutrition causes growth hormone insensitivity in chondrocytes by inducing leptin receptor overlapping transcript (LEPROT) and leptin receptor overlapping transcript-like 1 (LEPROTL1) expression. J Biol Chem. (2013) 288:27375–83. 10.1074/jbc.M113.462218
87.
TouvierTConte-AuriolFBriandOCudejkoCPaumelleRCaronSet al. LEPROT and LEPROTL1 cooperatively decrease hepatic growth hormone action in mice. J Clin Invest. (2009) 119:3830–8. 10.1172/JCI34997
88.
NdiayeKCastonguayABenoitGSilversidesDWLussierJG. Differential regulation of Janus kinase 3 (JAK3) in bovine preovulatory follicles and identification of JAK3 interacting proteins in granulosa cells. J Ovarian Res. (2016) 9:71. 10.1186/s13048-016-0280-5
89.
CoquelleFMLevyTBergmannSWolfSGBar-ElDSapirTet al. Common and divergent roles for members of the mouse DCX superfamily. Cell Cycle (2006) 5:976–83. 10.4161/cc.5.9.2715
90.
ZhangYLXiaYYuCRichardsJSLiuJFanHY. CBP-CITED4 is required for luteinizing hormone-triggered target gene expression during ovulation. Mol Hum Reprod. (2014) 20:850–60. 10.1093/molehr/gau040
91.
SohYQJunkerJPGillMEMuellerJLVan OudenaardenAPageDC. A gene regulatory program for meiotic prophase in the fetal ovary. PLoS Genet. (2015) 11:e1005531. 10.1371/journal.pgen.1005531
92.
EvsikovAVGraberJHBrockmanJMHamplAHolbrookAESinghPet al. Cracking the egg: molecular dynamics and evolutionary aspects of the transition from the fully grown oocyte to embryo. Genes Dev. (2006) 20:2713–27. 10.1101/gad.1471006
93.
TasakaKMasumotoNMizukiJIkebuchiYOhmichiMKurachiHet al. Rab3B is essential for GnRH-induced gonadotrophin release from anterior pituitary cells. J Endocrinol. (1998) 157:267–74. 10.1677/joe.0.1570267
94.
Al-MatubsiHYBreedWJenkinGFaircloughRJ. Co-localization of Rab3B and oxytocin to electron dense granules of the sheep corpus luteum during the estrous cycle. Anat Rec. (1999) 254:214–21.
95.
OpdamFJEchardACroesHJVan Den HurkJAVan De VorstenboschRAGinselLAet al. The small GTPase Rab6B, a novel Rab6 subfamily member, is cell-type specifically expressed and localised to the Golgi apparatus. J Cell Sci. (2000) 113(Pt 15):2725–35.
96.
HutagalungAHNovickPJ. Role of Rab GTPases in membrane traffic and cell physiology. Physiol Rev. (2011) 91:119–49. 10.1152/physrev.00059.2009
97.
Salilew-WondimDWangQTesfayeDSchellanderKHoelkerMHossainMMet al. Polycystic ovarian syndrome is accompanied by repression of gene signatures associated with biosynthesis and metabolism of steroids, cholesterol and lipids. J Ovarian Res. (2015) 8:24. 10.1186/s13048-015-0151-5
98.
LivakKJSchmittgenTD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (2001) 25:402–8. 10.1006/meth.2001.1262
99.
HayashidaTNakaharaKMondalMSDateYNakazatoMKojimaMet al. Ghrelin in neonatal rats: distribution in stomach and its possible role. J Endocrinol. (2002) 173:239–45. 10.1677/joe.0.1730239
100.
Fernandez-FernandezRNavarroVMBarreiroMLVigoEMTovarSSirotkinAVet al. Effects of chronic hyperghrelinemia on puberty onset and pregnancy outcome in the rat. Endocrinology (2005) 146:3018–25. 10.1210/en.2004-1622
101.
BrownREMathiesonWBStapletonJNeumannPE. Maternal behavior in female C57BL/6J and DBA/2J inbred mice. Physiol Behav. (1999) 67:599–605. 10.1016/S0031-9384(99)00109-2
102.
TingenCMBristol-GouldSKKiesewetterSEWellingtonJTSheaLWoodruffTK. Prepubertal primordial follicle loss in mice is not due to classical apoptotic pathways. Biol Reprod. (2009) 81:16–25. 10.1095/biolreprod.108.074898
103.
LiewSHNguyenQNStrasserAFindlayJKHuttKJ. The ovarian reserve is depleted during puberty in a hormonally driven process dependent on the pro-apoptotic protein BMF. Cell Death Dis. (2017) 8:e2971. 10.1038/cddis.2017.361
104.
OktayKSchenkenRSNelsonJF. Proliferating cell nuclear antigen marks the initiation of follicular growth in the rat. Biol Reprod. (1995) 53:295–301. 10.1095/biolreprod53.2.295
105.
WoodruffTKMatherJP. Inhibin, activin and the female reproductive axis. Annu Rev Physiol. (1995) 57:219–44. 10.1146/annurev.ph.57.030195.001251
106.
SpencerSJEmmerzaalTLKoziczTAndrewsZB. Ghrelin's role in the hypothalamic-pituitary-adrenal axis stress response: implications for mood disorders. Biol Psychiatry (2015) 78:19–27. 10.1016/j.biopsych.2014.10.021
107.
GillSSPulidoOM. Glutamate receptors in peripheral tissues: current knowledge, future research, and implications for toxicology. Toxicol Pathol. (2001) 29:208–23. 10.1080/019262301317052486
108.
DelhantyPJNeggersSJVan Der LelyAJ. Should we consider des-acyl ghrelin as a separate hormone and if so, what does it do?Front Horm Res. (2014) 42:163–74. 10.1159/000358345
109.
HerbisonAEPorteousRPapeJRMoraJMHurstPR. Gonadotropin-releasing hormone neuron requirements for puberty, ovulation, and fertility. Endocrinology (2008) 149:597–604. 10.1210/en.2007-1139
110.
PangasSAJorgezCJMatzukMM. Growth differentiation factor 9 regulates expression of the bone morphogenetic protein antagonist gremlin. J Biol Chem. (2004) 279:32281–6. 10.1074/jbc.M403212200
111.
BayneRADonnachieDJKinnellHLChildsAJAndersonRA. BMP signalling in human fetal ovary somatic cells is modulated in a gene-specific fashion by GREM1 and GREM2. Mol Hum Reprod. (2016) 22:622–33. 10.1093/molehr/gaw044
112.
ShinYHRenYSuzukiHGolnoskiKJAhnHWMicoVet al. Transcription factors SOHLH1 and SOHLH2 coordinate oocyte differentiation without affecting meiosis I. J Clin Invest. (2017) 127:2106–17. 10.1172/JCI90281
113.
GeorgesAAugusteABessiereLVanetATodeschiniALVeitiaRA. FOXL2: a central transcription factor of the ovary. J Mol Endocrinol. (2014) 52:R17–33. 10.1530/JME-13-0159
114.
WoodruffTKKrummenLMcCrayGMatherJP. In situ ligand binding of recombinant human [125I] activin-A and recombinant human [125I]inhibin-A to the adult rat ovary. Endocrinology (1993) 133:2998–3006. 10.1210/endo.133.6.8243328
115.
WoodruffTKBeseckeLMGroomeNDraperLBSchwartzNBWeissJ. Inhibin A and inhibin B are inversely correlated to follicle-stimulating hormone, yet are discordant during the follicular phase of the rat estrous cycle, and inhibin A is expressed in a sexually dimorphic manner. Endocrinology (1996) 137:5463–67. 10.1210/endo.137.12.8940372
116.
ParrishEMSiletzAXuMWoodruffTKSheaLD. Gene expression in mouse ovarian follicle development in vivo versus an ex vivo alginate culture system. Reproduction (2011) 142:309–18. 10.1530/REP-10-0481
117.
FindlayJKHuttKJHickeyMAndersonRA. How is the number of primordial follicles in the ovarian reserve established?Biol Reprod. (2015) 93:111. 10.1095/biolreprod.115.133652
118.
NelsonSMTelferEEAndersonRA. The ageing ovary and uterus: new biological insights. Hum Reprod Update (2013) 19:67–83. 10.1093/humupd/dms043
119.
RonD. Translational control in the endoplasmic reticulum stress response. J Clin Invest. (2002) 110:1383–8. 10.1172/JCI0216784
120.
ErnstEHGrondahlMLGrundSHardyKHeuckASundeLet al. Dormancy and activation of human oocytes from primordial and primary follicles: molecular clues to oocyte regulation. Hum Reprod. (2017) 32:1684–700. 10.1093/humrep/dex238
121.
SarmaJVWardPA. The complement system. Cell Tissue Res. (2011) 343:227–35. 10.1007/s00441-010-1034-0
122.
SobinoffAPMahonyMNixonBRomanSDMcLaughlinEA. Understanding the Villain: DMBA-induced preantral ovotoxicity involves selective follicular destruction and primordial follicle activation through PI3K/Akt and mTOR signaling. Toxicol Sci. (2011) 123:563–75. 10.1093/toxsci/kfr195
123.
PeplingMESpradlingAC. Mouse ovarian germ cell cysts undergo programmed breakdown to form primordial follicles. Dev Biol. (2001) 234:339–51. 10.1006/dbio.2001.0269
124.
Bristol-GouldSKKreegerPKSelkirkCGKilenSMMayoKESheaLDet al. Fate of the initial follicle pool: empirical and mathematical evidence supporting its sufficiency for adult fertility. Dev Biol. (2006) 298:149–54. 10.1016/j.ydbio.2006.06.023
125.
Meijs-RoelofsHMVan CappellenWAVan LeeuwenECKramerP. Short- and long-term effects of an LHRH antagonist given during the prepubertal period on follicle dynamics in the rat. J Endocrinol. (1990) 124:247–53. 10.1677/joe.0.1240247
126.
FlawsJAAbbudRMannRJNilsonJHHirshfieldAN. Chronically elevated luteinizing hormone depletes primordial follicles in the mouse ovary. Biol Reprod. (1997) 57:1233–7. 10.1095/biolreprod57.5.1233
127.
VaithiyanathanKLiewSHZerafaNGamageTCookMO'reillyLAet al. BCL2-modifying factor promotes germ cell loss during murine oogenesis. Reproduction (2016) 151:553–62. 10.1530/REP-15-0561
128.
ChapmanCCreeLShellingAN. The genetics of premature ovarian failure: current perspectives. Int J Womens Health (2015) 7:799–810. 10.2147/IJWH.S64024
129.
QinYJiaoXSimpsonJLChenZJ. Genetics of primary ovarian insufficiency: new developments and opportunities. Hum Reprod Update (2015) 21:787–808. 10.1093/humupd/dmv036
130.
RossettiRFerrariIBonomiMPersaniL. Genetics of primary ovarian insufficiency. Clin Genet. (2017) 91:183–98. 10.1111/cge.12921
Summary
Keywords
acyl ghrelin, des-acyl ghrelin, ovarian follicles, reproductive success, RNA-seq
Citation
Sominsky L, Goularte JF, Andrews ZB and Spencer SJ (2019) Acylated Ghrelin Supports the Ovarian Transcriptome and Follicles in the Mouse: Implications for Fertility. Front. Endocrinol. 9:815. doi: 10.3389/fendo.2018.00815
Received
29 August 2018
Accepted
27 December 2018
Published
15 January 2019
Volume
9 - 2018
Edited by
Katja Teerds, Human and Animal Physiology, Wageningen University, Netherlands
Reviewed by
Li Meng, South China Agricultural University, China; Francisco Gaytán, Universidad de Córdoba, Spain
Updates
Copyright
© 2019 Sominsky, Goularte, Andrews and Spencer.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luba Sominsky luba.sominsky@rmit.edu.au
This article was submitted to Reproduction, a section of the journal Frontiers in Endocrinology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.