Introduction: In Asian and Asian-American women, soy consumption is associated with reduced breast cancer risk, perhaps due to its effects on estrogen production or metabolism. In a sample of Asian-American women, we investigated the associations of usual adult soy intake with the urinary concentrations of 15 estrogens and estrogen metabolites (EM) measured using liquid chromatography-tandem mass spectrometry.

Methods: Participants included 430 Chinese-American, Japanese-American, and Filipino-American women, ages 20 to 55 years, and living in San Francisco-Oakland (California), Los Angeles (California), or Oahu (Hawaii). They were postmenopausal (n = 167) or premenopausal in luteal phase (n = 263) when 12-hour urine samples were collected. Robust linear regression was used to assess soy tertiles as predictors of log-transformed EM measures. Individual and grouped EM were considered as concentrations (pmol/mg creatinine) and as percentages of total EM (%EM).

Results: Factor analysis confirmed that EM groups defined by metabolic pathways appropriately captured covariation in EM profiles. Total EM concentrations were not significantly associated with soy in premenopausal or postmenopausal women. Among all women, %2-hydroxylated EM and %4-hydroxylation pathway EM were 16% higher (Ptrend = 0.02) and 19% higher (Ptrend = 0.03) in the highest versus lowest soy tertiles, respectively. In contrast, 16% hydroxylated EM were 11% lower (Ptrend < 0.01). Results were consistent across ethnic and menopausal groups and after adjustment for westernization measured by birthplace (Asia or United States).

Discussion: Findings suggest that regular soy intake is associated with increased ratios of 2:16-pathway EM and with higher relative levels of 4-hydroxylated EM. The observed variations in estrogen metabolism might modify breast cancer risk. (Cancer Epidemiol Biomarkers Prev 2009;18(10):2751–60)

Estrogens play important roles in the pathophysiology of breast tumors and are also recognized as causal factors in the etiology of this disease (1, 2). It has long been known that estrogens exert mitogenic effects on cells in the breast via receptor-mediated signaling (3). More recently, it has been recognized that some estrogen metabolites could be converted into reactive oxidative species which can damage DNA directly (4). Thus, estrogens and estrogen metabolites (jointly referred to as EM) might contribute to cancer risk by acting as cancer-promoting growth factors, or as cancer-initiating mutagens. Although prospective studies have shown consistent associations between circulating estrogen levels and risk of subsequent breast cancer in postmenopausal women, the relative importance of proposed pathways in breast carcinogenesis is not yet understood.

Bradlow et al. suggested that variability among women in susceptibility to breast cancer may result from interindividual variations in estrogen metabolism (5). Metabolism of estrogens occurs in the liver and kidneys as well as in breast and other target tissues, and includes oxidative metabolism (hydroxylation) and conjugative metabolism (glucuronidation, sulfation, glutathionation, and/or O-methylation; ref. 6). Hydroxylation of estrone and estradiol most commonly occurs at ring carbon positions 2, 4, or 16 to yield 2-hydroxylated EM, 4-hydroxylated EM, or 16-hydroxylated EM, respectively (7). EM may vary in binding affinity with estrogen receptors, susceptibility to oxidative conversion into reactive species, susceptibility to conjugative metabolism, and bioavailability in target tissues. Therefore, EM are likely to vary with respect to specific roles in cancer etiology or cancer prevention.

Soy foods are a dietary staple in East Asian countries and may play a role in this region's historically low incidence of breast cancer. Epidemiologic studies conducted to date among Asian women have found weak to moderate reductions in breast cancer risk associated with soy intake in adulthood (8, 9). In three studies which included measures of childhood and/or adolescent soy intake, results suggest stronger protective effects of soy consumed early in life (10-12). Isoflavones found in soy and soy-derived foods have structural similarities to estrogens and have been observed to have both estrogenic and antiestrogenic effects in animal models and in humans (13, 14). The health effects of isoflavones may be mediated through modulation of estrogen receptor signaling (15), but could also be mediated through their effects on enzymes with roles in production and/or metabolism of estrogens (16).

Various studies have tried to assess the effects of soy intake on endogenous levels of EM, including observational studies (17-19), short-term feeding studies (20-23), and randomized controlled trials (24). However, the studies have produced conflicting results.

We therefore decided to study the association between usual adult soy intake and urinary concentrations of 15 EM among Asian-American women. Study participants were of Chinese, Japanese, and Filipino ancestry and were living in San Francisco-Oakland (California), Los Angeles (California), and Oahu (Hawaii) when they were enrolled as controls in the population-based Asian-American Breast Cancer Study. In that study, we observed statistically significant inverse trends in breast cancer risk across soy tertiles; risk was 58% lower and 29% lower, respectively, in women in the highest versus lowest tertiles of childhood soy intake and usual adult soy intake (12, 25). In the present analysis, we test the hypothesis that adult intake of soy foods is associated with creatinine-adjusted urinary EM concentrations and/or relative proportions of urinary EM.

Study Population

Subjects for this cross-sectional analysis were drawn from a population-based breast cancer case-control study that has been described in detail previously (26). Briefly, cases in the parent study were Asian-American women of Chinese, Japanese, or Filipino ethnicity, living in San Francisco-Oakland (California), Los Angeles (California), and Oahu (Hawaii) who were diagnosed with breast cancer at 20 to 55 y of age during 1983 to 1987. Controls were frequency-matched to breast cancer cases by ethnicity, year of birth in 5-y age groups, and study center. These controls were the subjects for the current study. Controls completed in-home interviews which included questions about diet, anthropometry, menstrual and menopausal histories, and residential history. In addition, they were invited to provide overnight urine and/or fasting blood samples.

Of 966 controls in the parent study, 687 (71%) provided urine. Those who participated in the urine collection did not differ significantly from nonparticipants with respect to age, ethnicity, study center, or birthplace (Asia/West). Additionally, we excluded participants for whom urinary EM or urinary creatinine were not assayed (n = 16); participants who reported recent (<6 mo) exogenous hormone use (n = 83), pregnancy (n = 12), or lactation (n = 5); and participants with missing data on adult soy intake (n = 2).

Because menopausal status is an important determinant of circulating and urinary endogenous estrogen levels and so many of the participants in this study were ages 45 to 55 y and potentially perimenopausal, we assessed menopausal status carefully using data from multiple sources, including self-reported histories of menopause and surgical procedures from the baseline interview, surgical procedures reported in interviews at the time of urine collection, a postcard mailed by participants to study investigators recording the day of the first menstrual period following urine collection, and levels of follicle-stimulating hormone, progesterone, and estradiol in serum collected on the same day as urine. Decision rules were designed to identify women who were clearly premenopausal (having regular menstrual cycles or evidence of continuing ovulation) and women who were clearly postmenopausal (questionnaire data indicating cessation of menstrual cycles due to natural menopause at least 1 y prior to urine collection or circulating hormone levels indicating cessation of ovarian estrogen production). Postcard data and circulating progesterone levels were used to assess the phase of the menstrual cycle. Participants who reported continuing periods were given urine collection appointments that coincided with the midluteal phase (days 19-26) of the menstrual cycle. Because Xu et al. found that estrogen metabolism changes over the course of the menstrual cycle (27), premenopausal women thought to be in other phases of the menstrual cycle at the time of urine collection were excluded from most analyses.

Based on our criteria, of the 569 participants available for analysis, 167 women were classified as postmenopausal, 263 were classified as premenopausal in luteal phase, 98 were classified as premenopausal but not in luteal phase, and 41 women were assigned unknown/perimenopausal status due to missing or ambiguous information. Thus, we included 430 women, utilizing both premenopausal women in luteal phase and postmenopausal women in this study. Analyses were conducted separately for the two menopausal groups as well as combined. In sensitivity analyses, we also considered whether inclusion of premenopausal women who were not in luteal phase would modify study findings.

Soy Measures

During baseline interviews, participants were queried about their usual intake of approximately 60 food items/food groups, of which 10 to 15 were specific to their ethnic group. Frequencies of intake were recorded as times per day, week, month, or year, according to the most convenient time frame for the food item and respondent. Estimates of weekly soy intake for Chinese and Filipino women were based on reported frequency of intake of “a tofu dish made with any fresh, dried or deep-fried tofu product” whereas estimates for Japanese women were based on frequencies of intake of tofu, natto, and miso soup. Miso soup and natto were weighted as having 1/4 and 1/12, respectively, of the soy content of a tofu dish. Serving size was not queried. Participants were sorted into tertiles of soy intake using cutpoints generated from the distribution in the entire study sample (n = 430). Because soy intake was measured using more food items in Japanese compared with other ethnic groups, and because ethnic groups may differ both in EM profiles and in dietary patterns, all statistical models have been adjusted for ethnicity.

Urine Samples and Laboratory Methods

Participants were instructed to collect 12-h overnight urine samples using half-gallon containers kept at 4°C on ice or in the refrigerator. Boric acid was added as a preservative (28). After delivery to the study center, urine was mixed and aliquoted into conical tubes and sent to a repository for long-term storage at −80°C.

Aliquots were sent to Science Applications International Corporation (Frederick, MD), where 15 urinary EM were measured, including estrone (E1), estradiol (E2), estriol (E3), 2-hydroxyestrone (2-OHE1), 2-methoxyestrone (2-MeOE1), 2-hydroxyestradiol (2-OHE2), 2-methoxyestradiol(2-MeOE2), 2-hydroxyestrone-3-methyl ether (3-MeOE1), 4-hydroxyestrone (4-OHE1), 4-methoxyestrone (4-MeOE1), 4-methoxyestradiol (4-MeOE2), 16α-hydroxyestrone (16α-OHE1), 17-epiestriol (17-epiE3), 16-ketoestradiol (16-ketoE2), and 16-epiestriol (16-epiE3). The analytical method for measurement of urinary EM has been described previously (29). In brief, an internal standard solution containing four deuterium-labeled EM was added to each urine sample. Urinary EM and deuterium-labeled EM were hydrolyzed using β-glucuronidase/sulfatase from Helix pomatia, and then, following a 20-h incubation, extracted with dichloromethane. EM were quantitatively dansylated to improve their ionization efficiency. Liquid chromatography-tandem mass spectrometry analysis was done using a ThermoFinnigan TSQ Quantum-AM triple quadrupole mass spectrometer equipped with an electrospray ionization source and coupled directly to a Surveyor HPLC system (ThermoFinnigan). Both the chromatography system and the mass spectrometer were controlled using Xcalibur software (ThermoFinnigan).

The mass spectrometry conditions used in this study were as follows: source, positive mode ESI; spray voltage, 4,600 V; sheath and auxiliary gas, nitrogen; sheath gas pressure, 49 arbitrary units; auxiliary gas pressure, 23 arbitrary units; ion transfer capillary temperature, 350°C; scan type, selected reaction monitoring; collision gas, argon; collision gas pressure, 1.5 mTorr. The selected reaction monitoring conditions for the protonated molecules (MH+) of EM-dansyl and deuterium-labeled EM-dansyl were as follows: E1m/z 504 → 171 collision energy, 42 eV; E2m/z 506 → 171 collision energy, 43 eV; E3, 16-epi E3, and 17-epi E3m/z 522 → 171 collision energy, 43 eV; 16-ketoE2 and 16α-OHE1m/z 520 → 171 collision energy, 43 eV; 2-MeOE1, 4-MeOE1, and 3-MeOE1m/z 534 → 171 collision energy, 42 eV; 2-MeOE2 and 4-MeOE2m/z 536 → 171 collision energy, 43 eV; 2-OHE1 and 4-OHE1m/z 753 → 170 collision energy, 44 eV; 2-OHE2m/z 755 → 170 collision energy, 43 eV; d4-E2m/z 510 → 171 collision energy, 43 eV; d3-E3 and d3-16-epiE3m/z 525 → 171 collision energy, 43 eV; d5-2-MeOE2m/z 541 → 171 collision energy, 43 eV; d5-2-OHE2m/z 760 → 170 collision energy, 43 eV. The following mass spectrometry variables were used for all experiments: scan width, 0.7 u; scan time, 0.50 s; Q1 peak width, 0.70 u full width at half maximum; Q3 peak width, 0.70 u full width at half maximum.

Calibration curves for the 15 EM were constructed by plotting EM-dansyl/deuterium-labeled EM-dansyl peak area ratios obtained from calibration standards versus amounts of EM and fitting these data using linear regression with 1/X weighting. The amount of each EM in a urine sample was then interpolated using this linear function. Based on structural similarity and retention times, d4-E2 was used as the internal standard for E2 and E1; d3-E3 for E3, 16-ketoE2, and 16α-OHE1; d3-16-epiE3 for 16-epiE3 and 17-epiE3; d5-2-MeOE2 for 2-MeOE2, 4-MeOE2, 2-MeOE1, 4-MeOE1, and 3-MeOE1; d5-2-OHE2 for 2-OHE2, 2-OHE1, and 4-OHE1.

Quality control samples at three concentrations (0.12, 0.96, and 6.4 ng of each EM/mL) were included in each batch of assays. In addition, urine samples from two premenopausal and two postmenopausal women were used as blinded quality control samples. Four quality control samples, including two from the same subject, were randomly included in every batch of ∼40 samples. Total laboratory coefficients of variation were <10% for all EM except 4-methoxyestradiol (15%), and were ≤4% for estrone, estradiol, and estriol, and ∼1% for total EM. In general, laboratory coefficients of variation decreased as mean urinary EM concentration increased. Creatinine levels were obtained to adjust EM levels for differences in urine volume (30).

Statistical Methods

Statistical analyses were conducted using SAS v. 9.1 (SAS Institute, Cary, NC). Absolute levels of each urinary EM were expressed in picomoles per milligram of creatinine. Relative EM levels (%EM) were expressed as a percentage of total EM. Each of these measures was log-transformed to better approximate normal distributions.

To assess the effect of soy intake on EM measures, we first fit standard linear regression models and then robust linear regression models (PROC ROBUSTREG with option MM, which are less sensitive to outliers), to each EM measure as a dependent variable. We tested our hypotheses of linear trend across soy tertiles by assigning median frequency of soy intake to each tertile and treating the variable as a continuous covariate. We used the regression coefficient associated with the highest soy tertile to estimate the percent difference based on the formula: 100 [exp (β̂) − 1]. This represents a covariate-adjusted estimate of the difference between EM measures for the highest vs. the lowest soy tertile, expressed as a percent change.

Models of EM and %EM measures were fit separately for premenopausal and postmenopausal strata and were adjusted for age and ethnicity. In additional models we adjusted, in turn, for body mass index (BMI) and birthplace (Asia/West) to assess whether observed associations were independent of other aspects of acculturation. To assess effect modification, we fit models for EM grouped by metabolic pathways and expressed as a percentage of total, to sample subgroups defined by menopausal status, ethnicity, birthplace (Asia/United States), and BMI. The statistical significance of interactions between soy intake and ethnicity, BMI, and birthplace were evaluated by comparing models with and without interaction terms using a likelihood ratio test.

We used Spearman correlations and factor analysis to evaluate covariation among urinary EM (pmol/mg creatinine), separately in premenopausal-luteal and postmenopausal women. Because factor analysis confirmed that EM cluster according to metabolic pathways described in the literature, EM groups were defined to include EM deriving from specific metabolic events: hydroxylation (at carbon 2, 4, or 16 of the estrogen ring) and methylation. %EM groups were calculated by summing absolute levels of related metabolites and then dividing by total EM; these measures represent relative proportions of total urinary EM in a metabolic pathway. The composition of EM groups can be inferred from headings in Table 1, and are also described in the supplementary appendix.

Table 1.

Absolute and relative concentrations of urinary EM in premenopausal-luteal (n = 263) and postmenopausal (n = 167) Asian-American women

EM (pmol/mg creatinine)%EM
Premenopausal-lutealPostmenopausalPremenopausal-lutealPostmenopausal
MedianIQRMedianIQRMedianIQRMedianIQR
Parent estrogens 
    Estrone (E123.9 16.6-35.7 2.8 1.7-4.0 12.7% 10.0-15.8% 11.3% 7.3-15.2% 
    Estradiol (E210.4 7.8-15.9 1.6 0.9-2.7 5.9% 4.3-8.1% 6.3% 4.2-10.0% 
2-Hydroxylation pathway EM catechols 
    2-Hydroxyestrone (2-OHE126.5 16.2-40.9 3.2 1.8-5.0 14.3% 10.0-19.2% 12.5% 8.3-18.7% 
    2-Hydroxyestradiol (2-OHE22.8 1.6-4.8 0.6 0.3-1.3 1.5% 0.9-2.3% 2.3% 1.3-4.9% 
Methylated catechols 
    2-Methoxyestrone (2-MeOE15.2 3.3-7.9 0.7 0.4-1.0 2.7% 1.8-4.0% 2.7% 1.6-3.9% 
    2-Methoxyestradiol (2-MeOE20.5 0.3-0.8 0.1 0.1-0.2 0.3% 0.2-0.4% 0.5% 0.2-0.8% 
    2-Hydroxyestrone-3-methyl ether (3-MeOE10.9 0.7-1.7 0.3 0.2-0.6 0.6% 0.3-0.9% 1.3% 0.7-2.0% 
4-Hydroxylation Pathway EM catechols 
    4-Hydroxyestrone (4-OHE13.5 2.1-5.3 0.7 0.4-1.0 1.8% 1.2-2.9% 2.8% 1.7-3.9% 
Methylated catechols 
    4-Methoxyestrone (4-MeOE10.2 0.1-0.4 0.1 0.1-0.2 0.1% 0.1-0.2% 0.5% 0.2-0.8% 
    4-Methoxyestradiol (4-MeOE20.1 0.0-0.2 0.04 0.02-0.09 0.04% 0.02-0.08% 0.2% 0.1-0.4% 
16-Hydroxylation pathway EM 
    16α-Hydroxyestrone (16α-OHE111.7 7.3-17.3 1.3 0.8-2.1 5.9% 4.6-8.0% 5.3% 3.7-7.3% 
    Estriol (E362.3 36.2-90.7 5.6 3.6-9.4 32.2% 23.6-39.5% 23.6% 16.0-34.0% 
    17-Epiestriol (17-epiE31.1 0.6-2.4 0.2 0.1-0.3 0.6% 0.3-1.2% 0.6% 0.3-1.3% 
    16-Ketoestradiol (16-ketoE222.0 14.3-33.6 3.7 2.2-5.4 11.8% 9.1-15.0% 14.3% 10.5-20.3% 
    16-Epiestriol (16-epiE39.0 5.5-12.3 0.9 0.6-1.3 4.6% 3.6-5.5% 3.6% 2.6-4.6% 
Total EM 194.9 138.5-267.9 25.3 17.7-34.8     
EM (pmol/mg creatinine)%EM
Premenopausal-lutealPostmenopausalPremenopausal-lutealPostmenopausal
MedianIQRMedianIQRMedianIQRMedianIQR
Parent estrogens 
    Estrone (E123.9 16.6-35.7 2.8 1.7-4.0 12.7% 10.0-15.8% 11.3% 7.3-15.2% 
    Estradiol (E210.4 7.8-15.9 1.6 0.9-2.7 5.9% 4.3-8.1% 6.3% 4.2-10.0% 
2-Hydroxylation pathway EM catechols 
    2-Hydroxyestrone (2-OHE126.5 16.2-40.9 3.2 1.8-5.0 14.3% 10.0-19.2% 12.5% 8.3-18.7% 
    2-Hydroxyestradiol (2-OHE22.8 1.6-4.8 0.6 0.3-1.3 1.5% 0.9-2.3% 2.3% 1.3-4.9% 
Methylated catechols 
    2-Methoxyestrone (2-MeOE15.2 3.3-7.9 0.7 0.4-1.0 2.7% 1.8-4.0% 2.7% 1.6-3.9% 
    2-Methoxyestradiol (2-MeOE20.5 0.3-0.8 0.1 0.1-0.2 0.3% 0.2-0.4% 0.5% 0.2-0.8% 
    2-Hydroxyestrone-3-methyl ether (3-MeOE10.9 0.7-1.7 0.3 0.2-0.6 0.6% 0.3-0.9% 1.3% 0.7-2.0% 
4-Hydroxylation Pathway EM catechols 
    4-Hydroxyestrone (4-OHE13.5 2.1-5.3 0.7 0.4-1.0 1.8% 1.2-2.9% 2.8% 1.7-3.9% 
Methylated catechols 
    4-Methoxyestrone (4-MeOE10.2 0.1-0.4 0.1 0.1-0.2 0.1% 0.1-0.2% 0.5% 0.2-0.8% 
    4-Methoxyestradiol (4-MeOE20.1 0.0-0.2 0.04 0.02-0.09 0.04% 0.02-0.08% 0.2% 0.1-0.4% 
16-Hydroxylation pathway EM 
    16α-Hydroxyestrone (16α-OHE111.7 7.3-17.3 1.3 0.8-2.1 5.9% 4.6-8.0% 5.3% 3.7-7.3% 
    Estriol (E362.3 36.2-90.7 5.6 3.6-9.4 32.2% 23.6-39.5% 23.6% 16.0-34.0% 
    17-Epiestriol (17-epiE31.1 0.6-2.4 0.2 0.1-0.3 0.6% 0.3-1.2% 0.6% 0.3-1.3% 
    16-Ketoestradiol (16-ketoE222.0 14.3-33.6 3.7 2.2-5.4 11.8% 9.1-15.0% 14.3% 10.5-20.3% 
    16-Epiestriol (16-epiE39.0 5.5-12.3 0.9 0.6-1.3 4.6% 3.6-5.5% 3.6% 2.6-4.6% 
Total EM 194.9 138.5-267.9 25.3 17.7-34.8     

Abbreviation: IQR, interquartile range.

Finally, we calculated the ratios of EM which may be associated with breast cancer risk according to various hypotheses of estrogen-mediated carcinogenesis, including the ratio of 2-hydroxylated EM/16-hydroxylated EM, the ratio of 4-hydroxylated EM/2-hydroxylated EM, and the ratio of unmethylated/methylated catechols in 2-hydroxylation and 4-hydroxylation pathways.

%EM groups and EM ratios were each log-transformed for use as dependent variables in regression analyses. Regression models to assess the associations of soy intake with EM groups and ratios were fit to premenopausal and postmenopausal women combined and were adjusted for menopausal status, age, and ethnicity. In the combined population, we evaluated effect modification by other variables including menopausal status, using the approach outlined above. P ≤ 0.05 were considered statistically significant. All tests were two-sided.

In Table 1, medians and interquartile ranges are presented according to menopausal status for each urinary EM, expressed first as an absolute concentration (pmol/mg creatinine) and then as a proportion of total EM. Even within menopausal groups, urinary EM concentrations varied by factors of 10 to 100. Median EM concentrations were 2 to 10 times higher in premenopausal-luteal women compared with postmenopausal women. In contrast, interquartile ranges for %EM in premenopausal and postmenopausal women overlapped. Most median %EM differed by <2 percentage points between premenopausal-luteal and postmenopausal women. One exception is estriol, which had a markedly lower percentage of total EM in postmenopausal women compared with their premenopausal-luteal counterparts (23.6% versus 32.2%, respectively).

Participant characteristics by tertiles of soy intake are presented in Table 2 for premenopausal-luteal and postmenopausal women. Soy intake was nonsignificantly higher in postmenopausal compared with premenopausal-luteal women (P = 0.09), with postmenopausal women making up 32.7% of the lowest tertile, but 44% of the highest soy tertile. Soy intake was associated with ethnicity in both premenopausal and postmenopausal women. Filipino women had the lowest frequency of soy intake (with median intakes of 0.2 and 0.5 servings per week in premenopausal-luteal and postmenopausal women, respectively) and so were overrepresented in the lowest soy tertiles. In contrast, the median frequency of soy intake for each menopausal subgroup was approximately one serving per week in both Chinese-American and Japanese-American women. Soy intake was not significantly associated with age, study center, birthplace, or BMI in either menopausal group.

Table 2.

Participant characteristics by menopausal status and tertiles of usual adult soy intake

CharacteristicPremenopausal women in luteal phase (n = 263)Postmenopausal women (n = 167)
Soy tertiles*LowMediumHighP for trend/differenceLowMediumHighP for trend/difference
Weekly servings of soy<0.50.5-1.0≥1.1<0.50.5-1.0≥1.1
n1097579535064
Age (mean ± SD) 41.4 ± 4.3 42.4 ± 4.8 42.2 ± 5.4 0.30 54.2 ± 4.2 55.3 ± 3.5 54.3 ± 3.8 0.27 
BMI (mean ± SD) 20.8 ± 2.4 20.9 ± 2.6 21.2 ± 3.9 0.63 22.4 ± 3.7 22.5 ± 3.0 21.5 ± 2.6 0.15 
Ethnicity (%)    <0.001    <0.001 
    Filipino 40.4 16.0 7.6  56.6 26.0 9.4  
    Japanese 33.9 44.0 62.0  22.6 36.0 65.6  
    Chinese 25.7 40.0 30.4  20.8 38.0 25.0  
Study center (%)    0.72    0.66 
    Hawaii 41.3 48.0 38.0  39.6 42.0 46.9  
    Los Angeles 28.4 22.7 26.6  32.1 30.0 20.3  
    San Francisco 30.3 29.3 35.4  28.3 28.0 32.8  
Place of birth (%)    0.19    0.46 
    United States 55.6 45.3 59.5  54.7 54.0 50.0  
    Asia 44.4 54.7 40.5  45.3 46.0 50.0  
CharacteristicPremenopausal women in luteal phase (n = 263)Postmenopausal women (n = 167)
Soy tertiles*LowMediumHighP for trend/differenceLowMediumHighP for trend/difference
Weekly servings of soy<0.50.5-1.0≥1.1<0.50.5-1.0≥1.1
n1097579535064
Age (mean ± SD) 41.4 ± 4.3 42.4 ± 4.8 42.2 ± 5.4 0.30 54.2 ± 4.2 55.3 ± 3.5 54.3 ± 3.8 0.27 
BMI (mean ± SD) 20.8 ± 2.4 20.9 ± 2.6 21.2 ± 3.9 0.63 22.4 ± 3.7 22.5 ± 3.0 21.5 ± 2.6 0.15 
Ethnicity (%)    <0.001    <0.001 
    Filipino 40.4 16.0 7.6  56.6 26.0 9.4  
    Japanese 33.9 44.0 62.0  22.6 36.0 65.6  
    Chinese 25.7 40.0 30.4  20.8 38.0 25.0  
Study center (%)    0.72    0.66 
    Hawaii 41.3 48.0 38.0  39.6 42.0 46.9  
    Los Angeles 28.4 22.7 26.6  32.1 30.0 20.3  
    San Francisco 30.3 29.3 35.4  28.3 28.0 32.8  
Place of birth (%)    0.19    0.46 
    United States 55.6 45.3 59.5  54.7 54.0 50.0  
    Asia 44.4 54.7 40.5  45.3 46.0 50.0  

NOTE: For continuous variables, the P value is for trend and comes from linear regression. For continuous covariates, the P value is for difference and comes from a χ2 test.

*Cutpoints for soy tertiles were determined based on the distribution of soy intake in the combined sample (n = 430) of premenopausal-luteal and postmenopausal women.

We tested for associations of urinary concentrations of individual EM with soy intake by modeling each EM measure using robust linear regression and adjusting for age and ethnicity (Table 3). No statistically significant associations were seen between soy intake tertiles and urinary concentrations of E1, E2, E3, or total EM in either menopausal group. In postmenopausal women, urinary levels of 2-MeOE1 were directly associated with soy intake (Ptrend = 0.02). Also among postmenopausal women, urinary levels of 16-ketoE2 declined significantly across soy tertiles (Ptrend = 0.009). These associations remained statistically significant when birthplace or BMI were added to the model. Among premenopausal-luteal women, no significant associations were observed for absolute levels of any individual EM with soy intake.

Table 3.

Percent difference in concentrations of urinary EM among Asian-American women in the highest vs. lowest tertiles of soy intake, by menopausal status

EM measurePercent difference*
EM (pmol/mg creatinine)EM (% of total)
Premenopausal-luteal women (n = 263)Postmenopausal women (n = 167)Combined premenopausal-luteal and postmenopausal women (n = 430)
Highest vs. lowest tertile P trendHighest vs. lowest tertileP trendHighest vs. lowest tertileP trend
Parent EM 
    Estrone (E1−1.8 0.80 7.5 0.70 7.7 0.21 
    Estradiol (E2−7.7 0.45 5.9 0.68 0.6 0.94 
2-Hydroxylation pathway EM catechols 
    2-Hydroxyestrone (2-OHE11.7 0.93 7.8 0.85 14.7 0.06 
    2-Hydroxyestradiol (2-OHE2−4.1 0.79 14.5 0.51 15.9 0.15 
Methylated catechols 
    2-Methoxyestrone (2-MeOE114.7 0.31 51.8 0.03 34.7 <0.01 
    2-Methoxyestradiol (2-MeOE27.9 0.56 17.3 0.35 27.8 0.02 
    2-Hydroxyestrone-3-methyl ether (3-MeOE117.6 0.33 14.2 0.53 12.1 0.22 
4-Hydroxylation pathway EM catechols 
    4-Hydroxyestrone (4-OHE15.9 0.68 25.2 0.19 22.1 0.02 
Methylated catechols 
    4-Methoxyestrone (4-MeOE16.6 0.94 22.0 0.24 10.6 0.37 
    4-Methoxyestradiol (4-MeOE2−7.6 0.75     15.6 0.72 6.9 0.64 
16-Hydroxylation pathway EM 
    16α-Hydroxyestrone (16α-OHE1−7.5 0.33 −19.6 0.17 −6.6 0.28 
    Estriol (E3−7.5 0.30 −2.8 0.98 −5.2 0.36 
    17-Epiestriol (17-epiE32.6 0.90 −25.0 0.30 −8.4 0.51 
    16-Ketoestradiol (16-ketoE2−10.0 0.26 −29.7 0.01 −12.1 0.02 
    16-Epiestriol (16-epiE3−13.0 0.12 −16.6 0.25 −8.4 0.04 
Total EM −2.3 0.62 −2.3 0.62   
EM measurePercent difference*
EM (pmol/mg creatinine)EM (% of total)
Premenopausal-luteal women (n = 263)Postmenopausal women (n = 167)Combined premenopausal-luteal and postmenopausal women (n = 430)
Highest vs. lowest tertile P trendHighest vs. lowest tertileP trendHighest vs. lowest tertileP trend
Parent EM 
    Estrone (E1−1.8 0.80 7.5 0.70 7.7 0.21 
    Estradiol (E2−7.7 0.45 5.9 0.68 0.6 0.94 
2-Hydroxylation pathway EM catechols 
    2-Hydroxyestrone (2-OHE11.7 0.93 7.8 0.85 14.7 0.06 
    2-Hydroxyestradiol (2-OHE2−4.1 0.79 14.5 0.51 15.9 0.15 
Methylated catechols 
    2-Methoxyestrone (2-MeOE114.7 0.31 51.8 0.03 34.7 <0.01 
    2-Methoxyestradiol (2-MeOE27.9 0.56 17.3 0.35 27.8 0.02 
    2-Hydroxyestrone-3-methyl ether (3-MeOE117.6 0.33 14.2 0.53 12.1 0.22 
4-Hydroxylation pathway EM catechols 
    4-Hydroxyestrone (4-OHE15.9 0.68 25.2 0.19 22.1 0.02 
Methylated catechols 
    4-Methoxyestrone (4-MeOE16.6 0.94 22.0 0.24 10.6 0.37 
    4-Methoxyestradiol (4-MeOE2−7.6 0.75     15.6 0.72 6.9 0.64 
16-Hydroxylation pathway EM 
    16α-Hydroxyestrone (16α-OHE1−7.5 0.33 −19.6 0.17 −6.6 0.28 
    Estriol (E3−7.5 0.30 −2.8 0.98 −5.2 0.36 
    17-Epiestriol (17-epiE32.6 0.90 −25.0 0.30 −8.4 0.51 
    16-Ketoestradiol (16-ketoE2−10.0 0.26 −29.7 0.01 −12.1 0.02 
    16-Epiestriol (16-epiE3−13.0 0.12 −16.6 0.25 −8.4 0.04 
Total EM −2.3 0.62 −2.3 0.62   

*Percent difference = 100[exp(β̂) − 1], where β is the coefficient from robust linear regression corresponding to women in the highest compared with the lowest tertile of soy on log-transformed EM measures (pmol/mg creatinine or mol% of total EM). For menopause-specific strata, regression models were adjusted for continuous age and ethnicity (Chinese, Japanese, and Filipino); for the combined sample, regression models were also adjusted for menopausal status.

P for linear trend was tested by assigning median frequency of soy intake to each tertile and treating the variable as a continuous covariate.

P < 0.05, associated with tertile estimate.

In Table 3, we also present the results of robust linear regression models for individual EM expressed as a percentage of total EM, in a combined sample of premenopausal-luteal and postmenopausal women. 4%-OHE1, 2%-MeOE1, and 2%-MeOE2 increased significantly, whereas 16%-keto E2 and 16%-epiE3 each had statistically significant inverse trends across tertiles of soy intake. Other %EM in these pathways had similar but nonsignificant associations with soy. This observation led us to examine covariation among EM.

Factor analysis with rotation resulted in the extraction of four independent factors in each menopausal group (Fig. 1). The composition of factors was observed to be similar in premenopausal-luteal and postmenopausal women, in spite of the large differences between these groups in absolute levels of total EM. The predominant factor in both premenopausal and postmenopausal women, termed “catechols” and accounting for 46% and 30% of the variance in EM profiles, respectively, included all three catechol estrogens (2-OHE1, 2-OHE2, and 4-OHE1) and some of their five methylated metabolites (2-MeOE1 and 3-MeOE1 in both groups, and 2-MeOE2 in premenopausal women only). The second factor in both groups (termed “16-hydroxylated pathway EM”) included all five EM that were initially hydroxylated at the 16-C position (16-keto E2, 16a-OHE1, 16-epiE3, E3, and 17-epiE3), and accounted for 18% and 16%, respectively, of the variance in EM profiles in the two menopausal groups.

Figure 1.

Covariation among urinary EM in Asian-American women by menopausal status. Factor analysis with rotation resulted in the extraction of four independent factors with eigenvalues of >1 in premenopausal-luteal and postmenopausal Asian-American women. Factors were named according to predominant contributing EM. Each EM is graphed by loadings for the three factors accounting for the most variance in urinary EM concentrations. Solid shapes, statistically significant EM contributing to factors that correspond with axes in the graph; open shapes, EM contributing to the fourth extracted factor for each menopausal group (no axis shown).

Figure 1.

Covariation among urinary EM in Asian-American women by menopausal status. Factor analysis with rotation resulted in the extraction of four independent factors with eigenvalues of >1 in premenopausal-luteal and postmenopausal Asian-American women. Factors were named according to predominant contributing EM. Each EM is graphed by loadings for the three factors accounting for the most variance in urinary EM concentrations. Solid shapes, statistically significant EM contributing to factors that correspond with axes in the graph; open shapes, EM contributing to the fourth extracted factor for each menopausal group (no axis shown).

Close modal

The next factor, termed “parent estrogens”, was the third factor in premenopausal women, accounting for 7% of variation in EM profiles, and the fourth in postmenopausal women, accounting for 8% of variation in EM profiles. This factor included high loadings for E1 and E2 in both menopausal groups; in postmenopausal women only, two EM from the 16-pathway also contributed in a moderate fashion (E3 and 16-EpiE3). A fourth factor termed “methylated catechols” (4-MeOE1 and 4-MeOE2 in both groups, and 2-MeOE2 in postmenopausal women only) accounted for 6% and 14% of variation in premenopausal and postmenopausal women, respectively. The fact that methylated catechols in the 4-hydroxylation pathway did not cluster with those in the 2-hydroxylation pathway led us to define groups that examined these metabolites separately.

In Table 4, we present the results of robust linear regression models, fitted in each menopausal stratum and then in a combined sample, for EM grouped by metabolic pathway and expressed as a percentage of total EM. Patterns were similar in premenopausal-luteal and postmenopausal strata; accordingly, there was no statistically significant effect modification of the associations by menopausal status. In the combined sample, %2-hydroxylation pathway EM and %4-pathway EM each significantly increased across increasing tertiles of soy intake; these pathways were 16% higher (Ptrend = 0.02) and 19% higher (Ptrend = 0.03), respectively, in the highest compared with the lowest tertile of soy intake. The 16%-hydroxylated pathway EM were 11% lower in the highest versus the lowest tertile of soy intake (Ptrend < 0.01). %Parent estrogens showed no significant trend across soy tertiles. %Methylated catechols, like %catechols, increased significantly across the range of soy intakes (Ptrend < 0.01). Although %catechol estrogens in the 2-hydroxylation and 4-hydroxylation pathways each increased significantly with soy intake, trends in %methylated catechols differed by pathway. Methylated catechols of the 2-pathway increased as a proportion of total EM across soy tertiles (Ptrend < 0.01), whereas %methylated catechols in the 4-hydroxylation pathway showed no linear trend (Ptrend = 0.80).

Table 4.

Percent difference in urinary EM grouped by metabolic pathway in highest vs. lowest tertiles of soy intake among premenopausal (luteal) and postmenopausal Asian-American women

Percent difference*
Premenopausal-luteal (n = 263)Postmenopausal (n = 167)Combined (n = 430)
Highest vs. lowest tertileP trendHighest vs. lowest tertileP trendHighest vs. lowest tertileP trend
Grouped EM by pathway (% of total EM) 
    2-Hydroxylation pathway EM 13.5 0.06 20.0 0.16 16.3 0.02 
    4-Hydroxylation pathway EM 16.7 0.11 20.3 0.15 18.6 0.03 
    16-Hydroxylation pathway EM −7.1 0.05 −12.5 0.08 −10.6 <0.01 
    Parent estrogens 1.6 0.70 7.9 0.50 6.1 0.27 
    Catechol estrogens 10.9 0.10 15.3 0.35 15.0 0.05 
    Methylated catechols 22.8 0.06 35.8 0.02 23.9 0.01 
    Methylated catechols of the 2-hydroxylation pathway 25.8 0.04 38.9 0.02 27.4 <0.01 
    Methylated catechols of the 4-hydroxylation pathway 4.3 0.90 17.8 0.41 2.6 0.81 
EM ratios 
    4-Pathway/2-pathway −1.7 0.95 −1.3 0.98 −1.5 0.86 
    2-Pathway/16-hydroxylation pathway 21.9 0.05 42.0 0.10 31.4 0.01 
    4-Pathway/16-hydroxylation pathway 24.5 0.08 39.6 0.09 32.7 0.01 
    2-Pathway catechols/methylated 2-pathway catechols −8.74 0.49 −5.1 0.68 −2.3 0.80 
    4-Pathway catechols/methylated 4-pathway catechols 13.89 0.30 16.4 0.52 20.5 0.11 
Percent difference*
Premenopausal-luteal (n = 263)Postmenopausal (n = 167)Combined (n = 430)
Highest vs. lowest tertileP trendHighest vs. lowest tertileP trendHighest vs. lowest tertileP trend
Grouped EM by pathway (% of total EM) 
    2-Hydroxylation pathway EM 13.5 0.06 20.0 0.16 16.3 0.02 
    4-Hydroxylation pathway EM 16.7 0.11 20.3 0.15 18.6 0.03 
    16-Hydroxylation pathway EM −7.1 0.05 −12.5 0.08 −10.6 <0.01 
    Parent estrogens 1.6 0.70 7.9 0.50 6.1 0.27 
    Catechol estrogens 10.9 0.10 15.3 0.35 15.0 0.05 
    Methylated catechols 22.8 0.06 35.8 0.02 23.9 0.01 
    Methylated catechols of the 2-hydroxylation pathway 25.8 0.04 38.9 0.02 27.4 <0.01 
    Methylated catechols of the 4-hydroxylation pathway 4.3 0.90 17.8 0.41 2.6 0.81 
EM ratios 
    4-Pathway/2-pathway −1.7 0.95 −1.3 0.98 −1.5 0.86 
    2-Pathway/16-hydroxylation pathway 21.9 0.05 42.0 0.10 31.4 0.01 
    4-Pathway/16-hydroxylation pathway 24.5 0.08 39.6 0.09 32.7 0.01 
    2-Pathway catechols/methylated 2-pathway catechols −8.74 0.49 −5.1 0.68 −2.3 0.80 
    4-Pathway catechols/methylated 4-pathway catechols 13.89 0.30 16.4 0.52 20.5 0.11 

*Percent difference = 100[exp(β̂)−1], where β is the coefficient from robust linear regression corresponding to women in the highest compared to the lowest tertile of soy on log-transformed EM measures (pmol/mg creatinine or mol% of total EM). For menopause-specific strata, regression models were adjusted for continuous age and ethnicity (Chinese, Japanese, and Filipino); for the combined sample, regression models were also adjusted for menopausal status.

P for linear trend was tested by assigning median frequency of soy intake to each tertile and treating the variable as a continuous covariate.

P < 0.05, for tertile estimate.

Results for the ratios of EM were consistent with findings for EM groups. Ratios of 2-pathway EM and 4-pathway EM to 16α-pathway metabolites significantly increased across the range of soy intake (Ptrend = 0.01 for each). Among all women, the ratio of catechols/methylated catechols in the 4-pathway was 20% higher in the highest soy tertile compared with the lowest (Ptrend = 0.11); no trend was noted in the same ratio for the 2-pathway.

Finally, we considered several potential confounders (data not shown) of the EM/soy relationships. Birthplace is considered a marker of acculturation and of breast cancer risk (26). BMI might play a causal role in breast cancer etiology and has sometimes been observed to be inversely associated with soy intake (31). Additional adjustment of all models for birthplace and for BMI did not modify the direction or magnitude of observed associations.

Figure 2 shows the percentage of difference between the highest and lowest soy tertiles in relative levels of urinary EM for each stratum defined, in turn, by menopausal status, ethnicity, birthplace, and BMI. Results suggest that associations of %EM groups with soy intake were not significantly modified by these factors. One statistically significant interaction was noted: in women with BMI below the median, the %2-hydroxylation pathway showed a significant increasing trend across soy tertiles and was 16.6% higher in the highest versus lowest soy tertile; among heavier women, the %2-hydroxylation pathway also increased with soy intake but the magnitude of the effect was smaller, with a 10% increase seen in the highest versus lowest tertile (P for interaction = 0.04). This statistically significant finding may reflect the large number of statistical tests done.

Figure 2.

Percentage of difference between the highest and lowest tertiles of soy intake in relative levels of urinary EM (percentage of total) by menopausal status, ethnicity, and BMI. Percentage of difference was calculated based on regression coefficients associated with the highest tertile of soy intake in robust regression models fitted to predict log-transformed EM measures, using the formula 100[exp(β̂) − 1]. For menopause-specific strata, regression models were adjusted for continuous age and ethnicity (Chinese, Japanese, and Filipino); for ethnic-specific strata, regression models were adjusted for continuous age and menopausal status; for birthplace and BMI, models were adjusted for continuous age, ethnicity, and menopausal status. There was no statistically significant effect modification of the associations between soy and grouped %EM by menopausal status, by ethnicity or by birthplace. *, P < 0.05, for interaction.

Figure 2.

Percentage of difference between the highest and lowest tertiles of soy intake in relative levels of urinary EM (percentage of total) by menopausal status, ethnicity, and BMI. Percentage of difference was calculated based on regression coefficients associated with the highest tertile of soy intake in robust regression models fitted to predict log-transformed EM measures, using the formula 100[exp(β̂) − 1]. For menopause-specific strata, regression models were adjusted for continuous age and ethnicity (Chinese, Japanese, and Filipino); for ethnic-specific strata, regression models were adjusted for continuous age and menopausal status; for birthplace and BMI, models were adjusted for continuous age, ethnicity, and menopausal status. There was no statistically significant effect modification of the associations between soy and grouped %EM by menopausal status, by ethnicity or by birthplace. *, P < 0.05, for interaction.

Close modal

Inclusion of premenopausal women in other phases of the menstrual cycle also did not modify the direction or magnitude of the observed associations; however, it did result in wider confidence limits for many estimates.

In this study of premenopausal and postmenopausal Asian-American women with varying diets, lifestyles, and risks of breast cancer, we found no significant trends in absolute levels of estrone, estradiol, or total EM across tertiles of soy intake. This is consistent with the findings of a long-term intervention trial in which 220 premenopausal women were randomized to a soy intervention or control arm; in this study, investigators found no differences between groups in circulating estrone, estradiol, or free estradiol at four time points up to 24 months (24). In contrast, in two observational studies, investigators found that high soy intake was associated with reduced serum estrone in postmenopausal Chinese women (18) and reduced serum estradiol in Japanese premenopausal women (17). Evidence from short-term feeding trials (<6 months) has been inconsistent (16, 32, 33).

Although our findings do not support the hypothesis that soy intake influences endogenous production of estrogens, they do suggest that regular soy consumption modifies estrogen metabolism, perhaps through effects on enzymes involved in the hydroxylation of estrogens. In our study, women who reported more soy intake had relatively higher levels of 2- and 4-hydroxylated EM and relatively lower levels of 16-hydroxylated EM.

There are more than 10 human cytochrome P450 isoforms capable of hydroxylating estrone and estradiol to produce a wide range of potential metabolites (34). These phase I enzymes vary in distribution across target tissues, catalytic activity, and regiospecificity. The activities of several specific hepatic enzymes likely influence EM profiles; for example, CYP1A2 preferentially hydroxylates estrogens at the second carbon, whereas CYP3A4 and CYP3A5 have relatively higher activity for the production of 16-hydroxylated EM (35). The effects of polymorphic variants in extrahepatic enzymes, CYP1A1, which preferentially produces catechol estrogens, and CYP1B1, which selectively produces 4-hydroxylated EM, could also be detected in women's profiles of circulating or excreted EM (36, 37). Isoflavones can affect enzyme activity directly through competitive inhibition, or indirectly, through interactions with receptors that induce enzyme expression (38). In vitro studies suggest that exposure to soy isoflavones could modulate both the expression and activity of several of these enzymes; however, the net effects of soy exposure on estrogen metabolism in vivo remain unclear (39, 40).

Several observational studies have identified dietary or lifestyle components associated with levels of 2-OHE1, 16α-OHE1, or with their ratio (41, 42). The evidence regarding soy intake and levels of these EM is not conclusive. Of four intervention trials conducted in premenopausal women, two found an increased ratio of 2-hydroxyestrone/16α-hydroxyestrone during periods of exposure to high levels of soy isoflavones (16, 43) whereas two did not (44, 45). A feeding trial conducted with postmenopausal participants (33) showed a decreased ratio of 2-hydroxyestrone/16α-hydroxyestrone during a period of daily supplementation with 65 mg of soy isoflavones in comparison to measures at baseline and during use of the control supplement; but found no significant effect of 135 mg of soy isoflavones per day on the same parameter. This same study showed decreased levels of 4-OHE1 during use of isoflavone supplements at both doses. A more recent trial showed no effects of supplementation with soy flour on the ratio of 2-hydroxyestrone/16α-hydroxyestrone (46). Two studies found that soy interventions modified the ratio of 2-OHE1/16α-OHE1 only in a metabolic or ethnic subgroup (22, 47).

Methylation is an important pathway for detoxification and excretion of catechol estrogens (48). The enzyme catechol-O-methyl transferase is present in liver, kidney, and in breast tissue. Several in vitro studies have found that exposure of breast cells to soy isoflavones results in reduced expression and/or activity of catechol-O-methyl transferase (49, 50). Our findings suggest no significant difference in ratios of catechols/methylated catechols in the 2-hydroxylation pathway, but do show nonsignificant increases in this ratio for EM in the 4-pathway. Results of studies in human liver cells and in hamster kidney cells suggest that methylation by catechol-O-methyl transferase is more efficient for 2-hydroxylated EM than for 4-hydroxylated EM (51). Methylated 2-hydroxylated EM might actually inhibit the methylation of 4-hydroxylated catechols (52). Our findings suggest that soy isoflavones may directly or indirectly reduce the rate at which 4-hydroxyestrone is methylated.

Our study has several strengths. It includes a large population-based sample of Asian-American women with varying levels of acculturation, corresponding to a wide range of dietary and lifestyle choices, and a 6-fold gradient in breast cancer risk (26). Measures of usual adult soy intake represent long-term exposures at doses that are relevant to levels and patterns of intake in free-living populations. Another major strength of this study is the use of a highly sensitive, specific, and reliable assay to measure 15 of the most prevalent EM in urine. Many studies on the role of estrogen metabolism in breast cancer risk have used genetic variants as surrogates because these are easy to measure on a large scale; however, simulations suggest that genetic variants across metabolic pathways contribute in a complex fashion to the levels of EM (53). The EM profile is a phenotypic measure and thus provides a direct way to test hypotheses about the effects of dietary and lifestyle factors on estrogen metabolism and to study the effects of estrogen metabolism on disease risk.

The primary weakness of this observational study is the potential for confounding. We assessed two potential confounders and found that the observed associations were independent of birthplace (Asia/West), and of BMI, which are each associated with acculturation and breast cancer risk. Findings were also relatively consistent across ethnic subgroups, and are thus less likely to result from confounding by unmeasured dietary factors because dietary patterns are somewhat different in each ethnic group. Another study weakness is that assessment of soy intake was based on a limited set of queried soy foods, and on frequency of intake without reference to portion size. Results of a study of Asian-American women living in Los Angeles county (54), suggest that the queried foods may account for most of the variation in soy intake among Japanese participants, but only half of the variability among Chinese and Filipino participants.

Pooled analyses of prospective studies in postmenopausal women have provided definitive evidence that high levels of circulating estrogens (estrone, estradiol, and free estradiol) are associated with increased breast cancer risk (2). Our data suggests that if soy decreases breast cancer risk, its effects are not mediated by effects on absolute levels of estrogens or total EM. Increasing soy intake was associated with increased levels of 2-hydroxylated and 4-hydroxylated metabolites, and reduced levels of 16-hydroxylated metabolites. These variations in estrogen metabolism are consistent with a protective effect of soy on breast cancer risk based on Bradlow's theory of estrogen-mediated carcinogenesis (5). However, the “genotoxic metabolites” hypothesis, which posits oxidized products of catechol estrogens as primary causal agents in breast cancer (4), suggests that observed variations in estrogen metabolism could increase risk.

Methylated catechol estrogens have reduced estrogenicity, are readily cleared from circulation, and unlike the catechol estrogens from which they are derived, are not readily oxidized to form DNA-damaging quinones (6). The metabolite 2-methoxyestradiol has antiproliferative and antiangiogenic effects and therefore may play a role as an endogenous anticancer agent (55). Thus, observed trends of increasing methylated catechols in the 2-pathway across tertiles of soy intake could be consistent with the anticancer effects of soy.

Soy may also have effects on breast cancer risk mediated through other causal pathways; for example, through interactions of isoflavones with estrogen receptors. Further study of estrogen metabolism and EM phenotypes in relation to cancer risk will provide more information about the mechanisms by which soy can modify breast cancer risk.

No potential conflicts of interest were disclosed.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

We thank David Check of the Biostatistics Branch, Division of Cancer Epidemiology & Genetics, National Cancer Institute for assistance with figures.

1
Yager
JD
,
Davidson
NE
. 
Estrogen carcinogenesis in breast cancer
.
N Engl J Med
2006
;
354
:
270
82
.
2
Key
T
,
Appleby
P
,
Barnes
I
,
Reeves
G
. 
Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies
.
J Natl Cancer Inst
2002
;
94
:
606
16
.
3
Key
TJ
. 
Hormones and cancer in humans
.
Mutat Res
1995
;
333
:
59
67
.
4
Cavalieri
E
,
Chakravarti
D
,
Guttenplan
J
, et al
. 
Catechol estrogen quinones as initiators of breast and other human cancers: implications for biomarkers of susceptibility and cancer prevention
.
Biochim Biophys Acta
2006
;
1766
:
63
78
.
5
Bradlow
HL
,
Hershcopf
R
,
Martucci
C
,
Fishman
J
. 
16 α-Hydroxylation of estradiol: a possible risk marker for breast cancer
.
Ann N Y Acad Sci
1986
;
464
:
138
51
.
6
Raftogianis
R
,
Creveling
C
,
Weinshilboum
R
,
Weisz
J
. 
Estrogen metabolism by conjugation
.
J Natl Cancer Inst Monogr
2000
:
113
24
.
7
Jefcoate
CR
,
Liehr
JG
,
Santen
RJ
, et al
. 
Tissue-specific synthesis and oxidative metabolism of estrogens
.
J Natl Cancer Inst Monogr
2000
:
95
112
.
8
Wu
AH
,
Yu
MC
,
Tseng
CC
,
Pike
MC
. 
Epidemiology of soy exposures and breast cancer risk
.
Br J Cancer
2008
;
98
:
9
14
.
9
Trock
BJ
,
Hilakivi-Clarke
L
,
Clarke
R
. 
Meta analysis of soy intake and breast cancer risk
.
J Natl Cancer Inst
2006
;
98
:
459
71
.
10
Shu
XO
,
Jin
F
,
Dai
Q
, et al
. 
Soyfood intake during adolescence and subsequent risk of breast cancer among Chinese women
.
Cancer Epidemiol Biomarkers Prev
2001
;
10
:
483
8
.
11
Wu
AH
,
Wan
P
,
Hankin
J
, et al
. 
Adolescent and adult soy intake and risk of breast cancer in Asian-Americans
.
Carcinogenesis
2002
;
23
:
1491
6
.
12
Korde
LA
,
Wu
AH
,
Fears
T
, et al
. 
Childhood soy intake and breast cancer risk in Asian American women
.
Cancer Epidemiol Biomarkers Prev
2009
;
18
:
1050
9
.
13
Setchell
KD
,
Borriello
SP
,
Hulme
P
,
Kirk
DN
,
Axelson
M
. 
Nonsteroidal estrogens of dietary origin: possible roles in hormone-dependent disease
.
Am J Clin Nutr
1984
;
40
:
569
78
.
14
Kurzer
MS
,
Xu
X
. 
Dietary phytoestrogens
.
Annu Rev Nutr
1997
;
17
:
353
81
.
15
McCarty
MF
. 
Isoflavones made simple-genistein's agonist activity for the β-type estrogen receptor mediates their health benefits
.
Med Hypotheses
2006
;
66
:
1093
114
.
16
Xu
X
,
Duncan
AM
,
Merz
BE
,
Kurzer
MS
. 
Effects of soy isoflavones on estrogen and phytoestrogen metabolism in premenopausal women
.
Cancer Epidemiol Biomarkers Prev
1998
;
7
:
1101
8
.
17
Nagata
C
,
Kabuto
M
,
Kurisu
Y
,
Shimizu
H
. 
Decreased serum estradiol concentration associated with high dietary intake of soy products in premenopausal Japanese women
.
Nutr Cancer
1997
;
29
:
228
33
.
18
Wu
AH
,
Stanczyk
FZ
,
Seow
A
,
Lee
HP
,
Yu
MC
. 
Soy intake and other lifestyle determinants of serum estrogen levels among postmenopausal Chinese women in Singapore
.
Cancer Epidemiol Biomarkers Prev
2002
;
11
:
844
51
.
19
Verkasalo
PK
,
Appleby
PN
,
Davey
GK
,
Key
TJ
. 
Soy milk intake and plasma sex hormones: a cross-sectional study in pre- and postmenopausal women (EPIC-Oxford)
.
Nutr Cancer
2001
;
40
:
79
86
.
20
Cassidy
A
,
Bingham
S
,
Setchell
KD
. 
Biological effects of a diet of soy protein rich in isoflavones on the menstrual cycle of premenopausal women
.
Am J Clin Nutr
1994
;
60
:
333
40
.
21
Maskarinec
G
,
Morimoto
Y
,
Novotny
R
, et al
. 
Urinary sex steroid excretion levels during a soy intervention among young girls: a pilot study
.
Nutr Cancer
2005
;
52
:
22
8
.
22
Wu
AH
,
Stanczyk
FZ
,
Hendrich
S
, et al
. 
Effects of soy foods on ovarian function in premenopausal women
.
Br J Cancer
2000
;
82
:
1879
86
.
23
Wu
AH
,
Stanczyk
FZ
,
Martinez
C
, et al
. 
A controlled 2-mo dietary fat reduction and soy food supplementation study in postmenopausal women
.
Am J Clin Nutr
2005
;
81
:
1133
41
.
24
Maskarinec
G
,
Franke
AA
,
Williams
AE
, et al
. 
Effects of a 2-year randomized soy intervention on sex hormone levels in premenopausal women
.
Cancer Epidemiol Biomarkers Prev
2004
;
13
:
1736
44
.
25
Wu
AH
,
Ziegler
RG
,
Horn-Ross
PL
, et al
. 
Tofu and risk of breast cancer in Asian-Americans
.
Cancer Epidemiol Biomarkers Prev
1996
;
5
:
901
6
.
26
Ziegler
RG
,
Hoover
RN
,
Pike
MC
, et al
. 
Migration patterns and breast cancer risk in Asian-American women
.
J Natl Cancer Inst
1993
;
85
:
1819
27
.
27
Xu
X
,
Duncan
AM
,
Merz-Demlow
BE
,
Phipps
WR
,
Kurzer
MS
. 
Menstrual cycle effects on urinary estrogen metabolites
.
J Clin Endocrinol Metab
1999
;
84
:
3914
8
.
28
Falk
RT
,
Rossi
SC
,
Fears
TR
, et al
. 
A new ELISA kit for measuring urinary 2-hydroxyestrone, 16α-hydroxyestrone, and their ratio: reproducibility, validity, and assay performance after freeze-thaw cycling and preservation by boric acid
.
Cancer Epidemiol Biomarkers Prev
2000
;
9
:
81
7
.
29
Xu
X
,
Keefer
LK
,
Ziegler
RG
,
Veenstra
TD
. 
A liquid chromatography-mass spectrometry method for the quantitative analysis of urinary endogenous estrogen metabolites
.
Nat Protoc
2007
;
2
:
1350
5
.
30
Flores
OR
,
Sun
L
,
Vaziri
ND
,
Miyada
DS
. 
Colorimetric rate method for the determination of creatinine as implemented by the Beckman Creatinine Analyzer 2
.
Am J Med Technol
1980
;
46
:
792
8
.
31
Maskarinec
G
,
Aylward
AG
,
Erber
E
,
Takata
Y
,
Kolonel
LN
. 
Soy intake is related to a lower body mass index in adult women
.
Eur J Nutr
2008
;
47
:
138
44
.
32
Lu
LJ
,
Anderson
KE
,
Grady
JJ
,
Kohen
F
,
Nagamani
M
. 
Decreased ovarian hormones during a soya diet: implications for breast cancer prevention
.
Cancer Res
2000
;
60
:
4112
21
.
33
Xu
X
,
Duncan
AM
,
Wangen
KE
,
Kurzer
MS
. 
Soy consumption alters endogenous estrogen metabolism in postmenopausal women
.
Cancer Epidemiol Biomarkers Prev
2000
;
9
:
781
6
.
34
Lee
AJ
,
Cai
MX
,
Thomas
PE
,
Conney
AH
,
Zhu
BT
. 
Characterization of the oxidative metabolites of 17β-estradiol and estrone formed by 15 selectively expressed human cytochrome p450 isoforms
.
Endocrinology
2003
;
144
:
3382
98
.
35
Shou
M
,
Korzekwa
KR
,
Brooks
EN
, et al
. 
Role of human hepatic cytochrome P450 1A2 and 3A4 in the metabolic activation of estrone
.
Carcinogenesis
1997
;
18
:
207
14
.
36
Sowers
MR
,
Wilson
AL
,
Kardia
SR
,
Chu
J
,
McConnell
DS
. 
CYP1A1 and CYP1B1 polymorphisms and their association with estradiol and estrogen metabolites in women who are premenopausal and perimenopausal
.
Am J Med
2006
;
119
:
S44
51
.
37
Greenlee
H
,
Chen
Y
,
Kabat
GC
, et al
. 
Variants in estrogen metabolism and biosynthesis genes and urinary estrogen metabolites in women with a family history of breast cancer
.
Breast Cancer Res Treat
2007
;
102
:
111
7
.
38
Murray
M
. 
Altered CYP expression and function in response to dietary factors: potential roles in disease pathogenesis
.
Curr Drug Metab
2006
;
7
:
67
81
.
39
Scott
LM
,
Durant
P
,
Leone-Kabler
S
, et al
. 
Effects of prior oral contraceptive use and soy isoflavonoids on estrogen-metabolizing cytochrome P450 enzymes
.
J Steroid Biochem Mol Biol
2008
;
112
:
179
85
.
40
Anderson
GD
,
Rosito
G
,
Mohustsy
MA
,
Elmer
GW
. 
Drug interaction potential of soy extract and Panax ginseng
.
J Clin Pharmacol
2003
;
43
:
643
8
.
41
Wu
AH
,
Arakawa
K
,
Stanczyk
FZ
, et al
. 
Tea and circulating estrogen levels in postmenopausal Chinese women in Singapore
.
Carcinogenesis
2005
;
26
:
976
80
.
42
Sowers
MR
,
Crawford
S
,
McConnell
DS
, et al
. 
Selected diet and lifestyle factors are associated with estrogen metabolites in a multiracial/ethnic population of women
.
J Nutr
2006
;
136
:
1588
95
.
43
Lu
LJ
,
Cree
M
,
Josyula
S
, et al
. 
Increased urinary excretion of 2-hydroxyestrone but not 16α-hydroxyestrone in premenopausal women during a soya diet containing isoflavones
.
Cancer Res
2000
;
60
:
1299
305
.
44
Martini
MC
,
Dancisak
BB
,
Haggans
CJ
,
Thomas
W
,
Slavin
JL
. 
Effects of soy intake on sex hormone metabolism in premenopausal women
.
Nutr Cancer
1999
;
34
:
133
9
.
45
Cohen
LA
,
Crespin
JS
,
Wolper
C
, et al
. 
Soy isoflavone intake and estrogen excretion patterns in young women: effect of probiotic administration
.
In Vivo
2007
;
21
:
507
12
.
46
Brooks
JD
,
Ward
WE
,
Lewis
JE
, et al
. 
Supplementation with flaxseed alters estrogen metabolism in postmenopausal women to a greater extent than does supplementation with an equal amount of soy
.
Am J Clin Nutr
2004
;
79
:
318
25
.
47
Nettleton
JA
,
Greany
KA
,
Thomas
W
, et al
. 
The effect of soy consumption on the urinary 2:16-hydroxyestrone ratio in postmenopausal women depends on equol production status but is not influenced by probiotic consumption
.
J Nutr
2005
;
135
:
603
8
.
48
Goodman
JE
,
Jensen
LT
,
He
P
,
Yager
JD
. 
Characterization of human soluble high and low activity catechol-O-methyltransferase catalyzed catechol estrogen methylation
.
Pharmacogenetics
2002
;
12
:
517
28
.
49
Lehmann
L
,
Jiang
L
,
Wagner
J
. 
Soy isoflavones decrease the catechol-O-methyltransferase-mediated inactivation of 4-hydroxyestradiol in cultured MCF-7 cells
.
Carcinogenesis
2008
;
29
:
363
70
.
50
van Duursen
MB
,
Sanderson
JT
,
de Jong
PC
,
Kraaij
M
,
van den Berg
M
. 
Phytochemicals inhibit catechol-O-methyltransferase activity in cytosolic fractions from healthy human mammary tissues: implications for catechol estrogen-induced DNA damage
.
Toxicol Sci
2004
;
81
:
316
24
.
51
Hirose
Y
,
Tsutsui
TW
,
Ohno
M
,
Barrett
JC
,
Tsutsui
T
. 
Effects of a catechol-O-methyltransferase inhibitor on catechol estrogen-induced cellular transformation, chromosome aberrations and apoptosis in Syrian hamster embryo cells
.
Int J Cancer
2007
;
120
:
1627
33
.
52
Weisz
J
,
Clawson
GA
,
Creveling
CR
. 
Biogenesis and inactivation of catecholestrogens
.
Adv Pharmacol
1998
;
42
:
828
33
.
53
Crooke
PS
,
Ritchie
MD
,
Hachey
DL
, et al
. 
Estrogens, enzyme variants, and breast cancer: a risk model
.
Cancer Epidemiol Biomarkers Prev
2006
;
15
:
1620
9
.
54
Wu
AH
,
Yu
MC
,
Tseng
CC
,
Twaddle
NC
,
Doerge
DR
. 
Plasma isoflavone levels versus self-reported soy isoflavone levels in Asian-American women in Los Angeles County
.
Carcinogenesis
2004
;
25
:
77
81
.
55
Zhu
BT
,
Conney
AH
. 
Is 2-methoxyestradiol an endogenous estrogen metabolite that inhibits mammary carcinogenesis?
Cancer Res
1998
;
58
:
2269
77
.