Intro
Endometriosis (EM) affects approximately 10% of women of reproductive age ( 1 ) and is defined by the presence of endometrial glands or stoma outside the uterine cavity. It is a non-malignant disease nonetheless associated with dysmenorrhea, dyspareunia, pelvic pain and infertility due to the presence of ectopic tissue and inflammation ( 2 – 4 ). Endometriotic lesions can be superficial peritoneal, ovarian or deeply infiltrating ( 5 , 6 ). The etiology of endometriosis is complex ( 7 ) and multifactorial thus several theories have been proposed to explain the clinical manifestation of endometriosis; Sampson’s theory of retrograde menstruation, coelomic metaplasia theory, Mullerian rests theory, stem cell theory, impaired immune system theory, and others ( 8 ). Endometriosis is considered a chronic inflammatory disease with altered peritoneal environment in patients with endometriosis. The ectopic lesions recruit immune cells which leads to production of pro-inflammatory molecules and cytokines and also promote angiogenesis and innervation and thus contribute to survival of these lesions ( 9 ).
EM is an estrogen-dependent disease with molecular hallmarks of genetic predisposition, altered hormonal milieu (estrogen dependence and progesterone resistance) and inflammation ( 10 , 11 ). Changes in steroid biotransformation pathways have been reported leading to an increased local production of estrogens in endometriosis lesions ( 12 ). EM lesions respond to hormonal signals such as estradiol (E 2 ) that regulate uterine tissue growth and triggers inflammation, and that are linked to pain symptoms ( 13 ). Excessive inflammation also leads to changes in sex steroid receptors (ERα and ERβ) expression and enhanced estrogen biosynthesis in endometriotic lesions, involving aromatase, sulfatase and other pathways ( 14 – 21 ). Treatment of pelvic pain in EM thus includes the use of nonsteroidal anti-inflammatory drugs, oral contraceptives and progestins.
Estrogens comprise a vast array of hydroxylated (OH) and methoxylated (MeO) catechol estrogen (CE) metabolites with diverse biological activities. The synthesis of CE metabolites from E 2 and estrone (E 1 ) involves various metabolic routes, namely the 2-hydroxylation (2OH), 4-hydroxylation (4OH) and 16-hydroxylation (16OH) pathways and the action of the catechol-O-methyltransferase (COMT) to form 2- and 4- MeOCEs ( 22 , 23 ). Besides acting as ligands of ERα and ERβ, CEs also present non-estrogenic properties ( 22 , 24 ). MeOCEs have antiangiogenic and antiproliferative actions whereas 4-OHCEs have procarcinogenic properties ( 22 , 25 – 27 ). Both the 2OH and 4OH CE derivatives generally have reduced estrogenic effects ( 24 , 28 ) as opposed to the 16OH pathway that retains most of its estrogenic properties, with a preferential action on the ERβ ( 24 ). In addition, E 2 was previously found to be associated with pain due to its effects on nerves and inflammation ( 13 ). Some CEs present estrogenic activities resembling E 2 and they may be prone to cause pain.
Given the suspected biological roles of catechol estrogens, and their associations with several hormone-sensitive diseases including endometrial cancer ( 29 , 30 ), we hypothesized that circulating levels of E 2 and/or its biologically active metabolites were associated with an altered risk of EM (primary objective) and severity of pain symptoms (secondary objective).
Results
Characteristics of EM and controls are despicted in
Table 1
. A total of 16 estrogen derivatives were quantified by MS in the serums of 341 women. Most estrogens and their oxidative metabolites (11 out of 16) were above LLOQ, except for 4OH-E 2 , 17epi-E 3 , 2MeO-E 2 , 4MeO-E 1 and 4MeO-E 2 detected in less than 12% of the cohort. These five estrogens were thus excluded in subsequent statistical analyses (
Supplementary Table 2
). Levels of estrogens are displayed in
Table 2
with E 1 , E 2 , 2OH-E 1 , E 3 , 16αOH-E 1 and 16keto-E 2 displaying the highest levels. In addition, 2OH-3MeOE 1 levels were higher in cases compared to controls (by 123%; P = 0.02) (
Table 2
). The analysis of hormone levels according to anatomic sites of EM (
Table 3
), showed higher levels of 2OH-E 1 (by 18%; P = 0.009), 2MeO-E 1 (38%; P = 0.002), sum of 2OH (7%; P = 0.013), sum of MeO (23%; P = 0.036), ratio of 2OH/sum E 1 /E 2 (24%; P = 0.023), and ratio of 2OH-E 1/ 16αOH-E 1 (101%; P = 0.033) in cases diagnosed with ovarian EM. No evidence of an association was observed in relation to surgical stage (
Supplementary Table 3
).
Steroid levels for EM cases (n=209) and for controls (n=115) that included subjects with benign pathologies and healthy controls.
Steroids are reported in pg/mL. The sum of hydroxy derivatives (OHs) includes 2OH, 4OH and 16OH. The sum of all catechol estrogens (CEs) includes OHs and MeOs. P -values were obtained based on the one-way analysis of covariance F test with correction when necessary. Continuous means of hormone levels were log transformed and adjusted for age and BMI for statistical analysis. Significant P values <0.05 are highlighted with a star (*) and text in bold. The median values and the 10-90% intervals for E 2 were 123.00 (31.6 -417.0 pg/mL) for women in the proliferative phase and 169.00 (13.20-333.00 pg/mL) for women in the secretory phase.
Steroid levels according to anatomic sites of disease in 209 endometriosis cases.
The sum of hydroxy derivatives (OHs) includes 2OH, 4OH and 16OH. The sum of all catechol estrogens (CEs) includes OHs and MeOs. Data represent bivariate analysis corrected with Tukey and adjusted for age and BMI. P -values were obtained with the F test using log transformed data. Only trends ( P <0.1) and significant associations ( P <0.05; in bold) are displayed.
The first objective aimed to establish the potential association of estrogens with the risk of developing EM. Women with higher circulating levels of specific catechol estrogens were shown to be more predisposed to EM risk in multivariable analysis adjusted for confounders including age, BMI, tobacco status, contraception, hormonal therapy and endometrial phase, which differ in EM cases compared to controls (
Table 1
). More specifically, women with higher levels of 2OH-3MeO-E 1 had an adjusted OR of 1.91 (95%CI 1.09-3.34; P =0.025). This finding was also observed when restricted to healthy subjects and controls with benign pathologies (OR =2.61 (95%CI 0.84–8.09); P =0.097 and 1.56 (95%CI 0.84–2.91); P =0.164) but did not reached significance. A lower risk of EM was observed in association with elevated 16OH derivatives (with OR values of 0.22 for 16epi-E 3 (95%CI 0.06–0.83); P =0.025) and 0.22 for 16keto-E 2 (95%CI 0.06–0.84); P =0.027). Since 2OH-3MeO-E 1 was the main metabolite associated with the risk of EM, we evaluated whether the correlation of this metabolite with the other estrogens was different between controls and cases. We observed that 2OH-3MeO-E 1 was weakly but significantly correlated with 2-OH derivatives in EM cases at 0.20 ( P <0.05) but not in controls (
Table 4
).
Spearman correlation coefficients among endogenous hormone levels in endometriosis cases (EM), benign conditions (BC) and healthy (H).
Correlation values were similar with all controls (when BC and H were analyzed together). The P -values of <0.05 are identified with a (*) and highlighted in bold. Correlations were tested using the Spearman Rho statistic test.
A secondary objective aimed to explore the relationship between hormone levels and symptoms of pain in EM cases (
Table 5
). Higher levels of 2OH-3MeO-E 1 were associated with the risk of pain in the abdominal, pelvic and back regions (OR =1.86 (95%CI 1.06-3.27) P =0.032). Higher levels of 16αOH-E 1 were inversely associated with the risk of pain in the abdominal, pelvic and back regions (OR =0.55 (95%CI 0.31-0.97); P =0.038). More frequent menstrual pain was associated with elevated E 1 /E 2 (OR = 1.90 (95%CI 1.03-3.51); P =0.041) whereas more severe menstrual pain was linked to a higher metabolic ratio of 4OH/sum of E 1 /E 2 (OR =1.95; (95%CI 1.05-3.61); P =0.035) (
Table 5
). Higher levels of 16keto-E 2 were associated with the risk of more severe dyspareunia experienced in the last three months (OR = 2.40 (95%CI 1.02-5.62); P =0.045. Higher levels of 4OH-E 1 were associated with a reduced risk dysuria or dyschezia experienced in the last three months (OR =0.32 (95%CI 0.12-0.89); P =0.028) (
Supplementary Table 4
).
Significant associations between pain and steroid levels in EM cases (n=209).
The sum of hydroxy derivatives (OHs) includes 2OH, 4OH and 16OH. The sum of all catechol estrogens (CEs) includes OHs and MeOs. Odds ratios (OR) and their P-values were obtained using a logistic regression model adjusted for age and BMI. The number of cases with data on clinical outcomes are identified next to each outcome. There was no evidence of an association for the “score of dysmenorrhea” and the “dyspareunia (frequency)” outcomes. Significant P-values (<0.05) are shown in bold.
Discussion
EM is a complex estrogen-sensitive condition characterized by a chronic inflammation process for which the potential role of estrogen metabolites remains to be fully investigated. We report that higher levels of 2OH-3MeO-E 1 were associated with an increased risk of EM, with an approximately two-fold higher median level observed in circulation of EM cases compared to controls. An enrichment of the 2OH metabolic pathway, with significantly higher levels of 2OH-E 1 , 2MeO-E 1 , sum of MeO and ratio of 2OH-E 1/ 16αOH-E 1 , was also observed in ovarian EM cases compared to cases affected with lesions at other anatomical sites. A perturbation of estrogen metabolism (2OH-3MeO-E 1 and 16αOH-E 1 ) was further associated with pain symptoms.
Estrogens and their receptors play a key role in the pathophysiology of EM. Studies reported higher levels of systemic and locally synthesized estrogens in EM cases promoting the growth of lesions ( 38 ). This increase in estrogens was attributed to the secretion of estrogens by the ovaries as well as their autocrine and paracrine action, and an increased aromatase activity in EM lesions that supports local E 2 synthesis ( 15 , 17 , 39 – 41 ). Additional changes in estrogen synthesis, as well as their metabolic and receptor pathways, have also been reported in support of an enhanced local production and action in EM lesions, creating a hyperestrogenic environment that affects hormone receptor function ( 38 , 42 – 49 ). These changes may be reflected in circulation of EM cases with higher levels of E 2 and/or its metabolites. In our study, 2OH-3MeO-E 1 was associated with an increased risk of EM. This observation is consistent with elevated COMT expression in EM lesions ( 50 ) leading to the formation of 2OH-3MeO-E 1 from its precursor 2OH-E 1 ( 23 ), potentially contributing to higher systemic levels of this metabolite. In support, we showed that endometrial tissue can contribute to systemic estrogen levels in the context of endometrial cancer that significantly declined after surgery ( 29 ). In addition to significant higher circulating levels of 2OH-3MeO-E 1 metabolites in EM cases compared to controls, higher levels of 2MeO-E 1 and the sum of MeOs were observed but they did not reached significance. The biological properties of 2OH-3MeO-E 1 have been poorly studied. We further noted that this metabolite was less correlated in circulation with the other estrogen derivatives and particularly in control subjects, suggesting a dysregulation in the presence of EM lesions associated with its precursors such as 2OH-E 1 , with a higher correlation coefficient for this metabolite in EM cases at 0.20 ( P <0.05). In fact, 2OH-3MeO-E 1 was higher in cases compared to controls, supporting a potential EM origin. Consistent with our observation, a previous study evaluated a subset of estrogen metabolites in preoperative urine samples of 62 EM cases and 52 controls and found increased levels of the 2OH-3MeO-E 1 precursor 2OH-E 1 ( 51 ). Our findings that the 2OH pathway is significantly more elevated in ovarian EM cases is also consistent with a study that used proton nuclear magnetic resonance (H-NMR) spectroscopy to investigate potential non-invasive metabolomic markers in 31 infertile women with stage II and III EM cases and 15 healthy or control women ( 52 ). They found that levels of the antiangiogenic 2MeO-E 1 /E 2 metabolites were higher in EM cases compared to controls. The enrichment of the 2OH metabolic pathway in ovarian EM cases is supported by higher tissular levels of CYP1A1, involved in the conversion of E 1 and E 2 to 2OH-E 1 and 2OH-E 2, reported to be 4-fold higher in the ovarian EM group ( 50 , 53 ). Inversely, the 16OH pathway was inversely associated with EM, consistent with downregulation of the involved enzyme pathways (CYP3A) by inflammation ( 54 , 55 ).
An estrogenic environment may be associated with more severe pain symptoms ( 13 ). Hence, the association between estrogen levels and pain outcomes may not be related only to the effect of estrogens, as E 2 metabolites have been documented to present receptor-independent biological activities and may contribute to the maintenance of the inflammatory milieu ( 16 ). Previous reports revealed that elevated ERβ is associated with proliferation, inflammation and pain transmission ( 46 , 56 , 57 ), coherent with the positive correlation observed in this study between 16keto-E 2 and dyspareunia in EM cases. However, the negative association between 16αOH-E 1 and pelvic, abdominal and back pain suggests more complex relationships. A component of pain in EM was shown to be related to inflammatory damage of nerve fibers with neuroprotective roles for ERβ ( 58 – 61 ). Also, a dysregulation of both the ERα and ERβ expression pathways was observed in the ectopic endometrium in EM compared to normal endometrium in favor of a superior ERβ to ERα ratio ( 46 , 62 , 63 ). 16-hydroxylated derivatives are amongst E 2 metabolites known to bind the ERβ receptor ( 24 ), which may explain the observed association with pain. Other studies showed that the ERα was correlated with symptoms in deep infiltrating EM ( 64 ) and that it could favor hyperalgesia by altering calcium release ( 61 ). Since the 2OH metabolites are known to bind ERα ( 24 ), this could explain the association between pain outcomes and 2OH metabolites, such as 2OH-3MeO-E 1 . Additional studies are required to uncover the precise biological function of the 2OH-3MeO-E 1 metabolite.
This pilot study provides a comprehensive quantification of estrogens in the circulation of EM cases and controls based on a sensitive mass-spectrometry assay. It is comprised of a significant sample size, surgical and histologic confirmation of case and control status, adjustment for confounding factors and examination of pain symptoms. A limitation is the fact that the control group also included patients with gynecological conditions other than EM, which may influence the hormonal milieu ( 40 , 65 ). Although cases and controls differed in confounding factors such as menstrual phase, these variables were included in the multivariate model for EM risk. Exploratory analyses in relation to pain symptoms were adjusted for age and BMI. Additional studies could provide levels of progesterone and its metabolites, shown to be dysregulated in EM and recognized to counteract the effect of E 2 ( 48 , 66 ), whereas the endometriotic intratissue estrogen levels may not reflect the corresponding systemic levels. Due to the exploratory nature of the study, no correction for multiple testing was applied, but our initial findings warrant replication in other cohorts.
Conclusions
We conclude that the 2OH-3MeO-E 1 metabolite represents a potential adverse feature of EM and that the 2OH pathway is associated with the risk of ovarian endometriotic lesions. Data also suggest an association between the 2OH metabolic pathway and the risk of unfavorable pain outcomes.
Ethics Statement
The studies involving human participants were reviewed and approved by National Medical Ethics Committee in Slovenia (#0120-127/2016/6) and the Ethics committee of the CHUQc – Université Laval (#2012-993). The patients/participants provided their written informed consent to participate in this study.
Data Availability
The original contributions presented in the study are included in the article/
Supplementary Material
. Further inquiries can be directed to the corresponding authors.
Materials|Methods
The study design corresponded to a retrospective case-control study comprising cases and controls from the same type of population ( 31 ). Part of this cohort was described previously ( 32 – 34 ). Patients’ enrolment took place from March 2008 to June 2018 at the Departments of Obstetrics and Gynecology at the University Medical Centre Ljubljana, Slovenia. The study comprised patients who visited gynecologist with problems/symptoms that are indicative for laparoscopy surgery. The inclusion criteria were an indication for a diagnostic laparoscopy for symptoms suggestive of EM such as pain, infertility, ovarian cysts, other gynecological pathologies such as myomas and tubal sterilization. The exclusion criteria were pregnancy, age below 18 years, menopausal status, gynecological malignancies, cancelled surgery, previous hysterectomy, drug abuse and HIV infection ( 32 ). Of the 341 women, 17 participants were excluded to manage confounders: six because of missing data on mandatory age and/or BMI, three due to unknown case or control status, four with polycystic ovary syndrome (PCOS) that could impact hormone levels, one due to a prior unidentified menopausal state, and three with significant anomalies of the menstrual cycle of unknown etiology. The remaining cohort of 324 women underwent either diagnostic laparoscopy or laparoscopic tubal sterilization and were divided according to presence (n=209, cases) or absence (n=115, controls) of EM. Controls were further divided into two groups (patients with benign pathologies (n=79) and healthy controls (n=35). Patients with benign pathologies had symptoms suggestive for EM (infertility and/or pain) or other gynecological pathologies. Healthy patients underwent laparoscopic tubal sterilization and had no symptoms suggestive for EM (
Figure 1
). A post-hoc ANOVA power analysis test (power package and R Statistical Software v4.1.2; R Core Team 2021) estimated that a sample size per group of 36 was suffcient). The clinical characteristics presented in
Table 1
included age, body mass index (BMI), type of EM (ovarian, ovarian and peritoneal, peritoneal, and deep infiltrating), rASRM stage of disease ( 35 ), smoking status (current, former or never), use of hormonal therapy (last three months), use of oral contraception (last three months), and endometrial phase (secretory or proliferative). Patient-filled surveys using validated numeric rating scales documented the outcomes of “abdominal, pelvic and back pain”, “dysmenorrhea (frequency)”, “dysmenorrhea (intensity)”, “score of dysmenorrhea”, “dyspareunia (frequency)”, “dyspareunia (intensity)” and “dysuria or dyschezia (frequency)” ( 36 , 37 ). For EM cases, data for pain-related outcomes were available for 98.6% to 99.5% of participants, except for the “score of dysmenorrhea’’ (61.2%), “dysmenorrhea (intensity)” (35.4%), and “dyspareunia (intensity)” (59.3%) outcomes. In control cases, data was available for 86.7% to 97.4% of participants, except for “score of dysmenorrhea’’ (56.5%) and “dyspareunia (intensity)” (55.7%) outcomes. The patient’s characteristics related to pain symptoms are presented in
Supplementary Table 1
. Pain related outcomes were dichotomized for the statistical analysis. The dichotomization for “abdominal, pelvic and back pain” was “yes or no”. For “dysmenorrhea (frequency)”, “dyspareunia (frequency)” and “dysuria or dyschezia (frequency)”, the dichotomization was “infrequent (never, almost never or sometimes) or frequent (quite often or very often)”. For “dyspareunia (intensity)” and “dysmenorrhea (intensity)”, the dichotomization was “mild (no or slight pain) or moderate to severe (medium or strong pain)”. For the “score of dysmenorrhea”, the dichotomization was “scores of ≤5 or of >5”. All participants provided an informed consent prior to their enrolment. This study was conducted in accordance with the declaration of Helsinki. This study was approved by the National Medical Ethics Committee in Slovenia (#0120-127/2016/6) and the ethics committee of the CHUQc – Université Laval (#2012-993).
Flowchart of patient cohort.
Characteristics and clinical data of endometriosis cases (n=209) and controls (n=115).
Surgical stages were classified using the revised American Society for Reproductive Medicine score (rASRM) classification of endometriosis. 1
P -values <0.0001 obtained using Pearson’s chi square or Fisher’s exact test when appropriate. 2 Last 3 months.
Blood samples were collected two days prior to surgery as described ( 32 ), and following strict standard operating procedures for collection, processing and storage at -80˚C to preserve stability of metabolites such as steroids. Briefly, 4 ml of blood sample was collected by venipuncture from the median cubital vein using BD Vacutainer tubes (#369032; Becton Dickinson and Company, NJ, USA). The collected samples were incubated for no more than 1 h at room temperature and then centrifuged at 1400 × g for 10 min at room temperature. The seperated serum was collected, aliquoted, and stored at -80°C until analysis. Only samples that were frozen/thawed once were used for analysis. A specific set of 16 estrogen derivatives were quantified in 250 μL of serum using a liquid-chromatography tandem mass spectrometry assay (LC-MS/MS) as described ( 29 ). The lower limit of quantification (LLOQ) was 5 pg/mL. Sums including all analytes and metabolic ratios were calculated for the different metabolic pathways. Catechol estrogens at levels below LLOQ (even if detected above the limit of detection) were considered undetected.
Differences in estrogen hormone levels between cases and controls, anatomic sites and surgical stages were determined by bivariate analyses on means of log transformed continuous hormone levels. The relationship between hormone levels was assessed using Spearman’s rank-order correlation. Odd ratios (OR) were obtained using dichotomized hormone levels (independent variables) based on the median levels of controls as performed in previous studies ( 29 ) in a multivariate logistic regression model, adjusted for age, BMI, smoking status, oral contraception (last three months), hormonal therapy (last three months) and the endometrial phase (secretory or proliferative). Logistic models and Fisher’s scoring were used to determine the ORs for pain related outcomes in cases. P -values were obtained using Pearson’s chi square, Fisher’s exact test, the one-way analysis of covariance F-test corrected with Tukey when appropriate or the Spearman Rho statistic test in the appropriate contexts. Results were considered statistically significant when P < 0.05. Statistical analyses were performed by the statistician (DS) using the software SAS 9.4 by SAS Institute Inc. (Cary, NC, USA). Due to the exploratory nature of the study, no adjustment for multiple comparison was done.
Author Contributions
Study concept and design: TR, CG. Patient recruitment and clinical data: MP, AV, JO, TR. Conducted experiments and mass spectrometry: PC, VT. Statistical analyses: DS. Drafting of the manuscript: J-PE, CG. Critical revision of the manuscript for important intellectual content: All authors. Obtaining funding: TR, CG.
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.