Section 2
The NHSII is a prospective cohort study established in 1989 among 116,429 female registered nurses from 25 to 42 years of age. The enrolled women completed a questionnaire at study enrollment and are asked biennially to update with their medical and lifestyle factors and disease diagnoses. Between 1996 and 1999, blood samples were collected among a subgroup of participants in the NHSII. Details of the blood collection procedure are described in a prior publication [ 18 ]. Among the 29,611 participants who donated blood, 18,521 women who were not taking oral contraceptives, pregnant, or breastfed within 6 months providing timed samples collected in the follicular phase (3 to 5 days after the start of the menstrual cycle) and luteal phase (an estimated 7 to 9 days prior to the onset of the next menstrual cycle). Participants separated the plasma for the follicular sample and stored it in the freezer until the luteal sample was collected. Participants who provided untimed samples were included in analyses of dehydroepiandrosterone (DHEA), DHEA sulfate (DHEAS), testosterone, free testosterone, androstenedione, sex hormone–binding globulin (SHBG), and prolactin, as the concentrations of these hormones do not fluctuate much during menstrual cycles [ 19 ]. All blood samples were shipped overnight with an ice pack, processed by our laboratory into plasma, red blood cells, and white blood cells, and then stored in continuously monitored liquid nitrogen freezers (≤−130°). The original vials were thawed once to subaliquot a sample to send to each respective laboratory, and a second thaw occurred at the time of assay. The collection, processing, and storage methods of blood samples in NHSII have been shown to not materially alter levels of sex hormones. Furthermore, the between-person variation for most of the sex hormones far outweighed their within-person variation even with up to 2 days of delay in processing after phlebotomy [ 20 , 21 ].
Participants included in this analysis were part of a hormone reproducibility study [ 19 ], or were the controls in nested case-control studies of breast cancer [ 22 ], ovarian cancer [ 23 ], endometriosis [ 24 ], and rheumatoid arthritis [ 25 ]. A total of 2000 premenopausal women had measured concentrations of plasma sex hormones (estradiol, estrone, estrone sulfate, progesterone, DHEA, DHEAS, testosterone, androstenedione), prolactin, and SHBG concentrations. Among these participants, 1989 had flavonoid consumption data available in 1995, 1999, or both. The study protocol was approved by the Institutional Review Boards of the Brigham and Women’s Hospital and the Harvard T.H. Chan School of Public Health.
The laboratory assay methods used to obtain the concentrations of each hormone were documented in prior publications [ 5 , 26 ]. Estrone, estradiol, and estrone sulfate concentration were assayed in luteal and follicular samples; testosterone, androstenedione, and prolactin concentrations were assayed in luteal and/or follicular, as well as untimed samples; DHEA, DHEAS, and SHBG were assayed in luteal and untimed samples; and progesterone was assayed in luteal samples only. The assays were performed in different batches at different laboratories for several different sets of nested case-control studies. Three batches of estrogens, five batches of testosterone, two batches of androstenedione, and one batch of progesterone were assayed at Quest Diagnostics (San Juan Capistrano, CA, USA). At this location, estrogens and testosterone were measured by organic extraction, celite column chromatography, and then radioimmunoassay (RIA); estrone sulfate was measured by extraction of estrone, enzyme hydrolysis, column chromatography, and then RIA of estrone; progesterone were measured by organic extraction and RIA. Four batches of estrogens and testosterone were assayed at Mayo Medical Laboratories (Rochester, MN, USA) using liquid chromatography–tandem mass spectrometry (LC-MS/MS). Two batches of DHEA and androstenedione and four batches of DHEAS, SHBG, and progesterone were assayed at the Royal Marsden Hospital (London, United Kingdom). At this location, androstenedione was measured by RIA; and DHEAS, SHBG, and progesterone were measured by chemiluminescent enzyme immunoassay. One batch of progesterone, three batches of SHBG, and prolactin were assayed at Massachusetts General Hospital (Boston, MA, USA), where progesterone was measured by RIA, SHBG was measured by chemiluminescent enzyme immunoassay, and prolactin was measured by microparticle enzyme immunoassay (AxSYM Immunoassay System). One batch of SHBG and progesterone was assayed at Boston Children’s Hospital (Boston, MA, USA). Batch correction was performed as descibed in the statistical method section.
Free and percent free estradiol were calculated by the law of mass action according to the method described by Sodergard et al. [ 27 ]. The detection limits of the assays were as follows: 2 pg/mL for estradiol, 10 pg/mL for estrone, 40 pg/mL for estrone sulfate (in each laboratory), 3 ng/dL for androstenedione, 1 ng/dL for testosterone, 3 ng/dL for DHEA, 5 μg/dL for DHEAS, and 3 ng/dL for progesterone. When plasma hormone values were reported as less than the detection limit, we set the value to half this limit. All plasma hormone concentrations under the detection limits were set to half of the minimum detection value.
The reproducibility of plasma estrogens (intraclass correlation coefficients, ICC: 0.38–0.69), androgens (ICC: 0.56–0.94), progesterone (ICC: 0.29), prolactin (ICC: 0.64), and sex hormone binding globulin (ICC: 0.86) has been documented [ 19 ].
Flavonoid intake was calculated from semiquantitative food frequency questionnaires (FFQ) administered closest to the blood collection to represent the average intake around the time between 1995 and 1999. The participants were asked to estimate their usual intake of a list of ~131 food items, from which the intake of flavonoids were calculated using a database maintained by Harvard dietitians (available for request at: https://regepi.bwh.harvard.edu/health/nutrition.html ) that is largely based on the United States Department of Agriculture (USDA) food content database [ 28 ]. In the study that validated the FFQ against dietary records, the Spearman rank correlation coefficient is 0.55 for energy-adjusted total flavonoid intake, and from 0.41 to 0.58 for flavonoid subclasses [ 29 ]. For this study, we used average energy-adjusted intake of the nutrients from 1995 and 1999, or a single measure in 1995 or 1999 if there was only one available ( n = 192). We evaluated total flavonoids in our primary analyses. As secondary analyses, we examined 6 flavonoid subclasses-flavonols, flavanones, flavones, flavan-3-ols, anthocyanins, and isoflavones.
From the questionnaire completed at blood collection or the questionnaire immediately prior to blood draw, we obtained data on potential confounders for the association between flavonoid intake and plasma sex hormone concentrations: age at blood collection, smoking status, body mass index (BMI), height, age at menarche, usual menstrual cycle length, menstrual cycle pattern, parity, family history of breast cancer, history of biopsy-confirmed benign breast disease, physical activity, date and time of day of blood collection, fasting status, luteal day of blood collection, average alcohol consumption, weight change from age 18, and total energy intake.
Estradiol, free estradiol, estrone, and estrone sulfate were evaluated in follicular and luteal phases; progesterone was only evaluated in luteal phase; DHEA, DHEAS, testosterone, free testosterone, androstenedione, and prolactin were measured using the average of follicular and luteal phases or from a single luteal or untimed sample. We identified statistical outliers in the hormone data based on the generalized extreme studentized deviate (ESD) many-outlier detection approach [ 30 ]; the outliers were excluded from analyses. For each of the hormones, no more than 8 outliers were excluded. To address differences that might arise from assays being performed at multiple locations in multiple batches, all hormones were recalibrated to have a comparable distribution between all batches using the method of Rosner and colleagues [ 31 , 32 ]. In brief, we made statistical adjustment where the concentrations were combined to represent an average batch, then regressed on age, BMI (the strongest predictor of hormone levels), and indicator variables for each batch. Within each batch, the hormone levels were recalibrated according to the difference between the coefficients for that batch and the average of the batch coefficients. The adjusted concentrations thus accounted for between-batch variability independent of varying distributions of age and BMI across batches.
In our primary analyses, we calculated the geometric mean of each hormone concentration by quartile of total flavonoid intake using generalized linear models, with marginal covariate distribution in the study population. A contrast comparing the highest to the lowest quartile of flavonoid intake was also estimated. Confounders included in the models were based on subject matter knowledge and prior literature—age at blood collection in years (continuous), BMI at blood collection in kg/m 2 (<22.5, 22.5 to <25, 25 to <27.5, 27.5 to <30, ≥30), total energy intake in kcal/day (continuous), smoking (never smoker, past smoker with ≥5 years since quitting, past smoker with <5 years since quitting, current smoker of 20 g/day), physical activity (<3, 3–<9, 9–<18, 18–<27, ≥27 metabolic equivalent of task (MET)-hours/week), duration of past oral contraceptive use (never, <4 years, ≥4 years), age at menarche (13 years), age at first birth/parity (nulliparous, 1–2 children and age at first birth 3 children and age at first birth 3 children and age at first birth ≥25 years). Covariates related only to the outcomes—date of blood collection (month/year, 1/97, 2/97–1/98, 2/98–1/99, after 2/99), time of day of blood collection (1–8 a.m., 9 a.m. to noon, 1–4 p.m., 5 p.m. to midnight), and fasting status (<10, ≥10 h)—were also included in the models to reduce extraneous variation. Models that included luteal or average of timed samples also were adjusted for the difference between luteal blood draw date and date of next menstrual period (3–7, 8–12, 13–17, 18–21 days, unknown/untimed). We also evaluated the linear trend of flavonoid intake using median of each quartile modeled continuously and calculated a Wald test statistic. In secondary analyses, we repeated the main analyses for each subclass of flavonoids-flavonols, flavanones, flavones, flavan-3-ols, anthocyanins, and isoflavones. We also conducted a priori sensitivity analyses by excluding women who had luteal progesterone level <400 ng/dL ( n = 230, 11.6%) to restrict the analysis to women who donated blood samples during an ovulatory cycle.
Intro
Flavonoids, a diverse group of polyphenols, consist of more than 5000 different compounds produced in plants [ 1 , 2 ]. They have been considered one of the active components in fruits and vegetables that have antioxidant, anti-inflammatory, and anti-proliferative properties [ 3 ].
The initial interest in studying flavonoids stemmed from the observation that Asian populations have higher intake of soy products, which are high in flavonoids, particularly isoflavones, and lower breast cancer incidence rates compared to Western populations. As circulating estrogens and androgens are positively associated with the risk for breast cancer in pre- and postmenopausal women [ 4 , 5 , 6 , 7 ], flavonoids are also of special interest in relation to breast cancer due to the structural resemblance to estrogens [ 8 ], yet its consequent risk and benefit are still controversial. In vitro studies have shown promising anti-cancer effects of flavonoids [ 9 ]. In observational studies, however, the association is less clear, as most significant associations have been reported in Asian populations, in case-control studies, or for postmenopausal breast cancer, but generally not in Western populations, prospective cohort studies, or for premenopausal breast cancer [ 10 , 11 , 12 ].
The discrepancies between in vitro studies and human studies suggest that flavonoids are involved in more complex biological processes in vivo, yet the direction in which flavonoids affect sex steroid hormone-dependent activities is far from conclusive. Soymilk that is rich in isoflavone reduces circulating ovarian steroids and adrenal androgens and increases menstrual cycle length in premenopausal women [ 13 , 14 ]; while other subclasses of flavonoid have been shown to have binding affinities to estrogen receptor-α (ERα), and possible stimulation of ER-dependent transcriptional response that promotes growth of estrogen-dependent cells was implicated [ 15 , 16 , 17 ]. The effect of total flavonoid intake on risk of breast cancer is a combined result of the relative strengths in either direction, as well as the overall effect of all subclasses; it is important to elucidate whether any subcalsses of flavonoid modulate circulating estrogen and other sex hormones. To date, no study has characterized the relationship between dietary intake of specific subclasses of flavonoids and blood sex hormone levels as they vary over the menstrual cycle. In this study, we investigated the relationship among premenopausal women using timed samples in both follicular and luteal phases in the Nurses’ Health Study II (NHSII).
Results
The mean age at blood draw of the 1989 participants was 42.7 years, and the mean BMI was 26.1 kg/m 2 . Flavonoid intake did not vary substantially by height, BMI at age 18 or at blood draw, or reproductive factors. Those who had higher intake of total flavonoids tended to be more physically active, and less likely to be current smokers ( Table 1 ).
Table 2 shows the median concentration and 10th–90th percentile of the plasma sex hormones, prolactin, and SHBG. The Spearman rank correlation coefficients between the concentrations of the hormones, prolactin, and SHBG ranged from 0.0007 for follicular estradiol and androstenedione, to 0.88 for follicular estradiol and follicular-free estradiol.
In our study, the major sources of food for flavonoids were tea and apples, making up over 50% of the total flavonoid intake. The top contributing food items for the six flavonoid subclasses were flavonol (tea and onions), flavone (orange and orange juice), flavanone (orange and orange juice), flavan-3-ol (tea and apples), anthocyanin (blueberries and strawberries), isoflavone (tofu and soymilk) ( Table S3 ). Flavna-3-ol had foods with the highest intake among the six subclasses. The correlations between the subclasses of flavonoid are shown in Supplemental Table S1 . Flavonol and flavan-3-ol intakes were highly correlated, and since flavan-3-ol had the highest intake level, both were also highy correlated with total flavonoid intake (Spearman correlation coefficients 0.75–0.97, p < 0.0001). The correlations between other pairs of flavonoids subclasses were relatively weak (Spearman correlation coefficients <0.25) ( Supplementary Materials Table S1 ).
In this study, total flavonoid intake was not significantly associated with any of the sex hormone concentrations except for DHEAS ( p = 0.02); plasma DHEAS was 11.1% lower comparing the highest vs the lowest quartile of flavonoid intake (95% confidence interval (CI): −18.6%, −3.0%) ( Table 3 ). The percent difference comparing the highest versus the lowest quartile for the remaining sex hormones ranged from −7.4% to 5.1%.
In secondary analyses, we examined the six subclasses of flavonoids. Similar to total flavonoids, only DHEA and DHEAS showed significant differences between the extreme quartiles of the intake of the subclasses of flavonoid ( Table 4 ). Interestingly, DHEAS was significantly lower (−11.3% with 95% CI: −18.3, −3.7%, p -trend = 0.01) comparing the extreme quartiles of flavan-3-ol intake, however in contrast, DHEAS tended to increase significantly across quartiles of flavanone intake, although the difference between the highest and the lowest quartiles were only marginally significant (9.0% with 95% CI: −0.3%, 19.3%, p -trend = 0.04). Across the quartiles of flavan-3-ol intake, luteal progesterone concentration showed a significant increasing trend ( p = 0.03), although the contrast between the highest and the lowest quartile was not statistically significant (−7.6%, 95% CI: −14.78%, 0.20%). Anthocyanin intake was associated with significant inverse trends for DHEA (−18.0% with 95% CI: −27.9%, −6.7%, p -trend = 0.003) and follicular estrone sulfate, although the contrast between the extreme quartiles was not statistically significant for follicular estrone sulfate (−16.9% with 95% CI: −32.5%, 2.4%, p -trend = 0.02). Hormones with concentrations significantly higher comparing the highest and the lowest quartile of isoflavone were follicular estrone sulfate (25.7% with 95% CI: 3.6%, 52.6%) and DHEAS (12.7% with 95% CI: 3.5%, 22.8%), although neither had significant linear trend across quartiles of isoflavone intake.
In sensitivity analyses excluding anovulatory women (11.3%, n = 203), the results did not change materially. Total flavonoid intake remained significantly associated with DHEAS ( p -trend = 0.04). Among the six subclasses of flavonoids, associations that remained significant were between and flavan-3-ols and DHEAS ( p -trend = 0.03), and between anthocyanins and DHEA ( p -trend = 0.01).
Discussion
In this study of 1989 premenopausal women, we did not observe strong associations between intake of total flavonoids and/or its subclasses and most of the sex hormones examined. The only significant trends were that women with higher flavonoid intakes, especially that of flavan-3-ol and anthocyanins, had lower plasma DHEA and DHEAS concentrations; this has been associated with decreased risk of premenopausal breast cancer in the same cohort [ 33 ], as well as in a large pooled analysis [ 6 ].
Flavonoids in general have been considered as phytochemicals with the potential for antioxidant and anti-inflammatory activity [ 34 , 35 ], although the antioxidant capacity is limited after extensive metabolism that is sufficient to reduce plasma indices of oxidant status in vivo. In vitro studies suggest that they also have the promise as anti-tumor agents [ 3 ]-a number of flavonoids can induce apoptosis in cancer (including breast cancer) cell lines, arreste proliferation at the G2/M cell cycle checkpoint [ 36 , 37 ], and bound strongly with Estrogen Receptor β, (ER-β), which opposes the proliferative effects of ERα activation [ 38 ]. However, due to their estrogenic feature, some flavonoids also have binding affinity to ER-α, indicating possible promotion of proliferation [ 16 ]. The health benefit of even the most promising subclass-isoflavone-was challenged, given that soy isoflavone extracts may promote estrogen receptor positive breast cancer cell growth [ 39 ]. In addition, the majority of the evidence from lab-based experimental studies did not take into consideration the absorption and metabolic interactions within the human body. Epidemiologic studies that examined the relationship between flavonoids (dietary intake and circulating concentrations, as well as urinary levels) and breast cancer risk among humans have shown inconsistent results so far [ 11 ]. For example, short-term supplementation of soymilk rich in isoflavones decreased circulating estrogens among premenpausal women, and isoflavones were associated with lower premenopausal breast cancer risk in Asian and German populations [ 40 ]. However, studies also suggest that soy (the main source of isoflavone) supplementation may increase the expression of several tumor-promoting genes [ 41 ], and the breast cancer risk reduction by isoflavone was not seen in most Western countries.
A meta-analysis published in 2013 synthesized multiple prospective cohort and case-control studies on flavonoid intake and breast cancer risk. While total flavonoid intake, especially that of flavonols and flavones, was inversely associated with postmenopausal breast cancer, the results for premenopausal women were non-significant [ 42 ]. Among all proposed effects of flavonoid, most (e.g., antioxidation, anti-inflammation, inhibition of carcinogen-activating enzymes) are independent of menopausal status. What is most plausible to explain the different effects of flavonoid on pre and postmenopausal breast cancer is the complicated interaction between flavonoid and sex hormones. One may be concerned that flavonoids might alter sex hormones in premenopausal women in the undesired direction thus counteract the overall benefit. Our results are reassuring that flavonoids are not associated with higher concentrations of circulating sex hormones. In fact, DHEA and DHEAS demonstrated inverse linear trends with higher intake of anthocyanins and flavan-3-ols. In vitro studies have shown that DHEA increases proliferation in the human breast cancer cell line MCF-7 at physiological concentrations and has an antiproliferative effect at supraphysiologic concentrations [ 43 , 44 , 45 , 46 ]. DHEAS also has the potential to induce growth of ER-positive/AR-positive cells [ 47 ].
It is possible that a greater range of intake would result in further reduction of the hormones, and may lead to a decreased risk of breast cancer. In contrast, flavanones were associated with higher circulating DHEAS level. The higher intake of flavanones (median = 27.41 mg/d) compared to flavan-3-ols (median = 23.16 mg/d) and anthocyanins (median = 10.18 mg/d), may partially offset the lower concentrations shown with higher anthocyanin and flavan-3-ol intake, and more importantly emphasizes that the subclasses should be analyzed separately. Nevertheless, alternative mechanisms are needed to explain the heterogeneity due to menopausal status. Circulating sex hormone levels are relatively high before menopause, and therefore the degree to which dietary flavonoids can act against endogenous sex hormones may be limited; whereas after estrogen levels fall at menopause, the protective effect of flavonoids may become more evident.
Our study has several strengths. In a large sample of premenopausal women we assessed the association between multiple sex hormones and flavonoid intake by subclasses. The sex hormones that fluctuate during the menstrual cycle were timed, that is, measured in both follicular and luteal phases, adding accuracy to the analyses. We were also able to take advantage of high-quality dietary assessment based on validated FFQ, and had detailed data on a range of potentially important covariates. There are also limitations of our study. Measurement errors in the assessment of flavonoid intake could attenuate the associations, if any; meanwhile, the nutrient calculations for food items using the United States Department of Agriculture (USDA) food content database is imperfect. However, the relative accuracy of the exposures in our study is acceptable, since validation of FFQ versus 7-day dietary records suggested that the questionnaire was able to represent the average intake over the yearly period. The measurement error of the sex hormones examined should also be considered. The validity depends on both the assays performed at different labs and the reproducibility of a single measure. Notably, DHEAS, DHEA, and SHBG had the highest ICC among all hormones we assessed. The true associations between flavonoid intake and the other hormones (with lower ICCs) could have been undetectable because they were attenuated by random within-person variation. For samples that were untimed, the Spearman correlations between samples collected during the follicular and luteal phase were 0.60 for DHEA, 0.82 for DHEAS, 0.63 for testosterone, 0.55 for prolactin, and 0.86 for SHBG in this population. Therefore, any fluctuation might also have attenuated the observed association. We did not perform formal multiple testing adjustment, but given the number of statistical tests we conducted, some of the significant associations could be due to chance. If adjusted for multiple comparisons of the hormones, only the association between anthocyanins and DHEA would remain significant. We were also not able to control for all potential confounders, such as other nutrients rich in fruits and vegetables. However, in our study, the Pearson correlation coefficients between total flavonoid intake and intake of dietary fiber, carotenoids, total fruit, and total vegetables range from 0.06 to 0.11, suggesting that confounding by these nutrients/foods seems unlikely. Our subclass analyses, though, suggest that most of the significant association between flavonoids and DHEAS was driven by flavan-3-ol, which is enriched in tea [ 48 , 49 ]. However, in our study, the geometric mean concentrations of the biomarkers of interest did not change significantly by quartiles of tea consumption ( Supplementary Materials Table S2 ), indicating that there is little chance that the observed associations in our analyses was due to other substances in tea. Finally, this is a cross-sectional study by nature, therefore we cannot assess causality.
In conclusion, our study suggests that usual intake of total flavonoids at levels seen in this population had limited impact on circulating sex hormones in premenopausal women. Intake of anthocyanins and flavan-3-ols was associated with lower levels of DHEA and DHEAS, while it is important that the subclasses of flavonoid be considered separately as they may affect sex hormones differently.