Skip to main content
Download PDF

UHPLC-MS-MS analysis of oxylipins metabolomics components of follicular fluid in infertile individuals with diminished ovarian reserve

Reproductive Biology and Endocrinology volume 19, Article number: 143 (2021) Cite this article

Abstract

Background

Diminished ovarian reserve (DOR) refers to a decrease in the number and quality of oocytes in the ovary, which results in a lack of sex hormones and a decline of fertility in women. DOR can potentially progress to premature ovarian failure (POF), which has a negative impact on women's quality of life and is a major cause of female infertility. Oxidative stress is a major contributor to fertility decrease in DOR patients, affecting the follicular microenvironment, oocyte maturation, fertilization, and embryo development. Understanding intracellular signal transduction can be achieved by defining specific oxidized lipid components in follicular fluid (FF) of DOR infertile patients.

Methods

The oxylipins metabolic signatures in the FF of DOR patients and females with normal ovarian reserve (NOR) enrolled for the in vitro fertilization (IVF) cycle were analyzed using UHPLC-MS-MS technology. Principal component analysis (PCA) and orthogonal projections to latent structure discriminant analysis (OPLS-DA) were used to analyze the derived metabolomic profiles. Pathway enrichment analysis was carried out using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and MetaboAnalyst databases. Furthermore, the Spearman rank correlation coefficient was used to determine the correlation between age, FSH, AMH, AFC, oocytes retrieved, MII oocytes, fertilization, high-quality embryos, and the concentration of differential oxidized lipid metabolites in FF.

Results

Fifteen oxylipins metabolites were found to be lower in the FF of DOR patients than those in the NOR group, including ±20-HDoHE, ±5-iso PGF2α-VI, 12S-HHTrE, 15-deoxy-Δ12,14-PGJ2, 1a,1b-dihomo PGE2, 1a,1b-dihomo PGF2α, 20-COOH-AA, 20-HETE, 8S,15S-DiHETE, PGA2, PGD2, PGE1, PGF1α, PGF2α, and PGJ2. The pathway enrichment analysis revealed that the 15 differentially oxidized lipid metabolites were closely related to the arachidonic acid metabolic pathway. Correlation analysis revealed that the concentration of 8 different oxidized lipid metabolites in FF was negatively correlated to FSH and positively correlated with AFC. AMH, the number of oocytes retrieved, MII oocytes and fertilization, were all positively correlated with 9 different oxidized lipid metabolites, but only one metabolite was positively correlated with the number of high-quality embryos.

Conclusions

Metabolomic analysis of FF revealed that oxylipins metabolism disorders were closely related to ovarian reserve function. Among these oxylipins metabolites, arachidonic acid metabolism undergoes significant changes that may be related to oocyte development, resulting in decreased fertility in DOR patients.

Trial registration

ChiCTR, ChiCTR2000038182, Registered 12 September 2020-Retrospectively registered

Introduction

Diminished ovarian reserve (DOR) refers to a decrease in the number and quality of oocytes in the ovary, which leads to decreased female fertility and reproductive endocrine function disorders. The most common clinical manifestation of women under the age of 40 years is rare menstruation, amenorrhea, and infertility. If not treated promptly, it may progress to premature ovarian failure (POF) within 1~6 years [1]. The incidence of DOR in infertile women of different ages is estimated to range from 6% to 64% [2]. According to the US Center for Disease Control and Prevention’s national assisted reproductive technology data from 2014, DOR patients accounted for about 30% of pregnant women [3]. Due to the reduction in basic antral follicle count (AFC), the reactivity of DOR patients to ovulation-promoting drugs decreased during in vitro fertilization /intracytoplasmic sperm injection-embryo transfer (IVF/ICSI-ET), resulting in some phenomena such as poor ovarian responder (POR), high cycle cancellation rate, decreased number of oocytes, decreased number of high-quality embryos, and decreased cumulative pregnancy rate [4].

Oxidative stress (OS) is the phenomenon in which excessive reactive oxygen species (ROS) are produced after being stimulated in vivo and in vitro, exceeding the antioxidant capacity in the body, and there is a serious imbalance between oxidation and antioxidant system in the body, resulting in oxidative stress damage to cells and tissues [5]. ROS primarily consists of superoxide anion (O2-·), hydrogen peroxide (H2O2), hydroxyl radical (·OH), etc. Under physiological conditions in vivo, ROS acts as a key second messenger in intercellular signal transmission and modulates gene expression to maintain cell homeostasis. Excessive ROS accumulation in tissue cells, on the other hand, cause damage to proteins, nucleic acids, and lipids, resulting in the mutation or deletion of nuclear DNA and mitochondrial DNA, lipid peroxidation of cell membranes leading to changes in membrane fluidity, protein oxidative damage leading to the inactivation of important enzymes such as superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GSH-Px), and finally induce apoptosis and tissue structure damage [6]. The cell membrane, which is rich in polyunsaturated fatty acids, is a susceptible target of ROS attack. ROS reacts with phospholipids, enzymes, and membrane receptor-related macromolecules on the membrane surface to form malondialdehyde (MDA), 4-Hydroxynonenal, isomeric prostaglandins, and other active compounds, particularly lipid peroxides [7]. OS plays an important role in the physiological and pathological processes of the female reproductive system, which is one of the major reasons for the decline of female reproductive function. According to recent research, OS can cause ovarian endocrine dysfunction, oocyte quality decline, granulosa cell apoptosis, and then follicular atresia, which is an important reason for the decline in fertility of DOR patients, and has a significant impact on the follicular microenvironment, oocyte maturation, fertilization and embryo development [8,9,10].

Follicular fluid (FF) is the fluid in which oocytes exchange materials and metabolize energy with surrounding cells. Its metabolite composition can represent the level of the FF environment and, in turn, the level of oocyte metabolism. It is made up of secretions produced by peripheral granulosa cells and serum diffused by local capillaries, as well as plasma components that have crossed the blood follicular barrier, primarily hormones, growth factors, interleukins, anti-apoptosis factors, proteins, sugars, amino acids, active oxygen, and antioxidant enzymes [11, 12]. It has an impact on oocyte maturation, follicular wall rupture, fertilization, and the development of early embryos [13, 14]. Oxidative stress markers in FF are closely related to the growth, development, and maturation of oocytes. Infertility caused by DOR has long been an inescapable problem in reproductive medicine, and changes in associated metabolites in FF indicate the quality of oocytes. As a result, studying the oxylipins metabolomics in FF of infertile patients with DOR is useful.

Metabolomics, as an integral component of system biology, mimics the research ideas of genomics and proteomics, through quantitative analysis of metabolites in organisms, and determining the relative link between metabolites and physiological and pathological changes [15]. Metabolomics has recently been used to identify potential biomarkers in FF. The study of FF metabolomics is critical for evaluating and predicting the potential of in vitro fertilization and oocyte embryo development. Merhi et al. [16] used metabolomics to identify biomarkers of ovarian reserve function. The results showed that advanced glycation end-products (AGE) in FF were positively correlated with the level of anti-Mullerian hormone (AMH), and AGE could be used as an important biological index for predicting ovarian reserve in ART patients. Another metabolomics study found that the concentration of glucose in the FF of DOR patients decreased significantly, whereas the concentrations of lactate and progesterone increased significantly. At the same time, glucose uptake, lactic acid production, platelet-type phosphofructokinase gene expression in granulosa and cumulus cells increased significantly, while progesterone concentration decreased significantly, implying that decreased aerobic metabolism, increased anaerobic metabolism and high progesterone stimulation of DOR oocytes may be important reasons for their quality decline and early embryo dysplasia [17]. The metabolomics of oxylipins in FF of DOR patients, on the other hand, has not been documented.

In this study, UHPLC-MS-MS technology was used to detect oxylipins metabolites in FF of DOR infertile patients using target metabolomics, to identify the oxidized lipid metabolites in FF of DOR infertile patients, to screen potential biomarkers, and to investigate the possible mechanism of oocyte quality decline in DOR, which has far-reaching implications for understanding the pathogenesis of DOR.

Materials and methods

Sample collection and preparation

Patients with DOR (n= 20) and women with normal ovarian reserve (NOR) (n = 20, as a control group) were recruited from the Reproductive Medicine Center of Shanghai first maternity and infant hospital, Tongji university school of medicine, between May 2020 and January 2021. The study was approved by the Ethics Committees of Shuguang hospital which is affiliated with Shanghai University of traditional Chinese medicine (No.2020-833-40-01) and it was registered with the China clinical trial registration center (ChiCTR2000038182). Before the study, all subjects signed informed consent. DOR was diagnosed in patients based on the following diagnostic criteria [18]: (1) aged 20~40 years old. (2) requires one of the following three features to be met: 1) serum levels of basic follicle-stimulating hormone (FSH): 10~30 IU/L with or without FSH/ luteinizing hormone (LH) ≥ 3; 2) serum levels of AMH: 20~25 years old (AMH≤3.0 ng/ml), 26~30 years old (AMH≤2.5 ng/ml), 31~35 years old (AMH≤1.5 ng/ml), 36~40 years old (AMH≤1.0 ng/ml); 3) basic AFC on menstrual cycle (MC) 2~3: AFC≤ 5. DOR was diagnosed by meeting (1) and any two of (2). The control group included NOR patients who were infertile solely owing to tubal factors. The following were the detailed inclusion criteria: (1) aged 20~35 years old. (2) serum levels of basic FSH < 10IU/L, serum levels of basic LH<10 IU/L, with FSH/LH < 3. (3) the number of AFC: 8 < AFC < 24. (4) serum levels of AMH: 20~25 years old (3.0ng/ml<AMH≤8.4ng/ml), 26~30years old (2.5ng/ml< AMH≤8.4 ng/ml), 31~35 years old (1.5 ng/ml < AMH≤8.4ng/ml).

The study exclusion criteria were as follows: (1) endometriosis. (2) diagnosis of the polycystic ovarian syndrome. (3) presence of a functional ovarian cyst with E2 >100 pg/ml. (4) history of ovarian surgery such as ovarian cyst stripping and teratoma stripping. (5) receipt of hormone treatments within the preceding 3-month period. (6) any contraindications to ovarian stimulation treatment.

To avoid the influence of the ovulation induction scheme on the metabolites in FF in this study, all recruited patients, including DOR and NOR patients, adopted a micro-stimulation strategy for ovulation induction. Letrozole (Jiangsu Hengrui Medicine Co., China) was given orally for 5 days, beginning on MC 3, with a 5 mg daily dose. We began ovarian stimulation with a daily injection of gonadotropin (Gn) until the trigger day when the dominant follicular diameter of ≥ 10 mm was confirmed by transvaginal ultrasonography after MC 6. The protocol for human menopausal gonadotropin (hMG) (Zhuhai Lizhu Pharmaceutical Trading Co., China) was used as reported previously [19]. hMG was injected intramuscularly daily in alternate doses of 150 and 225 IU (150 IU of hMG given on the first day of ovarian stimulation, 225 IU on the second day, 150 IU on the third day, and so on). The hMG initiation dose was similar in the DOR and NOR groups. After 5~7 days, serum sex hormones [FSH, LH, estrogen (E2), and progesterone (P)] were detected, as well as a re-examination of the vaginal ultrasound. The hMG dosage was adjusted based on the level of sex hormones in the serum as well as the number and size of developing follicles. If more than one dominant follicle with a diameter≥18mm or more than three follicles with a diameter ≥16mm were present, 3000~5000IU of human chorionic gonadotropin (hCG) (Zhuhai Lizhu Pharmaceutical Trading Co., China) was used to induce oocyte maturation. We used transvaginal ultrasound-guided follicle aspiration to retrieve oocytes 36~38 h after hCG administration. Following oocyte separation, FF obtained from 3 dominant follicles was pooled in a 15 ml centrifuge tube and centrifuged at 3000 rpm for 15 min at 4°C to remove insoluble particles and cells. The supernatant was packaged in a 2ml freezing tube, and stored at -80°C in an ultra-low temperature storage box to avoid repeated freezing or thawing.

Serum hormone measurement and antral follicle calculation

Serum sex hormones (FSH, LH, E2, and P) were detected in all patients by radio-immunoassay and transvaginal ultrasonography for basic AFC on MC 2~3.

Metabolites extraction

A 1000 μL aliquot of each sample was accurately transferred to an Eppendorf tube. Following the addition of the isotopically-labeled internal standard mixture, the samples were vortexed for 30 s and sonicated for 5 min in an ice-water bath. SPE was used to further purify the sample. The SPE cartridges were equilibrated with 1 mL of methanol and 1 mL of water. After loading a sample (supernatant obtained as mentioned above), the cartridge was washed with 1 mL of 5% MeOH/H2O (v/v). The flow-through fraction was then discarded. Finally, the samples were eluted with 1 mL of MeOH, the eluent was evaporated to dryness under a gentle stream of nitrogen, and the eluent was reconstituted in 100 μL of 30% ACN/H2O (v/v). The clear supernatant was subjected to UHPLC-MS/MS analysis, after centrifugation (15 min, 12000 rpm, 4°C).

UHPLC-MRM-MS analysis

The UHPLC separation was performed using an EXIONLC System (Sciex), equipped with a Waters ACQUITY UPLC BEH C18 column (150 × 2.1 mm, 1.7 μm, Waters). The mobile phases A and B were 0.01% formic acid in the water, and 0.01% formic acid in acetonitrile, respectively. The temperature of the column was set to 50°C. The temperature of the auto-sampler was set to 4°C, and the injection volume was 10 μL.

For assay development, a SCIEX 6500 QTRAP+ triple quadrupole mass spectrometer (Sciex), equipped with an IonDrive Turbo V electrospray ionization (ESI) interface was used. Typical ion source parameters were as follows: Curtain Gas = 40 psi, IonSpray Voltage = -4500 V, temperature = 500°C, Ion Source Gas 1 = 30 psi, and Ion Source Gas 2 = 30 psi.

Flow injection analysis was used to analyze the MRM parameters for each of the targeted analytes by injecting the standard solutions of the individual analytes, into the API source of the mass spectrometer. In the MRM scan mode, several most sensitive transitions were used to optimize the collision energy for each Q1/Q3 pair. The Q1/Q3 pairs with the highest sensitivity and selectivity among the optimized MRM transitions per analyte were selected as a “quantifier” for quantitative monitoring. The additional transitions served as a “quantifier” to confirm the identity of the target analytes. SCIEX Analyst Work Station Software (Version 1.6.3) and Multiquant 3.03 software (Version 20.2) were used for MRM data acquisition and processing.

Data collection, processing and statistical analysis

The data were statistically analyzed using SPSS 25.0 software. The baseline characteristics of the study population are described. The continuous variables with a normal or near-normal distribution were expressed as the means ± standard deviations (SDs) and were analyzed by Student’s t-test, otherwise, a Mann–Whitney U test was used. The frequency (composition ratio) of the classification data was statistically described and evaluated using the chi-square test. For correlation analyses, the Spearman rank correlation coefficient was used. Statistical significance was defined as P < 0.05.

For the metabolomics analysis, the missing values in the raw data were filled by half of the minimum value. In addition, the overall normalization method was used for data analysis. SIMCA software (V16.0.2, Sartorius Stedim Data Analytics AB, Umea, Sweden) was used to perform multivariate analysis, principal component analysis (PCA), and orthogonal projections to latent structure discriminant analysis (OPLS-DA). The original data distribution was shown using PCA. OPLS-DA was used to examine the separation between two groups and to better understand the variables responsible for classification. The variable importance in the projection (VIP) of the first principal component obtained from the OPLS-DA analysis was determined. In the univariate analysis, metabolites with a VIP > 1.0 and P-value < 0.05 were considered significantly different. Moreover, the quality of the OPLS-DA model was assessed using standard parameters (R2Y and Q2). The Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.Genome.jp/kegg/) and the MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/) database were used for pathway enrichment analysis.

Results

Clinical characteristics and ovulation outcomes

The average age, body mass index (BMI), duration of infertility, type of infertility, bAFC, AMH, bFSH, bLH, bE2, P, the numbers of oocytes retrieved, MII oocytes, fertilization, and high-quality embryos of the participants are presented in Table 1. There was no significant difference between DOR and NOR groups in terms of BMI, duration of infertility, type of infertility, bLH, bE2, and P levels (P>0.05). The average age of the DOR group was significantly older than the NOR group (P <0.05), bFSH levels in patients with DOR were significantly higher than in NOR patients (P<0.05). However, bAFC, AMH, the number of retrieved oocytes, MII oocytes, fertilization, and high-quality embryos were significantly lower in DOR patients compared to NOR patients (P<0.05). This indicated that the ovarian function of the DOR group was significantly low.

Table 1 The demographic and clinical characteristics of patients with DOR and NOR
Full size table

Multivariate analysis of metabolites

After relative standard deviation de-noising, 103 peaks were detected and 68 metabolites were identified in this study. The median vale was used to fill in the missing variables. The final dataset, which included information on the peak number, sample name, and normalized peak area was imported into SIMCA software for multivariate analysis. To minimize the influence of both noise and high variance of the variables, the data were scaled and logarithmically transformed. Following these changes, PCA and OPLS-DA were used to comprehensively compare the FF metabolomic profiles and determine the degree of diversity between the DOR and NOR groups. Sixteen quality control (QC) samples were included to confirm the stability and repeatability of the system, and the QCs were clustered together and separated from the study subjects' samples, indicating that the PCA was correct (Fig. 1a). The PCA score plot revealed significant differences between DOR and NOR, with one outlier located beyond the 95% Hotelling’s T-squared ellipse (Fig. 1b). As shown in Fig. 1b, the majority of DOR samples clustered to the left, while all NOR samples clustered to the right. However, one DOR sample on the right was also visible in the PCA score plots.

Fig. 1
figure 1

PCA score plots, OPLS-DA score plots, and corresponding validation plot of OPLS-DA results derived from FF metabolomics profiles comparing DOR and NOR. a. PCA score plot with QC samples. b. PCA score plot without the QC samples. c. OPLS-DA score plot. d. The OPLS-DA model’s permutation test

Full size image

An OPLS-DA was performed to compare the metabolite between DOR and NOR. The OPLS-DA score plot revealed that the oxidative lipid metabolites in FF differed significantly between DOR and NOR samples (Fig.1c). As shown in Fig.1c, all NOR samples clustered to the right, while the majority of DOR clustered to the left. However, there was one DOR sample located on the right. This finding suggested that there were substantial changes in oxidative lipid metabolites between DOR and NOR patients. Furthermore, the permutation test yielded R2Y(cum) and Q2(cum) values of 0.36 and −0.8, respectively, indicating the lack of overfitting and the high predictive ability of the OPLS-DA model (Fig. 1d), and the suitability of the model for subsequent optimization analyses.

Identification of significantly different metabolites and pathways associated with DOR

In this study, 15 FF metabolites with VIP values > 1 in the OPLS-DA analysis and P-values < 0.05 in the univariate analysis were detected between DOR and NOR samples (Table 2), indicating that these 15 oxidized lipid metabolites were associated with ovarian reserve function. In particular, 15 metabolites were all reduced and no metabolites were increased in the DOR group compared to the NOR group. These 15 differentially oxidized lipid metabolites associated with ovarian reserve function were then summarized in a volcano plot (Fig. 2a): the 15 down-regulated metabolites (blue spots) were located on the left, while no up-regulated metabolites were found in the FF of DOR compared to the NOR group. Furthermore, a heatmap was used to distinguish the upregulated and downregulated metabolites in DOR patients compared with NOR patients (Fig. 2b). R software (Version 3.6.3) was also used to display the variations in 15 oxidized lipid metabolites in FF between the two groups. Among the differential metabolites, the levels of 15 differential metabolites were all lower in the DOR group than in the NOR control group (Fig.3).

Table 2 Differential metabolites in FF samples between DOR and NOR groups
Full size table
Fig. 2
figure 2

Identification of the differential metabolomics profiles of FF between DOR and NOR patients based on a volcano plot and hierarchical clustering analysis. a. Volcano plot, down-regulated and up-regulated metabolites in DOR compared to NOR are marked in blue and red, respectively. The X-axis represents the log2 fold change of metabolites, while the Y-axis represents the fold change of the –log10 P value determined by the Student’s t-test. The variable importance in the projection (VIP) value is represented by the dot size. b. Heatmap of the hierarchical clustering analysis. There are fifteen distinct metabolites presented

Full size image
Fig. 3
figure 3

a-o The violin diagram shows the concentration of 15 different metabolite profiles in the quantitative analysis of DOR and NOR groups. The y-axis represents integration data for mass spectral concentration. Each component group is represented by the x-axis. *P <0.05, **P <0.01, ***P<0.001

Full size image

KEGG pathway analysis was performed to investigate the enriched metabolic pathways associated with ovarian reserve function, and an overexpression analysis of the differentially oxidized lipid metabolites was performed to reveal the biological functional background of FF metabolic disorders in infertile patients with DOR. According to the results of the KEGG Mapper Search, differentially oxidized lipid metabolites were involved in metabolic pathways, including arachidonic acid (AA) metabolism, metabolic pathways, serotonergic synapse, neuroactive ligand-receptor interaction, asthma, phospholipase D signaling pathway, vascular smooth muscle contraction, oxytocin signaling pathway, FC epsilon RI signaling pathway, African trypanosomiasis, bile secretion, and ovarian steroidogenesis. In addition, we input 15 distinct metabolites into the MetaboAnalyst 5.0 database for metabolic pathway enrichment analysis. Only one major AA metabolic pathway was discovered using the Holm P-value, false discovery rate (FDR), and Impact value. As a result, when combined with the findings of the KEGG database analysis, it was hypothesized that the AA metabolic pathway was closely related to DOR (Fig.4).

Fig. 4
figure 4

KEGG’s metabolic pathway for arachidonic acid. The pink dots represent differentially expressed metabolites involved in the arachidonic acid metabolic pathway, all of which were reduced in DOR patients. KEGG, Kyoto Encyclopedia of Genes and Genomes; DOR, diminished ovarian reserve

Full size image

Predictive value of 15 oxidized lipid metabolites in follicular fluid on ovarian reserve function

Using ROC curve analysis, the predictive value of 15 oxidized lipid metabolites in ovarian reserve function was determined, and the AUC, 95%CI, specificity, and sensitivity were obtained. At the same time, the well-known ovarian reserve function markers such as age, FSH, AMH, and AFC, were assessed and used as a reference for lipid metabolite oxidation (Table 3). The levels of ±20-HDoHE (AUC=0.782), 12S-HHTrE (AUC=0.762), 20-COOH-AA (AUC=0.758), 8S,15S-DiHETE (AUC=0.800), PGA2 (AUC =0.850), and PGE1 (AUC=0.818) in FF showed high sensitivity and specificities for the ovarian reserve function, with predictive values comparable to age (AUC=0.805). Furthermore, we found that the AUC values of FSH, AMH, and AFC were all greater than 0.95, which was superior to the conventional study to assess their clinical value [20]. The reason for the analysis could be that the patients in this study adhered to FSH, AMH, and AFC standards, demonstrating that the patients included are credible.

Table 3 The ROC reports between DOR and NOR group
Full size table

Correlation analysis of oxidized lipid metabolite concentrations and related parameters

Spearman correlation analysis was used to investigate the relationships between different oxidized lipid metabolites in FF and age, FSH, AMH, AFC, retrieved oocytes, MII oocytes, fertilization, high-quality embryos. The concentrations of 20-HDoHE, ±5-iso PGF2α-VI, 15-deoxy-Δ12,14-PGJ2, 1a,1b-dihomo PGE2, 20-COOH-AA, PGA2, PGD2, and PGE1 were found to be negatively correlated with FSH. Significant positive correlations were found between 20-HDoHE, ±5-iso PGF2α-VI, 12S-HHTrE, 1a,1b-dihomo PGE2, 20-COOH-AA, 8S,15S-DiHETE, PGA2, PGE1, PGF2α, and AMH. The concentrations of 20-HDoHE, ±5-iso PGF2α-VI, 1a,1b-dihomo PGE2, 20-COOH-AA, 8S,15S-DiHETE, PGA2, PGE1, and PGF2α were positively correlated with AFC. The oocytes retrieval was positively correlated with 20-HDoHE, ±5-iso PGF2α-VI, 12S-HHTrE, 1a,1b-dihomo PGE2, 20-COOH-AA, 8S,15SDiHETE, PGA2, PGE1, and PGF2α. MII oocytes and fertilization were positively correlated with the concentrations of 20-HDoHE,12S-HHTrE,1a,1b-dihomo PGE2, 20-COOH-AA, 8S,15S-DiHETE, PGA2, PGD2, and PGF2α. Only 20-COOH-AA had a positive correlation with high-quality embryos (Table 4).

Table 4 Correlation analysis between the concentration of oxidized lipid metabolites and related parameters
Full size table

Discussion

Oxylipins, also known as lipid mediators, are a class of oxidative metabolites produced by the autooxidation of polyunsaturated fatty acids (AA, linoleic acid, alpha-linolenic acid, DHA, EPA, etc.) or by cyclooxygenase (COX), lipoxygenase (LOX), and cytochrome P450 (CYP) enzymes. Oxylipins are signal transduction molecules that participate in nearly all physiological functions of the body and play a critical regulatory role in the organism’s life activities, including inflammatory response, immune defense, endocrine regulation, oxidative stress, etc [21]. Besides, it is closely related to the occurrence and development of a variety of diseases, including tumor [22], cardiovascular disease [23], diabetes [24], lung disease [25], and Alzheimer's disease [26]. However, no studies have been conducted to investigate the relationship between oxylipins and the occurrence and development of DOR. Therefore, we performed a UHPLC-MS-MS analysis of the composition of FF oxylipins metabolites in patients with DOR and NOR. Our findings demonstrated changes in oxylipins metabolites in FF of infertile patients with DOR, as well as changes in ovarian reserve function in infertile patients with DOR. We identified 15 oxylipins metabolites that were associated with DOR, including ±20-HDoHE, ±5-iso PGF2α-VI, 12S-HHTrE, 15-deoxy-Δ12,14-PGJ2, 1a,1b-dihomo PGE2, 1a,1b-dihomo PGF2α, 20-COOH-AA, 20-HETE, 8S,15S-DiHETE, PGA2, PGD2, PGE1, PGF1α, PGF2α, and PGJ2. The levels of the 15 differentially oxidized lipid metabolites in the FF of DOR patients were all lower than in the NOR group, and the differences were statistically significant. Based on the KEGG and Metaboanalyst database, the results of pathway enrichment analysis showed that the differential oxylipins metabolites were closely related to the AA metabolic pathway. These distinct oxylipins metabolites and the AA metabolic pathway may play a non-negligible role in DOR infertility. Dysfunctions in oxylipins metabolism in FF of DOR patients related to follicular development may provide potential diagnostic and therapeutic targets for promoting oocyte maturation.

Polyunsaturated fatty acids (PUFAs) are long-chain fatty acids with two or more double bonds, the n-3 and n-6 series fatty acids having biological significance. Imbalanced dietary intake of n-6 to n-3 PUFAs has been associated with cardiovascular and cerebrovascular diseases, cancer, inflammation, and autoimmune diseases [27]. Wathes et al. [28] demonstrated that dietary n-6 and n-3PUFAs can alter reproductive processes by acting as precursors for prostaglandin synthesis and by regulating the expression patterns of many key enzymes involved in prostaglandin and steroid metabolism. AA is an n-6 PUFA essential fatty acid that is cleaved from phospholipids by phospholipase A2 (PLA2), phosphatidyl inositol specific phospholipase C, and the fatty acylglycerol lipase. AA is then metabolized to eicosanoids via COX, LOX, or P450-dependent cyclooxygenase pathways [29]. A previous study demonstrated that high concentrations of AA derivatives in human FF at the time of oocyte retrieval significantly decreased oocyte’s ability to form pronuclei after ICSI [30]. As cell signaling intermediates, AA and its derivatives regulate cAMP activation, Ca2+ influx, Ca2+/CaM, PKC, MAPK, and PI3K/Akt, resulting in regulation of cellar growth, proliferation, and differentiation [31, 32]. Another study reported that AA is one of the most abundant polyunsaturated fatty acids [33]. The concentration of AA in FF in bovine ovaries accounts for about 2.5% of total fatty acids, whereas the concentration in plasma is 1.2% [34]. It has been reported that follicular AA regulates cumulus granulosa cells in human [35] and non-human mammals [33]. Zhang et al. [36] demonstrated that AA regulates the survival, gene expression, lipid formation, steroid production, and intracellular signaling pathways in bovine granulosa cells. This supports the hypothesis that AA in FF can directly influence the function of granulosa cells, thereby modulating follicular development and ovulation. Khajeh et al. [37] investigated the biological effects of AA in human cumulus granulosa cells (CGCs) after exposure to acetylsalicylic acid (ASA). They found that ASA treatment reduced E2 production, Cyp19a1 expression, GSH-Px activity, and estradiol receptor expression in CGCs. Importantly, the addition of AA abolished the ASA-induced reduction in E2 levels and expression of Cyp19a1, increased the antioxidant capacity of CGCs exposed to ASA by increasing GSH-Px activity, synthesis, and secretion of PGE2, and increasing the expression of the estrogen receptors. Li et al. [38] found that the levels of AA and its metabolites were elevated in FF of patients with PCOS, and insulin increased the generation of AA metabolites by COX-2. This is likely to be a novel molecular pathophysiological mechanism of PCOS. In the present study, we found that the level of AA metabolites in FF of patients with DOR was lower compared to its level in the NOR control group, and these decreased AA metabolites were generated via COX-2 and CYP pathway, but not LOX pathway. CYP enzymes mediate AA transformation into 20-hydroxyeicosatetraenoic acid (20-HETE) [39]. A study revealed a higher serum concentration of 20-HETE in obese women as compared to that in women with normal weight [40]. Elsewhere, 20-HETE concentration in urine and plasma of patients with metabolic syndrome was found to be higher as compared to that of control [41]. There is more evidence of a significant correlation between the level of 20-HETE in blood with BMI, whereby higher BMI denote higher levels of 20-HETE [42]. 20-HETE significantly impacts obesity associated with a high-fat diet, impaired insulin signaling, and insulin resistance (IR) [43]. Li et al. [38] found higher levels of 20-HETE in the FF of non-obese PCOS patients that underwent IVF as compared to that in the non-PCOS group, but the difference was not significant [38]. We, herein, reported lower levels of 20-HETE in FF of DOR patients as compared to that of subjects in the NOR control group. However, the correlation and mechanism between 20-HETE and DOR remain elusive.

AA is a well-known precursor of prostaglandins (PGs), which are synthesized PGs by prostaglandin-endoperoxide synthase S1 (PTGS1) or PTGS-2 previously identified as cyclooxygenase enzymes COX1 and COX2 [44, 45]. PTGS-2(COX2) is the critical prostaglandin-endoperoxide synthase in the ovary. Mounting evidence indicates that PTGS-2(COX2) deficient female mice exhibit multiple failures in female reproductive processes, including ovulation, fertilization, implantation, and decidualization [46, 47]. Prostaglandins, especially PGE2 and PGF2α, have previously been shown to regulate oocyte maturation, ovulation, and cumulus expansion [48, 49]. Additional evidence indicates that PGE2 plays a crucial role in protecting oocytes against oxidative stress, therefore, holds promise as a novel autocrine/paracrine player in the mechanisms that potentially drive successful oocyte maturation and oocyte survival in cow [50]. Moreover, PGF2α is synthesized in response to LH surge, however, its underlying role in the ovulatory process remains to be determined [51]. Some researchers suggest that both PGF2α and PGE2 are important in the ovulatory process [52, 53]. Pereira de Moraes et al. [54], however, concluded that PGF2α alone is unable to induce ovulation in cattle. Sharma et al. [55] investigated PGE2 treated follicles and revealed a reduction in the attributes of apoptosis in granulosa cells, whereas PGF2α saw an increase in apoptotic characteristics. Furthermore, Kemiläinen et al. [56] found a broad expression of HSD17B12 enzyme in both human and mouse ovaries, demonstrating a role for HSD17B12 in the synthesis of PGD2, PGE2, PGF2α, and TXB2, which are well-known prostaglandins for their regulatory role in ovarian function. Another product of the AA cyclooxygenase pathway is 15-deoxy-Δ12,14-PGJ2, a specific endogenous ligand for a peroxisome proliferator-activated receptor γ (PPARγ) in vivo [57]. PPARγ was detected in primary and secondary follicles, but was expressed highly in the large follicles, and is thought to play a role in regulating the expression of genes involved with growth, development, and/or differentiation of the follicle. Treatment of granulosa cells with 15-deoxy-Δ12,14-PGJ2 stimulated basal progesterone secretion but elicited no significant effect on FSH-stimulated steroid production, suggesting a potential role for PPARγ in regulating genes related to follicular differentiation [58]. Also, the PPARγ decline in response to LH is important for ovulation and/or luteinization [58]. Significantly higher levels of AA metabolites (such as PGE2, PGD2, PGJ2, 15-deoxy-Δ12,14-PGJ2, and PGF2α) produced via the COX pathway were found in FF of PCOS patients as compared to those of the non-PCOS group [38]. In the present study, the levels of AA metabolites (such as PGE2, PGD2, PGJ2, 15-deoxy-Δ12,14-PGJ2, and PGF2α) produced via the COX pathway in FF of DOR patients were lower than those in the NOR control group. However, the specific mechanism of DOR and the level of metabolites produced by the arachidonic acid COX enzyme metabolic pathway remain elusive.

Our analysis revealed15 types of differentially oxidized lipid metabolites in FF, possessing certain predictive value for ovarian reserve function. In particular, the levels of ±20-HDoHE, 12S-HHTrE, 20-COOH-AA, 8S,15S-DiHETE, PGA2, and PGE1 showed high sensitivities and specificities. Through Spearman correlation analysis, we also explored the correlations between the concentration of differentially oxidized lipid metabolites in FF and age, FSH, AMH, AFC, oocytes retrieved, MII oocytes, fertilization, high-quality embryos. Results demonstrated a negative correlation between the concentrations of eight lipid oxidative metabolites in FF with FSH but a positive correlation with AFC. Nine types of oxidized lipid metabolites were positively correlated with AMH, the number of oocytes retrieved, MII oocytes, and fertilization. Notably, one metabolite was positively correlated with the number of high-quality embryos. The changes in the content of oxidized lipid metabolites in FF of DOR patients suggest their potential role in oocyte development, maturation, fertilization, and early embryo development. However, it remains elusive as to how the above lipid-oxidizing metabolites play a role in FF and which metabolic regulation channels influence ovarian reserve function and oocyte quality, which warrants further exploration.

Intriguingly, our analysis showed that the differential oxidized lipid metabolites in FF of DOR patients were lower than those of NOR patients, which contradicts our original hypothesis of the search for markers of oxidized lipid metabolism in FF of DOR patients. The possible reasons may be that: First, our reproductive center routinely employed a double-cavity needle to wash oocytes, allowing DOR patients to obtain maximum levels of oocytes, which is not the case for NOR patients. As such, we insinuate the potential influence of needle washing liquid on oxidized lipid metabolites in FF. Second, the structure of some oxidized lipid metabolites in FF may be unstable; however, we made our best attempt to ensure that the processes of collecting, transporting, centrifuging, and shipping FF for detection have no underlying impact. It is suggested that the above process may still impact oxidized lipid metabolites in FF. Furthermore, because the limited amount of oxidized lipid metabolites that can be detected by oxidized lipid metabolomics at present, we speculate that the different oxidative lipid metabolites in FF of DOR and NOR patients are way much limited, to an extent that the up-regulated metabolites may not be detected. Finally, it is possible that we analyzed a small number of samples and the individual differences of clinical samples were large. Although patients were screened with strict adherence to the inclusion criteria, the influence of individual differences cannot be ignored. Therefore, the difference between DOR and NOR patients in the up-regulated oxidized lipid metabolites in FF is not significant. Another finding of this study is that the down-regulated oxidized lipid metabolites in FF of DOR patients are primarily linked to AA metabolism. Mounting evidence indicates that AA potentially induces and inhibits oxidative stress, for instance, it can increase calcium influx by activating the calcium channel on the cell membrane [59, 60], which subsequently activates NADPH oxidase, increasing ROS levels and induces oxidative stress damage [61]. AA can also activate the mRNA and protein levels of PPAR, ameliorate the activities of SOD, CAT, and GSH-Px, decrease the production of mitochondrial ROS to arrest oxidative stress injury [62]. These data strongly suggest a dominant role for AA and its metabolites in anti-oxidative stress in FF of DOR patients, which will be explored in our future study.

Conclusions

Through UHPLC-MS/MS, oxylipins metabolites changes in FF of the DOR and NOR patients were investigated. Metabonomic analysis of FF showed 15 differentially expressed oxylipins metabolites, associated with ovarian reserve function. The differentially oxidized lipid metabolites were lower in FF of DOR patients as compared to those of the NOR group. The pathway enrichment analysis demonstrated that the differentially oxidized lipid metabolites are mainly concentrated in the AA metabolic pathway, while AA is involved in the regulation of oocyte development and maturation, and its complex changes are closely related to follicular development. The differential oxylipins metabolites and AA metabolic pathway may play a significant role in DOR infertility. Dysfunctions in the metabolism of oxylipins related to follicular development in FF of DOR patients may provide potential detection avenues and therapeutic targets for promoting oocyte maturation. The findings may also provide a scientific basis for understanding the microenvironment of oocyte development and improving oocyte quality in patients with DOR.

Availability of data and materials

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

Abbreviations

DOR:

Diminished ovarian reserve

POF:

Premature ovarian failure

POR:

Poor ovarian responder

NOR:

Normal ovarian reserve

IVF:

In vitro fertilization

ICSI-ET:

intracytoplasmic sperm injection-embryo transfer

OS:

Oxidative stress

ROS:

Reactive oxygen species

MDA:

Malondialdehyde

FF:

Follicular fluid

AGE:

Advanced glycation end-products

AMH:

Anti-Mullerian hormone

FSH:

Follicle stimulating hormone

LH:

Luteinizing hormone

bAFC:

Basic antral follicle count

E2 :

Estrogen

P:

Progesterone

SWATH:

Sequential window acquisition of all theoretical fragment-ion spectra

MRM:

Multiple reaction monitor

PCA:

Principal component analysis

OPLS-DA:

Orthogonal projections to latent structure discriminant analysis

QC:

Quality control

VIP:

Variable importance in the projection

KEGG:

Kyoto Encyclopedia of Genes and Genomes

BMI:

Body mass index

HMG:

Human menopausal gonadotropin

HCG:

Human chorionic gonadotropin

SOD:

Superoxide dismutase

CAT:

Catalase

GSH-PX:

Glutathione peroxidase

AA:

Arachidonic acid

FDR:

False discovery rate

COX:

Cyclooxygenase

LOX:

Lipoxygenase

CYP:

Cytochrome P450

PUFAs:

Polyunsaturated fatty acids

PLA2:

Phospholipase A2

CGCs:

Cumulus granulosa cells

ASA:

Acetylsalicylic acid

20-HETE:

20-hydroxyeicosatetraenoic acid

IR:

Insulin resistance

PGs:

Prostaglandins

PTGS1:

Prostaglandin-endoperoxide synthase S1

PTGS2:

Prostaglandin-endoperoxide synthase S2

PPARγ:

Peroxisome proliferator-activated receptor γ

ROC:

Receiver operating characteristic

AUC:

Area under the curve

References

  1. Lu Q, Shen H, Li Y, Zhang C, Wang C, Chen X, et al. Low testosterone levels in women with diminished ovarian reserve impair embryo implantation rate: a retrospective case-control study. J Assist Reprod Genet. 2014;31(4):485–91.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Mutlu MF, Erdem A. Evaluation of ovarian reserve in infertile patients. J Turk Ger Gynecol Assoc. 2012;13(3):196–203.

    Article  PubMed  PubMed Central  Google Scholar 

  3. European Society for Human Reproduction and Embryology (ESHRE) Guideline Group on POI, Webber L, Davies M, Anderson R, Bartlett J, Braat D, et al. ESHRE Guideline: management of women with premature ovarian insufficiency. Hum Reprod. 2016;31(5):926–37.

    Article  Google Scholar 

  4. Bishop LA, Richter KS, Patounakis G, Andriani L, Moon K, Devine K. Diminished ovarian reserve as measured by means of baseline follicle-stimulating hormone and antral follicle count is not associated with pregnancy loss in younger in vitro fertilization patients. Fertil Steril. 2017;108(6):980–7.

    Article  CAS  PubMed  Google Scholar 

  5. Nita M, Grzybowski A. The Role of the Reactive Oxygen Species and Oxidative Stress in the Pathomechanism of the Age-Related Ocular Diseases and Other Pathologies of the Anterior and Posterior Eye Segments in Adults. Oxid Med Cell Longev. 2016;2016:3164734.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Wohlgemuth SE, Calvani R, Marzetti E. The interplay between autophagy and mitochondrial dysfunction in oxidative stress-induced cardiac aging and pathology. J Mol Cell Cardiol. 2014;71:62–70.

    Article  CAS  PubMed  Google Scholar 

  7. Yang Y, Sharma R, Sharma A, Awasthi S, Awasthi YC. Lipid peroxidation and cell cycle signaling: 4-hydroxynonenal, a key molecule in stress mediated signaling. Acta Biochim Pol. 2003;50(2):319–36.

    Article  CAS  PubMed  Google Scholar 

  8. Luderer U. Ovarian toxicity from reactive oxygen species. Vitam Horm. 2014;94:99–127.

    Article  CAS  PubMed  Google Scholar 

  9. Shamilova NN, Marchenko LA, Dolgushina NV, Zaletaev DV, Sukhikh GT. The role of genetic and autoimmune factors in premature ovarian failure. J Assist Reprod Genet. 2013;30(5):617–22.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Razi Y, Eftekhar M, Fesahat F, Dehghani Firouzabadi R, Razi N, Sabour M, et al. Concentrations of homocysteine in follicular fluid and embryo quality and oocyte maturity in infertile women: a prospective cohort. J Obstet Gynaecol. 2020;4:1–6.

    Google Scholar 

  11. Agarwal A, Aponte-Mellado A, Premkumar BJ, Shaman A, Gupta S. The effects of oxidative stress on female reproduction: a review. Reprod Biol Endocrinol. 2012;10:49.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Jia ZZ, Zhang JW, Zhou D, Xu DQ, Feng XZ. Deltamethrin exposure induces oxidative stress and affects meiotic maturation in mouse oocyte. Chemosphere. 2019;223:704–13.

    Article  CAS  PubMed  Google Scholar 

  13. Lord T, Aitken RJ. Oxidative stress and ageing of the post-ovulatory oocyte. Reproduction. 2013;146(6):R217–27.

    Article  CAS  PubMed  Google Scholar 

  14. Tian H, Liu K, Zhang Y, Qi Q, Wang C, Guan H, et al. Adult follicular fluid supplementation during in vitro maturation improves the developmental competence of prepubertal lamb oocytes. Theriogenology. 2019;130:157–62.

    Article  CAS  PubMed  Google Scholar 

  15. Sreekumar A, Poisson LM, Rajendiran TM, Khan AP, Cao Q, Yu J, et al. Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression. Nature. 2009;457(7231):910–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Merhi Z, Irani M, Doswell AD, Ambroggio J. Follicular fluid soluble receptor for advanced glycation end-products (sRAGE): a potential indicator of ovarian reserve. J Clin Endocrinol Metab. 2014;99(2):E226–33.

    Article  CAS  PubMed  Google Scholar 

  17. Pacella L, Zander-Fox DL, Armstrong DT, Lane M. Women with reduced ovarian reserve or advanced maternal age have an altered follicular environment. Fertil Steril. 2012;98(4):986–94 e1-2.

    Article  PubMed  Google Scholar 

  18. Tal R, Seifer DB. Ovarian reserve testing: a user's guide. Am J Obstet Gynecol. 2017;217(2):129–40.

    Article  PubMed  Google Scholar 

  19. Zhu X, Ye J, Fu Y. Late Follicular Phase Ovarian Stimulation Without Exogenous Pituitary Modulators. Front Endocrinol (Lausanne). 2020;11:487.

    Article  Google Scholar 

  20. Zhou SJ, Zhao MJ, Li C, Su X. The comparison of evaluative effectiveness between antral follicle count/age ratio and ovarian response prediction index for the ovarian reserve and response functions in infertile women. Medicine (Baltimore). 2020;99(36):e21979.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Gabbs M, Leng S, Devassy JG, Monirujjaman M, Aukema HM. Advances in Our Understanding of Oxylipins Derived from Dietary PUFAs. Adv Nutr. 2015;6(5):513–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Edin ML, Duval C, Zhang G, Zeldin DC. Role of linoleic acid-derived oxylipins in cancer. Cancer Metastasis Rev. 2020;39(3):581–2.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Nayeem MA. Role of oxylipins in cardiovascular diseases. Acta Pharmacol Sin. 2018;39(7):1142–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Tans R, Bande R, van Rooij A, Molloy BJ, Stienstra R, Tack CJ, et al. Evaluation of cyclooxygenase oxylipins as potential biomarker for obesity-associated adipose tissue inflammation and type 2 diabetes using targeted multiple reaction monitoring mass spectrometry. Prostaglandins Leukot Essent Fatty Acids. 2020;160:102157.

    Article  CAS  PubMed  Google Scholar 

  25. Lundström SL, Balgoma D, Wheelock ÅM, Haeggström JZ, Dahlén SE, Wheelock CE. Lipid mediator profiling in pulmonary disease. Curr Pharm Biotechnol. 2011;12(7):1026–52.

    Article  PubMed  Google Scholar 

  26. Devassy JG, Leng S, Gabbs M, Monirujjaman M, Aukema HM. Omega-3 Polyunsaturated Fatty Acids and Oxylipins in Neuroinflammation and Management of Alzheimer Disease. Adv Nutr. 2016;7(5):905–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xu H, Dong X, Ai Q, Mai K, Xu W, Zhang Y, et al. Regulation of tissue LC-PUFA contents, Δ6 fatty acyl desaturase (FADS2) gene expression and the methylation of the putative FADS2 gene promoter by different dietary fatty acid profiles in Japanese seabass (Lateolabrax japonicus). PLoS One. 2014;9(1):e87726.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Wathes DC, Abayasekara DR, Aitken RJ. Polyunsaturated fatty acids in male and female reproduction. Biol Reprod. 2007;77(2):190–201.

    Article  CAS  PubMed  Google Scholar 

  29. Hanna VS, Hafez EAA. Synopsis of arachidonic acid metabolism: A review. J Adv Res. 2018;11:23–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Ciepiela P, Bączkowski T, Drozd A, Kazienko A, Stachowska E, Kurzawa R. Arachidonic and linoleic acid derivatives impact oocyte ICSI fertilization--a prospective analysis of follicular fluid and a matched oocyte in a 'one follicle--one retrieved oocyte--one resulting embryo' investigational setting. PLoS One. 2015;10(3):e0119087.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Cheng Z, Sheldrick EL, Marshall E, Wathes DC, Abayasekara DR, Flint AP. Control of cyclic AMP concentration in bovine endometrial stromal cells by arachidonic acid. Reproduction. 2007;133(5):1017–26.

    Article  CAS  PubMed  Google Scholar 

  32. Zuccolo E, Dragoni S, Poletto V, Catarsi P, Guido D, Rappa A, et al. Arachidonic acid-evoked Ca2+ signals promote nitric oxide release and proliferation in human endothelial colony forming cells. Vascul Pharmacol. 2016;87:159–71.

    Article  CAS  PubMed  Google Scholar 

  33. Lapa M, Marques CC, Alves SP, Vasques MI, Baptista MC, Carvalhais I, et al. Effect of trans-10 cis-12 conjugated linoleic acid on bovine oocyte competence and fatty acid composition. Reprod Domest Anim. 2011;46(5):904–10.

    Article  CAS  PubMed  Google Scholar 

  34. Renaville B, Bacciu N, Comin A, Motta M, Poli I, Vanini G, et al. Plasma and follicular fluid fatty acid profiles in dairy cows. Reprod Domest Anim. 2010;45(1):118–21.

    Article  CAS  PubMed  Google Scholar 

  35. Moran LJ, Tsagareli V, Noakes M, Norman R. Altered Preconception Fatty Acid Intake Is Associated with Improved Pregnancy Rates in Overweight and Obese Women Undertaking in Vitro Fertilisation. Nutrients. 2016;8(1):10.

    Article  PubMed Central  CAS  Google Scholar 

  36. Zhang N, Wang L, Luo G, Tang X, Ma L, Zheng Y, et al. Arachidonic Acid Regulation of Intracellular Signaling Pathways and Target Gene Expression in Bovine Ovarian Granulosa Cells. Animals (Basel). 2019;9(6):374.

    Article  PubMed Central  Google Scholar 

  37. Khajeh M, Nouri M, Ghasemzadeh A, Mehdizadeh A, Shanehbandi D, Yousefi S, et al. Arachidonic acid alleviates the detrimental effects of acetylsalicylic acid on human granulosa cells performance in vitro. Mol Reprod Dev. 2020;87(5):607–19.

    Article  CAS  PubMed  Google Scholar 

  38. Li S, Qi J, Tao Y, Zhu Q, Huang R, Liao Y, et al. Elevated Levels of Arachidonic Acid Metabolites in Follicular Fluid of PCOS Patients. Reproduction. 2019; REP-19-0136.R2.

  39. Kroetz DL, Xu F. Regulation and inhibition of arachidonic acid omega-hydroxylases and 20-HETE formation. Annu Rev Pharmacol Toxicol. 2005;45:413–38.

    Article  CAS  PubMed  Google Scholar 

  40. Peterson SJ, Vanella L, Gotlinger K, Jiang H, Singh SP, Sodhi K, et al. Oxidized HDL is a potent inducer of adipogenesis and causes activation of the Ang-II and 20-HETE systems in human obese females. Prostaglandins Other Lipid Mediat. 2016;123:68–77.

    Article  CAS  PubMed  Google Scholar 

  41. Tsai IJ, Croft KD, Mori TA, Falck JR, Beilin LJ, Puddey IB, et al. 20-HETE and F2-isoprostanes in the metabolic syndrome: the effect of weight reduction. Free Radic Biol Med. 2009;46(2):263–70.

    Article  CAS  PubMed  Google Scholar 

  42. Issan Y, Hochhauser E, Guo A, Gotlinger KH, Kornowski R, Leshem-Lev D, et al. Elevated level of pro-inflammatory eicosanoids and EPC dysfunction in diabetic patients with cardiac ischemia. Prostaglandins Other Lipid Mediat. 2013;100-101:15–21.

    Article  CAS  PubMed  Google Scholar 

  43. Gilani A, Pandey V, Garcia V, Agostinucci K, Singh SP, Schragenheim J, et al. High-fat diet-induced obesity and insulin resistance in CYP4a14-/- mice is mediated by 20-HETE. Am J Physiol Regul Integr Comp Physiol. 2018;315(5):R934–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kudo I, Murakami M. Prostaglandin E synthase, a terminal enzyme for prostaglandin E2 biosynthesis. J Biochem Mol Biol. 2005;38(6):633–8.

    CAS  PubMed  Google Scholar 

  45. Asselin E, Lacroix D, Fortier MA. IFN-tau increases PGE2 production and COX-2 gene expression in the bovine endometrium in vitro. Mol Cell Endocrinol. 1997;132(1-2):117–26.

    Article  CAS  PubMed  Google Scholar 

  46. Hester KE, Harper MJ, Duffy DM. Oral administration of the cyclooxygenase-2 (COX-2) inhibitor meloxicam blocks ovulation in non-human primates when administered to simulate emergency contraception. Hum Reprod. 2010;25(2):360–7.

    Article  CAS  PubMed  Google Scholar 

  47. Lim H, Paria BC, Das SK, Dinchuk JE, Langenbach R, Trzaskos JM, et al. Multiple female reproductive failures in cyclooxygenase 2-deficient mice. Cell. 1997;91(2):197–208.

    Article  CAS  PubMed  Google Scholar 

  48. Takahashi T, Morrow JD, Wang H, Dey SK. Cyclooxygenase-2-derived prostaglandin E(2) directs oocyte maturation by differentially influencing multiple signaling pathways. J Biol Chem. 2006;281(48):37117–29.

    Article  CAS  PubMed  Google Scholar 

  49. Marei WF, Abayasekara DR, Wathes DC, Fouladi-Nashta AA. Role of PTGS2-generated PGE2 during gonadotrophin-induced bovine oocyte maturation and cumulus cell expansion. Reprod Biomed Online. 2014;28(3):388–400.

    Article  CAS  PubMed  Google Scholar 

  50. Boruszewska D, Kowalczyk-Zieba I, Suwik K, Staszkiewicz-Chodor J, Jaworska J, Lukaszuk K, et al. Prostaglandin E2 affects in vitro maturation of bovine oocytes. Reprod Biol Endocrinol. 2020;18(1):40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Duffy DM, Ko C, Jo M, Brannstrom M, Curry TE. Ovulation: Parallels With Inflammatory Processes. Endocr Rev. 2019;40(2):369–416.

    Article  PubMed  Google Scholar 

  52. Algire JE, Srikandakumar A, Guilbault LA, Downey BR. Preovulatory changes in follicular prostaglandins and their role in ovulation in cattle. Can J Vet Res. 1992;56(1):67–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Berisha B, Rodler D, Schams D, Sinowatz F, Pfaffl MW. Prostaglandins in Superovulation Induced Bovine Follicles During the Preovulatory Period and Early Corpus Luteum. Front Endocrinol (Lausanne). 2019;10:467.

    Article  Google Scholar 

  54. Pereira de Moraes F, Amaral D’Avila C, Caetano de Oliveira F, Ávila de Castro N, Diniz Vieira A, Schneider A, et al. Prostaglandin F2α regulation and function during ovulation and luteinization in cows. Theriogenology. 2021;171:30–7.

    Article  CAS  PubMed  Google Scholar 

  55. Sharma AK, Sharma RK. Effect of prostaglandins E2 and F2α on granulosa cell apoptosis in goat ovarian follicles. Iran J Vet Res. 2020;21(2):97–102.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Kemiläinen H, Adam M, Mäki-Jouppila J, Damdimopoulou P, Damdimopoulos AE, Kere J, et al. The Hydroxysteroid (17β) Dehydrogenase Family Gene HSD17B12 Is Involved in the Prostaglandin Synthesis Pathway, the Ovarian Function, and Regulation of Fertility. Endocrinology. 2016;157(10):3719–30.

    Article  PubMed  CAS  Google Scholar 

  57. Forman BM, Tontonoz P, Chen J, Brun RP, Spiegelman BM, Evans RM. 15-Deoxy-delta 12, 14-prostaglandin J2 is a ligand for the adipocyte determination factor PPAR gamma. Cell. 1995;83(5):803–12.

    Article  CAS  PubMed  Google Scholar 

  58. Komar CM, Curry TE Jr. Inverse relationship between the expression of messenger ribonucleic acid for peroxisome proliferator-activated receptor gamma and P450 side chain cleavage in the rat ovary. Biol Reprod. 2003;69(2):549–55.

    Article  CAS  PubMed  Google Scholar 

  59. Fiorio Pla A, Genova T, Pupo E, Tomatis C, Genazzani A, Zaninetti R, et al. Multiple roles of protein kinase a in arachidonic acid-mediated Ca2+ entry and tumor-derived human endothelial cell migration. Mol Cancer Res. 2010;8(11):1466–76.

    Article  CAS  PubMed  Google Scholar 

  60. Thompson J, Mignen O, Shuttleworth TJ. The N-terminal domain of Orai3 determines selectivity for activation of the store-independent ARC channel by arachidonic acid. Channels (Austin). 2010;4(5):398–410.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Tseng CL, Wei JW. Homologous desensitization of histamine-mediated signal transduction system in C6 glioma cells. Chin J Physiol. 2013;56(2):90–100.

    CAS  PubMed  Google Scholar 

  62. Wang ZJ, Liang CL, Li GM, Yu CY, Yin M. Neuroprotective effects of arachidonic acid against oxidative stress on rat hippocampal slices. Chem Biol Interact. 2006;163(3):207–17.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are thankful to all the participants of the study. We would like to thank Shanghai Biotree Biotechnology for the UHPLC-MS/MS and data analyses.

Funding

This work was supported by the Key Project of the National Key Research and Development Program of "Modernization Research of Traditional Chinese Medicine" -Systematic Study on the Criteria for Differentiation of "Kidney-Yang Deficiency Syndrome" (No. SQ2018YFC170061)

Author information

Authors and Affiliations

  1. Shanghai University of Traditional Chinese Medicine, Shanghai, 201203, China

    Chengcheng Liang

  2. Department of Gynecology, Shuguang Hospital affiliated to Shanghai University of Traditional Chinese Medicine, 201203, Shanghai, China

    Chengcheng Liang, Xiaole Zhang, Cong Qi & Hui Hu

  3. Shanghai First Maternity and Infant Hospital, Tongji University School of Medicine, Shanghai, 201204, China

    Qinhua Zhang, Xiuxian Zhu & Yonglun Fu

Authors
  1. Chengcheng Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaole Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Cong Qi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Hui Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qinhua Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xiuxian Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yonglun Fu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LCC performed the literature search, data analysis and interpretation, and wrote the manuscript. ZXL and HH constructed the tables and figures. LCC and ZXX performed the sample collection and data analysis. ZQH, QC and FYL contributed to the design of the research study and provided critical review of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qinhua Zhang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Research Ethics Committees of Shuguang hospital affiliated to Shanghai university of traditional Chinese medicine (No.2020-833-40-01). All participants signed an informed consent form prior to the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons 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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Liang, C., Zhang, X., Qi, C. et al. UHPLC-MS-MS analysis of oxylipins metabolomics components of follicular fluid in infertile individuals with diminished ovarian reserve. Reprod Biol Endocrinol 19, 143 (2021). https://doi.org/10.1186/s12958-021-00825-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12958-021-00825-x

Keywords

Reproductive Biology and Endocrinology

ISSN: 1477-7827

Contact us