Skip to main content
  • Research article
  • Open access
  • Published:

Breast cancer risk factors, survival and recurrence, and tumor molecular subtype: analysis of 3012 women from an indigenous Asian population



Limited evidence, mostly from studies in Western populations, suggests that the prognostic effects of lifestyle-related risk factors may be molecular subtype-dependent. Here, we examined whether pre-diagnostic lifestyle-related risk factors for breast cancer are associated with clinical outcomes by molecular subtype among patients from an understudied Asian population.


In this population-based case series, we evaluated breast cancer risk factors in relation to 10-year all-cause mortality (ACM) and 5-year recurrence by molecular subtype among 3012 women with invasive breast cancer in Sarawak, Malaysia. A total of 579 deaths and 314 recurrence events occurred during a median follow-up period of ~ 24 months. Subtypes (luminal A-like, luminal B-like, HER2-enriched, triple-negative) were defined using immunohistochemical markers for hormone receptors and human epidermal growth factor receptor 2 (HER2) in conjunction with histologic grade. Hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between risk factors and ACM/recurrence were estimated in subtype-specific Cox regression models.


We observed heterogeneity in the relationships between parity/breastfeeding, age at first full-term pregnancy (FFP), family history, body mass index (BMI), and tumor subtype (p value < 0.05). Among luminal A-like patients only, older age at menarche [HR (95% CI) ≥15 vs ≤ 12 years = 2.28 (1.05, 4.95)] and being underweight [HRBMI < 18.5kg/m2 vs. 18.5–24.9kg/m2 = 3.46 (1.21, 9.89)] or overweight [HR25–29.9kg/m2 vs. 18.5–24.9kg/m2 = 3.14 (1.04, 9.50)] were associated with adverse prognosis, while parity/breastfeeding [HRbreastfeeding vs nulliparity = 0.48 (0.27, 0.85)] and older age at FFP [HR > 30 vs < 21 years = 0.20 (0.04, 0.90)] were associated with good prognosis. For these women, the addition of age at menarche, parity/breastfeeding, and BMI, provided significantly better fit to a prognostic model containing standard clinicopathological factors alone [LRχ2 (8df) = 21.78; p value = 0.005]. Overall, the results were similar in relation to recurrence.


Our finding that breastfeeding and BMI were associated with prognosis only among women with luminal A-like breast cancer is consistent with those from previously published data in Western populations. Further prospective studies will be needed to clarify the role of lifestyle modification, especially changes in BMI, in improving clinical outcomes for women with luminal A-like breast cancer.


In addition to impacting incidence, lifestyle and environmental risk factors for breast cancer may influence disease progression. Several studies have previously evaluated this question, with mixed results. While some studies have documented older age at menarche [1,2,3], early age at first full-term pregnancy (FFP) [4, 5] and nulliparity [6, 7] to be associated with adverse prognosis in breast cancer patients, others have reported better prognosis in relation to these risk factors [8,9,10,11,12]. Discrepancies in reported associations may be explained by differences in study populations, risk factor distributions, and potential confounders, but could also be due to heterogeneity inherent in breast cancer.

Findings from expression profiling studies have been used to classify breast cancers into intrinsic subtypes (i.e. luminal A, luminal B, human epidermal growth factor receptor 2 (HER2)-enriched, basal-like, and normal-like subtypes), which were associated with different prognoses [13] and can be corroborated by immunohistochemical (IHC) markers for hormone receptors (i.e. estrogen receptor (ER), progesterone receptor [PR]) and HER2. Recently, proxies of the extent of tumor proliferation have been endorsed to refine subgroups that recapitulate the intrinsic subtypes more accurately than using hormone receptors and HER2 alone [14, 15]. Epidemiological studies have shown that associations between breast cancer risk factors vary by tumor subtypes. For example, parity and early age at FFP are associated with decreased risk of luminal breast tumors, but they do not protect and may even increase the risk for ER-negative or triple-negative breast cancers [16,17,18,19].

Three previous studies have evaluated the relationship between breast cancer risk factors and survival according to molecular subtype, one among women in Seoul, South Korea [20] and the other two involving analyses of US-based prospective breast cancer cohorts [21, 22]. Results from these studies suggest that the associations between late age at menarche [20], breastfeeding [21], high body mass index (BMI) [22] and survival after breast cancer might differ according to molecular subtype. However, findings from these studies are yet to be validated in independent populations and, to our knowledge, no study has specifically examined risk factors in relation to survival according to subtypes defined by the recent IHC classification scheme accounting for proliferation in an Asian population.

Despite racial and geographic variations in the incidence, presentation, and outcome of breast cancer; so far, most investigations on risk factors in relation to tumor subtypes and survival have been conducted in European populations. This analysis, therefore, aims to evaluate the association between breast cancer risk factors and tumor molecular subtypes, defined by hormone receptors and HER2 in conjunction with histologic grade; and to examine the relationship between risk factors and survival by molecular subtype among women in Sarawak, Malaysia.


Study population

Sarawak is a Malaysian state on Borneo with a multiethnic composition, comprising of native Borneo populations (51%), Chinese (25%) and Malays (24%) [23]. Overall, 3355 women with invasive breast cancer diagnosed and treated between 2003 and 2016 in the Department of Radiotherapy, Oncology, and Palliative Care, Sarawak General Hospital where ~ 93% of all breast cancer cases diagnosed in Sarawak are treated, were recruited for this study. Of these, 106 (~ 3%) did not participate by not filling the questionnaire leading to a participating rate of ~ 97%. Of the 3249 who participated, 168 (~ 5%) were lost to follow-up and 69 did not have complete information on ER, PR, HER2, and grade that is needed to generate breast cancer subtypes hence were excluded from further analysis. Ultimately, 3012 women representing ~ 90% of the original population were included in the current analysis. Information on lifestyle and environmental risk factors were obtained from questionnaires that were administered to participants at enrollment, which was approximately 4 weeks after diagnosis, while information on tumor characteristics was obtained from clinical records. Weight and height measurements were obtained in the clinic as part of the clinical workup for the calculation of chemotherapy doses. Recordings were performed by a trained member of staff using a weighing scale. Patients were given follow-up appointments to the clinic during which recurrence was evaluated and clinically confirmed. For those living in the outskirts of the city, if recurrence was suspected, the patients were referred to our clinic for further evaluation. Furthermore, a research assistant made regular calls to check the patient’s status, whether alive or dead. The current analysis included a follow-up period of 153 months (median follow-up = 24 months). Ethical approval for this project was provided by the Ethics Committee of the National Institutes of Health, Malaysia. This study did not involve the use of personal identifying information; hence, it was exempted from review by the National Institutes of Health (NIH) Office of Human Subject Research Protections [23].

Breast cancer subtype definition

IHC staining for ER, PR, and HER2 was performed on formalin-fixed, paraffin-embedded tissue sections as has been previously described [24]. Molecular subtypes were defined using the St Gallen classification, proposed for the recapitulation of the intrinsic subtypes using IHC and proliferation markers [14, 15]. According to the St Gallen classification scheme, luminal breast cancers can be further distinguished into subgroups based on their level of proliferation (using KI67 or histologic grade) and hormone receptor expression patterns. Accordingly, luminal tumors that homogeneously express hormone receptors (i.e. ER+ and PR+) and low proliferation are classified as luminal A-like while those that heterogeneously express hormone receptors (i.e. ER+/PR or ER/PR+) and/or those that homogeneously or heterogeneously express hormone receptors (i.e. ER+ and/or PR+) but are also high proliferating (high KI67 or grade 3) and/or HER2+ are classified as luminal B. In keeping with this definition, we utilized ER, PR, and HER2 in addition to histologic grade [25], to define subtypes as follows: Luminal A-like: ER+ and PR+, HER2 and low grade (histologic grade 1 or 2); Luminal B-like: ER+ and/or PR+, HER2 and high-grade (histologic grade 3) or ER+ and/or PR+, HER2+ (regardless of levels of histologic grade); HER2-enriched: ER and PR and HER2+; and triple-negative: ER and PR and HER2.

Statistical analysis

Frequency tables were used to assess the distribution of risk factors and clinicopathological characteristics among the different subtypes. The chi-square test was used to assess differences for categorical variables and the Kruskal-Wallis test was used for continuous variables.

We categorized risk factors based on what is the convention for each variable and in accordance with what has been published in large-scale studies of breast cancer [18, 26]. We categorized age at menarche [≤12 years (early menarche), 13, 14 and ≥ 15 years (late menarche)]; family history of breast cancer in a first-degree relative [yes and no]; age at FFP [< 21, 21–24.9, 25–30, > 30 years] and age at menopause [≤ 50 and > 50 years] similarly as in previously published articles [18, 26]. For BMI, we adopted the World Health Organization classification [< 18.5 kg/m2 (underweight); 18.5–24.9 kg/m2 (normal weight); 25–29.9 kg/m2 (overweight) and ≥ 30 kg/m2 (obese)]. To test for associations between risk factors and molecular subtypes, we constructed a polytomous logistic regression model with tumor subtype as the outcome (luminal A-like subtype as the reference category) and risk factors (age at menarche [≤ 12 years (early menarche, reference category), 13, 14, and ≥ 15 years (late menarche)], parity and breastfeeding [nulliparity (reference category), parity but no breastfeeding, parity and breastfeeding], age at FFP [< 21 (reference category), 21–24.9, 25–30, > 30 years], family history [yes and no (reference category)], and BMI [< 18.5 kg/m2 (underweight); 18.5–24.9 kg/m2 (normal weight, reference category); 25–29.9 kg/m2 (overweight) and ≥ 30 kg/m2 (obese)] as explanatory variables, with adjustment for age at diagnosis (< 35, 35–45, 45–55, 55–65, 65–75, > 75 years) and ethnicity (Chinese, Malay, Native).

The association between breast cancer subtypes and all-cause mortality/recurrence was determined using Kaplan-Meier survival curves and Cox-proportional hazards regression models, which included adjustments for standard prognostic parameters including age at diagnosis, ethnicity, BMI, histologic grade, TNM stage I–IV [i.e. size (T), nodal status (N) and metastasis (M)], systemic therapy (endocrine (tamoxifen or aromatase inhibitor (AI)) and chemotherapy), radiotherapy and surgery. Follow-up started at diagnosis of breast cancer and ended at time of event (recurrence/death) or censoring (end of follow-up or, for the recurrence analysis, also death). For all-cause mortality, we censored at 10 years because this is the threshold at which most breast cancers are, by convention, considered cured in the absence of recurrence or death. We adopted a two-step approach in our survival analyses. In the first step, each of the above risk factors was modeled separately in basic models adjusted for standard prognostic factors separately for each tumor subtype. To test for heterogeneity in risk factor and survival relationships by subtype, we included an interaction term between each risk factor and tumor subtype. Violation of the proportionality assumption of the hazard model was tested by modeling each risk factor as a time-varying covariate. In the second step, it was decided, a priori, that factors that were associated with survival with P <  0.1 in the basic model were to be mutually adjusted for in a multivariable model that included the standard prognostic factors mentioned above. Using likelihood ratio (LR) test, we compared this model with one containing only the clinicopathological factors. For sensitivity analysis, we conducted survival analysis for women stratified into two age groups (< 50 yrs and ≥ 50 yrs). Also, we performed additional sensitivity analysis by excluding women with stage IV disease from our multivariate analyses for both all-cause mortality and recurrence. Results were very similar from these sensitivity analyses as compared to analyses including all women and we therefore presented results from all patients. All analyses were two-sided and performed using Stata statistical software version 14.0 (StataCorp, College Station, TX, USA).


In total, our analysis included 3012 invasive breast cancer cases, with a total of 579 deaths in 10 years and 314 recurrence events in 5 years. The mean age at diagnosis was 52 years and mean BMI was 25 kg/m2. The majority of the patients were Chinese (48%) and had early-stage (I and II, 56%) and HR-positive (66%) tumors (Table 1). Of the 3012 patients, 1016 (34%) were luminal A-like, 989 (33%) were luminal B-like, 387 (13%) were HER2-enriched, and 620 (20%) were triple-negative, respectively.

Table 1 Distribution of risk factors and clinicopathological characteristics by tumor subtype

Distribution of risk factors and clinicopathological characteristics by tumor subtype

As shown in Table 1, women with the luminal A-like subtype were slightly older than those with other subtypes. The distributions of ethnicity (P value < 0.001), age at FFP (P value = 0.019), breastfeeding practices (P value = 0.01), family history (P value = 0.05) and BMI (P value = 0.03) differed by subtype. No differences were observed in the distributions of age at menarche, age at menopause and parity according to subtype. The frequencies of all clinicopathological parameters differed by subtype, with low-grade, small, early-stage and node-negative tumors being more frequent for the luminal A-like subtype (Table 1).

Table 2 shows the associations between examined risk factors and molecular subtype in the multivariable polytomous regression model. Compared with women with the luminal A-like subtype, women with the luminal B-like, HER2-enriched and triple-negative tumors were significantly more likely to be Malay and Native than Chinese. Furthermore, women with other tumor subtypes were more likely to be parous and have breastfed [odds ratio (OR) (95% CI) parity and breastfeeding vs nulliparity = 1.44 (1.05, 1.98); 1.64 (1.06, 2.25); and 1.54 (1.07, 2.22) for luminal B, HER2-enriched and triple-negative subtypes, respectively] and less likely to experience their FFP after the age of 30 years [OR (95% CI) > 30 years vs < 21 years = 0.63 (0.42, 0.94); 0.57 (0.32, 1.02); and 0.58 (0.36, 0.93) for luminal B, HER2-enriched and triple-negative subtypes, respectively] than those with the luminal A-like subtype. Women with HER2-enriched [OR (95% CI) BMI > 30 kg/m2 vs 18.5–24.9 kg/m2 = 0.55 (0.34, 0.89); P value = 0.02] and triple-negative [OR (95% CI) BMI > 30 kg/m2 vs 18.5–24.9 kg/m2 = 0.59 (0.40, 0.88); P value = 0.01] tumors were significantly less likely to be obese than those with the luminal A-like subtype.

Table 2 OR and 95% CI from a polytomous logistic regression model testing the associations between breast cancer risk factors and tumor molecular subtype

Breast cancer risk factors in relation to all-cause mortality and recurrence by subtype

Overall, all-cause mortality and recurrence differed significantly by tumor subtype. In general, women with luminal A-like tumors had better survival outcomes than those with the other subtypes (Fig. 1). As shown in Table 3, in basic models for each risk factor (with adjustment for standard prognostic factors in addition to age, ethnicity and BMI), later age at menarche, parity/breastfeeding, and being underweight were significantly associated with 10-year all-cause mortality in the luminal A-like but not any of the other subtypes. Also, later age at FFP showed a suggestive association with mortality in luminal A-like patients (P trend = 0.08) but not the other subtypes. Results were similar in basic models for recurrence (Table 4).

Fig. 1
figure 1

Kaplan-Meier curves for the associations between breast cancer molecular subtypes and (a) 10-year all-cause mortality (b) 5-year recurrence-free survival among 3012 women diagnosed and treated in Sarawak General Hospital, Sarawak Malaysia

Table 3 HR and 95% CI for the associations between risk factors and 10-year all-cause mortality by tumor molecular subtype
Table 4 HR and 95% CI for the associations between risk factors and 5-year recurrence by tumor molecular subtype

In the multivariable model with the mutual adjustment for ethnicity, menarche, parity/breastfeeding, age at FFP, family history, and BMI in addition to standard clinicopathological factors and treatment variables, increasing age at menarche [hazard ratio (HR) (95% confidence interval (CI) ≥15 years vs ≤ 12 years = 2.28 (1.05, 4.95); P value for trend (P trend) = 0.06]; parity/breastfeeding [HR (95% CI) vs nulliparity = 0.48 (0.27, 0.85); P trend = 0.01]; older age at FFP [HR (95% CI) > 30 vs < 21 years = 0.20 (0.04, 0.90); P trend = 0.06]; and being underweight [HR (95% CI) vs normal weight = 3.46 (1.21, 9.89); P value = 0.02] or overweight [HR (95% CI) vs normal weight = 3.14 (1.04, 9.50); P value = 0.04] remained significantly associated with 10-year all-cause mortality in women with the luminal A-like subtype (Table 5 and Fig. 2). For these women, the addition of age at menarche, parity/breastfeeding, and BMI, provided significantly better fit to a model containing clinicopathological factors alone [LRχ2 (8df) = 21.78; P value = 0.005]. In general, the results were consistent in relation to recurrence (Table 5).

Table 5 Multivariate HR and 95% CI for the association between breast cancer risk factors and 10-year all-cause mortality and 5-year recurrence among women with luminal A-like subtype breast cancer
Fig. 2
figure 2

Survival curves for the multivariate association between (a) age at menarche, (b) parity and breastfeeding and (c) BMI and 10-year all-cause mortality among 1016 women with the luminal A-like subtype of breast cancer

When we examined the association between duration of breastfeeding and all-cause-mortality/recurrence for luminal A-like cases with complete information on breastfeeding duration (N = 719), we observed an inverse association between each breastfeeding duration category and all-cause mortality [HR (95% CI) vs nulliparity = 0.37 (0.18, 0.85), 0.86 (0.35, 2.11), 0.53 (0.24, 1.17) for < 6, 6–10, and > 10 months, respectively (P trend = 0.38)] and recurrence [HR (95% CI) vs nulliparity = 0.47 (0.19, 1.16), 0.73 (0.23, 2.37), 0.05 (0.01, 0.42) for < 6, 6–10, and > 10 months, respectively (P trend = 0.002)]. Among women who breastfed, all-cause mortality did not significantly vary by breastfeeding durations (comparing > 10 months to < 6 months, P value = 0.38) but women who breastfed for > 10 months tended to have better recurrence outcomes [HR (95% CI) vs < 6 months = 0.11 (0.01, 0.93); P value = 0.04].


In this study involving over 3000 invasive breast cancer cases from a population-based case series in Sarawak, Malaysia, with detailed demographic, risk factor, pathology, and follow-up data, we investigated several established breast cancer risk factors in relation to tumor subtypes and patient outcomes. We found differences in the prevalence of parity and breastfeeding, age at FFP, family history of breast cancer and obesity across different breast tumor subtypes. In general, traditional breast cancer risk factors (older age at FFP, higher BMI, lower parity, lack of breastfeeding) seem to show higher frequencies among women with the luminal A-like subtype compared with women with other subtypes. Further, we found that age at menarche, breastfeeding, and BMI were independent prognostic factors for both overall mortality and breast cancer recurrence but only for women with the luminal A-like subtype, which had better survival and recurrence outcomes than the other subtypes.

Our findings that nulliparity and older age at FFP were more prevalent in luminal A-like patients are consistent with those reported in studies in Western countries [16, 27]. However, unlike the observation of higher BMI and shorter breastfeeding duration in triple-negative patients among Western, especially African American, women [18, 19, 28, 29], we found lower frequencies of obesity and breastfeeding among HR (HER2-enriched and triple-negative) than luminal A-like tumors, which may be reflective of population/ethnic differences. In line with this hypothesis, a previous study conducted in South Korea [30] also showed a higher frequency of breastfeeding among women with luminal B or HER2-enriched than those with luminal A disease. In another study involving 730 Mexican women with breast cancer, Martinez and colleagues [31] reported the prevalence of breastfeeding to be higher among women with triple-negative than luminal A tumors. Similarly, results from a multiethnic study showed an inverse association between triple-negative tumors and breastfeeding in White, Hispanic and African American but, notably, not in Asian women, for whom breastfeeding for > 2 months was associated with an 86% increased likelihood of triple-negative tumors [32].

The prevalence of obesity is still much lower in most Asian populations compared with other race/ethnicity groups. In contrast to the reduced breast cancer risk associated with higher BMI among premenopausal Western women, obesity is associated with increased risk for both premenopausal and postmenopausal Asian women [33,34,35]. The heterogeneity of obesity by tumor subtype among Asian cases remains unclear. Results from both our study and the South Korean study [30] suggest that obesity was less frequent among patients with HER2-enriched tumors. In combination with our finding that women with the HER2-enriched subtype were more likely to be parous and to have breastfed than women with the luminal A-like subtype, our data suggest that these factors (parity, breastfeeding, and low BMI) may not protect against HER2-enriched breast cancers. The decreasing prevalence of these factors associated with the adoption of westernized lifestyles may, therefore, not affect the incidence of this subtype, which is known to be more prevalent among Asian women [36, 37]. More research to understand the risk factors associated with the HER2-enriched subtype is warranted.

Most epidemiological studies considering the prognostic significance of age at menarche in breast cancer have treated the disease as a homogeneous entity, and results from these studies are largely conflicting [1,2,3, 8, 38]. We found that older age at menarche was associated with worse prognosis but only among women with the luminal A-like subtype. This finding is consistent with results from a previous study of women in East Asia [20] that also evaluated relationships between risk factors and survival according to subtypes. It is unclear why late menarche and younger age at FFP lead to worse survival outcomes in women with the luminal A-like subtype since they are well known protective factors in terms of breast cancer risk. One possibility is that because early menarche and late age at FFP increase breast cancer risk through prolonged and sustained exposure of the mammary epithelium to the mitogenic effects of reproductive hormones [39, 40], these factors predispose more strongly to HR+ tumors which have better prognosis than HR tumors [18, 41]. Although this association was confined to the luminal A-like subtype, which, by definition, is HR+, expression of hormone receptors in tumors occurs in a spectrum. Whereas some tumors have very high expression levels, others have lower levels despite crossing the threshold for consideration as HR+. Due to differences in cumulative lifetime exposure to endogenous estrogens, luminal A-like tumors occurring among women with late menarche may have lower levels of hormone receptor expression, hence worse survival/recurrence outcomes, than those occurring among women with early menarche. Indeed, this is consistent with the finding by Song et al. that longer duration of endogenous estrogen exposure was associated with better survival [30]. On the other hand, late menarche and early age at FFP may be indicative of lower socioeconomic status (SES) which, in turn, may be reflective of less exposure to “westernized” environments/lifestyles. Nonetheless, adjustment of known surrogates for SES did not change the associations between these factors and survival. That late age at menarche and early age at FFP were not associated with survival/recurrence in the other subtypes may be due to the masking effect of other more aggressive tumor features, which are inherent in these subtypes.

Our findings of breastfeeding being associated with better outcomes are generally consistent with previous reports [2, 12, 42]. A distinctive feature of luminal A-like tumors defined in our analysis is the low levels of proliferation, indicated by low histologic grade. Findings from one previous study showed that the protective effect of breastfeeding on breast cancer mortality was stronger for tumors with low expression of proliferation genes [21]. In our study, breastfeeding was associated with a preponderance of lobular carcinoma and small size tumors, both of which are highly correlated with low levels of proliferation [43,44,45].

Results from several studies, summarized in two comprehensive reviews and meta-analyses [46, 47], are supportive of the prognostic value of BMI in breast cancer. The association between BMI and survival after breast cancer is thought to be U-shaped [47,48,49], with underweight and overweight/obese women more likely to suffer worse survival outcomes than their normal weight counterparts. In our study, being underweight and overweight, but not obese, were suggestively associated with worse prognosis, but this might be due to the low frequency of obesity in this population (13%). Compared with overweight/obesity, the effect of underweight in breast cancer survival is less well-studied. Overall, our finding of an association between underweight and worse survival outcomes in breast cancer is in line with those of several other studies involving Asian populations [49,50,51,52]. Whilst insulin resistance, chronic inflammation, and altered adipokine and cytokine production have been proposed to underlie the obesity-cancer link [53], the precise mechanisms underpinning the relationship between underweight BMI and disease progression are not well understood. Chronic, pre-diagnostic, malnutrition may contribute to weight loss in cancer patients and may independently influence outcomes in the disease. However, when we examined the impact of indicators of socioeconomic deprivation, as surrogates for chronic malnutrition, our estimates remained unchanged.

Strengths of this study include: a population-based breast cancer case series in an understudied Asian population and the collection of detailed questionnaire information, which allowed us to account for various confounding variables including sociodemographic factors. Several limitations of this study should be noted. First, we did not include controls in our study and therefore our case-case comparisons could not be translated into relative risk estimates. Second, relatively small sample sizes for some of the subtypes may have affected our power to detect significant associations. Third, data on the specific cause of death were not available and, therefore, we only evaluated risk factors in relation to all-cause, but not breast cancer-specific, mortality even though some of the risk factors, such as BMI, are important predictors of death from other causes [54]. Nonetheless, the consistency of our results for all-cause mortality and for breast cancer recurrence indicate that these factors may contribute to breast cancer-specific mortality in a similar manner.


In conclusion, our data indicate that risk factors for breast cancer are differentially associated with tumor subtypes and exert subtype-specific influence on survival/recurrence from the disease. Specifically, we observed that menarche after the age of 15 years, FFP after 30 years, underweight or overweight BMI, and breastfeeding practices were associated with survival/recurrence only among women with luminal A-like tumors. These findings are supportive of the prognostic value of reproductive and lifestyle-related factors in tumors with biologically favorable profiles, and could have implications for clinical counseling and for the development of subtype-specific prognostic tools. Future prospective studies are needed to delineate the role of lifestyle modification, especially changes in BMI, in improving clinical outcomes for women with luminal A-like breast cancer.



All-cause mortality


Body mass index


Confidence interval


Estrogen receptor


First full-term pregnancy


Human epidermal growth factor receptor 2


Hazard ratio






Likelihood ratio


Odds ratio


Progesterone receptor


  1. Caleffi M, Fentiman IS, Birkhead BG. Factors at presentation influencing the prognosis in breast cancer. Eur J Cancer Clin Oncol. 1989;25(1):51–6.

    Article  CAS  Google Scholar 

  2. Trivers KF, Gammon MD, Abrahamson PE, Lund MJ, Flagg EW, Kaufman JS, et al. Association between reproductive factors and breast cancer survival in younger women. Breast Cancer Res Treat. 2007;103(1):93–102.

    Article  Google Scholar 

  3. Orgéas CC, Hall P, Rosenberg LU, Czene K. The influence of menstrual risk factors on tumor characteristics and survival in postmenopausal breast cancer. Breast Cancer Res. 2008;10(6):R107.

    Article  Google Scholar 

  4. Schouten LJ, Hupperets PS, Jager JJ, Volovics L, Wils JA, Verbeek AL, et al. Prognostic significance of etiological risk factors in early breast cancer. Breast Cancer Res Treat. 1997;43(3):217–23.

    Article  CAS  Google Scholar 

  5. Kroman N, Wohlfahrt J, West Andersen K, Mouridsen HT, Westergaard T, Melbye M. Parity, age at first childbirth and the prognosis of primary breast cancer. Br J Cancer. 1998;78(11):1529–33.

    Article  CAS  Google Scholar 

  6. Papatestas AE, Mulvihill M, Josi C, Ioannovich J, Lesnick G, Aufses AH. Parity and prognosis in breast cancer. Cancer. 1980;45(1):191–4.

    Article  CAS  Google Scholar 

  7. Green A, Beral V, Moser K. Mortality in women in relation to their childbearing history. Br Med J. 1988;297(6645):391–5.

    Article  CAS  Google Scholar 

  8. Korzeniowski S, Dyba T. Reproductive history and prognosis in patients with operable breast cancer. Cancer. 1994;74(5):1591–4.

    Article  CAS  Google Scholar 

  9. Black MM, Hankey BF, Barclay THC. Parity as a prognostic factor in young breast Cancer Patients2. J Natl Cancer Inst. 1983;70(1):27–30.

    CAS  PubMed  Google Scholar 

  10. Mohle-Boetani JC, Grosser S, Whittemore AS, Malec M, Kampert JB, Paffenbarger RS. Body size, reproductive factors, and breast cancer survival. Prev Med. 1988;17(5):634–42.

    Article  CAS  Google Scholar 

  11. Olson SH, Zauber AG, Tang J, Harlap S. Relation of time since last birth and parity to survival of young women with breast Cancer. Epidemiology. 1998;9(6):669–71.

    Article  CAS  Google Scholar 

  12. Phillips K-A, Milne RL, Friedlander ML, Jenkins MA, McCredie MRE, Giles GG, et al. Prognosis of premenopausal breast Cancer and childbirth prior to diagnosis. J Clin Oncol. 2004;22(4):699–705.

    Article  Google Scholar 

  13. Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7.

    Article  Google Scholar 

  14. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thürlimann B, et al. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol. 2013;24(9):2206–23.

    Article  CAS  Google Scholar 

  15. Senkus E, Kyriakides S, Penault-Llorca F, Poortmans P, Thompson A, Zackrisson S, et al. Primary breast cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2013;24(6):vi7–vi23.

    Article  Google Scholar 

  16. Millikan RC, Newman B, Tse C-K, Moorman PG, Conway K, Smith LV, et al. Epidemiology of basal-like breast cancer. Breast Cancer Res Treat. 2008;109(1):123–39.

    Article  Google Scholar 

  17. Jatoi I, Anderson WF. Qualitative age interactions in breast cancer studies: a mini-review. Future Oncol. 2010;6(11):1781–8.

    Article  Google Scholar 

  18. Yang XR, Chang-Claude J, Goode EL, Couch FJ, Nevanlinna H, Milne RL, et al. Associations of breast cancer risk factors with tumor subtypes: a pooled analysis from the breast Cancer association consortium studies. J Natl Cancer Inst. 2011;103(3):250–63.

    Article  Google Scholar 

  19. Anderson WF, Pfeiffer RM, Wohlfahrt J, Ejlertsen B, Jensen M-B, Kroman N. Associations of parity-related reproductive histories with ER± and HER2± receptor-specific breast cancer aetiology. Int J Epidemiol. 2017;46(1):373.

    Article  Google Scholar 

  20. Song N, Choi J-Y, Sung H, Jeon S, Chung S, Song M, et al. Tumor subtype-specific associations of hormone-related reproductive factors on breast Cancer survival. PLoS One. 2015;10(4):e0123994.

    Article  Google Scholar 

  21. Kwan ML, Bernard PS, Kroenke CH, Factor RE, Habel LA, Weltzien EK, et al. Breastfeeding, PAM50 Tumor Subtype, and Breast Cancer Prognosis and Survival. J Natl Cancer Inst. 2015;107(7):djv087–djv.

    Article  Google Scholar 

  22. Cespedes Feliciano EM, Kwan ML, Kushi LH, Chen WY, Weltzien EK, Castillo AL, et al. Body mass index, PAM50 subtype, recurrence, and survival among patients with nonmetastatic breast cancer. Cancer. 2017;123(13):2535–42.

    Article  CAS  Google Scholar 

  23. Yang XR, Devi BCR, Sung H, Guida J, Mucaki EJ, Xiao Y, et al. Prevalence and spectrum of germline rare variants in BRCA1/2 and PALB2 among breast cancer cases in Sarawak, Malaysia. Breast Cancer Res Treat. 2017;165(3):687–97.

    Article  CAS  Google Scholar 

  24. Devi CRB, Tang TS, Corbex M. Incidence and risk factors for breast cancer subtypes in three distinct south-east Asian ethnic groups: Chinese, Malay and natives of Sarawak, Malaysia. Int J Cancer. 2012;131(12):2869–77.

    Article  CAS  Google Scholar 

  25. Tamimi RM, Colditz GA, Hazra A, Baer HJ, Hankinson SE, Rosner B, et al. Traditional breast cancer risk factors in relation to molecular subtypes of breast cancer. Breast Cancer Res Treat. 2012;131(1):159–67.

    Article  CAS  Google Scholar 

  26. Mustapha A, Jenny C-C, Raza AH, Nilanjan C, Penny C, Frances D, et al. Etiology of hormone receptor positive breast cancer differs by levels of histologic grade and proliferation. Int J Cancer. 2018;143(4):746–57.

    Article  Google Scholar 

  27. Anderson KN, Schwab RB, Martinez ME. Reproductive risk factors and breast cancer subtypes: a review of the literature. Breast Cancer Res Treat. 2014;144(1):1–10.

    Article  Google Scholar 

  28. Chen L, Li CI, Tang M-TC, Porter P, Hill DA, Wiggins CL, et al. Reproductive factors and risk of luminal, HER2-overexpressing, and triple-negative breast Cancer among multiethnic women. Cancer epidemiology biomarkers &amp. Prevention. 2016;25(9):1297–304.

    CAS  Google Scholar 

  29. Kwan ML, Kushi LH, Weltzien E, Maring B, Kutner SE, Fulton RS, et al. Epidemiology of breast cancer subtypes in two prospective cohort studies of breast cancer survivors. Breast Cancer Res. 2009;11(3):1–13.

    Article  Google Scholar 

  30. Song N, Choi J-Y, Sung H, Chung S, Song M, Park SK, et al. Heterogeneity of epidemiological factors by breast tumor subtypes in Korean women: a case–case study. Int J Cancer. 2014;135(3):669–81.

    Article  CAS  Google Scholar 

  31. Martinez ME, Wertheim BC, Natarajan L, Schwab R, Bondy ML, Daneri-Navarro A, et al. Reproductive factors, heterogeneity, and breast tumor subtypes in women of Mexican descent. Cancer Epidemiol Biomarkers Prev. 2013;22(10):1853–61.

    Article  Google Scholar 

  32. Shinde SS, Forman MR, Kuerer HM, Yan K, Peintinger F, Hunt KK, et al. Higher parity and shorter breastfeeding duration. Cancer. 2010;116(21):4933–43.

    Article  Google Scholar 

  33. Wada K, Nagata C, Tamakoshi A, Matsuo K, Oze I, Wakai K, et al. Body mass index and breast cancer risk in Japan: a pooled analysis of eight population-based cohort studies. Ann Oncol. 2014;25(2):519–24.

    Article  CAS  Google Scholar 

  34. Bandera EV, Maskarinec G, Romieu I, John EM. Racial and ethnic disparities in the impact of obesity on breast Cancer risk and survival: a global perspective. Adv Nutr. 2015;6(6):803–19.

    Article  CAS  Google Scholar 

  35. Amadou A, Ferrari P, Muwonge R, Moskal A, Biessy C, Romieu I, et al. Overweight, obesity and risk of premenopausal breast cancer according to ethnicity: a systematic review and dose-response meta-analysis. Obes Rev. 2013;14(8):665–78.

    Article  CAS  Google Scholar 

  36. Clarke CA, Keegan THM, Yang J, Press DJ, Kurian AW, Patel AH, et al. Age-specific incidence of breast Cancer subtypes: understanding the black–white crossover. J Natl Cancer Inst. 2012;104(14):1094–101.

    Article  Google Scholar 

  37. Horne HN, Beena Devi CR, Sung H, Tang TS, Rosenberg PS, Hewitt SM, et al. Greater absolute risk for all subtypes of breast cancer in the US than Malaysia. Breast Cancer Res Treat. 2015;149(1):285–91.

    Article  Google Scholar 

  38. Juret P, Couette JE, Mandard AM, Carre A, Delozier T, Brune D, et al. Age at menarche as a prognostic factor in human breast cancer. Eur J Cancer. 1976;12(9):701–4.

    Article  CAS  Google Scholar 

  39. 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.

    Article  CAS  Google Scholar 

  40. Colditz GA. Relationship between estrogen levels, use of hormone replacement therapy, and breast Cancer. J Natl Cancer Inst. 1998;90(11):814–23.

    Article  CAS  Google Scholar 

  41. Althuis MD, Fergenbaum JH, Garcia-Closas M, Brinton LA, Madigan MP, Sherman ME. Etiology of hormone receptor–defined breast Cancer: a systematic review of the literature. Cancer Epidemiol Biomarkers Prev. 2004;13(10):1558–68.

    CAS  PubMed  Google Scholar 

  42. Alsaker MDK, Opdahl S, Asvold BO, Romundstad PR, Vatten LJ. The association of reproductive factors and breastfeeding with long term survival from breast cancer. Breast Cancer Res Treat. 2011;130(1):175–82.

    Article  Google Scholar 

  43. Urruticoechea A, Smith IE, Dowsett M. Proliferation marker Ki-67 in early breast cancer. J Clin Oncol. 2005;23(28):7212–20.

    Article  CAS  Google Scholar 

  44. Reed A, Kutasovic J, Lakhani S, Simpson P. Invasive lobular carcinoma of the breast: morphology, biomarkers and 'omics. Breast Cancer Res. 2015;17(1):12.

    Article  Google Scholar 

  45. Arpino G, Bardou V, Clark G, Elledge R. Infiltrating lobular carcinoma of the breast: tumor characteristics and clinical outcome. Breast Cancer Res. 2004;6(3):R149–R56.

    Article  Google Scholar 

  46. Chlebowski RT, Aiello E, McTiernan A. Weight loss in breast Cancer patient management. J Clin Oncol. 2002;20(4):1128–43.

    Article  Google Scholar 

  47. Chan DSM, Vieira AR, Aune D, Bandera EV, Greenwood DC, McTiernan A, et al. Body mass index and survival in women with breast cancer—systematic literature review and meta-analysis of 82 follow-up studies. Ann Oncol. 2014;25(10):1901–14.

    Article  CAS  Google Scholar 

  48. Nechuta S, Chen WY, Cai H, Poole EM, Kwan ML, Flatt SW, et al. A pooled analysis of post-diagnosis lifestyle factors in association with late estrogen-receptor–positive breast cancer prognosis. Int J Cancer. 2016;138(9):2088–97.

    Article  CAS  Google Scholar 

  49. Zhang M, Cai H, Bao P, Xu W, Qin G, Shu XO, et al. Body mass index, waist-to-hip ratio and late outcomes: a report from the shanghai breast Cancer survival study. Sci Rep. 2017;7(1):6996.

    Article  Google Scholar 

  50. Moon H-G, Han W, Noh D-Y. Underweight and breast Cancer recurrence and death: a report from the Korean breast Cancer society. J Clin Oncol. 2009;27(35):5899–905.

    Article  Google Scholar 

  51. Kawai M, Minami Y, Nishino Y, Fukamachi K, Ohuchi N, Kakugawa Y. Body mass index and survival after breast cancer diagnosis in Japanese women. BMC Cancer. 2012;12(1):149.

    Article  Google Scholar 

  52. Kawai M, Tomotaki A, Miyata H, Iwamoto T, Niikura N, Anan K, et al. Body mass index and survival after diagnosis of invasive breast cancer: a study based on the Japanese National Clinical Database—Breast Cancer Registry. Cancer Med. 2016;5(6):1328–40.

    Article  Google Scholar 

  53. Iyengar NM, Hudis CA, Dannenberg AJ. Obesity and inflammation: new insights into breast Cancer development and progression. Am Soc Clin Oncol Educ Book. 2013;33:46–51.

    Article  Google Scholar 

  54. Thomas F, Bean K, Pannier B, Oppert J-M, Guize L, Benetos A. Cardiovascular mortality in overweight subjects. The Key Role of Associated Risk Factors. 2005;46(4):654–9.

    CAS  Google Scholar 

Download references


The authors acknowledge Julie Buckland and Jane Demuth at Information Management Systems for data management support.


This research was supported by the Intramural Research Program of the National Institutes of Health, National Cancer Institute, Division of Cancer Epidemiology and Genetics.

Availability of data and materials

The datasets used and analyzed for the current study will be available from the corresponding author upon reasonable request.

Author information

Authors and Affiliations



MA, HS, and XRY conceived the study design and analytical concept. MA performed statistical analyses and drafted the manuscript. HS, RMP, and XRY contributed to the interpretation of results and the critical revision of the manuscript. DBCR and TST contributed to the acquisition of patient information and pathology data. JG helped with data management. All authors participated in the revision of the manuscript and approved the final manuscript.

Corresponding author

Correspondence to Xiaohong R. Yang.

Ethics declarations

Ethics approval and consent to participate

All study participants provided written informed consent and the project was approved by the Ethics Committee of the National Institutes of Health, Malaysia. The study was also exempted from review by the Office of Human Subject Research Protections at the National Institutes of Health since NIH investigators do not have the access to the personal identifying information.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Abubakar, M., Sung, H., BCR, D. et al. Breast cancer risk factors, survival and recurrence, and tumor molecular subtype: analysis of 3012 women from an indigenous Asian population. Breast Cancer Res 20, 114 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: