Skip to main content
Download PDF

Urinary phytoestrogens and the risk of uterine leiomyomata in US women

BMC Women's Health volume 23, Article number: 261 (2023) Cite this article

Abstract

Background

Uterine leiomyomata (UL) is a common gynecological disease in women. Studied on the relationship between single metabolites of urinary phytoestrogens and UL, especially for the combined effects of mixed metabolites on UL still are insufficient.

Methods

In this cross-sectional study, we included 1,579 participants from the National Health and Nutrition Examination Survey. Urinary phytoestrogens were assessed by measuring urinary excretion of daidzein, genistein, equol, O-desmethylangolensin, enterodiol and enterolactone. The outcome was defined as UL. Weighted logistic regression was used to analyze the association between single metabolites of urinary phytoestrogens and UL. Notably, we adopted the weighted quantile sum (WQS) regression, Bayesian kernel machine regression (BKMR), and quantile g-computation (qgcomp) models, to investigate the combined effects of six mixed metabolites on UL.

Results

The prevalence of UL was approximately 12.92%. After adjusting age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers, total energy, daidzein, genistein, O-desmethylangolensin, enterodiol, and enterolactone, the association of equol with UL was significant [Odds ratio (OR) = 1.92, 95% confidence interval (CI): 1.09–3.38]. In the WQS model, mixed metabolites of urinary phytoestrogen had a positive association with UL (OR = 1.68, 95%CI: 1.12–2.51), with the highest weighted chemical of equol. In the gpcomp model, equol had the largest positive weight, followed by genistein and enterodiol. In the BKMR model, equol and enterodiol have positive correlation on UL risk, while enterolactone has negative correlation.

Conclusion

Our results implied a positive association between the mixed metabolites of urinary phytoestrogen and UL. This study provides evidence that urinary phytoestrogen-metabolite mixture was closely related to the risk of female UL.

Peer Review reports

Background

Uterine leiomyomata (UL) are the most common solid tumors in women [1, 2]. It is estimated that up to 80% of women will develop UL during their lifetime [1, 3], with 25–30% of them experiencing significant symptoms, including chronic pelvic pain, dysmenorrhea, abnormal vaginal discharge, and abnormal menstruation [3, 4]. UL continues to pose a serious disease burden for women.

Although the underlying pathology of UL is not particularly clear, it has been suggested to be an estrogen-dependent tumor [5]. Phytoestrogens are a group of plant compounds that are similar in chemical structure to mammalian estrogens, and they can be absorbed from food, circulate in the blood, and are excreted in the urine [6,7,8]. Previous studies have reported the effects of phytoestrogens on UL [9, 10]. For example, a case–control study included 328 eligible subjects from the Diagnostic Unit of the University Hospital of the West Indies, found that there was no association of urine daidzein, genistein, equol, enterolactone, total phytoestrogens and uterine fibroid (diagnosed by abdominal and/or vaginal ultrasonography) using binary logistic regression analysis [9]. A cross-sectional study contains 1,204 participants performed by Zhang Y, et al., implied that equol was significantly associated with the risk of UL after adjusting for age, race, pregnant status, ovary removed status, use of female hormones, body mass index (BMI), menopausal status and urinary creatinine levels [10]. There are still contradictions regarding the effects of phytoestrogens on uterine fibroids. Importantly, these studies on the association of phytoestrogens with UL have focused on the effects of single chemicals [9, 10]. Generally speaking, humans are often exposed to many chemicals simultaneously, and the cumulative effect of multiple chemicals is of concern [11]. Nevertheless, little is known about the mixed effects of multiple chemicals in phytoestrogens on UL.

Herein, this study aimed at investigating the relationship between single metabolites of urinary phytoestrogens and UL in US women, and the combined effects of mixed metabolites on UL risk.

Methods

Population selection

In this cross-sectional study, all data were drawn from the National Health and Nutrition Examination Survey (NHANES) database. The NHANES is a cross-sectional survey conducted by the National Center for Health Statistics (NCHS) of the Centers for Disease Control and Prevention using a multilayer probability sampling design, which aim to assess the health and nutritional status of adults and children in the United States [12]. The NHANES survey combines interviews and physical examinations [13]. The requirement of ethical approval and informed consent of the subjects for this was waived by the Institutional Review Board of The First Affiliated Hospital of Soochow University, because the data was accessed from NHANES (a publicly available database). All methods were carried out in accordance with relevant guidelines and regulations.

In this study, we used data from four cycles of the NHANES database (NHANES 1999–2000, 2001–2002, 2003–2004, 2005–2006). For participants in the NHANES database, only women aged 20–54 were asked diagnostic questions about UL (n = 6,508). Participants who met one of the following criteria were excluded: (1) Women without measurement of urinary phytoestrogen concentrations; (2) Women without assessment of UL; (3) Women with missing information of covariates related to UL. Ultimately, 1,579 participants were included in this study (Fig. 1).

Fig. 1
figure 1

Flowchart of population selection. NHANES = National Health and Nutrition Examination Survey; UL = uterine leiomyomata

Full size image

Assessment of urinary phytoestrogen

Urinary phytoestrogens were assessed by measuring urinary excretion of isoflavones (including daidzein, genistein, equol, and O-desmethylangolensin) and enterolignans (including enterodiol and enterolactone) [14]. The collection of urine specimens was carried out in the Mobile Examination Centers, and stored at -20 °C until analyzed [14]. The analyses of urinary excretion were accomplished by using the high-performance liquid chromatography (HPLC)-tandem mass spectrometric (MS) detection in the survey 1999–2004 and HPLC-atmospheric pressure photoionization- MS in the survey 2005–2006 [15]. For 1,579 participants of this study, 1 participant were below the lower limit of detection (LOD) for daidzein (0.40 ng/mL), 9 participants were below the lower LOD for genistein (0.20 ng/mL), 2 participants were below the lower LOD for equol (0.06 ng/mL), 29 participants were below the lower LOD for O-desmethylangolensin (0.20 ng/mL), 0 participants were below the lower LOD for enterodiol (0.04 ng/mL) and 0 participants were below the lower LOD for enterolactone (0.10 ng/mL) [16]. In the case of results below the LOD, the value of this variable is the LOD divided by the square root of two (https://wwwn.cdc.gov/Nchs/Nhanes/1999-2000/PHPYPA.htm#URXDAZ). The concentration of daidzein, genistein, equol, O-desmethylangolensin, enterodiol, and enterolactone in urinary phytoestrogens was corrected by creatinine in this study. Geometric mean and tertiles of each phytoestrogen metabolite (ug/g creatinine) were presented in Supplemental Table 1.

Assessment of uterine leiomyomata

The outcome was considered as UL. Participants in the NHANES database were classified as patients with UL when they answered “Yes” to the question “Has a doctor or other health professional ever told you that you had uterine fibroids?”.

Potential covariates

We extracted some characteristics of participants from the NHANES database, including age (years), race/ethnicity (non-Hispanic White/ non-Hispanic Black/ others), marital status (married/ never married/ others), education level [high school and below/ high school grad/ general educational development (GED) or equivalent/ some college or associate of arts (AA) degree/college graduate or above], poverty-to-income ratio (PIR, < 1.0/ ≥ 1.0), smoking status (yes/no), drinking status (yes/no), BMI (kg/m2), waist circumference (cm), cotinine (ng/mL), age at menarche (years), menopausal status (yes/no), ovary removed status (yes/no), hysterectomy (yes/no), use of female hormones (yes/no), hormones/hormone modifiers, pregnancy status (yes/no), number of gravidities, fiber (gm) and total energy (kcal). PIR was classified as in the NHANES database ≥ 1.0 (meaning household income was above the poverty line) and < 1.0 (meaning household income is at or below the poverty line). Smoking status and drinking status in the NHANES database was based on participants’ self‐report. BMI was calculated as weight (kg) divided by height squared (m2). Cotinine was assessed measured in serum using isotope dilution-high performance liquid chromatography/atmospheric pressure chemical ionization tandem mass spectrometry. Similarly, when the result is below the LOD, the value of cotinine is the LOD divided by the square root of two. Information on age at menarche, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers, pregnancy status and number of gravidities was obtained from the reproductive health questionnaire. Use of female hormones was judged by self-report " Have you/Has SP ever used female hormones such as estrogen and progesterone?" and drug code 97–101 in the NHANES database. Hormones/hormone modifiers was defined according to drug codes [97–98, 97–103, 97–288, 97–295, 97–377, 97–411, 97–413, 97–414, 97–416, 97–417, 97–418, 97–420, 97–422, 97–423, 97–426, 97–495].

Statistical analysis

Given the nature of the complex sampling of the NHANES database, we used a weighted analysis: weight variables for the urinary metabolites measurement (WTSB2YR and WTSPH2YR) and study design variables (SDMVPSU and SDMVSTRA). The measurement data were tested for normality using Kolmogorov–Smirnov, and normally distributed measurement data were described as mean (standard error) [Mean (SE)] and compared between two groups using independent samples t-test; non-normal data were described as median and quartiles [M (Q1, Q3)] and compared between groups using Mann–Whitney U rank sum test. Categorical data were described as number of cases and composition ratio N (%) and compared between groups using chi-square test and rank data using rank sum test. In the present study, we adopted chain equation multiple interpolation method based on random forest for some missing data of the variables. The miceforest package in python is used for interpolation processing (https://pypi.org/project/miceforest/). A sensitivity analysis was performed on the data before and after interpolation (Supplemental Table 2). SAS (version 9.4), Python (version 3.9) and R (version 4.0) software were used for statistical analyses. P < 0.05 was considered as statistically significant difference.

First, we performed weighted univariate logistic regression to screen covariates. Then, weighted logistic regression was used to analyze the association between single metabolites of urinary phytoestrogens and UL. Odds ratio (OR) and 95% confidence interval (CI) were calculated in the study. Last, we adopted three statistical models: weighted quantile sum (WQS) regression, Bayesian kernel machine regression (BKMR), and quantile g-computation (qgcomp) models, to investigate the effects of six mixed metabolites on UL.

Weighted quantile sum (WQS) regression

WQS regression was used to investigate the effects of six mixed metabolites on UL and identify the predominant metabolite. The study sample was randomly divided into training dataset (30%, n = 474) and validation dataset (70%, n = 1,105). Exposure to each metabolite in the training dataset was first divided into tertiles. The tertiles were then added together to generate an overall tertiles score for each metabolite. An empirical weight for each metabolite in the mixture was estimated using the bootstrapping method [17]. The WQS score is a combination of six mixed metabolites, representing the whole-body burden of six urinary phytoestrogens [10]. The weight of each metabolite in the WQS score indicates the contribution of each metabolite to the overall result [18]. Metabolites with an estimated weight greater than 0.333 (1/3) were considered to be significant contributors to the WQS score. Using 10,000 bootstrap samples from the training dataset (30%), we calculated the weights for WQS scores. Using the validation dataset (70%), we assess the statistical significance of WQS scores [19]. In addition, WQS regression requires that all exposure-outcome associations be focused in the same direction. Therefore, we estimated the positive and negative effects of the six metabolites on UL separately. R package gWQS was adopted to perform the analysis.

Quantile g-computation (qgcomp) model

gqcomp is a parameterized and generalized linear model based on application of g-computation, aimed to assess the effect of increasing all exposures in the mixture by one quatile simultaneously [20]. In this study, the gqcomp.noboot function was applied to estimate exposure effects, which divides six mixed metabolites into tertiles, assigns a positive or negative weight to each metabolite. If a metabolite has multiple effects in different directions, a positive or negative weight is interpreted as the proportion of exposure effects that have a negative (or positive) effect on UL, with a total weight of up to 2. The relationship of each metabolite endpoint and the mixed metabolites was assessed separately, and the finding models were used to estimate the scaled effect sizes, variable-specific coefficients, and overall model fit p-values. Metabolites with an estimated weight greater than 0.05 were considered to be significant contributors to the gqcomp scores. R package qgcomp was adopted to perform the analysis.

Bayesian kernel machine regression (BKMR)

BKMR is a supervised approach, which could identify nonlinear and nonadditive associations of exposure-outcome [21]. In this study, the BKMR model with 10,000 iterations was adopted. Genistein, equol and enterodiol were divided into two groups according to their positive correlation with UL, while daidzein, O-desmethylangolensin, and enterolactone were divided into one group according to their negative correlation with UL. The combined effect was calculated by comparing mixed metabolites at or above the 60th percentile with the 50th percentile. Group posterior inclusion probability (GroupPIP) and Conditional posterior inclusion probability (CondPIP) represent the probability of each group and metabolite in each group included in the model, representing their contribution to the overall effect. R package bkmr was adopted to perform the analysis.

Results

Population characteristics

Table 1 presents the general characteristics of 1,579 eligible participants. The average age was 37.81 years. Approximately 69.00% of participants reported a history of drinking, and 32.14% of participants indicated that they were menopausal. In addition, all participants were divided into UL group (n = 204) and non-UL group (n = 1,375). Age, race/ethnicity, marital status, drinking status, BMI, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers, number of gravidities and total energy were significantly different between UL group and non-UL group (P < 0.05).

Table 1 The general characteristics of included participants
Full size table

Correlation between single metabolites of urinary phytoestrogens and UL

As shown in Supplemental Table 3, the result of univariate logistic regression indicated that age, race/ethnicity, marital status, drinking status, BMI, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy might be covariates for this current study. The weighted logistic regression was used to assess the individual effect of each metabolite on UL (Table 2). After adjusting for age, race/ethnicity, marital status, drinking status, BMI, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy, equol in the tertile 3 showed significant association with UL (Model 1: OR = 1.92, 95%CI: 1.07–3.43, P = 0.029). After further adjusting for age, race/ethnicity, marital status, drinking status, BMI, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers, total energy, daidzein, genistein, O-desmethylangolensin, enterodiol, and enterolactone, the association of equol in the tertile 3 with UL remained significant (Model 2: OR = 1.92, 95%CI: 1.09–3.38, P = 0.024; Fig. 2).

Table 2 The individual effect of each metabolite on UL by using weighted logistic regression
Full size table
Fig. 2
figure 2

The association between single metabolite of urinary phytoestrogens and uterine leiomyomata in women in the multivariable logistic regression model. Other metabolites of urinary phytoestrogen were further adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy

Full size image

WQS, qgcomp and BKMR models to assess the combined association between six metabolites and UL

The WQS model was employed to estimate the combined effect of six metabolites of urinary phytoestrogen on UL. In the adjusted model (Table 3), mixed metabolites of urinary phytoestrogen had a positive association with UL (P = 0.011), and a tertile increase in the WQS index was related to a 68% increased risk of UL (95%CI: 1.12–2.51). We also calculated the estimated chemical weights of for each WQS index (Fig. 3). The highest weighted chemical in the WQS model was equol, followed by enterodiol and genistein.

Table 3 WQS model to estimate association between six metabolites and UL
Full size table
Fig. 3
figure 3

WQS model regression index weights for uterine leiomyomata. Model was adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy

Full size image

Similar to the WQS model, a tertile increase in the gpcomp index was associated with risk of UL in the adjusted model (Table 4, OR = 1.51, 95%CI: 1.05–2.18, P = 0.027). Figure 4 shows the estimated weight of each metabolite on the UL risk. Equol had the largest positive weight, followed by genistein and enterodiol, respectively.

Table 4 Qgcomp model to assess the combined association between six metabolites and UL
Full size table
Fig. 4
figure 4

gqcomp model regression index weights of the mixture on uterine leiomyomata risk. Model was adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy

Full size image

Supplemental Table 4 summarizes the GroupPIP and CondPIP derived from the BKMR model for six metabolites. The GroupPIP of two group (genistein, equol and enterodiol; 0.34) was higher than one group (daidzein, O-desmethylangolensin, and enterolactone; 0.04). Enterodiol (CondPIP = 0.89) contributed most to the model for the UL risk. Figure 5 indicates the overall associations between six metabolites and UL risk. Although the high concentrations of all metabolites were not statistically different compared to their 50th percentile, the overall effect on UL of the mixture of exposures at the 60th and above quantiles showed an upward trend. As all other metabolites were at their median levels, equol and enterodiol have positive correlation on UL risk, while enterolactone has negative correlation (Supplemental Fig. 1). In addition, we also found that there may be an interaction between enterodiol and enterolactone on UL risk (Supplemental Fig. 2).

Fig. 5
figure 5

Combined effects of six metabolites of urinary phytoestrogens on uterine leiomyomata risk. Model was adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy

Full size image

Discussion

In this study including 1,579 US women, we assessed the relationship of urinary phytoestrogens and UL risk by using a number of statistical models. Overall, the weighted multivariate logistic regression indicated a correlation between equol and UL risk. By the WQS and gpcomp models, we observed a positive association between mixed metabolites of urinary phytoestrogen and UL risk. WQS model further identified that equol made the most contribution in the association between metabolite mixture of urinary phytoestrogen and UL risk. In the BKMR model, there was no significant association between overall mixed metabolites and UL appeared, but there was a trend towards an increase. Additionally, equol and enterodiol also showed a positive correlation with UL risk in gpcomp and BKMR models.

Previous studies have focused on the relationship between individual chemicals and health outcomes, but in fact, humans are often exposed to mixtures of multiple pollutants/chemicals [19, 22]. In recent years, several novel statistical methods have been developed to assess the impact of exposure to chemical mixtures on health outcomes, including WQS regression [17,18,19], gpcomp [20] and BKMR [21]. A review assessed the relationship between exposure to mixtures of per- and polyfluoroalkyl substances and adverse health outcomes, and highlighted the importance of WQS and BKMR for assessment of the effects of exposure to mixtures [23]. In addition, a cross-sectional study performed in US population found a positive association between combined exposures to mercury, arsenic, cadmium and lead measured in urine and higher estimated glomerular filtration rate using WQS regression [24], and they also indicated that there might be influence for exposure to multiple metals on kidney function. In the study of Zhang Y, et al., they reported that mixed exposure of ten commonly exposed endocrine-disrupting chemicals had a significant positive association with UL in WQS and BKMR models, the weight distribution showed the highest weights for mercury (weight = 0.35) and equol (weight = 0.29) [10]. However, to our knowledge, the association between the mixed metabolites of urinary phytoestrogen and UL has not been studied so far.

Unlike previous study [5, 10, 23], this study considered the mixed effect of six metabolites of urinary phytoestrogen (daidzein, genistein, equol, O-desmethylangolensin, enterodiol, and enterolactone) on UL risk by three approaches (WQS regression, qgcomp, and BKMR). These results also indicated that mixed metabolites of urinary phytoestrogen were positively linked to the UL risk, with the greatest effect being from equol. Equol was related to an increased risk of UL. Our results are also consistent with previous study [10]. Equol, a metabolite of soy isoflavone daidzein, has estrogenic and antioxidant activity [25]. Several studies have showed that equol has a beneficial impact on metabolic diseases [26, 27]. But, estrogen-dependent diseases such as UL, are likely to be exacerbated by the estrogenic effects of equol. As described in an animal study, equol may trigger uterine tissue hyperplasia by increasing luminal epithelial cell height and myometrial and stromal thickness, which further lead to UL [28]. Our results agree with a previous study that estradiol could stimulate growth of UL, and was considered to be associated with increased risk of UL [29]. Although we found a combined effects of mixed metabolites on UL risk, the molecular mechanism related to the relationship of phytoestrogen and UL remains unclear. Further exploration is needed regarding the potential mechanisms in the association.

The main strength of this study was the use of WQS regression, qgcomp, and BKMR, which allowed us to assess the mixed metabolites of urinary phytoestrogen and UL risk. Some limitations for this study should be considered. First of all, because of the design of this cross-sectional study, there was a limitation in the causal relationship between urinary phytoestrogens and UL. Second, some possible confounders were lacking in this NHANES database, such as family history of UL. We did not adjust for history of hysterectomy because they may be a consequence of the outcome [30]. Third, for participants in the NHANES database, a single spot urine sample was only collected for metabolites analysis. The concentrations of metabolites of phytoestrogens may vary over time. Fourth, we excluded 4,587 women who were not measurement of urinary phytoestrogen concentrations. Urinary phytoestrogens were tested in 1/3 of the participants aged 6 years and older in the NHANES database. However, this study considered the weights in the analysis, so the bias was relatively small. Prospective studies with large sample size are warranted to further analyze the relationship of urinary phytoestrogens and UL, and the related mechanisms.

Conclusion

In summary, our results implied an association of equol and UL. Importantly, WQS regression, qgcomp, and BKMR models was adopted to analyze the combined effects of mixed metabolites on UL risk. A positive association between the mixed metabolites of urinary phytoestrogen and UL was also identified, with the greatest contribution from equol. This study provides evidence that urinary phytoestrogen-metabolite mixture was closely related to the risk of female UL and further research is needed to explore the detailed mechanism.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the NHANES database, https://wwwn.cdc.gov/nchs/nhanes/.

Abbreviations

UL:

Uterine leiomyomata

NHANES:

National Health and Nutrition Examination Survey

NCHS:

National Center for Health Statistics

HPLC:

High-performance liquid chromatography

MS:

Mass spectrometric

BMI:

Body mass index

SE:

Standard error

OR:

Odds ratio

CI:

Confidence interval

WQS:

Weighted quantile sum

References

  1. Baird DD, Dunson DB, Hill MC, Cousins D, Schectman JM. High cumulative incidence of uterine leiomyoma in black and white women: ultrasound evidence. Am J Obstet Gynecol. 2003;188:100–7.

    Article  PubMed  Google Scholar 

  2. Brasky TM, Bethea TN, Wesselink AK, Wegienka GR, Baird DD, Wise LA. Dietary fat intake and risk of uterine leiomyomata: a prospective ultrasound study. Am J Epidemiol. 2020;189:1538–46.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Lewis TD, Malik M, Britten J, San Pablo AM, Catherino WH. A comprehensive review of the pharmacologic management of uterine leiomyoma. Biomed Res Int. 2018;2018:2414609.

    Article  PubMed  PubMed Central  Google Scholar 

  4. McWilliams MM, Chennathukuzhi VM. recent advances in uterine fibroid etiology. Semin Reprod Med. 2017;35:181–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Alset D, Pokudina IO, Butenko EV, Shkurat TP. The effect of estrogen-related genetic variants on the development of uterine leiomyoma: meta-analysis. Reprod Sci. 2022;29:1921–9.

    Article  CAS  PubMed  Google Scholar 

  6. Atkinson C, Lampe JW, Scholes D, Chen C, Wähälä K, Schwartz SM. Lignan and isoflavone excretion in relation to uterine fibroids: a case-control study of young to middle-aged women in the United States. Am J Clin Nutr. 2006;84:587–93.

    Article  CAS  PubMed  Google Scholar 

  7. Liu T, Li N, Yan YQ, Liu Y, Xiong K, Liu Y, et al. Recent advances in the anti-aging effects of phytoestrogens on collagen, water content, and oxidative stress. Phytother Res. 2020;34:435–47.

    Article  PubMed  Google Scholar 

  8. Jackson MD, McFarlane-Anderson ND, Simon GA, Bennett FI, Walker SP. Urinary phytoestrogens and risk of prostate cancer in Jamaican men. Cancer Causes Control. 2010;21:2249–57.

    Article  PubMed  Google Scholar 

  9. Simon GA, Fletcher HM, Golden K, McFarlane-Anderson ND. Urinary isoflavone and lignan phytoestrogen levels and risk of uterine fibroid in Jamaican women. Maturitas. 2015;82:170–5.

    Article  CAS  PubMed  Google Scholar 

  10. Zhang Y, Lu Y, Ma H, Xu Q, Wu X. Combined Exposure to Multiple Endocrine Disruptors and Uterine Leiomyomata and Endometriosis in US Women. Front Endocrinol (Lausanne). 2021;12:726876.

    Article  PubMed  Google Scholar 

  11. Chiu YH, Bellavia A, James-Todd T, Correia KF, Valeri L, Messerlian C, et al. Evaluating effects of prenatal exposure to phthalate mixtures on birth weight: a comparison of three statistical approaches. Environ Int. 2018;113:231–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Karia PS, Joshu CE, Visvanathan K. Association of oophorectomy and fat and lean body mass: evidence from a population-based sample of U.S Women. Cancer Epidemiol Biomarkers Prev. 2021;30:1424–32.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Karanth SD, Washington C, Cheng TD, Zhou D, Leeuwenburgh C, Braithwaite D, et al. inflammation in relation to sarcopenia and sarcopenic obesity among older adults living with chronic comorbidities: results from the National Health and Nutrition Examination Survey 1999–2006. Nutrients. 2021;13:3957.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Reger MK, Zollinger TW, Liu Z, Jones J, Zhang J. Association between urinary phytoestrogens and C-reactive protein in the continuous national health and nutrition examination survey. J Am Coll Nutr. 2017;36:434–41.

    Article  CAS  PubMed  Google Scholar 

  15. Xiong G, Huang C, Zou Y, Tao Z, Zou J, Huang J. Associations of urinary phytoestrogen concentrations with nonalcoholic fatty liver disease among adults. J Healthc Eng. 2022;2022:4912961.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Martínez Steele E, Monteiro CA. Association between dietary share of ultra-processed foods and urinary concentrations of phytoestrogens in the US. Nutrients. 2017;9(3):209.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Nguyen HD, Kim MS. Effects of heavy metals on cardiovascular diseases in pre and post-menopausal women: from big data to molecular mechanism involved. Environ Sci Pollut Res Int. 2022;29:77635–55.

    Article  CAS  PubMed  Google Scholar 

  18. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of weighted quantile sum regression for highly correlated data in a risk analysis setting. J Agric Biol Environ Stat. 2015;20:100–20.

    Article  PubMed  Google Scholar 

  19. Duc HN, Oh H, Kim MS. The effect of mixture of heavy metals on obesity in individuals ≥50 years of age. Biol Trace Elem Res. 2022;200:3554–71.

    Article  CAS  PubMed  Google Scholar 

  20. Keil AP, Buckley JP, O’Brien KM, Ferguson KK, Zhao S, White AJ. A quantile-based g-computation approach to addressing the effects of exposure mixtures. Environ Health Perspect. 2020;128:47004.

    Article  PubMed  Google Scholar 

  21. Bobb JF, Claus Henn B, Valeri L, Coull BA. Statistical software for analyzing the health effects of multiple concurrent exposures via Bayesian kernel machine regression. Environ Health. 2018;17:67.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Yu L, Liu W, Wang X, Ye Z, Tan Q, Qiu W, Nie X, Li M, Wang B, Chen W. A review of practical statistical methods used in epidemiological studies to estimate the health effects of multi-pollutant mixture. Environ Pollut. 2022;306:119356.

    Article  CAS  PubMed  Google Scholar 

  23. Rosato I, ZareJeddi M, Ledda C, Gallo E, Fletcher T, Pitter G, Batzella E, Canova C. How to investigate human health effects related to exposure to mixtures of per- and polyfluoroalkyl substances: a systematic review of statistical methods. Environ Res. 2022;205:112565.

    Article  CAS  PubMed  Google Scholar 

  24. Sanders AP, Mazzella MJ, Malin AJ, Hair GM, Busgang SA, Saland JM, et al. Combined exposure to lead, cadmium, mercury, and arsenic and kidney health in adolescents age 12–19 in NHANES 2009–2014. Environ Int. 2019;131:104993.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Mayo B, Vázquez L, Flórez AB. Equol: a bacterial metabolite from the daidzein isoflavone and its presumed beneficial health effects. Nutrients. 2019;11:2231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Liu J, Mi S, Du L, Li X, Li P, Jia K, et al. The associations between plasma phytoestrogens concentration and metabolic syndrome risks in Chinese population. PLoS One. 2018;13:e0194639.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Takahashi A, Kokubun M, Anzai Y, Kogre A, Ogata T, Imaizumi H, et al. Association between equol production and metabolic syndrome in Japanese women in their 50s–60s. Menopause. 2022;29:1196–9.

    Article  PubMed  Google Scholar 

  28. Brown NM, Lindley SL, Witte DP, Setchell KD. Impact of perinatal exposure to equol enantiomers on reproductive development in rodents. Reprod Toxicol. 2011;32:33–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Nowak RA. Fibroids: pathophysiology and current medical treatment. Baillieres Best Pract Res Clin Obstet Gynaecol. 1999;13:223–38.

    Article  CAS  PubMed  Google Scholar 

  30. Marsh EE, Al-Hendy A, Kappus D, Galitsky A, Stewart EA, Kerolous M. Burden, prevalence, and treatment of uterine fibroids: a survey of U.S. women. J Womens Health (Larchmt). 2018;27:1359–67.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

  1. Department of Obstetrics and Gynecology, The First Affiliated Hospital of Soochow University, No.899 Pinghai Road, Suzhou, 215006, P. R. China

    Fang Yang & Youguo Chen

Authors
  1. Fang Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Youguo Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

FY and YC designed the study. FY wrote the manuscript. FY and YC collected, analyzed and interpreted the data. YC critically reviewed, edited and approved the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Youguo Chen.

Ethics declarations

Ethics approval and consent to participate

The requirement of ethical approval and informed consent of the subjects for this was waived by the Institutional Review Board of The First Affiliated Hospital of Soochow University, because the data was accessed from NHANES (a publicly available database). All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplemental Table 1.

Distribution of urinary phytoestrogen levels. Supplemental Table 2. Sensitivity analysis of data before and after interpolation. Supplemental Table 3. The Selection of covariates by univariate logistic regression. Supplemental Table 4. GroupPIP and CondPIP of six metabolites.

Additional file 2: Supplemental Fig. 1.

Univariate exposure–response functionbetween metabolite exposure and UL with fixing all the other metabolites at their median level. Model was adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy.

Additional file 3: Supplemental Fig. 2.

Bivariate exposure–response function for metabolites in UL, with exposure 1 metabolite at its 10%, 50%, and 90% levels and other metabolites fixed at their median levels. Model was adjusted for age, race/ethnicity, marital status, drinking status, body mass index, waist circumference, menopausal status, ovary removed status, use of female hormones, hormones/hormone modifiers and total energy.

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

Yang, F., Chen, Y. Urinary phytoestrogens and the risk of uterine leiomyomata in US women. BMC Women's Health 23, 261 (2023). https://doi.org/10.1186/s12905-023-02381-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12905-023-02381-5

Keywords

BMC Women's Health

ISSN: 1472-6874

Contact us