- Open Access
Transcriptomics of cumulus cells – a window into oocyte maturation in humans
Journal of Ovarian Research volume 13, Article number: 93 (2020)
Cumulus cells (CC) encapsulate growing oocytes and support their growth and development. Transcriptomic signatures of CC have the potential to serve as valuable non-invasive biomarkers for oocyte competency and potential. The present sibling cumulus-oocyte-complex (COC) cohort study aimed at defining functional variations between oocytes of different maturity exposed to the same stimulation conditions, by assessing the transcriptomic signatures of their corresponding CC. CC were collected from 18 patients with both germinal vesicle and metaphase II oocytes from the same cycle to keep the biological variability between samples to a minimum. RNA sequencing, differential expression, pathway analysis, and leading-edge were performed to highlight functional differences between CC encapsulating oocytes of different maturity.
Transcriptomic signatures representing CC encapsulating oocytes of different maturity clustered separately on principal component analysis with 1818 genes differentially expressed. CCs encapsulating mature oocytes were more transcriptionally synchronized when compared with CCs encapsulating immature oocytes. Moreover, the transcriptional activity was lower, albeit not absent, in CC encapsulating mature oocytes, with 2407 fewer transcripts detected than in CC encapsulating immature (germinal vesicle - GV) oocytes. Hallmark pathways and ovarian processes that were affected by oocyte maturity included cell cycle regulation, steroid metabolism, apoptosis, extracellular matrix remodeling, and inflammation.
Herein we review our findings and discuss how they align with previous literature addressing transcriptomic signatures of oocyte maturation. Our findings support the available literature and enhance it with several genes and pathways, which have not been previously implicated in promoting human oocyte maturation. This study lays the ground for future functional studies that can enhance our understanding of human oocyte maturation.
Cumulus cells (CC) provide somatic support to the maturing oocyte, and together they comprise the functional unit known as Cumulus-Oocyte-Complex (COC) [1,2,3,4]. Understanding oocyte maturation and associated pathologies can help improve fertility treatments, as well as culture media conditions in the lab. Cumulus cells can be collected without compromising the oocyte, and their transcriptomic signatures are valuable non-invasive biomarkers for processes within the oocyte [5, 6].
Oocyte maturation is contingent on rapid transcription and translation, governed by paracrine and autocrine signaling prior to ovulation [7,8,9]. Once matured, the MII oocyte is less transcriptionally active than its precursors, relying on stored mRNA transcripts that were acquired throughout its maturation, to undergo successful fertilization and support early embryo development until embryonic genome activation [10,11,12]. Moreover, in addition to the stored transcripts, there is active transportation of transcripts from cumulus cells to the growing oocyte through trans-zonal projections [13, 14]. These projections are critical for both oocyte and cumulus cell differentiation .
Currently, oocyte assessment relies mainly on morphological criteria that provide little insight on oocyte quality and competence . Furthermore, available techniques for maturing oocytes in-vitro are inefficient and do not provide good alternatives in cases were pathologies of oocyte maturation lead to retrieval of multiple immature eggs despite adequate controlled ovarian stimulation. This is why molecular investigation of processes in COCs responsible for nuclear and cytoplasmic maturation, as well as pathologies that could arise, are key in improving patient treatment and outcomes. The objectives of this study were: 1) to profile the transcriptome of CC from mature MII and immature GV oocytes, from the same treatment cycle, 2) to use the above profile to validate existing transcriptomic literature exploring human oocyte maturation, and 3) to provide a comprehensive list of genes impacted by in vitro maturation which can be used for future studies exploring human oocyte maturation.
Collected samples and patient demographics
A total of 22 cumulus-oocyte complexes (COCs) were collected from 11 patients, for RNA sequencing, with a mean age of 33.9 ± 1.6 years old, mean BMI of 25.4 ± 2.3 kg/m2, mean Anti-Mullerian Hormone (AMH) levels of 29.8 ± 12.5 pmol/L, mean Day 3 FSH levels of 7.9 ± 1.4 mIU/mL, and the mean number of collected oocytes was 10.6 ± 1.9 (Table 1). All samples had sufficient number of sequencing reads, high average quality scores, and high sequence alignment rates sufficient for differential gene expression analysis, as per guidelines previously published for quality control of RNAseq experiments  (Supplemental Table S2).
Samples clustered according to the degree of maturity of the encapsulated oocyte
A total of 6220 genes were detected in the MII-CC cohort, 202 of which were unique, and 8627 genes were detected in the GV-CC cohort, 2609 of which were unique (Fig. 1a). Unsupervised hierarchical clustering demonstrated a separation of the GV-CC cohort from the MII-CC cohort. Notably, the MII-CC cohort clustered more tightly than the GV-CC cohort, indicating that with maturation there is decreased inter-sample variability (Fig. 1b). This was further demonstrated by principal component analysis (PCA) (Fig. 1c) with 16.3% of variability in the dataset corresponding to oocyte maturity (PC1).
Differential expression reveals marked differences in gene expression according to degree of maturity of the encapsulated oocyte
In order to ensure the sample size and sequencing depth were sufficient to capture biologically significant differences between the MII and GV CCs, a post-hoc power analysis was conducted showing that the minimum required biological replicates is 6 and the minimum sequencing depth is 10 million reads per replicate, both of which were exceeded in this study . When comparing MII-CC with GV-CC cohorts, 1818 genes were differentially expressed (Supplemental Table S3), which comprise 10.3% of annotated RefSeq genes (2 < FC < − 2 and FDR < 0.05). Of these, 40 genes changed by 10-fold or more, 207 genes changed by 5 to 10-fold, and the rest (1571 genes) changed by 2 to 5-fold. When accounting for the direction of differential expression, 1031 genes were significantly downregulated (5.9% of annotated genes) and 787 genes were significantly upregulated (4.5% of annotated genes) in the MII-CC cohort compared with GV-CC (Fig. 1d). The top 20 genes enriched in the different maturation stages are reported in Table 2, all differentially expressed genes are reported in Supplemental Table S3.
Novel findings from this study enhance available literature exploring processes that lead to synchronized oocyte maturity
When comparing our differentially expressed genes to those previously reported in the literature, 42 genes were associated with oocyte maturation in both this study and in previous literature (Table 3). Forty-five genes previously correlated with oocyte maturation were not differentially expressed in the current study (Supplemental Table S4) . Three thousand five hundred and fifty-four genes were differentially expressed in our study and have not been annotated previously in studies exploring oocyte maturation. Of these, 129 are known to be regulated by at least one gonadotropin, 82 have been previously reported to be regulated by LH alone, 16 by FSH alone and 31 by both LH and FSH (Supplemental Table S5).
Pathway analysis and leading-edge analysis revealed significant differences between maturational stages
Gene set enrichment analysis (GSEA) was conducted to determine the pathways and cellular processes that are altered throughout oocyte maturity, thus allowing for the interpretation of the complex interactions between differentially expressed genes (Fig. 2) . Furthermore, leading edge analysis (LEA) helped determine which genes were driving pathway enrichment scores, and to better understand the major biological differences between the two cohorts.
Pathways and processes that were primarily enriched in downregulated genes in the MII-CC cohort included
Nuclear maturation, chromatin remodeling and replication initiation, faithful chromosome segregation, apoptosis and inflammation (Fig. 2). Furthermore, specific genes identified as biologically significant genes for oocyte maturation including nuclear maturation (Fig. 3a), chromatin remodeling and DNA replication initiation (Fig. 3b), and apoptosis and inflammation (Fig. 3c) are further highlighted in Fig. 3.
Pathways and processes that were primarily enriched in upregulated genes in the MII-CC cohort included
Extracellular matrix (ECM) components and its remodeling enzymes, and steroid metabolism and processing (Fig. 2). Genes identified as biologically significant genes for oocyte maturation including ECM remodeling (Fig. 3d), and steroid metabolism (Fig. 3e) are further highlighted in Fig. 3.
Leading edge analysis identified several genes involved in cell cycle control (CDK1, CCNB1, CCNB2, CCNA2, BUB1, and CDC20), DNA replication initiation (MCM2–7), and centromere assembly and organization (CENPF), among others, as the genes that were most significantly driving the gene set enrichment analysis (Supplemental Table S6). GSEA on DE genes known to be regulated by at least one gonadotropin revealed two major pathways; transcriptional regulation of tp53 (apoptosis), overall enriched in downregulated genes, and metabolic biosynthesis, overall enriched in upregulated genes. LEA identified several genes involved in cell cycle control and cell death indicating that the MII-CC cohort, in response to LH and FSH decreases cell death signaling and increases biosynthesis.
Validation of NGS results by qPCR
We selected 16 differentially expressed genes as determined by RNAseq and are known to be involved in various pathways of CC expansion and oocyte maturation. For all selected genes, similar fold changes were observed using qPCR as were observed using NGS (Fig. 4a). In addition, when assessing multiple CCs from 7 additional patients at both the MII and GV stage (from 2/3 GV-CC and 2/3 MII-CC per patient) (a total of 34 COCs), the expression of all tested genes was consistent within each patient. Interestingly, the expression of all genes used for validation across all MII-CC were consistent both within and between patients (∆Ct SEM (range) = 0.26(0.12–0.56)). GV-CC were also consistent within and between patients, however to a lesser extent (∆Ct SEM (range) = 0.40(0.21–0.57)) (Fig. 4b).
This was further validated using Pearson’s R correlation which also demonstrated that the cohort with the highest similarity was MII-CC within the same patient, with a R2 of 0.954 ± 0.023 followed by GV-CC within the same patient, with a R2 of 0.906 ± 0.036. The correlation coefficients and significance are outlined in Table 4. Taken together, this indicates that the heterogeneity of the population of CCs within the ovary is low, and sampling one mature MII-CC and one immature GV-CC for NGS is representative of the cohort of CCs.
This study is novel in the choice of cohorts for comparison. We included CC encapsulating oocytes arrested at the GV stage, despite being exposed to adequate stimulation, and compared them with CC encapsulating MII oocytes that matured in vivo, from the same patients during an IVF cycle. This comparison allowed us to fine-tune our understanding of oocyte maturation in-vivo. Previous human oocyte maturation studies analyzed COCs from in-vitro maturation cycles [16, 29,30,31,32,33]. This is why their findings are more relevant for processes in in-vitro maturation per se, but they do not tease out factors specifically associated with failed maturation despite adequate in-vivo exposure to controlled ovarian hyper stimulation (COH).
Forty-two differentially expressed genes in our study, have been previously associated with oocyte maturation and cumulus cell expansion in IVF treatments (Table 3). Following our extensive literature search, we identified 45 genes, which were previously correlated with human oocyte maturation, but not captured by our study design, possibly because they may impact IVM alone and may not reflect in vivo maturation (Supplemental Table S4). A third group of 3554 genes captured in the current study but not in previous studies which represents a novel group of genes that should be further explored as they have not been previously implicated in human oocyte maturation (Supplemental Table S3). Of these, 129 genes have been previously shown to be regulated by either LH, FSH, or both (Supplemental Table S5).
In this study, several factors and their regulators involved in nuclear maturation and cell cycle control were differentially expressed between cumulus cells encapsulating oocytes of different maturity, reiterating findings from previous studies [19, 20, 29, 34]. These include cell cycle regulators (BIRC5, BUB1, BUB1B, CCNA2, CCNB, CDK1, FBXO5 MAD2L1, and PTTG1) and components of the centromere (CENPA, CENPE, and CENPH) . In our MII-CC cohort we observed downregulation of MCM2–7, which form the hexameric pre-replication protein complex. This complex is involved in initiating replication forks and recruiting other DNA replication related proteins. We also observed downregulation of TOP2A, which relaxes supercoiled and circular DNA molecules. Reinforcing available literature that states that while crucial at the MI stage for chromatin remodeling [21, 22], its activity decreases in mature oocytes .
Apoptosis was also attenuated in the MII-CC cohort, further supporting decreased cell turnover with advanced maturity. Related pathways including Wnt pathway and Akt-pathway were affected, as demonstrated by downregulation of SFRP4, a potent inhibitor of Wnt signaling , and upregulation of OSMR, an activator of Akt-mediated proliferation . These findings corroborate previous literature reporting downregulation of SFRP4 during oocyte maturation [26, 35, 36], and upregulation of OSMR in bovine preovulatory follicles post-triggering by gonadotropins .
Extracellular matrix remodeling was also altered between the two maturity cohorts, as evident by members of the matrix metalloproteinases (MMP) family and their inducers (MMP11 and SPARC1L). Again, this supports previous literature showing significant decrease of MMP11 in granulosa cells following hCG administration . This effect is further demonstrated by increased expression of TNC, NID2, and SPOCK2 - all ECM proteins and MMP substrates [26, 39,40,41]. Notably, well characterized ECM remodeling enzymes, ADAMTS1 and SERPINE2, were also differentially expressed, aligning with previous studies [42, 43]. Both play critical roles in follicular remodeling during follicular growth and rupture , by metabolizing Versican and Hyaluronan which lead to cumulus cell matrix expansion and attenuation .
Another key process enhanced in follicular niche maturation is inflammation, which is crucial for ovulation. Upon gonadotropin stimulation, the follicle wall is weakened, thereby facilitating its eventual rupture . In our MII-CC cohort, we observed marked upregulation of genes associated with inflammation, including members of the Interleukin and TGF-beta families. Among the genes upregulated in our MII-CC cohort were IL18R1 which promotes cumulus cell expansion , and TGFBR3 which promotes cellular differentiation, migration, adhesion and extracellular matrix production [48, 49]. IL6ST which is part of the cytokine receptor complex (gp130) was also upregulated in the MII-CC cohort, consistent with previous studies in non-human primates and equine models [50, 51].
Key players that emerged in our cohort as being significant for cumulus cells to facilitate oocyte maturation are AREG, EREG, PTGS2, and STAR. Two factors at the heart of this complex process are AREG and EREG, which have been shown to mediate the LH signal driving cumulus expansion and oocyte maturation [19, 32, 52]. They also activate the EGF receptor (EGFR) which in turn releases matrix metalloproteinases (MMPs) and promotes cumulus expansion [52, 53]. Furthermore, in conjunction with progesterone, AREG and EREG enhance PTGS2 (also upregulated in our MII-CC cohort) via EGF to increase prostaglandin production and maintenance of chromosomal spindles [33, 54,55,56]. In addition, AREG mediates hCG-induced STAR expression (also upregulated in our MII-CC cohort), which plays a key role in steroid and progesterone production in human granulosa cells , and is a potential predictive biomarker for nuclear maturation  and oocyte quality . It is important to note, that despite being well defined as key in ovarian maturation [32, 52, 59], EREG has not been found to be differentially expressed in previous genomic signature studies addressing this question. This further highlights the importance of our study design in better refining the pathophysiology of oocyte maturation.
Other critical factors confirmed by this study are IL1, FSHR, BDNF, HSD11B1, and HSD17B1, all of which are implicated in the control of steroid synthesis and paracrine response to steroids.
IL1 (both alpha and beta subunits), which stimulates steroidogenesis, was upregulated in the MII-CC cohort with a concurrent decreased expression of FSHR in the same cohort, substantiating what was previously observed in rodents and humans [60, 61]. BDNF, which modulates granulosa cell function via FSHR-coupled signaling pathway, to affect aromatase-mediated steroidogenesis, was also downregulated in our MII-CC cohort .
HSD11B1, the enzyme responsible for cortisone production, an essential substrate for steroid hormone synthesis, was upregulated in our MII-CC cohort. A companion enzyme, HSD17B1, catalyzes the last step in estrogen metabolism converting E1 of low estrogenic activity to E2 of high activity using cortisone as a substrate . HSD17B1 has not been captured in previous human studies, but was downregulated in our MII-CC cohort, consistent with the results seen in a previous bovine study , and further highlighting the advantage of our study design.
Both LHCGR and FSHR were differentially expressed in the MII-CC cohort when compared to the GV-CC cohort. To further explore how this may impact the transcriptional profiles, we identified all genes differentially expressed between these two cohorts that are known to be regulated by LH and/or FSH. GSEA based on these 129 genes, identified two major pathways: regulation of apoptosis (via p53) and biosynthesis of various biomolecules.
Overall, apoptosis was enriched in downregulated genes. Interestingly, several major players in the regulation of apoptosis, including BIRC5, TP53, HMGB1, HMGB2, and SFRP4 are also known to be regulated by LH and/or FSH [24, 35, 65,66,67]. Taken together, these findings suggest that the inability of the GV-CC cohort to appropriately respond to FSH/LH may, in turn, lead to failure of the CCs to effectively dampen the apoptotic signals, causing the COC to enter a stage of maturation arrest and follicular atresia.
Overall, biosynthesis was enriched in upregulated genes among the MII-CC cohort. Notably, several members of the CYP family, which were upregulated, and are involved in the biosynthesis of estrogen and androgens, are known to be regulated by LH and/or FSH [68,69,70]. Taken together, these findings suggest that the MII-CCs are responding adequately to gonadotropin administration.
Collectively, it appears that the GV-CC cohort, failed to adequately synthesize estrogen, despite exposure to gonadotropins and, thus, began to upregulate genes involved in apoptosis. This may be due to insufficient LH and/or FSH receptors on these COCs or due to another underlying malady.
Finally, we show that PDE3A, known to improve nuclear-cytoplasmic synchrony , is significantly upregulated in our MII-CC cohort. While this gene has not been studied in cumulus cells in the context of oocyte maturation in humans, it has been shown that an increase in oocyte PDE3A activity causes delayed spontaneous meiotic maturation, coupled with extended gap junctional communication between the CC and the oocyte. Such a delay has a positive effect on oocyte cytoplasmic maturation, thereby improving oocyte developmental potential . The fact that upregulation of this gene was captured by our study design speaks once again to the strength of our study and to what it adds to current literature.
Methodological strengths of this study include (i) a sibling COC design allowing to minimize the biologic variability between cohorts, (ii) exploring transcriptomic dynamics in cumulus cells, which are considered valuable non-invasive markers for oocyte quality [73,74,75], and (iii) performing next generation sequencing (NGS), which is the most unbiased approach currently available for exploring transcriptomic signatures. A methodological weakness of this study is our inability to compare our findings with a third cohort of CCs encapsulating naïve GVs from the same stimulation cycle. Furthermore, to keep our sequencing sample size small, only 2 cumulus masses were selected from each patient; 1 to represent the immature GV COCs, and the other to represent the mature MII COCs. To lower the risk for selection bias due to this study design, we chose to perform validation studies that established the tight correlations within each cohort. Lastly, our small sample size and study design did not allow for tracking outcomes on an individual oocyte basis.
In conclusion, our findings enhance current literature on oocyte maturation by identifying CC genes not previously associated with oocyte maturation that may be involved in this process. Our novel list of genes can serve as a springboard for future studies. Our future studies will focus on determining the functional significance of these findings and on attempting to identify how different treatment options may favor a more synchronized mature/competent state. In addition, to further validate genes that are critical for oocyte maturation and competency, further large-scale studies correlating gene expression with clinical outcomes using a targeted transcriptome panel are needed.
Materials and methods
Patient recruitment, data collection and cumulus cell isolation
Cumulus cell samples were collected from eighteen patients undergoing IVF-ICSI cycles at the CReATe Fertility Centre (Toronto, ON, Canada), between August 2016 and June 2017. Exclusion criteria were patients diagnosed with PCOS, as per Rotterdam criteria, as well as patients with endometriosis diagnosed by laparoscopy. Samples from eleven patients were used for RNA-seq (22 COCs, 11 mature (MII), and 11 immature (GV)), and samples from seven additional patients were used for qPCR validation of the findings (a total of 18 mature (MII) and 16 immature (GV) COCs). Patients were treated using a standard antagonist protocol, with initial gonadotropin dosing and subsequent adjustments at the discretion of the treating physician.
Ultrasound guided oocyte retrieval was performed 35-36 h post hCG injection. COCs were identified under a stereomicroscope and only COCs completely and tightly enclosed by compact CCs were used for this study to minimize the potential collection of contaminating granulosa cells. Selected COCs were serially washed three times in Quinn’s Advantage Medium (Sage, USA) to remove cellular contaminants, and to further reduce the possibility of granulosa cell contamination. CCs were then mechanically separated from each oocyte individually in Quinn’s Advantage Medium (Sage, USA), under paraffin oil by one experienced embryologist within 1 h of oocyte retrieval. The oocytes corresponding to individually collected CC were separately exposed to hyaluronidase (80 IU/ml) immediately after mechanical separation of CCs, washed in Quinn’s Advantage Medium (Sage, USA). Maturational stage was assessed through the observation of the nucleus of the oocyte. Oocytes with an extruded polar body were deemed mature (MII), oocytes with an intact germinal vesicle was deemed immature (GV), oocytes without an observable germinal vesicle or an extruded polar body were deemed MI and excluded from further analysis. The CC were collected from single MII (n = 29) and GV oocytes (n = 27), frozen separately in 300ul of RNA lysis buffer RL (Norgen Biotek, Canada), and stored at − 80 °C until RNA extraction. Clinical data including patient demographics, medical history, and ovarian stimulation related parameters, were collected for all patients enrolled in this study.
RNA extraction and NGS library preparation
RNA extraction, cDNA conversion and NGS library construction and normalization were conducted as previously described . Briefly, for NGS one MII-CC and one GV-CC were chosen at random from each of the 11 patients for RNA extraction using the Total RNA Purification Kit Micro (Norgen Biotek, Canada). The quantity was assessed using Qubit RNA HS (ThermoFisher, USA) and RNA integrity assessed using 2100 Bioanalyser RNA 6000 Pico Total RNA Kit (Agilent Technologies, Canada). cDNA was synthesized using the SMART-seq v4 Ultra Low Input RNA Kit (Takara Bio Inc., Japan) according to the sample preparation guide and using 14 rounds of amplification. Sequencing libraries were constructed using Nextera XT (Illumina, USA) and 1 ng of amplified cDNA according to the sample preparation guide. Final sequencing libraries were assessed for quantity and quality using the KAPA Library Quantification Kit (KAPA Biosystems, Switzerland) and 2100 Bioanalyser High Sensitivity DNA Kit (Agilent Technologies), respectively. Normalized libraries were pooled, denatured, diluted to 1.4 pmol/l and loaded onto a High Output (300 cycle) flow cell (Illumina, USA) followed by sequencing (2 × 127 bp) on a NextSeq 550 (Illumina, USA).
The recommendations outlined by Ching et al. 2014 were followed when selecting the differential expression package, as well as using a paired-sample RNA-Seq design as suggested . FASTQ files were generated using bcl2fastq2 (v2.17) and the read quality was assessed. Sequences were trimmed based on quality (Phred > 28). Raw trimmed reads were aligned to Human Genome Assembly 38 (hg38) using STAR (v2.5.3a)  and quantified to RefSeq (Release 84). Low expressed transcripts were excluded (maximum counts < 10) and differential expression (DE) was conducted on the remaining counts using DESeq2 (v3.5) . Principal Component Analysis (PCA) and hierarchical clustering (HC) were conducted to assess the relationship between samples and determine covariates contributing to variation in the dataset. Principal component 1 (PC1) accounts for the largest proportion of the variability observed within the dataset, PC2 accounts for the second largest, and so on. Differentially expressed genes were identified by comparing all mature CC samples (MII-CC) to all immature CC samples (GV-CC), and were deemed to be differentially expressed if the gene had a Fold change (FC) of more than absolute value of 2, and a false discovery rate (FDR) < 0.05. This analysis was conducted in Partek Flow (version 8.0.19.0408).
Gene Set Enrichment Analysis (GSEA) was performed to determine the effect all differentially expressed genes have on cellular processes and functions . The resulting pathway list was cross referenced with a custom gene set created and supported by the Bader Lab (University of Toronto) which is comprised of all GO database, KEGG, and Reactome gene sets (v2018-12-01) (http://download.baderlab.org/EM_Genesets/) . Genes that could not be mapped to any gene-set term were excluded from the comparison. Gene sets with 10 or fewer genes and/or a q-value > 0.05 were excluded from further analysis. Following GSEA, leading edge analysis (LEA) was conducted to determine what genes were driving the gene set enrichment score, as well as to highlight genes that were shared between gene sets.
To further explore the impact FSH and/or LH may have on the transcriptome, we identified all differentially expressed genes that are known to be regulated by LH, FSH or both  and performed GSEA and LEA as described previously.
NGS validation by qPCR
Sixteen genes (and one reference gene) were chosen for validation from the list of differentially expressed genes. The choice of genes was based on previous annotations deeming these genes as biologically significant, identification by leading edge analysis, and/or participation in key ovarian gene pathways. Pre-designed and validated PrimeTime™ qPCR assays (IDT, USA) were used for validation of NGS results with RPLP0 as the reference gene. All targets were assayed in duplicate using PrimeTime™ Gene Expression MasterMix (IDT, USA) (polymerase activation at 95 °C for 3 min; 45 cycles of 15 s denaturation at 95 °C and 1 min annealing/extension at 60 °C). Relative fold change (ΔΔCt) was employed to quantify gene expression . Data analysis was performed using GraphPad Prism (version 5.02). The list of primers and probes used for validation are given in Supplemental Table S1. qPCR for the above validated genes was also performed on multiple COCs of the same maturational stages on samples from 7 additional patients. This was carried out to ensure the validity of a random choice of a single COC as a representative of all COCs at the same maturational stage from the same patient. Similarities between individual CCs was measured by creating an overall Pearson correlation between the vectors of variables (16 targets genes) using SPSS Proximities. This created an overall Pearson R2. Paired t-test was used to compare correlations within GV-CC and within MII-CC groups, separately. A non-paired test was used to compare overall GV-CC and overall MII-CC correlations due to an unequal number of observations for some individuals.
Gene annotation and literature search
To determine the clinical significance of our bioinformatic findings, differentially expressed genes were cross referenced with available datasets in the literature by searching the PubMed database for previous studies assessing the transcriptome of human CC using NGS, Microarray, or qPCR. Differentially expressed genes were further reviewed in depth using the Ovarian Kaleidoscope Database  and GeneCards Human Gene databases (http://www.genecards.org/), to correlate our bioinformatic findings with hallmark physiological and pathological processes in the ovary.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article (Supplementary Table S7).
Eppig JJ. Oocyte control of ovarian follicular development and function in mammals. Reproduction. 2001;122(6):829–38.
Ackert CL, Gittens JE, O'Brien MJ, Eppig JJ, Kidder GM. Intercellular communication via connexin43 gap junctions is required for ovarian folliculogenesis in the mouse. Dev Biol. 2001;233(2):258–70.
Matzuk MM, Burns KH, Viveiros MM, Eppig JJ. Intercellular communication in the mammalian ovary: oocytes carry the conversation. Science. 2002;296(5576):2178–80.
Gilchrist RB, Lane M, Thompson JG. Oocyte-secreted factors: regulators of cumulus cell function and oocyte quality. Hum Reprod Update. 2008;14(2):159–77.
McKenzie LJ, Pangas SA, Carson SA, Kovanci E, Cisneros P, Buster JE, et al. Human cumulus granulosa cell gene expression: a predictor of fertilization and embryo selection in women undergoing IVF. Hum Reprod. 2004;19(12):2869–74.
van Montfoort APA, Geraedts JPM, Dumoulin JCM, Stassen APM, Evers JLH, Ayoubi TAY. Differential gene expression in cumulus cells as a prognostic indicator of embryo viability: a microarray analysis. Mol Hum Reprod. 2008;14(3):157–68.
Su YQ, Sugiura K, Woo Y, Wigglesworth K, Kamdar S, Affourtit J, et al. Selective degradation of transcripts during meiotic maturation of mouse oocytes. Dev Biol. 2007;302(1):104–17.
Eichenlaub-Ritter U, Vogt E, Yin H, Gosden R. Spindles, mitochondria and redox potential in ageing oocytes. Reprod BioMed Online. 2004;8(1):45–58.
Eppig JJ. Reproduction: oocytes call, Granulosa Cells Connect. Curr Biol. 2018;28(8):R354–R6.
Bachvarova R. Gene expression during oogenesis and oocyte development in mammals. Dev Biol (N Y 1985). 1985;1:453–524.
Braude P, Bolton V, Moore S. Human gene expression first occurs between the four- and eight-cell stages of preimplantation development. Nature. 1988;332(6163):459–61.
De La Fuente R, Viveiros MM, Burns KH, Adashi EY, Matzuk MM, Eppig JJ. Major chromatin remodeling in the germinal vesicle (GV) of mammalian oocytes is dispensable for global transcriptional silencing but required for centromeric heterochromatin function. Dev Biol. 2004;275(2):447–58.
Macaulay AD, Gilbert I, Scantland S, Fournier E, Ashkar F, Bastien A, et al. Cumulus Cell Transcripts Transit to the Bovine Oocyte in Preparation for Maturation. Biol Reprod. 2016;94(1):16.
Macaulay AD, Gilbert I, Caballero J, Barreto R, Fournier E, Tossou P, et al. The gametic synapse: RNA transfer to the bovine oocyte. Biol Reprod. 2014;91(4):90.
Baena V, Terasaki M. Three-dimensional organization of transzonal projections and other cytoplasmic extensions in the mouse ovarian follicle. Sci Rep. 2019;9(1):1262.
Wells D, Patrizio P. Gene expression profiling of human oocytes at different maturational stages and after in vitro maturation. Am J Obstet Gynecol. 2008;198(4):455 e1–9 discussion e9–11.
Ching T, Huang S, Garmire LX. Power analysis and sample size estimation for RNA-Seq differential expression. RNA (New York). 2014;20(11):1684–96.
Busby MA, Stewart C, Miller CA, Grzeda KR, Marth GT. Scotty: a web tool for designing RNA-Seq experiments to measure differential gene expression. Bioinformatics. 2013;29(5):656–7.
Ouandaogo ZG, Frydman N, Hesters L, Assou S, Haouzi D, Dechaud H, 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. Hum Reprod. 2012;27(8):2438–47.
Ouandaogo ZG, Haouzi D, Assou S, Dechaud H, Kadoch IJ, De Vos J, et al. Human cumulus cells molecular signature in relation to oocyte nuclear maturity stage. PloS one. 2011;6(11):e27179-e.
Li X-M, Yu C, Wang Z-W, Zhang Y-L, Liu X-M, Zhou D, et al. DNA Topoisomerase II Is Dispensable for Oocyte Meiotic Resumption but Is Essential for Meiotic Chromosome Condensation and Separation in Mice1. Biology of Reproduction. 2013;89(5):1-11.
Assou S, Cerecedo D, Tondeur S, Pantesco V, Hovatta O, Klein B, et al. A gene expression signature shared by human mature oocytes and embryonic stem cells. BMC Genomics. 2009;10:10.
Yerushalmi GM, Salmon-Divon M, Ophir L, Yung Y, Baum M, Coticchio G, et al. Characterization of the miRNA regulators of the human ovulatory cascade. Sci Rep. 2018;8(1):15605.
M Drake J, R Friis R, Dharmarajan A. The role of sFRP4, a secreted frizzled-related protein, in ovulation 2003. 389–97 p.
Abir R, Ao A, Jin S, Barnett M, Van den Hurk R, Freimann S, et al. Immunocytochemical detection and reverse transcription polymerase chain reaction expression of oncostatin M (OSM) and its receptor (OSM-Rβ) in human fetal and adult ovaries. Fertil Steril. 2005;83(4):1188–96.
Devjak R, Fon Tacer K, Juvan P, Virant Klun I, Rozman D, Vrtačnik BE. Cumulus cells gene expression profiling in terms of oocyte maturity in controlled ovarian Hyperstimulation using GnRH agonist or GnRH antagonist. PLoS One. 2012;7(10):e47106.
Iager AE, Kocabas AM, Otu HH, Ruppel P, Langerveld A, Schnarr P, et al. Identification of a novel gene set in human cumulus cells predictive of an oocyte's pregnancy potential. Fertil Steril. 2013;99(3):745–52.e6.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.
Yerushalmi GM, Salmon-Divon M, Yung Y, Maman E, Kedem A, Ophir L, et al. Characterization of the human cumulus cell transcriptome during final follicular maturation and ovulation. Mol Hum Reprod. 2014;20(8):719–35.
Molinari E, Bar H, Pyle AM, Patrizio P. Transcriptome analysis of human cumulus cells reveals hypoxia as the main determinant of follicular senescence. Mol Hum Reprod. 2016;22(8):866–76.
Hannah MB, Kylie RD, Melanie S-M, Robert BG, Jeremy GT, Darryl LR. Failure to launch: aberrant cumulus gene expression during oocyte in vitro maturation. Reproduction. 2017;153(3):R109–R20.
Ben-Ami I, Komsky A, Bern O, Kasterstein E, Komarovsky D, Ron-El R. In vitro maturation of human germinal vesicle-stage oocytes: role of epidermal growth factor-like growth factors in the culture medium. Hum Reprod. 2011;26(1):76–81.
Dhali A, Javvaji PK, Kolte AP, Francis JR, Roy SC, Sejian V. Temporal expression of cumulus cell marker genes during in vitro maturation and oocyte developmental competence. J Assist Reprod Genet. 2017;34(11):1493–500.
Wei L, Liang X-W, Zhang Q-H, Li M, Yuan J, Li S, et al. BubR1 is a spindle assembly checkpoint protein regulating meiotic cell cycle progression of mouse oocyte. Cell Cycle. 2010;9(6):1112–21.
Maman E, Yung Y, Cohen B, Konopnicki S, Dal Canto M, Fadini R, et al. Expression and regulation of sFRP family members in human granulosa cells. Mol Hum Reprod. 2011;17(7):399–404.
Xin H, Cuifang H, Xiaofang S, Xiaoyan L, Yinghua S, Yuhua Z, et al. Differences in the transcriptional profiles of human cumulus cells isolated from MI and MII oocytes of patients with polycystic ovary syndrome. Reproduction. 2013;145(6):597–608.
Martins KR, Haas CS, Ferst JG, Rovani MT, Goetten ALF, Duggavathi R, et al. Oncostatin M and its receptors mRNA regulation in bovine granulosa and luteal cells. Theriogenology. 2019;125:324–30.
McCord LA, Li F, Rosewell KL, Brännström M, Curry TE. Ovarian expression and regulation of the stromelysins during the periovulatory period in the human and the rat. Biol Reprod. 2011;86(3):78.
Yasuda K, Hagiwara E, Takeuchi A, Mukai C, Matsui C, Sakai A, et al. Changes in the distribution of Tenascin and Fibronectin in the mouse ovary during Folliculogenesis, atresia, Corpus Luteum formation and Luteolysis. Zool Sci. 2005;22(2):237–45 9.
Bagavandoss P. Temporal expression of tenascin-C and type I collagen in response to gonadotropins in the immature rat ovary. Acta Histochem. 2014;116(7):1125–33.
Hernandez-Gonzalez I, Gonzalez-Robayna I, Shimada M, Wayne CM, Ochsner SA, White L, et al. Gene expression profiles of cumulus cell oocyte complexes during ovulation reveal cumulus cells express neuronal and immune-related genes: does this expand their role in the ovulation process? Mol Endocrinol. 2006;20(6):1300–21.
Lu CH, Lee RK, Hwu YM, Lin MH, Yeh LY, Chen YJ, et al. Involvement of the serine protease inhibitor, SERPINE2, and the urokinase plasminogen activator in cumulus expansion and oocyte maturation. PLoS One. 2013;8(8):e74602.
Pan B, Li J. MicroRNA-21 up-regulates metalloprotease by down-regulating TIMP3 during cumulus cell-oocyte complex in vitro maturation. Mol Cell Endocrinol. 2018;477:29–38.
Sayasith K, Lussier J, Sirois J. Molecular characterization and transcriptional regulation of a disintegrin and metalloproteinase with thrombospondin motif 1 (ADAMTS1) in bovine preovulatory follicles. Endocrinology. 2013;154(8):2857–69.
Young KA, Tumlinson B, Stouffer RL. ADAMTS-1/METH-1 and TIMP-3 expression in the primate corpus luteum: divergent patterns and stage-dependent regulation during the natural menstrual cycle. Mol Hum Reprod. 2004;10(8):559–65.
Russell DL, Robker RL. Molecular mechanisms of ovulation: co-ordination through the cumulus complex. Hum Reprod Update. 2007;13(3):289–312.
Tsuji Y, Tamaoki TH, Hasegawa A, Kashiwamura S, Iemoto A, Ueda H, et al. Expression of interleukin-18 and its receptor in mouse ovary. Am J Reprod Immunol. 2001;46(5):349–57.
Wang S, Liu W, Pang X, Dai S, Liu G. The Mechanism of Melatonin and Its Receptor MT2 Involved in the Development of Bovine Granulosa Cells. International Journal of Molecular Sciences. 2018;19(7). https://www.mdpi.com/about/announcements/784.
Myers M, van den Driesche S, McNeilly AS, Duncan WC. Activin a reduces luteinisation of human luteinised granulosa cells and has opposing effects to human chorionic gonadotropin in vitro. J Endocrinol. 2008;199(2):201–12.
Murphy MJ, Halow NG, Royer PA, Hennebold JD. Leukemia inhibitory factor is necessary for ovulation in female rhesus macaques. Endocrinology. 2016;157(11):4378–87.
Sessions-Bresnahan DR, Carnevale EM. Age-associated changes in granulosa cell transcript abundance in equine preovulatory follicles. Reprod Fertil Dev. 2015;27(6):906–13.
Park JY, Su YQ, Ariga M, Law E, Jin SL, Conti M. EGF-like growth factors as mediators of LH action in the ovulatory follicle. Science. 2004;303(5658):682–4.
Light A, Hammes SR. LH-induced Steroidogenesis in the mouse ovary, but not testis, requires matrix metalloproteinase 2- and 9-mediated cleavage of Upregulated EGF receptor ligands. Biol Reprod. 2015;93(3):65.
Choi Y, Wilson K, Hannon PR, Rosewell KL, Brännström M, Akin JW, et al. Coordinated regulation among progesterone, prostaglandins, and EGF-like factors in human ovulatory follicles. J Clin Endocrinol Metab. 2017;102(6):1971–82.
Nuttinck F, Gall L, Ruffini S, Laffont L, Clement L, Reinaud P, et al. PTGS2-related PGE2 affects oocyte MAPK phosphorylation and meiosis progression in cattle: late effects on early embryonic development. Biol Reprod. 2011;84(6):1248–57.
Marei WF, Abayasekara DR, Wathes DC, Fouladi-Nashta AA. Role of PTGS2-generated PGE2 during gonadotrophin-induced bovine oocyte maturation and cumulus cell expansion. Reprod BioMed Online. 2014;28(3):388–400.
Fang L, Yu Y, Zhang R, He J, Sun YP. Amphiregulin mediates hCG-induced StAR expression and progesterone production in human granulosa cells. Sci Rep. 2016;6:24917.
Feuerstein P, Cadoret V, Dalbies-Tran R, Guerif F, Bidault R, Royere D. Gene expression in human cumulus cells: one approach to oocyte competence. Hum Reprod. 2007;22(12):3069–77.
Ashkenazi H, Cao X, Motola S, Popliker M, Conti M, Tsafriri A. Epidermal growth factor family members: endogenous mediators of the ovulatory response. Endocrinology. 2005;146(1):77–84.
Calder MD, Caveney AN, Smith LC, Watson AJ. Responsiveness of bovine cumulus-oocyte-complexes (COC) to porcine and recombinant human FSH, and the effect of COC quality on gonadotropin receptor and Cx43 marker gene mRNAs during maturation in vitro. Reprod Biol Endocrinol. 2003;1:14.
Guzman L, Adriaenssens T, Ortega-Hrepich C, Albuz FK, Mateizel I, Devroey P, et al. Human antral follicles <6 mm: a comparison between in vivo maturation and in vitro maturation in non-hCG primed cycles using cumulus cell gene expression. Mol Hum Reprod. 2013;19(1):7–16.
Xie M, Li M, Zhou J, Ding X, Shao Y, Jing J, et al. Brain-derived neurotrophic factor promotes human granulosa-like tumor cell steroidogenesis and proliferation by activating the FSH receptor-mediated signaling pathway. Sci Rep. 2017;7(1):180.
Järvensivu P, Saloniemi-Heinonen T, Awosanya M, Koskimies P, Saarinen N, Poutanen M. HSD17B1 expression enhances estrogen signaling stimulated by the low active estrone, evidenced by an estrogen responsive element-driven reporter gene in vivo. Chem Biol Interact. 2015;234:126–34.
de Souza DK, Salles LP, Camargo R, Gulart LVM, Costa E, Silva S, de Lima BD, et al. Effects of PI3K and FSH on steroidogenesis, viability and embryo development of the cumulus-oocyte complex after in vitro culture. Zygote. 2018;26(1):50–61.
Leon PM, Campos VF, Kaefer C, Begnini KR, McBride AJ, Dellagostin OA, et al. Expression of apoptotic genes in immature and in vitro matured equine oocytes and cumulus cells. Zygote. 2013;21(3):279–85.
Kumazawa Y, Kawamura K, Sato T, Sato N, Konishi Y, Shimizu Y, et al. HCG up-regulates survivin mRNA in human granulosa cells. Mol Hum Reprod. 2005;11(3):161–6.
Ni XR, Sun ZJ, Hu GH, Wang RH. High concentration of insulin promotes apoptosis of primary cultured rat ovarian granulosa cells via its increase in extracellular HMGB1. Reprod Sci. 2015;22(3):271–7.
Du XH, Zhou XL, Cao R, Xiao P, Teng Y, Ning CB, et al. FSH-induced p38-MAPK-mediated dephosphorylation at serine 727 of the signal transducer and activator of transcription 1 decreases Cyp1b1 expression in mouse granulosa cells. Cell Signal. 2015;27(1):6–14.
Liu J, Tian Y, Ding Y, Heng D, Xu K, Liu W, et al. Role of CYP51 in the regulation of T3 and FSH-induced Steroidogenesis in female mice. Endocrinology. 2017;158(11):3974–87.
Patel SS, Beshay VE, Escobar JC, Suzuki T, Carr BR. Molecular mechanism for repression of 17alpha-hydroxylase expression and androstenedione production in granulosa cells. J Clin Endocrinol Metab. 2009;94(12):5163–8.
Nogueira D, Cortvrindt R, Everaerdt B, Smitz J. Effects of long-term in vitro exposure to phosphodiesterase type-3 inhibitors on follicle and oocyte development. Reproduction. 2005;130(2):177–86.
Thomas RE, Thompson JG, Armstrong DT, Gilchrist RB. Effect of specific phosphodiesterase isoenzyme inhibitors during in vitro maturation of bovine oocytes on meiotic and developmental capacity. Biol Reprod. 2004;71(4):1142–9.
Assou S, Haouzi D, Mahmoud K, Aouacheria A, Guillemin Y, Pantesco V, et al. A non-invasive test for assessing embryo potential by gene expression profiles of human cumulus cells: a proof of concept study. Mol Hum Reprod. 2008;14(12):711–9.
Hamel M, Dufort I, Robert C, Gravel C, Leveille MC, Leader A, et al. Identification of differentially expressed markers in human follicular cells associated with competent oocytes. Hum Reprod. 2008;23(5):1118–27.
Racowsky C, Needleman DJ. Cumulus cell gene expression as a potential biomarker for oocyte quality. Fertil Steril. 2018;109(3):438–9.
Fuchs Weizman N, Wyse BA, Gat I, Balakier H, Sangaralingam M, Caballero J, et al. Triggering method in assisted reproduction alters the cumulus cell transcriptome. Reprod Biomed Online. 2019;39(2):211–24.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS One. 2010;5(11):e13984.
Ovarian Kaleidoscope Database. 1999. Available from: http://okdb.appliedbioinfo.net/. Accessed 17 Jul 2020.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.
The authors would like to thank all CReATe Fertility Centre patients for the donation of their material to this study, Janice Montbriand for her help with statistical analysis, all embryologists, nurses, and staff for helping with data collection.
This project was supported through reinvestment of clinical earnings at CReATe Fertility Clinic.
Ethics approval and consent to participate
This study was approved by the University of Toronto Research Ethics Board (Approval #29237). All subjects provided written informed consent for the donation of their biological waste material, which included collection of cumulus cells, and extraction of de-identified clinical information including age, body mass index, treatment regime and laboratory outcomes.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of primers and probes used for qPCR validation of RNAseq results. Supplemental Table S2: RNAseq quality control metrics for each cohort presented as mean ± SEM (Range). Supplemental Table S3: All differentially expressed genes between MII-CC and GV-CC cohorts. Supplemental Table S4: Differentially expressed genes known to be regulated by LH, FSH, or both with their associated fold change, and number of gene sets that are found in following leading edge analysis. Supplemental Table S5: Previously identified genes implicated in oocyte maturation in vitro that were not differentially expressed in this study. Supplemental Table S6: Top 50 genes identified by leading edge analysis ranked by the number of unique gene-sets the gene was identified in. Supplemental Table S7: Normalized read counts of all samples.
About this article
Cite this article
Wyse, B.A., Fuchs Weizman, N., Kadish, S. et al. Transcriptomics of cumulus cells – a window into oocyte maturation in humans. J Ovarian Res 13, 93 (2020). https://doi.org/10.1186/s13048-020-00696-7
- Cumulus cells
- Cumulus-oocyte complex
- Gene expression
- Oocyte maturation
- Assisted reproductive technology