Total Effective Xenoestrogen Burden in Serum Samples and Risk of Endometrial Cancer in the Spanish Screenwide Case-Control Study.

Compared with controls, cases were more likely to be obese, nulliparous, have a history of cancer and were less likely to have ever used hormonal contraceptives ( Table 1 ). History of cancer included breast cancer in 45.4% ( n = 5 / 11 ) of controls and 45.8% ( n = 13 / 24 ) of cases, whereas the rest of the past cancers included a variety of types, such as colon, ovary, and skin melanoma and nonmelanoma. Cases had similar concentrations of serum TEXB- α and TEXB- β levels, PCB-138, PCB-180, and HCB and higher p , p ′ -DDE and β -HCH concentrations than controls ( Table 3 ). Descriptive characteristics of controls according to serum TEXB- α and TEXB- β levels are shown in Table 2 . Significant differences were only observed in serum TEXB- α according to BMI among the controls. Serum levels of TEXB- α and TEXB- β were highly correlated among controls [Pearson correlation coefficient for continuous log-transformed variables: 0.84 (95% CI: 0.78, 0.88), and among cases: 0.69 (95% CI: 0.60, 0.77)]. None of the individual contaminants were positively correlated with serum TEXB levels among the controls ( Table 2 ). Serum levels of TEXB and individual organohalogenated compounds in Spanish participants diagnosed with endometrial cancer and controls frequency-matched to cases by age ( N = 306 ), Screenwide case–control study, 2017–2021. Note: LOD for TEXB- α and TEXB- β were set at 0.1  Eeq pM / mL , and LOD for individual organohalogenated compounds were set at 0.05  ng / mL . HCB, hexachlorobenzene; HCH, hexachlorocyclohexane; IQR, interquartile range; LOD, limit of detection; PCB, polychlorinated biphenyl ; p , p ′ -DDE, p , p ′ -dichlorodiphenyldichloroethylene; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Two-sample Wilcoxon rank-sum (Mann–Whitney) test comparing cases and controls. Median TEXB- α and TEXB- β levels did not significantly differ among cases ( 3.02 × 10 − 1 and 1.25  Eeq pM / mL , for TEXB- α and TEXB- β , respectively) and controls ( 4.17 × 10 − 1 and 1.28  Eeq pM / mL ; p = 0.653 and 0.933, respectively, Table 3 ). In adjusted models, the risk for endometrial cancer did not linearly increase with increasing category levels of both TEXB- α and TEXB- β ( p linear  trend = 0.566 and 0.679, respectively, Table 4 ). However, positive associations were observed among those in the second category of exposure compared with the lowest [ OR = 2.11 (95% CI: 1.13, 3.94) for TEXB- α , and OR = 3.32 (95% CI: 1.62, 6.81) for TEXB- β ], whereas no significant associations were observed among those in the third category of exposure [ OR = 1.22 (95% CI: 0.64, 2.31) for TEXB- α , and OR = 1.58 (95% CI: 0.75, 3.33) for TEXB- β ]. In mutually adjusted models for TEXB- α and TEXB- β levels, the association of TEXB- α with endometrial cancer risk was attenuated [ OR = 1.45 (95% CI: 0.61, 3.47) for the second category of exposure compared with the lowest], whereas estimates for TEXB- β were slightly lower but still statistically significant [ OR = 2.68 (95% CI: 1.03, 6.99); Table 4 ]. Generalized additive models revealed an inverted-U risk trend across serum TEXB- α and TEXB- β levels, with a sharp increase in risk for exposures < 0.5 and 1  Eeq pM / mL , respectively, and a downturn at higher levels ( Figure 1 ). Associations between TEXB and endometrial cancer risk among participants diagnosed with endometrial cancer and controls frequency-matched to cases by age ( N = 303 ), Screenwide case–control study, 2017–2021. Note: BMI, body mass index; CI, confidence interval; LOD, limit of detection; OR, odds ratio; Ref, reference; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Three participants were dropped from the model owing to missing data given that there were no missing values in the case group. Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y) and education level (high school or below, some college, college or above). Obtained using Wald test from logistic models, representing whether the estimate of the category differ significantly from the reference category. Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), and past history of cancer (no, yes). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), past history of cancer (no, yes), and the other fraction of TEXB ( < LOD , below median, above median). Participants were grouped into three groups of exposure, those below the LOD were classified as the lowest exposure category ( < LOD ; < 0.1 Eeq pM/mL), and those with detectable values were grouped considering the median of serum TEXB- α ( median = 5.72 × 10 − 1 Eeq pM/mL) and TEXB- β levels ( median = 1.79 Eeq pM/mL) based on their distributions among the controls. Obtained using Wald test from logistic models treating ordinal variables as continuous variables, representing whether there is a statistically significant linear trend in the estimates across the ordered categories. Generalized additive models for serum levels of TEXB of (A) α and (B) β fractions and endometrial cancer with 2 degrees of freedom among participants in the Screenwide case–control study ( n = 303 ), 2017–2021. Generalized additive models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), past history of cancer (no, yes), and the other fraction of TEXB ( < LOD , below median, above median). Median for TEXB- α = 5.72 × 10 − 1 and median for TEXB- β = 1.79 Eeq pM/mL. Smooth estimates are represented in solid lines and 95% confidence intervals in dashed lines. Numeric data can be found in Excel Tables S1(alpha) and S2 (beta). Note: BMI, body mass index; LOD, limit of detection; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Associations of TEXB- α by cancer histology in multivariate models revealed significant associations for endometrioid cancers between the second category of exposure compared with the lowest ( OR = 2.37 ; 95% CI: 1.19, 4.72), whereas no statistically significant associations were found for non-endometrioid cancers ( OR = 1.36 ; 95% CI: 0.50, 3.73). For TEXB- β , significant associations were observed in both cancer groups [ OR = 2.74 (95% CI: 1.23, 6.09), and OR = 4.44 (95% CI: 1.39, 14.24) for endometrioid and non-endometrioid, respectively, comparing the second category of exposure with the lowest; Table 5 ]. When serum TEXB- α and TEXB- β levels were mutually adjusted, associations for endometrioid cancers were no longer statistically significant [ OR = 1.68 (95% CI: 0.63, 4.50), and 2.07 (95% CI: 0.70, 6.14) for TEXB- α and TEXB- β , respectively], and the association between TEXB- β and non-endometrioid cancers remained similar [ OR = 4.41 (95% CI: 1.06, 18.34), comparing the second category of exposure with the lowest]. Associations for endometrial cancer and TEXB, by histological subtype ( N = 303 ), Screenwide case–control study, 2017–2021. Note: Three participants were dropped from the model owing to missing data given that there were no missing values in the case group. BMI, body mass index; CI, confidence interval; LOD, limit of detection; OR, odds ratio; Ref, reference; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), and education level (high school or below, some college, college or above). Obtained using Wald test from logistic models, representing whether the estimate of the category differ significantly from the reference category. Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), and past history of cancer (no, yes). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), past history of cancer (no, yes), and the other fraction of TEXB ( < LOD , below median, above median). Participants were grouped into three groups of exposure, those below the LOD were classified as the lowest exposure category ( < LOD ; < 0.1 Eeq pM/mL), and those with detectable values were grouped considering the median of serum TEXB- α ( median = 5.72 × 10 − 1 Eeq pM/mL) and TEXB- β levels ( median = 1.79 Eeq pM/mL) based on their distributions among the controls. Obtained using Wald test from logistic models treating ordinal variables as continuous variables, representing whether there is a statistically significant linear trend in the estimates across the ordered categories. Individual xenoestrogens contained in the α -fraction generally showed null associations with endometrial cancer, except for a few positive associations ( Table 6 ). In particular, the ORs for endometrial cancer comparing the third with the first tertile in multivariate models were 1.68 (95% CI: 0.94, 3.00) for p , p ′ -DDE, 0.95 (95% CI: 0.50, 1.80) for HCB, 1.11 (95% CI: 0.62, 2.01) for α -HCH , 2.01 (95% CI: 0.99, 4.09) for β -HCH , 1.25 (95% CI: 0.70, 2.25) for γ -HCH , 0.96 (95% CI: 0.50, 1.88) for PCB-138, and 1.03 (95% CI: 0.56, 1.9) for PCB-180. Models further adjusted for TEXB- β yielded similar but attenuated associations [1.62 (95% CI: 0.89, 2.94) for p , p ′ -DDE, 0.98 (95% CI: 0.51, 1.89) for HCB, 1.16 (95% CI: 0.63, 2.12) for α -HCH , 1.82 (95% CI: 0.88, 3.76) for β -HCH , 1.38 (95% CI: 0.76, 2.53) for γ -HCH , 0.81 (95% CI: 0.40, 1.61) for PCB-138, and 1.04 (95% CI: 0.55, 1.95) for PCB-180, comparing the third with the first tertile; Table 5 ]. Associations between individual organohalogenated compounds and endometrial cancer risk among participants in the Screenwide case–control study ( N = 303 ), 2017–2021. Note: Three participants were dropped from the model owing to missing data given that there were no missing values in the case group. BMI, body mass index; CI, confidence interval; HCB, hexachlorobenzene; HCH, hexachlorocyclohexane; LOD, limit of detection; OR, odds ratio; PCBs, polychlorinated biphenyls; p , p ′ -DDE, p , p ′ -dichlorodiphenyldichloroethylene; Ref, reference; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), and education level (high school or below, some college, college or above). Obtained using Wald test from logistic models, representing whether the estimate of the category differ significantly from the reference category. Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), and past history of cancer (no, yes). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), past history of cancer (no, yes), and fraction of TEXB- β [ < LOD ( < 0.1 Eeq pM/mL,) below median, above median; TEXB- β median = 1.79 Eeq pM/mL]. Obtained using Wald test from logistic models treating ordinal variables as continuous variables, representing whether there is a statistically significant linear trend in the estimates across the ordered categories. Stratified analyses by type of control and menopausal status yielded no significant interactions, although estimates for the β   fraction were higher among nongynecologic controls ( Figure 2 ). Results by BMI revealed lower estimates among those participants with obesity ( BMI ≥ 30   kg / m 2 ), although interactions were not statistically significant. Similarly, estimates were generally higher among hormonal contraceptive users than among nonusers, and a 6-fold increase in risk was observed for TEXB- β comparing the second category of exposure with the lowest, although the interactions were not statistically significant. Odds ratios (ORs) for endometrial cancer comparing categories (above and below the median compared with below the LOD) of TEXB of α and β fractions by subgroups, among participants in the Screenwide case–control study ( n = 303 ), 2017–2021. Numbers of participants in each category are provided for controls (Co) and cases (Ca). Logistic models adjusted for age ( < 60 , 60–69, ≥ 70 y), education (high school or below, some college, college or above), BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), past history of cancer (no, yes), and the other fraction of TEXB ( < LOD , below median, above median). Median for TEXB- α = 5.72 × 10 − 1 and median for TEXB- β = 1.79 Eeq pM/mL. Results for normal BMI are omitted owing to the limited number of participants in the case categories. All p -interactions > 0.05 . Note: BMI, body mass index; CI, confidence interval; HC, hormonal contraception; PMH, postmenopausal hormone use; LOD, limit of detection; Ref, reference; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens).

Cases were recruited within the Screenwide study, a Spanish case–control study (2017–2021), in which 182 consecutive incident endometrial cancer cases were enrolled, as well as 266 hospital controls (including 190 women with benign gynecological pathology and 76 women with nongynecological diseases) frequency-matched to cases by age. Benign gynecological pathology among controls included conditions such as endometriosis, fibroids, benign cysts, prolapse, and polyps, as well as regular gynecologic controls in healthy women. Inclusion criteria were having an intact uterus and, in addition for cases, having an incident diagnosis of endometrial cancer. Endometrial cancer cases were classified into the classic two broad subtypes, the so-called type I, which included the endometrioid histological subtype, and type II, which included the rest of histological subtypes. 37 Exclusion criteria for all participants included pregnancy, puerperium (8 wk), communication difficulties, intellectual disability, and treatment with chemotherapy or radiotherapy during the previous 6 months. Response rates among the eligible participants were 89.7% for cases, 80.4% for controls with benign gynecological pathology, and 76.8% for hospital controls without benign gynecological pathology. Further details have been published elsewhere. 38 Participants were identified and enrolled through outpatient visits and visits for anesthesia prior to surgery. Peripheral blood serum was drawn from participants upon enrollment before any cancer treatment, using SST II Advance BD Vacutainers, and was aliquoted and stored at − 80 ° C . Data on all covariates included in model adjustments were self-collected through face-to-face interviews conducted by trained personnel. Reference dates were age at diagnosis for cases and age at enrollment for controls. In cases, the age at diagnosis closely matched the age at enrollment, and the questions referred to typical exposures before the cancer diagnosis. Anthropometric questions included the usual weight 1 y prior to the diagnosis, which was used to compute body mass index (BMI). Lifetime occupational history was assessed for all jobs held for at least 1 y. Prior to the interview, all participants were informed about the purpose of the study and signed an informed consent form. Among those cases who agreed to donate blood ( N = 174 ), we did not include samples from 6 cases with < 2 mL of serum available, and 1 sample from a case enrolled after the shipment to the site of analyses. We further randomly discarded 11 samples owing to budget constraints, yielding 156 samples from cases for analysis. We included 150 controls (75 controls with benign gynecological pathology, and 75 women attending hospital for nongynecological diseases), who were frequency-matched to cases by 5-y age intervals. The study followed all the requirements established by the Ethics Committee for Clinical Research and was approved by the Ethics Committee for Clinical Research from the Bellvitge University Hospital (references: PR128/16 and PR348/19), as previously described in Peremiquel-Trillas et al. 38 Participation in the study was voluntary, and all eligible participants signed an informed consent form. The Screenwide study followed the national and international directives on ethics and data protection (Declaration of Helsinki and subsequent amendments; EU Reglament 2016/679) and the Spanish laws on data protection (Organic Law 3/2018; Law 14/2007 Biomedical Research). The study was registered in the National Register of Biobanks/Collections (C.0004389). The total effective xenoestrogen burden (TEXB), which assesses the combined estrogenic effect of chemical mixtures, was assessed in serum samples from cases and controls, collected when the participants signed the informed consent form—in cases before starting oncological treatment—using a previously standardized methodology described by Fernández et al. 35 and Pastor-Barriuso et al. 36 Briefly, 3 mL of methanol were added to each serum sample ( 3 mL ), and the mixture was extracted with 7.5 mL of hexane:ethyl ether (1:1 vol/vol), using a Bond Elut polychlorinated biphenyl (PCB) cartridge (Varian), previously prepared with 1.5 mL of hexane. The dried eluates were reconstituted with 200 μ L of hexane and injected in duplicate ( 100 + 100 μ L ) into a high-performance liquid chromatograph (HPLC). The semipreparative HPLC method was designed and validated to effectively separate persistent halogenated xenoestrogens from endogenous hormones. The method collects two HPLC fractions: the α -fraction , containing the more lipophilic compounds (PCBs, organochlorine pesticides and metabolites, and halogenated bisphenols, among others), and the β -fraction , containing the endogenous hormones together with more polar xenoestrogens (nonhalogenated bisphenols, phytoestrogens, polyphenols, and mycoestrogens) other than those eluted in the α -fraction . Of the first 11 mL collected (11 min) corresponding to the α -fraction , 1 mL was reserved for subsequent analysis of the presence and concentration of some persistent compounds. Organohalogenated compounds typically found within the α -fraction , including the organochlorine pesticides p , p ′ -dichlorodiphenyldichloroethylene ( p , p ′ -DDE), hexachlorobenzene (HCB), hexachlorocyclohexane (HCH), and three PCBs congeners (138, 153, and 180), were quantified by gas chromatography-tandem mass spectrometry (GC-MS/MS), using an Agilent 7890 GC with split-splitless inlet and 7693 ALS autosampler. 36 These environmental pollutants were selected, not only for their activity as xenoestrogenic endocrine disruptors through different pathways, but also for their associations with other hormone-dependent tumors such as breast cancer. 39 p -Chlorodibenzophenone was used as the internal standard and a limit of detection (LOD) for all of the chemicals was set at 0.05  ng / mL . The remaining 10 -mL duplicates of each fraction were combined, dried, resuspended in 1 mL of experimental steroid-free medium (phenol red-free medium supplemented with 1 mL of charcoal–dextran fetal bovine serum), and tested for estrogenic activity in the E-SCREEN bioassay. 40 The E-SCREEN bioassay is an in vitro assay based on the enhanced proliferation of human breast cancer cells (MCF-7) in the presence of estrogen-active substances. The MCF-7 BUS subtype, characterized as providing the highest proliferative response to estradiol, was used. The α and β fractions were assayed in triplicate, at three different dilutions (1:1, 1:5, and 1:10), along with negative (cells treated with experimental steroid-free mediums) and positive controls (cells treated with 100 pM of estradiol). The proliferative effect (PE) of α and β fractions were expressed as MCF-7 cell proliferation–fold over control, transformed into estradiol equivalent units (Eeq) by reading from a dose–response curve with estradiol; and expressed as the estradiol equivalent concentration in picomolar per milliliter of serum (Eeq pM/mL). 36 The LODs for TEXB- α and TEXB- β were set up at 0.1  Eeq pM / mL , which was the lowest concentration required to induce a significantly different PE from that observed in steroid-free control cells. For quality control, 10 serum samples were analyzed in triplicate through independent extraction, HPLC fractionation, and E-SCREEN bioassay. The interassay coefficients of variation for TEXB- α and TEXB- β were 15.9% and 14.4%, respectively. Researchers at the University of Granada and Biosanitary Institute of Granada (ibs.GRANADA) labs, responsible for the extraction, fractionation, and cell culture analysis (i.e., E-SCREEN) were blinded to the characteristics of the study population. Descriptive analyses were performed using medians and interquartile ranges (IQRs) per participant for continuous data and counts and percentages for categorical data. Chi-square tests were used for categorical variables, and Mann–Whitney tests for continuous variables. Pearson coefficients were employed to assess the correlation between continuous variables. A level equal to half of the first value above the LOD was imputed (single imputation 41 ) for 30.1% ( n = 92 ) and 18.6% ( n = 57 ) of women with TEXB- α and TEXB- β determinations below the LOD, respectively, as well as for the individual contaminants (47.4% for p , p ′ -DDE, 26.8% for HCB, 0.0% for α -HCH , 0.0% for β -HCH , 0.3% for γ -HCH , 22.9% for PCB-138, and 0.0% for PCB-180). Participants were grouped into three groups of exposure: Those below the LOD ( < 0.1 Eeq pM/mL) were classified as the lowest exposure category, and those with detectable values were grouped considering the median of serum TEXB- α and TEXB- β levels based on their distributions among the controls. We used unconditional logistic regression models adjusted for potential confounders to estimate odds ratios (ORs) and 95% confidence intervals (CIs). Estimated OR for endometrial cancer risk comparing tertiles of specific organohalogenated compounds [PCB-138, PCB-180, HCB, HCH (isomers α , β , and γ ), and p , p ′ -DDE] based on their control distributions, were also calculated. Basic models included adjustment for age at reference date ( < 60 , 60–69, ≥ 70 y) and education level (high school or below, some college, college or above). We further adjusted multivariable models for BMI ( < 25 , 25–29.9, ≥ 30   kg / m 2 ), use of hormonal contraceptives (ever, never), parity (0, 1–2, ≥ 3 ), and past history of other nonendometrial cancers (yes, no). These variables were selected based on their association with the outcome (endometrial cancer) or the exposure (TEXB) in univariate analyses, following the definition of a confounder. Variables associated with endometrial cancer in the bivariate analyses were BMI ( p < 0.001 ), use of hormonal contraceptives ( p = 0.002 ), parity ( p = 0.036 ), and past history of cancer ( p = 0.039 ; Table 1 ). BMI was also associated with TEXB- α ( p = 0.021 ; Table 2 ). Similarly, a third model was further mutually adjusted for serum TEXB- α and TEXB- β levels ( < LOD , below median, above median), that is, models evaluating TEXB- α levels were adjusted for TEXB- β levels and vice versa. Given that during premenopause, estrogen levels are higher, we performed sensitivity analyses restricted to postmenopausal women. Similarly, analyses restricted to never users of postmenopausal hormone therapy, and stratified analyses by hormonal contraception (ever, never), were performed to exclude the effect of exogenous hormone use. We conducted analyses exclusively on women with no prior history of cancer to investigate the potential influence of past cancer treatments on hormonal burden. Selection bias due to the inclusion of controls with benign gynecologic conditions was also inspected using stratified analyses by type of control (gynecologic, nongynecologic). Given that xenoestrogens may accumulate in adipose tissue, stratified analyses were additionally performed by BMI (normal, overweight, obese). We performed likelihood-ratio tests for the interaction terms between main exposure variables and the specified variables in logistic models. Results with fewer than five participants per category were omitted, given that they were based on a limited number of participants. In addition, generalized additive models were used to inspect the linearity of associations of log-transformed serum TEXB- α and TEXB- β levels with endometrial cancer risk. To test for linear trend, ordinal variables were treated as continuous variables. Data were missing for < 5 % of participants for all variables, except for use of hormone replacement therapy, where data for < 10 % were missing. In models, missing covariate values were treated as distinct categories. All tests were two-tailed with a significance level of 0.05. All analyses were conducted by using Stata (version 16.0; StataCorp). Descriptive characteristics [ n or n (%)] of Spanish participants diagnosed with endometrial cancer and controls frequency-matched to cases by age ( N = 306 ), Screenwide case–control study, 2017–2021. Note: BMI, body mass index; NA, not applicable; PMH, postmenopausal hormone use. Reference dates were age at diagnosis for cases and age at enrollment for controls. Chi-squared calculated without missing values comparing cases with controls. Descriptive characteristics [ n (%) or median (IQR range) by serum TEXB- α and TEXB- β levels among Spanish controls ( N = 150 ), in the Screenwide case–control study, 2017–2021. Note: BMI, body mass index; HCB, hexachlorobenzene; HCH, hexachlorocyclohexane; LOD, limit of detection; PCB, polychlorinated biphenyl; PMH, postmenopausal hormone use; p , p ′ -DDE, p , p ′ -dichlorodiphenyldichloroethylene; TEXB, total effective xenoestrogen burden; TEXB- α , TEXB of the alpha fraction (represents the combined estrogenic effect of mixtures of lipophilic organohalogenated xenoestrogens); TEXB- β , TEXB of the beta fraction (represents the combined estrogenic activity of endogenous hormones and more polar xenoestrogens). Participants were grouped into three groups of exposure, those below the LOD were classified as the lowest exposure category ( < LOD ; < 0.1 Eeq pM/mL), and those with detectable values were grouped considering the median of serum TEXB- α ( median = 5.72 × 10 − 1 Eeq pM/mL) and TEXB- β levels ( median = 1.79 Eeq pM/mL). Reference dates were age at diagnosis for cases and age at enrollment for controls. Chi-squared calculated without missing values. Nonparametric equality-of-medians test.

This study has identified for the first time a positive association between serum TEXB levels and endometrial cancer risk. In particular, an inverted-U risk trend across serum TEXB- α levels was observed, and women in the second category of exposure had a 45% increase in risk, but that was not seen among those in the highest exposure categories compared with the lowest exposure category. TEXB- β values were also positively associated with endometrial cancer, with an almost 3-fold increase in risk among women in the second category of exposure compared with the lowest. Positive associations were somewhat expected, given that circulating endogenous estrogens were included in the β   fraction, and that exposure to circulating endogenous estrogens is an established risk factor for endometrial cancer. 13 , 16 , 18 Associations by endometrial type revealed that estimates for TEXB- β were higher among non-endometrioid cancers than for endometrioid cancers, whereas typically non-endometrioid cancers have been considered less driven by estrogens. 8 , 19 However, our results were limited by small sample sizes for type II cancers. Other stratified analyses yielded estimates generally > 1 and all p -interactions were not statistically significant. An inverted-U risk trend for endometrial cancer in relation to combined exposure to serum lipophilic xenoestrogens ( TEXB- α ) is biologically plausible, given that both endogenous hormones and some EDCs have shown to exert nonclassical dose–response curves, with lower doses exerting more potent effects than higher doses. 42 Nonmonotonic dose–response curves are common in endocrinology. Thus, in vivo studies in rats have shown that estrogen’s impact on mammary gland development varies with dosage, with low-to-moderate doses promoting terminal end bud formation and ductal elongation, whereas higher doses hinder these processes. 43 For example, in the US Food and Drug Administration Clarity study evaluating Sprague–Dawley rats, perinatal exposure to bisphenol A (BPA) was found to induce breast cancer at the lowest tested dose, but not at the highest doses. 44 Similarly, related studies revealed nonmonotonic effects of BPA on the developing rat mammary gland. 45 Other nonmonotonic dose–response curves have also been observed in epidemiological studies for different EDCs and several outcomes, including diabetes 46 – 48 and BMI, 49 , 50 among others. 51 These nonmonotonic dose–response curves suggest intricate interactions in which estrogen can act as both agonist and antagonist on multiple targets. 43 , 44 , 52 The nonmonotonic dose response observed in our study could be attributed to competition to specific receptors between endogenous hormones and xenoestrogens present in serum samples. 51 In situations where endogenous hormones do not fully saturate these receptors, xenoestrogens may bind to the available ones, thereby enhancing the overall cellular response. Under these circumstances, at low-to-moderate TEXB levels, xenoestrogens mimic the actions of endogenous estrogen by promoting endometrial proliferation, as suggested using mathematical models. 51 , 53 Tumor cell proliferation effects in endometrial cancer are mainly mediated through ER α , whereas the role of ER β remains inconclusive. 54 However, when TEXB levels are high, xenoestrogens may surpass natural ligands in ER α binding competition, potentially leading to an attenuation of the overall biological response owing to their lower estrogenic activity and/or potency. This competition could lead to partial antagonistic behavior, reducing the stimulatory impact of xenoestrogen in endometrial tissues and, in turn, potentially diminishing the risk of endometrial cancer. 51 , 53 Positive estimates were observed for β -HCH and p , p ′ -DDE, whereas for the rest of identified chemical compounds, no clear associations with endometrial cancer were found. The individual chemicals, analyzed in this study, were also not associated with TEXB- α , consistent with what has been previously described in the literature. 35 , 36 , 55 Several factors, including possible interactions established between the chemical compounds involving additive, synergistic or antagonistic mechanisms, and/or unmeasured substances, could explain the absence of agreement. Humans are exposed to both lipophilic and more polar xenoestrogens EDCs, such as PCB, brominated flame retardants, alkylphenolic compounds, polycyclic aromatic hydrocarbons compounds, parabens, phthalates, BPA, pesticides, and certain metals, among others. 21 , 42 , 56 Epidemiology literature evaluating xenoestrogens in endometrial cancer is generally scarce. Two human studies evaluated exposure to individual organochlorines and PCBs using serum exposure biomarkers, yielding null associations in accordance with results on individual contaminants. 27 , 28 Other epidemiological studies have evaluated organochlorines and PCBs using dietary variables, 23 as well as adipose tissue exposure biomarkers among reduced samples sizes, 25 which also yielded null results. Dietary assessments have been also used to assess personal exposure to cadmium, with positive results in a meta-analyses of two studies, 29 but not in following large prospective studies. 22 , 24 Urinary concentrations of BPA, parabens, and phthalate metabolites have been recently evaluated in a study among 139 case–control sets, yielding null results for all single chemicals evaluated, except for mono- n -butyl phthalate (MnBP) and dibutyl phthalate (DBP) excretion. These had statistically significant associations for the second tertiles of exposure, but not the third tertiles, 26 also suggesting an inverted-U risk shape pattern consistent with our results. The null results observed for parabens by Sarink et al. 26 contrasted with the positive associations observed by Dogan et al., 30 although Dogan et al. assessed individual chemical exposure in endometrial and myometrial tissue samples in a smaller study sample size (33 cases and 49 controls). Urinary biomarkers have been also used to evaluate exposure to few alkylphenols in 49 endometrial cancer cases, showing positive associations for chemical exposures above the median. 31 Interestingly, most of the above individual chemical compounds tested elutes in the β -HPLC fraction. In addition to the heterogeneity of exposure designs and assessments, most epidemiological studies assessing EDCs and endometrial cancer have focused on the concentration of individual chemicals or a limited number of chemicals. Historically, this one-chemical-at-a-time approach to exposure assessment has led to insufficient knowledge of the human health effects caused by exposure to mixtures. There was a need to implement improved tools, such as in vitro bioassays, that allow the analysis of both exposure and biological activity resulting from the mixture of chemicals present in human biological samples. 57 TEXB of the α -fraction represents the combined estrogenic effect of organohalogenated lipophilic xenoestrogens mixtures, whereas the TEXB of the β fraction represents the combined estrogenic activity of endogenous hormones plus more polar xenoestrogens. Thus, the TEXB bioassay represents an effective way to explore the combined impact of these compounds and to overcome the unpredictability of xenoestrogen interactions derived from possible additive, synergistic, or antagonistic mechanisms present in complex mixtures. These interactions would explain the null results found when EDCs are analyzed individually. 58 The E-SCREEN assay is a cell proliferation assay based on the enhanced proliferation of human breast cancer cells (MCF-7) in the presence of estrogen-active substances; thus, it primarily measures estrogenic activity but no other hormonal activities that could be relevant to the development of endometrial cancer. The assay’s specificity to estrogenic activity may overlook the opposing effect of progesterone or the influence of other hormones, such as androgens. 12 , 14 To attain a more comprehensive understanding of the hormonal context of endometrial cancer, future research could explore complementary assays or studies targeting other hormonal pathways and their potential interactions. Similar to findings shown in a previously reported breast cancer case–control study, 36 the combined estrogenic activity in serum extracts from cases and controls in this study was not associated with age, or with any other participant characteristic, except BMI. Interestingly, the associations were attenuated among obese women in stratified analyses. Obesity could act as a relevant confounder, given that it is associated with both case–control status and the levels of environmental lipophilic chemicals. The conversion of androgens to estrogens in adipose tissue could contribute, among others, to the role of obesity in endometrial cancer. 10 Aromatase levels increase with age and adiposity, and therefore this conversion results in estrogen-induced endometrial proliferation, as shown in some epidemiological studies. 59 Obesity, on the other hand, could lead to a lower serum concentration of persistent environmental chemicals, given that these compounds would be diluted in human adipose tissue. 48 However, obese controls in our study were more likely to have higher serum TEXB- α than nonobese controls, and adjusting estimates for BMI yielded stronger associations. Other studies, however, have reported no uniform relationship between chemical levels in serum and adipose tissue, highlighting that the adipose-to-serum ratios of these chemicals vary considerably among women. 60 Furthermore, obesity could potentially serve as a mediator in the present study, given the influence of certain persistent organic pollutants as obesogens leading to weight gain, 48 adding complexity to the analysis of the association between TEXB and endometrial cancer. We observed a positive correlation between serum TEXB- α and TEXB- β levels, and when we accounted for TEXB- β levels, we noticed a reduction in the strength of the association between TEXB- α and endometrial cancer risk. This association suggests that both TEXB- α and TEXB- β levels could be affected by an unspecified common exposure to both lipophilic and more polar xenoestrogens. 35 However, previous extensive assessment corroborated that lipophilic organohalogenated xenoestrogens that elute in the α -fraction are distinct from more polar xenoestrogens that elute, together with endogenous sex steroids, in the β -fraction . 35 TEXB- β levels are also influenced by the endogenous hormone levels of the participants at the time of blood collection, which could directly modulate the risk of endometrial cancer. Although the adjustment of TEXB- α estimates for TEXB- β can help mitigate the influence of endogenous hormones and correlated xenoestrogens found in the β fraction, this correction can also introduce M-bias or collider bias, resulting in a misleading attenuation of the results. 61 Recent literature using simulated data indicates that when a collider acts as a confounder, giving priority to confounding control should supersede the avoidance of M-bias. 62 Nonetheless, the effect of TEXB- α on endometrial cancer is expected to fall within the estimates with and without adjustment for TEXB- β . 61 Estrogenicity was below the LOD for TEXB- α in one-third of the serum samples and in 18.6% for TEXB- β . As stated above, the HPLC methodology was developed to separate lipophilic xenoestrogens ( α -fraction ) from sex steroid hormones and more polar xenoestrogens ( β -fraction ). 35 Many of the chemicals that elute in HPLC ( α -fraction ) have been subject to restrictive measures on their production and use owing to their persistence and suspected health effects. 63 Accordingly, TEXB- α levels in controls of the present study were lower than those in a previous similar population that was enrolled about a decade earlier and evaluated using the same methodology. 36 Consequently, the general population is exposed to relatively low but frequent doses of these compounds, which would explain the lower sensitivity of the bioassay to obtain a different signal than the negative control. This may justify why the serum TEXB values obtained in these participants are lower than those found in the multicase–control study, MCC-Spain. In that study, healthy women and women with breast cancer were enrolled between 2008 and 2013, and geometric means of 8.32 and 9.94  Eeq pM / mL were observed among cases for TEXB- α and TEXB- β levels, respectively, and 2.99 and 5.96  Eeq pM / mL , respectively, among controls. 36 The strengths of this study include the detailed assessment of covariates that allowed several sensitivity analyses and the first-time use of a reliable biomarker for the combined estrogenic effect of EDC mixtures in endometrial cancer risk. Our sample size was relatively large, considering the previous literature on xenoestrogens, as well as the serum-demanding and time-consuming characteristics of the TEXB bioassay. However, our sample size was too small to evaluate certain associations, especially for stratified analyses. Serum samples were collected at diagnosis for cases owing to the case–control design, which could lead to reverse causation bias if hormones or xenoestrogens levels changed after the onset of the disease. However, we expect this bias to be relatively low given that all serum samples were collected before any surgery, chemotherapy, or hormone therapy for cancer. Moreover, if serum estrogen depletion occurs as a result of tumors’ reliance on estrogen, this would lead to inverse estimates, further bolstering our conclusions regarding the positive associations. However, we used a single serum specimen as a surrogate for usual xenoestrogen exposures, and consequently, these exposures may not accurately represent the etiologic window for a cancer with a lengthy latent period. Endogenous estrogens and hydrophilic xenoestrogens primarily account for the hormonal activity observed in the β -fraction . Sexual hormones exhibit higher potency compared with nonpersistent xenoestrogen components that are detected within the β -fraction , such as BPA, parabens, and phthalates. Human exposure to these nonpersistent substances is usually assessed through repeated urine samples because they all have relatively short half-lives and are swiftly eliminated as metabolites in this biological matrix. Obtaining repeated samples from participants would have been highly beneficial in characterizing this exposure, particularly considering the long latency period of endometrial cancer. Adjustment variables were based on self-reported data, and therefore recall bias can result in residual confounding. The inclusion of controls with benign gynecologic conditions could lead to selection bias; these controls may have a different distribution of risk factors compared with the general population, especially those related to hormonal exposures (because they are often prescribed hormonal compounds), which could lead to attenuated estimates. Restricted analyses to controls without these benign gynecologic conditions, however, yielded consistent estimates. Hospital controls also may not be representative of the population hormonal burden, and further evaluation using prospective designs may yield additional insights.

We evaluated a reliable marker of the combined estrogenic effect of EDC mixtures in relation to endometrial cancer etiology in a Spanish case–control study. Both the combined estrogenic activity of mixtures of persistent xenoestrogens and endogenous hormones together with nonpersistent xenoestrogens in serum samples were positively associated with endometrial cancer risk with an inverted-U dose–response pattern. The use of in vitro bioassays on human samples may provide relevant insights in the way we understand the negative impact of chemical mixtures on human health effects.

Endometrial cancer ranks as the most common gynecological tumor in regions with a very high human development index. 1 The burden of this cancer is expected to increase worldwide owing to an increased aging population, given that it is typically a postmenopausal cancer, among others. 2 Endometrial cancer has been traditionally classified into two subtypes: type I (endometrioid) and type II (non-endometrioid). Type I tumors, the most common subtype, are typically low-grade and good-prognosis cancers. 3 Type II tumors constitute 10%–15% of endometrial cancers, are high-grade carcinomas, and show a poor prognosis. 3 Surrogates of the classification system developed by The Cancer Genome Atlas (TCGA) Consortium, which is based on molecular features, are now being integrated into clinical practice. 4 , 5 Endometrial cancer is a hormone-dependent cancer, and an estrogen imbalance plays a crucial role in its development. 6 Nulliparity, age at last birth, and infertility have repeatedly been associated with endometrial cancer. 7 – 9 Obesity can also increase the risk of endometrial cancer via the conversion of androgens to estrogens via aromatase activity in adipose tissue and other mechanisms. 10 A past history of breast cancer has also been associated with endometrial cancer, especially in estrogen receptor (ER)-positive breast cancer patients treated with tamoxifen. 11 Circulating estrogens and estrogen metabolites have been repeatedly associated with endometrial cancer, 12 – 14 as have genetic variants in sex hormone-related genes, such as CYP19. 15 , 16 Mendelian randomization studies have supported the hypothesis that endometrial cancer risk is causally driven by estrogens. 16 – 18 Both types of endometrial cancer share common risk factors, and to a lesser extent, the etiology of type II tumors has also been shown to have an estrogenic influence. 8 , 19 Although exposure to estrogens is an established risk factor for endometrial cancer, little is known regarding the effect of endocrine-disrupting chemicals (EDCs); that is, environmental pollutants that have the potential to act as xenoestrogens or influence estrogenic activity through different pathways. 20 Most of the literature originates from basic in vitro and in vivo data, 21 and the scarce observational literature has shown mostly null results, 22 – 28 with some exceptions. 29 – 31 However, the available literature suffers from limited sample sizes, has heterogeneous designs and exposure assessments, and focuses on a limited number of pollutants. Exposure assessments were performed using dietary variables, 22 – 24 , 32 urinary markers, 26 , 31 serum biomarkers, 27 , 28 adipose tissue biomarkers, 25 or endometrial and myometrial tissue biomarkers. 30 Functional biomarkers measuring the joint estrogenic activity of EDC mixtures offer a validated approach for combined exposure assessment. 33 – 36 Thus, a case–control study of 186 patients with breast cancer and 196 controls reported a positive association between the joint effect of environmental xenoestrogens extracted from serum samples and breast cancer risk, 36 which share certain etiologic similarities with endometrial cancer. However, to our knowledge no study has evaluated the combined effects of xenoestrogens on endometrial cancer. In the present study, we aimed to assess the association between the combined estrogenic activity of mixtures of xenoestrogens present in serum samples and the risk of endometrial cancer in the Screenwide Spanish case–control study.