Skip to main content

The effect of growth hormone on the metabolome of follicular fluid in patients with diminished ovarian reserve



Increasing evidence supports that the co-treatment with growth hormone (GH) enhances ovarian response and oocyte quality during controlled ovarian stimulation (COS) in patients with diminished ovarian reserve (DOR). The composition of follicular fluid (FF) plays an essential role in oocyte development and mirrors the communication occurring between the oocyte and follicular microenvironment. However, the effect of GH on the FF metabolome remains unclear.


This prospective observational study recruited DOR patients undergoing in vitro fertilization (IVF) cycles with minimal stimulation protocol for COS. Each patient receiving GH co-treatment was matched to a patient without GH co-treatment by propensity score matching. The FF was collected after isolating oocytes and assayed by gas chromatograph-mass spectrometry (GC-MS) metabolomics. The Pearson correlation was performed to evaluate the relationship between the number of oocytes retrieved and the levels of differential metabolites. The KEGG database was used to map differential metabolites onto various metabolic pathways.


One hundred thirty-four FF metabolites were identified by GC-MS metabolomics. Twenty-four metabolites, including glutathione, itaconic acid and S-adenosylmethionin (SAM) showed significant differences between the GH and control groups (p-value < 0.05 and q-value < 0.1). In addition, the number of oocytes retrieved was significantly higher in the GH group compared to the control group (3 vs 2, p = 0.04) and correlated with the levels of five differential metabolites. Among them, the levels of antioxidant metabolite itaconic acid were upregulated by GH administration, while SAM levels were downregulated.


The co-treatment with GH during COS may improve oocyte development by altering FF metabolite profiles in DOR patients. However, given the downregulation of SAM, a regulator of genomic imprinting, the potential risk of imprinting disturbances should not be neglected.


Diminished ovarian reserve (DOR), also called as poor ovarian reserve (POR), refers to reduced oocyte quantity or quality in women of reproductive age [1]. Women with DOR exhibit decreased response to ovarian stimulation or reduced fecundity compared to women of comparable age [1]. The prevalence of DOR is estimated at 10% among infertile women. Patients with DOR usually seek assisted reproductive technology as fertility treatment. Unfortunately, DOR is associated with ovarian hypo-response, a higher chance of cycle cancellation, reduced numbers of oocytes retrieved and embryos, lower rates of embryo implantation, clinical pregnancy and live birth, and higher rates of miscarriage and aneuploid in in vitro fertilization (IVF) cycles [2,3,4]. Various adjuvant therapies, including growth hormone (GH), testosterone, dehydroepiandrosterone (DHEA), L-arginine, and coenzyme Q10, have been developed to improve IVF outcomes in DOR patients [5,6,7,8,9].

GH, a single-chain peptide secreted by somatotrophic cells within the anterior pituitary in a pulsatile pattern, is involved in regulating protein synthesis, cell proliferation and metabolism [10]. The expression of GH and its receptor (GHR) was detected in human oocytes, granulosa and stromal cells, implying that GH may play a critical role in human reproduction [11]. Thirty years ago, GH was first prescribed to improve the ovarian response in patients who had previously responded sub-optimally to standard ovarian stimulation regimens for IVF [12]. Several recent studies showed that GH treatment improved IVF outcomes of patients with DOR or poor ovarian response, including clinical pregnancy and live birth, probably through the beneficial effect on the number of oocytes retrieved and good-quality embryo formation [13,14,15]. The beneficial role of GH treatment in DOR patients undergoing IVF was further supported by two systematic reviews including randomized controlled trials [16, 17]. GH was demonstrated to improve ovarian response or oocyte quality through upregulating the density of GHR in human oocytes and granulosa cells [18, 19], stimulating the secretion of insulin like growth factor 1 [20], and promoting follicle-stimulating hormone (FSH)-mediated ovarian steroidogenesis [21]. However, the underlying mechanisms have not been fully elucidated.

Follicular fluid (FF) is deemed as a transudate of serum components and the secretions of theca and granulosa cells in ovarian follicles [22], constituted by a variety of molecules, including proteins, steroid hormones, metabolites, and polysaccharides. FF surrounding the oocyte forms an essential element for the follicular microenvironment and plays a crucial role in oocyte development [23]. Links have been established between FF metabolites and the clinical outcomes of patients undergoing IVF [24,25,26]. Notably, FF metabolite profiles of DOR patients are significantly different from that of women with a normal ovarian reserve [27, 28]. These suggested that GH may improve IVF outcomes of DOR patients by changing FF metabolite profiles.

This study aimed to investigate the effect of GH adjuvant treatment on the FF metabolome of DOR patients undergoing IVF and analyze the potential metabolic pathways in which the differential metabolites are involved.

Materials and methods

Study design and participants

This prospective observational study was conducted in the Center for Reproductive Medicine at the Second Affiliated Hospital of Chongqing Medical University. The consecutive IVF cycles with minimal stimulation protocol for controlled ovarian stimulation (COS) between May 2019 and November 2019 were screened for patients with DOR, defined as serum anti-Müllerian hormone (AMH) level below 1.1 ng/ml or antral follicle count (AFC) below seven [29, 30]. We excluded those DOR patients with: 1) body mass index (BMI) ≥24 kg/m2; 2) a history of endometriosis, ovarian surgery, polycystic ovary syndrome, or other endocrine or autoimmune dysfunction; 3) other adjuvant treatment including DHEA and melatonin; 4) previous IVF attempts more than once; or 5) less than 3 months since last IVF attempt. The AFC and basal sex hormone levels were evaluated on days 2-4 of the menstrual cycle. The propensity score matching (PSM) was performed to reduce the potential bias from confounding variables including age and basal FSH levels. Each patient receiving GH co-treatment was matched to a patient without GH co-treatment by PSM. The patients who were not matched were excluded.

IVF treatment

Clomiphene citrate 50 mg per day was administered in conjunction with human menopausal gonadotropin 75 IU per day starting on menstrual cycle day 2 or 3. Patients in the GH group were co-treated with human recombinant GH (Jintropin, JenSci, China) 3 IU per day during COS. When the leading follicle reached a diameter of 18 mm or greater, ovulation triggering was performed with human chorionic gonadotropin (hCG) 10,000 IU, and the number of follicles ≥14 mm was recorded. Oocytes were retrieved 34–36 h later through follicular flushing.

FF collection and metabolite extraction

To avoid blood contamination, FF from only the follicle with the largest diameter was aspirated first and collected during oocyte retrieval. FF samples were immediately centrifuged at 1000 g, 4 °C for 10 min to remove cellular components. The supernatant was collected, divided into 300 μl aliquots, and stored at − 80 °C until sample preparation. Two hundred fifty microliters aliquots of thawed FF were transferred into 1.5 ml Eppendorf tubes. Subsequently, 20 μl internal standard (2,3,3,3-d4-alanine, 10 mM) was added to each aliquot. To precipitate protein from the follicular sample, 730 μl of cold methanol was added to the follicular aliquots, followed by cooling at − 20 °C for 40 min. Then the supernatant was collected by centrifugation at 1500 g for 20 min and dried in a SpeedVac (LABCONCO, #7810041) at 1000 g for 8 hr. The extracted metabolites were stored at − 80 °C freezer prior to derivatization.

Methyl chloroformate (MCF) derivatization and gas chromatography-mass spectrometry (GC-MS)

Dried follicular extracts were resuspended in 200 μl of sodium hydroxide (1 M), and MCF derivatization was added to reduce the boiling points of compounds for GC-MS analysis, according to the methodology published by Smart et al. [31]. All follicular samples were analyzed in a single randomized order, and volatile compounds were separated by An Agilent 7890B GC system using a ZB-1701 GC capillary column (30 m × 250 μm id × 0.15 μm with a 5 m guard column, Phenomenex) and analyzed m/z of ions by an Agilent 5977As mass spectrometer (Agilent, California, USA) with electron impact ionization via electron emission at 70 eV. The GC temperature ramps and MS settings were operated according to the protocol published in Han et al. [32].

Data extraction and normalization

The GC peaks were first deconvoluted by Automated Mass Spectral Deconvolution & Identification System software (AMIDS). The compounds were identified by Ion fragmentation patterns and GC retention time to our in-house mass spectral library built by chemical standards. The MassOmics XCMS R-based software was implemented to extrapolate the relative concentration of the metabolites through the peak height of the most abundant fragmented ion mass [33]. To achieve the reproducibility robustness along with minimizing sample preparation and instrumental variabilities, the relative concentration of the identified metabolites were normalized in the sequence of internal standards (D4-alanine, D5-phenylalanine, or D2-tyrosine), batch correction by quality control of pooled samples, and total ion concentration (TIC) of the follicular metabolome. Lastly, negative blanks were used to eliminate carryover contaminations.

Statistical analysis

Categorical data were expressed as frequency or percentage, and statistical comparison between groups was carried out using Pearson’s chi-squared test or Fisher’s exact test. The Shapiro–Wilk test was used to test the normal distribution of continuous variables, which were expressed as mean ± standard deviation if normally distributed and as median (25th percentile, 75th percentile) if not normally distributed. Statistical comparison was performed using the Student t-test or Mann–Whitney U test for continuous variables, where appropriate. Prior to statistical analysis of the FF metabolome, metabolite concentration was transformed by log2 scale and Pareto scaling to establish Gaussian distribution for this dataset. Partial least squares discriminant analysis (PLS-DA) and model validation were performed by MetaboAnalyst 5.0 ( Since clinical characteristics associated with IVF outcomes were matched between the GH and control groups, Student t-test and false discovery rate were implemented to calculate the significance of FF metabolites between two groups using R software. Only both two-tailed p-values and q-values less than 0.05 and 0.1 respectively were considered statistically significant. The areas under the receiver operating characteristic (ROC) curve were performed using pROC R package [34]. The forest plot displaying the Pearson correlation between the number of oocytes retrieved and the levels of differential metabolites was illustrated using the ggplot2 R package [35]. The Sankey diagram, which links differential metabolites into their KEGG metabolic pathways, was performed by the Online website The metabolic network was illustrated based on the KEGG global metabolism map using MetaboAbalyst 5.0.


Clinical characteristics of participants

A total of 64 DOR patients, comprising 32 in the GH group and 32 in the control group, were included in this study. The clinical characteristics of the participants are presented in Table 1. The baseline characteristics were comparable between the GH and control groups, including age, BMI, duration of infertility, previous conception, previous IVF attempt, AMH, AFC, in addition to the levels of basal FSH, luteinizing hormone (LH) and estradiol (E2). There was no significant difference in the days of stimulation, FSH dose, or LH dose between both groups. It was noteworthy that the number of follicles above 14 mm on the day of hCG was significantly higher in the GH group compared to the control group (3 vs 2, p = 0.03; Fig. 1A). More importantly, this difference was also found for the number of oocytes retrieved between the GH and control groups (3 vs 2, p = 0.04; Fig. 1B). Moreover, there were trends toward higher E2 levels on the day of hCG (1380.8 ± 638.2 vs 1160.3 ± 582.6, p = 0.15) and normal fertilization rate (82.4% vs 69.1%, p = 0.05) in the GH group compared with the control group.

Table 1 Clinical characteristics of participants
Fig. 1
figure 1

The numbers of follicles ≥14 mm on the day of hCG (A) and oocytes retrieved (B)

The difference of FF metabolite profiles between the GH and control groups

Based on all 134 metabolites identified from FF, the PLS-DA displayed a distinct separation between the GH and control groups via principal components PC1, PC2 and PC3, representing 7.7, 10.5, and 5.7% of the variance respectively (Fig. 2). The leave-one-out validation model showed the best PLS-DA performance using three accumulative PC (Accuracy = 0.83, R2 = 0.88, Q2 = 0.54). Among identified metabolites, 24 of them appeared to be significantly different in concentration between two groups with p-value and q-value less than 0.05 and 0.1, respectively (Fig. 3 and Table S1). In FF from patients with GH administration, the levels of itaconic acid, glutathione, cis-aconitic acid, N-alpha-acetyllysine, stearic acid, tridecane, and the majority of organic acids were significantly elevated, while those of S-adenosylmethionine (SAM), 2-oxobutyric acid, citramalic acid, and butylated hydroxytoluene were significantly reduced compared with controls. Furthermore, most of the unsaturated-chain fatty acids displayed reduced levels in the GH group, including linolelaidic acid, 9-heptadecenoic acid, and palmitelaidic acid.

Fig. 2
figure 2

PLS-DA of the FF metabolome between the GH and control groups. The right bar graphs evaluate the prediction model performance via leave-one-out validation. R2 indicates how well the model represents the data. Q2 indicates how reproducible is the PLS-DA model. PLS-DA partial least squares discriminant analysis, FF follicular fluid, GH growth hormone

Fig. 3
figure 3

The differences in the FF metabolome between the GH and control groups. The relative concentrations of FF metabolites are illustrated via a log2 scale. Positive values represent higher metabolite levels in the dividend group (GH group) than the divisor group (control group), whereas negative values represent lower metabolite levels in the dividend group than the divisor group. Only the metabolites with a p-value less than 0.05 (Student t-test) and a q-value less than 0.01 (false discovery rate) are displayed. FF follicular fluid, GH growth hormone, TCA tricarboxylic acid

The correlation between the number of oocytes retrieved and the levels of FF metabolites

The number of oocytes retrieved is positively correlated with the levels of conjugated linoleic acid, itaconic acid, and tridecane, while negatively correlated with the levels of D − norleucine and SAM (Fig. 4).

Fig. 4
figure 4

The correlation between the number of oocytes retrieved and the levels of differential metabolites. Statistical analyses were conducted using Pearson correlation. TCA tricarboxylic acid, CI confidence interval

The metabolic pathway enrichment analysis of differential metabolites

As shown in Fig. 5, these differential metabolites were mapped into their corresponding metabolic pathways in order to investigate the potential biological processes altered by GH. FF metabolites, including conjugated linoleic acid, cis-aconitic acid, stearic acid and glutathione, were involved in various downregulated center carbon metabolism. Especially, glutathione was branched into the most diverse metabolic pathways, namely ferroptosis, glutathione metabolism, thyroid hormone synthesis, ATP-binding cassette (ABC) transporters, biosynthesis of cofactors, as well as cysteine and methionine metabolism. And the other differential metabolites, 2-oxobutyric acid, lysine, SAM, and linolelaidic acid were involved in the upregulation of amino acid or fatty acid metabolism. In addition, SAM is also linked to DNA methylation. The significant metabolic pathways were connected to their common metabolites and constructed in silico into a metabolic network using KEGG global metabolic framework, as displayed in Fig. 6. The metabolic network demonstrated that some differential metabolites could be linked with the shortest distance by the following metabolic pathways: glutathione metabolism, cysteine and methionine metabolism, linoleic acid metabolism, biosynthesis of unsaturated fatty acid, tricarboxylic acid (TCA) cycle, glyoxylate metabolism, and lysine degradation.

Fig. 5
figure 5

A Sankey diagram displays how differential metabolites participated in upregulated or downregulated KEGG metabolic pathways. ABC ATP-binding cassette, TCA tricarboxylic acid

Fig. 6
figure 6

The KEGG metabolic network in which FF metabolites of DOR patients were involved. The metabolic network was based on KEGG global metabolic pathways. The red dots and red letters are the metabolites that were statistically significant between the GH and control groups. The yellow dots and black letters are identified metabolites that displayed no significance. The red lines are metabolic pathways that connect differential metabolites with the shortest distance. FF follicular fluid, DOR diminished ovarian reserve, GH growth hormone


This research was the first to employ metabolomics approaches to analyze the effect of GH administration on the FF metabolome in DOR patients undergoing IVF. Our results showed that the co-treatment with GH during COS is involved in altering FF metabolite profiles. Twenty-four differential metabolites were identified between the GH and control groups. Among them, the concentrations of antioxidant metabolites itaconic acid and glutathione were increased by GH administration, while SAM concentration was reduced. In addition, the number of oocytes retrieved increased with GH administration and was correlated with the levels of five differential metabolites, including itaconic acid and SAM. These findings may reveal the potential therapeutic mechanism of GH in improving oocyte development of DOR patients.

All essential intermediates of the TCA cycle were detected in FF. Interestingly, a significantly higher level of itaconic acid was observed in the GH group. Itaconic acid and itaconate generate from the decarboxylation of cis-aconitate, a TCA cycle intermediate. The derivative of itaconic acid, four-octyl itaconate, has been demonstrated to attenuate H2O2-induced reactive oxygen species (ROS) production, lipid peroxidation and DNA damage, in addition to neuronal cell death and apoptosis [36]. Another study found that itaconate slowed down TCA cycle metabolism and reduced ROS levels to improve brain function by inhibiting succinate dehydrogenase [37]. Oxidative stress due to ROS accumulation in the follicular microenvironment plays a critical role in ovarian aging or DOR development [38,39,40]. Furthermore, reduced ovarian response to COS is also associated with increased oxidative stress in the follicular microenvironment [41]. Antioxidants, such as coenzyme Q10 and melatonin, improve ovarian response and embryo quality in DOR patients undergoing IVF [42, 43]. Recently, GH was reported to decrease ROS levels in granulosa cells and increase the number of high-quality embryos in patients with poor ovarian response, but the underlying mechanism is unclear [14]. In the present study, the level of itaconic acid was positively correlated with the number of oocytes retrieved. These findings suggest that GH treatment improves ovarian response probably by reducing ROS generation through elevating the levels of itaconic acid in the follicular microenvironment.

The glutathione levels in the GH group were also higher than those in the control group. Likewise, GH treatment led to increased glutathione levels in rat ovary tissues [44]. Glutathione is a powerful antioxidant removing ROS and protects oocytes against oxidative damage and exhibits a beneficial effect on the quality of ovine and bovine oocytes [45,46,47]. Most importantly, higher glutathione levels in granulosa cells and FF have been found to be associated with increased fertilization potential of oocytes in IVF patients [48, 49]. Moreover, the KEGG metabolic pathway analysis revealed a downregulation of ferroptosis and glutathione catabolism, the latter of which might be owing to the suppression of glutathione degradation enzyme by GH administration [50]. Ferroptosis, an iron-dependent cell death, can be initiated by glutathione depletion [51]. Nowadays, there is growing evidence that ferroptosis may also be a key driver of pathology in age-related diseases [52]. Although little has been reported on the direct association between shortlisted metabolic pathways (glutathione metabolism, ferroptosis, and fatty acid metabolism; Fig. 5) and patients with DOR, Liang et al. demonstrated that oxylipins metabolism was significantly altered in FF of DOR patients [28]. Oxylipins are a class of lipid metabolites that derive from the oxidation of unsaturated fatty acids, and their production can be affected by ROS and fatty acid levels, which are related to our shortlisted metabolic pathways. Taken together, GH may rejuvenate the oocytes of DOR patients by raising glutathione levels and then attenuating oxidative stress in the follicular microenvironment.

Last but not least, there was a significant decrease of SAM in the GH group compared to the control group. This finding is supported by the fact that daily injections of GH inhibited the activity and mRNA transcription of SAM synthetase [53]. Besides a major methyl donor for the methylation of DNA, RNA and histone, SAM is one of the glutathione precursors [54], and our KEGG metabolic network also reveals the connection between them. Recently it was proved that the enhanced biosynthesis of SAM in the germline of Drosophila ovaries led to aging-related defects in oogenesis [55]. Although it has not been evaluated in human ovaries, the SAM biosynthesis is regarded as evolutionarily conserved. In addition, a variety of evidence suggested a strong association between methylation and cell aging [56]. Compared to women with normal ovarian reserve, a distinctive DNA methylation profile was found in mural granulosa cells of DOR patients [57]. In our study, the SAM levels in FF were negatively correlated with the number of oocytes retrieved. Therefore, GH administration may improve ovarian response by reducing SAM biosynthesis in follicles. In concert with the decreased SAM levels induced by GH in this study, GH was reported to result in age- and sex-dependent DNA hypomethylation [58]. Given these findings, it is necessary to notice the possibility of aberrant genomic imprinting, which gives rise to the onset of imprinting disorders in the offspring.


The findings in this study suggest that the co-treatment with GH during COS alters FF metabolite profiles and in turn increases the number of oocytes retrieved in DOR patients. The novel data reported here show the elevating effect of GH on the antioxidant metabolites itaconic acid and glutathione. In addition, the levels of SAM, a regulator of genomic imprinting, were downregulated by GH treatment. In future work, the antioxidant capability of itaconic acid/glutathione and SAM synthetase activity regulated by GH should be validated. Combined with previous findings, GH has been proved to improve IVF outcomes, but we should be aware of the potential risk of imprinting disturbances. And large multi-center cohort studies are warranted to evaluate this possible harmful effect of GH on offsprings.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.



ATP-binding cassette


Antral follicle count


Anti-Müllerian hormone


Automated Mass Spectral Deconvolution & Identification System software


Area under the receiver operating characteristic curve


Body mass index


Controlled ovarian stimulation




Diminished ovarian reserve




Follicular fluid


Follicle-stimulating hormone


Gas chromatography-mass spectrometry


Growth hormone


Growth hormone receptor


Human chorionic gonadotropin


In vitro fertilization


Luteinizing hormone


Methyl chloroformate


Receiver operating characteristic


Partial least squares discriminant analysis


Poor ovarian reserve


Principal components


Reactive oxygen species




Tricarboxylic acid


Total ion concentration


  1. Practice Committee of the American Society for Reproductive M. Testing and interpreting measures of ovarian reserve: a committee opinion. Fertil Steril. 2020;114:6.

    Google Scholar 

  2. Chang Y, Li J, Li X, Liu H, Liang X. Egg quality and pregnancy outcome in young infertile women with diminished ovarian reserve. Med Sci Monit. 2018;24:7279–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Devine K, Mumford SL, Wu M, DeCherney AH, Hill MJ, Propst A. Diminished ovarian reserve in the United States assisted reproductive technology population: diagnostic trends among 181,536 cycles from the Society for Assisted Reproductive Technology Clinic Outcomes Reporting System. Fertil Steril. 2015;104:3.

    Article  Google Scholar 

  4. Jaswa EG, McCulloch CE, Simbulan R, Cedars MI, Rosen MP. Diminished ovarian reserve is associated with reduced euploid rates via preimplantation genetic testing for aneuploidy independently from age: evidence for concomitant reduction in oocyte quality with quantity. Fertil Steril. 2021;115:4.

    Article  Google Scholar 

  5. Li XL, Wang L, Lv F, Huang XM, Wang LP, Pan Y, et al. The influence of different growth hormone addition protocols to poor ovarian responders on clinical outcomes in controlled ovary stimulation cycles: a systematic review and meta-analysis. Medicine (Baltimore). 2017;96:12.

    Google Scholar 

  6. Nagels HE, Rishworth JR, Siristatidis CS, Kroon B. Androgens (dehydroepiandrosterone or testosterone) for women undergoing assisted reproduction. Cochrane Database Syst Rev. 2015;(11):CD009749.

  7. Li J, Yuan H, Chen Y, Wu H, Wu H, Li L. A meta-analysis of dehydroepiandrosterone supplementation among women with diminished ovarian reserve undergoing in vitro fertilization or intracytoplasmic sperm injection. Int J Gynaecol Obstet. 2015;131:3.

    Article  Google Scholar 

  8. Battaglia C, Salvatori M, Maxia N, Petraglia F, Facchinetti F, Volpe A. Adjuvant L-arginine treatment for in-vitro fertilization in poor responder patients. Hum Reprod. 1999;14:7.

    Article  Google Scholar 

  9. Giannubilo SR, Orlando P, Silvestri S, Cirilli I, Marcheggiani F, Ciavattini A, et al. CoQ10 supplementation in patients undergoing IVF-ET: the relationship with follicular fluid content and oocyte maturity. Antioxidants (Basel). 2018;7:10.

    Google Scholar 

  10. Donato J Jr, Wasinski F, Furigo IC, Metzger M, Frazao R. Central regulation of metabolism by growth hormone. Cells. 2021;10:1.

    Article  Google Scholar 

  11. Abir R, Garor R, Felz C, Nitke S, Krissi H, Fisch B. Growth hormone and its receptor in human ovaries from fetuses and adults. Fertil Steril. 2008;90(4 Suppl):1333–9.

    Article  CAS  PubMed  Google Scholar 

  12. Owen EJ, West C, Mason BA, Jacobs HS. Co-treatment with growth hormone of sub-optimal responders in IVF-ET. Hum Reprod. 1991;6:4.

    Article  Google Scholar 

  13. Cai MH, Gao LZ, Liang XY, Fang C, Wu YQ, Yang X. The effect of growth hormone on the clinical outcomes of poor ovarian reserve patients undergoing in vitro fertilization/intracytoplasmic sperm injection treatment: a retrospective study based on POSEIDON criteria. Front Endocrinol (Lausanne). 2019;10:775.

    Article  PubMed  Google Scholar 

  14. Gong Y, Zhang K, Xiong D, Wei J, Tan H, Qin S. Growth hormone alleviates oxidative stress and improves the IVF outcomes of poor ovarian responders: a randomized controlled trial. Reprod Biol Endocrinol. 2020;18:1.

    Article  Google Scholar 

  15. Dogan S, Cicek OSY, Demir M, Yalcinkaya L, Sertel E. The effect of growth hormone adjuvant therapy on assisted reproductive technologies outcomes in patients with diminished ovarian reserve or poor ovarian response. J Gynecol Obstet Hum Reprod. 2021;50:2.

    Article  Google Scholar 

  16. Liu FT, Hu KL, Li R. Effects of growth hormone supplementation on poor ovarian responders in assisted reproductive technology: a systematic review and meta-analysis. Reprod Sci. 2021;28:4.

    Google Scholar 

  17. Zhang Y, Zhang C, Shu J, Guo J, Chang HM, Leung PCK, et al. Adjuvant treatment strategies in ovarian stimulation for poor responders undergoing IVF: a systematic review and network meta-analysis. Hum Reprod Update. 2020;26:2.

    Article  Google Scholar 

  18. Weall BM, Al-Samerria S, Conceicao J, Yovich JL, Almahbobi G. A direct action for GH in improvement of oocyte quality in poor-responder patients. Reproduction. 2015;149:2.

    Article  Google Scholar 

  19. Regan SLP, Knight PG, Yovich JL, Arfuso F, Dharmarajan A. Growth hormone during in vitro fertilization in older women modulates the density of receptors in granulosa cells, with improved pregnancy outcomes. Fertil Steril. 2018;110:7.

    Article  Google Scholar 

  20. Dosouto C, Calaf J, Polo A, Haahr T, Humaidan P. Growth hormone and reproduction: lessons learned from animal models and clinical trials. Front Endocrinol (Lausanne). 2019;10:404.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Xu YM, Hao GM, Gao BL. Application of growth hormone in in vitro fertilization. Front Endocrinol (Lausanne). 2019;10:502.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Fuentes A, Sequeira K, Tapia-Pizarro A, Munoz A, Salinas A, Cespedes P, et al. Androgens profile in blood serum and follicular fluid of women with poor ovarian response during controlled ovarian stimulation reveals differences amongst POSEIDON stratification groups: a pilot study. Front Endocrinol (Lausanne). 2019;10:458.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Revelli A, Delle Piane L, Casano S, Molinari E, Massobrio M, Rinaudo P. Follicular fluid content and oocyte quality: from single biochemical markers to metabolomics. Reprod Biol Endocrinol. 2009;7:40.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wallace M, Cottell E, Gibney MJ, McAuliffe FM, Wingfield M, Brennan L. An investigation into the relationship between the metabolic profile of follicular fluid, oocyte developmental potential, and implantation outcome. Fertil Steril. 2012;97:5.

    Article  Google Scholar 

  25. Shehadeh A, Bruck-Haimson R, Saidemberg D, Zacharia A, Herzberg S, Ben-Meir A, et al. A shift in follicular fluid from triacylglycerols to membrane lipids is associated with positive pregnancy outcome. FASEB J. 2019;33:9.

    Article  Google Scholar 

  26. Nagy RA, Homminga I, Jia C, Liu F, Anderson JLC, Hoek A, et al. Trimethylamine-N-oxide is present in human follicular fluid and is a negative predictor of embryo quality. Hum Reprod. 2020;35:1.

    Article  Google Scholar 

  27. de la Barca JMC, Boueilh T, Simard G, Boucret L, Ferré-L'Hotellier V, Tessier L, et al. Targeted metabolomics reveals reduced levels of polyunsaturated choline plasmalogens and a smaller dimethylarginine/arginine ratio in the follicular fluid of patients with a diminished ovarian reserve. Hum Reprod. 2017;32:11.

    Google Scholar 

  28. Liang C, Zhang X, Qi C, Hu H, Zhang Q, Zhu X, et al. UHPLC-MS-MS analysis of oxylipins metabolomics components of follicular fluid in infertile individuals with diminished ovarian reserve. Reprod Biol Endocrinol. 2021;19:1.

    Article  Google Scholar 

  29. Busnelli A, Somigliana E, Cirillo F, Levi-Setti PE. Is diminished ovarian reserve a risk factor for miscarriage? Results of a systematic review and meta-analysis. Hum Reprod Update. 2021;27:6.

    Article  Google Scholar 

  30. Cohen J, Chabbert-Buffet N, Darai E. Diminished ovarian reserve, premature ovarian failure, poor ovarian responder--a plea for universal definitions. J Assist Reprod Genet. 2015;32:12.

    Article  Google Scholar 

  31. Smart KF, Aggio RBM, Van Houtte JR, Villas-Bôas SG. Analytical platform for metabolome analysis of microbial cells using methyl chloroformate derivatization followed by gas chromatography-mass spectrometry. Nat Protoc. 2010;5:10.

    Article  Google Scholar 

  32. Han TL, Cannon RD, Gallo SM, Villas-Boas SG. A metabolomic study of the effect of Candida albicans glutamate dehydrogenase deletion on growth and morphogenesis. NPJ Biofilms Microbiomes. 2019;5(1):13.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Guo G, McKenzie EJ, Jones B, Zarate E, Seymour J, Baker PN, et al. MassOmics: an R package of a cross-platform data processing pipeline for large-scale GC-MS untargeted metabolomics datasets. Zendo. 2021.

  34. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.

    Book  Google Scholar 

  36. Liu H, Feng Y, Xu M, Yang J, Wang Z, Di G. Four-octyl itaconate activates Keap1-Nrf2 signaling to protect neuronal cells from hydrogen peroxide. Cell Commun Signal. 2018;16:1.

    Article  Google Scholar 

  37. Cordes T, Lucas A, Divakaruni AS, Murphy AN, Cabrales P, Metallo CM. Itaconate modulates tricarboxylic acid and redox metabolism to mitigate reperfusion injury. Mol Metab. 2020;32:122–35.

    Article  CAS  PubMed  Google Scholar 

  38. Wang S, Zheng Y, Li J, Yu Y, Zhang W, Song M, et al. Single-cell transcriptomic atlas of primate ovarian aging. Cell. 2020;180:3.

    Article  Google Scholar 

  39. Park SU, Walsh L, Berkowitz KM. Mechanisms of ovarian aging. Reproduction. 2021;162:2.

    Article  Google Scholar 

  40. Li XH, Wang HP, Tan J, Wu YD, Yang M, Mao CZ, et al. Loss of pigment epithelium-derived factor leads to ovarian oxidative damage accompanied by diminished ovarian reserve in mice. Life Sci. 2019;216:129–39.

    Article  CAS  PubMed  Google Scholar 

  41. Uppangala S, Fernandes G, Salian SR, Kumar P, Talevi R, Kalthur G, et al. Reduced ovarian response to controlled ovarian stimulation is associated with increased oxidative stress in the follicular environment. Reprod Biol. 2020;20:3.

    Article  Google Scholar 

  42. Xu Y, Nisenblat V, Lu C, Li R, Qiao J, Zhen X, et al. Pretreatment with coenzyme Q10 improves ovarian response and embryo quality in low-prognosis young women with decreased ovarian reserve: a randomized controlled trial. Reprod Biol Endocrinol. 2018;16:1.

    Article  Google Scholar 

  43. Jahromi BN, Sadeghi S, Alipour S, Parsanezhad ME, Alamdarloo SM. Effect of melatonin on the outcome of assisted reproductive technique cycles in women with diminished ovarian reserve: a double-blinded randomized clinical trial. Iran J Med Sci. 2017;42:1.

    Google Scholar 

  44. Yigiter M, Halici Z, Odabasoglu F, Keles ON, Atalay F, Unal B, et al. Growth hormone reduces tissue damage in rat ovaries subjected to torsion and detorsion: biochemical and histopathologic evaluation. Eur J Obstet Gynecol Reprod Biol. 2011;157:1.

    Article  Google Scholar 

  45. Luberda Z. The role of glutathione in mammalian gametes. Reprod Biol. 2005;5:1.

    Google Scholar 

  46. Li HJ, Sutton-McDowall ML, Wang X, Sugimura S, Thompson JG, Gilchrist RB. Extending prematuration with cAMP modulators enhances the cumulus contribution to oocyte antioxidant defence and oocyte quality via gap junctions. Hum Reprod. 2016;31:4.

    Article  Google Scholar 

  47. He J, Yao G, He Q, Zhang T, Fan H, Bai Y, et al. Theaflavin 3, 3′-digallate delays ovarian aging by improving oocyte quality and regulating granulosa cell function. Oxidative Med Cell Longev. 2021;2021:7064179.

    Article  Google Scholar 

  48. Tola EN, Mungan MT, Uguz AC, Naziroglu M. Intracellular Ca2+ and antioxidant values induced positive effect on fertilisation ratio and oocyte quality of granulosa cells in patients undergoing in vitro fertilisation. Reprod Fertil Dev. 2013;25:5.

    Article  Google Scholar 

  49. Nishihara T, Matsumoto K, Hosoi Y, Morimoto Y. Evaluation of antioxidant status and oxidative stress markers in follicular fluid for human in vitro fertilization outcome. Reprod Med Biol. 2018;17:4.

    Article  Google Scholar 

  50. Brown-Borg HM, Rakoczy SG, Uthus EO. Growth hormone alters components of the glutathione metabolic pathway in Ames dwarf mice. Ann N Y Acad Sci. 2004;1019:317–20.

    Article  CAS  PubMed  Google Scholar 

  51. Yang WS, SriRamaratnam R, Welsch ME, Shimada K, Skouta R, Viswanathan VS, et al. Regulation of ferroptotic cancer cell death by GPX4. Cell. 2014;156:1–2.

    Article  Google Scholar 

  52. Stockwell BR, Friedmann Angeli JP, Bayir H, Bush AI, Conrad M, Dixon SJ, et al. Ferroptosis: a regulated cell death nexus linking metabolism, redox biology, and disease. Cell. 2017;171:2.

    Article  Google Scholar 

  53. Oscarsson J, Gardmo C, Eden S, Mode A. Pulsatile growth hormone secretion decreases S-adenosylmethionine synthetase in rat liver. Am J Physiol Endocrinol Metab. 2001;280:2.

    Article  Google Scholar 

  54. Lu SC. Regulation of glutathione synthesis. Mol Asp Med. 2009;30:1–2.

    Article  Google Scholar 

  55. Hayashi Y, Kashio S, Murotomi K, Hino S, Kang W, Miyado K, et al. Biosynthesis of S-adenosyl-methionine enhances aging-related defects in Drosophila oogenesis. Sci Rep. 2022;12:1.

    Article  Google Scholar 

  56. Sen P, Shah PP, Nativio R, Berger SL. Epigenetic mechanisms of longevity and aging. Cell. 2016;166:4.

    Article  Google Scholar 

  57. Olsen KW, Castillo-Fernandez J, Chan AC, la Cour Freiesleben N, Zedeler A, Bungum M, et al. Identification of a unique epigenetic profile in women with diminished ovarian reserve. Fertil Steril. 2021;115:3.

    Article  Google Scholar 

  58. Takasugi M, Hayakawa K, Arai D, Shiota K. Age- and sex-dependent DNA hypomethylation controlled by growth hormone in mouse liver. Mech Ageing Dev. 2013;134:7–8.

    Article  Google Scholar 

Download references


We thank our colleagues at the Center for Reproductive Medicine for their assistance in sample collection. We appreciate all patients who participated in this study.


This study was funded by the National Natural Science Foundation of China (81871185, 82001554, 82101710), Medical Research Foundation of Sichuan Medical Association (S18071), Chongqing Municipal Education Commission (KJZD-K202100407), Natural Science Foundation of Chongqing (cstc2021jcyj-msxmX0213, cstc2021jcyj-msxmX0097), the Kuanren Talents Program of the Second Affiliated Hospital of Chongqing Medical University and the Program for Youth Innovation in Future Medicine, Chongqing Medical University.

Author information

Authors and Affiliations



FH, FW, TH, and LH designed the study. HZ and HC included patients. SZ and FH collected data. FH, ZY, CS and HY collected samples. YY performed the experiments. TH, SZ, and SG analyzed data. FH, SZ and FW wrote the manuscript. TH and LH revised the manuscript. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Shen Zhang, Lina Hu or Ting-li Han.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the Second Affiliated Hospital of Chongqing Medical University (No. 2018-79). All methods were carried out in accordance with the Guiding Principles for Research with Human Samples issued by the Second Affiliated Hospital of Chongqing Medical University. All participants signed their written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

The area under the receiver operating characteristic curve (AUC), p-value, q-value and peak intensity of identified metabolites.

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

He, F., Wang, F., Yang, Y. et al. The effect of growth hormone on the metabolome of follicular fluid in patients with diminished ovarian reserve. Reprod Biol Endocrinol 21, 21 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: