Section 2
This prospective observational study was carried out in the Division of Gynecology and Obstetrics, Department of Surgical Sciences, University of Cagliari, Cagliari, Italy, between November 2015 and June 2016. The trial was registered with the number NCT02337816 and approved by the Duilio Casula University Hospital Ethics Committee, Monserrato, Cagliari (ENDOMETAB01-Prot. 2015/3649).
The manuscript was prepared following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement [ 44 ].
Fifty-four consecutive women, ranging in age from 18 to 50 years (mean age 58.7), who needed surgery for suspected endometriosis or other clinical conditions were enrolled in the study. Exclusion criteria were assessed as follows: hormonal treatment two months before surgery, menopausal status, pregnancy, gynecological cancers, and pelvic inflammatory disease. Before joining the study, each woman provided written informed consent, and the patient’s anonymity was preserved. No study advertising was made, and no remuneration was offered. None of the enrolled patients suffered from autoimmune diseases, nor did they have symptoms referable to potential autoimmune diseases. Clinical data were collected, including socio-demographic data from patients, symptoms, comorbidities, previous therapies, and ultrasound data.
8–12 h of fasting were observed before surgery, and all women were treated with an osmotic laxative (Isocolan ® , Giuliani S.p.a., Milano, Italy). Before surgery, each patient’s blood sample was collected in the middle of the follicular phase and sent to the laboratory within twenty minutes. Serum samples were stored at −80 ° C until use. After surgery, we also collected a surgical report to verify the diagnosis. Based on the diagnosis, ten patients were excluded due to the presence of neoplastic or inflammatory pathologies, and the remaining patients were divided into four groups: (1) 22 patients with symptoms and with a histological diagnosis of endometriosis (group E); (2) 10 asymptomatic patients without a diagnosis of endometriosis with uterine fibroids or pelvic organ prolapse (group NE); (3) 4 patients with symptoms but without a diagnosis of endometriosis (4 patients); (4) 3 asymptomatic patients with a diagnosis of endometriosis. The third and fourth groups were not considered for metabolomic analyses due to the low number of patients.
The American Society for Reproductive Medicine classification was used to assess the severity of endometriosis [ 45 ].
Considering group E, all patients were classified at stage IV as severe endometriosis and presented the following symptoms before surgery: dysmenorrhea, chronic pelvic pain, dyspareunia, dyschezia, and dysuria. Metabolomic analysis was performed, including groups E and NE, so 37 patients were recruited. Table 1 and Table 2 show all demographic data.
A total of nineteen patients (11 E, 8 NE) showed the presence of comorbidity, nineteen patients (14E and 5 NE) reported previous intake of hormone therapies, thirteen patients (12 E, 1 NE) had ultrasound suggestive of endometriosis, and 22 patients (all belonging to the E class) had a surgical diagnosis of endometriosis.
The samples were prepared as described in a previous study [ 22 ]. Sera were centrifuged at 4500 rpm for 10 min at 4 °C, and 400 µL of supernatant were transferred to an Eppendorf tube. Hydrophilic and lipophilic metabolites were extracted and separated using a modified Folch method [ 46 ]. Four hundred µL of each serum sample was mixed with 600 µL of methanol containing an internal standard of Succinic acid-2,2,3,3-d4 (Sigma-Aldrich, St. Louis, MO, USA), 600 µL of chloroform, and 175 µL of distilled water. Then, they were centrifuged at 4500 rpm for 20 min at 4 °C. The lipid chloroform phase and the water/methanol phase have been separated.
A total of 150 µL of hydrophilic phase was dried with speed-vac ON and derivatized as follows: 50 μL of a methoxyamine solution in pyridine (10 mg/mL) (Sigma-Aldrich, St. Louis, MO, USA) was added, and after one hour at 70 °C, 100 μL of N-Methyl-N- (trimethylsilyl)-trifluoroacetamide, MSTFA (Sigma-Aldrich, St. Louis, MO, USA) was also added and left at room temperature for 1 h. The samples were resuspended with 150 μL of hexane (Sigma-Aldrich, St. Louis, MO, USA). A microliter of the derivatized sample was injected split-free into a 7890A gas chromatograph coupled with a 5975C Network mass spectrometer (Agilent Technologies, Santa Clara, CA, USA) equipped with a 30 m × 0.25 mm ID, fused silica capillary column of 30 m 0.25 mm ID, which was chemically bound with 0.25 μM TG-5MS stationary phase 0.25 M TG-5MS (Thermo Fisher Scientific, Waltham, MA, USA). The injector temperature was 250 °C. The gas flow through the column was 1 mL/min. The transfer line temperature was 280 °C. The initial temperature of the column was kept at 50 °C for 3 min, then increased from 50 °C to 250 °C at 10 °C/min and maintained at 250 °C for 12 min. After instrumental analysis, the identification of the metabolites was carried out using the standard NIST 08 ( http://www.nist.gov/srd/mslist.cfm (accessed on 16 November 2022)), Fiehn 2013 ( http://fiehnlab.ucdavis.edu/Metabolite-Library-2007 (accessed on 16 November 2022)) and GMD ( http://gmd.mpimp-golm.mpg.de (accessed on 16 November 2022)) mass spectra libraries (match ≥ 40%) and when available, by comparison with authentic standards. Data processing was performed using a Knime pipeline [ 47 ]. In summary, peak detection and deconvolution were performed in an R-XCMS package, filtering using blank samples and keeping features present in ≥50% of the samples. Missing value imputation was performed with the random forest algorithm. Concentrations of discriminant metabolites were obtained by the chromatogram area and then normalized by the total area (=100).
SIMCA-P software (ver. 16.0, Umetrics, Sweden) was used to perform multivariate statistical analysis. First, variables were UV-scaled, and principal component analysis (PCA) was applied. It is essential to explore sample distributions without a priori classification and identify potential outliers through the DmodX and Hotelling T 2 tests.
Then, Partial Least Squares Discriminant Analysis (PLS-DA) was applied. It maximizes the discrimination between samples assigned to different classes [ 48 ]. Statistical parameters such as variance and predictive ability (R 2 X, R 2 Y, Q 2 ) were established to evaluate the suitability of the models. R 2 and Q 2 can vary between 0 and 1 values. R 2 indicates the proportion of variance in the data explained by the model and indicates goodness to fit. Q 2 indicates the proportion of variance in the data predictable by the model and indicates the goodness of prediction. A good prediction model is achieved when Q 2 > 0.5 [ 48 ].
Furthermore, a permutation test (n = 400) was performed to validate the model. The PLS-DA model scores were tested using CV-ANOVA (ANalysis Of Variance), testing of Cross-Validated predictive residuals, to verify significance ( p < 0.05). CV-ANOVA is a diagnostic tool for assessing the reliability of the models [ 49 ]. The discriminant metabolites were selected by evaluating the loading plot and the VIP value. VIPs of more than 1 are the most relevant terms to explain Y (two classes assignment).
Univariate statistical analyses of the data were performed using GraphPad Prism software (version 7.01, GraphPad Software, Inc., San Diego, CA, USA). The Mann–Whitney U test was conducted to verify the significance of the metabolic results obtained from multivariate statistical analysis. ROC curves are designed to test the sensitivity and specificity of a selected metabolite pool.
MetaboAnalyst 5.0 ( www.metaboanalyst.ca (accessed on 16 November 2022)), a web server designed to obtain comprehensive analysis, visualization, and interpretation of metabolomic data [ 50 ], was used to generate metabolic pathways. This method allows a correlation between metabolite changes and metabolic networks.
Intro
Endometriosis is a chronic inflammatory disease characterized by the implantation of the stromal and/or endometrial glandular epithelium outside the uterus [ 1 ]. The ectopic endometrium is receptive to the cyclical fluctuation of ovarian steroids in terms of proliferation, differentiation, and bleeding [ 2 , 3 , 4 ].
It is estimated that endometriosis affects between 10 and 20 percent of women of reproductive age [ 5 ], causing dysmenorrhea, dyspareunia, noncyclical pelvic pain, infertility, and consequently, impairment of quality of life, sexual health, and psychological well-being [ 6 , 7 , 8 , 9 , 10 ].
Several hypotheses on its pathogenesis have been postulated: genetic, hormonal, environmental, and immunological factors have been contemplated; however, the exact role of each etiological factor remains unclear [ 11 , 12 , 13 ].
This complex scenario makes endometriosis a highly heterogeneous disease characterized by infertility and pain with a strong immuno-inflammatory component [ 14 , 15 ]. Indeed, a dysregulated immune response seems to be one of the most critical aetiologic factors related to increased production of pro-inflammatory cytokines, autoantibodies, growth factors, oxidative stress, increased number and activation of macrophages, which are also responsible for adhesion, anatomy changes, scars, and fibrosis [ 16 ].
Until now, endometriosis has been underdiagnosed because symptoms are often confused with many other medical and surgical conditions [ 17 ]. In addition, the severity of the symptomatology is not related to the extension of the disease [ 18 ].
On the other hand, no imaging method can accurately detect pelvic endometriosis with the accuracy needed to replace surgery [ 2 ]. Therefore, a thorough understanding of the pathophysiology, implementing non-invasive tools, and discovering new biomarkers are essential for developing novel diagnostic approaches for this debilitating condition [ 19 , 20 ].
In recent years, the introduction of omics sciences employing modern technologies has improved the knowledge of several diseases, and it continues to gain importance both in clinical and research settings [ 21 , 22 ]. Genomics and proteomics provided a more complete understanding of some biological events related to gene and protein expression. At the same time, metabolomics is defined as the “quantitative measurement of the dynamic multi-parametric metabolic response of living systems to pathophysiologic stimuli or genetic modification” [ 23 ] and offers the possibility to quantify and identify a broad range of endogenous and exogenous low molecular weight metabolites (lipids, amino acids, sugars, biogenic amines, etc.), providing a comprehensive view of metabolic changes during various physiological and pathological processes [ 24 , 25 ]. Small metabolites play a fundamental role in biological systems, representing potential candidates to understand disease phenotypes and pathological mechanisms and impacting the improvement of diagnosis of disease, establishment of therapeutic approach, response to the therapy, and classification of the patients [ 26 ].
Several analytical techniques can acquire the metabolic profile and simultaneously measure many metabolites or potential biomarkers in a single experiment. Among these, mass spectrometry (MS) and nuclear magnetic resonance (NMR) are the most widely used techniques [ 27 , 28 , 29 , 30 ]. NMR is a rapid and nondestructive analytical method that requires simple sample preparation, and spectra can be used to identify several metabolites. However, only the most abundant metabolites are detected because it is relatively insensitive, and millimolar to high micromolar concentrations are often required [ 31 ]. Conversely, MS approaches offer combined sensitivity and selectivity platforms for metabolomics analysis [ 32 ]. Selectivity among similar metabolites can be ensured by various MS/MS techniques and by combination with separation modules, such as gas chromatography (GC) or liquid chromatography (LC). Several pieces of evidence suggested that a specific metabolomics profile of endometriosis patients resulted from the analysis of the samples with MS techniques, which allowed the measurement of several classes of metabolites, including lipids, in both serum and peritoneal fluid [ 33 , 34 , 35 ].
The metabolic approach is a valuable and promising way to discover new biomarkers for endometriosis and reveal metabolic pathways associated with disease pathogenesis [ 20 , 22 , 36 , 37 , 38 , 39 , 40 , 41 ]. Because of the unique strength of this field, our study focused on metabolic changes in endometriosis patients (E) using gas chromatography–mass spectrometry (GC–MS) to improve understanding of metabolic pathways involved in endometrial pathogenesis and to provide new functional biomarkers. GC–MS is a high-throughput analytical platform that has been used for untargeted studies of altered metabolism in a variety of applications [ 28 , 42 ]. GC–MS-based metabolomics allow the identification and quantification of small-molecule metabolites (<650 Da), drugs, and toxins, using chemical derivatization to make these compounds volatile for the analysis. Several database annotations using large spectral libraries identify over 200 compounds from human body fluids [ 43 ]. As a second point of our study, we aimed to compare the patterns observed in the endometriosis group with those detected in a cohort of asymptomatic women without endometriosis (NE).
Results
This study aimed to investigate serum metabolic changes associated with endometriosis. Therefore, thirty-two serum samples from women were analyzed with GC–MS: twenty-two patients were affected by symptomatic endometriosis, and ten were unaffected.
Twenty-two metabolites, including organic, amino, fatty, and sugars, were identified in the samples.
As a first step, PCA non-supervised analysis was performed, and through the application of the T 2 -Hotelling test, eight outliers (six E and two NE) were identified and excluded from the subsequent analysis. The PLS-DA analysis was then performed to explore the presence of any differences between the NE and pathological E samples ( Figure 1 A).
The resulting model showed good statistical parameters: R 2 X = 0.45, R 2 Y = 0.731, Q 2 = 0.453, and p = 0.003. Then, the permutation test was used to validate the model: R 2 intercept = 0.42; Q 2 intercept = −0.24 ( Figure 1 B).
The VIP list obtained from the PLS-DA model showed that seven metabolites were the most responsible for separating the two groups, having a VIP value greater than 1. Metabolites were threonic acid, 3-hydroxybutyric acid, proline, alanine, valine, lactate, and phenylalanine.
These seven metabolites underwent a Mann–Whitney U test to find significant changes. Only five metabolites showed a p -value < 0.05 ( Figure 2 A) and were used to build the ROC curves to test their sensitivity and specificity as diagnostic biomarkers ( Figure 2 B).
Moreover, we performed the combined ROC-curve with the metabolites, which significantly changed their concentration, and we found that the best ROC-curve was obtained combined with 3-OH butyrate, threonic acid, and alanine (AUC = 0.91 CI = 0.8–1, p -value = 0.002).
To better understand which pathways were involved in the disease, all metabolites with a VIP value greater than 1 were used for the pathway analysis ( Figure 3 ). As shown in the figure, the altered pathways were the synthesis and degradation of ketone bodies, the biosynthesis of phenylalanine, tyrosine, tryptophan, and phenylalanine metabolism.
Discussion
Endometriosis is a well-characterized disease and represents an example of a complex disease in which different factors are involved, such as inflammation, cell invasion and angiogenesis, genetic mutations, steroid hormone metabolism, and oxidative stress [ 12 , 14 , 51 , 52 ]. It is well accepted that identifying precise serum biomarkers for endometriosis could avoid the need for laparoscopy, making diagnosing the disease less invasive and more accessible. A non-invasive early diagnosis could have an impact on the management and medical or surgical treatment of the disease [ 53 , 54 ]. The biomarker most consistently studied in endometriosis has been CA125 [ 55 , 56 , 57 , 58 ]. Mol and colleagues found that CA125 may be more beneficial in diagnosing advanced stages (III–IV) than stages I and II of the disease [ 59 ]. There is increasing evidence that global metabolic profiles are useful for identifying diagnostic biomarkers in easily accessible biofluids [ 22 , 60 , 61 , 62 ]. Metabolomics has emerged as an effective tool for monitoring disease progression [ 63 ] and distinguishing between diseased and non-diseased states [ 64 ].
The present prospective observational translational study is an attempt to identify differently expressed metabolites in serum samples from endometriosis patients using GC–MS-based metabolomics, trying to clarify the pathogenesis of this disease and find new biomarkers for a non-invasive diagnosis. We start from a previous experience in the Division of Gynecology and Obstetrics, University of Cagliari, Italy, to study the identification and quantification of a wide range of metabolites in serum samples from women affected and not affected by endometriosis using 1 H-NMR [ 22 ].
Partially in line with our previous study, the results showed an increase in the concentrations of proline, β-hydroxybutyric acid (a common feature with our past study), and threonic acid, and a decrease in valine and alanine, as well as an alteration of pathways such as biosynthesis and synthesis of phenylalanine, tyrosine, and tryptophan, and degradation of ketone bodies.
The development of fibrosis in endometrial lesions is a complex phenomenon whose fundamental mechanisms are not yet fully understood. Fibrosis is consistently present in all forms of endometriosis, leading to classic symptoms related to the disease, such as pain and infertility [ 65 ]. Proline is an important and essential component of collagen and has been shown to fuel protein production, which is needed for cell proliferation [ 66 ]. Some previous studies suggest that collagen concentration is significantly higher in chronic endometriosis-associated fibrosis and may play an important role in patients with endometriosis [ 67 ]. In addition, Plaka et al. [ 68 ] have shown an increase in serum proline levels in these patients. It could be related to increased proline biosynthesis and its increased use in collagen synthesis as support. Consistent with this evidence, our results showed an increase in serum proline levels in patients affected by endometriosis, similar to Dutta et al. [ 38 ], who also found a positive correlation between serum proline levels with early (minimal/mild) disease. However, as demonstrated by this latter mentioned study, proline’s trend in endometriosis patients could be controversial, especially considering the different degrees of severity of the disease. Moreover, proline level was decreased in follicular fluids of patients affected by endometriosis compared to control subjects [ 69 ].
Oxidative stress, defined as an imbalance between reactive oxygen species (ROS) and antioxidants, can cause general inflammation reactions in the peritoneal cavity and may be involved in the pathophysiology of endometriosis [ 51 , 70 ]. In our study, we found a significant increase in threonic acid in patients with endometriosis: it is a primary product of the oxidative breakdown of ascorbic acid [ 71 ]. The elevated serum levels of this sugar acid shown in our study, as well as the altered balance in the concentration of ketone bodies (which are by-products of the glutathione oxidation pathway), further support the fundamental role that oxidative stress plays in the pathogenesis of this disease [ 72 , 73 ]. Furthermore, changes in the ketone body concentration balance may be indirect signs of oxidative stress, as ketone is a by-product of the glutathione oxidation pathway, and the strong relationship between endometriosis and oxidative stress has already been demonstrated [ 74 , 75 ].
Ketone bodies are also involved in energy production. Glucose metabolism fuels 95% or more ATP production, but in case of an energetic challenge, the sources of alternative fuels are ketones. In fact, they provide acetyl coenzyme A (acetyl-CoA), which enters the tricarboxylic acid cycle (TCA) to generate ATP [ 76 ]. Bioenergetic metabolism is essential for cells of the immune system to perform its specific functions [ 77 ]. A situation of chronic immune activation, such as endometriosis, can exceed the physiological bioenergetic metabolism consuming a considerable amount of energy (up to 2000 kJ/day and more) [ 78 ].
Furthermore, elevated ketone body β-hydroxybutyric acid levels impaired bone marrow cell proliferation, lymphocyte proliferation in vitro, and in vitro chemotactic differentials of leukocytes in animals [ 79 ]. This mechanism could also characterize women with endometriosis. In addition, serum concentrations of β-hydroxybutyric acid were also high in patients with ovarian cancer [ 80 ]. Some endometriotic tissue behaviors may be similar to those of cancerous cells, so that some metabolic changes may be shared between the two diseases [ 81 ].
We also found a decrease in alanine in the serum of patients affected by endometriosis compared to those not affected. This result is in line with Jana et al.’s study [ 82 ], where an increased utilization of this amino acid as a significant gluconeogenic precursor was suggested to meet the high glucose demand and consumption by fast-growing endometrial cells outside of uterine tissue. In the presence of energetic imbalance, alanine blood levels decrease more than any other amino acid, mainly due to the fuel for gluconeogenesis [ 83 ]. This phenomenon is driven by the glucose-alanine cycle, where pyruvate is indirectly transformed into alanine through transamination in skeletal muscle. Interestingly, some studies suggested that ketone body infusion lowers L-alanine blood levels more than any other amino acid [ 84 , 85 ]. Despite the fact that we found a decreased level of alanine, other studies found an increased level of it in patients affected by endometriosis compared to control subjects, as well as valine, which we found decreased in endometriosis patients compared to NE [ 86 ]. It is a branched-chain amino acid (BCAA) and, together with leucine and isoleucine, can provide carbon precursors useful for synthesizing glucose and other molecules to fuel TCA metabolism and oxidative phosphorylation and provide energy for cells [ 87 ]. Furthermore, BCAAs can supply nitrogen nucleotide and protein synthesis [ 88 ]. However, the interactions between BCAAs and systemic metabolism are not entirely understood. Still, interesting evidence suggested a role in maintaining immune function and homeostasis [ 89 ] with the regulation of several biological processes, such as autophagy and cell growth and proliferation [ 90 ]. Furthermore, in vitro studies reported that immune cells use BCAA as substrates, showing an increase in the uptake of these amino acids [ 91 ]. As endometriosis is strongly associated with an inflammatory state [ 14 ], and the immune system is believed to play a central role in its etiology, pathophysiology [ 92 ], and associated morbidities, it is reasonable to expect a decrease in BCAA valine. Our findings are consistent with Santonastaso et al.’s metabolomic study, where a valine level reduction was found in the follicular fluid of endometriosis patients [ 69 ], although other studies showed a reverse trend [ 38 , 93 ].
Considering the results obtained through the pathway analysis, the biosynthesis of phenylalanine, tyrosine, and tryptophan, together with the metabolism of ketone bodies previously discussed, emerged as the most altered.
The biosynthesis of phenylalanine, tyrosine, and tryptophan is strongly linked to an inflammatory environment [ 86 ]. The relationship between the metabolic pathways of phenylalanine and inflammation has been reported in the literature, and it has been observed that during body inflammation, the expression of phenylalanine levels increases. Inflammatory cytokines induce muscle breakdown with the consequent release of phenylalanine as a substrate for gluconeogenesis to supply the increased metabolic demand during inflammation processes [ 94 ]. Therefore, a higher phenylalanine concentration could be associated with the catabolic state.
Tryptophan is involved in two major metabolic pathways: kynurenine and serotonin (5-HT) formation [ 95 ]. The kynurenine metabolic pathway (KP) represents about > 90% of peripheral tryptophan metabolism in mammals and plays a crucial role in many biological processes, including immunoregulation [ 96 , 97 , 98 ]. Dysfunction of immune responses is suggested to promote endometriosis [ 99 ] which, together with the presence of inflammation, is correlated with alteration in the kynurenine metabolic pathway with significant changes in the concentrations of its metabolic products [ 100 , 101 ]. The activity of indolamine 2,3-dioxygenase (IDO), the enzyme that promotes tryptophan degradation along the KP, is strongly induced by pro-inflammatory stimuli [ 102 ]. It has been proposed that, during inflammation, tryptophan metabolism occurs mainly in immune cells, while in homeostatic/healthy conditions, most TRP degradation occurs in the liver [ 103 ]. This can be explained by the increased demand for energy metabolism, redox balance in immune cells, and the requirement of KP metabolites for immune regulatory functions [ 104 , 105 ]. Tryptophan metabolites, 4,6-dihydroxyquinoline, and 5-hydroxy-L-tryptophan were markedly altered in the urine of endometriosis model rats [ 106 ].
Conclusions
The results of the present study demonstrated that the application of the GC–MS metabolomics strategy to discover new potential biomarkers has been useful in the endometriosis diagnostic process and in clarifying the disease’s pathophysiology. Indeed, a specific pathological metabolic fingerprint, in line with our previous work, was established. Furthermore, increased concentrations of proline, β-hydroxybutyric acid, and threonic acid and a decrease in valine and alanine were found in the endometriosis patients compared to the non-endometriosis group. This paper describes a translational study where the common clinical practice “marries” the basic research with innovative analytical and statistical techniques.
However, despite the innovative nature of this study, it shows some limitations. First, the number of enrolled patients is relatively small, and further studies are needed with a larger cohort of patients, possibly with the inclusion of a group of healthy control subjects to better define the pathological metabolic profile of patients affected by endometriosis. Second, the analysis of the mean age of the two classes of patients showed a p -value < 0.05. Despite the fact that we didn’t find a correlation between age and the metabolomics profile, this could be considered a confounding factor. Furthermore, it could be desirable for the combination of different analytical platforms to amplify the panel of the metabolites, which could be candidates as biomarkers to improve the clinical management of the patients affected by endometriosis. In this light, NMR spectroscopy, LC–MS, and GC–MS can give complementary information.
Finally, considering that fibrosis represents a hallmark in the endometriosis scenario, the lack of the evaluation of its level in the patients and the correlation with the metabolic changes represent a weak point of this study, as well as the lack of the measure of the enzymatic activity relative to the pathways that we found altered, which could confirm their involvement in the pathological processes.
In conclusion, literature data about the metabolic changes in patients affected by endometriosis are still conflicting. For this reason, further studies are needed to better clarify this interesting topic.
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.