JBRA Assist. Reprod. 2019;23(4):323-327
ORIGINAL ARTICLE
doi: 10.5935/1518-0557.20190008
Serum metabolites as predictive molecular markers of ovarian response to controlled stimulation: a pilot study
1Fertility Medical Group. São Paulo, SP - Brazil
2Instituto Sapientiae - Centro de Estudos e Pesquisa em Reprodução Assistida. São Paulo, SP - Brazil
3Grupo de Bio-orgânica e Bioanalítica. Departamento de Química - Instituto de Ciências Ambientais, Químicas e Farmacêuticas - UNIFESP Diadema Campus, Diadema, SP - Brazil
CONFLICT OF INTEREST
The authors have no conflict of interest to declare.
ABSTRACT
Objective: This study aimed to look into the use of serum metabolites as potential biomarkers of response to controlled ovarian stimulation (COS) in patients undergoing intracytoplasmic sperm injection (ICSI) cycles.
Methods: This case-control study analyzed serum samples from 30 patients aged <36 years undergoing COS for ICSI in a university-affiliated assisted reproduction center from January 2017 to August 2017. The samples were split into three groups based on response to COS as follows: poor responders: <4 retrieved oocytes (PR group, n=10); normal responders: ≥ 8 and ≤ 12 retrieved oocytes (NR group, n=10); and hyper-responders: >25 retrieved oocytes (HR, n=10). The metabolic profiles of the serum samples were compared between the groups through Principal Component Analysis (PCA). Receiver Operating Characteristic (ROC) curves were built to assess the power of the model at predicting response to COS.
Results: PCA clearly distinguished between PR, NR and HR, and 10 ions were chosen as potential biomarkers of response to COS. These ions were more specific for PR than for NR. The ROC curve considering PR and NR had an area under the curve of 99.6% (95% CI: 88.9 - 100%).
Conclusion: The preliminary evidence discussed in this study suggests that serum metabolites may be used as predictive molecular markers of ovarian response to controlled stimulation. The integration of clinical and “omics” findings may allow the migration toward an era of personalized treatment in reproductive medicine.
Keywords: biomarker, controlled ovarian stimulation, ICSI, IVF, metabolomics
INTRODUCTION
Controlled ovarian stimulation (COS) is a key step of in vitro fertilization used to stimulate the development of multiple follicles, enable the retrieval of multiple oocytes, and thus allow the selection of the best embryo for transfer. The retrieval of 5 -14 oocytes has been defined as a clinically appropriate ovarian response to COS (Kyrou et al., 2009), but poor and excessive response to COS are not rare.
The balance between retrieving too few or too many oocytes is difficult to strike, and may affect the chances of taking a baby back home and even increase the risk of patients suffering negative health outcomes through events such as the ovarian hyperstimulation syndrome (OHSS) (Kumbak et al., 2009).
Nevertheless, it is nearly impossible to accurately predict ovarian response and tailor stimulation protocols based on parameters such as maternal age, antral follicle count, anti-Müllerian hormone level, or prior poor/excessive response to COS. The development of noninvasive techniques to predict response to COS would allow individualized treatment, significantly increase treatment success rates, and alleviate the physical, emotional, and economic burden faced by patients.
Metabolomics is a great tool for the comprehensive study of the dynamic changes of the metabolome, and provides a powerful platform to discover biomarkers and improve diagnostic and therapeutic monitoring (Eckhart et al., 2012; Wang et al., 2013; Hertel et al., 2016). In the field of assisted reproduction, most of the studies have focused on the embryonic metabolome through the analysis of spent culture media (Botros et al., 2008; Cortezzi et al., 2013). In addition, follicular fluid metabolites have also been studied to better understand oocyte competence (O'Gorman et al., 2013), ovarian aging (de la Barca et al., 2017), endometriosis (Santonastaso et al., 2017), and polycystic ovarian syndrome (PCOS) (Zhang et al., 2017). To date, there are no studies investigating the correlation between metabolomic profiling and response to COS.
This pilot study aimed to look into the use of serum metabolites as potential biomarkers of response to controlled ovarian stimulation (COS) in patients undergoing intracytoplasmic sperm injection (ICSI) cycles.
MATERIAL AND METHODS
Experimental Design, Inclusion and Exclusion criteria
This case-control study analyzed serum samples from 30 patients aged <36 years undergoing COS for ICSI in a university-affiliated assisted reproduction center from January 2017 to August 2017. The samples were split into three groups based on response to COS as follows: poor responders: <4 retrieved oocytes (PR group, n=10); normal responders: ≥8 and ≤12 retrieved oocytes (NR group, n=10); and hyper-responders: >25 retrieved oocytes (HR, n=10). The metabolic profiles of the serum samples were compared between the groups. The patients gave informed consent prior to joining the study. The local institutional review board approved the study.
Controlled Ovarian Stimulation and Laboratory Procedures
The patients were prescribed recombinant FSH (Gonal-F®, Merck KGaA, Darmstadt, Germany) for controlled ovarian stimulation and a GnRH antagonist (GnRH - Cetrotide® Merck KGaA, Darmstadt, Germany) for pituitary suppression. Follicular growth was monitored using transvaginal ultrasound examination starting on day 4 of gonadotropin administration. When adequate follicular growth and serum E2 levels were observed, recombinant hCG (Ovidrel®, Merck KGaA, Darmstadt, Germany) was administered to trigger the final follicular maturation. The oocytes were collected 35 hours after hCG administration through transvaginal ultrasound ovum pick-up. The retrieved oocytes were assessed to determine their nuclear status, and the ones in metaphase II were submitted to ICSI following routine procedures (Palermo et al., 1997).
Sample Preparation, Metabolite Extraction, and Mass Spectrometry
Metabolites were extracted on ice for protein precipitation using Methanol/Chloroform (2:1 v/v). The mixture was centrifuged at 8,000×g for 15 min. After centrifugation, the supernatant was collected and transferred to a 96-well plate, which was placed inside a microTOF-QII™ mass spectrometer equipped with an Apollo II electrospray ion source (Bruker, Billerica, USA) coupled to a UFLC Prominence binary liquid chromatograph (Shimadzu, Kyoto, Japan). The plate was stored into the SIL-30AC autosampler at 10ºC prior to analysis. The samples were directly injected into the analyzer in a 1 uL volume, by a 20 mmol/L ammonium formate solution in acetonitrile/2-propanol (4:1, v/v) at a flow rate of 200 µL/min. Spectra were acquired in the positive mode using a range of m/z 50-1200 Da. Sodium formate clusters in isopropyl alcohol within the m/z 50-1200 Da range were used as the calibration standard.
Data processing and statistical analysis
Patient and cycle characteristics were analyzed on SPSS Statistics 21 (IBM, New York, NY, USA). Variables were tested for normality and group homogeneity using the Shapiro-Wilk and Levenne tests, respectively. When needed, the samples were standardized using the z-score. The variables were compared between groups through one-way ANOVA, followed by the Bonferroni post-hoc test. Variables were described as mean values ± standard deviation and significance was attributed when α was 5%. Mass spectrometry results were obtained using the DataAnalysis 4.1 sofware (Bruker Daltonics Bremen, Germany) and data analyses were performed on MetaboAnalyst 3.0 (http://www.metaboanalyst.ca). Log2 was used to normalize intensity values followed by self-scaling. Principal component analysis (PCA), an unsupervised method, was applied to the data set to detect intrinsic clusters based on metabolic profiles. From these analyses, a list of ions responsible for group discrimination was obtained. These ions were used to build a receiver operating characteristic (ROC) curve and to evaluate the strength of the model at predicting response to COS. Metabolite attribution was performed based on the Human metabolites database (http://www.hmdb.ca/), with a maximum mass tolerance of 0.01 Da. The maximum mass error was 30 ppm. For the attribution, only molecules containing hydrogen (M+H+), sodium (M+Na+) and potassium (M+K+) as adducts were considered.
RESULTS
Patient and cycle characteristics
The patient and cycle characteristics are described in Table 1. As expected, the level of estradiol on hCG trigger day, the number of aspirated follicles, the number of retrieved follicles, and the number of mature follicles were higher in the HR group followed by the NR group, while the PR group presented the lowest results.
Table 1. Patient and cycle characteristics for the poor, normal, and hyper-responder groups
Metabolomic analysis
Considering components 1 and 2, PCA clearly distinguished between the PR, NR, and HR groups (Figure 1).
Figure 1. Variance among groups according to principal component analysis (PCA score plot)
The ion masses associated to the separation of the groups were obtained by the PCA loading plot, from which 10 ions were chosen as potential biomarkers (Table 2).
Table 2. Tentative biomarker attribution, loading values, and predictive values for high and poor responders
The relative abundances of individual biomarkers are demonstrated in the box plot charts (Figure 2). The PCA loading values for the selected ions based on principal component 1 are presented in Table 2.
Figure 2. Boxes plot charts for the relative abundances of individual biomarkers
Our evidence demonstrated that all ions selected in the present study were more specific for the PR group when compared to the NR group. The ROC curve considering the PR and NR groups presented an area under the curve (AUC) of 99.6% (95% CI: 88.9 - 100%, Figure 3).
Figure 3. ROC curve and predicted class probabilities for normal and poor responders
DISCUSSION
The success of assisted reproductive technology (ART) treatments is highly dependent on the response to COS. Moreover, considering the physical, emotional, and economic burden faced by patients undergoing ART treatment, finding a method to predict ovarian response to COS would be a major step forward in the field of reproductive medicine.
This study was able to identify metabolites that might be used as predictive molecular markers of response to COS. PCA clearly distinguished between PR, NR and HR groups, and 10 ions were chosen as potential biomarkers of response to COS. To date, this is the first study to investigate whether blood plasma metabolites might predict ovarian response to COS. The advantage of the method used in this study is having knowledge of the patients’ ability to respond to gonadotropins before the start of stimulation. Previous metabolomic studies in assisted reproduction investigated the metabolic profile on seminal plasma (Deepinder et al., 2007) and fluid (Gupta et al., 2011), embryo spent culture media (Cortezzi et al., 2013), follicular (O'Gorman et al., 2013) and endometrial fluid (Braga et al., 2017).
OMICS technologies study cellular events and interactions from deoxyribonucleic acid (DNA) and genes to metabolites in a global way. Although metabolomics was recognized as a separate area of science much later than the other “omics” - such as genomics, transcriptomics, and proteomics - it provides a powerful platform for the discovery of novel biomarkers. The main advantages of metabolomics are its familiarity to the actual phenotype and the number of possible low molecular weight bio-compounds (Yoshida et al., 2012).
In conclusion, our preliminary evidence suggests that serum metabolites might work as predictive molecular markers of ovarian response to controlled stimulation. We believe that the quality of emerging data will eventually allow the individualization of COS, and further studies will validate new biomarkers of ovarian response. To date, the technology and software around metabolomics are still developing as the human metabolome is being mapped. The integration of clinical and “omics” findings will eventually allow the migration toward an era of personalized treatment in the field of reproductive medicine.
REFERENCES
Botros L, Sakkas D, Seli E. Metabolomics and its application for non-invasive embryo assessment in IVF. Mol Hum Reprod. 2008;14:679-90.
Medline Crossref
Braga DP, Montani DA, Setti AS, Pilli G, Godoy A, Eberlin MN, Iaconelli A Jr., Borges E Jr., Lo Turco EG. Uterine fluid lipidomic as an endometrial receptivity predictive tool. Fertil Steril. 2017;108:e358.
Crossref
Cortezzi SS, Cabral EC, Trevisan MG, Ferreira CR, Setti AS, Braga DP, Figueira Rde C, Iaconelli A Jr., Eberlin MN, Borges E Jr. Prediction of embryo implantation potential by mass spectrometry fingerprinting of the culture medium. Reproduction. 2013;145:453-62.
Medline Crossref
de la Barca JMC, Boueilh T, Simard G, Boucret L, Ferre-L'Hotellier V, Tessier L, Gadras C, Bouet PE, Descamps P, Procaccio V, Reynier P, May-Panloup P. Targeted metabolomics reveals reduced levels of polyunsaturated choline plasmalogens and a smaller dimethylarginine/arginine ratio in the follicular fluid of patients with a diminished ovarian reserve. Hum Reprod. 2017;32:2269-78.
Medline Crossref
Deepinder F, Chowdary HT, Agarwal A. Role of metabolomic analysis of biomarkers in the management of male infertility. Expert Rev Mol Diagn. 2007;7:351-8.
Medline Crossref
Eckhart AD, Beebe K, Milburn M. Metabolomics as a key integrator for "omic" advancement of personalized medicine and future therapies. Clin Transl Sci. 2012;5:285-8.
Medline Crossref
Gupta A, Mahdi AA, Ahmad MK, Shukla KK, Jaiswer SP, Shankhwar SN. 1H NMR spectroscopic studies on human seminal plasma: a probative discriminant function analysis classification model. J Pharm Biomed Anal. 2011;54:106-13.
Medline Crossref
Hertel J, Friedrich N, Wittfeld K, Pietzner M, Budde K, Van der Auwera S, Lohmann T, Teumer A, Völzke H, Nauck M, Grabe HJ. Measuring Biological Age via Metabonomics: The Metabolic Age Score. J Proteome Res. 2016;15:400-10.
Medline Crossref
Kumbak B, Ulug U, Erzik B, Akbas H, Bahceci M. Early clinical pregnancy loss rate in poor responder patients does not change compared to age-matched normoresponders. Fertil Steril. 2009;91:106-9.
Medline Crossref
Kyrou D, Popovic-Todorovic B, Fatemi HM, Bourgain C, Haentjens P, Van Landuyt L, Devroey P. Does the estradiol level on the day of human chorionic gonadotrophin administration have an impact on pregnancy rates in patients treated with rec-FSH/GnRH antagonist? Hum Reprod. 2009;24:2902-9.
Medline Crossref
O'Gorman A, Wallace M, Cottell E, Gibney MJ, McAuliffe FM, Wingfield M, Brennan L. Metabolic profiling of human follicular fluid identifies potential biomarkers of oocyte developmental competence. Reproduction. 2013;146:389-95.
Medline Crossref
Palermo GD, Colombero LT, Rosenwaks Z. The human sperm centrosome is responsible for normal syngamy and early embryonic development. Rev Reprod. 1997;2:19-27.
Medline Crossref
Santonastaso M, Pucciarelli A, Costantini S, Caprio F, Sorice A, Capone F, Natella A, Iardino P, Colacurci N, Chiosi E. Metabolomic profiling and biochemical evaluation of the follicular fluid of endometriosis patients. Mol Biosyst. 2017;13:1213-22.
Medline Crossref
Wang X, Zhang A, Sun H. Power of metabolomics in diagnosis and biomarker discovery of hepatocellular carcinoma. Hepatology. 2013;57:2072-7.
Medline Crossref
Yoshida M, Hatano N, Nishiumi S, Irino Y, Izumi Y, Takenawa T, Azuma T. Diagnosis of gastroenterological diseases by metabolome analysis using gas chromatography-mass spectrometry. J Gastroenterol. 2012;47:9-20.
Medline Crossref
Zhang Y, Liu L, Yin TL, Yang J, Xiong CL. Follicular metabolic changes and effects on oocyte quality in polycystic ovary syndrome patients. Oncotarget. 2017;8:80472-80.
Medline Crossref