Association of reproductive history with breast tissue characteristics and receptor status in the normal breast

Introduction Reproductive history has been associated with breast cancer risk, but more knowledge of the underlying biological mechanisms is needed. Because of limited data on normal breast tissue from healthy women, we examined associations of reproductive history and established breast cancer risk factors with breast tissue composition and markers of hormone receptors and proliferation in a nested study within the Karolinska Mammography project for risk prediction for breast cancer (Karma). Materials and methods Tissues from 153 women were obtained by ultrasound-guided core needle biopsy as part of the Karma project. Immunohistochemical staining was used to assessed histological composition of epithelial, stromal and adipose tissue, epithelial and stromal oestrogen receptor (ER) and progesterone receptor (PR) status, and Ki-67 proliferation status. An individualised reproductive score including parity, number of pregnancies without birth, number of births, age at first birth, and duration of breastfeeding, was calculated based on self-reported reproductive history at the time of the Karma study entry. All analyses were adjusted for age and BMI. Results Cumulated reproductive score was associated with increased total epithelial content and greater expression of epithelial ER. Parity was associated with greater epithelial area, increased epithelial–stromal ratio, greater epithelial ER expression and a lower extent of stromal proliferation. Increasing numbers of pregnancies and births were associated with a greater epithelial area in the entire study set, which remained significant among postmenopausal women. Increasing numbers of pregnancies and births were also associated with a greater expression of epithelial ER among postmenopausal women. Longer duration of breastfeeding was associated with greater epithelial area and greater expression of epithelial PR both in the entire study set and among postmenopausal women. Breastfeeding was also positively associated with greater epithelial ER expression among postmenopausal women. Prior use of oral contraceptives was associated with lower epithelial–stromal ratio amongst all participants and among pre- and postmenopausal women separately. Conclusion Reproductive risk factors significantly influence the epithelial tissue compartment and expression of hormone receptors in later life. These changes remain after menopause. This study provides deeper insights of the biological mechanisms by which reproductive history influences epithelial area and expression of hormone receptors, and as a consequence the risk of breast cancer. Electronic supplementary material The online version of this article (10.1007/s10549-018-4768-0) contains supplementary material, which is available to authorized users.


Introduction
Reproductive history, like parity, age at first birth and number of births, has consistently been shown to be associated with breast cancer risk [1]. Women who have undergone a full time pregnancy before 20 years of age, for example, have a 50% reduced lifetime risk of developing breast cancer when compared to nulliparous women [2]. Following sex and age, mammographic density is considered one of the strongest risk factors for female breast cancer [3][4][5][6]. Mammographic density reflects the composition of supporting connective-(stromal), glandular epithelial-and adipose tissue of the breast. Reproductive history influences the mammographic density [7], thus suggesting a potential mechanism for the effect on risk by breast tissue composition alteration.
During puberty, pregnancy and postmenopausal involution, important shifts occur in normal breast tissue composition in response to changing hormone levels [8][9][10][11][12]. Epithelial proportions increase with hormonal exposure during pregnancy and lactation, whereas menopausal involution of the breast is associated with regression of lobules and relative increase of adipose tissue [8,13]. Although the underlying biology behind reproductive risk factors and breast cancer risk is not fully understood, proposed mechanisms include increased maturation of the breast parenchyma, greater differentiation of epithelial cells, and change in oestrogen responsiveness of the mammary gland [14]. The frequency of epithelial cells expressing the proliferation marker Ki-67 has also been positively associated with breast cancer risk among premenopausal women [15]. It remains inconclusive however, if parity and reproductive history alter epithelial proliferation [16][17][18].
A few studies, with conflicting results, have investigated how reproductive behaviour influences oestrogen receptor (ER) and progesterone receptor (PR) expression in normal mammary tissue [19][20][21][22]. Changing hormone levels caused by reproductive history may also influence the stromal composition, but previous studies show inconsistent results [8,[23][24][25]. We have previously found associations between epithelial and stromal hormone receptors and breast tissue characteristics, thus supporting the idea of hormonal regulation of tissue composition [26].
To understand the mechanisms by which reproductive history influences breast cancer risk, better knowledge of how these factors influence the breast tissue composition is needed. In this study, we examined associations of reproductive history and established breast cancer risk factors with breast tissue composition and markers of hormone receptors and proliferation in a nested study within the Karolinska Mammography project for risk prediction for breast cancer (Karma) [27].

Study population and tissue samples
Karma (Karolinska Mammography project for risk prediction for breast cancer) is a population-based prospective cohort study initiated in January 2011 and comprises 70,877 women attending mammography screening or clinical mammography at four hospitals in Sweden [27,28]. Women participating in the Karma project were recruited for this study. The methods used to collect this material have been described in details elsewhere by Gabrielson et al. [26,28] and will be given only in brief here. In total, 153 healthy women without prior history of breast cancer, other cancer or breast surgery, were included in the study. Characteristics of all participants are found in Table 1. Information on risk factors and exposures were collected by questionnaire at study enrolment [28]. BMI was calculated based on self-reported height and weight, and was missing for two participants. Areas of dense breast tissue were located by ultrasound and core needle biopsies were extracted from the densest part of the left breast. Tissues were formalin fixed and paraffin embedded. All participants signed an informed consent and the ethical review board at Karolinska Institutet approved the study.

Image analysis
Preparation of slides and image analysis has been described previously [26]. In brief, for each block a single haematoxylin and eosin section was prepared, on which a pathologist examined the tissue to confirm normal histology. All IHC sections were scanned in an Aperio ScanScope XT slide scanning system (Aperio Technologies, USA) at × 40 magnification and blinded before manually read using the Image-Scope viewing software. Digital images of the sections were captured using the ScanScope photo tool. Stromal, epithelial and adipose areas were manually selected using the Annotations tool. The total areas of analysis (comprising stromal, epithelial and adipose tissues) generated through the markup images were used to calculate the proportion of each tissue type. Calculations of epithelial and stromal protein expressions were done without knowledge of the different exposures of the study and have been described previously [26]. In brief, per cent epithelial nuclear expression of ER, PR and Ki-67 was assessed manually for each section separately. Likewise, stromal nuclear expression of ER, PR and Ki-67 was manually categorised as positive (≥ 1% positive cells) or negative for each section. Detailed information on tissue and protein characteristics of the cohort is found in Table 1.

Statistical analysis
Breast cancer risk factors and lifestyle factors included in the study were age at biopsy, BMI, use of oral menopausal hormone therapy, menopausal status (pre-or postmenopausal), previous benign breast disease history (defined as ever having any biopsy or partial removal procedure, with no cancer detected), age at menarche, parity, number of pregnancies, and number of births. Age at first birth and minimum duration of breastfeeding was assessed among parous women. Biopsy-based outcomes were total percentage epithelial and stromal content in the biopsy, total percentage epithelial immunohistological expressions of ER, PR and Ki-67, and positive or negative stromal expression of ER, PR and Ki-67. Linear regression was used to examine the association between established breast cancer risk factors and reproductive history with epithelial, stromal and adipose tissue measures expressed as percentages of the total area of the section, and percentage expressions of epithelial ER, PR and Ki-67 made from the histological sections. To study the balance in tissue composition between the epithelial and the stromal compartment, the ratios of epithelial-to-stromal proportions were evaluated by linear regression.
The assumption of normal distribution was inspected in all analyses and log transformations were made when necessary. All linear regression analyses of epithelial and stromal tissue compartment distributions and epithelial protein expressions were adjusted for age and BMI unless otherwise indicated.
Logistic regression was used to examine the association between established breast cancer risk factors and reproductive history with stromal protein expression. The analyses of stromal protein expressions were adjusted for age.
We created a reproductive score reflecting the reproductive activity of each study participant. The individual score was calculated by summarising the number of reproductive events (range total score 0-13) and it was defined by including the variables parous (no = 0, yes = 1), number of pregnancies without births (range 0-4), number of births (range 0-5), age at first birth (categories: nulliparous = 0; < 20 = 1; 20-25 = 2; > 25 = 3) and minimum duration of breastfeeding (categories: 0 = 0; < group mean = 1; > group mean = 2). Analyses of proportional growth by quartile groups of epithelial area (%) and epithelial ER (%) with reproductive score were assessed using univariate analysis of variance. Linear associations between percentage epithelial area or epithelial ER and reproductive score (continuous) were analysed by linear regression. Odds ratios (ORs) were evaluated using multinomial logistic regression comparing quartile groups of epithelial area or epithelial ER (cut points based on the overall distribution and the lowest 25% as the reference group in all analyses) with reproductive score for all women and stratified by menopausal status.
In secondary analyses, we stratified by menopausal status at biopsy. In total, 55 women were premenopausal at biopsy collection, and 95 were postmenopausal. Menopausal status was missing for three participants.
Two-tailed p values were used for all analyses with a p value of less than 0.05 considered to be statistically significant. Statistical analyses were performed using SPSS IBM, version 23.

Association of risk factors with epithelial area
Distributions of demographics, menstrual and reproductive history, and selected risk factors for breast cancer are described in Table 1. Histological sections from a total of 153 individuals were examined for association of risk factors with the epithelial tissue compartment ( Table 2).
In the total dataset, women at higher age showed smaller epithelial area and the epithelial area seemed to decrease by 2.2% per increasing age at biopsy (p = 0.049). Greater BMI was associated with a 5.3% smaller epithelial area per increasing BMI unit (p = 0.026), which remained significant among postmenopausal women (p = 0.038). Postmenopausal parous women had a 70.5% greater epithelial area compared to nulliparous postmenopausal women (p = 0.037) ( Table 2). Increasing number of pregnancies, number of births and longer duration of breastfeeding were all associated with a significantly greater epithelial area (17.5% per pregnancy; 21.0% per birth; 2.8% per month;  Fig. 1).
To better capture the combined effects of reproductive factors on the epithelial tissue compartment, we created a reproductive score reflecting the reproductive activity of each study participant. Total per cent epithelial area increased by 43.4% (p < 0.001) per quartile increase in reproductive score when adjusted for age and BMI in the total data set, which remained similar among postmenopausal women (39.9%; p = 0.001) (Supplemental Table S1).
Among all women, a higher reproductive score was positively associated with a 11.1% greater epithelial content per unit increase (p for trend = 0.001) (Table 3), with the OR in the highest quartile compared with the lowest quartile being 1.35 (95% CI 1.13-1.62) when adjusting for age and BMI. The findings remained significant for postmenopausal women (11.0% increase per unit reproductive score) (p for trend = 0.002) ( Table 3, Fig. 1), with the OR in the highest quartile compared with the lowest quartile being 1.28 (95% CI 1.04-1.57).

Association of risk factors with epithelial protein expression
In the total data set, higher age was significantly associated with a 1.0% higher expression of ER per year increase at biopsy (p = 0.033) ( Table 4) and a 2.2% lower expression of PR per year (p = 0.049) (Supplemental Table S2). Furthermore, higher BMI was associated with a significantly higher expression of PR (6.4%; p = 0.041). Increasing number of births was significantly associated with greater expression of ER in the total data set (10.0% per birth; p = 0.017), which remained significant among postmenopausal women (13.2% per birth; p = 0.013) ( Table 4, Fig. 2). Parous postmenopausal women also had 40.8% higher ER expression (p = 0.026), and greater ER by increasing number of pregnancies (7.4% per pregnancy; p = 0.035). Minimum duration of breastfeeding was positively associated with a 3.2% higher expression of PR in all women per additional month of breastfeeding (p = 0.003), with similar results among postmenopausal women (4.4%; p = 0.004). Among postmenopausal women, breastfeeding was also associated with greater ER (1.1%; p = 0.040).
In the total data set, a higher reproductive score was positively associated with a 3.0% greater epithelial ER expression per unit increase (p for trend = 0.043), with the OR in the highest quartile of reproductive score (median ER 43.0%) compared with the lowest quartile (median ER 16.1%) being 1.29 (95% CI 1.04-1.61), when adjusting for age and BMI. When stratifying be menopausal status, a higher reproductive score was positively associated with 4.5% greater epithelial ER expression per unit increase (p for trend = 0.017) (OR 1.239; 95% CI 0.953-1.612 for the highest quartile) when adjusting for age and BMI (Table 5, Fig. 2) among postmenopausal women. We did not see any significant associations between reproductive score and epithelial ER among premenopausal women.
Amongst premenopausal women, we found no significant associations between any of the risk factors and epithelial hormone receptor expressions (Table 4 and Supplemental  Table S2). Likewise, we found no significant associations between any of the risk factors and epithelial Ki-67 (Supplemental Table S3).

Association of risk factors with stromal compartment and stromal protein expression
The only risk factor significantly associated with stromal area in our study was BMI, which was associated with a smaller stromal area among all women (p < 0.001) (Supplemental Table S4). This finding remained significant after stratifying into premenopausal and postmenopausal women (p = 0.002) and (p = 0.035), respectively. We also investigated the epithelial-stromal ratio in our cohort. In the total dataset, increasing age was associated with a reduction of the epithelial-stromal ratio (p = 0.008) (Supplemental Table S5). Higher age at first birth was associated with a greater epithelial-stromal ratio in the total cohort (p = 0.004) and among premenopausal women (p = 0.012). Among postmenopausal women, parous women had a greater epithelial-stromal ratio (p = 0.036); furthermore, greater epithelial-stromal ratio was associated with higher number of births (p = 0.021). Prior use of oral contraceptives was associated with lower epithelial-stromal ratio amongst all participants and after stratification by menopausal status (total, p = 0.020; premenopausal, p = 0.034; and postmenopausal, p = 0.028). Table 3 Multinomial logistic regression of quartiles per cent epithelial area expressed as percentage of the total area of the section and reproductive score for all women and stratified by menopausal status Adjusted for age and BMI  In the total dataset, nulliparous women more frequently exhibited positive stromal expression of Ki-67 (p = 0.040). Likewise, higher BMI was associated with expression of stromal Ki-67 among postmenopausal women (p = 0.030). Lack of stromal PR was associated with higher age at first birth among all women (p = 0.012), and in postmenopausal women (p = 0.042).

Discussion
We found that reproductive history significantly influences the epithelial tissue content and expression of hormone receptors in later life and that these changes remain after menopause. Specifically, our results show that parity, number of pregnancies, number of births and duration of breastfeeding all were individually associated with a greater proportion of epithelial tissue after adjusting for age and BMI. Consistently, we found a strong positive relationship between epithelial content and the reproductive score we created reflecting the combined effects of reproductive history.
The human breast is defined by the degree of complexity of the secretory lobules, generally categorised as type 1-4, in the order of increasing complexity [10,12,29], although the cut points defining lobule types remain under discussion. The non-cancerous breast of parous women contains more type 3 lobules with a concomitant reduction in type 1 lobules [10]. Parous women also have more terminal duct lobular units (TDLU), and an increasing number of live births is associated with increasing TDLU counts [30]. Intuitively, the dramatically increased exposure to female sex hormones during each pregnancy could increase epithelial proliferation and thus the total amount of epithelial glandular tissue. To further support our findings, greater epithelial proportions are associated with parity and number of births in age-adjusted models [8]. Other studies, however, have not found any associations between reproductive history and glandular area from forensic autopsies [23] or epithelial content in non-neoplastic adjacent breast tissues from invasive breast tumours [25]. We did not investigate lobule types in this study, although we would expect that the epithelial proportion would increase with more complex branching lobule types, which by definition contain greater amounts of epithelia. A few studies, with conflicting results, have evaluated associations between parity and/or reproductive risk factors and expression of ER in human breast tissues [19][20][21][22]. The heterogeneity between results in the previous studies and the present one may reflect different statistical approaches and number of observations, menopausal status of subjects, and the type of tissue samples used. Expression of ER in TDLUs in near proximity to breast tumours has been positively associated with ER expression in the tumour, as compared with TDLUs further away [31]. Field effects surrounding the breast cancer may thus influence ER expression, and should be taken into consideration when analysing benign epithelium adjacent to the tumour. We only included non-diseased women without prior history of breast cancer or other cancer in our study. Further studies are needed to determine the relationship between parity and reproductive history with epithelial expression of ER in the normal human epithelium.
While expression of ER in benign epithelium was initially associated with breast cancer risk [20], recent studies have failed to confirm this [15,[31][32][33]. Circulating oestradiol levels and ER expression in normal epithelium is inversely associated [31,34], and parity and number of births are associated with lower concentrations of circulating oestrogens [35,36]. These findings suggest that the positive association between reproductive history and epithelial ER expression found in our study more likely reflects lower levels of circulating oestrogens rather than a direct causal relationship between parity and epithelial ER. Interestingly, although epithelial ER expression was associated with reproductive history, we found no association between reproductive history and epithelial proliferation after adjusting for age and BMI. In fact, epithelial proliferation is inversely associated with ER [34], and the frequency of Ki-67-expressing cells has been positively associated with increased risk of breast cancer, where high Ki-67+/low ER+ cell frequency was significantly associated with a 4.5-fold higher risk of breast cancer compared to low Ki-67+/high ER+ cell frequencies [15]. Collectively, this suggests that epithelial expression of ER without an increased mitotic activity does not mediate increased risk; rather it is a cumulative effect reflecting lower levels of circulating oestrogens as a consequence of parity and reproductive history. Likewise, and in agreement with previous findings, the positive association between reproductive factors protective of breast cancer and epithelial area without a corresponding increase of epithelial proliferation suggests that a greater epithelial area per se does not mediate an increased risk. Parityassociated protection against breast cancer may thus be mediated through non-breast epithelial tissue specific mechanisms, such as reduced endogenous hormone levels through exhausted ovaries.
Duration of breastfeeding was associated with increased epithelial area later in life in our overall study population and among postmenopausal women. These findings are in agreement with that of Figueroa et al. [30] who found a positive association between ever having breastfed and the number of TDLUs in postmenopausal women, and others who found positive associations between fibroglandular tissue and duration of lactation [24,37]. On the other hand, our results conflict with those studies that found no association [38,39]. We did not observe any significant changes in stromal compartment associated with duration of breastfeeding, however; amount of adipose tissue was inversely associated with breastfeeding (Supplemental information). Duration of breastfeeding was also positively associated with epithelial expression of PR in our overall study population and among postmenopausal women. Breastfeeding is protective for the luminal A, B and basal-like subtypes of breast cancer, with the largest protective effect of breastfeeding on the risk of the latter [40]. Basal-like cancers are thought to originate from undifferentiated luminal progenitor cells [41]. Fully differentiated type 4 lobules do not form until the end of pregnancy and during lactation when the breast attains its maximum development and increase in lobule numbers and size [42]. Progesterone is necessary for the lobuloalveolar mammary gland development [43,44], and the luminal population of epithelial cells express more PR compared to the stem cell-enriched fraction [45,46]. After menopause, the type 4 lobules regress to predominantly type 1-lobules, although more differentiated [47]. Data also suggest that breastfeeding is associated with less involution in postmenopausal women [30]. Although we did not determine lobule types and differentiation in this study, we hypothesise that epithelial expression of PR would increase with more complex and differentiated branching lobule types. Mammographic density is inversely associated with reproductive history [48], suggesting that reproductive factors also influence the stromal compartment. Although we did not detect any significant associations between stromal content and reproductive history in this study, parous women tended to have lesser stroma compared to nulliparous women. Similarly, parity and a higher reproductive score were non-significantly associated with lower mammographic density, with the greatest effect observed among postmenopausal women (data not shown). Parity was, however, significantly associated with lower expression of Ki-67 in stromal cells, and increased epithelial-stromal ratio in postmenopausal women only, further substantiating the hypothesis of reproductive history-associated regulation of the stromal content. Notably, although we observed stromal expression of ER, we did not find any associations between stromal ER and reproductive history after adjusting for age.
Collectively, our results suggest that while both the epithelial and stromal compartment express ER they are differentially responsive to oestrogen. More studies are needed in this field to clarify through which mechanisms parity regulate the stromal compartment. To the best of our knowledge, this is the first study investigating how parity influences stromal proliferation.
Strengths of our study include analysis of normal breast tissues from women without any history of breast cancer, quantitative measures of breast tissue composition, ER, PR and Ki-67 expression, use of standardised and optimised sample collection procedures, and a rich collection of epidemiologic data through the Karma cohort. Although our study is relatively small, it is among the largest exploring relationships between breast cancer risk factors, tissue characteristics and marker expression in normal breast tissue. However, the restricted number of samples in stratified analyses, particularly among premenopausal women, may bias our results. Also, all epidemiological data were self-reported at the Karma study entry, which may lead to misreporting and missing information. All biopsies included in this study were ultrasound-guided and collected in the densest part of the breast, which may limit the generalizability of findings. Nonetheless, we have previously shown that these tissues are representative of breast tissue composition in the normal breast [26].
In conclusion, this study provides deeper insights of the biological mechanisms by which reproductive history influences the risk of breast cancer. Specifically, our results show that several reproductive risk factors significantly influence breast tissue composition through the epithelial tissue compartment and expression of hormone receptors in later life, and that these changes remain after menopause. Association of risk reducing reproductive risk factors with greater epithelial area and increased expression of ER, but without increased mitotic activity, suggest that these morphologic changes alone do not contribute to an increased risk of breast cancer.