Skip to main content

Dehydroepiandrosterone supplementation and the impact of follicular fluid metabolome and cytokinome profiles in poor ovarian responders



Poor ovarian responders (POR) are women undergoing in-vitro fertilization who respond poorly to ovarian stimulation, resulting in the retrieval of lower number of oocytes, and subsequently lower pregnancy rates. The follicular fluid (FF) provides a crucial microenvironment for the proper development of follicles and oocytes through tightly controlled metabolism and cell signaling. Androgens such as dehydroepiandrosterone (DHEA) have been proposed to alter the POR follicular microenvironment, but the impact DHEA imposes on the FF metabolome and cytokine profiles is unknown. Therefore, the objective of this study is to profile and identify metabolomic changes in the FF with DHEA supplementation in POR patients.


FF samples collected from 52 POR patients who underwent IVF with DHEA supplementation (DHEA +) and without (DHEA-; controls) were analyzed using untargeted liquid chromatography-tandem mass spectrometry (LC–MS/MS) metabolomics and a large-scale multiplex suspension immunoassay covering 65 cytokines, chemokines and growth factors. Multivariate statistical modelling by partial least squares-discriminant regression (PLSR) analysis was performed for revealing metabolome-scale differences. Further, differential metabolite analysis between the two groups was performed by PLSR β-coefficient regression analysis and Student’s t-test.


Untargeted metabolomics identified 118 FF metabolites of diverse chemistries and concentrations which spanned three orders of magnitude. They include metabolic products highly associated with ovarian function – amino acids for regulating pH and osmolarity, lipids such fatty acids and cholesterols for oocyte maturation, and glucocorticoids for ovarian steroidogenesis. Four metabolites, namely, glycerophosphocholine, linoleic acid, progesterone, and valine were significantly lower in DHEA + relative to DHEA- (p < 0.05–0.005). The area under the curves of progesterone glycerophosphocholine, linoleic acid and valine are 0.711, 0.730, 0.785 and 0.818 (p < 0.05–0.01). In DHEA + patients, progesterone positively correlated with IGF-1 (Pearson r: 0.6757, p < 0.01); glycerophosphocholine negatively correlated with AMH (Pearson r: -0.5815; p < 0.05); linoleic acid correlated with estradiol and IGF-1 (Pearson r: 0.7016 and 0.8203, respectively; p < 0.01 for both).

In DHEA- patients, valine negatively correlated with serum-free testosterone (Pearson r: -0.8774; p < 0.0001). Using the large-scale immunoassay of 45 cytokines, we observed significantly lower MCP1, IFNγ, LIF and VEGF-D levels in DHEA + relative to DHEA.


In POR patients, DHEA supplementation altered the FF metabolome and cytokine profile. The identified four FF metabolites that significantly changed with DHEA may provide information for titrating and monitoring individual DHEA supplementation.


Poor ovarian responders (POR) are a sub-group of infertile women that account for 9–26% of in vitro fertilization (IVF) indications [1, 2]. In patients designated as “poor responders,” so-called due to poor response to ovarian stimulation given during IVF workup, the limited number of obtained oocytes remains the major problem in optimizing the live birth rates [3]. As a result of a lower number of oocytes retrieved, there are fewer embryos to select and transfer, and subsequently these patients have lower pregnancy rates per transfer and lower cumulative pregnancy rates per started cycle compared with normal responders. In PORs the mechanism of ovarian insufficiency can be multifactorial with causes such as ovarian surgery especially in case of endometrioma [4, 5], uterine artery embolization for the treatment of uterine leiomyoma [6, 7], genetic defects, chemotherapy, radiotherapy, autoimmune disorders, single ovary, chronic smoking [8, 9], or linked to diseases such as diabetes mellitus Type I [10]. However, in most cases, follicular depletion plausibly reflecting premature ovarian aging [11], clinically translates into a reduction of implantation rates, an increase of early pregnancy loss, and disappointingly low IVF success [12, 13].

In each menstrual cycle, human ovaries produce a single dominant follicle. Growth of the dominant follicle encompasses enlargement of the oocyte, replication of follicular cells, and formation and expansion of a fluid-filled follicular antrum or cavity, providing a specialized microenvironment for the development of oocytes. Follicular fluid (FF) that fills the antrum cavity is derived from the surrounding theca capillaries, abundant and easily accessible during IVF procedures due to ample volume being produced during follicle maturation [14]. FF are rich in metabolites, notably hormones, amino acids and lipids that are critical for oocyte growth and development, which determines subsequent potential to achieve fertilization and embryo development [15]. As such, constituents of the FF surrounding the oocyte provides a unique biochemical window to the growth and differentiation of the oocyte [16]. To unravel the biochemical composition of human FF and its impact on oocyte development, metabolomic analyses using gas chromatography–mass spectrometric and proton nuclear magnetic resonance metabolomic analyses have been conducted [17,18,19,20,21,22], as were proteomic analyses [23,24,25,26,27,28]. These studies mainly report the FF profiles of IVF patients, whereas the FF metabolome of POR remains poorly characterized. Furthermore, the effect of dehydroepiandrosterone (DHEA) on the FF metabolome has not been previously studied.

DHEA is a steroid produced in the adrenal cortex and the ovarian theca cells in women that is converted into more active forms of androgens such as testosterone [5]. It has been suggested that DHEA supplementation may increase the number of available follicles in PORs through an increased serum level of insulin-like growth factor, increased follicular response to follicle stimulating hormone (FSH), shifts to aerobic metabolism [29] and improved quality of oocytes [30]. However, the efficacy of DHEA pre-treatment has been controversial, with partial to reasonable clinical evidence being observed [31,32,33,34,35,36,37]. Based on these findings, we conducted a study to evaluate whether the FF metabolome differed in POR patients treated with DHEA or not, and whether the FF metabolome may be predictive of IVF outcome. Furthermore, DHEA has immunoregulatory functions [38], and a large-scale study of FF cytokines was conducted in parallel to reveal DHEA immunomodulatory targets.

Material and methods

Ethical approval and study population

The local Institutional Review Board approved the study (CIRB/2011/404/D) and written informed consent was obtained from each participant. A prospective case–control study was conducted to evaluate the metabolic and cytokine effect of DHEA administration in women, below the age of 42 starting their IVF treatment who met one of the two following features of POR (an abnormal ovarian reserve test and/or a previous poor response to ovarian stimulation in an IVF cycle) were assessed for eligibility [3]. The ESHRE working group on Poor Ovarian Response Definition of diminished ovarian reserves (AMH < 1.0 ng/mL or Day 2 or 3 FSH > 10 IU/L), or women with fewer than four oocytes retrieved with either standard long or antagonist protocols was used to defined POR in this study [3]. Inclusion criteria included women with diminished ovarian reserves (anti-müllerian hormone < 1.0 ng/mL or D2/3 follicle stimulating hormone > 10 IU/L), or women with fewer than four oocytes retrieved with either standard long or antagonist protocols. The study excluded women with previous or current DHEA supplementation, use of corticosteroids within the past three months, major systemic illnesses, and allergy to DHEA, women with women BMI > 37.5. A total of 60 subjects was enrolled into the study. 30 eligible patients received DHEA (Pharma Natural, USA) at the dose of 75 mg/day for three to eight months prior starting their controlled ovarian stimulation (COS), herein known as DHEA + , and 30 patients who received no treatment enrolled into the study (DHEA-; controls), but two DHEA + and six control patients did not complete the study. In DHEA + , two patients were withdrawn for stopping their DHEA treatment 4 and 5 months before their IVF treatment. In DHEA-, six subjects were withdrawn as they have postponed or stop further infertility treatment. Therefore, the study was conducted for 28 patients in the DHEA + group and 24 in the control DHEA- group. Table 1 summarizes the baseline characteristics of the patients in this study. Average age (DHEA-: 36 years; DHEA + , 37 years) and body mass index (DHEA-: 24 kg/m2; DHEA + , 23 kg/m2) were similar in both groups (p > 0.05). Baseline hormones between the two groups were comparable.

Table 1 Baseline characteristics of poor ovarian responder patients in this study

IVF/ICSI protocol

All individuals received the same stimulation protocol, same starting dose of gonadotropin, and fertilization technique. Briefly, the IVF/ICSI treatment cycle involved an antagonist-based COS protocol consisting of daily sub-cutaneous injections of recombinant-FSH (Puregon, Follitropin beta, 300iu; MSD, USA) and highly-purified human menopausal gonadotropin (Menopur; Menotropin, 150 IU; Ferring Pharmaceuticals, Germany) with initiation of gonadotropin releasing hormone antagonist (Ganirelix, Orgalutan, 0.25 mg s/c; MSD, USA) on day 5 of COS. The dose of Menopur and Puregon could be further increased depending on individual ovarian response. All patients had this standardized antagonist (short) protocol: no agonist (long) protocol was used. Human chorionic gonadotropin (i.m 10,000 IU hCG; Pregnyl; MSD, USA) was administered when at least one follicle measured ≥ 17 mm in diameter (averaged orthogonal measurements). The endometrial thickness, peak estradiol and progesterone levels were assessed on the day of human chorionic gonadotropin (hCG) trigger. Ultrasound-guided trans-vaginal oocyte retrieval was performed 36 h after hCG administration. The effect of DHEA supplementation on the markers of ovarian reserve (anti-müllerian hormone; AMH), follicular function (IGF-1), ovarian follicular levels of estradiol, testosterone, and DHEA, collected from the lead follicle at the time of OPU were assessed through ELISA as previously described [39].

Embryo transfer was performed on day 2 or day 3 of embryo-culture, and luteal phase support was achieved with vaginal progesterone (micronized progesterone, Utrogestan, 200 mg three times a day, Besins-International, France). Pregnancy was established by serum beta-hCG seventeen days post embryo transfer. Clinical pregnancy will be established by a transvaginal ultrasound four weeks after embryo transfer. IVF/ICSI clinical and hormonal outcomes are shown in Table 2. There was no significant difference in clinical pregnancy, number of oocytes retrieved, metaphase II oocytes or number of embryos transferred.

Table 2 Primary and secondary outcomes between DHEA+ and DHEA - control groups

Sample preparation

FF (DHEA + , N = 18 and DHEA-, N = 16) were divided for metabolomics and cytokine analyses. For untargeted metabolomics analysis, sample preparation followed previously published reports with some modifications [40, 41]. A volume of 50 µL from each FF sample was thawed at 4 °C, and FF proteins were precipitated with 200 µL ice-cold methanol. After vortexing, the mixture was centrifuged at 16,000 rpm for 10 min at 4 °C and the supernatant was collected and evaporated to dryness in a speedvacuum evaporator. The dry extracts were then re-dissolved in 200 µL of water/methanol (98:2; v/v) for liquid chromatography-tandem mass spectrometry (LC–MS/MS) analysis.

A pooled quality control (QC) sample was generated to allow comparison of analytic behavior over long periods of time. The pooled reference samples were for the purposes of quality control (i.e., to ensure relative consistency among identical samples within days) and for quality assurance (i.e., to ensure consistent results between days). They did not contribute data to downstream statistical analysis.

Liquid Chromatography-Tandem Mass Spectrometry-based Metabolomics

The supernatant fraction from sample preparation step was analyzed using Agilent 1290 ultra-high pressure (performance) liquid chromatography system (Waldbronn, Germany) equipped with a Bruker impact II Q-TOF mass spectrometer with its normal electrospray ionization (ESI) ion source (Bruker Daltonics). 2.5 μL of samples was injected and were separated using Waters Acquity HSS T3 (2.1 mm i.d. × 100 mm, 1.8 µm) at a flow rate of 0.2 mL/min. The oven temperature was set at 50 °C. The gradient elution involved a mobile phase consisting of (A) 0.1% formic acid in water and (B) 0.1% formic acid in methanol. The initial condition was set at 5% B. A 5.5 min linear gradient to 60% B was applied, followed by a 13 min gradient to 98% B (total 24 min including wash and re-equilibration) at a flow rate of 0.4 ml/min. The ion spray voltage was set at 4,500 V, and the Dry Temperature was maintained at 150 °C. The drying nitrogen gas flow rate and the nebulizer gas pressure were set at 8.0 L/min and 26 psi, respectively. Calibration of the system was performed using sodium formate clusters before data acquisition. The stability of the LC–MS method was examined and evaluated by sodium formate clusters (1 mM NaOH, 0.1% formic acid, 50% 2-propanol) infused into the system.

The ESI mass spectra were acquired in positive ion mode. Mass data were collected between m/z 100 and 1000 at a rate of three scans per second. Auto MS/MS was triggered at 8 Hz with duty cycle of 1.5 s. Threshold was set at 1500 counts, with active exclusion activated after 3 spectra, released after 0.3 min and overwritten if the current or previous intensity changes. MS/MS spectra were acquired at collision energy of 20–50 eV automatically varied by the charge states and the intensities of the selected precursors. Fragment spectra acquisition was carried out at a scan rate dependent on the MS precursor intensities—MS/MS spectra for high-intensity precursors were acquired for a shorter time (90,000 counts, 12 Hz) than low-intensity precursor ions (10,000 counts, 6 Hz) thus allowing for a balancing of maximal scan time and MS/MS spectral quality. As shown in Figure S1, the eight pooled quality control samples clustered in Principal Component Analysis (PCA) scores plots, and together with retention time CV% < 0.1 min, peak m/z values 3 mDa, and relative standard deviations of peak areas < 20%, there was good system stability, mass accuracy and reproducibility of the chromatographic separation during the whole LC–MS/MS sequence. In addition, intensity CV% of the identified compounds in pooled quality control samples are low (average 6%). PCA hotelling (T2) revealed one DHEA + subject as an outlier (D4) and was removed from further analysis (Figure S1).

Compound identification

Structure identification was achieved via the following in MetaboScape (version 2.0): elemental composition was predicted via isotopic pattern following the rules (i) mSigma of MS1: 20 with tolerance of 5 ppm and (ii) MS2: 50 with tolerance of 2 mDa of the differential metabolites was searched against Bruker HMDB (Human Metabolome Database) using a precursor match of ± 10 mDa, minimum score of 400 and minimum match score of 250. Progesterone, glycerophosphocholine, linoleic acid and valine were structurally confirmed using chemical standards.

Multiplex immunoassay analysis

Forty five  cytokines were detected and measured using ProCartaplex (EBioscience, CA, USA) as previously reported [BDNF, EGF, Eotaxin (CCL11), FGF-2 (FGF basic), GM-CSF, CXCL1 (GROα), HGF, IFNγ, IFNα, IL-1RA, IL-1β, IL-1α, IL-2, IL-4, IL-5, IL-6, IL-7, CXCL8 (IL-8), IL-9, IL-10, IL-12 p70, IL-13, IL-15, IL-17A, IL-18, IL-21, IL-22, IL-23, IL-27, IL-31, CXCL10 (IP-10), LIF, CCL2 (MCP-1), CCL3 (MIP-1α), CCL4 (MIP-1β), βNGF, PDGF-BB, PLGF, CCL5 (RANTES), SCF, CXCL12 (SDF1α), TNFα, LTA (TNFβ), VEGF-A, VEGF-D] [42]. Briefly, 5 μL of FFs were diluted with 5 μL Universal Dilution Buffer, and mixed with 50 μL of antibody-conjugated, magnetic beads in a 96 well DropArray plate (Curiox Biosystems, Singapore) and rotated at 450 rpm for 120 min at 25 °C while protected from light. Beads were internally dyed with different concentrations of two spectrally distinct fluorophores and covalently conjugated to antibodies against the 45 cytokines, chemokines and growth factors. The plate was washed three times with wash buffer (PBS, 0.05% Tween-20) on the LT210 Washing Station (Curiox) before adding 10 μL of secondary antibody and rotating at 450 rpm for 30 min at 25 °C protected from light. Subsequently, the plate was washed three times with wash buffer, and 10 μL of streptavidin–phycoerythrin added and rotated at 450 rpm for 30 min at 25 °C protected from light. The plate was again washed thrice with wash buffer; 60 μL of reading buffer was then added and the samples read using the Bio-Plex Luminex 200 (BioRad). The beads were classified by the red classification laser (635 nm) into its distinct sets, while a green reporter laser (532 nm) excites the phycoerythrin, a fluorescent reporter tag bound to the detection antibody. Quantitation of the 45 cytokines was then determined by extrapolation to a six or seven-point standard curve using five-parameter logistic regression modelling. Calibrations and validations were performed prior to runs and on a monthly basis respectively.

Statistical analysis

GraphPad Prism 6 (GraphPad Software Inc.) was used for performing all statistical analyses. Data were checked for normal distribution using Kolmogorov–Smirnov test. Unpaired or paired t-test was performed, as appropriate, to determine statistical significance between groups form normally distributed data. Mann–Whitney U test was used for non-normally distributed data. For comparing more than three groups, the data were analyzed using ANOVA test, followed by the t-test with Bonferroni adjustment. P < 0.05 was considered significant. Metabolomic data was further analyzed by Principal Component Analysis (PCA) and Partial Least Squares Regression (PLSR) modelling (Unscrambler X version 10.1) after the normalization of data by first centering the data to the median and scaling it by division with the standard deviation. Full cross-validation was applied in PLSR to increase model performance and for the calculation of β-coefficient regression values [43]. Metabolites with β-coefficient regression values ≥ 1 are considered to have contribute significantly to the PLSR model. In this study, metabolites fulfilling both PLSR β-coefficient regression values > 1.2 and Student t-test p < 0.05 was considered as differential.


Follicular fluid metabolomic analysis of POR patients

Significant increases were observed for clinical hormonal markers DHEA-S and free testosterone in the DHEA + patients compared to the DHEA- controls with treatment (p < 0.001; Figure S2). Estradiol, free testosterone and DHEA-S were significantly higher 4 months after treatment in the DHEA + patients (p < 0.001; Table 3).

Table 3 Comparison of AFC, ovarian volume and hormonal test at baseline and after 4 months of treatment in the DHEA group

From a total 2717 time-aligned features, an average of 903 features was chosen for auto MS/MS mode. From these, a total of 100 metabolites were identified via chemical standard confirmed HMDB [44]. An average of 65 MS/MS confirmed metabolites was identified per patient, which was similar in terms of metabolite identified in either DHEA- or DHEA + subjects (p = 0.8; range:59–76 metabolites), providing a global metabolome view of the FF metabolome in POR patients. The FF metabolome spanned three orders of magnitude, and was composed of a range of chemically diverse metabolites including, glycerophospholipids and derivatives (glycerophosphocholine, phosphatidylcholines), fatty acids (heptadecanoic acid, linoleic acid, vaccenic acid, myristic acid), cholesterols (isocaproic acid, 7-ketocholesterol), glucocorticoids (11-deoxycortisol or cortexolone, cortisol, corticosterone), hormones (17-hydroxyprogesterone, deoxycorticosterone, 11α-hydroxyprogesterone, 16-dehydroprogesterone, androstenedione, epitestosterone, progesterone, pregnenolone). Other metabolites included bile acids (3b-hydroxy-5-cholenoic acid, 3-oxocholic acid, glycocholic acid), peptides and derivatives (3-indolepropionic acid), lactones (delta-hexanolactone/caprolactone), lactic acid, vitamin D3 and sphingosine (Table 4).

Table 4 List of identified follicular fluid metabolites

Altered POR follicular fluid metabolome in response to DHEA

Next, partial least squares-discriminant regression analysis (PLSR) was performed based on the overall features to compare the FF metabolomic profiles between the DHEA + and DHEA- groups. At metabolome-scale, the PLSR score plot showed that the FF metabolome clearly distinguished DHEA + patients from DHEA- patients (Fig. 1). Further analysis revealed that in DHEA- controls, progesterone was the most abundant FF metabolite, followed by L-alanine, L-phenylalanine, pyridine, L-leucine. The top five metabolites in terms of abundance collectively made up close to half (48%) of the DHEA- FF metabolome (Fig. 2A). In DHEA + , the FF metabolome profile of highly abundant metabolites was different, with cortisol as the most abundance metabolite, followed by L-alanine, L-phenylalanine, pyridine, L-isoleucine and L-leucine. These top six metabolites collectively made up ~ 49.5% of the DHEA + FF metabolome (Fig. 2B). Interestingly, pyridine, considered a non-endogenous metabolite (HMDB0000926) was found in such high abundance suggesting it came from the synthesis of DHEA [45]. The observed MS/MS spectra of pyridine at various eV matched very well with HMDB database (Figure S3), which suggested its correct identification. As a precursor to testosterone and estrogen, DHEA could be converted to testosterone, and aromatized to estrogen; in the case of POR, exogenous DHEA was proposed to increase androgens in promoting folliculogenesis and potentiate the effects of gonadotropins [8, 46, 47]. FF testosterone was detected in our metabolomics profiling, although the differences between DHEA + and DHEA- group were small [DHEA-: mean signal intensity = 2294.5; DHEA + : mean signal intensity = 2267.75 (testosterone), p > 0.05; Figure S3B]. Next, we screened for prominent metabolites that are differential with DHEA treatment in POR using the criteria of β-coefficient regression values > 1.2 and p < 0.05. Among the FF metabolites, glycerophosphocholine, linoleic acid, progesterone, and L-valine fulfilled the screening criteria and were significantly lower in DHEA + relative to DHEA- (Student’s t-test, p < 0.05–0.005; Fig. 3A-D). Although not significant, pregnenolone, a cholesterol metabolite and steroid that is upstream of DHEA metabolism, was detected only in DHEA + (6/18 subjects), and not DHEA- (0/16 subjects). Receiver operating characteristic (ROC) analyses of the four metabolites revealed area under the curve (AUC) ranging from 0.711 (progesterone), 0.730 (glycerophosphocholine), 0.785 (linoleic acid) and 0.818 (L-valine) (p < 0.05–0.01; Fig. 3E-H), suggesting the plausible utility of these FF metabolites in monitoring DHEA treatment. Additionally, linoleic acid and L-valine remained significantly lower in DHEA + (p < 0.05, p < 0.001 for both) when women with endometriosis (N = 5) were removed from analysis, strongly suggesting the significant effect of DHEA on these metabolites (Figure S4).

Fig. 1
figure 1

Partial Least Squares Scores plot of DHEA- and DHEA+ follicular fluid metabolome. Metabolomic data was median centred and scaled by division with the standard deviation. The follicular fluid metabolome distinguished POR subjects on DHEA supplementation (DHEA+, red) and without DHEA supplementation (DHEA-, blue)

Fig. 2
figure 2

Histogram of follicular fluid metabolites in poor ovarian responders with and without DHEA supplementation. Follicular metabolome coverage and metabolite abundance as quantified by untargeted LC–MS/MS metabolomics in (A) DHEA- controls and (B) DHEA + poor ovarian response subjects. Metabolites were ranked according to their intensity counts

Fig. 3
figure 3

Significantly changed follicular fluid metabolites in DHEA+ and DHEA- patients. a-d Dot plots of significantly changed metabolites in poor ovarian responder patients. e-f Corresponding receiver operating curve (ROC) analyses of the metabolites. Area under curve (AUC) of the metabolites and their P-values are reported

Correlation of FF metabolites to biochemical hormones

In DHEA + patients, progesterone positively correlated with IGF-1 (Pearson r: 0.6757, p < 0.01); glycerophosphocholine negatively correlated with AMH (Pearson r: -0.5815; p < 0.05); linoleic acid correlated with estradiol and IGF-1 (Pearson r: 0.7016 and 0.8203, respectively; p < 0.01 for both; Figure S5A-D). Valine did not correlate with any biochemical hormones.

In DHEA- patients, few metabolites correlated weakly with estradiol, AMH,DHEA sulphate, serum-free testosterone, and IGF-1, with the exception of valine with serum-free testosterone (Pearson r: -0.8774; p < 0.0001; Figure S5E).

FF cytokine profile in response to DHEA

Of the 45 cytokines, chemokines and growth factors investigated, 22 were detected in human FF, comprising of 10 cytokines (IFNγ, IL12p70, IL13, IL1b, TNFα, IL1Ra, IL5, IL7, IL10, IL18), 6 chemokines (eotaxin, IP-10, MCP1, MIP1β, SCF, SDF-1α) and 8 growth factors (bNGF, BDNF, EGF, HGF, LIF, PIGF, VEGF-A, VEGF-D). Among them, FF MCP1, IFNγ, LIF and VEGF-D were significant lower in DHEA + compared to DHEA- (p = 0.03, 0.014, 0.031, 0.0161 respectively; Fig. 4). No correlation was found between the significant cytokines and metabolites.

Fig. 4
figure 4

Significantly changed follicular fluid cytokines in DHEA+ and DHEA- patients. Among 45 cytokines, chemokines and growth factors measured by multiplex immunoassay, A MCP-1, B IFNg, C LIF and (D) VEGF-D were significantly lower in POR subjects with DHEA supplementation. Student’s t-tests were performed and p<0.05 is considered statistically significant.


In this study, we report the first, and largest metabolome and large-scale cytokine study to-date of the FF of POR/DOR patients with DHEA supplementation. A total of 118 metabolites that included lipids, fatty acids, glucocorticoids, hormones, bile acids, peptides, and 22 cytokines were identified in the FF of POR/DOR patients. Among these, four FF metabolites, namely glycerophosphocholine, linoleic acid, progesterone, and L-valine varied significantly between DHEA + and DHEA- control subjects. This LC–MS/MS metabolomics study extends the human FF metabolome space in terms of characterization of its constituents, providing new insights into the complexities of oocyte development especially in POR women [48], as well as with DHEA supplementation [17,18,19, 49, 50].

The identification of four differential metabolites, glycerophosphocholine, linoleic acid, progesterone, and valine in this study may alert us to the metabolic effects of exogenous DHEA supplementation and plausibly using them to achieve improved outcomes [36, 51, 52]. Choline and derivatives are an emerging class of metabolites critical in developmental competence of fertilized oocytes [19], and glycerophosphorylcholine was found increased in the DHEA + group. Glycerophosphorylcholine is formed from the breakdown of phosphatidylcholine, and is an organic osmolyte, plausibly affecting concentrations of other constitutes of FF [16], and regulation of the diffusion of compounds into FF necessary for folliculogenesis and oogenesis [53]. PORs are known to exhibit a low diffusion of exogenous gonadotropin into FF, which is correlated with poor IVF outcomes [54]. It is conceivable that DHEA induced the metabolism of phosphatidylcholine to glycerophosphorylcholine. Progesterone is one of the key hormones for the progress of the first meiotic division in oocyte maturation, but changes to progesterone levels with DHEA supplementation has been controversial [55]. Our metabolomics study revealed for the first time that DHEA supplementation led to a decrease in FF progesterone levels, but to what impact lower progesterone induced by DHEA supplementation has on PORs remains to be investigated. Valine degradation has been previously reported in a proteomics study comparing competent versus incompetent buffalo oocyte proteome [56]. In a metabolomics study, valine metabolism was also identified in bovine cumulus and cumulus-oocyte-complex-conditioned media that undergo oocyte maturation [57]; although in both studies, valine was not directly detected in the omics profiling. In humans, degenerate oocytes or germinal vesicles that failed meiotically to reach metaphase II deplete valine more than competent oocytes. In other words, lower valine levels in culture media which is consistent with our results, and suggest plausible biological roles of valine in oocyte maturation. Interestingly, we noted a segregation of DHEA + patients with low and high level of valine, with the high valine group approaching concentrations of the DHEA- control group. Together with valine’s high AUC value in DHEA + , and that valine negatively correlated with testosterone, it is tempting to speculate that valine can be used as biomarker for monitoring individual DHEA supplementation. Linoleic acid is the most abundant polyunsaturated fatty acid in bovine [58] and human FF (Fig. 2), and varying concentrations of linoleic acid have reportedly different effects on oocyte maturation. At a concentration of 100 µM, linoleic acid added to maturation media inhibits bovine oocyte maturation and subsequent blastocyst development through increasing prostaglandin E2 concentration in the medium, decreasing intracellular cAMP, decreasing phosphorylation of the MAPK1 and AKT and inhibited germinal vesicle breakdown [58, 59]. Conversely, at concentrations at 50 µM or below, linoleic acid improved oocyte quality by increasing the content of neutral lipids stored in lipid droplets [59]. FF linoleic acid’s high AUC value and its correlation with estradiol suggest that can be another biomarker for titrating and monitoring individual DHEA supplementation.

The elevated DHEA-sulphate levels coupled with a lack of difference in FF testosterone with DHEA supplementation suggest the following possibilities in POR patients: (i) inadequate DHEA conversion to testosterone due to polymorphism in SULT2A1, CYP19A1 and FMR1 genes [60], or (ii) long CAG repeats in androgen receptor gene which is linked to its lower transcriptional activity at the promoters of genes involved in the metabolism of DHEA to testosterone [61]. The former is unlikely: in a case–control study involving 94 subjects, androgen secretion was not impaired in pre-ovulatory follicles of POR compared to normal responders, and similar levels of follicular testosterone levels was reported [62]. However, ethnicity and genetic predispositions might play a role as Chinese women are reported to have higher free androgens and African American women lower [63], which might explain their differences in pregnancy rates in association with IVF than those observed among other ethnic groups. Conversely, long CAG repeats is associated with risk of POR and oocyte insensitivity to androgenic stimulation [64], thus hinting a tenable rationale on the observed similar FF androgen levels between the DHEA- controls and DHEA + subjects in this study and others [36]. The abundance of cortisol in DHEA + subjects is interesting, in particular that DHEA reduces circulating cortisol [65], indicating follicular versus systemic difference in how DHEA affects cortisol levels. In vitro, it was noted that DHEA suppresses cortisol activity [66], including the antagonist effects of DHEA on the anti-inflammatory responses induced by cortisol via glucocorticoid receptor-mediated pathways [67]. It is noteworthy that high FF cortisol levels found in fertilized IVF individuals compared to unfertilized individuals led to the postulation that oocyte exposure to cortisol is required with oocyte maturation [68]. The higher levels of FF cortisol observed in DHEA + subjects therefore argues for a compensatory response to modulate the ratio of the two hormones in an attempt to maintain a favourable FF response to mature oocytes [67].

In mouse models of polycystic ovary syndrome, treatment with DHEA resulted in increased production of cytokines such as serum TNFα, IL-6, IL12p70, and IFNγ [69, 70]. In this study, DHEA supplementation led to the reduction of FF IFNγ, LIF, MCP-1, and VEGF-D levels. It appears that DHEA modulates chemokines and growth factors in POR FF without a clear Th1 or Th2 immune response as proposed [52]. LIF or leukemia inhibitory factor is expressed in the ovary and controls follicular growth [71]. It was reported that LIF suppressed the growth of primary, secondary, and early antral follicles in cultured ovarian tissues [72]. The authors postulated that LIF produced in the late antral or graafian follicles is secreted to suppress the growth of the neighbouring primary, secondary, and early antral follicles as part of follicular growth [16]. Interestingly, when hCG is administered in rhesus macaques, at 12 h follicular LIF levels increase and induce follicle rupture and ovulation and decrease at 24 h [73]. In our study, the number of MII oocytes and oocytes trended higher in the DHEA + group, suggesting that the biological roles of LIF might have been achieved (follicular maturation and rupture) but inadequate to generate a clinically significant outcome. In vitro results suggested that follicles produce VEGF-A, with VEGF-A inducing the expanding vasculature to support the increased needs growing follicles [74]. The decrease in VEGF-A in DHEA + individuals is intriguing. Fisher et al., described that in cultured follicles, the rise in VEGF-A levels in faster-growing follicles are dependent on FSH dose and oxygen tension [75]. There have been reports that DHEA inhibits oxygen consumption in neurons [76], tempting the postulation that DHEA inhibited oxygen consumption in follicle that subsequently led to lower production of VEGF-A in DHEA + individuals. Further, the lack of correlation between the significant cytokines and metabolites suggests that DHEA converting to steroids which subsequently modulate cytokine production within the follicular microenvironment is more complex than originally thought.

We note various strengths of this study. Firstly, due to the highly confident identification based on MS/MS, and mass accuracy of LC–MS/MS-based metabolomics, we were able to distinguish progesterone from DHEA, an advantage over interference-prone immunoassays that face a cross-reactivity bioanalytical problem [55]. Similarly, LC–MS/MS-based determination of androgens was preferred over immunoassays due to strong interference from DHEA [77]. We did not detect E1 and E2; because for phenolic hydroxyl group of estrogens to act as proton donors, the signal would be more sensitive in the negative ion mode electrospray ionization [78] than in the positive ion mode which was used in this study. Aside from previously reported constituents of FF such as linoleic acid [17, 22], amino acids [18], and steroids including progesterone, testosterone [79], this study also captured metabolic products of ovarian steroidogenesis, cholesterols and glucocorticoids in the FF. Secondly, this study is the largest-to-date, providing a global view, specifically of the effect of DHEA on the FF metabolome and cytokine profile in POR patients. Thirdly, the recruitment of women of relatively advanced age, low mean number of recovered oocytes (4.9) and low clinical pregnancy rate (5/52 or 9.6%), which contrast with other studies [22, 33] suggest the appropriate inclusion of POR patients.

Many studies of POR, including this study, suffer from important limitations. Firstly, POR patients represent a heterogeneous group of patients of different prognosis and with a range of patient and biochemical characteristics. To-date, no definition has been able to correlate their presentation with pregnancy or live-birth prognosis [80, 81], contributing to the difficulties in designing studies and trials to evaluate therapeutic modalities. Secondly, with respect to the similar clinical pregnancy rate between DHEA + and DHEA- subjects, this might be construed as a limitation in terms of sample size. However, various studies, including randomized clinical trials have not led to a clear outcome in terms of an improvement in clinical pregnancy rate, the number of oocytes retrieved, and embryos formed [31,32,33,34,35,36,37].

In conclusion, our study provided new insights to POR FF at the metabolome level, and as indicated from the FF metabolome analysis, exogenous DHEA to these patients altered the overall metabolome coverage and abundance to four metabolites. LC–MS global (untargeted) metabolomics analysis has provided the ability to reveal biologically relevant changes within a system, even at sensitive ranges before the precedence of gross morphological or phenotypical changes [82]. Hypotheses generated from this study included plausible mechanisms underlying DHEA metabolism, and the potential utility of glycerophosphocholine, linoleic acid, progesterone, and L-valine as markers to assess DHEA supplementation. Therefore, future directions include targeted quantitative LC–MS/MS approaches to be developed to detect and quantify four “responder” metabolites in approaches similar to those previously conducted on human peritoneal fluids and sera [83,84,85] to design treatment based on metabolomics profiles. Steroid hormones including testosterone should also be quantified via LC–MS/MS to establish baseline levels before commencing DHEA supplementation. Further, comparing POR and normal responders will provide further insights to the alteration of the FF metabolome, and reach a deeper understanding of underpinning pathophysiology to PORDisclosure of conflict-of-interest statement.

Availability of data and materials

Available upon request.





Electrospray ionization


Follicular fluid


Follicle stimulating hormone


Human chorionic gonadotropin


in vitro Fertilization


Liquid chromatography-tandem mass spectrometry


Poor ovarian responders


Partial least squares-discriminant regression


Quality control (QC); Principal Component Analysis


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

    Article  CAS  PubMed  Google Scholar 

  2. 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;1109.

  3. Ferraretti AP, La Marca A, Fauser BCJM, 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:1616–24.

    Article  CAS  PubMed  Google Scholar 

  4. Kolibianakis EM, Venetis CA, Diedrich K, Tarlatzis BC, Griesinger G. Addition of growth hormone to gonadotrophins in ovarian stimulation of poor responders treated by in-vitro fertilization: A systematic review and meta-analysis. Hum. Reprod. Update. 2009. p. 613–22.

  5. Haning RV, Flood CA, Hackett RJ, Loughlin JS, McClure N, Longcope C. Metabolic clearance rate of dehydroepiandrosterone sulfate, its metabolism to testosterone, and its intrafollicular metabolism to dehydroepiandrosterone, androstenedione, testosterone, and dihydrotestosterone in vivo. J Clin Endocrinol Metab. 1991;72:1088–95.

    Article  CAS  PubMed  Google Scholar 

  6. Casson PR, Santoro N, Elkind-Hirsch K, Carson SA, Hornsby PJ, Abraham G, et al. Postmenopausal dehydroepiandrosterone administration increases free insulin-like growth factor-I and decreases high-density lipoprotein: A six- month trial. Fertil Steril. 1998;70:107–10.

    Article  CAS  PubMed  Google Scholar 

  7. Barad D, Gleicher N. Effect of dehydroepiandrosterone on oocyte and embryo yields, embryo grade and cell number in IVF. Hum Reprod. 2006;21:2845–9.

    Article  CAS  PubMed  Google Scholar 

  8. Casson PR. Dehydroepiandrosterone supplementation augments ovarian stimulation in poor responders: a case series. Hum Reprod. 2000;15:2129–32.

    Article  CAS  PubMed  Google Scholar 

  9. Burger HG. Androgen production in women. Fertil Steril. 2002;77:S3-5.

    Article  PubMed  Google Scholar 

  10. Roy S, Mahesh VB, Greenblatt RB. Effect of dehydroepiandrosterone and δ4-androstenedione on the reproductive organs of female rats: Production of cystic changes in the ovary. Nature. 1962;196:42–3.

    Article  CAS  PubMed  Google Scholar 

  11. Broekmans FJ, Soules MR, Fauser BC. Ovarian Aging: Mechanisms and Clinical Consequences. Endocr Rev. 2009;30:465–93.

    Article  CAS  PubMed  Google Scholar 

  12. Barad D, Brill H, Gleicher N. Update on the use of dehydroepiandrosterone supplementation among women with diminished ovarian function. J Assist Reprod Genet. 2007. p. 629–34.

  13. Kyrou D, Kolibianakis EM, Venetis CA, Papanikolaou EG, Bontis J, Tarlatzis BC. How to improve the probability of pregnancy in poor responders undergoing in vitro fertilization: a systematic review and meta-analysis. Fertil Steril. 2009;91:749–66.

    Article  CAS  PubMed  Google Scholar 

  14. Bächler M, Menshykau D, De Geyter C, Iber D. Species-specific differences in follicular antral sizes result from diffusion-based limitations on the thickness of the granulosa cell layer. MHR Basic Sci Reprod Med. 2014;20:208–21.

    Article  Google Scholar 

  15. Seli E, Babayev E, Collins SC, Nemeth G, Horvath TL. Minireview: Metabolism of Female Reproduction: Regulatory Mechanisms and Clinical Implications. Mol Endocrinol. 2014;28:790–804.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Rodgers RJ, Irving-Rodgers HF. Formation of the ovarian follicular antrum and follicular fluid. Biol Reprod. 2010;82:1021–9.

    Article  CAS  PubMed  Google Scholar 

  17. O’Gorman A, Wallace M, Cottell E, Gibney MJ, McAuliffe FM, Wingfield M, et al. Metabolic profiling of human follicular fluid identifies potential biomarkers of oocyte developmental competence. Reproduction. 2013;146:389–95.

    Article  PubMed  Google Scholar 

  18. Lan Xia, Zhao X, Sun Y, Hong Y, Yuping Gao SH. Metabolomic profiling of human follicular fluid from patients with repeated failure of in vitro fertilization using gas chromatography/mass spectrometry. Int J Clin Experimenal Pathol. 2014;7:7220–9.

  19. 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:1078-1084.e8.

    Article  CAS  PubMed  Google Scholar 

  20. Piñero-Sagredo E, Nunes S, de los Santos MJ, Celda B, Esteve V. NMR metabolic profile of human follicular fluid. NMR Biomed. 2010;23:485–95.

  21. Ban Y, Ran H, Chen Y, Ma L. Lipidomics analysis of human follicular fluid form normal-weight patients with polycystic ovary syndrome: a pilot study. J Ovarian Res. 2021;14:135.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. He F, Wang F, Yang Y, Yuan Z, Sun C, Zou H, et al. The effect of growth hormone on the metabolome of follicular fluid in patients with diminished ovarian reserve. Reprod Biol Endocrinol. 2023;21:21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zamah AM, Hassis ME, Albertolle ME, Williams KE. Proteomic analysis of human follicular fluid from fertile women. Clin Proteomics. 2015;12:5.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Angelucci S, Ciavardelli D, Di Giuseppe F, Eleuterio E, Sulpizio M, Tiboni GM, et al. Proteome analysis of human follicular fluid. Biochim Biophys Acta - Proteins Proteomics. 2006;1764:1775–85.

    Article  CAS  Google Scholar 

  25. Jarkovska K, Martinkova J, Liskova L, Halada P, Moos J, Rezabek K, et al. Proteome mining of human follicular fluid reveals a crucial role of complement cascade and key biological pathways in women undergoing in vitro fertilization. J Proteome Res. 2010;9:1289–301.

    Article  CAS  PubMed  Google Scholar 

  26. Hanrieder J, Nyakas A, Naessén T, Bergquist J. Proteomic analysis of human follicular fluid using an alternative bottom-up approach. J Proteome Res. 2008;7:443–9.

    Article  CAS  PubMed  Google Scholar 

  27. Estes SJ, Ye B, Qiu W, Cramer D, Hornstein MD, Missmer SA. A proteomic analysis of IVF follicular fluid in women ≤32 years old. Fertil Steril. 2009;92:1569–78.

    Article  CAS  PubMed  Google Scholar 

  28. Chen F, Spiessens C, D’Hooghe T, Peeraer K, Carpentier S. Follicular fluid biomarkers for human in vitro fertilization outcome: Proof of principle. Proteome Sci. 2016;14:17.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Li C-J, Lin L-T, Tsui K-H. Dehydroepiandrosterone Shifts Energy Metabolism to Increase Mitochondrial Biogenesis in Female Fertility with Advancing Age. Nutrients. 2021;13:2449.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yakin K, Urman B. DHEA as a miracle drug in the treatment of poor responders; hype or hope? Hum Reprod. 2011;26:1941–4.

    Article  CAS  PubMed  Google Scholar 

  31. PG A, Simi G, Ruggiero M, Pinelli S, OM DB, Papini F, et al. DHEA supplementation improves follicular microenviroment in poor responder patients. Gynecol. Endocrinol. 2012. p. 669–73.

  32. Yeung TWY, Chai J, Li RHW, Lee VCY, Ho PC, Ng EHY. A randomized, controlled, pilot trial on the effect of dehydroepiandrosterone on ovarian response markers, ovarian response, and in vitro fertilization outcomes in poor responders. Fertil Steril. Elsevier Inc.; 2014;102:108–115.e1.

  33. Kara M, Aydin T, Aran T, Turktekin N, Ozdemir B. Does dehydroepiandrosterone supplementation really affect IVF-ICSI outcome in women with poor ovarian reserve? Eur J Obstet Gynecol Reprod Biol. 2014;173:63–5.

    Article  CAS  PubMed  Google Scholar 

  34. Chern C-U, Tsui K-H, Vitale SG, Chen S-N, Wang P-H, Cianci A, et al. Dehydroepiandrosterone (DHEA) supplementation improves in vitro fertilization outcomes of poor ovarian responders, especially in women with low serum concentration of DHEA-S: a retrospective cohort study. Reprod Biol Endocrinol. 2018;16:90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Chen S-N, Tsui K-H, Wang P-H, Chern C-U, Wen Z-H, Lin L-T. Dehydroepiandrosterone Supplementation Improves the Outcomes of in vitro Fertilization Cycles in Older Patients With Diminished Ovarian Reserve. Front Endocrinol (Lausanne). 2019;10.

  36. Hu Q, Hong L, Nie M, Wang Q, Fang Y, Dai Y, et al. The effect of dehydroepiandrosterone supplementation on ovarian response is associated with androgen receptor in diminished ovarian reserve women. J Ovarian Res. 2017;10:32.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Wang Z, Yang A, Bao H, Wang A, Deng X, Xue D, et al. Effect of dehydroepiandrosterone administration before in vitro fertilization on the live birth rate in poor ovarian responders according to the Bologna criteria: A randomised controlled trial. BJOG An Int J Obstet Gynaecol. 2022;129:1030–8.

    Article  CAS  Google Scholar 

  38. Hazeldine J, Arlt W, Lord JM. Dehydroepiandrosterone as a regulator of immune cell function. J Steroid Biochem Mol Biol. 2010;120:127–36.

    Article  CAS  PubMed  Google Scholar 

  39. Lee YH, Goh WW Bin, Ng CK, Raida M, Wong L, Lin Q, et al. Integrative toxicoproteomics implicates impaired mitochondrial glutathione import as an off-target effect of troglitazone. J Proteome Res. 2013;12:2933–45.

  40. Harden SL, Zhou J, Gharanei S, Diniz-da-Costa M, Lucas ES, Cui L, et al. Exometabolomic Analysis of Decidualizing Human Endometrial Stromal and Perivascular Cells. Front Cell Dev Biol. 2021;9.

  41. Loy SL, Zhou J, Cui L, Tan TY, Ee TX, Chern BSM, et al. Discovery and validation of peritoneal endometriosis biomarkers in peritoneal fluid and serum. Reprod Biomed Online. 2021;

  42. Zhou J, Chern BSM, Barton-Smith P, Phoon JWL, Tan TY, Viardot-Foucault V, et al. Peritoneal Fluid Cytokines Reveal New Insights of Endometriosis Subphenotypes. Int J Mol Sci. 2020;21:3515.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ku CW, Ong LS, Goh JP, Allen J, Low LW, Zhou J, et al. Defects in protective cytokine profiles in spontaneous miscarriage in the first trimester. F&S Sci. 2023;4:36–46.

    Article  Google Scholar 

  44. Wishart DS, Tzur D, Knox C, Eisner R, Guo AC, Young N, et al. HMDB: the Human Metabolome Database. Nucleic Acids Res. 2007;35:D521–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Williams JR, Boehm JC. Studies on the synthesis of dehydroepiandrosterone (DHEA) phosphatide. Steroids. 1995;60:333–6.

    Article  CAS  PubMed  Google Scholar 

  46. Nielsen ME, Rasmussen IA, Kristensen SG, Christensen ST, Mollgard K, Wreford Andersen E, et al. In human granulosa cells from small antral follicles, androgen receptor mRNA and androgen levels in follicular fluid correlate with FSH receptor mRNA. Mol Hum Reprod. 2011;17:63–70.

    Article  CAS  PubMed  Google Scholar 

  47. Gleicher N, Weghofer A, Barad DH. The role of androgens in follicle maturation and ovulation induction: friend or foe of infertility treatment? Reprod Biol Endocrinol. 2011;9:116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Emori MM, Drapkin R. The hormonal composition of follicular fluid and its implications for ovarian cancer pathogenesis. Reprod Biol Endocrinol. 2014;12:60.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Bedaiwy M, Shahin AY, AbulHassan AM, Goldberg JM, Sharma RK, Agarwal A, et al. Differential expression of follicular fluid cytokines: relationship to subsequent pregnancy in IVF cycles. Reprod Biomed Online. 2007;15:321–5.

    Article  CAS  PubMed  Google Scholar 

  50. Baskind NE, Orsi NM, Sharma V. Follicular-phase ovarian follicular fluid and plasma cytokine profiling of natural cycle in vitro fertilization patients. Fertil Steril. 2014;102:410–8.

    Article  CAS  PubMed  Google Scholar 

  51. Wiser A, Gonen O, Ghetler Y, Shavit T, Berkovitz A, Shulman A. Addition of dehydroepiandrosterone (DHEA) for poor-responder patients before and during IVF treatment improves the pregnancy rate: A randomized prospective study. Hum Reprod. 2010;25:2496–500.

    Article  CAS  PubMed  Google Scholar 

  52. Zhang J, Qiu X, Gui Y, Xu Y, Li D, Wang L. Dehydroepiandrosterone improves the ovarian reserve of women with diminished ovarian reserve and is a potential regulator of the immune response in the ovaries. Biosci Trends. 2015;9:350–9.

    Article  CAS  PubMed  Google Scholar 

  53. Hsueh AJW, Kawamura K, Cheng Y, Fauser BCJM. Intraovarian control of early folliculogenesis. Endocr Rev. 2015;36:1–24.

    Article  CAS  PubMed  Google Scholar 

  54. Nagata Y, Honjou K, Sonoda M, Sumii Y, Inoue Y, Kawarabayashi T. Pharmacokinetics of exogenous gonadotropin and ovarian response in in vitro fertilization. Fertil Steril. 1999;72:235–9.

    Article  CAS  PubMed  Google Scholar 

  55. Forman EJ, Franasiak JM, Scott RT. Elevated progesterone levels in women on DHEA supplementation likely represent assay interference. J Assist Reprod Genet. 2015;32:661–661.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Chen L, Zhai L, Qu C, Zhang C, Li S, Wu F, et al. Comparative Proteomic Analysis of Buffalo Oocytes Matured in vitro Using iTRAQ Technique. Sci Rep. 2016;6:31795.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Uhde K, van Tol HTA, Stout TAE, Roelen BAJ. Metabolomic profiles of bovine cumulus cells and cumulus-oocyte-complex-conditioned medium during maturation in vitro. Sci Rep. 2018;8:9477.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Marei WF, Wathes DC, Fouladi-Nashta AA. Impact of linoleic acid on bovine oocyte maturation and embryo development. Reproduction. 2010;139:979–88.

    Article  CAS  PubMed  Google Scholar 

  59. Carro M, Buschiazzo J, Ríos GL, Oresti GM, Alberio RH. Linoleic acid stimulates neutral lipid accumulation in lipid droplets of maturing bovine oocytes. Theriogenology. 2013;79:687–94.

    Article  CAS  PubMed  Google Scholar 

  60. Shohat-Tal A, Sen A, Barad DH, Kushnir V, Gleicher N. Genetics of androgen metabolism in women with infertility and hypoandrogenism. Nat Rev Endocrinol. 2015;11:429–41.

    Article  CAS  PubMed  Google Scholar 

  61. Borgbo T, Macek M, Chrudimska J, Jeppesen JV, Hansen LL, Andersen CY. Size matters: Associations between the androgen receptor CAG repeat length and the intrafollicular hormone milieu. Mol Cell Endocrinol. 2016;419:12–7.

    Article  CAS  PubMed  Google Scholar 

  62. de los Santos MJ, García-Laez V, Beltrán D, Labarta E, Zuzuarregui JL, Alamá P, et al. The follicular hormonal profile in low-responder patients undergoing unstimulated cycles: Is it hypoandrogenic? Hum Reprod. 2013;28:224–9.

  63. Sutton-Tyrrell K. sex hormone-binding globulin and the free androgen index are related to cardiovascular risk factors in multiethnic premenopausal and perimenopausal women enrolled in the Study of Women Across the Nation (SWAN). Circulation. 2005;111:1242–9.

    Article  CAS  PubMed  Google Scholar 

  64. Lledo B, Llácer J, Ortiz JA, Martinez B, Morales R, Bernabeu R. A pharmacogenetic approach to improve low ovarian response: The role of CAG repeats length in the androgen receptor gene. Eur J Obstet Gynecol Reprod Biol. 2018;227:41–5.

    Article  CAS  PubMed  Google Scholar 

  65. Kroboth PD, Amico JA, Stone RA, Folan M, Frye RF, Kroboth FJ, et al. Influence of DHEA Administration on 24-Hour Cortisol Concentrations. J Clin Psychopharmacol. 2003;23:96–9.

    Article  CAS  PubMed  Google Scholar 

  66. Buoso E, Lanni C, Molteni E, Rousset F, Corsini E, Racchi M. Opposing effects of cortisol and dehydroepiandrosterone on the expression of the receptor for Activated C Kinase 1: Implications in immunosenescence. Exp Gerontol. 2011;46:877–83.

    Article  CAS  PubMed  Google Scholar 

  67. Kamin HS, Kertes DA. Cortisol and DHEA in development and psychopathology. Horm Behav. 2017;89:69–85.

    Article  CAS  PubMed  Google Scholar 

  68. Keay SD. Higher cortisol:cortisone ratios in the preovulatory follicle of completely unstimulated IVF cycles indicate oocytes with increased pregnancy potential. Hum Reprod. 2002;17:2410–4.

    Article  CAS  PubMed  Google Scholar 

  69. Sander V, Luchetti CG, Solano ME, Elia E, Di Girolamo G, Gonzalez C, et al. Role of the N, N’-dimethylbiguanide metformin in the treatment of female prepuberal BALB/c mice hyperandrogenized with dehydroepiandrosterone. Reproduction. 2006;131:591–602.

    Article  CAS  PubMed  Google Scholar 

  70. Solano ME, Sander VA, Ho H, Motta AB, Arck PC. Systemic inflammation, cellular influx and up-regulation of ovarian VCAM-1 expression in a mouse model of polycystic ovary syndrome (PCOS). J Reprod Immunol. 2011;92:33–44.

    Article  CAS  PubMed  Google Scholar 

  71. Abir R, Fisch B, Jin S, Barnnet M, Freimann S, Van den Hurk R, et al. Immunocytochemical detection and RT-PCR expression of leukaemia inhibitory factor and its receptor in human fetal and adult ovaries. Mol Hum Reprod. 2004;

  72. Komatsu K, Koya T, Wang J, Yamashita M, Kikkawa F, Iwase A. Analysis of the Effect of Leukemia Inhibitory Factor on Follicular Growth in Cultured Murine Ovarian Tissue. Biol Reprod. 2015;93.

  73. Murphy MJ, Halow NG, Royer PA, Hennebold JD. Leukemia Inhibitory Factor Is Necessary for Ovulation in Female Rhesus Macaques. Endocrinology. 2016;157:4378–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Xu J, Bernuci MP, Lawson MS, Yeoman RR, Fisher TE, Zelinski MB, et al. Survival, growth, and maturation of secondary follicles from prepubertal, young, and older adult rhesus monkeys during encapsulated three-dimensional culture: Effects of gonadotropins and insulin. Reproduction. 2010;

  75. Fisher TE, Molskness TA, Villeda A, Zelinski MB, Stouffer RL, Xu J. Vascular endothelial growth factor and angiopoietin production by primate follicles during culture is a function of growth rate, gonadotrophin exposure and oxygen milieu. Hum Reprod. 2013;28:3263–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Safiulina D, Peet N, Seppet E, Zharkovsky A, Kaasik A. Dehydroepiandrosterone Inhibits Complex I of the Mitochondrial Respiratory Chain and is Neurotoxic In Vitro and In Vivo at High Concentrations. Toxicol Sci. 2006;93:348–56.

    Article  CAS  PubMed  Google Scholar 

  77. Heald AH, Butterworth A, Kane JW, Borzomato J, Taylor NF, Layton T, et al. Investigation into possible causes of interference in serum testosterone measurement in women. Ann Clin Biochem. 2006;43:189–95.

    Article  CAS  PubMed  Google Scholar 

  78. Guo T, Gu J, Soldin OP, Singh RJ, Soldin SJ. Rapid measurement of estrogens and their metabolites in human serum by liquid chromatography-tandem mass spectrometry without derivatization. Clin Biochem. 2008;41:736–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Wen X, Li D, Tozer AJ, Docherty SM, Iles RK. Estradiol, progesterone, testosterone profiles in human follicular fluid and cultured granulosa cells from luteinized pre-ovulatory follicles. Reprod Biol Endocrinol. 2010;8:117.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Oudendijk JF, Yarde F, Eijkemans MJC, Broekmans FJM, Broer SL. The poor responder in IVF: is the prognosis always poor? A systematic review Hum Reprod Update. 2012;18:1–11.

    Article  CAS  PubMed  Google Scholar 

  81. Esteves SC, Yarali H, Vuong LN, Carvalho JF, Özbek İY, Polat M, et al. Cumulative delivery rate per aspiration IVF/ICSI cycle in POSEIDON patients: a real-world evidence study of 9073 patients. Hum Reprod. 2021;36:2157–69.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Cui L, Lu H, Lee YH. Challenges and emergent solutions for LC-MS/MS based untargeted metabolomics in diseases. Mass Spectrom Rev. 2018;37:772–92.

    Article  CAS  PubMed  Google Scholar 

  83. Lee YH, Yang JX, Allen JC, Tan CS, Chern BSM, Tan TY, et al. Elevated peritoneal fluid ceramides in human endometriosis-associated infertility and their effects on mouse oocyte maturation. Fertil Steril. 2018;110:767-777.e5.

    Article  CAS  PubMed  Google Scholar 

  84. Lee YH, Tan CW, Venkatratnam A, Tan CS, Cui L, Loh SF, et al. Dysregulated sphingolipid metabolism in endometriosis. J Clin Endocrinol Metab. 2014;99:E1913–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Lee YH, Cui L, Fang J, Chern BSM, Tan HH, Chan JKY. Limited value of pro-inflammatory oxylipins and cytokines as circulating biomarkers in endometriosis – a targeted ‘omics study. Sci Rep. 2016;6:26117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors would like to thank the patients for their participation and nurses for their frontline service.

Disclosure of conflict-of-interest statement

The authors have nothing to disclose.


This work was supported by SingHealth Foundation (SHF/CTG034/2010) to VVF, and the National Medical Research Council Centre Grant Programme (NMRC/CG/M003/2017) to YHL. JKYC received salary support from Singapore’s Ministry of Health’s National Medical Research Council (NMRC/CSA(SI)/008/2016).

Author information

Authors and Affiliations



VVF and JCKY recruited patients. JZ performed metabolomics and cytokine analyses. DB performed statistical analysis. YT provided expert help on metabolomics. LYH and JCKY supervised the work and resources. All authors agreed to submission of the manuscript and have agreed to the order of authors.

Corresponding authors

Correspondence to Jerry. K. Y. Chan or Yie Hou Lee.

Ethics declarations

Ethics approval and consent to participate

The local Institutional Review Board approved the study (CIRB/2011/404/D) and written informed consent was obtained from each participant.

Consent for publication


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: Supplementary Figure 1.

Principle component analysis reveals DHEA+4 (arrow) as a potential outlier and was removed from subsequent analysis. Supplementary Figure 2. (a) MS/MS spectra of pyridine at increasing eV. (b) Follicular fluid testerosterone levels as measured by metabolomics. DHEA+, POR subjects on DHEA supplementation and DHEA- control without DHEA supplementation. Supplementary Figure 3. (a) Dot Plots of Linoleic acid and L-Valine after removal of women with endometriosis (N=5), (b) ROC curves of Linoleic acid and L-Valine after removal of women with endometriosis (N=5). Supplementary Figure 4. Histograms of estradiol, anti-müllerian hormone (AMH), DHEA-sulphate and insulin Growth Factor-1 (IGFBP-1) concentrations as determined by immunoassay. NS, not significant. Supplementary Figure 5. Scatter plots of (a) progesterone with IGF-1 (Pearson r: 0.6757, p<0.01), (b) linoleic acid with estradiol (Pearson r: 0.7016, p<0.01), (c) linoleic acid with IGF-1 (Pearson r: 0.8203, p<0.01), (d) glycerophosphocholine negatively correlated with AMH (Pearson r: -0.5815; p<0.05), (e) valine with serum-free testosterone (Pearson r: -0.8774; p<0.0001). Linear regression lines are shown.

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

Viardot-Foucault, V., Zhou, J., Bi, D. et al. Dehydroepiandrosterone supplementation and the impact of follicular fluid metabolome and cytokinome profiles in poor ovarian responders. J Ovarian Res 16, 107 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: