Abstract
Ovarian follicular development and ovulation are complex and tightly regulated processes that involve regulation by microRNAs (miRNAs). We previously identified differentially expressed mRNAs between human cumulus granulosa cells (CGCs) from immature early antral follicles (germinal vesicle - GV) and mature preovulatory follicles (metaphase II - M2). In this study, we performed an integrated analysis of the transcriptome and miRNome in CGCs obtained from the GV cumulus-oocyte complex (COC) obtained from IVM and M2 COC obtained from IVF. A total of 43 differentially expressed miRNAs were identified. Using Ingenuity IPA analysis, we identified 7288 potential miRNA-regulated target genes. Two hundred thirty-four of these target genes were also found in our previously generated ovulatory gene library while exhibiting anti-correlated expression to the identified miRNAs. IPA pathway analysis suggested that miR-21 and FOXM1 cooperatively inhibit CDC25A, TOP2A and PRC1. We identified a mechanism for the temporary inhibition of VEGF during ovulation by TGFB1, miR-16-5p and miR-34a-5p. The linkage bioinformatics analysis between the libraries of the coding genes from our preliminary study with the newly generated library of regulatory miRNAs provides us a comprehensive, integrated overview of the miRNA-mRNA co-regulatory networks that may play a key role in controlling post-transcriptomic regulation of the ovulatory process.
Similar content being viewed by others
Introduction
Ovarian follicular development and ovulation in mammals are highly complex and tightly regulated processes that involve the selection of a dominant follicle, reactivation of oocyte meiosis, rupture of the follicle wall, cumulus-oocyte expansion and tissue remodeling to form the corpus luteum. Oocyte growth and maturation are linked to the developmental properties of somatic cell populations in the expanding follicle as bi-directional signaling that occurs between oocytes and granulosa cells1. The genes involved in follicle maturation and ovulation are the subject of growing interest and have been well-studied by several laboratories, including ours2,3.
MicroRNAs (miRNAs) are small (~22-nucleotide), non-coding regulatory RNA molecules encoded by plants, animals and some viruses that regulate a variety of developmental and physiological processes (reviewed in4,5).
In the last few years, by virtue of the application of high-throughput sequencing techniques, it became clear that eukaryotes transcribe up to 90% of their genomic DNA. However, only 1-2% of these transcripts are protein-coding genes6.
Several studies have found specific miRNAs that target important ovarian molecules such as progesterone receptor (PGR)7, aromatase (CYP19A1), follicle stimulating hormone receptor (FSHR)8 and the FYN pathway9. miRNAs were suggested to function as RNA-based paracrine factors in the crosstalk between oocytes and support cells of the follicle10. Furthermore, the expression of a subset of miRNAs has been shown to vary during follicle/luteal development11,12,13. Effects on granulosa cell function and/or gene expression, largely in mice, have been reported for a handful of miRNAs, and these include effects on cell survival (miR-21 and miR-23a), proliferation (miR-145, miR-503, and miR-224), estradiol production (miR-224, miR-383, and miR-378) and terminal differentiation (miR-132 and miR-212) of cultured cells (reviewed in14,15).
In humans, miRNAs have already been profiled in cumulus versus mural granulosa cells8, Polycystic Ovary Syndrome (PCOS) patients16 and cumulus cells according to oocyte nuclear maturity17. miRNAs have been suggested as potential biomarkers in IVF18.
In the past several years, miRNA and mRNA integrated analysis (MMIA) has become a tool for examining the biological functions of miRNA expression. As the biological effects of miRNAs are due to their modulation of target RNA expression, accurate miRNA target prediction is essential to any study of miRNA function.
To understand the molecular and regulatory mechanisms necessary for follicular maturation and ovulation, we aimed to uncover the miRNA transcriptome involved in these processes. Recent advances in genomics allow a systematic approach to identify critical genes and regulated pathways involved in oocyte maturation and ovulation. Using whole transcriptome sequencing, we recently identified mRNAs that are differentially expressed between immature mid-antral follicles and mature pre-ovulatory follicles3. The resulting database provides unprecedented insight into the processes and pathways involved in follicular maturation and ovulation. To complete the identification of factors involved in the ovulatory process, the aim of this work was to generate a library of global miRNAs involved in this process and to link the new ovulatory miRNA library with the previously described mRNA library. This enabled us to identify new regulatory mechanisms responsible for the final follicular maturation and ovulatory processes.
Results
miRNA profile differences in compact and expanded cumulus cells
To elucidate the roles of miRNAs in human folliculogenesis and ovulation, we generated a library of regulated miRNAs during the final stage of follicular maturation and ovulation. We used NanoString on RNA extracted from two groups of cumulus granulosa cells: (1) Compact cumulus cells surrounding GV oocytes that were acquired during IVM treatment (CCGV n = 4) and (2) Expanded cumulus cells enclosing M2 oocytes that were acquired during IVF treatment (CCM2, n = 3). One CCGV sample was an outlier and was removed from further downstream analysis. Thus, the final number of samples in each group was three.
A Venn diagram of these groups is shown in Fig. 1B, illustrating the numbers of selective and co-expressed miRNAs in the two conditions, CCGV and CCM2. A total of 45 miRNAs were detected in CCM2 and 11 miRNAs were detected in CCGV. Most notably, 35 miRNAs were exclusively detected in the CCM2 group; only 1 was exclusively detected in the CCGV group and 10 miRNAs were shared between the two conditions (Table 1). hsa-miR-378e was exclusively expressed in the CCGV group, although the differential expression level was not statistically significant.
With NanoString data inspected, all quality control parameters were within the acceptable ranges. Unsupervised hierarchical cluster analysis based on the expression of all detected miRNAs shows that the biological replicates cluster to their corresponding groups (Fig. 1A).
Almost all of the expressed miRNAs were significantly differentially expressed (43 of the 46 identified miRNAs) between CCGV and CCM2, with a fold change >2 and a false discovery rate (FDR) < 5%. Of these, 25 represent unique seed sequences. Our results indicate that all differentially expressed miRNAs were up-regulated in CCM2 cells.
Validation of selected miRNA expression levels as estimated by qPCR
We used qPCR analysis to validate the results obtained by NanoString. We selected three miRNAs that were identified as differentially expressed by NanoString. qPCR was performed in three different experiments, as shown as fold induction in Fig. 1C, and demonstrated agreement with NanoString relative expression data.
Identification of differentially expressed miRNA target genes – in silico analysis
Because miRNAs act by regulating the expression of target genes, precise miRNA target prediction is important for elucidating miRNA function. The identification of miRNA target genes was performed using the microRNA target filter tool of QIAGEN’s Ingenuity Pathway Analysis software, taking into account only experimentally validated and highly confident predictions of miRNA target interactions. We found 7288 putative miRNA target genes of which 727 (294 unique) were experimentally validated (Supp. Tables S3 and S4).
Identification of differentially expressed miRNAs that negatively correlated with ovulatory target genes - experimental analysis
To further dissect the putative ovulatory miRNA targets, we combined our previously generated mRNA library3 and miRNA expression data using the microRNA target filter tool of QIAGEN’s Ingenuity Pathway Analysis19. A Venn diagram depicting the relationships between the putative miRNA targets and experimentally differentially expressed mRNAs is depicted in Fig. 1D.
We chose to focus on negatively correlated targets. Since all of the differentially expressed CCM2 miRNAs were upregulated, all chosen paired mRNA targets were downregulated. The 696 downregulated mRNAs in the CCM2 group (compared to CCGV group) were paired with the 7288 putative miRNA targets. Of the 696 downregulated mRNAs, 234 were targets of 22 differentially expressed miRNAs (of the 25 having unique seed sequences). In total, these 234 genes formed 479 miRNA-target gene pairs with an inverse correlation of expression (Table 2).
To understand the putative roles of these miRNA target genes, two lists of putative regulated genes were generated. The first list included all of the putative miRNA targets (in silico analysis) and the second list included only the negatively correlated differentially expressed targets (experimental analysis). This contributes to both a broad and targeted view of the ontologies and pathways related to the roles of miRNAs in ovulation. Gene Ontology (GO) and pathway analyses were performed using the GeneAnalytics tool20.
Gene Ontology and pathway classification of the in silico analysis of miRNA target genes
For the unfiltered miRNA targets (7288 genes), we identified 78 high-score GO terms (Fig. 2A and Supp. Table S5). The GO terms included 9 “transcription regulation” ontologies and 15 “signaling” ontologies, including neurotropin TRK signaling, FGF receptor signaling, EGF receptor signaling, TGFB receptor signaling, WNT signaling, VEGF signaling, insulin receptor signaling and RAS signaling. Other important processes included axon guidance, nervous system development, negative regulation of cell proliferation, in utero embryonic development and angiogenesis and cell migration.
We identified 144 high-score pathways that included 67 signaling pathways (Fig. 2B and Supp. Table S6). The top scoring signaling pathways were ERK signaling, HGF signaling, GPCR signaling, MAPK signaling, P70-S6K signaling and RAS signaling. Other pathways included NFAT and cardiac hypertrophy, proteoglycans in cancer, focal adhesion, endocytosis and axon guidance.
Gene Ontology and pathway classification of the negatively correlated miRNA target genes
For the subset of negatively correlated miRNA targets (234 genes), only four high scoring and four medium scoring GO terms were identified. The GO terms included DNA replication, ascending aorta morphogenesis, extracellular matrix organization and mitotic cell cycle. The top medium scoring GO terms included collagen fibril, basement membrane organization and signal transduction (Fig. 2C and Supp. Table S7). The nine high scoring pathways included RB in cancer, integrated pancreatic cancer pathway, diseases associated with O-glycosylation, cell cycle, E2F-mediated DNA replication, degradation of extracellular matrix and ERK signaling (Fig. 2D and Supp. Table S8).
Upstream regulator analysis of negatively correlated miRNA target genes
We used Ingenuity’s IPA upstream regulator analytic feature to identify a cascade of upstream transcriptional regulators that can explain the observed gene expression changes of the negatively correlated genes beyond the identified miRNAs. These regulators can be transcription factors (TFs) and any gene or small molecule that has been observed experimentally to affect gene expression. TF and miRNA may mutually regulate each another or regulate a shared target gene and, as such, enhance the robustness of gene regulation. Hence, TF-miRNA regulatory network analysis will be helpful to decipher gene expression regulatory mechanisms of the ovulatory process.
A total of 68 significant upstream regulators were identified (Z-score > ±2, Supp. Table S9). The top ten inhibited regulators and activated regulators are presented in Tables 3 and 4, respectively.
As expected, among the activated regulators, we identified 5 miRNAs, 3 of which were also differentially expressed in our analysis (Supp. Table S9).
We plotted several of the upstream regulators and their targets, as shown in Fig. 3. We identified miR-21 to be involved in the regulation of CDC25A, DDAH1, ECT2, MSH2, NUSAP1, PBK, PRC1, RAD51AP1, STMN1 and TOP2A (Fig. 3A). Interestingly, most genes are involved in reproductive disease (highlighted in pink).
We identified let-7 to be involved in the regulation of AURKB, BRCA2, CDC25A, CDC6, CHEK1, EZH2, FANCD2, MCM3, PPP1R12B, RRM2 and THBS1 (Fig. 3B).
Among the notable interactions found was the reciprocal crosstalk between miR-21 (activated) and FOXM1 (inhibited) signaling. Both transcriptional regulators inhibit STMN1, TOP2A, CDC25A and PRC1 (Fig. 3C).
Also noted was the reciprocal crosstalk between TGFB1 (inhibited), miR-16-5p and miR-34a-5p (both activated) signaling (Fig. 4). All three transcription regulators inhibit VEGFA expression observed during the ovulatory process. TGFB1 (inhibited) activates E2F1, while miR-34a (activated) inhibits E2F1 signaling to the net effect of E2F1 inhibition. TGFB1 also inhibits mir-34a.
Interestingly, both upstream regulators, miR-16-5p (seed sequence of miR-424-5p) and miR-34a-5p, are also upregulated in granulosa cells during ovulation (Table 2 and Fig. 1C).
Functional validation of upstream regulator analysis results
We chose several of the putative miRNA target genes and regulators that were identified by the upstream regulator analysis and performed qPCR analysis. In CCM2 cells, FOXM1 and TOP2A were downregulated, whereas miR-21 was upregulated (Figs 3C and 4B, fold induction). Additionally, miR-34a-5p was upregulated, while TOP2A and CD47 were downregulated (Fig. 4, fold induction).
Discussion
A number of previous studies have examined mRNA and/or miRNA expression in cumulus granulosa cells. However, no study has examined global miRNA expression in human cumulus granulosa cells during final follicular maturation and ovulation. Moreover, this is the first study to link the miRNA-target gene pairs with an inverse correlation of expression. This approach allows the identification of the regulated mRNA-miRNA networks during final follicular maturation and ovulation.
Differential expression analysis revealed dramatic changes in miRNA expression during the CC maturation process; all differentially expressed miRNAs were up-regulated in CCM2 cells.
Compared to published miRNA expression in CCs obtained from M2 COC during IVF treatment8, we observed that 8 of the 10 most abundant miRNAs were also present in the CCM2 sample and 43 of the 46 expressed miRNAs were identified by Velthut-Meikas et al.
Forty-two of the differentially expressed miRNAs were previously identified in granulosa cells. Among them, 14 are involved in steroidogenesis or PCOS, 14 are involved in apoptosis (eight of which belong to the let-7 family), 12 have been implicated in bovine follicular development and 20 were described as differentially expressed between human cumulus and mural granulosa cells. miR-4454 was not previously described in granulosa cells and was recently implicated in the pathogenesis of cartilage degeneration by promoting inflammatory, catabolic, and cell death activity in chondrocytes21. We believe that miR-4454 is implicated in the control of final follicular maturation perhaps by promoting similar proinflammatory-like activities in granulosa cells.
The differentially expressed miRNAs were predicted to regulate 696 target genes that were differentially expressed in our previously generated ovulatory cDNA library, of which, 234 displayed an inverse correlation.
During the final stages of ovulation, several upregulated processes were previously identified by the analysis of global mRNA expression, such as cellular movement, inflammatory response, immune cell trafficking, tissue development and lipid metabolism3,22,23. However, several processes, such as tissue morphology, DNA replication and cell cycle, are downregulated3,22,23. These downregulated processes were also represented in our analyses of miRNA targets. We observed that a large number of the putative miRNA target genes are related to cell cycle and DNA replication processes. These are probably downstream to the FSH blocking effect following the LH surge, thus leading to decreased proliferation24. Interestingly, all of the differentially expressed miRNAs were upregulated in our analysis. This result suggests that the main role of miRNA in the late ovulatory process is more toward control and inhibition of granulosa cell proliferation and less toward other processes found when analyzing gene, rather than miRNA, expression.
The top identified pathways among the miRNA negatively correlated targets include several pathways that pertain to cellular proliferation and survival, including “RB in cancer”, “E2F regulation of DNA replication” and “cell cycle”. Among them, the retinoblastoma protein (RB) was identified as an important factor in ovarian physiology. Conditional deletion of the RB gene in murine ovarian granulosa cells leads to increased follicular recruitment and premature ovarian failure25.
We also identified the O-glycosylation of proteins pathway, which might pertain to cumulus matrix maintenance since cumulus complexes from O-glycan-deficient oocytes were smaller and contained fewer CCs, although fertility was not impaired26.
Several pathways involved in ovulation were identified in our analyses, including the ERK signaling pathway27, HGF signaling28, MAPK pathway27, P70-S6K signaling29, neurotropin TRK signaling30, FGF receptor signaling31, EGF receptor signaling31, TGFB receptor signaling32, WNT signaling33, VEGF signaling34, insulin receptor signaling35, and RAS signaling36. These findings strengthen and validate our miRNA results.
As expected, more potential novel pathways, including NFAT (nuclear factor of activated T-cells) and PAK (p21-activated kinase) signaling pathways, were found in the broader analysis, which included all miRNA targets (and not only the anti-correlated ones).
NFAT has not been described in granulosa cells, although it has been implicated in GnRH signaling37. NFAT was also described as a regulator of COX2 in endometrial stromal cells38 and other tissues and may play a role in prostaglandin regulation during ovulation.
PAK signaling was implicated in cell migration by altering actin cytoskeletal dynamics downstream to Rac/Cdc42, although this activity was not described in granulosa cells and may be important during cumulus expansion and ovulation39.
miR-21 is highly expressed and upregulated by LH in murine granulosa cells40 and ovine follicles12. miR-21 suppression results in granulosa cell apoptosis40. Our results further corroborate miR-21 upregulation in humans during the ovulatory process and, by combining our negatively correlated ovulatory gene expression data, we suggest a mechanism of action for miR-21 function in granulosa cells. We found 10 different transcripts regulated by miR-21 that were downregulated in our library, all of which are involved in the cell cycle and apoptosis. These findings expand our knowledge of the processes that lead to the observed inhibition of proliferation on granulosa cells during the final stages of the ovulatory process. Some of the targets were already identified in granulosa cells, including CDC25a (apoptosis, cell cycle41), MSH2 (apoptosis42), PRC1 (apoptosis43) and STMN1 (cell migration, apoptosis44). Other targets that were not previously described in granulosa cells include TOP2A (apoptosis, cell cycle), DDAH1 (apoptosis), ECT2 (cell cycle), RAD51AP1 (cell cycle), PBK (apoptosis) and NUSAP1 (cell migration, apoptosis). The known and previously unknown targets provide insight into miR-21 activity in granulosa cells. The negative correlation between TOP2A and miR-21 was experimentally validated (Figs 1C and 4B).
Forkhead box M1 (FOXM1), a member of the large family of Forkhead box transcription factors, is highly expressed in proliferating cells and plays pivotal roles in embryonic and fetal development, DNA replication and mitosis45. The role of FOXM1 in ovulation and granulosa cells has not yet been established. FOXM1 is downregulated in granulosa cells obtained from obese patients undergoing IVF compared to lean patients46 and was also found as an inhibited upstream regulator of negatively correlated miRNA target genes in our analysis. Our observation suggests that the reciprocal activation of miR-21 and suppression of FOXM1 synergize to cause the observed inhibition of granulosa cell proliferation during the final follicular maturation just prior to ovulation. Interestingly, high expression of miR-21 is implicated in increased proliferation in cancer and inflammation47. However, the inhibition of breast cancer cell growth by 3, 3′-diindolylmethane is mediated by the upregulation of miR-21 and CDC25A and FOXM1 suppression48, consistent with that observed in cumulus cells during ovulation3. These in silico findings between FOXM1, TOP2A, and miR-21 were experimentally validated (Figs 1C and 4B).
The let-7 family of miRNAs is believed to act as tumor suppressors49 and has been implicated in granulosa cell atresia in porcine50. This family was among the most upregulated miRNAs during final follicular maturation. Several of its negatively correlated targets were previously identified in granulosa cells, including AURKB (cell cycle51), CDC25A (apoptosis, cell cycle41), EZH2 (histone methylation52), and THBS1 (adhesion, angiogenesis53). Other targets that were not described in granulosa cells include BRCA2 (DNA repair), CDC6 (DNA replication), FANCD2 (DNA repair), MCM3 (DNA replication), PPP1R12B (myosin phosphatase), RRM2 (DNA synthesis) and CHEK1 (cell cycle). Most negatively correlated let-7 targets are related to DNA replication/repair and cell cycle progression. Thus, the putative role of let-7 upregulation during final follicular maturation may be in mediating the observed inhibition of granulosa cell proliferation.
We determined that miR-34a is upregulated in CCM2 compared to CCGV. miR-34a is a tumor suppressor gene that suppresses ovarian cancer proliferation and motility by targeting AXL receptor tyrosine kinase54. miR-34a levels are increased in sheep granulosa cells during the follicular to luteal transition12. In addition, miR-34a may be involved in granulosa cell apoptosis in pig ovaries by targeting the inhibin B gene55. The integrated mRNA-miRNA analysis suggests that miR-34a acts through the regulation of E2F1 (transcription), CD47 (adhesion), MCM3 (proliferation), VEGFA (angiogenesis) and TGFB1 signaling during final follicular maturation and luteinization.
Furthermore, the roles of the TGFB1 signaling pathways in folliculogenesis have been extensively studied56. The Ingenuity’s IPA upstream regulator analytic identified 38 inhibited TGFB1 targets of the 234 putative miRNA negatively correlated targets, making it the most significant miRNA-regulated pathway. When TGFB1 targets were plotted along with miR-16-5p (seed of miRNA-424-3p) and miR-34a targets, we determined that VEGFA expression was affected by all of these regulators. It was previously shown by us3 and others2 that VEGFA expression is reduced just prior to ovulation and increased in the luteal phase. This identified regulatory network suggests a novel cooperative mechanism for VEGF transient downregulation.
In conclusion, integrated analysis between the expression of coding genes from our preliminary study with the newly generated library of regulatory miRNAs enabled us to better understand the regulation of coding ovulatory genes by non-coding transcripts and their function from the single-molecule level to whole pathways.
Methods
Ethical approval
The current study was approved by the Sheba Medical Center Institutional Review Board (IRB) Committee (ethical approval number 8707-11-SMC), and written informed consent was obtained from each patient.
All experiments were performed in accordance with relevant guidelines and regulations.
IVF protocol
Ovarian stimulation was carried out as previously described57,58 according to the “short antagonist” protocol. The protocol consisted of controlled ovarian hyperstimulation with recombinant FSH (rFSH), either Gonal-F; Merck Serono or Puregon Pen; Schering Plough) and human menopausal gonadotropin (HMG; Menopur; Ferring) followed by the addition of ovarian suppression with GnRH antagonists (0.25 mg/day, Cetrorelix, Cetrotide; Serono International, SR) when the leading follicle was more than 12 mm in diameter. When three or more follicles exceeded 18 mm in diameter, 250 µg of hCG (Ovitrelle; Merck Serono) was administered to trigger ovulation. Transvaginal follicular aspiration was performed 35 hours later with ultrasound guidance.
For NanoString analysis, three women (ages 25–35 years) undergoing IVF donated cumulus cells from one M2 COC (see Supp. Table S1 for details).
For the miRNA NanoString validation experiments by qPCR, cumulus cells were obtained from 3–4 different women and pooled to generate a single replicate. Each woman donated CGCs from one M2 COC. The women’s ages and infertility etiologies were similar to those described in Table S1.
For functional validation of upstream regulator analysis, cumulus cells were obtained from 3–4 different women and pooled to generate a single replicate. Each woman (total −21 women) donated CGCs of one M2 or one GV COC. Due to the limited availability of GV COC obtained from IVM, we used cumulus cells of GV COC obtained during IVF. The age of the women was 28–40 years (average 36.5 years). The etiology of infertility included male factor and PGD.
IVM protocol
Four women (ages 28–40 years) undergoing IVM were selected for this study (see Supp. Table S1 for details). In the final analysis, only 3 women were included; one woman (Patient code 59 A in Table S1) was excluded because she was an outlier. Each woman donated one single COC. IVM cycles were carried out as previously described59. Briefly, sonographic assessment of the antral follicle count and of endometrial thickness was carried out on day 3 of a spontaneous menstrual cycle. The serum concentrations of estradiol and progesterone were also determined. Next, 150 IU/day rFSH were administered to the patients for 3 days. A second evaluation was performed on day 6. An injection of 10000 IU hCG (Pregnyl; Organon, Oss, Holland) was administered subcutaneously when the endometrial thickness was ≥5 mm and the leading follicle was at least 12 mm. Oocyte retrieval was carried out 36 hours later. Compact CCs were obtained from surplus germinal vesicle (GV) oocytes that were acquired during IVM treatment (CCGV group, provided by Dr. Ruben Fadini from the Biogenesi Reproductive Medicine Center, Istituti Clinici Zucchi, Monza, Italy).
For NanoString analysis, four women (ages 28–40 years) undergoing IVM were selected (see Supp. Table S1 for details). In the final analysis, only 3 women were included; one woman (Patient code 59 A in Table S1) was excluded because she was an outlier. Each woman donated one single COC. For the miRNA NanoString validation experiments, cumulus cells were obtained from 3–4 different women and pooled to generate a single replicate. Each woman donated CGCs from one GV COC. The women’s ages and infertility etiologies were similar to those described in Table S1.
Cumulus granulosa cell collection
CGCs were obtained through oocyte denudation in the course of intracytoplasmic sperm injection (ICSI) procedures. After oocyte retrieval, CGCs of each oocyte were removed using hyaluronidase (SAGE) and a glass denudation pipette (Swemed). The CGCs were washed in PBS and centrifuged at 5000 × g for 5 minutes at room temperature. The resulting pellets were stored at −80 °C until RNA isolation.
RNA extraction
Total RNA was extracted from seven independent CCs samples of seven different women obtained from compact GV (4 samples) and expanded M2 of single COC (3 samples) using a Micro RNA Isolation Kit (Zymo Research Corp., CA USA) according to the manufacturer’s instructions. RNA purity and concentration were assessed using a NanoDrop spectrophotometer (NanoDrop 2000C, Thermo Fisher Scientific).
Quantitative real-time PCR
cDNA synthesis for the detection of mature miRNAs was performed with the miScript Reverse Transcription Kit (QIAGEN) using a blend of oligo(dT) and random primers (iScript™ cDNA Synthesis Kit). cDNA synthesis for the quantification of mRNA expression was performed with the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). miRNA expression was determined using the miScript SYBR Green Kit (QIAGEN) and mRNA expression was determined using Fast SYBR Green Master Mix (Applied Biosystems). The StepOnePlus Real-Time PCR System (Applied Biosystems) was used to detect amplification.
Expression levels were normalized to ACTB (for mRNA) and RNU6B (for miRNA). qPCR results were analyzed with StepOne software. Relative gene expression was calculated using the delta-delta Ct method60. Details of the primers are shown in Table S2.
NanoString and bioinformatics sample analysis
Evaluation of miRNA expression was carried out using NanoString technology, which is not based on sequencing but rather based on a digital molecular barcoding system. Each barcode is attached to a single target-specific probe corresponding to a specific miRNA. The output from this technology is the digital count of 800 tested miRNAs.
A total of 3 μl (20–120 ng) of each RNA sample was prepared according to the manufacturer’s instructions. Mature miRNAs were ligated to a species-specific tag sequence (miRtag) via a thermally controlled splinted ligation. miRNA assays were performed according to the NanoString miRNA Assay Manual. Hybridizations were carried out by combining 5 µl of each miRNA assay with 20 µl of nCounter Reporter probes in hybridization buffer and 5 µl of nCounter Capture probes for a total reaction volume of 30 µl. miRNA expression was assessed using the NanoString nCounter system (NanoString Technologies, Seattle, WA, USA), which enables multiplexed direct digital counting of 800 human miRNA molecules. The raw count data from the 6 samples were analyzed using the edgeR61 and limma R62 packages. First, TMM normalization was applied to the data followed by voom transformation. To assess the number of miRNAs expressed under each condition, we required that probes be detectable in all replicates of the tested condition. As detectable, we defined probes having a normalized expression level that exceeded that of the “highest expressed” negative control probe. This resulted in 54 expressed miRNAs. Next, linear models were used to assess differential expression using the limma method. Pairwise comparisons were performed using moderated t statistics. To perform multiple group comparisons, one-way ANOVA was applied, except the residual mean squares were moderated between genes (see the limma package for details62). miRNAs were considered differentially expressed if their FDR was <5%.
The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.
Bioinformatics analysis of the library
The resultant differentially expressed miRNAs were analyzed through the use of IPA (Ingenuity Pathway Analysis QIAGEN Inc. https://www.qiagenbioinformatics.com/products/ingenuitypathway-analysis)19. Putative miRNA targets found using Ingenuity miRNA target filter for experimentally validated and putative predicted targets (high confidence level) – in silico analysis.
Using whole transcriptome sequencing, we previously generated an mRNA ovulatory genes library3. The library was comprised of the same groups as in the current study (CCGV and CCM2), encompassing all differentially expressed mRNAs between immature mid-antral follicles and mature pre-ovulatory follicles.
Our previously generated mRNA library and the miRNA putative targets were combined using the microRNA target filter to generate a subset of the negatively correlated miRNA targets – experimental analysis.
Both mRNA target lists (in silico analysis target list and experimental analysis target list) were analyzed for GO terms and pathways using the GeneAnalytics tool20.
The DE miRNAs and negatively correlated targets were also analyzed using the Ingenuity Upstream Regulator Analysis tool to predict upstream molecules, including miRNA and transcription factors, that may have caused the observed gene expression changes.
Statistics
Each experiment was performed at least three times. Data are expressed as the mean ± standard error of the mean (SEM) and were evaluated using Student’s t-test with a two-tailed distribution, with two samples equaling variance, or with ANOVA (ANalysis Of VAriance) for more than two variances using the post hoc Tukey test assuming equal variance, or the Games–Howell test for unequal variance. When appropriate, the Kruskal–Wallis non-parametric comparison test was used. For all statistical analyses, SPSS 22 software (IBM, Armonk, NY, USA) was used. P-value < 0.05 was considered statistically significant.
References
Li, Q., McKenzie, L. J. & Matzuk, M. M. Revisiting oocyte-somatic cell interactions: in search of novel intrafollicular predictors and regulators of oocyte developmental competence. Mol Hum Reprod 14, 673–678, https://doi.org/10.1093/molehr/gan064 (2008).
Wissing, M. L. et al. Identification of new ovulation-related genes in humans by comparing the transcriptome of granulosa cells before and after ovulation triggering in the same controlled ovarian stimulation cycle. Human reproduction (Oxford, England) 29, 997–1010, https://doi.org/10.1093/humrep/deu008 (2014).
Yerushalmi, G. M. et al. Characterization of the human cumulus cell transcriptome during final follicular maturation and ovulation. Mol Hum Reprod. https://doi.org/10.1093/molehr/gau031 (2014).
Bartel, D. P. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297 (2004).
Plasterk, R. H. Micro RNAs in animal development. Cell 124, 877–881, https://doi.org/10.1016/j.cell.2006.02.030 (2006).
Kaikkonen, M. U., Lam, M. T. Y. & Glass, C. K. Non-coding RNAs as regulators of gene expression and epigenetics. Cardiovascular Research 90, 430–440, https://doi.org/10.1093/cvr/cvr097 (2011).
Toms, D., Xu, S., Pan, B., Wu, D. & Li, J. Progesterone receptor expression in granulosa cells is suppressed by microRNA-378-3p. Mol Cell Endocrinol 399, 95–102, https://doi.org/10.1016/j.mce.2014.07.022 (2015).
Velthut-Meikas, A. et al. Research resource: small RNA-seq of human granulosa cells reveals miRNAs in FSHR and aromatase genes. Molecular endocrinology 27, 1128–1141, https://doi.org/10.1210/me.2013-1058 (2013).
Ninio-Many, L., Grossman, H., Shomron, N., Chuderland, D. & Shalgi, R. microRNA-125a-3p reduces cell proliferation and migration by targeting Fyn. J Cell Sci 126, 2867–2876, https://doi.org/10.1242/jcs.123414 (2013).
Hawkins, S. M. & Matzuk, M. M. Oocyte-somatic cell communication and microRNA function in the ovary. Ann Endocrinol (Paris) 71, 144–148, https://doi.org/10.1016/j.ando.2010.02.020 (2010).
Fiedler, S. D., Carletti, M. Z., Hong, X. & Christenson, L. K. Hormonal regulation of MicroRNA expression in periovulatory mouse mural granulosa cells. Biol Reprod 79, 1030–1037, https://doi.org/10.1095/biolreprod.108.069690 (2008).
McBride, D. et al. Identification of miRNAs associated with the follicular-luteal transition in the ruminant ovary. Reproduction (Cambridge, England) 144, 221–233, https://doi.org/10.1530/rep-12-0025 (2012).
Schauer, S. N., Sontakke, S. D., Watson, E. D., Esteves, C. L. & Donadeu, F. X. Involvement of miRNAs in equine follicle development. Reproduction (Cambridge, England) 146, 273–282, https://doi.org/10.1530/rep-13-0107 (2013).
Donadeu, F. X., Schauer, S. N. & Sontakke, S. D. Involvement of miRNAs in ovarian follicular and luteal development. The Journal of endocrinology 215, 323–334, https://doi.org/10.1530/JOE-12-0252 (2012).
Imbar, T. & Eisenberg, I. Regulatory role of microRNAs in ovarian function. Fertility and sterility 101, 1524–1530, https://doi.org/10.1016/j.fertnstert.2014.04.024 (2014).
Liu, S. et al. Altered microRNAs expression profiling in cumulus cells from patients with polycystic ovary syndrome. J Transl Med 13, 238, https://doi.org/10.1186/s12967-015-0605-y (2015).
Moreno, J. M. et al. Follicular fluid and mural granulosa cells microRNA profiles vary in in vitro fertilization patients depending on their age and oocyte maturation stage. Fertility and sterility, https://doi.org/10.1016/j.fertnstert.2015.07.001 (2015).
Scalici, E. et al. Circulating microRNAs in follicular fluid, powerful tools to explore in vitro fertilization process. Scientific reports 6, 24976, https://doi.org/10.1038/srep24976 (2016).
Kramer, A., Green, J., Pollard, J. Jr & Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 30, 523–530, https://doi.org/10.1093/bioinformatics/btt703 (2014).
Ben-Ari Fuchs, S. et al. GeneAnalytics: An Integrative Gene Set Analysis Tool for Next Generation Sequencing, RNAseq and Microarray Data. OMICS 20, 139–151, https://doi.org/10.1089/omi.2015.0168 (2016).
Nakamura, A. et al. Identification of microRNA-181a-5p and microRNA-4454 as mediators of facet cartilage degeneration. JCI Insight 1, https://doi.org/10.1172/jci.insight.86820 (2016).
Lee, Y. S. et al. Extensive effects of in vitro oocyte maturation on rhesus monkey cumulus cell transcriptome. Am J Physiol Endocrinol Metab 301, E196–209, https://doi.org/10.1152/ajpendo.00686.2010 (2011).
Ouandaogo, Z. G. et al. Differences in transcriptomic profiles of human cumulus cells isolated from oocytes at GV, MI and MII stages after in vivo and in vitro oocyte maturation. Human reproduction (Oxford, England), https://doi.org/10.1093/humrep/des172 (2012).
Yong, E. L., Baird, D. T. & Hillier, S. G. Mediation of gonadotrophin-stimulated growth and differentiation of human granulosa cells by adenosine-3′,5′-monophosphate: one molecule, two messages. Clin Endocrinol (Oxf) 37, 51–58 (1992).
Andreu-Vieyra, C., Chen, R. & Matzuk, M. M. Conditional deletion of the retinoblastoma (Rb) gene in ovarian granulosa cells leads to premature ovarian failure. Molecular endocrinology 22, 2141–2161, https://doi.org/10.1210/me.2008-0033 (2008).
Ploutarchou, P., Melo, P., Day, A. J., Milner, C. M. & Williams, S. A. Molecular analysis of the cumulus matrix: insights from mice with O-glycan-deficient oocytes. Reproduction (Cambridge, England) 149, 533–543, https://doi.org/10.1530/REP-14-0503 (2015).
Fan, H. Y. et al. MAPK3/1 (ERK1/2) in ovarian granulosa cells are essential for female fertility. Science 324, 938–941, https://doi.org/10.1126/science.1171396 (2009).
Zachow, R. & Uzumcu, M. The hepatocyte growth factor system as a regulator of female and male gonadal function. The Journal of endocrinology 195, 359–371, https://doi.org/10.1677/JOE-07-0466 (2007).
Alam, H. et al. Follicle-stimulating hormone activation of hypoxia-inducible factor-1 by the phosphatidylinositol 3-kinase/AKT/Ras homolog enriched in brain (Rheb)/mammalian target of rapamycin (mTOR) pathway is necessary for induction of select protein markers of follicular differentiation. The Journal of biological chemistry 279, 19431–19440, https://doi.org/10.1074/jbc.M401235200 (2004).
Vera, C., Tapia, V., Vega, M. & Romero, C. Role of nerve growth factor and its TRKA receptor in normal ovarian and epithelial ovarian cancer angiogenesis. Journal of ovarian research 7, 82, https://doi.org/10.1186/s13048-014-0082-6 (2014).
Chaves, R. N., de Matos, M. H., Buratini, J. Jr & de Figueiredo, J. R. The fibroblast growth factor family: involvement in the regulation of folliculogenesis. Reprod Fertil Dev 24, 905–915, https://doi.org/10.1071/rd11318 (2012).
Juengel, J. L. & McNatty, K. P. The role of proteins of the transforming growth factor-beta superfamily in the intraovarian regulation of follicular development. Hum Reprod Update 11, 143–160, https://doi.org/10.1093/humupd/dmh061 (2005).
Fan, H. Y. et al. Beta-catenin (CTNNB1) promotes preovulatory follicular development but represses LH-mediated ovulation and luteinization. Molecular endocrinology 24, 1529–1542, https://doi.org/10.1210/me.2010-0141 (2010).
Otani, N. et al. The vascular endothelial growth factor/fms-like tyrosine kinase system in human ovary during the menstrual cycle and early pregnancy. The Journal of clinical endocrinology and metabolism 84, 3845–3851, https://doi.org/10.1210/jcem.84.10.6025 (1999).
Brogan, R. S., Mix, S., Puttabyatappa, M., VandeVoort, C. A. & Chaffin, C. L. Expression of the insulin-like growth factor and insulin systems in the luteinizing macaque ovarian follicle. Fertility and sterility 93, 1421–1429, https://doi.org/10.1016/j.fertnstert.2008.12.096 (2010).
Wayne, C. M., Fan, H. Y., Cheng, X. & Richards, J. S. Follicle-stimulating hormone induces multiple signaling cascades: evidence that activation of Rous sarcoma oncogene, RAS, and the epidermal growth factor receptor are critical for granulosa cell differentiation. Molecular endocrinology 21, 1940–1957, https://doi.org/10.1210/me.2007-0020 (2007).
Binder, A. K., Grammer, J. C., Herndon, M. K., Stanton, J. D. & Nilson, J. H. GnRH regulation of Jun and Atf3 requires calcium, calcineurin, and NFAT. Molecular endocrinology 26, 873–886, https://doi.org/10.1210/me.2012-1045 (2012).
Abraham, F., Sacerdoti, F., De León, R., Gentile, T. & Canellada, A. Angiotensin II Activates the Calcineurin/NFAT Signaling Pathway and Induces Cyclooxygenase-2 Expression in Rat Endometrial Stromal Cells. PLoS ONE 7, e37750, https://doi.org/10.1371/journal.pone.0037750 (2012).
Edwards, D. C., Sanders, L. C., Bokoch, G. M. & Gill, G. N. Activation of LIM-kinase by Pak1 couples Rac/Cdc42 GTPase signalling to actin cytoskeletal dynamics. Nature cell biology 1, 253–259, https://doi.org/10.1038/12963 (1999).
Carletti, M. Z., Fiedler, S. D. & Christenson, L. K. MicroRNA 21 blocks apoptosis in mouse periovulatory granulosa cells. Biol Reprod 83, 286–295, https://doi.org/10.1095/biolreprod.109.081448 (2010).
Pan, H. A. et al. CDC25 protein expression and interaction with DAZL in human corpus luteum. Fertility and sterility 92, 1997–2003, https://doi.org/10.1016/j.fertnstert.2008.09.025 (2009).
Donadeu, F. X., Mohammed, B. T. & Ioannidis, J. A miRNA target network putatively involved in follicular atresia. Domest Anim Endocrinol 58, 76–83, https://doi.org/10.1016/j.domaniend.2016.08.002 (2017).
Ndiaye, K., Fayad, T., Silversides, D. W., Sirois, J. & Lussier, J. G. Identification of downregulated messenger RNAs in bovine granulosa cells of dominant follicles following stimulation with human chorionic gonadotropin. Biol Reprod 73, 324–333, https://doi.org/10.1095/biolreprod.104.038026 (2005).
Garcia-Rudaz, C. et al. Excessive ovarian production of nerve growth factor elicits granulosa cell apoptosis by setting in motion a tumor necrosis factor alpha/stathmin-mediated death signaling pathway. Reproduction (Cambridge, England) 142, 319–331, https://doi.org/10.1530/rep-11-0134 (2011).
Bella, L., Zona, S., Nestal de Moraes, G. & Lam, E. W. FOXM1: A key oncofoetal transcription factor in health and disease. Semin Cancer Biol 29, 32–39, https://doi.org/10.1016/j.semcancer.2014.07.008 (2014).
Merhi, Z. et al. Adiposity Alters Genes Important in Inflammation and Cell Cycle Division in Human Cumulus Granulosa Cell. Reprod Sci 22, 1220–1228, https://doi.org/10.1177/1933719115572484 (2015).
Sheedy, F. J. Turning 21: Induction of miR-21 as a Key Switch in the Inflammatory Response. Front Immunol 6, 19, https://doi.org/10.3389/fimmu.2015.00019 (2015).
Jin, Y. 3,3′-Diindolylmethane inhibits breast cancer cell growth via miR-21-mediated Cdc25A degradation. Mol Cell Biochem 358, 345–354, https://doi.org/10.1007/s11010-011-0985-0 (2011).
Shimizu, S. et al. The let-7 family of microRNAs inhibits Bcl-xL expression and potentiates sorafenib-induced apoptosis in human hepatocellular carcinoma. J Hepatol 52, 698–704, https://doi.org/10.1016/j.jhep.2009.12.024 (2010).
Cao, R. et al. Expression and preliminary functional profiling of the let-7 family during porcine ovary follicle atresia. Mol Cells 38, 304–311, https://doi.org/10.14348/molcells.2015.2122 (2015).
Oktay, K., Buyuk, E., Oktem, O., Oktay, M. & Giancotti, F. G. The c-Jun N-terminal kinase JNK functions upstream of Aurora B to promote entry into mitosis. Cell Cycle 7, 533–541, https://doi.org/10.4161/cc.7.4.5660 (2008).
Maekawa, R. et al. Changes in gene expression of histone modification enzymes in rat granulosa cells undergoing luteinization during ovulation. Journal of ovarian research 9, 15, https://doi.org/10.1186/s13048-016-0225-z (2016).
Berisha, B., Schams, D., Rodler, D., Sinowatz, F. & Pfaffl, M. W. Expression and localization of members of the thrombospondin family during final follicle maturation and corpus luteum formation and function in the bovine ovary. J Reprod Dev 62, 501–510, https://doi.org/10.1262/jrd.2016-056 (2016).
Li, R. et al. MiR-34a suppresses ovarian cancer proliferation and motility by targeting AXL. Tumour Biol 36, 7277–7283, https://doi.org/10.1007/s13277-015-3445-8 (2015).
Tu, F. et al. miR-34a targets the inhibin beta B gene, promoting granulosa cell apoptosis in the porcine ovary. Genet Mol Res 13, 2504–2512, https://doi.org/10.4238/2014.January.14.6 (2014).
Myers, M. & Pangas, S. A. Regulatory roles of transforming growth factor beta family members in folliculogenesis. Wiley Interdiscip Rev Syst Biol Med 2, 117–125, https://doi.org/10.1002/wsbm.21 (2010).
Elizur, S. E. et al. Factors predicting IVF treatment outcome: a multivariate analysis of 5310 cycles. Reprod Biomed Online 10, 645–649 (2005).
Hourvitz, A. et al. Role of embryo quality in predicting early pregnancy loss following assisted reproductive technology. Reprod Biomed Online 13, 504–509 (2006).
Fadini, R. et al. Effect of different gonadotrophin priming on IVM of oocytes from women with normal ovaries: a prospective randomized study. Reprod Biomed Online 19, 343–351 (2009).
Yuan, J. S., Wang, D. & Stewart, C. N. Jr. Statistical methods for efficiency adjusted real-time PCR quantification. Biotechnology journal 3, 112–123, https://doi.org/10.1002/biot.200700169 (2008).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140, https://doi.org/10.1093/bioinformatics/btp616 (2010).
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43, e47, https://doi.org/10.1093/nar/gkv007 (2015).
Zhou, J. et al. The let-7g microRNA promotes follicular granulosa cell apoptosis by targeting transforming growth factor-beta type 1 receptor. Mol Cell Endocrinol 409, 103–112, https://doi.org/10.1016/j.mce.2015.03.012 (2015).
Gebremedhn, S. et al. MicroRNA Expression Profile in Bovine Granulosa Cells of Preovulatory Dominant and Subordinate Follicles during the Late Follicular Phase of the Estrous Cycle. PLoS One 10, e0125912, https://doi.org/10.1371/journal.pone.0125912 (2015).
Sirotkin, A. V., Ovcharenko, D., Grossmann, R., Laukova, M. & Mlyncek, M. Identification of microRNAs controlling human ovarian cell steroidogenesis via a genome-scale screen. J Cell Physiol 219, 415–420, https://doi.org/10.1002/jcp.21689 (2009).
Navakanitworakul, R. et al. Characterization and Small RNA Content of Extracellular Vesicles in Follicular Fluid of Developing Bovine Antral Follicles. Scientific reports 6, 25486, https://doi.org/10.1038/srep25486 (2016).
Wang, C. et al. MicroRNA-125a-5p induces mouse granulosa cell apoptosis by targeting signal transducer and activator of transcription 3. Menopause 23, 100–107, https://doi.org/10.1097/GME.0000000000000507 (2016).
Sang, Q. et al. Identification of microRNAs in human follicular fluid: characterization of microRNAs that govern steroidogenesis in vitro and are associated with polycystic ovary syndrome in vivo. The Journal of clinical endocrinology and metabolism 98, 3068–3079, https://doi.org/10.1210/jc.2013-1715 (2013).
Zhang, Q. et al. MicroRNA-181a suppresses mouse granulosa cell proliferation by targeting activin receptor IIA. PLoS One 8, e59667, https://doi.org/10.1371/journal.pone.0059667 (2013).
Salilew-Wondim, D. et al. The expression pattern of microRNAs in granulosa cells of subordinate and dominant follicles during the early luteal phase of the bovine estrous cycle. PLoS One 9, e106795, https://doi.org/10.1371/journal.pone.0106795 (2014).
Xie, S., Batnasan, E., Zhang, Q. & Li, Y. MicroRNA Expression is Altered in Granulosa Cells of Ovarian Hyperresponders. Reproductive Sciences 23, 1001–1010, https://doi.org/10.1177/1933719115625849 (2016).
Xu, B., Zhang, Y. W., Tong, X. H. & Liu, Y. S. Characterization of microRNA profile in human cumulus granulosa cells: Identification of microRNAs that regulate Notch signaling and are associated with PCOS. Mol Cell Endocrinol 404, 26–36, https://doi.org/10.1016/j.mce.2015.01.030 (2015).
Hossain, M. M. et al. Identification and characterization of miRNAs expressed in the bovine ovary. BMC Genomics 10, 443, https://doi.org/10.1186/1471-2164-10-443 (2009).
Yao, N. et al. Follicle-stimulating hormone regulation of microRNA expression on progesterone production in cultured rat granulosa cells. Endocrine 38, 158–166, https://doi.org/10.1007/s12020-010-9345-1 (2010).
Sontakke, S. D., Mohammed, B. T., McNeilly, A. S. & Donadeu, F. X. Characterization of microRNAs differentially expressed during bovine follicle development. Reproduction (Cambridge, England) 148, 271–283, https://doi.org/10.1530/REP-14-0140 (2014).
Xu, S., Linher-Melville, K., Yang, B. B., Wu, D. & Li, J. Micro-RNA378 (miR-378) regulates ovarian estradiol production by targeting aromatase. Endocrinology 152, 3941–3951, https://doi.org/10.1210/en.2011-1147 (2011).
Lei, L., Jin, S., Gonzalez, G., Behringer, R. R. & Woodruff, T. K. The regulatory role of Dicer in folliculogenesis in mice. Mol Cell Endocrinol 315, 63–73, https://doi.org/10.1016/j.mce.2009.09.021 (2010).
Sathyapalan, T., David, R., Gooderham, N. J. & Atkin, S. L. Increased expression of circulating miRNA-93 in women with polycystic ovary syndrome may represent a novel, non-invasive biomarker for diagnosis. Scientific reports 5, 16890, https://doi.org/10.1038/srep16890 (2015).
Jiang, L. et al. MicroRNA-93 promotes ovarian granulosa cells proliferation through targeting CDKN1A in polycystic ovarian syndrome. The Journal of clinical endocrinology and metabolism 100, E729–738, https://doi.org/10.1210/jc.2014-3827 (2015).
Kim, Y. J. et al. MicroRNA Profile of Granulosa Cells after Ovarian Stimulation Differs According to Maturity of Retrieved Oocytes. Geburtshilfe Frauenheilkd 76, 704–708, https://doi.org/10.1055/s-0041-111173 (2016).
Diez-Fraile, A. et al. Age-associated differential microRNA levels in human follicular fluid reveal pathways potentially determining fertility and success of in vitro fertilization. Hum Fertil (Camb) 17, 90–98, https://doi.org/10.3109/14647273.2014.897006 (2014).
Acknowledgements
This study was supported by grants from the Ministry of Health, Israel (3-00000-7410) and from the Chaim Sheba Medical Center (R&D Sheba grant 5/2014).
Author information
Authors and Affiliations
Contributions
A.H., G.M.Y., E.M. and M.D. designed the study, performed the analysis and interpretation of the data, prepared figures and tables, and wrote and finalized the manuscript. A.H., G.M.Y., M.D. and L.S. prepared figures and tables. Y.Y., M.B., R.M. and L.O. contributed to the interpretation of the data and manuscript preparation. L.O., Y.Y. and M.D.C. were in charge of tissue and oocyte processing. A.H., G.C., R.F., M.M.R. and M.D.C. were involved in patient recruitment. All of the authors contributed to data analysis and reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yerushalmi, G.M., Salmon-Divon, M., Ophir, L. et al. Characterization of the miRNA regulators of the human ovulatory cascade. Sci Rep 8, 15605 (2018). https://doi.org/10.1038/s41598-018-33807-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-33807-y
Keywords
This article is cited by
-
Epistasis Detection via the Joint Cumulant
Statistics in Biosciences (2022)
-
New biomarkers in peripheral blood of patients with ovarian cancer: high expression levels of miR-16-5p, miR-17-5p, and miR-638
Archives of Gynecology and Obstetrics (2022)
-
Transcriptomics of cumulus cells – a window into oocyte maturation in humans
Journal of Ovarian Research (2020)
-
Target prediction and validation of microRNAs expressed from FSHR and aromatase genes in human ovarian granulosa cells
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.