Skip to main content

Comparative analyses in transcriptome of human granulosa cells and follicular fluid micro-environment between poor ovarian responders with conventional controlled ovarian or mild ovarian stimulations



Both mild and conventional controlled ovarian stimulation are the frequently used protocols for poor ovarian responders. However, there are some debates about which treatment is better. Moreover, little is known about the follicular physiology after the two ovarian stimulation protocols. This study was intended to investigate the features in granulosa cells and follicular fluid micro-environment after the two different ovarian stimulation protocols in poor responders.


Granulosa cells RNA were sequenced using Illumina Hiseq technology. Specific differently expressed genes and proteins were verified by real-time quantitative PCR and Western blot analysis. Moreover, hormone and cytokine concentrations in the follicular fluid were measured by electrochemiluminescence immunoassay and enzyme-linked immunoabsorbent assay. The correlation between the results of molecular experiments and the laboratory outcomes were analyzed by Spearman correlation analysis.


The differentially expressed genes between the two groups were involved in 4 signaling pathways related to the follicular development; three proteins pertinent to the TGF-β signaling pathway were expressed differently in granulosa cells between the two, and the constituents in the follicular fluid were also different. Further, a correlation between the TGF-β signaling pathway and the good-quality embryo was observed.


The present study made a comparison for the first time in the transcriptome of human granulosa cells and the follicular fluid micro-environment between poor responders with the conventional controlled ovarian stimulation or the mild ovarian stimulation, showing that the TGF-β signaling pathway may correlate with the good-quality of embryos in the mild group, which may be instrumental to the choice of optimal management for IVF patients.


In in vitro fertilization (IVF) practice, a noticeable proportion of women receiving the conventional controlled ovarian stimulation (COS) undergo poor ovarian response (POR), suffering low numbers of retrieved oocytes and transferrable embryos, and low pregnancy rate [1]. The situation of POR patients is critically challenging the ART clinicians although many adjunctive strategies for the COS protocol have been developed, e.g., increasing gonadotropin (Gn) dose, supplementing with exogenous luteinizing hormone (LH), decreasing gonadotropin-releasing hormone agonist (GnRH-a) dose and using adjunctive growth hormone [2,3,4,5]. Recently, the mild ovarian stimulation highlights in clinicians’ focus and is favored by some owing to its lower economic burden on patients. However, it remains debatable which protocol is clinically much better, thus, it is necessary to investigate the different influence on the follicular development under the conditions of the conventional COS and the mild ovarian stimulation. The major difference between the conventional COS and the mild ovarian stimulation lies on the daily used dosage of Gn, a key regulator in the follicular development. In vitro experiments showed that an appropriate dose of Gn stimulated follicle growth while a very high dose of Gn led to a decreased follicle survival and the contemporary accumulation of hazardous metabolic products [6], which even made detrimental effects on oocytes and embryos [7,8,9]. Yet, the mild ovarian stimulation, aiming at a more physiological response, may attenuate such kind of damage.

A series of architectural transformations in follicles during folliculogenesis depend on a complicated and subtle cross-talking mechanism among the three of granulosa cells, the follicular fluid (FF) micro-environment and oocytes. Some studies reported that granulosa cells affected the growth, development [10], energy uptake [11] and transcription process of the oocyte [12]. Besides, the FF micro-environment, the most immediate niche rich in various hormones and cytokines, acts as important communication media in the antral follicle and transports nutrients [13] for the developing oocytes. Thus, different stimulation protocols, influencing the gene expression in granulosa cells and the FF micro-environment differently [14, 15], may have a significant effect on the cell signaling and on the oocyte quality potentially [16]. Because it is infeasible to consume oocytes freely from patients, the granulosa cell and the FF represent two indispensable factors to fathom folliculogenesis. Theoretically, the transcription information from granulosa cells, which are the somatic cell most adjacent to the oocyte, may mirror the developmental competence of the associated oocyte. Therefore, comparing the FF micro-environment and the gene expression of granulosa cells in POR patients with different ovarian stimulation protocols is of great significance for understanding the effects of different protocols on the follicular physiology.

Here, aiming to provide reliable physiological data for selecting optimal management for POR patients, we compared the pertinent data or measurements of granulosa cells, the FF micro-environment and oocytes from POR patients treated with the conventional COS and the mild ovarian stimulation.


Study design and participants

The prospective study, conducted at the Reproductive Medicine Research Center, the Sixth Affiliated Hospital of Sun Yat-sen University, was approved by the Ethics Committee of Reproductive Medicine and Prenatal Diagnosis, the Sixth Affiliated Hospital of Sun Yat-sen University (2019ZSLYEC-002S), and all participants rendered written informed consent.

The POR patients, recruited according to the Bologna criteria [17], received either the mild ovarian stimulation (the mild group) or the conventional COS (the COS group) for the ART treatment. Patients in the mild group, starting with 5 mg of letrozole (Letrozole Tablets, Hengrui Medicine, China) daily from the 3rd to 7th day of the menstrual cycle in conjunction with the injection of 150 IU of recombinant follicle stimulation hormone (FSH) (Puregon, MSD, Germany) on the 4th and 6th day, were administered 150 IU of recombinant FSH since the 8th day of the menstrual cycle until the day of human chorionic gonadotropin (hCG) administration. And when the serum estradiol (E2) level was ≥ 200 pg/ml, a gonadotropin-releasing hormone antagonist (GnRH-ant, Cetrotide, Merck-Serono, Switzerland, or Orgalutran, MSD, Germany) was given subcutaneously (0.25 mg/day). The COS group was instituted in terms of a GnRH-a stop protocol with a 14-day mid-luteal down-regulation (0.05 mg/day; Triptorelin, Ferring, Germany), followed by daily stimulation with 300 IU of recombinant FSH (Puregon, MSD, Germany). The final maturation of oocytes was induced with 0.25 mg of hCG (Ovidrel, Merck Serono S.p.A., Italy) when at least one leading follicle reached 18 mm, or 3 follicles reached 17 mm in diameter. Oocyte retrieval was performed transvaginally under ultrasound guidance 36–37 h past the hCG injection. Standard IVF or ICSI was performed as appropriate. The patients’ demographic and clinical characteristics, including age, body mass index (BMI) and the serum basal sexual hormone levels, were recorded for following applications.

Embryo evaluation

Embryo evaluation methods are as described previously [18]. Day-3 embryo morphology was scored according to blastomere number, the degree of fragmentation and extent of asymmetry [19]. Briefly, each embryo was assigned a fragmentation score of 0%, 1–9%, 10–25%, 26–50%, or > 50%, and an asymmetry score of perfect, moderate or severe asymmetry. A good-quality embryo was defined as having ≥ 6 cells with < 10% fragmentation and either no asymmetry or moderate asymmetry. A transferrable embryo was defined as having ≥ 4 cells with < 26% fragmentation and either no asymmetry or moderate asymmetry.

The quality of blastocyst stage embryos was assessed according to the criteria of Gardner and Schoolcraft [20] based on the degree of expansion and hatching status of the blastocoel cavity (1-6), the size of the inner cell mass (A–C) and the development of the trophectoderm (A–C). A good-quality blastocyst was defined as a quality score ≥ 3BB.

Collection of follicular fluid (FF) and granulosa cells

The FF and granulosa cells were obtained via ultrasound-guided transvaginal oocyte retrieval. The FF, collected from the largest follicle, was transferred to a 2-ml cryogenic vial (Corning, New York, USA) and then stored at -80℃ for further analyses of hormones and cytokines. The FF from the remaining follicles, pooled in a 50-ml conical centrifuge tube (Corning, New York, USA) after removal of oocyte-cumulus complexes, was centrifuged for 5 min at 500 g for the granulosa cell isolation. The granulosa cells were washed three times with cold sterile phosphate buffered saline (PBS) at 500 g for 5 min at room temperature (RT). Three ml of PBS were added for resuspending the cells and then 3 ml of Ficoll (Tianjin HaoYang Biological Manufacture Co., Ltd, Tianjin, China) were carefully added along the tube wall for gradient centrifugation. The samples were centrifuged at 1000 g for 25 min at RT for the recollection of granulosa cells, which, appearing in the middle interface, were pipetted and transferred to a 1.5-ml Eppendorf tube sitting on ice. Then, 200 μl of cold sterile PBS and 600 μl of red blood cell lysis buffer (CoWin Biosciences, Beijing, China) were added and mixed by inversion, and placed on ice for 15 min. After centrifugation at 1000 g for another 10 min at RT, the granulosa cells were washed again with cold sterile PBS for final centrifugation, and the cell precipitates were stored at -80℃ until RNA extraction [21,22,23]. Preceding the study, we performed the identification experiments on the extracted cells (Fig. S1).

RNA extraction

For independent extraction of the total RNA, the granulosa cells sampled from each individual patient were treated and incubated with 500 μl of TRIzon reagent according to the instruction (CoWin Biosciences, Beijing, China). Then, 100 μl of chloroform were added to the solution and mixed thoroughly. After phase separation and centrifugation, the aqueous layer was carefully pipetted into another RNase-free tube, and 500 μl of isopropanol were added to the aqueous phase and incubated for 10 min at RT for the following centrifugation and precipitation of the total RNA. The precipitate was resuspended in 500 μl of 75% ethanol and centrifuged to get RNA pellet. The RNA purity and integrity were measured by using the Nanodrop ND-1000 (NanoDrop Technologies, Wilmington, USA) and the Agilent 2200 TapeStation (Agilent Technologies, Santa Clara, USA), respectively.

cDNA Library construction and transcriptome sequencing

Fragmented RNAs, approximately 200 bp in average, were subjected to the first strand and second strand cDNA synthesis followed by adaptor ligation and PCR enrichment with a low-cycle according to instructions of NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (New England Biolabs, Ipswich, USA). The purified library products were evaluated using the Agilent 2200 TapeStation and Qubit3.0 (Life Technologies, Carlsbad, USA). The libraries were paired-end sequenced (PE150, sequencing reads = 150 bp) at Guangzhou RiboBio Co., Ltd (Guangzhou, China) using the Illumina Hiseq X-ten platform.

Gene expression analysis

The differential expression analysis of the two groups was performed using the DEGSeq2 R package (V1.18.1). DEGSeq analysis was used to provide statistical routines for determining differentially expressed genes (DEGs) using a model based on the negative binomial distribution. P-values were adjusted (Padj), which was calculated via the Benjamini–Hochberg method to exclude false-positive results. Genes with Padj < 0.05 were defined as being differentially expressed. All differentially expressed genes were selected for the GO (Gene Ontology) and the KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis. The GO enrichment analysis was performed using the KOBAS3.0 software. GO terms with corrected P-values less than 0.05 were considered to be significantly enriched. Enrichment of DEGs in the KEGG was analyzed using the KOBAS3.0 software.

Real-time quantitative PCR (RT-qPCR) analysis

The significant DEGs were confirmed by RT-qPCR. The primers (Table 1) were designed and synthesized by Sangon Biotech (Shanghai, China) using sequence information obtained from the NCBI Gene Database ( and blasted ( The total RNA from granulosa cells was reverse-transcribed into cDNA using PrimeScript™ RT Master Mix (Perfect Real Time) (Takara, Kusatsu, Japan). RT-qPCR was performed on LightCycler® 480 System (Roche, Basel, Switzerland) using 2 × RealStar Green Fast Mixture (GenStar, Beijing, China). All the genes were normalized to the gene of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and results expressed as relative fold changes using the 2–ΔΔCT method.

Table 1 Primers used for RT-PCR

Western blot analysis

According to the sequencing data and functional analysis, we analyzed the proteins associated with the transforming growth factor beta (TGF-β) signaling pathway in granulosa cells using the Western blot. Briefly, the whole protein samples were extracted, and their concentrations measured using the BCA Protein Assay Kit (GenStar, Beijing, China). Then, equal amounts of proteins (20 μg per lane) were loaded and separated on a 12% sodium dodecyl sulfate (SDS) polyacrylamide gel. Primary antibodies were used after proteins transferred to the 0.45-μm Immobilon-P PVDF Membrane, and GAPDH used as loading control. Signals were obtained in the linear range of detection and quantified with the Bio-Rad ChemiDoc Imaging System. The grey values of target bands were analyzed by ImageJ software and normalized to the reference band. The primary antibodies used were TGF-β2 (19,999–1-AP; 1:2000 diluted; Proteintech), TGF-βR2 (66,636–1-lg; 1:2000 diluted; Proteintech), Smad2 (12,570–1-AP; 1:2000 diluted; Proteintech), Smad3 (25,494–1-AP; 1:2000 diluted; Proteintech).

Hormone assay of follicular fluid and correlation analysis

The FF was thawed to RT and fully mixed before analysis. The FF anti-müllerian hormone (AMH), FSH, LH, E2, progesterone (P), testosterone (T) and prolactin (PRL) were measured by electrochemiluminescence immunoassay (ECLIA) following the instructions (Roche Diagnostics GmbH, Mannheim, Germany). The TGF-β2, bone morphogenetic protein-15 (BMP-15) and growth differentiation factor-9 (GDF-9) concentrations were detected using enzyme-linked immunosorbent assay (ELISA). And the correlation analysis between the cytokine concentrations and the laboratory outcomes was performed.

Statistical analysis

All clinical parameters and intrafollicular hormonal concentration values were expressed as means ± standard deviation (SD). Differences in the laboratory outcomes and the intrafollicular hormone concentrations between the two groups were compared via Student’s t-test. The statistical analyses of the Western blot and RT-qPCR results were performed using Student’s t-test. Spearman’s correlation coefficient was calculated to detect the correlation between the TGF-β2 level in the FF and the rate of the good-quality embryo. Analyses were performed using SPSS (version 23, Chicago, USA), and significance defined as a P value < 0.05.


Clinical characteristics and outcomes of patients

Table 2 shows the clinical characteristics of the recruited 48 POR patients, who were categorized into two groups, i.e., the mild group (n = 27) and the COS group (n = 21), with no statistical differences in age, BMI, serum AMH and basal hormonal levels. Compared with the COS group, both the Gn stimulation duration and the total Gn dosage were significantly shorter and lower in the mild group, and also the serum E2 level on trigger day statistically lower. Although a higher number of retrieved oocytes were obtained in the COS group, no statistical differences observed regarding the numbers of the metaphase II oocyte, 2 pronuclei (2PN), transferrable embryo and good-quality embryo. Noteworthily, significantly higher rates of the metaphase II oocyte and the good-quality embryo were observed in the mild group (Table 2).

Table 2 Clinical characteristics and outcomes of the two groups

Transcriptional profiles between samples from the two groups

We obtained 55–67 million 150 bp reads for each sample, and 92 Gb of raw data in total were obtained for all the samples. Clean data, which were used for subsequent analyses, were ~ 8 Gb per sample (Supplementary data Table S1). We calculated the Pearson correlation coefficient of every two samples by using the gene expression. As shown in Fig. 1, the expression patterns were generally homogeneous between every two samples, all R2 > 0.9 except for COS_3 and COS_4 (R2 = 0.896).

Fig. 1
figure 1

Pearson correlation analysis between samples. The color key from blue to white indicates the correlation from high to low

Analyses of DEGs and function enrichment

The clustering analysis was used to determine the DEGs (Fig. 2A) between the two groups, and totally 425 genes were found to be differentially expressed, with 192 being up-regulated and 233 down-regulated in the mild group (P < 0.05, absolute value of log2(Foldchange) > 1) (Fig. 2B).

Fig. 2
figure 2

Cluster and filter analyses of DEGs. A Heatmap of the DEGs between the two groups. The color key from blue to red indicates the relative gene expression level from low to high. B Volcano map showing the DEGs. The x-axis shows the foldchange in gene expression, and the y-axis shows significant statistical differences. Red dots represent up-regulated genes, and green dots represent down-regulated genes

The GO enrichment analysis was performed to explore the biological functions of the 425 DEGs, and three categories of biological functions were identified, including biological process, cellular component, and molecular function. Five representative GO enrichment results in each category are presented in Fig. 3. Fig. 3A shows the up-regulated DEGs in the mild group, which are related to cytokine activity and regulation activity, and Fig. 3B shows the down-regulated DEGs in the mild group, which are involved in the function of immune response.

Fig. 3
figure 3

GO enrichment analysis of DEGs between the two groups. The x-axis shows the five representative enrichment terms in three categories, including biological process in green bars, cellular component in blue bars and molecular function in red bars. The y-axis represents -log10(P value). Compared with COS group, A genes up-regulated in mild group, and B genes down-regulated in mild group

The 425 DEGs were also assigned to the KEGG enrichment analysis, and Table 3 shows the top ten pathways (P < 0.05), among which are some associated with the oocyte development, including the cytokine-cytokine receptor interaction, TGF-β signaling pathway, cGMP-PKG signaling pathway, and metabolic pathway.

Table 3 Top ten KEGG pathways of DEGs between the two groups

RT-qPCR validation

Based on the sequencing data and functional analyses, 20 genes were selected for validation with RT-qPCR not only according to the most significant fold change variations but also in terms of the role they play in the pathway associated with the oocyte formation and development (Table 1). Of the 20 genes, 6 genes (TGF-β2, CLCF1, NPR1, ADRA1D, DAN, CCL26) exhibited by RNA sequencing in the mild group were confirmed in upregulation and 1 gene (IL21R) confirmed downregulation (Fig. 4).

Fig. 4
figure 4

Validation of the RNA sequencing results using RT-qPCR. A Bars extending to the right or left of zero coordinate represent up or down regulated genes, respectively. Black bars represent expression of transcript in RNA sequencing while grey bars represent expression of the same transcript presented by RT-qPCR. B Comparison of expression of 7 selected genes between the two groups. Asterisk (*) above the bar represents statistical difference in RT-qPCR experiments. *P < 0.05; ***P < 0.001

Protein expressions in granulosa cells

To determine the different protein expression levels of the TGF-β signaling pathway in granulosa cells between the two groups, such protein products as TGF-β2, TGF-βR2, Smad2 and Smad3 were detected. The expression of TGF-β2, TGF-βR2 and Smad3 in the mild group were significantly upregulated compared with the COS group (Fig. 5).

Fig. 5
figure 5

Western blot analysis of TGF-β signaling pathway specific proteins TGF-β2, TGF-βR2, Smad2 and Smad3 in the mild and COS groups. GAPDH used as a control. A Western blot results. B Quantity analysis based on the Western blot results. Mild group (n = 10), and COS group (n = 8). *P < 0.05; **P < 0.01; ***P < 0.001; NS No significant

Follicular fluid analysis and correlation analysis

Different hormone and cytokine concentrations in the follicular fluid (FF) of the two groups are listed in Table 4. Compared with the COS group, the FSH, PRL and P levels in the FF were significantly lower in the mild group while the LH and T levels significantly higher. As for the cytokine concentrations, a significantly higher TGF-β2 concentration and a significantly lower GDF-9 level were detected in the mild group while the BMP-15 levels were comparable between the two groups. Correlation analysis revealed that there was a positive correlation between the TGF-β2 concentration and the rate of the good-quality embryo (Fig. 6).

Table 4 Hormone and cytokine concentrations in FF between the two groups
Fig. 6
figure 6

Spearman correlation analysis between the TGF-β2 level in the FF and the rate of the good-quality embryo. The scatter diagrams of the two groups show the monotonic relationship between the TGF-β2 level and the rate of the good-quality embryo, r is correlation coefficient and P < 0.05 indicates that the correlation is statistically significant


The follicle, underlying the follicular development, oocyte maturation, cumulus expansion and ovulation, is a histological and functional architecture supporting a series of complicated, subtle and continuous cross-talkings among the oocyte, the somatic follicular cells (particularly the granulosa cell) and the FF micro-environment [24] (Fig. 7). Granulosa cells have close relationship with the oocyte, affecting the oocyte development and contributing to the FF micro-environment through a series of dynamic regulations in levels of the transcription, protein and metabolism [25]. The FF components may be directly altered by the hormonal, paracrine and autocrine signaling pathways [26], which, further, has been suggested to influence the oocyte quality, early embryo development and subsequent pregnancy [27, 28]. Nowadays, as an indispensable medicine for ovarian stimulation in the ART practice, exogenous Gn is used to promote follicular growth via its influence on the cross-talking among the oocyte, granulosa cell and FF. Gn dose is the main difference in our two protocols, but the effect of this difference on the oocyte development and the mechanism underlying which it affects the oocyte development remain unclear although Lu et al. showed that Gn stimulation may potentially induce meiotic errors of human oocytes through analyzing the transcriptome of the granulosa cell after natural and Gn stimulation cycles [29].

Fig. 7
figure 7

Schematic diagram of the results in the background of cross-talkings within the follicle. The red and blue arrows represent the COS group and the mild group, respectively. The up and down arrows represent up-regulation and down-regulation, respectively. a Histologic and functional architecture of a mature follicle. Panels b, c, and d summarize the main findings of this study. b The KEGG analysis of DEGs between the two groups. After sequencing and analyzing the DEGs of granulosa cells from the two groups, the KEGG analysis showed four signaling pathways associated with the oocyte development. c Comparative analyses of the rate of good-quality embryos as well as the hormones and cytokines in the FF. Significantly higher rates of good-quality embryos were observed in the mild group. Compared with the COS group, the FSH and P levels in the FF were significantly lower in the mild group while the LH and T levels significantly higher. As for the cytokine concentrations, a significantly higher TGF-β2 concentration and a significantly lower GDF-9 level were detected in the mild group. d Comparative analyses of TGF-β2, TGF-βR2 and Smad3 in granulosa cells as well as TGF-β2 in the FF. The expressions of TGF-β2, TGF-βR2 and Smad3 in the mild group were significantly upregulated compared with the COS group

The results of this study showed that the DEGs between the two protocols were mainly involved in the TGF-β signaling pathway, cytokine-cytokine receptor interactions, metabolic pathways and cGMP-PKG signaling pathway, all of which are involved in the developmental processes of the follicle and oocyte. The RNA sequencing result, qRT-PCR result and Western blot analysis revealed the up-regulation of the genes and proteins related to the TGF-β signaling pathway in the mild group, and a higher rate of the good-quality embryo also observed in the mild group. The correlation analysis further confirmed a positive correlation between the TGF-β2 level in the FF and the rate of the good-quality embryo. These results imply that the TGF-β signaling pathway may be closely associated with the embryo quality, which is consistent to literatures reporting that the TGF-β signaling pathway was directly related to early embryonic development and blastocyst formation in bovine, and affected the oocyte by mediating the proliferation of granulosa cells [30,31,32,33].

As an important link within the cross-talking among the oocyte, granulosa cell and FF, the TGF-β signaling pathway is under complicated regulations of the FF components, which has not been completely addressed. In this study, lower intrafollicular FSH was found in the mild group due to the use of fewer doses of Gn. Although the clinical significance of lower FSH remains uncertain, previous study pointed out that high-dose Gn stimulation could impair embryo development potential [34], which is supported by our study, showing a higher rate of the good-quality embryo in the mild group. And, in our experimental system, through detecting the granulosa cell and the FF from the same individual sample, observed in the mild group were a significant up-regulation of the TGF-β signaling pathway-related genes and proteins in the cells, and a decrease of FSH along with the increase of LH and T in the FF. However, Gueripel et al. showed that the increase of FSH and LH could up-regulate the TGF-β signaling pathway in immature mice after Gn stimulation [35]. This inconsistence can be explained in the following aspects. First, a major difference between us lies in that we used letrozole in the mild group. Letrozole, an aromatase inhibitor blocking the transformation of androstenedione into estrogen, will cause different T level in the FF, and T could affect the expression of TGF-β signaling [36]. Second, the granulosa cell and the FF in our study were sampled during ovulation while their samples were from the follicular phase, and the different time windows will present different observation facts because the sexual hormones relevant to the follicular development appear and effect in terms of a periodical and dynamic regulating axis. Thirdly, the different experimenting systems, such as measurement methods and experiment species, may also contribute to the inconsistence.

We also measured the cytokines relevant to the oocyte development in the FF, i.e., TGF-β2, BMP-15 and GDF-9. Comparable BMP-15 levels were found in the two groups while a significantly lower GDF-9 and a significantly higher TGF-β2 in the FF were found in the mild group, whose oocyte maturity rate and good-quality embryo rate were both higher. This observation is inconsistent to the previous reports [37, 38] showing that higher GDF-9 and BMP-15 levels in the FF were significantly correlated with a higher oocyte maturation and a better embryo quality, which can be explained as follows: 1) The increase of GDF-9 is positively correlated with the increase of FSH [39], and the higher level of FSH in the COS group accounted for the higher level of GDF-9; 2) Only if the heterodimer between GDF-9 and BMP-15 has been structured was there any biological activity, such as activating downstream molecules, promoting the proliferation of granulosa cells and regulating the growth and development of oocytes [40]. Therefore, the measurement of GDF-9 or BMP-15 levels alone does not accurately reflect the quality of embryos.

Here are three limitations in this study. First, only the mural granulosa cells were sampled for the transcriptome data analysis while the cumulus granulosa cells were spared because of the difficulty in obtaining enough cells for sequencing. Second, morphologic assessment rather than preimplantation genetic testing (PGT) was applied to rating the quality of embryos because PGT can only be used for such specific patients as with genetic diseases, recurrent pregnancy loss, repeated IVF failure, etc., according to the regulations of Chinese health administration. A slight concern is that GnRH agonist and GnRH antagonist were used in our study. But, according to literatures [41, 42], different GnRH analogs make small difference in the granulosa cell at the transcriptome level, meaning that they could not undermine our conclusion. In the future, we would like to validate biological functions of the DEGs and reveal the underlying molecular mechanisms.


Collectively, based on the transcriptome sequencing technique, after comparing the gene and protein expression in the granulosa cell and the alteration in the FF components from the samples of POR patients with the mild ovarian stimulation or the conventional COS, we found that the TGF-β2 signaling pathway was correlated with the good quality of oocytes, and played a crucial role in the cross-talking among the three of the granulosa cell, FF and oocyte, which implicates that the mild ovarian stimulation protocol is more beneficial to POR patients.

Availability of data and materials

Sequence data from this article have been deposited with the GEO under accession number GSE191322.



In in Vitro Fertilization


Controlled Ovarian Stimulation


Poor Ovarian Response




Luteinizing Hormone


Gonadotropin-Releasing Hormone agonist


Follicular Fluid


Follicle Stimulation Hormone


Human Chorionic Gonadotropin

E2 :



Body Mass Index


Phosphate Buffered Saline


Room Temperature


Differentially Expressed Genes

P adj :

Adjusted P-values


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Real-Time Quantitative PCR


Glyceraldehyde-3-Phosphate Dehydrogenase


Transforming Growth Factor Beta


Sodium Dodecyl Sulfate


Anti-Müllerian Hormone








Electrochemiluminescence Immunoassay


Bone Morphogenetic Protein-15


Growth Differentiation Factor-9


Enzyme-Linked Immunosorbent Assay


Standard Deviation


2 Pronuclei


Preimplantation Genetic Testing


  1. Biljan MM, Buckett WM, Dean N, Phillips SJ, Tan SL. The outcome of IVF-embryo transfer treatment in patients who develop three follicles or less. Hum Reprod. 2000;15(10):2140–4.

    Article  CAS  Google Scholar 

  2. Schoolcraft WB, Surrey ES, Minjarez DA, Stevens JM, Gardner DK. Management of poor responders: can outcomes be improved with a novel gonadotropin-releasing hormone antagonist/letrozole protocol? Fertil Steril. 2008;89(1):151–6.

    Article  CAS  Google Scholar 

  3. Surrey ES, Schoolcraft WB. Evaluating strategies for improving ovarian response of the poor responder undergoing assisted reproductive techniques. Fertil Steril. 2000;73(4):667–76.

    Article  CAS  Google Scholar 

  4. Pandian Z, McTavish AR, Aucott L, Hamilton MPR, Bhattacharya S. Interventions for 'poor responders' to controlled ovarian hyper stimulation (COH) in in‐vitro fertilisation (IVF). Cochrane Database of Systematic Reviews. 2010;(1):CD004379.

  5. Papathanasiou A, Searle BJ, King NM, Bhattacharya S. Trends in “poor responder” research: lessons learned from RCTs in assisted conception. Hum Reprod Update. 2016;22(3):306–19.

  6. Kreeger PK, Fernandes NN, Woodruff TK, Shea LD. Regulation of mouse follicle development by follicle-stimulating hormone in a three-dimensional in vitro culture system is dependent on follicle stage and dose. Biol Reprod. 2005;73(5):942–50.

    Article  CAS  Google Scholar 

  7. Van Blerkom J, Davis P. Differential effects of repeated ovarian stimulation on cytoplasmic and spindle organization in metaphase II mouse oocytes matured in vivo and in vitro. Hum Reprod. 2001;16(4):757–64.

    Article  Google Scholar 

  8. Hodges CA, Ilagan A, Jennings D, Keri R, Nilson J, Hunt PA. Experimental evidence that changes in oocyte growth influence meiotic chromosome segregation. Hum Reprod. 2002;17(5):1171–80.

    Article  CAS  Google Scholar 

  9. Katz-Jaffe MG, Trounson AO, Cram DS. Chromosome 21 mosaic human preimplantation embryos predominantly arise from diploid conceptions. Fertil Steril. 2005;84(3):634–43.

    Article  Google Scholar 

  10. Thomas FH, Vanderhyden BC. Oocyte-granulosa cell interactions during mouse follicular development: regulation of kit ligand expression and its role in oocyte growth. Reprod Biol Endocrinol. 2006;4:19.

    Article  Google Scholar 

  11. Eppig JJ. Analysis of mouse oogenesis in vitro. Oocyte isolation and the utilization of exogenous energy sources by growing oocytes. J Exp Zool. 1976;198(3):375–82.

    Article  CAS  Google Scholar 

  12. De La Fuente R, Eppig JJ. Transcriptional activity of the mouse oocyte genome: companion granulosa cells modulate transcription and chromatin remodeling. Dev Biol. 2001;229(1):224–36.

    Article  Google Scholar 

  13. Basuino L, Silveira CF Jr. Human follicular fluid and effects on reproduction. JBRA Assist Reprod. 2016;20(1):38–40.

    Article  Google Scholar 

  14. Baskind NE, Orsi NM, Sharma V. Impact of exogenous gonadotropin stimulation on circulatory and follicular fluid cytokine profiles. Int J Reprod Med. 2014;2014:218769.

    Article  Google Scholar 

  15. Papler TB, Bokal EV, Tacer KF, Juvan P, Virant Klun I, Devjak R. Differences in cumulus cells gene expression between modified natural and stimulated in vitro fertilization cycles. J Assist Reprod Genet. 2014;31(1):79–88.

    Article  Google Scholar 

  16. Adriaenssens T, Wathlet S, Segers I, Verheyen G, De Vos A, Van der Elst J, et al. Cumulus cell gene expression is associated with oocyte developmental quality and influenced by patient and treatment characteristics. Hum Reprod. 2010;25(5):1259–70.

    Article  CAS  Google Scholar 

  17. Ferraretti AP, La Marca A, Fauser BC, Tarlatzis B, Nargund G, Gianaroli L, et al. ESHRE consensus on the definition of “poor response” to ovarian stimulation for in vitro fertilization: the Bologna criteria. Hum Reprod. 2011;26(7):1616–24.

  18. Liu X, Li T, Wang B, Xiao X, Liang X, Huang R. Mild stimulation protocol vs conventional controlled ovarian stimulation protocol in poor ovarian response patients: a prospective randomized controlled trial. Arch Gynecol Obstet. 2020;301(5):1331–9.

    Article  CAS  Google Scholar 

  19. Racowsky C, Ohno-Machado L, Kim J, Biggers JD. Is there an advantage in scoring early embryos on more than one day? Hum Reprod. 2009;24(9):2104–13.

    Article  Google Scholar 

  20. Gardner DK, Lane M, Stevens J, Schlenker T, Schoolcraft WB. Blastocyst score affects implantation and pregnancy outcome: towards a single blastocyst transfer. Fertil Steril. 2000;73(6):1155–8.

    Article  CAS  Google Scholar 

  21. Fan Y, Chang Y, Wei L, Chen J, Li J, Goldsmith S, et al. Apoptosis of mural granulosa cells is increased in women with diminished ovarian reserve. J Assist Reprod Genet. 2019;36(6):1225–35.

    Article  Google Scholar 

  22. Kwintkiewicz J, Nishi Y, Yanase T, Giudice LC. Peroxisome proliferator-activated receptor-gamma mediates bisphenol A inhibition of FSH-stimulated IGF-1, aromatase, and estradiol in human granulosa cells. Environ Health Perspect. 2010;118(3):400–6.

    Article  CAS  Google Scholar 

  23. Muhammad F, Yivgi-Ohana N, Shveiky D, Orly J, Alexander S, Laufer N. Levels of steroidogenic acute regulatory protein and mitochondrial membrane potential in granulosa cells of older poor-responder women. Fertil Steril. 2009;91(1):220–5.

    Article  Google Scholar 

  24. Uyar A, Torrealday S, Seli E. Cumulus and granulosa cell markers of oocyte and embryo quality. Fertil Steril. 2013;99(4):979–97.

    Article  CAS  Google Scholar 

  25. Dumesic DA, Meldrum DR, Katz-Jaffe MG, Krisher RL, Schoolcraft WB. Oocyte environment: follicular fluid and cumulus cells are critical for oocyte health. Fertil Steril. 2015;103(2):303–16.

    Article  Google Scholar 

  26. Von Wald T, Monisova Y, Hacker MR, Yoo SW, Penzias AS, Reindollar RR, et al. Age-related variations in follicular apolipoproteins may influence human oocyte maturation and fertility potential. Fertil Steril. 2010;93(7):2354–61.

    Article  Google Scholar 

  27. Liu N, Ma Y, Li R, Jin H, Li M, Huang X, et al. Comparison of follicular fluid amphiregulin and EGF concentrations in patients undergoing IVF with different stimulation protocols. Endocrine. 2012;42(3):708–16.

    Article  CAS  Google Scholar 

  28. Yang WJ, Liu FC, Hsieh JS, Chen CH, Hsiao SY, Lin CS. Matrix metalloproteinase 2 level in human follicular fluid is a reliable marker of human oocyte maturation in in vitro fertilization and intracytoplasmic sperm injection cycles. Reprod Biol Endocrinol. 2015;13:102.

    Article  Google Scholar 

  29. Lu CL, Yan ZQ, Song XL, Xu YY, Zheng XY, Li R, et al. Effect of exogenous gonadotropin on the transcriptome of human granulosa cells and follicular fluid hormone profiles. Reprod Biol Endocrinol. 2019;17(1):49.

    Article  Google Scholar 

  30. Zhang K, Rajput SK, Lee KB, Wang D, Huang J, Folger JK, et al. Evidence supporting a role for SMAD2/3 in bovine early embryonic development: potential implications for embryotropic actions of follistatin. Biol Reprod. 2015;93(4):86.

    Article  Google Scholar 

  31. Ashry M, Lee K, Mondal M, Datta TK, Folger JK, Rajput SK, et al. Expression of TGFbeta superfamily components and other markers of oocyte quality in oocytes selected by brilliant cresyl blue staining: relevance to early embryonic development. Mol Reprod Dev. 2015;82(3):251–64.

    Article  CAS  Google Scholar 

  32. Hardy K, Mora JM, Dunlop C, Carzaniga R, Franks S, Fenwick MA. Nuclear exclusion of SMAD2/3 in granulosa cells is associated with primordial follicle activation in the mouse ovary. J Cell Sci. 2018;131(17):jcs218123.

    Article  Google Scholar 

  33. Bondestam J, Huotari MA, Moren A, Ustinov J, Kaivo-Oja N, Kallio J, et al. cDNA cloning, expression studies and chromosome mapping of human type I serine/threonine kinase receptor ALK7 (ACVR1C). Cytogenet Cell Genet. 2001;95(3–4):157–62.

    Article  CAS  Google Scholar 

  34. Xu YW, Peng YT, Wang B, Zeng YH, Zhuang GL, Zhou CQ. High follicle-stimulating hormone increases aneuploidy in human oocytes matured in vitro. Fertil Steril. 2011;95(1):99–104.

    Article  CAS  Google Scholar 

  35. Gueripel X, Benahmed M, Gougeon A. Sequential gonadotropin treatment of immature mice leads to amplification of transforming growth factor beta action, via upregulation of receptor-type 1, Smad 2 and 4, and downregulation of Smad 6. Biol Reprod. 2004;70(3):640–8.

    Article  CAS  Google Scholar 

  36. Abdel-Aziz AM, Gamal El-Tahawy NF, Salah Abdel Haleem MA, Mohammed MM, Ali AI, Ibrahim YF. Amelioration of testosterone-induced benign prostatic hyperplasia using febuxostat in rats The role of VEGF/TGFbeta and iNOS/COX-2. Eur J Pharmacol. 2020;889:173631.

    Article  CAS  Google Scholar 

  37. Gode F, Gulekli B, Dogan E, Korhan P, Dogan S, Bige O, et al. Influence of follicular fluid GDF9 and BMP15 on embryo quality. Fertil Steril. 2011;95(7):2274–8.

    Article  CAS  Google Scholar 

  38. Wu YT, Tang L, Cai J, Lu XE, Xu J, Zhu XM, et al. High bone morphogenetic protein-15 level in follicular fluid is associated with high quality oocyte and subsequent embryonic development. Hum Reprod. 2007;22(6):1526–31.

    Article  CAS  Google Scholar 

  39. Kobayashi N, Orisaka M, Cao M, Kotsuji F, Leader A, Sakuragi N, et al. Growth differentiation factor-9 mediates follicle-stimulating hormone-thyroid hormone interaction in the regulation of rat preantral follicular development. Endocrinology. 2009;150(12):5566–74.

    Article  CAS  Google Scholar 

  40. Mottershead DG, Sugimura S, Al-Musawi SL, Li JJ, Richani D, White MA, et al. Cumulin, an Oocyte-secreted Heterodimer of the Transforming Growth Factor-beta Family, Is a Potent Activator of Granulosa Cells and Improves Oocyte Quality. J Biol Chem. 2015;290(39):24007–20.

    Article  CAS  Google Scholar 

  41. Devjak R, FonTacer K, Juvan P, VirantKlun I, Rozman D, BokalVrtacnik E. 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.

    Article  CAS  Google Scholar 

  42. Kolibianakis EM, Collins J, Tarlatzis BC, Devroey P, Diedrich K, Griesinger G. Among patients treated for IVF with gonadotrophins and GnRH analogues, is the probability of live birth dependent on the type of analogue used? A systematic review and meta-analysis. Hum Reprod Update. 2006;12(6):651–71.

    Article  CAS  Google Scholar 

Download references


The authors thank the participating infertile women, doctors, laboratory staff, nurses, and clerks at the Reproductive Medicine Research Center of the Sixth Affiliated Hospital of Sun Yat-sen University, and we also thank Dr. Shen Chen for his writing assistance. The author(s) read and approved the final manuscript.


This research was supported by the Natural Science Foundation of Guangdong Province, No.2019A1515011764 and Merck Serono China Research Fund for Fertility Agreement (2015).

Author information

Authors and Affiliations



RH and XYL conceived study conception and study design. RH undertook the recruitment of patients. XPL, PYC, HSM, ZQZ and PS undertook the collection of experimental materials. XPL, CCZ, TBW and JHC conducted the experiments. XPL performed data analyses and drafted the manuscript, and RH revised the manuscript. RH and XPL interpreted the data.

Corresponding author

Correspondence to Rui Huang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Reproductive Medicine and Prenatal Diagnosis, the Sixth Affiliated Hospital of Sun Yat-sen University (2019ZSLYEC-002S), and all participants rendered written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, X., Mai, H., Chen, P. et al. Comparative analyses in transcriptome of human granulosa cells and follicular fluid micro-environment between poor ovarian responders with conventional controlled ovarian or mild ovarian stimulations. Reprod Biol Endocrinol 20, 54 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: