- Research article
- Open Access
Associations of reproductive breast cancer risk factors with breast tissue composition
Breast Cancer Research volume 23, Article number: 70 (2021)
We investigated the associations of reproductive factors with the percentage of epithelium, stroma, and fat tissue in benign breast biopsy samples.
This study included 983 cancer-free women with biopsy-confirmed benign breast disease (BBD) within the Nurses’ Health Study and Nurses’ Health Study II cohorts. The percentage of each tissue type (epithelium, stroma, and fat) was measured on whole-section images with a deep-learning technique. All tissue measures were log-transformed in all the analyses to improve normality. The data on reproductive variables and other breast cancer risk factors were obtained from biennial questionnaires. Generalized linear regression was used to examine the associations of reproductive factors with the percentage of tissue types, while adjusting for known breast cancer risk factors.
As compared to parous women, nulliparous women had a smaller percentage of epithelium (β = − 0.26, 95% confidence interval [CI] − 0.41, − 0.11) and fat (β = − 0.34, 95% CI − 0.54, − 0.13) and a greater percentage of stroma (β = 0.04, 95% CI 0.01, 0.08). Among parous women, the number of children was inversely associated with the percentage of stroma (β per child = − 0.01, 95% CI − 0.02, − 0.00). The duration of breastfeeding of ≥ 24 months was associated with a reduced proportion of fat (β = − 0.30, 95% CI − 0.54, − 0.06; p-trend = 0.04). In a separate analysis restricted to premenopausal women, older age at first birth was associated with a greater proportion of epithelium and a smaller proportion of stroma.
Our findings suggest that being nulliparous as well as having a fewer number of children (both positively associated with breast cancer risk) is associated with a smaller proportion of epithelium and a greater proportion of stroma, potentially suggesting the importance of epithelial-stromal interactions. Future studies are warranted to confirm our findings and to elucidate the underlying biological mechanisms.
Breast cancer remains the most commonly diagnosed cancer in women in the USA and worldwide . The vast majority of breast tumors are carcinomas that arise from the breast epithelium. Sarcomas of the breast are exceedingly rare and are thought to originate from stromal components of the breast (< 1% of all breast tumors) . It has long been recognized that women with a greater proportion of fibroglandular breast tissue (combined epithelium and stroma) as reflected on a mammogram (also referred to as breast density) are at a greater risk of breast cancer . While several breast cancer risk factors are suggested to influence breast tissue composition and thus subsequent breast cancer risk, the epidemiological evidence on these relationships remains very limited.
Reproductive factors related to childbearing are also recognized as breast cancer risk factors. Parity, younger age at first birth, and breastfeeding are associated with reduced breast cancer risk [4,5,6,7,8]. A longer period between menarche and first pregnancy, on the other hand, is associated with increased breast cancer risk [9,10,11,12]. Whether any of these factors could influence adult breast tissue composition is unclear.
Some previous studies of associations between reproductive factors and mammographic breast density, a well-established strong risk factor reflective of relative amounts of fibroglandular vs. fatty tissue content on the mammogram, found inverse associations of parity and positive associations of age at first birth and duration of breastfeeding with breast density [13,14,15,16,17,18]. In our recent study of reproductive factors and breast density, parous women as compared to nulliparous women had lower percent breast density (proportion of fibroglandular tissue out of the total breast area), smaller absolute dense area (area of fibroglandular tissue), and greater non-dense area (area of adipose tissue). The positive associations of breastfeeding with absolute dense and non-dense areas were limited to premenopausal women, while the positive association of the age at first child’s birth with percent density and the inverse association with the non-dense area were limited to postmenopausal women . Despite this evidence on the associations of reproductive factors with tissue composition on mammograms, only a few studies have examined these associations using direct measurement of tissue components in the non-malignant breast tissue of cancer-free women. Gabrielson et al. found positive associations of parity and duration of breastfeeding with epithelial area , while an earlier study on associations of reproductive factors with the proportion of epithelium or stroma found no such associations .
In this study, we aimed to assess the associations of several reproductive variables (parity, age at first birth, breastfeeding, time since last pregnancy, and duration of the time between menarche and first birth) with the extent of epithelial, stromal, fibroglandular, and fat tissue in non-malignant breast tissue from benign breast biopsy samples using prospective data in cancer-free women from the Nurses’ Health Study (NHS) and Nurses’ Health Study II (NHSII) and a deep-learning computational pathology method for tissue composition assessment.
Materials and methods
Our analysis included cancer-free women (controls) from the nested case-control study of breast cancer conducted among the subcohort of women with biopsy-confirmed benign breast disease (BBD) in the NHS and NHSII cohorts [22, 23]. These prospective cohorts followed registered nurses in the USA who were 30–55 years (NHS) or 25–42 years old (NHSII) at enrollment. After the administration of the initial questionnaire, the information on breast cancer risk factors (body mass index [BMI], reproductive history, postmenopausal hormone [PMH] use, and alcohol use) and any diagnoses of cancer or other diseases (including BBD) was updated through biennial questionnaires which were then confirmed via medical record review [13, 24]. Details of this nested case-control study and the BBD assessment have been previously described [22, 23].
Early NHS questionnaires (1976, 1978, and 1980) asked whether the participant had ever been diagnosed with “fibrocystic disease” or “other BBD” and whether she had been hospitalized in relation to this diagnosis. Beginning in 1982, the NHS questionnaires specifically asked about a history of biopsy-confirmed BBD (fibrocystic disease or other BBD). The initial 1989 NHS II questionnaire and all subsequent biennial questionnaires also asked the participants to report any diagnosis of BBD and to indicate whether it was confirmed by biopsy or aspiration.
Cases were women with biopsy-confirmed BBD who reported a diagnosis of breast cancer during 1976–1998 for the NHS and 1989–1999 for the NHSII following their BBD diagnosis. Using incidence density sampling, four women with biopsy-confirmed BBD who were free of breast cancer at the time of the matching case’s diagnosis (controls) were matched to the respective case on the year of birth and year of benign breast biopsy . We attempted to obtain BBD pathology records and archived biopsy specimens for all cases and controls from their hospital pathology departments; our ability to obtain biopsy blocks did not significantly differ by case and control status. Women were excluded if they had evidence of in situ or invasive carcinoma or unknown lesion type at the time of benign breast biopsy (22 cases and 12 controls). Only controls from this nested case-control study were used to examine the associations of reproductive factors with the extent of different tissue types. Out of 1907 controls, 983 had tissue readings, and information on reproductive factors was included in this analysis. Women with and without available tissue readings (as explained under the “Whole-slide image acquisition” section) had similar distributions of breast cancer risk factors.
The study protocol was approved by the institutional review boards of the Brigham and Women’s Hospital and Harvard T.H. Chan School of Public Health and those of participating registries as required. Consent was obtained or implied by the return of questionnaires.
Benign breast biopsy confirmation and BBD subtypes
Hematoxylin and eosin (H&E) breast tissue slides were retrieved for biopsy-confirmed BBD patients who gave permission to review their biopsy records. The slides were independently reviewed by one of three pathologists in a blinded fashion, i.e., the evaluating pathologists were blinded to the type of BBD noted on the original diagnosis [26, 27]. Any slide identified as having either questionable atypia or atypia was jointly reviewed by two pathologists. For each set of slides, a detailed worksheet was completed, and the benign breast biopsy was classified according to the categories of Page et al.  as non-proliferative, proliferative without atypia, or atypical hyperplasia (ductal or lobular hyperplasia) .
Whole-slide image acquisition
H&E slides were digitized into whole-slide images at × 20 (n = 93) or × 40 (n = 890) using the Panoramic SCAN 150 (3DHISTECH Ltd., Budapest, Hungary). For women with good-quality slides, up to six slides from different tissue blocks were digitized. H&E slides that were not digitized were due to poor quality, slides too thick to fit into the scanner, and plastic mounting coverslips. Attempts to create new H&E slides were not always possible due to missing (or returned to hospital) blocks, old style blocks not created using tissue cassettes, or poor-quality blocks . Slides were successfully digitized for approximately 80% of all control women in the original nested case-control study.
Quantification of epithelium, stroma, and fat
Whole-slide images were processed using a deep-learning computational pathology method to segment BBD tissues into epithelial, stroma, and fat regions. Tissue image analysis included normal terminal duct lobular units (TDLUs) and BBD lesions (referred to as “non-malignant” throughout this manuscript). Details of the image analysis method and its performance are described elsewhere . Briefly, to evaluate the tissue segmentation network, precision, recall, and Dice similarity coefficient were calculated using the held-out test set (n = 48). Dice similarity coefficient is the harmonic mean of precision (i.e., sensitivity) and recall (i.e., positive predictive value) and assesses how accurate the automated segmentation compares with ground truth on a pixel-wise basis. The range for Dice similarity coefficient is from 0 to 1, with 1 indicating perfect overlap. The majority of the precision, recall, and Dice similarity coefficient values of the tissue segmentation network and nuclei detection were > 0.75 . For more details about the nuclear segmentation network, please refer to the previously published methods paper by Vellal et al. .
For each whole-slide image, our method computed total, epithelial, stromal, and adipose tissue areas in pixels. We next calculated the average percent of each tissue type out of the total area across all available slides for each woman (median = 3, range 1–4), weighted by the total tissue area of the slides. We examined the associations of reproductive factors with the percentage of each of these individual tissue regions as well as combined epithelial and stromal tissue (fibroglandular area).
The data on age at menarche, parity, age at first birth, and breastfeeding were available from baseline and biennial questionnaires, completed closest to the date of the biopsy. Among all eligible controls with tissue readings, the completeness of the data on parity was 98.3%. Among parous women with tissue readings, information on age at first birth and breastfeeding was available for 98.4% and 94.1% of the sample, respectively. Age at menarche was available for 99.6% of the sample. For 4 women with missing age at menarche, a median value in the study sample was imputed, as done in previous studies [31,32,33]. Time since the last pregnancy was available for 86.3% of the study sample.
Age at first birth was categorized as < 25, 25–29, and ≥ 30 years. Parity was defined both as a binary variable (nulliparous, parous) as well as categorical (1, 2, 3, and ≥ 4 children). Additionally, the number of children among parous women was modeled as a continuous variable. Age at first birth was modeled both as a categorical (< 25, 25–29, and ≥ 30 years) and as a continuous variable. Information on breastfeeding was collected one time (on the 1986 questionnaire), and women were asked to report the total months of breastfeeding. Lifetime duration of breastfeeding (sum of breastfeeding duration across all births) was classified as none to < 1, 1 to < 12, 12 to < 24, and ≥ 24 months. Age at menarche was modeled both as a categorical (< 12, 12, 13, and > 13 years) and as a continuous variable. The time interval between menarche and first birth as well as the time since the last pregnancy was modeled as continuous variables.
Information on breast cancer risk factors was obtained from the biennial questionnaires closest to the date of the biopsy. Women were considered to be postmenopausal if they reported (1) no menstrual periods within the 12 months before biopsy with natural menopause, (2) bilateral oophorectomy, or (3) hysterectomy with one or both ovaries retained and were 54 years or older for ever smokers or 56 years or older for never smokers [34, 35].
We used multivariate linear regression to examine the associations of parity, age at first birth, breastfeeding, and the interval between menarche and first birth with the proportion of epithelial, stromal, fibroglandular, and fat tissues. Because tissue type measures were non-normally distributed, we used log-transformed values in all the regression analyses to improve normality. The risk estimates were adjusted for age (continuous), body mass index (BMI, continuous), a family history of breast cancer (yes vs. no), alcohol use (none, > 0 to < 5, ≥ 5 g/day), age at menarche (< 12, 12, 13, > 13), menopausal status/postmenopausal hormone use (pre-, post-/no hormones, post-/past hormone use, post-/current hormone use, post-/unknown hormone use status), and study cohort (NHS, NHSII). Additionally, in the analysis of the association of breastfeeding, the estimates were adjusted for parity and age at first birth. In the analysis of the associations of parity and age at first birth, the risk estimates were mutually adjusted for these two variables. In the analysis for the interval between menarche and first birth, the estimates were adjusted for parity.
The analyses of all reproductive variables except nulliparity and age at menarche were limited to parous women only. Parity, age at first birth, and age at menarche were modeled both as continuous and categorical, and breastfeeding was modeled as categorical. The lowest category for parity (1 child), age at first birth (< 25 years), and breastfeeding (0 to < 1 month) were used as the reference. To assess the overall trend for each of the categorical reproductive variables, we used respective medians within each category. The time since the last pregnancy as well as the duration of the interval between menarche and first birth were modeled as continuous variables. As in our prior study on reproductive factors and mammographic breast density, some of the associations were limited to either pre- or postmenopausal women; in the secondary analysis, we examined these associations separately in premenopausal women; the small sample of postmenopausal women (n = 290) in our study did not allow us to draw meaningful conclusions for this stratum. Finally, to account for the potential influence of BBD lesions on the study findings, we additionally adjusted all models for the type of the BBD.
In addition to the main approach, we used the SAS Proc Glimmix procedure which accounts for non-normal data distributions to examine the associations with tissue types using their original continuous scale. The analyses were performed using the SAS software (version 9.4, SAS Institute, Cary, NC, USA).
In this study of 983 cancer-free women, 299 (30.4%) had non-proliferative disease, 559 (56.9%) had proliferative disease without atypia, and 125 (12.7%) had atypical hyperplasia, consistent with previously reported distributions of these BBD subtypes . The distribution of different tissue types across these subtypes in our study is presented in Supplementary Table 1. The average proportion of epithelium, stroma, and fat in our study sample was 9.1% (range 0.5–52.2%), 72.4% (range 23.6–99.0%), and 18.5% (range 0–71.3%), respectively.
In our study sample, the average age at the biopsy was 42 years (range 19–58 years). A majority of the women were premenopausal at the biopsy (62.3%). The majority of women were parous (89.9% for premenopausal and 92.8% for postmenopausal), and the majority of parous women had at least two children (87.0% for premenopausal and 93.7% for postmenopausal) and breastfed for at least 1 month (58.8% for premenopausal and 50.4% for postmenopausal). The average age at first birth was 25 years (range 15–40 years) for premenopausal women and 25 years (range 16–37 years) for postmenopausal women. Age-adjusted characteristics of pre- and postmenopausal women in the study by nulliparous status are presented in Table 1.
In multivariate analysis (Table 2), being nulliparous was significantly associated with a reduced proportion of epithelium (nulliparous vs. parous β = − 0.26, 95% confidence interval [CI] − 0.41, − 0.11) and fat tissue (nulliparous vs. parous β = − 0.34, 95% CI − 0.54, − 0.13) and an increased proportion of stroma (nulliparous vs. parous β = 0.04, 95% CI 0.01, 0.08). The duration of breastfeeding of 24 months or longer was associated with a reduced proportion of fat (breastfeeding ≥ 24 vs. 0 to < 1 month, β = − 0.30, 95% CI − 0.54, − 0.06; p-trend = 0.04). As during pregnancy and lactation, terminal duct-lobular units undergo differentiation with a simultaneous reduction in surrounding fat cells , we explored the potential effect of postpartum involution on the observed associations of breastfeeding with the proportion of fat by additionally adjusting these models for time since the last pregnancy. We observed a slight attenuation of the effects (β = − 0.22, 95% CI − 0.47; 0.03 with adjustment vs. β = − 0.30, 95% CI − 0.5; − 0.06 without adjustment). The findings also did not reach statistical significance (p-trend = 0.10), likely due to the smaller number of observations in these models (790 vs. 816) due to missing data on the time since the last pregnancy. Finally, we found no interaction of time since the last pregnancy with breastfeeding (p for interaction = 0.40).
Parity was associated with a reduced proportion of stroma (β per one child = − 0.01, 95% CI − 0.02, − 0.00; p-trend 0.02), and having a first child at age 25–29 years was associated with a larger proportion of epithelial tissue (age at first birth 25–29 vs. < 25 years, β = 0.12, 95% CI 0.03, 0.21), though there was no clear pattern for this association. The time since the last pregnancy as well as the duration of the interval between age at menarche and first birth was not associated with any of the tissue measures. None of the reproductive factors was associated with the proportion of fibroglandular tissue. The patterns of associations of the reproductive factors with tissue types were similar in the statistical analyses with the secondary modeling approach (Supplementary Table 2).
Among premenopausal women (Table 3), being nulliparous was associated with a greater proportion of stroma (nulliparous vs. parous β = 0.06, 95% CI 0.02, 0.10) and a smaller proportion of epithelium (nulliparous vs. parous β = − 0.22, 95% CI − 0.38, − 0.06) and fat (nulliparous vs. parous β = − 0.32, 95% CI − 0.56, − 0.08). Greater parity and older age at first birth were both associated with a greater proportion of epithelium and a smaller proportion of stroma. The time since the last pregnancy as well as the duration of the interval between age at menarche and first birth was not associated with the proportion of any of the tissue types. These patterns of associations were similar with the secondary modeling approach (Supplementary Table 3). Finally, the findings did not change after additional adjustment for the BBD subtype (Supplementary tables 4 and 5).
In this study of 983 cancer-free women, being nulliparous was associated with having a smaller percentage of epithelium and fat and a greater percentage of stroma, while the greater number of children was associated with a smaller percentage of stroma. Breastfeeding for 24 months or longer was associated with a decreased percentage of fat. No associations were observed for the interval between menarche and age at first birth and the positive associations of age at first birth with the percentage of epithelium and inverse associations with the percentage of stroma were seen in premenopausal women.
Consistent with our findings on the associations of nulliparity with epithelium, a recent study by Gabrielson et al. of core biopsy samples of non-malignant breast tissue from 153 cancer-free women found a greater proportion of epithelium in parous as compared to nulliparous women (parous vs. nulliparous β = 0.56, p = 0.07) . In contrast, an earlier study by Gertig et al. within NHS found no associations of parity with the proportion of epithelium or stroma. Compared to our study, this study was small (n = 300) and used a different method of computer-assisted image analysis . Finally, even though neither of the previous studies found associations of parity with the proportion of stroma, Gabrielson et al. observed inverse associations with stromal proliferation .
Breast tissue changes during pregnancy have been suggested as possible reasons for the long-term protective effect on breast cancer risk . The influence of full-term pregnancy on the breast tissue appears to be complex and some of the suggested mechanisms include changes in hormonal signaling in the breast, gene methylation and expression changes, long-term reduction in the levels of circulating hormones, and life-long reduction in the number of mammary stem cells [4, 10, 38,39,40]. Previous studies also suggest that hormonal changes during pregnancy may also influence the stromal composition, but the evidence remains inconsistent [16, 41, 42]. Interestingly, in nulliparous women in our study, we observed a smaller proportion of epithelium but a larger proportion of stroma. These findings suggest that the increased risk of breast cancer in nulliparous women might potentially be driven by the dominating stroma and epithelial-stromal interactions that play a pivotal role in normal mammary gland function by controlling and regulating normal processes in the breast and suppressing the expression of preneoplastic phenotypes [43, 44]. Further, despite these findings for stroma and epithelium, we did not observe associations with the percentage of fibroglandular tissue which could further indicate that the relative contributions of stroma and epithelium might be more important than their absolute amount. Importantly, in our previous investigation of associations of reproductive factors with breast density, we observed greater percent breast density in nulliparous as compared to parous women. A few previous studies suggested that the degree of mammographic breast density is driven predominantly by the changes in the extent and composition of stroma rather than epithelium [45,46,47]. Thus, these relationships appear complex and warrant further investigation. Finally, as previous studies show heterogeneity in the association of reproductive factors with molecular breast cancer subtypes, our findings might help to explain eventually the underlying mechanisms behind these associations. However, as our study did not look at the associations with breast cancer outcomes and was limited to cancer-free women only, we are unable to discuss this further in the context of our findings.
In our study, older age at first birth was associated with a greater proportion of epithelium and a smaller proportion of stroma, but these findings were limited to premenopausal women. We did not find any associations of the length of the time period between menarche and first birth with any of the tissue types. Gabrielson et al. found a marginally positive association of age at first birth with the proportion of epithelium, but the results were not significant in the small stratum (n = 55) of premenopausal women , and Gertig et al. did not find any associations . Younger age at first birth may reduce subsequent breast cancer risk by earlier induction of cellular differentiation in the breast . However, it remains unclear if the observed associations may represent the result of the long-lasting effects of this differentiation.
We report inverse associations of breastfeeding duration with the proportion of fat tissue. Consistent with our findings, Gabrielson et al. found inverse associations, though only marginally significant associations of breastfeeding with the percentage of adipose tissue in the breast (β = − 0.55, p = 0.05) . Studies on breast tissue remodeling after lactation in humans are very limited though some studies suggest that the protective effect of breastfeeding on breast cancer risk may be related to the increased cellular proliferation and epithelial exfoliation of breast tissue during lactation with subsequent apoptosis after discontinuation of breastfeeding that could result in the elimination of cells which may have DNA damage [4, 10, 38,39,40]. Animal models suggest that with discontinuation of lactation, the breast tissue undergoes postpartum involution and remodeling as the result of apoptosis, regression of alveoli, and adipocyte repopulation [36, 49]. During lactation, mammary gland adipose tissue undergoes significant remodeling with replacement of adipocytes by mammary alveolar structures and their subsequent re-differentiation (“reversion”) back into adipocytes after weaning . Whether this reversion completely restores the previous tissue structure and whether this mechanism could be applied to humans is unknown. Finally, some recent studies also suggest that women who breastfed may have lower adiposity as reflected in their BMI which could also affect the amount of adipose tissue in the breast . In our study, however, we observed no correlation between BMI and breastfeeding (correlation coefficient = − 0.01, p = 0.75), and the estimates for the duration of breastfeeding in all models were adjusted for BMI.
To our knowledge, this is the largest study to date exploring the associations of several reproductive variables with the proportion of epithelium, stroma, fibroglandular, and fat tissues. The analysis used data from the Nurses’ Health Study and Nurses’ Health Study II, established cohorts with more than 30 years of follow-up, confirmed benign breast disease status and comprehensive information on breast cancer risk factors. Our study has a few limitations. Despite the prospective nature of the cohort, a measurement error for a few select reproductive variables especially in postmenopausal women is possible. For example, previous studies had conflicting findings on the accuracy of recall for age at menarche [52,53,54] which could potentially influence the results for associations of the interval between menarche and first birth with tissue measures. Some reports suggest that recall inaccuracy for breastfeeding in older women can affect the estimated associations between breastfeeding and health outcomes . Next, the proportion of nulliparous women in our study (10%) is comparable to the proportion among cancer-free women from other studies (range 4–14% [56,57,58,59]), including those in NHS/NHSII [60, 61]. As our study includes only cancer-free women with clinically indicated biopsy resulting in BBD diagnosis and since the analysis of the whole-slide images included both the background normal tissue as well as benign lesions, the findings are expected to be generalizable to cancer-free women with BBD. Next, as whole-slide images included both normal tissue and lesions, we additionally adjusted our models for BBD subtypes to explore their influence on the study findings, and the results remained unchanged. Finally, due to the small proportion of postmenopausal women in our study sample (n = 290), we were unable to perform an analysis within this stratum.
We investigated the associations of several reproductive variables related to childbearing with the extent of epithelial, stromal, fibroglandular, and fat tissue. Our findings suggest that nulliparous women are more likely to have a lower percentage of epithelium and fat and a greater percentage of stroma as compared to parous women, potentially suggesting the importance of epithelial-stromal interactions. Parous women with a greater number of children appeared to have a smaller proportion of stroma. In premenopausal women, younger age at first birth was associated with a larger proportion of epithelium and a smaller proportion of stroma. Future studies are warranted to confirm our findings and to elucidate the underlying biological mechanisms.
Availability of data and materials
Benign breast disease
Body mass index
Hematoxylin and eosin
Nurses’ Health Study
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Feng Y, Spezia M, Huang S, Yuan C, Zeng Z, Zhang L, et al. Breast cancer development and progression: risk factors, cancer stem cells, signaling pathways, genomics, and molecular pathogenesis. Genes Dis. 2018;5(2):77–106.
McCormack VA, dos Santos SI. Breast density and parenchymal patterns as markers of breast cancer risk: a meta-analysis. Cancer Epidemiol Biomarkers Prev. 2006;15(6):1159–69.
Britt K, Ashworth A, Smalley M. Pregnancy and the risk of breast cancer. Endocr Relat Cancer. 2007;14(4):907–33.
Morris GJ. Breastfeeding, parity, and reduction of breast cancer risk. Breast J. 2009;15(5):562–3.
MacMahon B, Cole P, Lin TM, Lowe CR, Mirra AP, Ravnihar B, et al. Age at first birth and breast cancer risk. Bull World Health Organ. 1970;43:209–21.
Trichopoulos D, Hsieh CC, MacMahon B, Lin TM, Lowe CR, Mirra AP, et al. Age at any birth and breast cancer risk. Int J Cancer. 1983;31:701–4.
Kelsey JL, Gammon MD, John EM. Reproductive factors and breast cancer. Epidemiol Rev. 1993;15(1):36–47.
Li CI, Malone KE, Daling JR, Potter JD, Bernstein L, Marchbanks PA, et al. Timing of menarche and first full-term birth in relation to breast cancer risk. Am J Epidemiol. 2008;167(2):230–9.
Ritte R, Tikk K, Lukanova A, Tjonneland A, Olsen A, Overvad K, et al. Reproductive factors and risk of hormone receptor positive and negative breast cancer: a cohort study. BMC Cancer. 2013;13:584.
Rosner B, Colditz GA. Nurses’ Health Study: log-incidence mathematical model of breast cancer incidence. J Natl Cancer Inst. 1996;88(6):359–64.
Pike MC, Krailo MD, Henderson BE, Casagrande JT, Hoel DG. ‘Hormonal’ risk factors, ‘breast tissue age’ and the age-incidence of breast cancer. Nature. 1983;303(5920):767–70.
Vachon CM, Kuni CC, Anderson K, Anderson VE, Sellers TA. Association of mammographically defined percent breast density with epidemiologic risk factors for breast cancer (United States). Cancer Causes Control. 2000;11(7):653–62.
Yaghjyan L, Mahoney MC, Succop P, Wones R, Buckholz J, Pinney SM: Relationship between breast cancer risk factors and mammographic breast density in the Fernald Community Cohort. Br J Cancer. 2012;106:996–1003.
Loehberg CR, Heusinger K, Jud SM, Haeberle L, Hein A, Rauh C, et al. Assessment of mammographic density before and after first full-term pregnancy. Eur J Cancer Prev. 2010;19(6):405–12.
Prebil L, Ereman R, Powell M, Jamshidian F, Kerlikowske K, Shepherd J, et al. First pregnancy events and future breast density: modification by age at first pregnancy and specific VEGF and IGF1R gene variants. Cancer Causes Control. 2014;25(7):859–68.
Sung H, Ren J, Li J, Pfeiffer RM, Wang Y, Guida JL, et al. Breast cancer risk factors and mammographic density among high-risk women in urban China. npj Breast Cancer. 2018;4(1):3.
Nakajima E, Iwase T, Miyagi Y, Fujita T, Ikeda N, Ishikawa T, et al. Association of parity and infant feeding method with breast density on mammography. Acad Radiol. 2020;27(2):e24–6.
Yaghjyan L, Colditz GA, Rosner B, Bertrand KA, Tamimi RM. Reproductive factors related to childbearing and mammographic breast density. Breast Cancer Res Treat. 2016;158(2):351–9.
Gabrielson M, Chiesa F, Behmer C, Rönnow K, Czene K, Hall P. Association of reproductive history with breast tissue characteristics and receptor status in the normal breast. Breast Cancer Res Treat. 2018;170(3):487–97.
Gertig DM, Stillman IE, Byrne C, Spiegelman D, Schnitt SJ, Connolly JL, et al. Association of age and reproductive factors with benign breast tissue composition. Cancer Epidemiol Biomark Prev. 1999;8(10):873–9.
Tamimi RM, Colditz GA, Wang Y, Collins LC, Hu R, Rosner B, et al. Expression of IGF1R in normal breast tissue and subsequent risk of breast cancer. Breast Cancer Res Treat. 2011;128(1):243–50.
Collins LC, Baer HJ, Tamimi RM, Connolly JL, Colditz GA, Schnitt SJ. The influence of family history on breast cancer risk in women with biopsy-confirmed benign breast disease: results from the Nurses’ Health Study. Cancer. 2006;107(6):1240–7.
Tamimi RM, Byrne C, Colditz GA, Hankinson SE. Endogenous hormone levels, mammographic density, and subsequent risk of breast cancer in postmenopausal women. J Natl Cancer Inst. 2007;99(15):1178–87.
Tamimi RM, Rosner B, Colditz GA. Evaluation of a breast cancer risk prediction model expanded to include category of prior benign breast disease lesion. Cancer. 2010;116(21):4944–53.
Jacobs TW, Byrne C, Colditz G, Connolly JL, Schnitt SJ. Pathologic features of breast cancers in women with previous benign breast disease. Am J Clin Pathol. 2001;115(3):362–9.
Tamimi RM, Byrne C, Baer HJ, Rosner B, Schnitt SJ, Connolly JL, et al. Benign breast disease, recent alcohol consumption, and risk of breast cancer: a nested case–control study. Breast Cancer Res. 2005;7(4):R555.
Page DL, Dupont WD, Rogers LW, Rados MS. Atypical hyperplastic lesions of the female breast. A long-term follow-up study. Cancer. 1985;55(11):2698–708.
Kensler KH, Liu EZF, Wetstein SC, Onken AM, Luffman CI, Baker GM, et al. Automated quantitative measures of terminal duct lobular unit involution and breast cancer risk. Cancer Epidemiol Biomark Prev. 2020;29(11):2358–68.
Vellal AD, Sirinukunwattan K, Kensler KH, Baker GM, Stancu AL, Pyle ME, et al. Deep learning image analysis of benign breast disease to identify subsequent risk of breast cancer. JNCI Cancer Spectrum. 2021;5(1):pkaa119.
Yaghjyan L, Mahoney MC, Succop P, Wones R, Buckholz J, Pinney SM. Relationship between breast cancer risk factors and mammographic breast density in the Fernald Community Cohort. British J Cancer. 2012;106(5):996–1003.
Moran O, Eisen A, Demsky R, Blackmore K, Knight JA, Panchal S, et al. Predictors of mammographic density among women with a strong family history of breast cancer. BMC Cancer. 2019;19(1):631.
Yaghjyan L, Colditz GA, Collins LC, Schnitt SJ, Rosner B, Vachon C, et al. Mammographic breast density and subsequent risk of breast cancer in postmenopausal women according to tumor characteristics. J Natl Cancer Inst. 2011;103(15):1179–89.
Willett W, Stampfer MJ, Bain C, Lipnick R, Speizer FE, Rosner B, et al. Cigarette smoking, relative weight, and menopause. Am J Epidemiol. 1983;117(6):651–8.
Stampfer MJ, Willett WC, Colditz GA, Rosner B, Speizer FE, Hennekens CH. A prospective study of postmenopausal estrogen therapy and coronary heart disease. N Engl J Med. 1985;313(17):1044–9.
Radisky DC, Hartmann LC. Mammary involution and breast cancer risk: transgenic models and clinical studies. J Mammary Gland Biol Neoplasia. 2009;14(2):181–91.
Meier-Abt F, Bentires-Alj M, Rochlitz C. Breast cancer prevention: lessons to be learned from mechanisms of early pregnancy–mediated breast cancer protection. Cancer Res. 2015;75(5):803–7.
Murrow LM, Weber RJ, Caruso J, McGinnis CS, Gascard P, Borowsky AD, Desai TA, Thomson M, Tlsty T, Gartner ZJ: Pregnancy and obesity modify the epithelial composition and hormone signaling state of the human breast. bioRxiv 2020:430611.
Ghosh S, Gu F, Wang C-M, Lin C-L, Liu J, Wang H, et al. Genome-wide DNA methylation profiling reveals parity-associated hypermethylation of FOXA1. Breast Cancer Res Treat. 2014;147(3):653–9.
World Cancer Research Fund and American Institute for Cancer Research. Food, nutrition, physical activity, and the prevention of cancer: a global perspective. Washington, DC: American Institute for Cancer Research; 2007.
Li T, Sun L, Miller N, Nicklee T, Woo J, Hulse-Smith L, et al. The association of measured breast tissue characteristics with mammographic density and other risk factors for breast cancer. Cancer Epidemiol Biomarkers Prev. 2005;14(2):343–9.
Sun X, Sandhu R, Figueroa JD, Gierach GL, Sherman ME, Troester MA. Benign breast tissue composition in breast cancer patients: association with risk factors, clinical variables, and gene expression. Cancer Epidemiol Biomark Prev. 2014;23(12):2810–8.
Barcellos-Hoff MH, Medina D. New highlights on stroma-epithelial interactions in breast cancer. Breast Cancer Res. 2005;7(1):33–6.
Kim JB, Stein R, O’Hare MJ. Tumour-stromal interactions in breast cancer: the role of stroma in tumourigenesis. Tumor Biol. 2005;26(4):173–85.
Alowami S, Troup S, Al-Haddad S, Kirkpatrick I, Watson PH. Mammographic density is related to stroma and stromal proteoglycan expression. Breast Cancer Res. 2003;5(5):R129–35.
Lin SJ, Cawson J, Hill P, Haviv I, Jenkins M, Hopper JL, et al. Image-guided sampling reveals increased stroma and lower glandular complexity in mammographically dense breast tissue. Breast Cancer Res Treat. 2011;128(2):505–16.
Warren R, Lakhani SR. Can the stroma provide the clue to the cellular basis for mammographic density? Breast Cancer Res. 2003;5(5):225.
Lambe M, Hsieh CC, Chan HW, Ekbom A, Trichopoulos D, Adami HO. Parity, age at first and last birth, and risk of breast cancer: a population-based study in Sweden. Breast Cancer Res Treat. 1996;38(3):305–11.
O’Brien J, Martinson H, Durand-Rougely C, Schedin P. Macrophages are crucial for epithelial cell death and adipocyte repopulation during mammary gland involution. Development. 2012;139(2):269–75.
Wang QA, Song A, Chen W, Schwalie PC, Zhang F, Vishvanath L, et al. Reversible de-differentiation of mature white adipocytes into preadipocyte-like precursors during lactation. Cell Metab. 2018;28(2):282–288.e283.
Armenta RF, Kritz-Silverstein D, Wingard D, Laughlin GA, Wooten W, Barrett-Connor E, et al. Association of breastfeeding with postmenopausal visceral adiposity among three racial/ethnic groups. Obesity. 2015;23(2):475–80.
Must A, Phillips SM, Naumova EN, Blum M, Harris S, Dawson-Hughes B, et al. Recall of early menstrual history and menarcheal body size: after 30 years, how well do women remember? Am J Epidemiol. 2002;155(7):672–9.
Koprowski C, Coates RJ, Bernstein L. Ability of young women to recall past body size and age at menarche. Obes Res. 2001;9(8):478–85.
Cooper R, Blell M, Hardy R, Black S, Pollard TM, Wadsworth MEJ, et al. Validity of age at menarche self-reported in adulthood. J Epidemiol Commun Health. 2006;60(11):993–7.
Promislow JH, Gladen BC, Sandler DP. Maternal recall of breastfeeding duration by elderly women. Am J Epidemiol. 2005;161(3):289–96.
Aurin J, Thorlacius H, Butt ST. Age at first childbirth and breast cancer survival: a prospective cohort study. BMC Res Notes. 2020;13(1):9.
Minami Y, Ohuchi N, Taeda Y, Fukao A, Hisamichi S. Risk factors for benign breast disease according to histopathological type: comparisons with risk factors for breast cancer. Jpn J Cancer Res. 1998;89(2):116–23.
Tan M-M, Ho W-K, Yoon S-Y, Mariapun S, Hasan SN, Lee DS-C, et al. A case-control study of breast cancer risk factors in 7,663 women in Malaysia. Plos One. 2018;13(9):e0203469.
Butt S, Borgquist S, Anagnostaki L, Landberg G, Manjer J. Parity and age at first childbirth in relation to the risk of different breast cancer subgroups. Int J Cancer. 2009;125(8):1926–34.
Yaghjyan L, Wijayabahu A, Eliassen AH, Colditz G, Rosner B, Tamimi RM. Associations of aspirin and other anti-inflammatory medications with mammographic breast density and breast cancer risk. Cancer Causes Control. 2020;31(9):827–37.
Chen H, Yaghjyan L, Li C, Peters U, Rosner B, Lindström S, Tamimi RM: Association of interactions between mammographic density phenotypes and established risk factors with breast cancer risk, by tumor subtype and menopausal status. Am J Epidemiol. 2020;190(1):44–58.
This work was supported by the National Cancer Institute at the National Institutes of Health [CA240341 to L.Y., CA131332, CA175080, P01 CA087969 to R.M.T.; UM1 CA186107 to M.S.; U01 CA176726 to W.W.], Avon Foundation for Women, Susan G. Komen for the Cure®, Breast Cancer Research Foundation, and the BIDMC High School Summer Research Program to A.D.V.
We would like to thank the participants and staff of the NHS and NHSII 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. The authors assume full responsibility for the analyses and interpretation of these data.
This work was supported by the National Cancer Institute at the National Institutes of Health [CA240341 to L.Y.; CA131332, CA175080, P01 CA087969 to R.M.T.; UM1 CA186107 to M.S.; U01 CA176726 to W.W.], Avon Foundation for Women, Susan G. Komen for the Cure®, Breast Cancer Research Foundation, and the BIDMC High School Summer Research Program to A.D.V.
Ethics approval and consent to participate
The study protocol was approved by the institutional review boards of the Brigham and Women’s Hospital and Harvard T.H. Chan School of Public Health and those of participating registries as required. Consent was obtained or implied by the return of questionnaires.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Distribution of various tissue elements by the type of benign breast disease. Table S2. Associations of reproductive factors with tissue types in benign breast biopsy samples using proc glimmix procedure (Odds Ratios and 95% Confidence Intervals). Table S3. Associations of reproductive variables with tissue types in benign breast biopsy samples of premenopausal women using proc glimmix procedure (Odds Ratios and 95% Confidence Intervals). Table S4. Associations of reproductive variables with percentage of different tissue types (log-transformed) in benign breast biopsy samples (β coefficients and 95% Confidence Intervals), additionally adjusted for BBD subtype. Table S5. Associations of reproductive variables with percentage of different tissue types (log-transformed) in benign breast biopsy samples of premenopausal women (β coefficients and 95% Confidence Intervals), additionally adjusted for BBD subtype.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yaghjyan, L., Austin-Datta, R.J., Oh, H. et al. Associations of reproductive breast cancer risk factors with breast tissue composition. Breast Cancer Res 23, 70 (2021). https://doi.org/10.1186/s13058-021-01447-2
- Benign breast disease
- Age at first child
- Breast cancer risk