Skip to main content


Association of high obesity with PAM50 breast cancer intrinsic subtypes and gene expression

Article metrics



Invasive breast cancers are now commonly classified using gene expression into biologically and clinically distinct tumor subtypes. However, the role of obesity in breast tumor gene expression and intrinsic subtype is unknown.


Early-stage breast cancer (BC) patients (n = 1,676) were sampled from two prospective cohorts. The PAM50 qRT-PCR assay was used to: a) assess tumor gene expression levels for ESR1, PGR, ERBB2, and 10 proliferation genes and b) classify tumors into intrinsic subtype (Luminal A, Luminal B, Basal-like, HER2-enriched, Normal-like). Body mass index (BMI) around BC diagnosis (kg/m2) was categorized as: underweight (<18.5), normal (18.5-24), overweight (25–29), mildly obese (30–34), and highly obese (≥35). In a cross-sectional analysis, we evaluated associations of BMI with gene expression using linear regression models, and associations of BMI with non-Luminal A intrinsic subtypes, compared with Luminal A subtype, using multinomial logistic regression. Statistical significance tests were two-sided.


Highly obese women had tumors with higher expression of proliferation genes compared with normal weight women (adjusted mean difference = 0.44; 95% CI: 0.18, 0.71), yet mildly obese (adjusted mean difference = 0.16; 95% CI: −0.06, 0.38) and overweight (adjusted mean difference = 0.18; 95% CI: −0.01, 0.36) women did not. This association was stronger in postmenopausal women (p for interaction = 0.06). Being highly obese, however, was inversely associated with ESR1 expression (adjusted mean difference = −0.95; 95% CI: −1.47, −0.42) compared with being normal weight, whereas being mildly obese and overweight were not. In addition, women with Basal-like and Luminal B subtypes, relative to those with Luminal A subtype, were more likely to be highly obese, compared with normal-weight.


ER expression may not increase correspondingly with increasing degree of obesity. Highly obese patients are more likely to have tumor subtypes associated with high proliferation and poorer prognosis.


Invasive breast cancers are now commonly classified using gene expression into biologically and clinically distinct tumor subtypes known as Luminal A, Luminal B, Basal-like, and HER2-Enriched (HER2-E) [1,2]. Subtype information has been shown to be an independent predictor of breast cancer survival when used in multivariate analyses including standard clinicopathologic variables [3-6]. In 2009, Parker et al. derived a minimal gene set (PAM50) for classifying “intrinsic” subtypes of breast cancer [6,7]. The PAM50 gene set has high classification agreement with larger “intrinsic” gene sets previously used for subtyping [1,2,4,6], and is a feasible assay for application in clinical and epidemiologic studies that routinely use processed tumor tissue [8].

Substantial evidence suggests that obese women are at greater risk of postmenopausal breast cancer and have poorer breast cancer survival compared with normal-weight women [9-16]. Obesity might impact breast tumor development via increased estradiol production in adipose tissue in postmenopausal women, higher insulin levels, cellular interaction of leptin with insulin, and a constant pro-inflammatory state [17,18]. However, to our knowledge, the relationship between obesity before cancer diagnosis and likelihood of a specific tumor gene profile has not been examined. Therefore, in a cohort of 1,676 breast cancer survivors derived from two large prospective cohort studies, we explored cross-sectional associations of body mass index (BMI) around breast cancer diagnosis with PAM50-derived tumor expression of selected genes (ESR1, PGR, ERBB2, and proliferation) and intrinsic subtype.


Study population

The underlying study population was women from the LACE (PI: BJ Caan, [19]) and Pathways (PI: LH Kushi, [20]) prospective cohort studies of breast cancer survivors. A total of 2,135 LACE participants were 18–79 years old when diagnosed with early-stage breast cancer from 1997–2000 (AJCC stage I with tumor size ≥1 cm, stage II, or stage IIIA) and were identified primarily from the KPNC Cancer Registry (83%) or the Utah Cancer Registry (12%). Additional eligibility criteria included: being within 39 months of diagnosis to study enrollment (mean time = 23 months, 61% between 12 and 24 months), completion of chemotherapy or radiotherapy, no prior history of breast cancer or other cancer in the last 5 years.

The Pathways Study enrolled 4,505 women diagnosed with AJCC Stage I-IV breast cancer from 2006–2013 at KPNC with no previous diagnosis of other invasive cancer, at least 21 years of age at diagnosis, and spoke English, Spanish, or Chinese. Most women were approached for enrollment within two months of diagnosis (mean time = 1.8 months, range = 0.3-7.2 months).

Participants provided informed consent under human subjects’ protocols approved by the institutional review boards (IRB) at KPNC (CN-98BCaan-04-H) and the University of Utah (IRB_00038002). All human subjects’ research carried out in this study was in compliance with the Helsinki declaration (

Clinicopathologic characteristics

Clinicopathologic characteristics at cancer diagnosis, including disease stage, tumor size, nodal status, grade, estrogen receptor (ER) status, progesterone receptor (PR) status, and human epidermal growth factor receptor 2 (Her2) overexpression or amplification in the primary tumor, were abstracted from cancer registry data and medical record review.

Obesity and other covariates

Demographic and breast cancer risk factor data were collected at study enrollment on a mailed questionnaire (LACE) or in-person interview (Pathways),and included age at breast cancer diagnosis, race/ethnicity, education, menopausal status, smoking, and moderate-vigorous physical activity (metabolic equivalent (MET)-hours/week).

LACE women were asked to self-report their weight and height at 12 months before breast cancer diagnosis on the baseline questionnaire, which as completed on average 2 years post-diagnosis. Pathways women were asked to self-report their weight and height at the time of the baseline interview, conducted on average 2 months after breast cancer diagnosis. Then a BMI value which represents the period around breast cancer diagnosis was computed from these self-reported weight and height values, and categorized according to WHO international guidelines as a 5-level variable based on our previous work in obesity and breast cancer survival [21,22]. These categories (kg/m2) are: underweight (<18.5), normal (18.5-24), overweight (25–29), mildly obese (30–34), and highly obese (≥35).

Sampling strategy for PAM50 assay

For the PAM50 ancillary study (PI: BJ Caan), the LACE and Pathways cohorts were pooled with the overall goal to evaluate the performance of the PAM50 assay in a population-based study where patient characteristics, treatment patterns, and time of initial follow-up varied [23]. All LACE women from KPNC and Utah were eligible for the sub-study (n = 2,135), whereas Pathways women diagnosed from 2006–2008 were eligible (n = 2,172).

To further select eligible women for the PAM50 assay given limited study resources, we used a stratified case-cohort study design [24], with strata defined by clinical subtype based on immunohistochemistry (IHC) results for ER, PR, and Her2 [25]. The subcohort consisted of a random sample of women with the most common IHC subtype (ER+ or PR+, Her2-) (sampling fraction = 18%), and all women with the remaining less common subtypes having worse prognosis (ER+ or PR+, Her2+; ER-, PR-, Her2-; and ER-, PR-, Her2+) (sampling fraction = 100%). The cohort was followed for recurrence and survival through August 2013. Women who were not part of the subcohort but had an outcome of interest during this time were included. Out of 2,087 women selected for the case-cohort, 1,691 had tumor tissue successfully assayed by the PAM50. For this analysis, an additional 15 women were excluded due to missing BMI values, thus a total of n = 1,676 women comprised the final analytic sample with PAM50 data.

Tissue samples

For those selected into the case-cohort study, we obtained formalin-fixed, paraffin-embedded (FFPE) tissue blocks and corresponding slides from the surgical hospital or pathology storage facility. Slides were reviewed by one pathologist (R.E.F.). If the area of invasive tumor was observed to be smaller than 0.5 cm in diameter, the case was classified as ineligible. Tissue punches 1 mm in diameter were obtained from an area of the FFPE tissue block corresponding to the marked slide.

Gene expression assay and PAM50 intrinsic subtypes

Real-time reverse-transcription PCR (qRT-PCR) was conducted for the 50 target genes that comprise the PAM50 intrinsic subtype classifier [6], and details have been provided elsewhere [26,27]. Laboratory personnel (I.J.S.) were blinded to clinical information and received only a study identification number to track the sample.

The PAM50 assay yields an expression value for each gene that is relative to a reference gene. The raw data include both positive and negative values. For ease of interpretation of summary statistics, we transformed the values by adding 10 to all scores prior to analysis to make all values positive while preserving rank order. To determine intrinsic subtypes from the gene expression data, we applied centroid-based algorithms to the calibrated log-expression ratio for the 50 genes in the entire PAM50 assay. For each sample, this process generates five continuous-scale normalized subtype scores representing degree of Spearman correlation of gene expression with that of prototype Luminal A, Luminal B, Basal-like, HER2-E, and Normal-like breast tumors [6,26]. The subtype with the highest score became the predicted intrinsic subtype for that case.

Quantitative expression of individual genes that are standard breast cancer prognostic biomarkers (ESR1, PGR, ERBB2) were selected as variables of interest. Expression of 10 cell cycle regulation genes was averaged into a cell proliferation value (CENPF, ANLN, CDC20, CCNB1, CEP55, MYBL2, MKI67, UBE2C, RRM2, KIF2C).

Statistical analysis

All analyses incorporated sampling weights and the stratified sampling design for unbiased estimation of population parameters and valid estimates of standard errors [28,29]. This includes estimates of frequency distributions and chi-square tests of baseline characteristics and regression analyses using the ‘svy’ commands in Stata software, StataCorp, College Station, TX. Statistical significance tests were two-sided.

We described associations of obesity with intrinsic subtypes by fitting a multinomial logistic regression model. This method is similar to the case-case analysis approach widely used for dichotomous tumor characteristics [30], extended to the five subtype categories via the multinomial model. Treating the most prevalent subtype, Luminal A, as the base comparator outcome, we estimated odds ratios (OR) and 95% confidence intervals (CI) associated with BMI categories for each of the non-Luminal A subtypes. We also used multiple linear regression for point and interval estimation of adjusted differences in mean gene expression levels across BMI categories.

All models were adjusted for age at diagnosis, race/ethnicity, moderate-vigorous physical activity, and AJCC stage. Given that associations between obesity and breast cancer risk differ by menopausal status [31], models were also stratified by menopausal status (premenopausal vs. postmenopausal). In the non-stratified models, effect modification was evaluated by calculating p values for interaction via cross-product terms of BMI as a continuous variable and menopausal status.

All models were also run individually by cohort, and results were largely consistent with the combined cohort (Additional files 1 and 2), thus we present findings for the combined cohort below.


The distributions of demographic, clinical, and PAM50 intrinsic subtype by BMI are given in Table 1. Among the highly obese (≥35 kg/m2) and underweight (<18.5 kg/m2) groups, Basal-like subtypes were more common (18.9% and 22.6%, respectively), compared with the other BMI groups (8%-10%), whereas Luminal A was less common (36.3% and 33.9%, respectively), compared with the other groups (49%-56%). African American women were more likely to be highly obese (14.9%) and less likely to be normal weight, overweight, or mildly obese (3.7%-8.7%), whereas White women were more likely to be normal weight, overweight, or mildly obese (74.2%-75.1%) and less likely to be highly obese (59.4%). Both the highly obese and underweight groups had noticeably lower levels of moderate-vigorous physical activity (67.6% and 66.7% below median 18.9 MET-hours/week, respectively), in contrast to the other BMI groups that had higher (normal weight) or similar (overweight and mildly obese) levels relative to the median level. The normal-weight women had less comorbidities (9.1%) whereas the underweight women had more comorbidities (22.6%), compared with the other BMI groups (9.1%-18.2%).

Table 1 Demographic and clinical characteristics by BMI around breast cancer diagnosis, LACE and Pathways cohorts

ESR1 and proliferation unadjusted gene expression levels varied by BMI category, whereas there were no significant differences for PGR and ERBB2 levels (Table 2). Women who were highly obese or underweight had lower ESR1 gene expression (highly obese mean = 11.54, underweight mean = 10.95; p = 0.03) yet higher expression of proliferation genes (highly obese mean = 9.12, underweight mean = 9.07; p = 0.02), relative to women in the other BMI categories. When stratifying by menopausal status, lower ESR1 expression in the highly obese (mean = 10.78) and underweight (mean = 10.61) continued to be observed in premenopausal women only (p = 0.01), whereas higher proliferation expression was seen in postmenopausal women who were highly obese (mean = 9.09) but not underweight (mean = 8.40, p = 0.01).

Table 2 PAM50 gene expression levels by BMI around breast cancer diagnosis and menopausal status, LACE and Pathways cohorts

The adjusted mean differences in expression of ESR1, PGR, ERBB2, and proliferation genes by BMI category are given in Table 3. In models adjusted for age, race/ethnicity, moderate-vigorous physical activity, AJCC stage, and study, women who were highly obese (≥35 kg/m2) had tumors with 0.37 standard deviation (SD) higher expression of proliferation genes vs. normal-weight women (adjusted mean difference = 0.44; 95% CI: 0.18, 0.71), yet expression levels were fairly similar at 0.13 SD for the mildly obese (30.0-34.9 kg/m2) vs. normal-weight women (adjusted mean difference = 0.16; 95% CI: −0.06, 0.38). By contrast, the highly obese group had 0.34 SD lower expression of ESR1 compared with normal-weight women (adjusted mean difference = −0.95; 95% CI: −1.47, −0.42). Furthermore, compared with normal-weight women, the underweight women not only had 0.80 SD lower ESR1 expression (adjusted mean difference = −2.24; 95% CI: −3.81, −0.67 vs. 0), but 0.77 SD lower PGR expression as well (adjusted mean difference = −2.31; 95% CI: −3.98, −0.64). No associations of BMI with ERBB2 expression were observed.

Table 3 Adjusted mean difference in gene expression levels by BMI, overall and by menopausal status

When stratifying by menopausal status, the postmenopausal group had similar patterns of association to the overall cohort (Table 3). Highly obese women had tumors with 0.47 SD higher expression of proliferation genes compared with normal-weight women (adjusted mean difference = 0.54; 95% CI: 0.21, 0.86), yet mildly obese (adjusted mean difference = 0.17; 95% CI: −0.09, 0.42) and overweight (adjusted mean difference = 0.18; 95% CI: −0.04, 0.40) women did not (0.15 SD and 0.16 SD, respectively). In contrast, being highly obese was associated with 0.34 lower ESR1 expression (adjusted mean difference = −0.97; 95% CI: −1.61, −0.32) compared with normal weight, whereas being mildly obese (adjusted mean difference = 0.10; 95% CI: −0.38, 0.57) and overweight (adjusted mean difference = −0.10; 95% CI: −0.47, 0.26) were not (0.03 SD and 0.03 SD, respectively). Finally, being underweight compared with normal weight was associated with 0.83 SD lower PGR expression (adjusted mean difference = −2.46; 95% CI: −4.29, −0.62). We also examined expression levels in the very highly obese (≥40 kg/m2) postmenopausal women. Lower ESR1 expression and higher proliferation gene expression were observed at similar levels in both the highly obese (35–40 kg/m2) and very highly obese (>40 kg/m2) women compared with normal-weight women (data not shown).

In premenopausal women (Table 3), we also observed 0.37 SD lower, yet borderline significant, ESR1 expression among the highly obese (adjusted mean difference = −0.96; 95% CI: −1.97, 0.05), while ESR1 expression was 1.5 SD significantly lower among the underweight (adjusted mean difference = −3.90; 95% CI: −6.44, −1.36). There was no association between BMI and expression of proliferation genes. Effect modification of the proliferation associations by menopausal status was borderline statistically significant (p for interaction = 0.06).

The associations of BMI with intrinsic subtype are given in Table 4. In models adjusted for age, race/ethnicity, moderate-vigorous physical activity, AJCC stage, and study, breast cancer patients with Basal-like tumors had over triple the odds of being highly obese (≥35 kg/m2) vs. normal weight compared to those with Luminal A tumors (OR = 3.75; 95% CI: 1.97, 7.12). Women with Luminal B tumors also had increased odds of being highly obese compared to those with Luminal A tumors (OR = 2.44; 95% CI: 1.24, 4.79). There was little evidence for increased odds of being mildly obese in women with Basal-like (OR = 1.38; 95% CI: 0.81, 2.36) and Luminal B (OR = 1.14; 95% CI: 0.67, 1.95) tumors. While women with Basal-like and Luminal B tumors had increased odds of being overweight, the odds were of lower magnitude compared with being highly obese. Similarly, cases with Basal-like tumors had elevated odds of being underweight compared with those with Luminal A tumors (OR = 7.44; 95% CI: 2.03, 27.32). However, the number of underweight women was small (n = 13) and made the confidence interval wide for this association, thus limiting the interpretation of this result.

Table 4 Association of BMI around breast cancer diagnosis with PAM50 intrinsic subtype, overall and by menopausal status

In stratified analyses of BMI and intrinsic subtype by menopausal status (Table 4), the postmenopausal group had patterns of association similar to the overall cohort. These same relationships were not present in the premenopausal group, although no effect modification by menopausal status was observed (p for interaction = 0.52).


In this cohort study of 1,676 breast cancer survivors, we found that extreme obesity around breast cancer diagnosis was positively associated with poorer prognostic gene expression profiles and tumor subtypes. Highly obese women were more likely to have tumors with greater expression of proliferation genes, and lower expression of ESR1, which are characteristics of the Basal-like subtype. Correspondingly, compared to women with the Luminal A intrinsic subtype, those with the Basal-like subtype and Luminal B subtype had increased odds of being highly obese (≥35 kg/m2) around breast cancer diagnosis. We did not find comparable associations among women who were overweight or mildly obese. The association with proliferation gene expression was observed in postmenopausal but not premenopausal women. To our knowledge, ours is the first study to examine the association of different levels of obesity with tumor gene expression.

Interestingly, we observed that a very high BMI (≥35 kg/m2) around breast cancer diagnosis, while strongly and positively associated with both proliferation-related gene expression and associated intrinsic subtypes (Basal-like and Luminal B), was negatively associated with estrogen-related gene expression and associated intrinsic subtypes (Luminal A). Surprisingly, the latter finding is potentially inconsistent with the notion that among postmenopausal women, obesity is generally associated with higher plasma levels of estradiol from adipose tissue [31] and greater risk of ER+ (primarily Luminal A) breast cancer [32,33], and suggests that even if higher circulating estradiol is available in postmenopausal women, the tumor itself may be less responsive to endogenous estrogen depending on level of obesity. In a study of adipose gene expression and weight loss changes in postmenopausal women, greater weight loss was associated with borderline increased ESR1 expression (p for trend = 0.08) that could be possibly attributed to reduced adipose tissue inflammation [34]. Thus, perhaps at some threshold level of increasing obesity, tumor growth could be fueled by heightened inflammatory processes, rather than estrogen exposure, thus leading to decreased ESR1 expression and lower likelihood of developing the associated intrinsic subtypes. Finally, while we found that the associations for subtype might also be present in the overweight but not the mildly obese group, the magnitude of association was much smaller and could be due in part to chance.

Our observation of a possible threshold effect at high BMI being associated with a proliferative tumor gene expression profile, and development of Basal-like and Luminal B tumors, is consistent with our previous work which identified highly obese (≥40 kg/m2) breast cancer patients being at greatest risk for poorer prognosis and survival [22,35]. While the underlying biological mechanism linking obesity to tumor etiology is unclear, there are some intriguing and plausible hypotheses. Obesity can influence cancer risk by increased production of inflammatory factors, insulin and insulin-like growth factors (IGFs), and altered adipokines, resulting in a state of low-grade chronic inflammation [36]. Higher activity of the phosphatidylinositol 3-kinase (PI3K)-Akt pathway, which primarily regulates cellular proliferation, migration, and survival [37], has been well-described in Basal-like tumors [38-40]. Thus, perhaps higher levels of insulin and IGFs in the highly obese can drive the growth of Basal-like tumors, but not Luminal tumors, through this pathway [41,42]. Furthermore, higher circulating glucose levels in the highly obese could potentially support biosynthesis of Basal-like tumors, which have been shown to be more glycolytic than other tumor subtypes [43].

While we observed some intriguing associations among the underweight, including lower expression of ESR1 and PGR and being more likely to have a Basal-like tumor, we were limited by the small number of underweight women in our cohort (n = 13) to draw any definitive conclusions about this subgroup. To date, the role of underweight and tumor gene expression is largely unknown.

Strengths of this study include being the first to examine the relationship between obesity around breast cancer diagnosis and tumor gene expression, thus investigating potential molecular mechanisms of obesity on tumorigenesis. Given recent findings from large epidemiologic studies on high obesity and underweight, but not overweight or mild obesity being associated with poorer prognosis and survival [11,22,35], we were also able to examine the association of varying degrees of obesity in relation to subtype and gene expression. Finally, we used the PAM50 assay, which is a classification tool that has been shown to have better prognostic ability than surrogate IHC classification methods [23,26].

Several limitations should be noted. Weight and height were self-reported, yet substantial agreement between BMI based on self-reported, compared with measured, weight and height has been shown [44]. Also, the number of underweight women in our cohort was small (n = 13). However, we chose to keep the underweight group as a separate category, as we have previously observed elevated risks of breast cancer mortality in underweight women [22,35] and considered this analysis of obesity and gene expression exploratory. This was a cross-sectional analysis, thus causality could not be inferred between obesity around breast cancer diagnosis and tumor intrinsic subtype and gene expression. In addition, BMI reflects the relationship of weight to height and thus does not reflect between-individual variation in total adiposity [45] and metabolic risk profiles [46].

One should also consider that perhaps it is not weight at one timepoint around breast cancer diagnosis, but rather weight trajectories over the life course that may act on gene expression or tumor subtype [47]. Furthermore, compared with non-Hispanic Whites, African Americans and Hispanics are more likely to be obese [48,49], and African Americans are more likely to be diagnosed with poor prognosis subtypes [25,50,51]. Given these associations, as a sensitivity analysis, we restricted the statistical models to Whites only, and both subtype and gene expression results were essentially unchanged. Finally, LACE women were enrolled on average two years post-diagnosis, thus women with better prognosis subtypes (Luminal A) could have been more likely to survive to enrollment whereas those with poorer prognosis subtypes (Basal-like) were not. However, we found in other analyses that this potential survival bias was minimal [23], and again when we restricted the analyses by individual cohort, the results were similar (Additional files 1 and 2).


Among women with breast cancer, particularly postmenopausal women, those who were highly obese, but not mildly obese, around breast cancer diagnosis were more likely to have breast tumors with greater expression of proliferation genes and lesser expression of ESR1, and possible increased odds of being diagnosed with Basal-like and Luminal B tumor subtypes. These findings suggest that etiology of tumor subtypes may vary by degree of pre-existing obesity of the patient and propose novel insights into molecular mechanisms linking obesity, ER expression, and proliferation to breast tumor development.



Life After Cancer Epidemiology

ESR1 :

Estrogen receptor 1


Progesterone receptor


Epidermal growth factor receptor


Estrogen receptor


Progesterone receptor


Hazard ratio


Confidence interval


Luminal A


Luminal B








Real-time reverse-transcription PCR


  1. 1.

    Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.

  2. 2.

    Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001;98(19):10869–74.

  3. 3.

    Cheang MC, Voduc KD, Tu D, Jiang S, Leung S, Chia SK, et al. Responsiveness of intrinsic subtypes to adjuvant anthracycline substitution in the NCIC.CTG MA.5 randomized trial. Clin Cancer Res. 2012;18(8):2402–12.

  4. 4.

    Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, et al. The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006;7:96.

  5. 5.

    Nielsen TO, Parker JS, Leung S, Voduc D, Ebbert M, Vickery T, et al. A comparison of PAM50 intrinsic subtyping with immunohistochemistry and clinical prognostic factors in tamoxifen-treated estrogen receptor-positive breast cancer. Clin Cancer Res. 2010;16(21):5222–32.

  6. 6.

    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.

  7. 7.

    Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A. 2002;99(10):6567–72.

  8. 8.

    Sweeney C, Bernard P, Factor RE, Kwan ML, Habel LA, Quesenberry Jr CP, et al. Intrinsic subtypes from PAM50 gene expression assay in a population-based breast cancer cohort: Differences by age, race, and tumor characteristics. Cancer Epidemiol Biomarkers Prev. 2014;23(5):714–24.

  9. 9.

    Goodwin PJ, Boyd NF. Body size and breast cancer prognosis: a critical review of the evidence. Breast Cancer Res Treat. 1990;16(3):205–14.

  10. 10.

    Ryu SY, Kim CB, Nam CM, Park JK, Kim KS, Park J, et al. Is body mass index the prognostic factor in breast cancer?: a meta-analysis. J Korean Med Sci. 2001;16(5):610–4.

  11. 11.

    Conroy SM, Maskarinec G, Wilkens LR, White KK, Henderson BE, Kolonel LN. Obesity and breast cancer survival in ethnically diverse postmenopausal women: the Multiethnic Cohort Study. Breast Cancer Res Treat. 2011;129(2):565–74.

  12. 12.

    Ewertz M, Jensen MB, Gunnarsdottir KA, Hojris I, Jakobsen EH, Nielsen D, et al. Effect of obesity on prognosis after early-stage breast cancer. J Clin Oncol. 2011;29(1):25–31.

  13. 13.

    Protani M, Coory M, Martin JH. Effect of obesity on survival of women with breast cancer: systematic review and meta-analysis. Breast Cancer Res Treat. 2010;123(3):627–35.

  14. 14.

    Chen X, Lu W, Zheng W, Gu K, Chen Z, Zheng Y, et al. Obesity and weight change in relation to breast cancer survival. Breast Cancer Res Treat. 2010;122(3):823–33.

  15. 15.

    Patterson RE, Cadmus LA, Emond JA, Pierce JP. Physical activity, diet, adiposity and female breast cancer prognosis: a review of the epidemiologic literature. Maturitas. 2010;66(1):5–15.

  16. 16.

    Connolly BS, Barnett C, Vogt KN, Li T, Stone J, Boyd NF. A meta-analysis of published literature on waist-to-hip ratio and risk of breast cancer. Nutr Cancer. 2002;44(2):127–38.

  17. 17.

    Hursting SD, Lashinger LM, Wheatley KW, Rogers CJ, Colbert LH, Nunez NP, et al. Reducing the weight of cancer: mechanistic targets for breaking the obesity-carcinogenesis link. Best Pract Res Clin Endocrinol Metab. 2008;22(4):659–69.

  18. 18.

    Renehan AG, Roberts DL, Dive C. Obesity and cancer: pathophysiological and biological mechanisms. Arch Physiol Biochem. 2008;114(1):71–83.

  19. 19.

    Caan B, Sternfeld B, Gunderson E, Coates A, Quesenberry C, Slattery ML. Life After Cancer Epidemiology (LACE) Study: a cohort of early stage breast cancer survivors (United States). Cancer Causes Control. 2005;16(5):545–56.

  20. 20.

    Kwan ML, Ambrosone CB, Lee MM, Barlow J, Krathwohl SE, Ergas IJ, et al. The Pathways Study: a prospective study of breast cancer survivorship within Kaiser Permanente Northern California. Cancer Causes Control. 2008;19(10):1065–76.

  21. 21.

    Caan BJ, Kwan ML, Shu XO, Pierce JP, Patterson RE, Nechuta SJ, et al. Weight change and survival after breast cancer in the after breast cancer pooling project. Cancer Epidemiol Biomarkers Prev. 2012;21(8):1260–71.

  22. 22.

    Kwan ML, Chen WY, Kroenke CH, Weltzien EK, Beasley JM, Nechuta SJ, et al. Pre-diagnosis body mass index and survival after breast cancer in the After Breast Cancer Pooling Project. Breast Cancer Res Treat. 2012;132(2):729–39.

  23. 23.

    Caan BJ, Sweeney C, Habel LA, Kwan ML, Kroenke CH, Weltzien E, et al. Intrinsic subtypes from the PAM50 gene expression assay in a population-based breast cancer survivor cohort: Prognostication of short and long term outcomes. Cancer Epidemiol Biomarkers Prev. 2014;23(5):725–34.

  24. 24.

    Wacholder S. Practical considerations in choosing between the case-cohort and nested case–control designs. Epidemiology. 1991;2(2):155–8.

  25. 25.

    Carey LA, Perou CM, Livasy CA, Dressler LG, Cowan D, Conway K, et al. Race, breast cancer subtypes, and survival in the Carolina Breast Cancer Study. Jama. 2006;295(21):2492–502.

  26. 26.

    Bastien RR, Rodriguez-Lescure A, Ebbert MT, Prat A, Munarriz B, Rowe L, et al. PAM50 breast cancer subtyping by RT-qPCR and concordance with standard clinical molecular markers. BMC Med Genomics. 2012;5:44.

  27. 27.

    Ebbert MT, Bastien RR, Boucher KM, Martin M, Carrasco E, Caballero R, et al. Characterization of uncertainty in the classification of multivariate assays: application to PAM50 centroid-based genomic predictors for breast cancer treatment plans. J Clin Bioinforma. 2011;1:37.

  28. 28.

    Chambers RL, Skinner CJ. Analysis of survey data. Chichester, U.K.: Wiley; 2003.

  29. 29.

    Cochran WG. Sampling techniques. 3rd ed. New York: John Wiley & Sons; 1977.

  30. 30.

    Begg CB, Zhang ZF. Statistical analysis of molecular epidemiology studies employing case-series. Cancer Epidemiol Biomarkers Prev. 1994;3(2):173–5.

  31. 31.

    Key TJ, Appleby PN, Reeves GK, Roddam A, Dorgan JF, Longcope C, et al. Body mass index, serum sex hormones, and breast cancer risk in postmenopausal women. J Natl Cancer Inst. 2003;95(16):1218–26.

  32. 32.

    Canchola AJ, Anton-Culver H, Bernstein L, Clarke CA, Henderson K, Ma H, et al. Body size and the risk of postmenopausal breast cancer subtypes in the California Teachers Study cohort. Cancer Causes Control. 2012;23:473–85.

  33. 33.

    Munsell MF, Sprague BL, Berry DA, Chisholm G, Trentham-Dietz A. Body mass index and breast cancer risk according to postmenopausal estrogen-progestin use and hormone receptor status. Epidemiol Rev. 2014;36(1):114–36.

  34. 34.

    Campbell KL, Foster-Schubert KE, Makar KW, Kratz M, Hagman D, Schur EA, et al. Gene expression changes in adipose tissue with diet- and/or exercise-induced weight loss. Cancer Prev Res (Phila). 2013;6(3):217–31.

  35. 35.

    Kwan ML, John EM, Caan BJ, Lee VS, Bernstein L, Cheng I, et al. Obesity and mortality after breast cancer by race/ethnicity: the california breast cancer survivorship consortium. Am J Epidemiol. 2014;179(1):95–111.

  36. 36.

    Lumeng CN, Saltiel AR. Inflammatory links between obesity and metabolic disease. J Clin Invest. 2011;121(6):2111–7.

  37. 37.

    Bader AG, Kang S, Zhao L, Vogt PK. Oncogenic PI3K deregulates transcription and translation. Nat Rev Cancer. 2005;5(12):921–9.

  38. 38.

    Hoeflich KP, O’Brien C, Boyd Z, Cavet G, Guerrero S, Jung K, et al. In vivo antitumor activity of MEK and phosphatidylinositol 3-kinase inhibitors in basal-like breast cancer models. Clin Cancer Res. 2009;15(14):4649–64.

  39. 39.

    Marty B, Maire V, Gravier E, Rigaill G, Vincent-Salomon A, Kappler M, et al. Frequent PTEN genomic alterations and activated phosphatidylinositol 3-kinase pathway in basal-like breast cancer cells. Breast Cancer Res. 2008;10(6):R101.

  40. 40.

    Moulder SL. Does the PI3K pathway play a role in basal breast cancer? Clin Breast Cancer. 2010;10 Suppl 3:S66–71.

  41. 41.

    Lopez-Knowles E, O’Toole SA, McNeil CM, Millar EK, Qiu MR, Crea P, et al. PI3K pathway activation in breast cancer is associated with the basal-like phenotype and cancer-specific mortality. Int J Cancer. 2010;126(5):1121–31.

  42. 42.

    Moestue SA, Dam CG, Gorad SS, Kristian A, Bofin A, Maelandsmo GM, et al. Metabolic biomarkers for response to PI3K inhibition in basal-like breast cancer. Breast Cancer Res. 2013;15(1):R16.

  43. 43.

    Palaskas N, Larson SM, Schultz N, Komisopoulou E, Wong J, Rohle D, et al. 18 F-fluorodeoxy-glucose positron emission tomography marks MYC-overexpressing human basal-like breast cancers. Cancer Res. 2011;71(15):5164–74.

  44. 44.

    Craig BM, Adams AK. Accuracy of body mass index categories based on self-reported height and weight among women in the United States. Matern Child Health J. 2009;13(4):489–96.

  45. 45.

    Keys A, Fidanza F, Karvonen MJ, Kimura N, Taylor HL. Indices of relative weight and obesity. J Chronic Dis. 1972;25(6):329–43.

  46. 46.

    Florez H, Castillo-Florez S. Beyond the obesity paradox in diabetes: fitness, fatness, and mortality. Jama. 2012;308(6):619–20.

  47. 47.

    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.

  48. 48.

    Flegal KM, Carroll MD, Kit BK, Ogden CL. Prevalence of obesity and trends in the distribution of body mass index among US adults, 1999–2010. Jama. 2012;307(5):491–7.

  49. 49.

    Stommel M, Schoenborn CA. Variations in BMI and prevalence of health risks in diverse racial and ethnic populations. Obesity (Silver Spring). 2010;18(9):1821–6.

  50. 50.

    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):R31.

  51. 51.

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

Download references


This work was supported by the National Institutes of Health (R01 CA129059 to B.J.C. and R01 CA105274 to L.H.K.). Additional support was from the Bioinformatics and Biostatistics core resources of the Huntsman Cancer Institute (P30 CA042014). The Utah Cancer Registry is funded by Contract No. HHSN261201000026C from the National Cancer Institute’s SEER Program along with additional support from the Utah State Department of Health and the University of Utah.

We thank Bryan M. Langholz, Ph.D. at University of Southern California for biostatistical consultation on case-cohort methodology, and Ms. Carole Davis at University of Utah for laboratory assay support.

The content of this manuscript is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health.

Author information

Correspondence to Marilyn L Kwan.

Additional information

Competing interests

P.S.B. is named on the patent for PAM50 which is licensed to Bioclassifier LLC. The other authors declare that they have no competing interests.

Authors’ contributions

MLK contributed to conception and design, analysis, and interpretation of data, and drafted the manuscript. CHK contributed to analysis and interpretation of data. CS contributed to analysis and interpretation of data. PSB contributed to laboratory data acquisition and interpretation of data. EKW carried out the statistical analysis and interpretation of data. AC contributed to cohort data acquisition and interpretation of data. REF contributed to laboratory data acquisition and interpretation of data. KSM contributed to cohort data acquisition. IJS contributed to laboratory data acquisition. CPQ contributed to conception and design, and the statistical analysis. LAH contributed to analysis and interpretation of data. LHK had primary responsibility for study conception and data acquisition for the Pathways cohort, and interpretation of data for this manuscript. BJC had primary responsibility for study conception and data acquisition for the LACE cohort and for conception of the intrinsic subtype study. She contributed to analysis and interpretation of data for this manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1:

LACE cohort results.

Additional file 2:

Pathways cohort results.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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

Verify currency and authenticity via CrossMark

Cite this article

Kwan, M.L., Kroenke, C.H., Sweeney, C. et al. Association of high obesity with PAM50 breast cancer intrinsic subtypes and gene expression. BMC Cancer 15, 278 (2015) doi:10.1186/s12885-015-1263-4

Download citation


  • Breast cancer
  • Obesity
  • Body mass index
  • PAM50 intrinsic subtype classifier
  • Tumor subtype
  • Gene expression
  • ESR1
  • PGR
  • ERBB2
  • Proliferation