- Research article
- Open Access
Analysis of high fat diet induced genes during mammary gland development: identifying role players in poor prognosis of breast cancer
BMC Research Notesvolume 7, Article number: 543 (2014)
Epidemiological studies have shown that consumption of a high-fat diet (HFD) increases the risk of developing breast cancer (BC). Studies in rodents have shown HFD causes changes in the genetic programming of the maturing mammary gland (MG) increasing the susceptibility of developing the disease. Less is known about how HFD induced genes impact BC development. HFD exposure two weeks before conception to six weeks of age was previously shown to dramatically change MG gene expression in 10 week old mice. Therefore, we investigated these differentially expressed HFD-induced genes for their expression in BC using the NKI 295 breast tumor dataset.
To examine the potential role of HFD induced genes in BC, we first investigated whether these HFD-induced genes in mouse MGs were differentially expressed in different types of human BC. Of the 28 HFD induced genes that were differentially expressed between BC subtypes in the NKI set, 79% were significantly higher in basal-like BC. Next, we analyzed whether HFD induced genes were associated with BC prognosis utilizing gene expression and survival data for each HFD induced gene from the NKI data and constructed Kaplan Meier survival plots. Significantly, 93% of the prognosis associated genes (13/14) were associated with poor prognosis (P = 0.002). Kaplan Meier analysis with 249 non-basal-like BC found that all but one of the genes examined were still significantly associated with poor prognosis. Furthermore, gene set enrichment analysis (GSEA) with HFD microarray data revealed that invasive BC genes where enriched in HFD samples that also had lost expression of luminal genes.
HFD exposed mouse MGs maintain differential expression of genes that are found highly expressed in basal-like breast cancer. These HFD-induced genes associate with poor survival in numerous BC subtypes, making them more likely to directly impact prognosis. Furthermore, HFD exposure leads to a loss in the expression of luminal genes and a gain in expression of mesenchymal and BC invasion genes in MGs. Collectively, our study suggests that HFD exposure during development induces genes associated with poor prognosis, thus identifying how HFD diet may regulate BC development.
Recent evidence suggests that dietary content is one causal lifestyle factor that may contribute to breast cancer (BC) development. Studies focusing on migrant populations have shown that women migrating from low (Asian, Latin American) to high BC incidence rate countries (United States and other Western countries) acquire a higher risk for BC similar to those in the new country [1–6]. These studies suggest that dietary changes might be a strong contributing factor to the increased incidence rates observed in these migrating populations. A study analyzing the effects of dietary patterns and the risk of BC in women of different ethnicities concluded that women consuming a Western diet had a higher risk of BC . The Western diet consists of foods primarily high in fat. A higher intake of dietary fat has been shown to increase BC risk in adults [8, 9].
A high-fat diet (HFD) may also determine the type of BC that will develop translating to a more precise treatment strategy increasing patient survival. Breast cancer tumors have been previously classified into subtypes based on their unique molecular/genetic expression profiles that correlate with phenotypic characteristics and clinical outcome [10–13]. The expression of estrogen receptor (ER) is a distinguishing marker between these tumors subtypes in combination with other molecular cues. Luminal subtype A (luminal A) and luminal subtype B (luminal B) are tumors that have high expression levels of ER (ER+) and are associated with a favorable prognosis. ERBB2 is overexpressed in HER2 positive breast cancers whereas basal epithelial-like (basal-like) tumors are ER negative (ER-). HER2 and basal-like tumors are associated with a poor prognosis [10, 11, 13, 14]. Strikingly, the five year survival rate for patients with luminal A is 90% compared to as low as a 30% five year survival rate for patients with HER2 and basal-like BC tumors [14–16]. This difference is due, in part, to the difficulty in predicting a clinical course, tumor stage at diagnosis and availability of targeted therapy [14, 17]. Obesity significantly increases the risk of developing basal-like BC in premenopausal women, predicting a poor outcome in these individuals . Commensurately, obese individuals are much less likely to develop luminal BC. Therefore, identifying the impact a HFD has on the type of BC a women develops may be critical in identifying the disease early to help predict a clinical course of action that would increase that patient’s survival.
The impact that HFD has on developing BC in humans has yet to be fully elucidated, in part, due to the limited number of human studies (and ability) to adequately assess the effect diet has on genetic regulation of the developing mammary gland (MG). Studies using rodent models have provided the most compelling evidence linking developmental exposure to a HFD to breast carcinogenesis [19–21]. These studies culminate to identify several theories to how consumption of a HFD during pre and/or postnatal development may cause BC including genetic reprogramming of the MG [22–25]. However, it is unclear as to how HFD might contribute to the development of specific types of BC correlating to the severity of the disease and patient survival. Thus, identifying MG HFD-induced transcriptional programs that are pertinent to specific types of human BC could have a tremendous impact on understanding the etiology of diet induced BC.
The purpose of this study was to examine genes induced by exposure to a HFD during MG development in the mouse for the role they may play in human BC. Luitjen et. al. previously identified genes that were differentially expressed in 10 week old MGs after developmental HFD exposure. Strikingly, these changes were sustained even after HFD exposure, suggesting a long-term effect. To investigate how HFD may influence the development of BC, we examined the experimentally identified HFD-induced genes from the study by Luitjen et al. in human breast tumors. The expression of HFD- induced genes was examined in BC by utilizing the widely studied NKI BC dataset from the Netherlands Cancer Institute, which contains gene expression and patient follow up/survival data from 295 women between the ages of 26 and 53 . The NKI BC data set had been previously classified into five breast cancer subtypes . Using Kaplan Meier analysis and gene set enrichment analysis (GSEA) we found that HFD induced genes were associated with poor prognosis and invasive BC. We also found that HFD leads to a loss in luminal gene expression and a gain in basal-like gene expression. Our data shows for the first time that HFD-induced genes are highly expressed in BCs that are of the basal-like subtype and are associated with poor prognosis.
Collection of HFD-induced genes and gene expression analysis in BC subtypes
Microarray data of differentially expressed genes was obtained from a study analyzing the effects of early diet on the genetic programming of the MG in wild-type mice . As described previously , FVB wild-type mice were fed a HFD diet consisting of either a 24% high fat content of either n-6 (corn oils) or n-3 (flaxseed oil) PUFAs (polyunsaturated fatty acids) from two weeks prior conception to 6 weeks of age. Control animals were fed standard rodent chow (5% fat) during the same period. Gene expression profiling of the MG was performed at 10 weeks of age. Luitjen et al. utilized the data from both types of fat (n-3 and n-6) in ANOVA analysis to identify differentially expressed genes. We examined the top one hundred differentially expressed genes that were induced by HFD with a p-value of 0.0001 or less based on the Luitjen ANOVA analysis for representation in the NKI dataset. The NKI dataset contains gene expression and survival data for 295 human breast tumors; these samples have also been previously classified based on breast cancer subtypes (basal-like = 46 tumors; HER2 positive = 49 tumors; luminal A = 88 tumors; luminal B = 81 tumors; normal-like = 31 tumors) [10, 11]. Our analysis focused on the HFD-induced genes found in the NKI data set (41 total). The average expression of each HFD-induced gene was determined for each BC sub-type. In order to determine whether HFD genes were differentially expressed between breast cancer subtypes, we performed the two tiered, type two T-test to determine whether each gene was expressed signficantly different between the basal-like subtype and the other sub-types of BC (HER2, luminal A, luminal B and normal-like). The T-test gave the same p-value as one way ANOVA, however a T-test anaylsis was used given we were only examining a single factor between two groups at a time. The false discovery rate (FDR) was calculated for differentially expressed genes to correct for multiple hypothesis testing. Changes in gene expression between subtypes that had an FDR of 0.01 or less were considered to be signifcant. It is improtant to note that the data sets utilized in this article were from experiemnts previously approved by their respective institutial ethics commitees and follow their national legeslation [10, 25].
A Kaplan Meier survival analysis for each HFD induced gene was performed using the NKI dataset and MedCalc software (Mariakerke, Belgium) for each gene to be analyzed. The association with prognosis of human breast cancer for each gene was obtained from this analysis: poor prognosis; good prognosis; no association. 295 breast tumors were divided into two equal-sized groups based on the expression of each gene; group one represented low expression of the gene and group two represented high expression of the gene. All availabe time points were utilized. P-values of less than 0.05 were considered significant associations. Based on the Chang et al. classification of the NKI database, we utilized 249 breast cancer samples for the non-basal like breast tumor analysis.
Additional statistical analysis of poor prognosis genes
HFD-induced genes that were associated with poor prognosis were analyzed as a group for overall differential expression between BC subtypes using NKI data. We performed ONE-WAY ANOVA with a Dunnett’s post test to determine the significance of mean gene expression between BC subtypes using the basal-like subtype as the reference group. Probability values of less than 0.05 were considered significant. The PC IBM INSTAT and PRISM software were used for the statistical analyses and graphs, respectively (GraphPad, San Diego, CA, USA).
Gene set enrichment analysis
Gene set enrichment analysis (GSEA) was performed using the Broad Institute platform [26, 27]. Samples were analyzed with weighted, Signal2Noise default settings. Microarray gene expression data from Luitjen et al. was grouped based on diet (control low fat diet or HFD samples) . Both types of HFD (n-3 and n-6 as indicated) were independently studied by GSEA with similar results. (Additional file 1: Table S1) contains the list of 361 breast cancer gene sets from the Molecular Signatures database, Broad Institute, that were utilized for analysis. The luminal gene sets utilized in GSEA are in (Additional file 2: Table S2).
HFD during MG development induces basal-like BC genes
Murine MGs exposed to HFD during development (two weeks before conception to 6 weeks post birth) had dramatic changes in gene expression at ten weeks of age . In order to investigate how HFD may influence BC development, we examined genes from the Luitjen et al. study that were induced by HFD with a p-value of 0.0001 or less for expression in BC subtypes by utilizing data for these genes in the NKI BC dataset. Gene expression signatures classify breast cancer into 5 basic subtypes: basal-like, HER2, luminal A, luminal B and normal-like. The basal-like subtype is associated with the worst prognosis whereas the luminal A is associated with the best prognosis [12, 13]. The average expression of each HFD induced gene was determined for each BC subtype (Table 1). Changes with an FDR of 0.01 or less were considered significantly different. Twenty-eight out of the 41 HFD induced genes showed statistically different levels of gene expression between the luminal A (associated with good prognosis) and basal-like (associated with poor prognosis) subtypes of BC; 79% of these genes were significantly higher in basal-like BC compared to luminal A BC. Therefore, HFD exposed glands preferentially maintain induced expression of genes that are found highly expressed in basal-like BC.
HFD induced genes are strongly associated with poor prognosis in BC
In order to examine whether HFD induced genes were associated with BC prognosis, we performed Kaplan Meier survival analysis using HFD induced genes from 10-week old mammary fat pads and the NKI dataset, as described in the Methods. A total of fourteen genes out of forty-one were associated with prognosis in BC (p < 0.05). Only one HFD-induced gene was associated with good prognosis; the high expression of Esterase D (ESD), showed a higher survival probability when compared to the low gene expression (Figure 1, Table 1). Significantly, 93% of the prognosis associated genes (13 out of 14 HFD-induced genes) were associated with poor prognosis; the high expression of these genes showed a lower survival probability when compared to their low expression (Figure 1, Table 1) (P = 0.002). The AmiGO database  was utilized to identify gene ontologies that were associated with the HFD-induced poor prognosis genes. We found that the poor prognosis genes were associated with cell adhesion, proteolysis, immunological response, ion transport and DNA methylation (Table 2). Altogether, HFD-induced genes during mammary development are highly associated with poor survival in BC.
The strong association between HFD-induced genes and poor prognosis in human breast cancer prompted us to further investigate this phenomena. We first examined the expression distribution of these genes in the 295 breast cancer tumors by plotting the expression of individual genes for each breast tumor sample, (Additional file 3: Figure S1) and (Additional file 4: Figure S2). One can see that there is a high expression of each poor prognosis associated gene in basal-like breast cancer samples whereas there is lower expression in luminal A breast cancers. Each of the poor prognosis genes is significantly less well expressed in luminal BC compared to basal-like BC, Table 1. Genes associated with poor prognosis were analyzed as a group for differential expression in human breast cancer subtypes. The average expression of these genes showed significant statistical difference in HER2, luminal A, luminal B, and normal-like subtypes with respect to the basal-like subtype (ANOVA p < 0.01, Figure 2). Of note, the differences in expression were more striking between basal-like BC samples and luminal A, luminal B and normal-like tumors, whereas a less dramatic, but still significant difference was found between basal-like tumors and HER2 positive tumors. This result is consistent with the fact that HER2 positive cancer are associated with poor prognosis as well [12, 14].
Although it is clear that HFD-induced genes are more strongly expressed in basal-like breast cancer, we were interested to examine whether these genes were associated with prognosis in other subtypes of breast cancer. To investigate this, we performed Kaplan Meier analysis with 249 non-basal-like NKI samples and found that all but 1 (92%) of the genes examined (Matrix metalloproteinase 12, MMP12) were still significantly associated with poor prognosis (Table 3). Therefore, HFD-induced genes associate with poor survival in a broad range of human breast cancer subtypes, making them likely to directly impact prognosis and have an influence on survival in many types of breast cancer.
HFD exposure induces BC mesenchymal and invasion genes
HFD exposure during development in mice leads to dramatic changes in mammary fat pad gene expression. These changes strongly suggest that the breast tissue exposed to HFD are developmentally altered and may sustainably express programs that parallel those found in certain types of BC. To investigate this, we performed GSEA with 361 curated BC gene sets from the Broad Institute Molecular Signatures Database. Each of the 361 examined gene sets contained genes that were previously determined to be differentially expressed under defined condition in breast cancer. In GSEA, Kolmogorov-Smirnov-style statistics are employed to investigate whether gene sets are more highly expressed (enriched) between two conditions in microarray analyses. We compared microarray gene expression data from control low fat MGs to data from HFD MGs to determine whether any of the BC gene sets were significantly enriched in HFD samples. We found that genes associated with BC invasion were significantly enriched in HFD samples, Figure 3 and Table 4. Enrichment plots are depicted in Figure 3 in which genes are ranked based on their association with given phenotypes, in this case the phenotypes are HFD (on the right) or control diet (on the left). The lines underneath the enrichment plot depict genes that are contained on the gene list being investigated. One can see a striking enrichment of invasion genes in HFD treated samples, Figure 3. Therefore, HFD leads to an increase in the expression of genes involved in BC invasion.
To further examine HFD mediated changes in gene expression in the MG with regard to BC development, we performed GSEA with 53 luminal BC gene sets from the Molecular Signatures database, Broad Institute. We found that genes that are normally down-regulated in the luminal subtype of BC (and are highly expressed in mesenchymal BC) were enriched in HFD treated samples, suggesting that these mice may have altered breast development to a more basal phenotype (Figure 4 A-B, Table 5). Enrichment plots for analysis with gene sets that are down in luminal BC and high in mesenchymal BC (Charafe_luminal_versus_mesenchymal_DN) indicate a clear increase in mesenchymal BC genes upon HFD treatment, Figure 4. Furthermore, we found that stromal genes from poor prognosis BC (gene-set FINAK_BREAST_CANCER_SDPP_SIGNATURE, SDPP is an abbreviation for stroma-derived prognostic predictor of BC disease outcome), were strongly associated with HFD MGs with a Normalized Enrichment Score (NES) of 1.87, p = <0.001 (Table 4). One of the top enriched genes in our GSEA was the mesenchymal marker Vimentin. We assessed the average Vimentin expression in control MGs versus HFD fed MGs using the Luitjen et al. microarray dataset and found that HFD exposed MGs expressed significantly higher levels of the basal-like BC marker Vimentin. The average expression of Vimentin was strikingly tenfold higher with n-3 HFD and 20 fold higher with n-6 HFD (Figure 4C). In sum, HFD exposure leads to a loss in the expression of luminal genes and a gain in expression of mesenchymal and BC invasion genes in MGs.
BC is a heterogeneous disease. Identifying dietary factors that contribute to the development of different BC subtypes is important for the prevention and treatment of the disease. A HFD has been shown to increase BC risk in humans [8, 9] and several studies using rodent models have shown that exposure to a HFD in utero and/or during postnatal development significantly increase carcinogen stimulated MG tumorigenesis [19–21]. However, the developmental processes in the MG impacted by a HFD that mediate the heterogeneity of BC are unclear. Recently, work by Luitjen et al. revealed that HFD exposure during MG development leads to sustained changes in gene expression. In the present study, we found that genes induced by a HFD in the mouse developing MG were strongly expressed in the human basal-like subtype of BC, suggesting that HFD exposed MGs may have abnormal basal-like characteristics. In line with this, we found increased expression of mesenchymal genes by GSEA including the well described marker Vimentin in HFD gland samples. Our analysis reveals that HFD exposure during development leads to an increase in basal-like characteristics in MGs along with a concomitant loss of luminal gene expression (Figures 3–4), which may promote an increase in basal-like BC development. This pattern of gene expression parallels the increase in basal-like BC/decrease in luminal BC that is found in obese women . It remains to be determined as to which HFD exposed cells express the basal-like genes and whether the MG architecture is dramatically altered in a sustained manner in these models and in humans. Nonetheless, HFD profoundly affects MG gene expression in a sustained way, inducing genes that are expressed in basal-like BC.
Here we show that of the HFD-induced genes that were significantly associated with prognosis in BC, a striking 93% was significantly associated with poor prognosis. All but one these genes were still strongly associated with poor prognosis in the non-basal-like set of BCs (249 tumors examined). Therefore, the association between HFD-induced genes and prognosis in BC is highly compelling and spans many subtypes. While some skepticism could be drawn due to the relevance of cross-species analysis, a strong body of evidence has demonstrated the parallels between rodent and human MG structures, MG development and carcinogenesis [29–33]. Furthermore, mouse MG tumors have been shown to replicate the diversity of human BCs including tumor initiation, hormone dependency and histopathology across multiple MG endpoints [31–33]. Thus, our data suggests a model in which HFD impacts MG development to promote poor prognosis cancer.
Interestingly, we provide evidence that a HFD, regardless of content (n-3 or n-6 PUFAs), may preferentially promote basal-like BC. Epidemiological evidence supports the idea that a HFD increases BC risk. Whether this risk is a consequence of overall percentage of fat in an individual’s diet or type of PUFAs in the diet has long been debated. Numerous rodent studies have investigated the dietary impact of (n-3) and (n-6) PUFAs on BC risk. While rodent studies support the notion that diets high in (n-6) PUFAs increase breast carcinogenesis [34–37], diets primarily consisting of n-3 PUFAs have been suggested to protect against BC [38, 39]. However, Hilakivi-Clarke et al., showed that prepubertal rats exposed to a high (n-3) PUFA diet resulted in key biological changes within the MG reflecting an increased susceptibility to BC. Furthermore, studies in adult rats have shown that exposure to a diet high in (n-3) PUFAs do not inhibit carcinogen induced mammary tumorigenesis, but may promote it [40, 41]. While the timing of exposure could play a factor (fetal, adolescence, adult), human studies have yet to provide conclusive evidence that diets consisting of primarily (n-3) or (n-6) PUFAs play a profound role in increasing or decreasing BC risk. [42–44]. Thus, the influence of a HFD on BC development maybe more a consequence of the amount fat consumed regardless of the source and our data supports this hypothesis.
Given the striking association between HFD-induced genes during development and prognosis in BC, it is imperative to examine these phenomena in humans. The greatest strides in decreasing cancer mortality have been made by prevention as with the link between smoking and lung cancer. It will be important to determine whether developmental HFD exposure as in this study and/or exposure in adults lead to the induction of poor prognosis genes. Also, it will be critical to determine whether these changes are permanent or require the presence of HFD or are maintained for only a certain period of time. It will also be crucial to determine how any of the maintained changes in gene expression are propagated. Based on the HFD mediated induction of DNA (cytosine-5-)-methyltransferase 3 alpha (DNMT3A), perhaps an epigenetic change is responsible for these sustained shifts in gene expression. DNMT3A functions in de novo methylation, is important in development and altered expression levels have been found in several different types of human cancers . If the link between poor prognosis in breast cancer and developmental HFD exposure bare out in humans, efforts can be made to prevent such exposures. In addition, if this link holds true in human breast cancer, it may be possible to identify women who are more susceptible to poor prognosis cancers based on the expression of these newly identified HFD-induced basal-like characteristics such as the mesenchymal marker Vimentin. Importantly, vimentin has been shown to play a significant role in the epithelial-mesenchymal transition (EMT) process in BC [46, 47]. EMT is a cellular reprograming process that changes the shape of epithelial cells to exhibit a more motile mesenchymal phenotype (reviewed in ). In vitro Vimentin over expression in human BC cells has been shown to contribute to this process, thus increasing BC cell motility and invasive properties . In human breast cancer samples, vimentin expression is found in high-grade ductal carcinomas with low ER expression levels . Thus, identifying early changes in vimentin expression in breast epithelium due to HFD exposure might be used to predict clinical outcomes translating to better preventive treatment strategies.
In sum, we show that HFD induces the expression of genes that are associated with poor prognosis in BC. HFD exposed MGs showed high expression of the mesenchymal marker Vimentin as well as a loss in luminal markers. We propose that developmental HFD exposure leads to sustained changes in gene expression that promote the development of poor prognosis cancers including basal-like BC.
Availability of supporting data
The data supporting the results of the article is included within the article and in supplementary materials.
High fat diet
Epidermal growth factor receptor 2
Gene set enrichment analysis
Polyunsaturated fatty acids
False discovery rate
Matrix metalloproteinase 12
DNA (cytosine-5-)-methyltransferase 3 alpha
Netherlands cancer institute
Normalized enrichment scores.
Ziegler RG, Hoover RN, Pike MC, Hildesheim A, Nomura AM, West DW, Wu-Williams AH, Kolonel LN, Horn-Ross PL, Rosenthal JF, Hyer MB: Migration patterns and breast cancer risk in Asian-American women. J Natl Cancer Inst. 1993, 85: 1819-1827.
Thomas D, Karagas M: Migrant Studies. Cancer Epidemiology and Prevention. Edited by: Schottenfeld D, Fraumeni JF. 1996, New York: Oxford University Press, 236-254. 3
McCredie M: Cancer epidemiology in migrant populations. Recent Results Cancer Res. 1998, 154: 298-305.
Deapen D, Liu L, Perkins C, Bernstein L, Ross RK: Rapidly rising breast cancer incidence rates among Asian-American women. Int J Cancer. 2002, 99: 747-750.
Pike MC, Kolonel LN, Henderson BE, Wilkens LR, Hankin JH, Feigelson HS, Wan PC, Stram DO, Nomura AM: Breast cancer in a multiethnic cohort in Hawaii and Los Angeles: risk factor-adjusted incidence in Japanese equals and in Hawaiians exceeds that in whites. Canc Epidemiol Biomarker Prev. 2002, 11: 795-800.
John EM, Phipps AI, Davis A, Koo J: Migration history, acculturation, and breast cancer risk in Hispanic women. Canc Epidemiol Biomarker Prev. 2005, 14: 2905-2913.
Murtaugh MA, Sweeney C, Giuliano AR, Herrick JS, Hines L, Byers T, Baumgartner KB, Slattery ML: Diet patterns and breast cancer risk in Hispanic and non-Hispanic white women: the Four-Corners Breast Cancer Study. Am J Clin Nutr. 2008, 87: 978-984.
Boyd NF, Stone J, Vogt KN, Connelly BS, Martin LJ, Minkin S: Dietary fat and breast cancer risk revisited: a meta-analysis of the published literature. Br J Cancer. 2003, 89: 1672-1685.
Thiebaut AC, Kipnis V, Chang SC, Subar AF, Thompson FE, Rosenberg PS, Hollenbeck AR, Leitzmann M, Schatzkin A: Dietary fat and postmenopausal invasive breast cancer in the National Institutes of Health-AARP Diet and Health Study cohort. J Natl Cancer Inst. 2007, 99: 451-462.
van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347: 1999-2009.
Chang HY, Nuyten DS, Sneddon JB, Hastie T, Tibshirani R, Sorlie T, Dai H, He YD, van't Veer LJ, Bartelink H, van de Rijn M, Brown PO, van de Vijver MJ: Robustness, scalability, and integration of a wound-response gene expression signature in predicting breast cancer survival. Proc Natl Acad Sci USA. 2005, 102: 3738-3743.
Perou CM, Jeffrey SS, van de Rijn M, Rees CA, Eisen MB, Ross DT, Pergamenschikov A, Williams CF, Zhu SX, Lee JC, Lashkari D, Shalon D, Brown PO, Botstein D: Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci USA. 1999, 96: 9212-9217.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffry SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Lonning PE, Borresen-Dale A: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874.
Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, Deng S, Johnsen H, Pesich R, Geisler S, Demeter J, Perou CM, Lonning PE, Brown PO, Borresen-Dale AL, Botstein D: Repeated Observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA. 2003, 14: 8418-8423.
Carey LA: Race, breast cancer subtypes, and survival in the Carolina Breast Cancer Study. JAMA. 2006, 21: 2492-2502.
Rutqvist LE, Johansson H: Long-term follow-up of the randomized stockholm trial on adjuvant tamoxifen among postmenopausal patients with early stage breast cancer. Acta Oncologica. 2007, 2: 133-145.
Turkoz FP, Solak M, Petekkaya I, Keskin O, Kertmen N, Sarici F, Arik Z, Babacan T, Ozisik Y, Altundag K: The prognostic impact of obesity on molecular subtypes of breast cancer in premenopausal women. J BUON. 2013, 18: 335-341.
Hilakivi-Clarke L, Clarke R, Onojafe I, Raygada M, Cho E, Lippman M: A maternal diet high in n-6 polyunsaturated fats alters mammary gland development, puberty onset, and breast cancer risk among female rat offspring. Proc Natl Acad Sci USA. 1997, 94: 9372-9377.
Hilakivi-Clarke L, Clarke R, Lippman M: The influence of maternal diet on breast cancer risk among female offspring. Nutrition. 1999, 15: 392-401.
Moral R, Escrich R, Solanas M, Vela E, Costa I, de Villa MC, Escrich E: Diets high in corn oil or extra-virgin olive oil provided from weaning advance sexual maturation and differentially modify susceptibility to mammary carcinogenesis in female rats. Nutr Cancer. 2011, 63: 410-420.
Hilakivi-Clarke L, Olivo SE, Shajahan A, Khan G, Zhu Y, Zwart A, Cho E, Clarke R: Mechanisms mediating the effects of prepubertal (n-3) polyunsaturated fatty acid diet on breast cancer risk in rats. J Nutr. 2005, 135: 2946-2952.
Olivo-Marston SE, Zhu Y, Lee RY, Cabanes A, Khan G, Zwart A, Wang Y, Clarke R, Hilakivi-Clarke L: Gene signaling pathways mediating the opposite effects of prepubertal low-fat and high-fat n-3 polyunsaturated fatty acid diets on mammary cancer risk. Cancer Prev Res (Phila). 2008, 1: 532-545.
Medvedovic M, Gear R, Freudenberg JM, Schneider J, Bornschein R, Yan M, Mistry MJ, Hendrix H, Karyala S, Halbleib D, Heffelfinger D, Clegg DJ, Anderson MW: Influence of fatty acid diets on gene expression in rat mammary epithelial cells. Physiol Genomics. 2009, 38: 80-88.
Luijten M, Singh AV, Bastian CA, Westerman A, Pisano MM, Pennings JL, Verhoef A, Green ML, Piersma AH, de Vries A, Knudsen TB: Lasting effects on body weight and mammary gland gene expression in female mice upon early life exposure to n-3 but not n-6 high-fat diets. PLoS One. 2003, 8: e55603-
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34: 267-273.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550.
Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: AmiGO: online access to ontology and annotation data. Bioinformatics. 2009, 25: 288-289.
Rudland PS, Barraclough R, Fernig DG, Smith JA: Growth and differentiation of the normal mammary gland and its tumours. Biochem Soc Symp. 1998, 63: 1-20.
Medina D: Chemical carcinogenesis of rat and mouse mammary glands. Breast Dis. 2007, 28: 63-68.
Russo J, Russo IH: Development of the human breast. Maturitas. 2004, 49: 2-15.
Russo J, Russo I: Molecular basis of breast cancer. 2004, New York: Springer-Verlag
Singh M, McGinley JN, Thompson HJ: A comparison of the histopathology of premalignant and malignant mammary gland lesions induced in sexually immature rats with those occurring in the human. Lab Invest. 2000, 80: 221-231.
Carroll KK: Dietary fats and cancer. Am J Clin Nutr. 1991, 53: 1064S-1067S.
Fay MP, Freedman LS, Clifford CK, Midthune DN: Effect of different types and amounts of fat on the development of mammary tumors in rodents: a review. Cancer Res. 1997, 57: 3979-3988.
Hopkins GJ, Carroll KK: Relationship between amount and type of dietary fat in promotion of mammary carcinogenesis induced by 7,12-dimethylbenz[a]anthracene. J Natl Cancer Inst. 1979, 72: 1009-1012.
Rose DP, Connolly JM: Effects of dietary omega-3 fatty acids on human breast cancer growth and metastases in nude mice. J Natl Cancer Inst. 1993, 85: 1743-1747.
Maillard V, Bougnoux P, Ferrari P, Jourdan ML, Pinault M, Lavillonniere F, Body G, Le Floch O, Chajes V: N-3 and N-6 fatty acids in breast adipose tissue and relative risk of breast cancer in a case–control study in Tours, France. Int J Cancer. 2002, 98: 78-83.
Kaizer L, Boyd NF, Kriukov V, Tritchler D: Fish consumption and breast cancer risk: an ecological study. Nutr Cancer. 1989, 12: 61-68.
Sasaki T, Kobayashi Y, Shimizu J, Wada M, In'nami S, Kanke Y, Takita T: Effects of dietary n-3-to-n-6 polyunsaturated fatty acid ratio on mammary carcinogenesis in rats. Nutr Cancer. 1998, 30: 137-143.
Cohen LA, Chen-Backlund JY, Sepkovic DW, Sugie S: Effect of varying proportions of dietary menhaden and corn oil on experimental rat mammary tumor promotion. Lipids. 1993, 28: 449-456.
MacLennan M, Ma DW: Role of dietary fatty acids in mammary gland development and breast cancer. Breast Cancer Res. 2010, 12: 211-
Hunter DJ, Spiegelman D, Adami HO, Beeson L, van den Brandt PA, Folsom AR, Fraser GE, Goldbohm RA, Graham S, Howe GR, Kushi LH, Marshall JR, McDermott A, Miller AB, Speizer FE, Wolk A, Yaun S, Willett W: Cohort studies of fat intake and the risk of breast cancer–a pooled analysis. N Engl J Med. 1996, 334: 356-361.
Smith-Warner SA, Spiegelman D, Adami HO, Beeson WL, van den Brandt PA, Folsom AR, Fraser GE, Freudenheim JL, Goldbohm RA, Graham S, Kushi LH, Miller AB, Rohan TE, Speizer FE, Toniolo P, Willett WC, Wolk A, Zeleniuch-Jacquotte A, Hunter DJ: Types of dietary fat and breast cancer: a pooled analysis of cohort studies. Int J Cancer. 2001, 92: 767-774.
Kim MS, Kim YR, Yoo NJ, Lee SH: Mutational analysis of DNMT3A gene in acute leukemias and common solid cancers. APMIS. 2013, 121: 85-94.
Vuoriluoto K, Haugen H, Kiviluoto S, Mpindi JP, Nevo J, Gjerdrum C, Tiron C, Lorens JB, Ivaska J: Vimentin regulates EMT induction by Slug and ocogenic H-Ras and migration by governing Axl expression in breast cancer. Oncogene. 2011, 12: 1436-1448.
Korsching E, Packeisen J, Liedtke C, Hungermann D, Wulfing P, van Diest PJ, Brandt B, Boecker W, Buerger H: The origin of vimentin expression in invasive breast cancer: epithelial-mesenchymal transition, myoepithelial histogenesis or histogenesis from progenitor cells with bilinear differentiation potential?. J Pathol. 2005, 206: 451-457.
Thiery JP: Epithelial-mesenchymal transitions in tumor progression. Nat Rev Cancer. 2002, 2: 442-454.
Domagala W, Lasota J, Bartkowiak J, Weber K, Osborn M: Vimentin contributes to human mammary epithelial cell migration. J Cell Sci. 1999, 112 (Pt 24): 4615-4625.
RCMCh was supported by an undergraduate training grant from the Howard Hughes Medical Institute Precollege and Undergraduate Science Education Program- award number 52007568. RKD was supported by a Faculty Research Council (FRC) grant from the University of Texas-Pan American (UTPA). The content of this publication is the responsibility of the authors and does not represent the views of the HHMI.
The authors declare that they have no competing interests.
RKD conceived and designed the study. RCMCh and MK participated in the design of the study and conducted the experimental analysis. All authors participated in the analysis of data and review of the manuscript. RCMCh did the primary writing of this manuscript, but all authors contributed in this regard. All authors read and approved the final manuscript.