Background:

Differential associations between ovarian cancer risk factors and estrogen receptor-α (ERα) ovarian tumor expression have been noted; however, no research has assessed estrogen receptor-β (ERβ) expression. Thus, in exploratory analyses, we assessed the association of several factors with ovarian cancer risk by ERβ tumor status.

Methods:

We conducted a nested case–control study within the prospective Nurses' Health Study cohorts (NHS/NHSII), with exposures collected through biennial questionnaires. Paraffin-embedded tumor blocks were requested for cases diagnosed from 1976 to 2006 (NHS) and 1989 to 2005 (NHSII) and tissue microarrays were stained for nuclear ERβ (ERβ-nuc) and cytoplasmic ERβ (ERβ-cyto), with any staining considered positive (+). We obtained odds ratios (OR) and 95% confidence intervals (CI) using multivariate polytomous logistic regression.

Results:

We included 245 cases [43% ERβ-cyto (+) and 71% ERβ-nuc (+)] and 1,050 matched controls. An inverse association was observed between parity and risk of ERβ-nuc (+) (OR, parous vs. nulliparous: 0.46; 95% CI, 0.26–0.81), but not ERβ-nuc (–) tumors (OR, parous vs. nulliparous: 1.51; 95% CI, 0.45–5.04; Pheterogeneity = 0.04). Conversely, parity was inversely associated with ERβ-cyto (–) tumors (OR, parous vs. nulliparous: 0.42; 95% CI, 0.23–0.78), but was not associated with ERβ-cyto (+) tumors (OR, parous vs. nulliparous: 1.08; 95% CI, 0.45–2.63; Pheterogeneity = 0.05). Associations for other exposures, including hormone therapy, did not differ by ERβ-nuc or ERβ-cyto status.

Conclusions:

Our results suggest that parity may influence ovarian cancer risk, in part, through alterations in ERβ localization within tumor cells.

Impact:

Alterations in ERβ expression and localization appear to be important for ovarian cancer etiology. Future research should confirm our results and assess potential biologic mechanisms for the observed associations.

This article is featured in Highlights of This Issue, p. 2105

Extensive epidemiologic research has demonstrated that hormonal and reproductive exposures are associated with risk of ovarian cancer. However, the mechanisms by which these factors influence ovarian carcinogenesis and resulting tumor development remain poorly understood. Epithelial ovarian cancers are heterogeneous in their natural history, morphology, and gene and protein expression (1–3). These differences may be crucial to understand the etiology, prognosis, and effective treatment of ovarian cancer. It is notable that nearly all established hormonal and reproductive exposures demonstrated differential associations with ovarian cancer histologic subtypes (4); however, limited research is available on whether risk factor profiles also differ by hormonal receptor expression patterns in ovarian tumors.

While previous research has shown differential associations between ovarian cancer risk factors and estrogen receptor-α (ERα) and progesterone receptor (PR) status (5, 6), no previous studies have examined associations of ovarian cancer risk factors with estrogen receptor-β (ERβ) expression. Several studies suggest that ERβ acts as a tumor suppressor in ovarian tissue and ERβ expression may be lost with malignant transformation; however, the loss of ERβ may be dependent on cell localization (nuclear vs. cytoplasmic; refs. 7–10). De Stefano and colleagues noted that ERβ staining was more likely to be localized in the nucleus in normal ovarian tissue and conversely, was more likely to be localized in the cytoplasm in ovarian cancer cells (11). In addition, higher levels of ERβ protein expression have been associated with improved progression-free and overall survival, as well as decreased probability of lymph node metastasis in serous tumors, although the cellular localization of ERβ was not described in these studies (12, 13). Conversely, two recent studies reported less favorable prognosis for ovarian tumors expressing ERβ; however, ERβ localization differed in these two studies, with one study noting worse prognosis with higher ERβ nuclear staining and the other with higher ERβ cytoplasmic staining (14, 15). These conflicting results, potentially due to differences in ERβ cell localization (nuclear vs. cytoplasmic), as well as the limited data on associations between ovarian cancer risk factors and hormone receptor subtypes, suggest that additional research on ERβ expression, both nuclear and cytoplasmic, in ovarian tumors is needed.

Therefore, in exploratory analyses, we examined associations between known risk factors for ovarian cancer and risk of ovarian cancer by expression of nuclear and cytoplasmic ERβ in paraffin-embedded tumor tissue from ovarian cancers diagnosed among women in the Nurses' Health Studies (NHS and NHSII).

Study population

In 1976, 121,700 U.S. female nurses, ages 30–55, were enrolled into the NHS after completing a baseline questionnaire collecting data on disease diagnoses and exposures. The NHSII began in 1989 when 116,429 U.S. female nurses, ages 25–42, completed baseline questionnaires. Information on exposures and disease diagnoses has been updated through mailed biennial questionnaires in NHS and NHSII, with a follow-up rate of 85%–90% for both cohorts. We identified incident cases of ovarian cancer through self-reports on the biennial questionnaires, linkage to the National Death Index (16), or report from family members. The study protocol was approved by the Institutional Review Boards of the Brigham and Women's Hospital (Boston, MA) and Harvard T.H. Chan School of Public Health (Boston, MA), and those of participating registries as required.

Ovarian tumor block collection

Cases diagnosed from 1976 to 2006 for NHS and 1989 to 2005 for NHSII, and for whom we obtained a pathology report of the primary surgery, were included in this study. Case pathology reports were reviewed by a gynecologic pathologist, blinded to exposure status, to confirm the epithelial ovarian cancer diagnosis and classify the tumors by histologic type (serous/poorly differentiated, mucinous, endometrioid, clear cell, other), behavior (invasive and borderline), and stage (I/II and III/IV). For confirmed cases, we requested representative paraffin-embedded tissue blocks of the ovarian tumors. Through 2006, there were 1,083 confirmed NHS cases, of which 217 were included on tissue microarrays (TMA). In the NHSII, 46 of the 201 confirmed cases through 2005 were included on TMAs. Details of case exclusions have been described previously (5). Briefly, cases were excluded because of incomplete information to request tumor tissue (NHS, n = 290 and NHSII, n = 77), inability to obtain tumor blocks (NHS, n = 479 and NHSII, n = 65), or no primary tumor blocks available (NHS, n = 97 and NHSII, n = 13). For these analyses, we additionally excluded cases on the TMAs with no tumor in the core slice (NHS, n = 2), noninterpretable stains (NHS, n = 13 and NHSII, n = 2), or who withdrew consent (NHS, n = 1). Separately within each cohort, up to four controls (n = 1,050) were matched to each case on year of birth. The selected controls did not have a prior bilateral oophorectomy, menopause due to irradiation, or prior diagnosis of cancer (except nonmelanoma skin cancer), and were alive at the time of case diagnosis. Matched controls of excluded cases were included in the analyses.

Assays

Details of the TMA construction and IHC process have been described previously (6). Briefly, the primary tumor block for each case was selected by the study pathologist (J.L. Hecht), who circled the tumor on the slide, confirmed histology and behavior, assessed grade (1–3), and sent the slide/block to the Dana-Farber/Harvard Cancer Center Specialized Histopathology Services Core for TMA construction. Three cores per case were transferred to a recipient paraffin-embedded block and sections were cut to create array slides. Using IHC, the slides were processed and stained within 2 weeks of creation using the primary antibody ERβ (anti–estrogen receptor beta 1, clone: PPG5/10, Thermo Fisher Scientific Pierce). For nuclear ERβ (ERβ-nuc), the pathologist scored the number of reactive versus total cells (0%, 1%–10%, 11%–50%, 51%–90%, and >90%) for each of the three cores per case independently. For cytoplasmic ERβ (ERβ-cyto), the pathologist scored the cores as having any cytoplasmic versus no cytoplasmic ERβ staining. ERα (primary antibody: rabbit monoclonal; clone SP1; Neomarkers; dilution 1:40) and PR (primary antibody: mouse monoclonal; clone PgR 636; Dako; dilution 1:150) staining were conducted similarly to ERβ-nuc with the pathologist scoring the number of reactive versus total cells (0%, 1%–10%, 11%–50%, 51%–90%, and >90%) for each of the three cores per case independently. If tissue was missing from the slide or only a few cell clusters (<20 cells) were present, then the core was designated as noninterpretable.

Assessment of exposure and covariate information

For each case and their matched controls, exposure and covariate information were obtained from the biennial questionnaire cycle prior to case diagnosis. We assessed well-established ovarian cancer risk factors including age (years), oral contraceptive (OC) use (ever vs. never; never, >0–<1 year, 1–5 years, and >5 years; continuous), tubal ligation (ever vs. never), parity (parous vs. nulliparous; nulliparous, 1–2 children, and 3+ children; continuous), menopausal status (premenopausal and postmenopausal), hormone therapy (HT; ever vs. never), and family history of breast/ovarian cancer (yes vs. no). HT use was further divided into estrogen-only HT use (ever vs. never) and estrogen plus progesterone HT use (ever vs. never), for which these categories were not mutually exclusive. In analyses on parity, pregnancies lasting 6 months or longer were included. Ovulatory years were calculated as age at natural menopause (current age for premenopausal women) minus age at menarche, 1 year for each pregnancy, and OC duration, and was categorized into quartiles. In addition, we assessed age at menarche (<13 years vs. ≥13 years), age at natural menopause (<50 years vs. ≥50 years), height (continuous), body mass index (BMI; <25, 25–<30, and ≥30 kg/m2; continuous), and breastfeeding among parous women (ever vs. never; 0–<1 month, 1–6 months, 7–11 months, and ≥12 months; continuous). Age at natural menopause excluded women reporting a hysterectomy before menopause.

Statistical analysis

We pooled the NHS (201 cases and 866 controls) and NHSII (44 cases and 184 controls) data for all analyses. ERβ-nuc and ERβ-cyto were considered positive (+) if ≥1% of cells stained positive in any core and negative (–) if 0% of cells stained positive. Differential expression of ERβ-nuc and ERβ-cyto staining was assessed using a χ2 test. In addition, the χ2 test was used to assess the distribution of ERβ-nuc and ERβ-cyto staining positivity by histologic subtype, morphology, stage, and ERα, PR, and joint ERα/PR staining positivity. χ2 test for trend was used to assess the distribution of ERβ-nuc and ERβ-cyto staining positivity by grade. Spearman correlation coefficients were calculated between the maximum ERβ-nuc, ERα, and PR scores, which were measured as the percentage of cells staining positive (0%, 1%–10%, 11%–50%, 51%–90%, and >90%).

We used polytomous logistic regression (PLR) to assess the relationship between ovarian cancer risk factors and ovarian cancer risk by ERβ-nuc and ERβ-cyto staining positivity using a three category outcome [ERβ-nuc (+), ERβ-nuc (–), and controls, or ERβ-cyto (+), ERβ-cyto (–), and controls]. Women with missing exposure information were excluded for the specific analysis with missing data. For each exposure, we used a likelihood ratio (LR) test to compare an unconstrained model (i.e., the association between the exposure and ovarian cancer risk by hormone receptor expression was allowed to vary) with a constrained model (i.e., the association was not allowed to vary) to assess the heterogeneity in the ORs by ERβ-nuc and ERβ-cyto expression versus controls (17). Analyses were adjusted for age at diagnosis (years) and cohort (NHS/NHSII), with estimates allowed to vary by receptor status, as well as family history of breast and/or ovarian cancer (yes/no), OC use (continuous in months), number of pregnancies, and menopausal/HT status (premenopausal/postmenopausal + no HT use/postmenopausal + HT use/missing), for which estimates were constrained. In models of estrogen-only HT use, we further adjusted for estrogen plus progesterone HT use (ever/never) and other HT use (any HT use that did not include estrogen-only or estrogen plus progesterone; ever/never), for which estimates were constrained. In models of estrogen plus progesterone HT use, we further adjusted for estrogen-only HT use (ever/never) and other HT use (ever/never), for which estimates were constrained. Women missing OC duration (n = 25) or number of pregnancies (n = 19) were set to the cohort-specific median. In secondary analyses, we assessed joint ERβ-nuc and ERβ-cyto staining positivity using PLR with a five category outcome [ERβ-nuc (+)/ERβ-cyto (+), ERβ-nuc (+)/ERβ-cyto (–), ERβ-nuc (–)/ERβ-cyto (+), ERβ-nuc (–)/ERβ-cyto (–), and controls] adjusting for the same covariates as above.

In sensitivity analyses, we used case–case analyses to explore the influence of histology on the differential risk factor associations by ERβ-cyto status, as ERβ-cyto staining differed by histology. “Cases” were defined as receptor-positive tumors, whereas “controls” were defined as receptor-negative tumors. We compared an unconditional logistic regression model without adjustment for histology with a model with adjustment for histology using an LR test, with the exception of breastfeeding and OC use among ever OC users due to the small sample size of nonserous tumors. We adjusted for the same covariates as above in addition to histology (serous, mucinous, endometrioid, clear cell, and other). In addition, as ERβ-cyto staining differed by ERα expression, we assessed joint ERα and ERβ-cyto staining; because of the small sample, we excluded women with ERα (–) and ERβ-cyto (+) staining (n = 15).

We performed analyses using SAS version 9.4 (SAS Institute Inc.) and Stata version 15.1 (Stata Corp LP). All P values were two-sided and considered statistically significant if less than 0.05.

Among the 245 ovarian cancer cases included in the analyses, 43% stained positive for ERβ-cyto receptors and 71% stained positive for ERβ-nuc receptors (Table 1). The distribution of histology differed between ERβ-cyto (+) and ERβ-cyto (–) tumors, with ERβ-cyto (+) tumors more likely to be serous (75% vs. 59%) and less likely to be mucinous (1% vs. 8%) or clear cell (5% vs. 11%; P = 0.008). In addition, ERβ-cyto (+) tumors were more likely to be invasive (P = 0.02) and higher stage (P = 0.09) compared with ERβ-cyto (–) tumors. ERβ-nuc (–) tumors were more likely to be invasive compared with ERβ-nuc (+) tumors (P = 0.01). No other associations with tumor characteristics were observed for ERβ-nuc staining (P ≥ 0.60).

Table 1.

Distribution of ovarian tumor characteristics according to cytoplasmic and nuclear ERβ staining expression.

ERβ-cytoaERβ-nuca
All cases(+)(–)Pb(+)(–)Pb
Total, N (%) 245 (100%) 105 (43%) 140 (57%)  175 (71%) 70 (29%)  
Histology, n (%) 
 Serous 161 (66%) 79 (75%) 82 (59%) 0.008 115 (66%) 46 (66%) 0.62 
 Mucinous 12 (5%) 1 (1%) 11 (8%)  9 (5%) 3 (4%)  
 Endometrioid 41 (17%) 18 (17%) 23 (16%)  26 (15%) 15 (21%)  
 Clear cell 20 (8%) 5 (5%) 15 (11%)  16 (9%) 4 (6%)  
 Other 11 (4%) 2 (2%) 9 (6%)  9 (5%) 2 (3%)  
Morphology, n (%) 
 Borderline 41 (17%) 11 (10%) 30 (21%) 0.02 36 (21%) 5 (7%) 0.01 
 Invasive 204 (83%) 94 (90%) 110 (79%)  139 (79%) 65 (93%)  
Grade, n (%)c 
 1 22 (11%) 10 (11%) 12 (11%) 0.74 14 (10%) 8 (12%) 0.62 
 2 9 (4%) 3 (3%) 6 (6%)  6 (4%) 3 (5%)  
 3 170 (85%) 80 (86%) 90 (83%)  117 (86%) 53 (83%)  
Stage, n (%) 
 I/II 99 (40%) 36 (34%) 63 (45%) 0.09 72 (41%) 27 (39%) 0.71 
 III/IV 146 (60%) 69 (66%) 77 (55%)  103 (59%) 43 (61%)  
ERβ-cyto 
 (+) 105 (43%) N/A N/A  61 (35%) 44 (63%) <0.0001 
 (–) 140 (57%)    114 (65%) 26 (37%)  
ERβ-nuc 
 (+) 175 (71%) 61 (58%) 114 (81%) <0.0001 N/A N/A  
 (–) 70 (29%) 44 (42%) 26 (19%)     
ERα 
 (+) 185 (76%) 90 (86%) 95 (68%) 0.001 129 (74%) 56 (80%) 0.30 
 (–) 60 (24%) 15 (14%) 45 (32%)  46 (26%) 14 (20%)  
PR 
 (+) 119 (49%) 60 (57%) 59 (42%) 0.02 80 (46%) 39 (56%) 0.16 
 (–) 126 (51%) 45 (43%) 81 (58%)  95 (54%) 31 (44%)  
ERα/PR joint statusd 
 ERα(+)/PR(+) 110 (47%) 60 (57%) 50 (38%) 0.007 74 (44%) 36 (54%) 0.32 
 ERα(+)/PR(–) 75 (32%) 30 (29%) 45 (34%)  55 (33%) 20 (30%)  
 ERα(–)/PR(–) 51 (21%) 15 (14%) 36 (28%)  40 (24%) 11 (16%)  
ERβ-cytoaERβ-nuca
All cases(+)(–)Pb(+)(–)Pb
Total, N (%) 245 (100%) 105 (43%) 140 (57%)  175 (71%) 70 (29%)  
Histology, n (%) 
 Serous 161 (66%) 79 (75%) 82 (59%) 0.008 115 (66%) 46 (66%) 0.62 
 Mucinous 12 (5%) 1 (1%) 11 (8%)  9 (5%) 3 (4%)  
 Endometrioid 41 (17%) 18 (17%) 23 (16%)  26 (15%) 15 (21%)  
 Clear cell 20 (8%) 5 (5%) 15 (11%)  16 (9%) 4 (6%)  
 Other 11 (4%) 2 (2%) 9 (6%)  9 (5%) 2 (3%)  
Morphology, n (%) 
 Borderline 41 (17%) 11 (10%) 30 (21%) 0.02 36 (21%) 5 (7%) 0.01 
 Invasive 204 (83%) 94 (90%) 110 (79%)  139 (79%) 65 (93%)  
Grade, n (%)c 
 1 22 (11%) 10 (11%) 12 (11%) 0.74 14 (10%) 8 (12%) 0.62 
 2 9 (4%) 3 (3%) 6 (6%)  6 (4%) 3 (5%)  
 3 170 (85%) 80 (86%) 90 (83%)  117 (86%) 53 (83%)  
Stage, n (%) 
 I/II 99 (40%) 36 (34%) 63 (45%) 0.09 72 (41%) 27 (39%) 0.71 
 III/IV 146 (60%) 69 (66%) 77 (55%)  103 (59%) 43 (61%)  
ERβ-cyto 
 (+) 105 (43%) N/A N/A  61 (35%) 44 (63%) <0.0001 
 (–) 140 (57%)    114 (65%) 26 (37%)  
ERβ-nuc 
 (+) 175 (71%) 61 (58%) 114 (81%) <0.0001 N/A N/A  
 (–) 70 (29%) 44 (42%) 26 (19%)     
ERα 
 (+) 185 (76%) 90 (86%) 95 (68%) 0.001 129 (74%) 56 (80%) 0.30 
 (–) 60 (24%) 15 (14%) 45 (32%)  46 (26%) 14 (20%)  
PR 
 (+) 119 (49%) 60 (57%) 59 (42%) 0.02 80 (46%) 39 (56%) 0.16 
 (–) 126 (51%) 45 (43%) 81 (58%)  95 (54%) 31 (44%)  
ERα/PR joint statusd 
 ERα(+)/PR(+) 110 (47%) 60 (57%) 50 (38%) 0.007 74 (44%) 36 (54%) 0.32 
 ERα(+)/PR(–) 75 (32%) 30 (29%) 45 (34%)  55 (33%) 20 (30%)  
 ERα(–)/PR(–) 51 (21%) 15 (14%) 36 (28%)  40 (24%) 11 (16%)  

aOvarian tumors were classified as ERβ-cyto (+) and ERβ-nuc (+) if 1% or greater cells stained positive.

bP values were calculated using χ2 test, and χ2 test for trend (grade only) comparing ERβ-cyto (+) versus ERβ-cyto () and ERβ-nuc (+) versus ERβ-nuc ().

cA total of 201 cases had data on grade.

dExcluding nine cases with ERα(–)/PR(+) ovarian tumors.

Tumors categorized as ERβ-nuc (–) were more likely to be categorized as ERβ-cyto (+) compared with ERβ-cyto (–), with 63% of ERβ-nuc (–) tumors staining positive for ERβ-cyto and 65% of ERβ-nuc (+) tumors staining negative for ERβ-cyto (P < 0.0001; Table 1). ERβ-cyto (+) tumors were more likely to stain positive for ERα compared with ERβ-cyto (–) tumors (86% vs. 68% stained positive for ERα, respectively; P = 0.001). ERβ-cyto staining was also significantly associated with PR (P = 0.02) and joint ERα/PR staining status (P = 0.007). ERβ-nuc staining was not associated with ERα or PR staining in the dichotomous classification (P > 0.15) or in correlations between the percentage of cells staining positive (Spearman correlation for ERβ-nuc and ERα = 0.03, P = 0.69; Spearman correlation for ERβ-nuc and PR = −0.05, P = 0.47).

ERβ-cyto staining

An inverse association was observed between parity and risk of ERβ-cyto (–) tumors [OR, 0.42; 95% confidence interval (CI), 0.23–0.78 for parous vs. nulliparous), but no association was observed for ERβ-cyto (+) tumors (OR, 1.08; 95% CI, 0.45–2.63 for parous vs. nulliparous; Pheterogeneity = 0.05; Table 2). Interestingly, among parous women, there was a 16% (OR, 0.84; 95% CI, 0.71–0.99) decreased risk of ERβ-cyto (+) tumors per child, but no association with ERβ-cyto (–) tumors (OR, 1.01; 95% CI, 0.89–1.14 per child; Pheterogeneity = 0.07). In addition, there was a suggestion of a reduced risk of ERβ-cyto (–) tumors with increasing BMI (OR, 0.96; 95% CI, 0.92–1.00 per kg/m2), but not for ERβ-cyto (+) tumors (OR, 1.00; 95% CI, 0.96–1.04; Pheterogeneity = 0.09). The other ovarian cancer risk factors examined, including OC use, tubal ligation, and menopausal status, did not appear to differ by ERβ-cyto staining positivity (Pheterogeneity > 0.10).

Table 2.

Association between ovarian cancer risk factors and risk of ovarian cancer by ERβ-cyto staininga,b.

ERβ-cyto (+)ERβ-cyto (–)
ControlsCasesOR (95% CI)CasesOR (95% CI)Pheterogeneity
Height 
 Per 5 cm 1,049 105 1.12 (0.95–1.32) 140 1.19 (1.03–1.38) 0.54 
BMI (kg/m2
 <25 481 46 1.00 (Ref) 80 1.00 (Ref) 0.10 
 25–<30 329 27 0.82 (0.50–1.36) 31 0.55 (0.35–0.85)  
 ≥30 184 21 1.19 (0.69–2.06) 17 0.57 (0.32–0.98)  
 per kg/m2 994 94 1.00 (0.96–1.04) 128 0.96 (0.92–1.00) 0.09 
Family history of breast or ovarian cancer 
 No 896 85 1.00 (Ref) 118 1.00 (Ref) 0.49 
 Yes 154 20 1.40 (0.83–2.36) 22 1.11 (0.68–1.81)  
Parity 
 Nulliparous 59 1.00 (Ref) 16 1.00 (Ref) 0.05 
 Parous 978 98 1.08 (0.45–2.63) 119 0.42 (0.23–0.78)  
 Nulliparous 59 1.00 (Ref) 16 1.00 (Ref) 0.16 
 1–2 children 400 51 1.33 (0.54–3.27) 60 0.52 (0.27–0.98)  
 3+ children 578 47 0.85 (0.34–2.12) 59 0.34 (0.18–0.65)  
 Per childc 1,037 104 0.88 (0.76–1.01) 135 0.94 (0.83–1.06) 0.45 
 Per childd 978 98 0.84 (0.71–0.99) 119 1.01 (0.89–1.14) 0.07 
OC use 
 Never 487 47 1.00 (Ref) 74 1.00 (Ref) 0.33 
 Ever 544 54 0.90 (0.58–1.42) 64 0.69 (0.46–1.01)  
 Never 487 47 1.00 (Ref) 74 1.00 (Ref) 0.80 
 <1 year 120 10 0.79 (0.38–1.63) 13 0.63 (0.33–1.18)  
 1–5 years 217 24 1.03 (0.59–1.79) 29 0.80 (0.49–1.30)  
 >5 years 207 20 0.86 (0.48–1.54) 22 0.61 (0.36–1.04)  
 Per yearc 1,031 101 1.00 (0.98–1.03) 138 1.01 (0.99–1.03) 0.56 
 Per yeare 544 54 1.00 (0.98–1.03) 64 1.02 (1.00–1.04) 0.37 
Breast feedingd 
 Never 397 48 1.00 (Ref) 56 1.00 (Ref) 0.96 
 Ever 418 35 0.77 (0.48–1.24) 45 0.76 (0.50–1.17)  
 0–1 months 397 48 1.00 (Ref) 56 1.00 (Ref) 0.11 
 1–6 months 203 23 1.04 (0.60–1.79) 18 0.60 (0.34–1.06)  
 7–11 months 73 0.39 (0.12–1.29) 0.84 (0.40–1.80)  
 12+ months 142 0.59 (0.28–1.26) 18 0.98 (0.55–1.75)  
Ovulatory yearsf 
 Quartile 1 224 15 1.00 (Ref) 13 1.00 (Ref) 0.59 
 Quartile 2 224 22 1.46 (0.71–3.04) 34 2.70 (1.33–5.50)  
 Quartile 3 241 19 1.23 (0.56–2.74) 29 2.24 (1.04–4.83)  
 Quartile 4 208 20 1.38 (0.60–3.14) 27 2.21 (1.00–4.91)  
 Per year 897 76 1.01 (0.96–1.06) 103 1.02 (0.98–1.07) 0.55 
Tubal ligation 
 Never 837 92 1.00 (Ref) 123 1.00 (Ref) 0.94 
 Ever 213 13 0.54 (0.30–1.00) 17 0.53 (0.31–0.90)  
Age at menarche 
 <13 years 538 48 1.00 (Ref) 71 1.00 (Ref) 0.42 
 ≥13 years 512 57 1.16 (0.77–1.74) 69 0.94 (0.66–1.35)  
Menopause 
 Premenopausal 256 27 1.00 (Ref) 28 1.00 (Ref) 0.12 
 Postmenopausal 756 73 0.78 (0.38–1.60) 106 1.61 (0.86–3.03)  
Age at menopauseg 
 <50 years 239 16 1.00 (Ref) 28 1.00 (Ref) 0.63 
 ≥50 years 404 34 1.32 (0.71–2.46) 49 1.09 (0.66–1.80)  
HT useg 
 Never 304 17 1.00 (Ref) 27 1.00 (Ref) 0.70 
 Ever 404 51 2.32 (1.31–4.13) 68 2.02 (1.26–3.26)  
Estrogen-only HT useg,h 
 Never 495 35 1.00 (Ref) 52 1.00 (Ref) 0.85 
 Ever 171 29 2.50 (1.46–4.30) 38 2.35 (1.47–3.77)  
Estrogen + progesterone HT useg,i 
 Never 459 39 1.00 (Ref) 63 1.00 (Ref) 0.31 
 Ever 215 23 1.36 (0.78–2.35) 27 0.95 (0.58–1.54)  
ERβ-cyto (+)ERβ-cyto (–)
ControlsCasesOR (95% CI)CasesOR (95% CI)Pheterogeneity
Height 
 Per 5 cm 1,049 105 1.12 (0.95–1.32) 140 1.19 (1.03–1.38) 0.54 
BMI (kg/m2
 <25 481 46 1.00 (Ref) 80 1.00 (Ref) 0.10 
 25–<30 329 27 0.82 (0.50–1.36) 31 0.55 (0.35–0.85)  
 ≥30 184 21 1.19 (0.69–2.06) 17 0.57 (0.32–0.98)  
 per kg/m2 994 94 1.00 (0.96–1.04) 128 0.96 (0.92–1.00) 0.09 
Family history of breast or ovarian cancer 
 No 896 85 1.00 (Ref) 118 1.00 (Ref) 0.49 
 Yes 154 20 1.40 (0.83–2.36) 22 1.11 (0.68–1.81)  
Parity 
 Nulliparous 59 1.00 (Ref) 16 1.00 (Ref) 0.05 
 Parous 978 98 1.08 (0.45–2.63) 119 0.42 (0.23–0.78)  
 Nulliparous 59 1.00 (Ref) 16 1.00 (Ref) 0.16 
 1–2 children 400 51 1.33 (0.54–3.27) 60 0.52 (0.27–0.98)  
 3+ children 578 47 0.85 (0.34–2.12) 59 0.34 (0.18–0.65)  
 Per childc 1,037 104 0.88 (0.76–1.01) 135 0.94 (0.83–1.06) 0.45 
 Per childd 978 98 0.84 (0.71–0.99) 119 1.01 (0.89–1.14) 0.07 
OC use 
 Never 487 47 1.00 (Ref) 74 1.00 (Ref) 0.33 
 Ever 544 54 0.90 (0.58–1.42) 64 0.69 (0.46–1.01)  
 Never 487 47 1.00 (Ref) 74 1.00 (Ref) 0.80 
 <1 year 120 10 0.79 (0.38–1.63) 13 0.63 (0.33–1.18)  
 1–5 years 217 24 1.03 (0.59–1.79) 29 0.80 (0.49–1.30)  
 >5 years 207 20 0.86 (0.48–1.54) 22 0.61 (0.36–1.04)  
 Per yearc 1,031 101 1.00 (0.98–1.03) 138 1.01 (0.99–1.03) 0.56 
 Per yeare 544 54 1.00 (0.98–1.03) 64 1.02 (1.00–1.04) 0.37 
Breast feedingd 
 Never 397 48 1.00 (Ref) 56 1.00 (Ref) 0.96 
 Ever 418 35 0.77 (0.48–1.24) 45 0.76 (0.50–1.17)  
 0–1 months 397 48 1.00 (Ref) 56 1.00 (Ref) 0.11 
 1–6 months 203 23 1.04 (0.60–1.79) 18 0.60 (0.34–1.06)  
 7–11 months 73 0.39 (0.12–1.29) 0.84 (0.40–1.80)  
 12+ months 142 0.59 (0.28–1.26) 18 0.98 (0.55–1.75)  
Ovulatory yearsf 
 Quartile 1 224 15 1.00 (Ref) 13 1.00 (Ref) 0.59 
 Quartile 2 224 22 1.46 (0.71–3.04) 34 2.70 (1.33–5.50)  
 Quartile 3 241 19 1.23 (0.56–2.74) 29 2.24 (1.04–4.83)  
 Quartile 4 208 20 1.38 (0.60–3.14) 27 2.21 (1.00–4.91)  
 Per year 897 76 1.01 (0.96–1.06) 103 1.02 (0.98–1.07) 0.55 
Tubal ligation 
 Never 837 92 1.00 (Ref) 123 1.00 (Ref) 0.94 
 Ever 213 13 0.54 (0.30–1.00) 17 0.53 (0.31–0.90)  
Age at menarche 
 <13 years 538 48 1.00 (Ref) 71 1.00 (Ref) 0.42 
 ≥13 years 512 57 1.16 (0.77–1.74) 69 0.94 (0.66–1.35)  
Menopause 
 Premenopausal 256 27 1.00 (Ref) 28 1.00 (Ref) 0.12 
 Postmenopausal 756 73 0.78 (0.38–1.60) 106 1.61 (0.86–3.03)  
Age at menopauseg 
 <50 years 239 16 1.00 (Ref) 28 1.00 (Ref) 0.63 
 ≥50 years 404 34 1.32 (0.71–2.46) 49 1.09 (0.66–1.80)  
HT useg 
 Never 304 17 1.00 (Ref) 27 1.00 (Ref) 0.70 
 Ever 404 51 2.32 (1.31–4.13) 68 2.02 (1.26–3.26)  
Estrogen-only HT useg,h 
 Never 495 35 1.00 (Ref) 52 1.00 (Ref) 0.85 
 Ever 171 29 2.50 (1.46–4.30) 38 2.35 (1.47–3.77)  
Estrogen + progesterone HT useg,i 
 Never 459 39 1.00 (Ref) 63 1.00 (Ref) 0.31 
 Ever 215 23 1.36 (0.78–2.35) 27 0.95 (0.58–1.54)  

Abbreviation: Ref, reference.

aAnalyses were adjusted for cohort (NHS/NHSII), age at diagnosis (per year), duration of OC use (continuous in months), family history of breast/ovarian cancer (yes/no), menopausal and HT status (premenopausal/postmenopausal + ever HT/postmenopausal + never HT/missing HT), and number of pregnancies (continuous).

bTotal N does not add up to 1,295 due to missingness in the exposure (height missing = 1; BMI missing = 79; parity missing = 19; OC missing = 25; breastfeeding missing = 212; ovulatory years missing = 219; age at menopause missing = 165; menopause missing = 49; HT use missing = 64; estrogen-only HT use missing = 115; and estrogen + progesterone HT use missing = 109).

cAmong all women.

dAmong parous women.

eAmong women who ever used OCs.

fOvulatory years were calculated as age at natural menopause (or age at diagnosis for premenopausal women) minus age at menarche with additional subtraction for OC use duration and 1 year for each pregnancy.

gAmong postmenopausal women.

hAdditionally adjusted for estrogen plus progesterone HT use (ever/never) and other HT use (ever/never).

iAdditionally adjusted for estrogen-only HT use (ever/never) and other HT use (ever/never).

The ERβ-cyto associations above were essentially unchanged after adjusting for histology in case–case analyses (Supplementary Table S1). Compared with nulliparous cases, parous cases were twice as likely to develop ERβ-cyto (+) tumors as opposed to ERβ-cyto (–) tumors without adjustment for histology (OR, 2.52; 95% CI, 0.88–7.25). With adjustment for histology, the risk estimate decreased slightly to 2.42 (95% CI, 0.82–7.12). In analyses of joint ERα/ERβ-cyto staining among all women, there were decreased risks per child for ERα (–)/ERβ-cyto (–) tumors (OR, 0.77; 95% CI, 0.62–0.97 per child) and ERα (+)/ERβ-cyto (+) tumors (OR, 0.87; 95% CI, 0.75–1.01 per child), but no association was observed for ERα (+)/ERβ-cyto (–) tumors (OR, 1.01; 95% CI, 0.88–1.16 per child; Pheterogeneity = 0.08; Supplementary Table S2). Similar associations between number of children and ERα/ERβ-cyto status of the tumor were observed among parous women (Pheterogeneity = 0.03). The association between decreased risk of ERβ-cyto (–) tumors and increasing BMI was fairly similar between ERα (+)/ERβ-cyto (–) tumors (OR, 0.95; 95% CI, 0.91–1.00) and ERα (–)/ERβ-cyto (–) tumors (OR, 0.97; 95% CI, 0.91–1.03), but no association was observed for ERα (+)/ERβ-cyto (+) tumors (OR, 1.01; 95% CI, 0.97–1.05; Pheterogeneity = 0.18).

ERβ-nuc staining

A decreased risk of ERβ-nuc (+) tumors was observed among parous women compared with nulliparous women (OR, 0.46; 95% CI, 0.26–0.81), whereas no association was observed for ERβ-nuc (–) tumors (OR, 1.51; 95% CI, 0.45–5.04; Pheterogeneity = 0.04; Table 3). The other ovarian cancer risk factors examined, including OC use, breastfeeding, BMI, or ovulatory years, did not appear to differ by ERβ-nuc staining positivity (Pheterogeneity ≥ 0.09).

Table 3.

Association between ovarian cancer risk factors and risk of ovarian cancer by ERβ-nuc staininga,b.

ERβ-nuc (+)ERβ-nuc (–)
ControlsCasesOR (95% CI)CasesOR (95% CI)Pheterogeneity
Height 
 Per 5 cm 1,049 175 1.20 (1.05–1.38) 70 1.06 (0.86–1.30) 0.28 
BMI 
 <25 481 92 1.00 (Ref) 34 1.00 (Ref) 0.52 
 25–<30 329 42 0.65 (0.44–0.96) 16 0.65 (0.35–1.21)  
 ≥30 184 24 0.69 (0.43–1.13) 14 1.07 (0.56–2.05)  
 per kg/m2 994 158 0.97 (0.93–1.00) 64 1.01 (0.96–1.06) 0.13 
Family history of breast or ovarian cancer 
 No 896 144 1.00 (Ref) 59 1.00 (Ref) 0.66 
 Yes 154 31 1.29 (0.84–1.98) 11 1.09 (0.56–2.13)  
Parity 
 Nulliparous 59 19 1.00 (Ref) 1.00 (Ref) 0.04 
 Parous 978 150 0.46 (0.26–0.81) 67 1.51 (0.45–5.04)  
 Nulliparous 59 19 1.00 (Ref) 1.00 (Ref) 0.13 
 1–2 children 400 76 0.56 (0.31–1.01) 35 1.86 (0.55–6.31)  
 3+ children 578 74 0.37 (0.20–0.67) 32 1.16 (0.34–4.03)  
 Per childc 1,037 169 0.89 (0.80–1.00) 70 0.95 (0.82–1.12) 0.48 
 Per childd 978 150 0.94 (0.83–1.06) 67 0.92 (0.77–1.11) 0.84 
OC use 
 Never 487 85 1.00 (Ref) 36 1.00 (Ref) 0.31 
 Ever 544 89 0.84 (0.59–1.19) 29 0.61 (0.35–1.06)  
 Never 487 85 1.00 (Ref) 36 1.00 (Ref) 0.60 
 <1 year 120 17 0.72 (0.40–1.27) 0.61 (0.25–1.51)  
 1–5 years 217 42 1.02 (0.66–1.57) 11 0.58 (0.28–1.22)  
 >5 years 207 30 0.73 (0.46–1.18) 12 0.63 (0.31–1.30)  
 Per yearc 1,031 174 1.01 (0.99–1.03) 65 1.00 (0.97–1.03) 0.58 
 Per yeare 544 89 1.01 (0.99–1.03) 29 1.01 (0.98–1.04) 0.76 
Breast feedingd 
 Never 397 70 1.00 (Ref) 34 1.00 (Ref) 0.28 
 Ever 418 61 0.85 (0.58–1.24) 19 0.59 (0.33–1.07)  
 0–1 months 397 70 1.00 (Ref) 34 1.00 (Ref) 0.57 
 1–6 months 203 33 0.92 (0.58–1.45) 0.50 (0.22–1.11)  
 7–11 months 73 0.70 (0.33–1.49) 0.54 (0.16–1.82)  
 12+ months 142 19 0.83 (0.47–1.44) 0.76 (0.34–1.71)  
Ovulatory yearsf 
 Quartile 1 224 21 1.00 (Ref) 1.00 (Ref) 0.99 
 Quartile 2 224 41 2.05 (1.11–3.77) 15 2.03 (0.78–5.27)  
 Quartile 3 241 36 1.78 (0.91–3.48) 12 1.51 (0.54–4.23)  
 Quartile 4 208 34 1.79 (0.88–3.62) 13 1.71 (0.60–4.87)  
 Per year 897 132 1.02 (0.98–1.06) 47 1.02 (0.96–1.08) 0.98 
Tubal ligation 
 Never 837 152 1.00 (Ref) 63 1.00 (Ref) 0.57 
 Ever 213 23 0.57 (0.35–0.91) 0.44 (0.20–0.99)  
Age at menarche 
 <13 years 512 85 1.00 (Ref) 34 1.00 (Ref) 0.96 
 ≥13 years 438 90 1.04 (0.75–1.43) 36 1.02 (0.63–1.67)  
Menopause status 
 Premenopausal 256 40 1.00 (Ref) 15 1.00 (Ref) 0.91 
 Postmenopausal 756 126 1.21 (0.69–2.13) 53 1.14 (0.48–2.73)  
Age at menopauseg 
 <50 years 239 32 1.00 (Ref) 12 1.00 (Ref) 0.97 
 ≥50 years 404 60 1.18 (0.74–1.88) 23 1.16 (0.56–2.39)  
HT useg 
 Never 304 29 1.00 (Ref) 15 1.00 (Ref) 0.41 
 Ever 404 86 2.36 (1.50–3.71) 33 1.72 (0.91–3.24)  
Estrogen-only HT useg,h 
 Never 495 65 1.00 (Ref) 22 1.00 (Ref) 0.16 
 Ever 171 43 2.06 (1.50–3.71) 24 3.46 (1.85–6.50)  
Estrogen + progesterone HT useg,i 
 Never 459 68 1.00 (Ref) 34 1.00 (Ref) 0.09 
 Ever 215 40 1.31 (0.85–2.02) 10 0.67 (0.32–1.39)  
ERβ-nuc (+)ERβ-nuc (–)
ControlsCasesOR (95% CI)CasesOR (95% CI)Pheterogeneity
Height 
 Per 5 cm 1,049 175 1.20 (1.05–1.38) 70 1.06 (0.86–1.30) 0.28 
BMI 
 <25 481 92 1.00 (Ref) 34 1.00 (Ref) 0.52 
 25–<30 329 42 0.65 (0.44–0.96) 16 0.65 (0.35–1.21)  
 ≥30 184 24 0.69 (0.43–1.13) 14 1.07 (0.56–2.05)  
 per kg/m2 994 158 0.97 (0.93–1.00) 64 1.01 (0.96–1.06) 0.13 
Family history of breast or ovarian cancer 
 No 896 144 1.00 (Ref) 59 1.00 (Ref) 0.66 
 Yes 154 31 1.29 (0.84–1.98) 11 1.09 (0.56–2.13)  
Parity 
 Nulliparous 59 19 1.00 (Ref) 1.00 (Ref) 0.04 
 Parous 978 150 0.46 (0.26–0.81) 67 1.51 (0.45–5.04)  
 Nulliparous 59 19 1.00 (Ref) 1.00 (Ref) 0.13 
 1–2 children 400 76 0.56 (0.31–1.01) 35 1.86 (0.55–6.31)  
 3+ children 578 74 0.37 (0.20–0.67) 32 1.16 (0.34–4.03)  
 Per childc 1,037 169 0.89 (0.80–1.00) 70 0.95 (0.82–1.12) 0.48 
 Per childd 978 150 0.94 (0.83–1.06) 67 0.92 (0.77–1.11) 0.84 
OC use 
 Never 487 85 1.00 (Ref) 36 1.00 (Ref) 0.31 
 Ever 544 89 0.84 (0.59–1.19) 29 0.61 (0.35–1.06)  
 Never 487 85 1.00 (Ref) 36 1.00 (Ref) 0.60 
 <1 year 120 17 0.72 (0.40–1.27) 0.61 (0.25–1.51)  
 1–5 years 217 42 1.02 (0.66–1.57) 11 0.58 (0.28–1.22)  
 >5 years 207 30 0.73 (0.46–1.18) 12 0.63 (0.31–1.30)  
 Per yearc 1,031 174 1.01 (0.99–1.03) 65 1.00 (0.97–1.03) 0.58 
 Per yeare 544 89 1.01 (0.99–1.03) 29 1.01 (0.98–1.04) 0.76 
Breast feedingd 
 Never 397 70 1.00 (Ref) 34 1.00 (Ref) 0.28 
 Ever 418 61 0.85 (0.58–1.24) 19 0.59 (0.33–1.07)  
 0–1 months 397 70 1.00 (Ref) 34 1.00 (Ref) 0.57 
 1–6 months 203 33 0.92 (0.58–1.45) 0.50 (0.22–1.11)  
 7–11 months 73 0.70 (0.33–1.49) 0.54 (0.16–1.82)  
 12+ months 142 19 0.83 (0.47–1.44) 0.76 (0.34–1.71)  
Ovulatory yearsf 
 Quartile 1 224 21 1.00 (Ref) 1.00 (Ref) 0.99 
 Quartile 2 224 41 2.05 (1.11–3.77) 15 2.03 (0.78–5.27)  
 Quartile 3 241 36 1.78 (0.91–3.48) 12 1.51 (0.54–4.23)  
 Quartile 4 208 34 1.79 (0.88–3.62) 13 1.71 (0.60–4.87)  
 Per year 897 132 1.02 (0.98–1.06) 47 1.02 (0.96–1.08) 0.98 
Tubal ligation 
 Never 837 152 1.00 (Ref) 63 1.00 (Ref) 0.57 
 Ever 213 23 0.57 (0.35–0.91) 0.44 (0.20–0.99)  
Age at menarche 
 <13 years 512 85 1.00 (Ref) 34 1.00 (Ref) 0.96 
 ≥13 years 438 90 1.04 (0.75–1.43) 36 1.02 (0.63–1.67)  
Menopause status 
 Premenopausal 256 40 1.00 (Ref) 15 1.00 (Ref) 0.91 
 Postmenopausal 756 126 1.21 (0.69–2.13) 53 1.14 (0.48–2.73)  
Age at menopauseg 
 <50 years 239 32 1.00 (Ref) 12 1.00 (Ref) 0.97 
 ≥50 years 404 60 1.18 (0.74–1.88) 23 1.16 (0.56–2.39)  
HT useg 
 Never 304 29 1.00 (Ref) 15 1.00 (Ref) 0.41 
 Ever 404 86 2.36 (1.50–3.71) 33 1.72 (0.91–3.24)  
Estrogen-only HT useg,h 
 Never 495 65 1.00 (Ref) 22 1.00 (Ref) 0.16 
 Ever 171 43 2.06 (1.50–3.71) 24 3.46 (1.85–6.50)  
Estrogen + progesterone HT useg,i 
 Never 459 68 1.00 (Ref) 34 1.00 (Ref) 0.09 
 Ever 215 40 1.31 (0.85–2.02) 10 0.67 (0.32–1.39)  

Abbreviation: Ref, reference.

aAnalyses were adjusted for cohort (NHS/NHSII), age at diagnosis (per year), duration of OC use (continuous in months), family history of breast/ovarian cancer (yes/no), menopausal and HT status (premenopausal/postmenopausal + ever HT/postmenopausal + never HT/missing HT), and number of pregnancies (continuous).

bTotal N does not add up to 1,295 due to missingness in the exposure (height missing = 1; BMI missing = 79; parity missing = 19; OC use missing = 25; breastfeeding missing = 212; ovulatory years missing = 219; age at menopause missing = 165; menopause missing = 49; HT use missing = 64; estrogen-only HT use missing = 115; and estrogen + progesterone HT use missing = 109).

cAmong all women.

dAmong parous women.

eAmong women who ever used OCs.

fOvulatory years were calculated as age at natural menopause (or age at diagnosis for premenopausal women) minus age at menarche with additional subtraction for OC use duration and 1 year for each pregnancy.

gAmong postmenopausal women.

hAdditionally adjusted for estrogen plus progesterone HT use (ever/never) and other HT use (ever/never).

iAdditionally adjusted for estrogen-only HT use (ever/never) and other HT use (ever/never).

Joint ERβ-cyto and ERβ-nuc staining

The decreased risk for parity with ERβ-cyto (–) tumors appeared to be slightly stronger for ERβ-cyto (–)/ERβ-nuc (+) tumors (OR, 0.38; 95% CI, 0.20–0.74 for parous vs. nulliparous) compared with ERβ-cyto (–)/ERβ-nuc (–) tumors (OR, 0.69; 95% CI, 0.15–3.10 for parous vs. nulliparous; Pheterogeneity = 0.08; Supplementary Table S3). There was a suggestion of heterogeneity across hormone receptor status for duration of OC use, with a stronger decreased risk of ovarian cancer for ERβ-cyto (–)/ERβ-nuc (–) tumors (OR, 0.85; 95% CI, 0.70–1.02 per year) compared with the other tumor subtypes (ORs ranged from 0.98 to 1.02; Pheterogeneity = 0.06) and a similar association was observed among ever OC users.

This is the first study to assess the differential associations of ovarian cancer risk factors by ERβ-cyto and ERβ-nuc status of the ovarian tumors. ERβ-nuc was inversely related to ERβ-cyto. Notably, ERβ-cyto expression differed by tumor characteristics, particularly histology and morphology; however, there appeared to be a similar level of ERβ-nuc expression by tumor characteristics with the exception of morphology. Parity was associated with an inverse association for ERβ-cyto (–) as well as ERβ-nuc (+) ovarian tumors. Although the P value for heterogeneity was not significant, there was a stronger inverse association of parity with ERβ-cyto (–)/ERβ-nuc (+) tumors compared with ERβ-cyto (+)/ERβ-nuc (+) and ERβ-cyto (–)/ERβ-nuc (–) tumors. In addition, there was a suggestion of a stronger inverse association between ERβ-cyto (–) versus ERβ-cyto (+) tumors and BMI. Other risk factors, including HT use, did not demonstrate differential associations by ERβ staining.

We observed that 71% of ovarian tumors included in this study expressed nuclear staining and 43% expressed cytoplasmic staining for ERβ. Previous reports of ERβ staining have varied (53%–93% staining positive), with the majority of studies reporting that approximately 60%–80% of ovarian tumors stained positive for ERβ (10, 11, 13–15, 18), which is in-line with our observed expression of ERβ-nuc staining. Differences in ERβ staining between our study and previous studies may be due to the localization of ERβ staining (nuclear vs. cytoplasmic), differences in antibodies used to detect ERβ, and differences in techniques to determine ERβ levels (mRNA vs. IHC). Interestingly, Chan and colleagues observed in a small pilot study that mRNA levels of ERβ did not correlate with IHC staining for nuclear ERβ (18). Similar to our study, the majority of previous studies observed no significant differences of ERβ-nuc staining with histology or grade (10, 11, 13, 19). Contrary to two previous studies, we observed significant differences of ERβ-cyto staining with histology (11, 19). Interestingly, previous studies have noted a shift in ERβ staining with malignant transformation of ovarian cells (9, 11, 18). Lower ERβ-nuc staining was observed in ovarian cancer cells compared with normal ovarian cells, but higher ERβ-cyto staining was observed in ovarian cancer cells compared with normal ovarian cells. This is consistent with our finding of a strong inverse association of staining within these different cellular compartments, suggesting that malignant transformation may lead to a reduction in localization of ERβ in the nucleus and supporting examination of the ratio of ERβ-nuc to ERβ-cyto in future epidemiologic studies.

Our most notable finding was that parity was only inversely related to risk of ovarian cancer for ERβ-nuc (+) and ERβ-cyto (–) tumors. The latter finding is consistent with our prior work demonstrating that increasing number of pregnancies was only inversely associated with ERα (–) tumors, but not ERα (+); ERα was strongly positively related to ERβ-cyto in our data. Interestingly, in a mouse model, parous versus nulliparous mice had reduced nuclear ERα expression in mammary tissue, but increased ERβ-nuc expression; ERβ-cyto expression did not vary after pregnancy (20). In addition, low nuclear ERα and high ERβ-nuc expression in parous mice was associated with lower Ki67, a marker of proliferation (20, 21). Furthermore, in rats, parity was associated with reduced responsiveness to estradiol as demonstrated by lower uterine weight (22). Similar studies of the effect of parity on ERα and ERβ expression and localization in ovarian and fallopian tube tissue would provide potential insight into how parity may lower risk of ovarian cancer.

In addition, increasing BMI was associated with a decreased risk of ERβ-cyto (–) tumors, but was not associated with ERβ-cyto (+) tumors. Interestingly, genetic variation in ERβ is associated with risk of type 2 diabetes (23). Furthermore, in skin punch biopsies from postmenopausal women, higher adiposity was associated with lower expression of ERβ (24). However, these results should be interpreted with caution as most studies suggest that BMI is associated with a modest increased risk of ovarian cancer (25, 26). The other ovarian cancer risk factors did not appear to differ by ERβ staining.

As this is the first study that has assessed ovarian cancer risk factors by ERβ-nuc and -cyto status and due to our small sample size, future consortium studies with more ovarian tumor samples will be needed to confirm our findings. In addition, as ovarian tumors were only available for a subset of the ovarian cancer cases in NHS and NHSII, our study could have selection bias. To assess for selection bias, we compared the distribution of reproductive and hormonal factors and tumor characteristics between confirmed ovarian cancer cases, cases eligible for tumor block collection, cases with collected tumor blocks, and cases included on the TMA (5). Cases on the TMA were similar to all confirmed ovarian cancer cases with the exception that TMA cases were older at diagnosis, had lower stage tumors, and were slightly more likely to be parous (NHS II only), suggesting that the potential for selection bias within our study is low. As with all epidemiologic studies, there is the chance that our results may be due to residual confounding; however, we adjusted for known ovarian cancer risk factors. Our study also had important strengths. It is the first study to assess lifestyle and reproductive factors with the risk of ovarian cancer by ERβ staining, and our sample of 245 tumor blocks is large compared with other similar studies. In addition, exposure and covariate information were collected prospectively, eliminating the possibility of differential recall in cases versus controls. Finally, the use of TMAs limited variability in assay results as many cases could be stained at the same time.

Our results suggest that parity may be differentially associated with ovarian tumors based on both ERβ-cyto and ERβ-nuc expression. Specifically, parity may influence ovarian cancer risk, in part, through alterations in ERβ expression and localization within tumor cells. These results should be confirmed in future studies with a larger sample size. In addition, future studies should assess the biologic mechanisms that may underlie the observed associations.

A.L. Shafrir reports grants from the NIH and the Marriott Family Foundations outside the submitted work. M.S. Rice reports employment with Sanofi (all work on this article was completed prior to employment at Sanofi). J.L. Hecht reports grants from the NIH/NCI (2P01CA87969, principal investigator: M. Stampfer) during the conduct of the study. S.S. Tworoger reports grants from the NCI during the conduct of the study. No potential conflicts of interest were disclosed by the other authors.

The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. The authors assume full responsibility for analyses and interpretation of these data.

A.L. Shafrir: Formal analysis, methodology, writing–original draft. A. Babic: Formal analysis, methodology, writing–review and editing. M. Gates Kuliszewski: Writing–original draft, writing–review and editing. M.S. Rice: Data curation, writing–review and editing. M.K. Townsend: Writing–review and editing. J.L. Hecht: Conceptualization, data curation, visualization, writing–review and editing. S.S. Tworoger: Conceptualization, data curation, supervision, methodology, writing–review and editing.

The authors acknowledge the Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, as the home of the Nurses' Health Study. The authors thank the participants and staff of the Nurses' Health Study and Nurses' Health Study II for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, and WY. This project was supported by the NIH (P01 CA87969, to M. Stampfer and S.S Tworoger; UM1 CA186107, to M. Stampfer; U01 CA176726, to W.C. Willet; and R01CA54419 and P50CA105009, to D.W. Cramer) and the U.S. Department of Defense (W81XWH-10-1-02802, to K.L. Terry).

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.

1.
Kurman
RJ
,
Shih
IM
. 
The dualistic model of ovarian carcinogenesis revisited, revised, and expanded
.
Am J Pathol
2016
;
186
:
733
47
.
2.
Kroeger
PT
,
Drapkin
R
. 
Pathogenesis and heterogeneity of ovarian cancer
.
Curr Opin Obstet Gynecol
2017
;
29
:
26
34
.
3.
Köbel
M
,
Kalloger
SE
,
Boyd
N
,
McKinney
S
,
Mehl
E
,
Palmer
C
, et al
Ovarian carcinoma subtypes are different diseases: implications for biomarker studies
.
PLoS Med
2008
;
5
:
1749
60
.
4.
Wentzensen
N
,
Poole
EM
,
Trabert
B
,
White
E
,
Arslan
AA
,
Patel
AV
, et al
Ovarian cancer risk factors by histologic subtype: an analysis from the Ovarian Cancer Cohort Consortium
.
J Clin Oncol
2016
;
34
:
2888
98
.
5.
Shafrir
AL
,
Rice
MS
,
Gupta
M
,
Terry
KL
,
Rosner
BA
,
Tamimi
RM
, et al
The association between reproductive and hormonal factors and ovarian cancer by estrogen-α and progesterone receptor status
.
Gynecol Oncol
2016
;
143
:
628
35
.
6.
Hecht
JL
,
Kotsopoulos
J
,
Hankinson
SE
,
Tworoger
SS
. 
Relationship between epidemiologic risk factors and hormone receptor expression in ovarian cancer: results from the Nurses' Health Study
.
Cancer Epidemiol Biomarkers Prev
2009
;
18
:
1624
30
.
7.
Bossard
C
,
Busson
M
,
Vindrieux
D
,
Gaudin
F
,
Machelon
V
,
Brigitte
M
, et al
Potential role of estrogen receptor beta as a tumor suppressor of epithelial ovarian cancer
.
PLoS One
2012
;
7
:
e44787
.
8.
Lazennec
G
. 
Estrogen receptor beta, a possible tumor suppressor involved in ovarian carcinogenesis
.
Cancer Lett
2006
;
231
:
151
7
.
9.
Rutherford
T
,
Brown
WD
,
Sapi
E
,
Aschkenazi
S
,
Muñoz
A
,
Mor
G
. 
Absence of estrogen receptor-beta expression in metastatic ovarian cancer
.
Obstet Gynecol
2000
;
96
:
417
21
.
10.
Chan
KKL
,
Siu
MKY
,
Jiang
YX
,
Wang
JJ
,
Wang
Y
,
Leung
THY
, et al
Differential expression of estrogen receptor subtypes and variants in ovarian cancer: effects on cell invasion, proliferation and prognosis
.
BMC Cancer
2017
;
17
:
606
.
11.
De Stefano
I
,
Zannoni
GF
,
Prisco
MG
,
Fagotti
A
,
Tortorella
L
,
Vizzielli
G
, et al
Cytoplasmic expression of estrogen receptor beta (ERβ) predicts poor clinical outcome in advanced serous ovarian cancer
.
Gynecol Oncol
2011
;
122
:
573
9
.
12.
Halon
A
,
Nowak-Markwitz
E
,
Maciejczyk
A
,
Pudelko
M
,
Gansukh
T
,
Györffy
B
, et al
Loss of estrogen receptor beta expression correlates with shorter overall survival and lack of clinical response to chemotherapy in ovarian cancer patients
.
Anticancer Res
2011
;
31
:
711
8
.
13.
Burges
A
,
Brüning
A
,
Dannenmann
C
,
Blankenstein
T
,
Jeschke
U
,
Shabani
N
, et al
Prognostic significance of estrogen receptor alpha and beta expression in human serous carcinomas of the ovary
.
Arch Gynecol Obstet
2010
;
281
:
511
7
.
14.
Lenhard
M
,
Tereza
L
,
Heublein
S
,
Ditsch
N
,
Himsl
I
,
Mayr
D
, et al
Steroid hormone receptor expression in ovarian cancer: progesterone receptor B as prognostic marker for patient survival
.
BMC Cancer
2012
;
12
:
553
.
15.
Van Kruchten
M
,
Van Der Marel
P
,
De Munck
L
,
Hollema
H
,
Arts
H
,
Timmer-Bosscha
H
, et al
Hormone receptors as a marker of poor survival in epithelial ovarian cancer
.
Gynecol Oncol
2015
;
138
:
634
9
.
16.
Stampfer
MJ
,
Willett
WC
,
Speizer
FE
,
Dysert
DC
,
Lipnick
R
,
Rosner
B
, et al
Test of the National Death Index
.
Am J Epidemiol
1984
;
119
:
837
9
.
17.
Glynn
RJ
,
Rosner
B
. 
Methods to evaluate risks for composite end points and their individual components
.
J Clin Epidemiol
2004
;
57
:
113
22
.
18.
Chan
KKL
,
Wei
N
,
Liu
SS
,
Xiao-Yun
L
,
Cheung
AN
,
Ngan
HYS
. 
Estrogen receptor subtypes in ovarian cancer: a clinical correlation
.
Obstet Gynecol
2008
;
111
:
144
51
.
19.
Ciucci
A
,
Zannoni
GF
,
Travaglia
D
,
Petrillo
M
,
Scambia
G
,
Gallo
D
. 
Prognostic significance of the estrogen receptor beta (ERβ) isoforms ERβ1, ERβ2, and ERβ5 in advanced serous ovarian cancer
.
Gynecol Oncol
2014
;
132
:
351
9
.
20.
Dall
GV
,
Hawthorne
S
,
Seyed-Razavi
Y
,
Vieusseux
J
,
Wu
W
,
Gustafsson
J-A
, et al
Estrogen receptor subtypes dictate the proliferative nature of the mammary gland
.
J Endocrinol
2018
;
237
:
323
36
.
21.
Kass
L
,
Durando
M
,
Ramos
JG
,
Varayoud
J
,
Powell
CE
,
Luque
EH
, et al
Association of increased estrogen receptor beta2 expression with parity-induced alterations in the rat mammary gland
.
J Steroid Biochem Mol Biol
2004
;
91
:
29
39
.
22.
Keiler
AM
,
Bernhardt
R
,
Scharnweber
D
,
Jarry
H
,
Vollmer
G
,
Zierau
O
. 
Comparison of estrogenic responses in bone and uterus depending on the parity status in Lewis rats
.
J Steroid Biochem Mol Biol
2013
;
133
:
101
9
.
23.
Herrera-Lopez
EE
,
Castelan-Martinez
OD
,
Suarez-Sanchez
F
,
Gomez-Zamudio
JH
,
Peralta-Romero
JJ
,
Cruz
M
, et al
The rs1256031 of estrogen receptor β gene is associated with type 2 diabetes
.
Diabetes Metab Syndr
2018
;
12
:
631
3
.
24.
Meza-Muñoz
DE
,
Fajardo
ME
,
Pérez-Luque
EL
,
Malacara
JM
. 
Factors associated with estrogen receptors-alpha (ER-alpha) and -beta (ER-beta) and progesterone receptor abundance in obese and non obese pre- and post-menopausal women
.
Steroids
2006
;
71
:
498
503
.
25.
Beral
V
,
Hermon
C
,
Peto
R
,
Reeves
G
,
Brinton
L
,
Marchbanks
P
, et al
Ovarian cancer and body size: individual participant meta-analysis including 25,157 women with ovarian cancer from 47 epidemiological studies
.
PLoS Med
2012
;
9
:
e1001200
.
26.
Huang
T
,
Tworoger
SS
,
Willett
WC
,
Stampfer
MJ
,
Rosner
BA
. 
Associations of early life and adulthood adiposity with risk of epithelial ovarian cancer
.
Ann Oncol
2019
;
30
:
303
9
.