- Research article
- Open Open Peer Review
microRNA alterations in ALDH positive mammary epithelial cells: a crucial contributing factor towards breast cancer risk reduction in case of early pregnancy
BMC Cancervolume 14, Article number: 644 (2014)
microRNAs have recently succeeded in grabbing the center stage in cancer research for their potential to regulate vital cellular process like cell cycle, stem cell renewal and epithelial mesenchymal transition. Breast cancer is the second most leading cause of cancer related mortality in women. The main reason for mortality is chemoresistance and metastasis for which remnant stem cells are believed to be the cause. One of the natural ways to reduce the risk of breast cancer in women is early pregnancy. Unraveling the mechanism behind it would add to our knowledge and help in evolving newer paradigms for breast cancer prevention.
The current study deals with investigating transcriptomic differences in putative stem cells in mammary epithelial cell population (MECs) in terms of genes and microRNAs. In silico tools were used to identify potential mechanisms. ALDH positive MECs represent a putative stem cell population in the mammary gland.
MECs were extracted from the mammary gland of virgin and parous (one time pregnant) rats. ALDH positive MECs were sorted and used for transcriptional and translational analysis for genes and microRNAs. In silico analysis for target prediction and networking was performed through online portals of Target Scan and Metacore.
A total of 35 and 49 genes and microRNAs respectively were found to be differentially expressed within the two groups. Among the important genes were Lifr, Acvr1c, and Pparγ which were found to be targeted by microRNAs in our dataset like miR-143, miR-30, miR-140, miR-27b, miR-125a, miR-128ab, miR-342, miR-26ab, miR-181, miR-150, miR-23ab and miR-425. In silico data mining and networking also demonstrates that genes and microRNA interaction can have profound effects on stem cell renewal, cell cycle dynamics and EMT processes of the MEC population.
Our data clearly shows that certain microRNAs play crucial role in the regulation of ALDH positive MECs and favor an anti-carcinogenic environment in the post-partum gland. Some of the potential interplaying mechanisms in the ALDH positive MEC population identified through this study are p21, Lifr and Pparγ mediated cell cycle regulation, regulation of metastasis and expansion of stem cell pool respectively.
Breast cancer is one of the primary causes of cancer-related deaths and the most common malignancy in women worldwide . Despite different existing and potential strategies to treat breast cancer, its global incidence is predicted to increase in the coming years at a rate of 3.2 million new cases every year by 2050 . Currently, while there are good treatment options for patients with breast cancer, a high percentage of patients develop resistance to these treatments over time, and these therapies often have undesirable and harmful side effects.
It is well known that early, full-term pregnancy reduces the risk of breast cancer. A completed full-term pregnancy before the age of 20 years reduces breast cancer risk by 50% compared to nulliparous women . While this fact has been known for a long time, different reasons have been brought forward repeatedly to explain this phenomenon. Unfortunately, the mechanism behind this protective phenomenon is not well defined. One possible explanation is that parous women often have a different hormonal profile compared to nulliparous women . Thus, it is thought that alterations in the hormonal milieu both during and after pregnancy may contribute to the phenomenon of parity-induced protection against breast cancer. Animal studies have proven that short-term treatment with pregnancy levels of estrogen can be effective in reducing mammary cancer incidence [5, 6]. Also of special interest is the hormone prolactin, which has been found at reduced levels in the sera of parous women; coincidentally, prolactin-suppressing drugs have been shown to reduce mammary tumors . Further, growth hormone has also been demonstrated to be vital for breast cancer development and parity reduces the levels of growth hormone in circulation [8–10]. Thus, strong evidence suggests a definitive role for hormones in parity-induced protection against breast cancer. Further, some researchers have suggested that pregnancy results in the terminal differentiation of the mammary gland, resulting in the loss of a particular cell population that is prone to malignant transformation [11, 12]. However, other studies indicate that differentiation of the mammary gland per se is not sufficient to explain the phenomenon of parity-induced protection against breast cancer [5, 6, 13].
It has been well established that the mammary gland is partly comprised of a population of epithelial stem cells that are capable of self-renewal and are responsible for the generation of newer cell types specific to the gland. Therefore, a third theory was proposed that breast cancer arises primarily from the stem cell compartment and pregnancy may lead to protective changes in the stem cell population of the mammary gland. However, it remains highly debatable whether the mammary epithelial stem cell population is a primary contributing factor to the phenomenon of parity-induced protection [14–17], and additional work in this area is therefore needed. A recent report by Siwko et al. suggested that there is a persistent decrease in the number of mammary-repopulating units (mammary epithelial stem cells) after parity. In contrast, there are reports, including ours, which demonstrate that parity-induced protection is not due to changes in the number of cells in the mammary epithelium itself but is the result of systemic changes in the whole organism [18–20]. Thus, it is imperative to understand how the systemic environment influences mammary epithelial stem cells and how this may contribute to the protective phenomenon of parity. The significance of this study lies in the notion that stem cells are the initiators of carcinogenesis, according to the cancer stem cell theory.
Over the last couple of years, research in the field of breast cancer, has added strong lines of evidence, supporting the fact that microRNAs have a significant role to play in the regulation of the signaling pathways involved in oncogenesis. They have also been implicated in the maintenance of cancer stem cells via their ability to affect multiple pathways including cell proliferation, cell death [21–23], cell- cell communication and cell adhesion .
In this study, we demonstrate that pregnancy alters molecular processes in ALDH positive MECs (putative mammary epithelial stem cells), leading to a decreased risk of mammary cancer. Here, we primarily focused on the genetic differences of ALDH positive MECs from both virgin and parous animals, through gene and microRNA profiling. To the best of our knowledge, this is the first study to identify the parity-induced microRNA signature in the ALDH positive MECs that is associated with the reduced risk of breast cancer.
Virgin Lewis rats were purchased from Harlan Sprague–Dawley (Indianapolis and San Diego). The rats were housed in temperature controlled room with 12-h light/dark schedule. They were fed (Teklad 8640; Teklad, Madison, WI) and water ad libitum. To generate parous animals, seven week old virgin rats were mated with similar aged male rats. The pups were removed from the cage right after parturition. The mammary glands from the parous rats were removed six weeks after parturition to allow involution of the gland. Mammary glands were removed from age-matched virgin rats and used as controls. All procedures performed were approved and conducted in accordance with the Texas Tech University Health Sciences Center Institutional Animal Care and Use Committee’s guidelines.
Mammary epithelial cell isolation and stem cell enrichment
Isolation of MECs from the mammary gland was carried out using collagenase assisted cell dissociation. Briefly, all six pairs of mammary glands from both virigin and uniparous (one-time pregnant) Lewis rats were processed by mechanical and enzymatic dissociation to prepare a single-cell suspension. MECs were cultured overnight. They were then stained with ALDEFLUOR (Stem Cell Technologies), which enabled the selection of putative stem cells from MECs with strong aldehyde dehydrogenase (ALDH) activity. An inhibitor of ALDH was used as the negative control. The cells were analyzed on a flow cytometer in the green fluorescence channel (520–540 nm) and were then subjected to sorting. Sorting was performed using a flow cytometric cell sorter (BD FACS Aria) to collect the ALDHbright population. These sorted ALDH positive MECs were then used to prepare RNA and protein lysates for transcriptional and translational analysis, respectively.
The mammosphere formation assay was performed as previously described . Briefly, MECs from both virgin and parous animals were plated after preparation of a single-cell suspension using a 23-G needle. Cells were plated in ultralow attachment 6-well plates (Corning) with mammosphere media containing B27 supplement and LONZA Single Quot supplements (hydrocortisone, insulin, beta-mercaptoethanol, EGF, and gentamycin) in phenol red-free DMEM/F12 media (GIBCO). The cell density of this assay was optimized to 500 cells/cm2. The cells were not disturbed for 5 days before any change in media. After 7 days, any sphere larger than 50 μm was considered for counting and further analysis, using a sample size of six.
The mammospheres were collected by centrifugation at 115 × g for 5 min and were gently suspended in 200 μl mammosphere media. They were then plated on poly-lysine coated; 8-well chambered slides with mammosphere assay containing 1% fetal bovine serum and incubated at 37°C, with 5% CO2 for 3–4 hrs for attachment. These mammospheres were then fixed using 5% formaldehyde and blocked with 5% bovine serum albumin for 1 hr. They were then stained for stem cell markers using primary antibodies against SOX2 (goat IgG clone Y-17, 1:100 dilutions, Santa Cruz Biotechnology) and OCT3/4 (mouse IgG2b clone C-10, 1:50 dilution, Santa Cruz Biotechnology). Alexafluor 488 and 594 were used as secondary antibodies raised in species appropriate for the primary antibody. The spheres were washed and counterstained with DAPI and mounted. All slides were examined using a Nikon confocal microscope (Eclipse Ti, Nikon, Japan). Multicolor images were collected sequentially in three channels.
An EdU (5-ethnyl-2′-deoxyuridine) based kit; Click-iT EdU Imaging kit was used to perform the assay (Molecular probes, Life technologies). The sorted ALDH positive MECs were plated in the 8 well chamber slide with 1 × 104 cells/ well and incubated overnight at 37°C/ 5% CO2. 10 μM of EdU was incubated with the cells for 2 hrs at 37°C/5% CO2. The cells were then fixed with 3.7% formaldehyde for 15 min and permeabilized with 0.5% Triton-X-100 for 20 min. It was then incubated with Alexa fluor azide for 30 min to enable the detection of EdU. They were finally counterstained with DAPI and mounted for examination using a Nikon confocal microscope (Eclipse Ti, Nikon, Japan).
Gene and microRNA arrays and in silicoanalysis
Gene and global microRNA profiles were generated using the SABiosciences PCR (Cat No. PARN-405Z and PARN-047Z) and miRNome array (Cat No. MIRN-216Z), respectively. Briefly, RNA was extracted from ALDHpositive MECs of both virgin and parous animals using Trizol (Invitrogen). Replicates of mammary tissues from at least six animals from each group were used for gene and miRNome array analyses. For the gene array, 160 genes associated with stemness and stem cell development were analyzed. For the microRNA array, 653 of the most abundantly expressed and well-characterized microRNAs in the rat microRNA genome as annotated by miRBase Release 16 were profiled. Genes or microRNAs with at least a 2-fold increase or decrease in expression were considered significantly up- or downregulated, respectively.
For microRNA target prediction, the main in silico approach used depends on sequence complementarity. To predict the targets for the microRNA array, we used the online portal of TargetScan. The targets in which the paired sites were highly conserved were considered for further analysis. TargetScan (http://www.targetscan.org/) predicts the biological targets of microRNAs by searching for the presence of conserved 8mer and 7mer sites that match the seed region of each microRNA. To increase the signal-to-noise ratio, TargetScan requires strict complementarity between the seed region of the microRNA and the predicted target. TargetScan Human considers matches to annotate human UTRs and their orthologs, as defined by UCSC whole-genome alignments (http://genome.ucsc.edu/). Conserved targeting has also been detected within open reading frames (ORFs). MetaCore from Thomson Reuters was used to perform data mining and pathway analysis for the differentially regulated microRNAs and genes.
Western blot analysis
Protein lysates of ALDH positive MECs were prepared from glands of virgin and parous animals. Protein concentration was measured using the Pierce BCA (bicinchoninic acid) protein assay (Thermo Scientific) according to the manufacturer’s specifications. Equivalent amounts of protein (1 – 5 μg) were resolved using reducing SDS-PAGE in 4-20% gradient pre-cast Mini-Protean TGX polyacrylamide gels (Bio-Rad). Resolved proteins were transferred onto PVDF immunoblotting membranes and probed with the following antibodies as per manufacturer’s instructions: LIFR, PPAR γ, CYCLIN D, ACVR1C, LIN 28 (Santa Cruz), NOTCH 2, CYCLIN E1, SNAIL, SLUG, VIMENTIN, N-CADHERIN, E-CADHERIN, ZEB (Cell Signaling) and P21 (Abcam). Each membrane was also stripped and reprobed for β-actin as protein loading. For chemiluminescent detection of proteins, SuperSignal West Femto Chemiluminescent Substrate (ThermoScientific) was used according to the manufacturer’s instructions. Immunoblots were processed digitally on the ImageQuant LAS4000 biomolecular imager (GE Healthcare). The signal intensities for each antibody was densitometrically analyzed and normalized to actin bands.
Data are expressed as means with standard deviation or standard error. Student’s t-test was used to determine statistical significance between 2 groups. A value of p < 0.05 was considered a statistically significant difference.
The main aim of this study was to determine the molecular differences in ALDH positive MECs as a result of pregnancy. The model system that we used for this study is FACS-enriched MECs positive for ALDH from virgin and parous animals.
Parity does not influence the proportion or stemness of ALDH positive MECs
First, we quantified ALDH positive MECs from virgin and parous animals; we did not observe a statistically significant change in the number of ALDH positive MECs between the two groups (p = 0.22) (Figure 1A & 1B). Additionally, there were no statistically significant differences in their mammosphere-forming capacities (p = 0.35) (Figure 1C). The presence of stem cells in mammospheres derived from virgin and parous MECs was confirmed using stemness markers SOX2 and OCT4 (Figure 1D). Overall, the data indicate that parity does not influence the stem cell population by significantly altering the percentage of ALDH positive MECs compared to a virgin mammary gland. On contrary, our cell proliferation experiment indicated that virgin ALDH positive MECs had a slightly higher proliferation rate than parous ALDH positive MECs (Additional file 1: Figure S1).
Parity induces alterations in the genetic environment of the ALDH positive MEC population
We next investigated potential differences in the expression of genes that influence stemness and stem cell properties of ALDH positive MECs in virgin and parous animals. We identified a total of 21% (35/168) of genes that were differentially expressed in parous compared to virgin ALDH positive MECs. We observed that leukemia inhibitory factor receptor (Lifr), RGM domain family member A (Rgma), fibroblast growth factor 3 (Fgf3), and activin A receptor, type IC (Acvr1c) were all strongly upregulated in parous ALDH positive MECs (Figure 2 ; Table 1). Further, Lifr was the most upregulated gene, exhibiting a 7.6-fold increase in parous compared to virgin ALDH positive MECs. We investigated the protein levels of some of these differentially regulated genes, including LIFR, NOTCH2, ACVR1C, p21, CYCLIN D, CYCLIN E1, LIN 28, PPAR γ, SNAIL, SLUG, VIMENTIN, N-CADHERIN, E-CADHERIN and ZEB2. Interestingly, we found the same trend was reflected at the protein level as well (Figure 3). The downregulated genes in parous ALDH positive MECs included the metastasis-promoting gene, Zeb2, cell cycle-associated E2F transcription factor 5 (E2f5) cyclin E1 (Ccne1) and retinoblastoma-like protein 1 (Rbl1) indicating inhibitory activity at the G1–S phase of the cell cycle (Figure 2 ; Table 2).
Several of the differentially regulated genes were found to be involved in cell cycle regulation, stem cell self renewal, and differentiation, including Ccne1, Fgf1, Fgf3, Notch2, Myst histone acetyltransferase 1 (Myst1), Neurogenin 2 (Neurog2), Forkhead box A2 (Foxa2), and ISL Lim homeobox 1 (Isl1). In addition, our stem cell gene array data indicates that TGF-β, NOTCH, and WNT pathways which are involved in stem cell regulation were most influenced by parity. Interestingly, genes associated with identification of stemness, such as Abcg2 and Aldh1a1, were down regulated in parous ALDH positive MECs.
MicroRNA profiles in virgin and parous ALDH positive MECs
With the aim of finding novel mechanisms for parity-induced protection against breast cancer, we next investigated global microRNA expression in parous and virgin ALDH positive MECs. The global microRNA profiling data revealed differential regulation of 7.5% (49/653) of the total microRNAs between virigin and parous ALDH positive MECs. Among the differentially regulated microRNAs 12.2% (6/49) and 87.8% (43/49), were up and down regulated respectively; in parous compared to virgin ALDH positive MECs (Figure 4, Table 3 & 4). This screening confirmed the differential expression of several microRNAs in parous ALDH positive MECs that have been reported to be associated with breast and other cancers (Table 5). Many of the downregulated microRNAs in parous ALDH + MECs have been previously reported to be pro-carcinogenic (Table 5). In contrast, a few of the upregulated microRNAs in this case have been previously reported as anti-carcinogenic in nature (Table 5).
miR-497, miR-218a, miR-378*, miR-503, miR-7a/7c and miR-221 were upregulated in parous ALDH positive MECs. Most of these microRNAs have been shown to have tumor suppressive properties in breast and other cancers. miR-497, miR-218a, miR-378 and miR-503 were reported to have anti-carcinogenic effect in breast, glioma, liver, endometrial and gastric cancers respectively. These microRNAs suppress tumor growth by affecting cell proliferation, stem cell renewal, angiogenesis and cancer cell metabolism [26–31]. miR-221 has been shown to be involved in the promotion of epithelial-mesenchymal transition (EMT) in breast cancer cell lines, where it is regulated by Slug . However, in normal human umbilical vein endothelial cells, miR-221 has been shown to directly target and repress Zeb2 .
Both miR-27b and miR-181b were among the most downregulated microRNAs in parous compared to nulliparous ALDH positive MECs. Positive expressions of both these microRNAs have been correlated with poor clinical outcomes as they target genes associated with tumor suppression and neoplastic transformation (high mobility group A proteins) [34, 35].
Thus, the microRNA prolife of parous ALDH positive MECs appears to favor anti-carcinogenesis, at least considering the differentially expressed microRNAs observed in this study. However, functional studies to further support this statement are warranted.
In silicoanalysis reveals interconnections between differentially regulated microRNAs and genes
The results of the in silico analysis utilizing TargetScan revealed that many of the differentially regulated genes in our gene array can possibly be regulated by some of the differentially regulated microRNAs from the global microRNA profile (Table 3 & 4). Among the most targeted genes was Lifr, which drew our attention. It was found to be targeted by multiple microRNAs in our analysis, including miR-143, miR-30 family members, miR-140, miR-27b, miR-125a-5p, miR-128ab, and miR-342-3p. Interestingly, all these microRNAs were downregulated in the miRnome analysis, and correspondingly, Lifr was upregulated at the transcription and translational levels in parous ALDH positive MECs (Figure 3). We also looked at genes other than the ones differentially expressed in our gene array analysis Figure 5 and found that miR-125a targets many cell cycle genes, like Ccnd1, p21, and Cdk2, which may prove to be instrumental in unraveling the mechanism of parity-induced protection against breast cancer. Previously, miR-125a-5p was also reported to target Lin28, which is again very interesting in the context of stem cells. Here, we have demonstrated that parous ALDH positive MECs demonstrate a higher expression of Lin28 at the protein level than nulliparous ALDH positive MECs (Figure 3B).
Another highly targeted gene in our in silico analysis was Acvr1c or Alk7, a receptor involved in the nodal-activin pathway. It was found to be a putative target for let-7 family members, miR-26ab, miR-181 family, miR-150, miR-27b, miR-23ab, miR-425, miR-125a-5p, and miR-128ab. However, miR-125a-5p is the only microRNA that was found to target both Acvr1c and Lifr. Thus, we believe that Lifr and Acvr1c could be the major players in the modulation of ALDH positive MECs in response to pregnancy. Furthermore, because these genes are regulated by miR-125a-5p, it is possible that miR-125a-5p plays a key regulatory role in putative stem cells of MECs and is involved in parity-induced protection against breast cancer.
It is well known that early parity leads to a reduced risk of breast cancer. In this regard, numerous theories trying to explain the underlying mechanism have been postulated. However, at present, it has not been possible to completely unravel the mechanism behind this phenomenon. Previously, investigators have studied the differences between MECs of virgin and parous females of both humans and rodents. However, in this study, we have directed our efforts to look for the transcriptomic changes in the ALDH positive MEC population in virgin and parous rats. Our aim was to determine the underlying mechanism of parity-induced reduction in breast cancer risk in terms of changes introduced during the process of pregnancy in the ALDH positive MEC population. Therefore, in contrast to previous studies where changes in the entire mammary epithelium were analyzed, our study specifically scrutinized changes in putative stem cells of the MEC population.
To the best of our knowledge, this is the first report to describe changes in the global microRNA profiles between virgin and parous ALDH positive MECs. Our data identified the gene and microRNA signature induced by pregnancy in the ALDH positive MECs. Our in vitro data and in silico analysis led us to speculate about some of the possible phenomenon that can be attributed to the parity-induced reduction in breast cancer risk. We show here that ALDHpositive MECs of parous might differ from virgin in terms of regulation of the cell cycle, EMT processes, and tumor suppressors. The mammary gland is considered as an actively cycling tissue, which means that putative stem cells are under constant pressure to make decisions on cell cycle progression, proliferation, cell cycle arrest, and differentiation. There are several genes (Ccne1, Cdk2, E2f5 and Rbl1) in our dataset that demonstrate the possibility of an alteration in the regulation of cell cycle of ALDH positive MECs in response to pregnancy. Down regulation of microRNAs, like miR-28, miR-125a, and miR-503, could possibly lead to up regulation of p21 in parous ALDH positive MECs. This increase in p21 levels may lead to cell cycle arrest in the putative stem cells. This is further supported by our data, which show a down regulation of other significant cell cycle-associated genes, including Ccne1, Cdk2, E2f5, and Rbl1, in parous ALDH positive MECs. p21 is known to be suppressed in many human cancers; it was recently shown in a preclinical study that the up regulation of p21 leads to an anti-proliferative effect in breast cancer cell lines with cell cycle arrest in G1 phase . Further, our cell proliferation data also indicates a higher proliferation rate in virgin ALDH positive MECs. These finding suggests that faster cycling ALDH positive MECs in virgin mammary glands would have a much greater chance of accumulating genetic mutations after carcinogenic insult compared to ALDH positive MECs in parous mammary gland, which may cycle at a much lower rate.
However, this speculation does not completely explain the phenomenon of parity-induced reduction in breast cancer risk. It has been observed in animal models that upon carcinogenic insult, the frequency of occurrence of neoplastic lesions is almost similar in both virgin and parous mammary glands . This observation suggests that lesions in virgin glands are able to progress to breast cancer much more rapidly over time than in parous glands. Therefore, there must be additional factors other than cell cycle regulators that govern the progression of these neoplastic lesions. In this context, it is interesting to note that Lifr is highly upregulated in parous ALDH positive MECs. A recent report identified Lifr as a metastasis suppressor that executes its effects through the Hippo signaling pathway and prevents EMT . In addition, we demonstrated the down regulation of Zeb2 and the up regulation of miR-218a in parous ALDH positive MECs. miR-218a in this study has been predicted to target Zeb2, which is another factor known to suppress EMT. Interestingly, miR-218 has been proved to regulate HMGB1 and suppress migration and invasion. This further supports the idea that ALDH positive MECs from parous mammary glands may not undergo EMT as efficiently as those from virgin glands. If this is indeed the case, then the regulatory environment of ALDH positive MECs in the parous gland does not favor progression of the disease. Thus, we speculate that in response to pregnancy, ALDHpositive MECs are altered in ways that not only favor slower cell cycle rates, but also lower efficiency of EMT like events, thus inhibiting progression of the disease. In this regard, our current data demonstrates that microRNAs likely play a key role in the regulation of these processes and therefore represent an important factor to be considered for their contribution to parity-induced reduction in breast cancer risk.
Another important point worth mentioning is Pparγ (upregulated in parous ALDH positive MECs) is a highly targeted gene by most microRNAs in our dataset. Pparγ has been shown to promote the transcription of tumor-suppressor genes like Pten  and Brca1 . With regards to stem cells, inhibition of Pparγ leads to the expansion of the stem cell pool as represented by the increased percentage of CK5+ and CD29+/CD24+ cells . Interestingly, PPARγ-deficient mice also show an increased occurrence of hormone-dependent mammary cancers . Taken together, we speculate that increased expression of Pparγ in parous ALDH positive MECs contributes to the regulation of stem cell fate determination, thus controlling the incidence of tumorigenesis. According to our in silico analysis, Pparγ is likely regulated by microRNAs like let-7 family members, miR-30 family members, miR-27b, miR-23ab, miR-93, miR-25, miR-128ab, miR-320, and miR-135. In summary, the present study elucidates the differences in the ALDH positive MEC population in both parous and virgin glands. It was generally observed that there are changes at both the gene and microRNA levels in ALDH positive MECs of both groups, and these changes could contribute to reduced risk of breast cancer in parous mammary glands.
Overall our data indicates that pregnancy alters microRNAs and genes in the ALDH positive MECs. These changes in the putative epithelial stem cell population render it resistant to mammary carcinogenesis. Also, our study brings forward a gene and microRNA signature specific to the ALDH positive MECs. These can serve as potential biomarkers and lead to the development of novel therapeutic targets for breast cancer.
Mammary epithelial cells
Epithelial Mesenchymal Transition.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA Cancer J Clin. 2011, 61 (2): 69-90.
Hortobagyi G, de la Garza SJ, Pritchard K, Amadori D, Haidinger R, Hudis CA, Khaled H, Liu MC, Martin M, Namer M, O’Shaughnessy JA, Shen ZZ, Albain KS, ABREAST Investigators: The global breast cancer burden: variations in epidemiology and survival. Clin Breast Cancer. 2005, 6 (5): 391-
MacMahon B, Cole P, Lin TM, Lowe CR, Mirra AP, Ravnihar B, Salber EJ, Valaoras VG, Yuasa S: Age at first birth and breast cancer risk. Bull World Health Organ. 1970, 43 (2): 209-
M N: The Physiology of Reproduction. 1994, New York: Raven
Guzman RC, Yang J, Rajkumar L, Thordarson G, Chen X, Nandi S: Hormonal prevention of breast cancer: mimicking the protective effect of pregnancy. Proc Natl Acad Sci U S A. 1999, 96 (5): 2520-2525.
Rajkumar L, Guzman R, Yang J, Thordarson G, Talamantes F, Nandi S: Short-term exposure to pregnancy levels of estrogen prevents mammary carcinogenesis. Proc Natl Acad Sci U S A. 2001, 98 (20): 11755-11759.
Welsch C, Nagasawa H: Prolactin and mammary tumorigenesis: a review. Cancer Res. 1977, 37 (4): 951-963.
Thordarson G, Jin E, Guzman RC, Swanson SM, Nandi S, Talamantes F: Refractoriness to mammary tumorigenesis in parous rats: is it caused by persistent changes in the hormonal environment or permanent biochemical alterations in the mammary epithelia?. Carcinogenesis. 1995, 16 (11): 2847-2853.
Thordarson G, Semaan S, Low C, Ochoa D, Leong H, Rajkumar L, Guzman RC, Nandi S, Talamantes F: Mammary tumorigenesis in growth hormone deficient spontaneous dwarf rats; effects of hormonal treatments. Breast Cancer Res Treat. 2004, 87 (3): 277-290.
Dearth RK, Delgado DA, Hiney JK, Pathiraja T, Oesterreich S, Medina D, Dees WL, Lee AV: Parity-induced decrease in systemic growth hormone alters mammary gland signaling: a potential role in pregnancy protection from breast cancer. Cancer Prev Res (Phila). 2010, 3 (3): 312-321.
Russo J, Russo I: Biological and molecular bases of mammary carcinogenesis. Lab Invest. 1987, 57 (2): 112-137.
Russo J, Russo I: Toward a unified concept of mammary carcinogenesis. Prog Clin Biol Res. 1997, 396: 1-16.
Sivaraman L, Stephens L, Markaverich B, Clark J, Krnacik S, Conneely O, O’Malley B, Medina D: Hormone-induced refractoriness to mammary carcinogenesis in Wistar-Furth rats. Carcinogenesis. 1998, 19 (9): 1573-1581.
Siwko S, Dong J, Lewis M, Liu H, Hilsenbeck S, Li Y: Evidence that an early pregnancy causes a persistent decrease in the number of functional mammary epithelial stem cells–implications for pregnancy-induced protection against breast cancer. Stem Cells. 2008, 26 (12): 3205-3209.
Britt K, Ashworth A, Smalley M: Pregnancy and the risk of breast cancer. Endocr Relat Cancer. 2007, 14 (4): 907-933.
Russo J, Moral R, Balogh G, Mailo D, Russo I: The protective role of pregnancy in breast cancer. Breast Cancer Res. 2005, 7 (3): 131-142.
Trichopoulos D, Lagiou P, Adami H: Towards an integrated model for breast cancer etiology: the crucial role of the number of mammary tissue-specific stem cells. Breast Cancer Res. 2005, 7 (1): 13-17.
Smith GH, Medina D: A morphologically distinct candidate for an epithelial stem cell in mouse mammary gland. J Cell Sci. 1988, 90 (Pt 1)(Pt 1): 173-183.
Abrams T, Guzman R, Swanson S, Thordarson G, Talamantes F, Nandi S: Changes in the parous rat mammary gland environment are involved in parity-associated protection against mammary carcinogenesis. Anticancer Res. 1998, 18 (6A): 4115-4121.
Rajkumar L, Kittrell F, Guzman R, Brown P, Nandi S, Medina D: Hormone-induced protection of mammary tumorigenesis in genetically engineered mouse models. Breast Cancer Res. 2007, 9 (1): R12-
Yu F, Deng H, Yao H, Liu Q, Su F, Song E: Mir-30 reduction maintains self-renewal and inhibits apoptosis in breast tumor-initiating cells. Oncogene. 2010, 29 (29): 4194-4204.
Howe EN, Cochrane DR, Richer JK: Targets of miR-200c mediate suppression of cell motility and anoikis resistance. Breast Cancer Res. 2011, 13 (2): R45-
Radisky DC: miR-200c at the nexus of epithelial-mesenchymal transition, resistance to apoptosis, and the breast cancer stem cell phenotype. Breast Cancer Res. 2011, 13 (3): 110-
Liu C, Kelnar K, Liu B, Chen X, Calhoun-Davis T, Li H, Patrawala L, Yan H, Jeter C, Honorio S, Wiggins JF, Bader AG, Fagin R, Brown D, Tang DG: The microRNA miR-34a inhibits prostate cancer stem cells and metastasis by directly repressing CD44. Nat Med. 2011, 17 (2): 211-215.
Shaw F, Harrison H, Spence K, Ablett M, Simões B, Farnie G, Clarke R: A detailed mammosphere assay protocol for the quantification of breast stem cell activity. J Mammary Gland Biol Neoplasia. 2012, 17 (2): 111-117.
Li D, Zhao Y, Liu C, Chen X, Qi Y, Jiang Y, Zou C, Zhang X, Liu S, Wang X, Zhao D, Sun Q, Zeng Z, Dress A, Lin M, Kung H, Rui H, Liu L, Mao F, Jiang B, Lai L: Analysis of MiR-195 and MiR-497 expression, regulation and role in breast cancer. Clin Cancer Res. 2011, 17 (7): 1722-1730.
Tu Y, Gao X, Li G, Fu H, Cui D, Liu H, Jin W, Zhang Y: MicroRNA-218 Inhibits Glioma Invasion, Migration, Proliferation, and Cancer Stem-like Cell Self-Renewal by Targeting the Polycomb Group Gene Bmi1. Cancer Res. 2013, 73 (19): 6046-6055.
Xiao F, Zhang W, Chen L, Chen F, Xie H, Xing C, Yu X, Ding S, Chen K, Guo H, Cheng J, Zheng S, Zhou L: MicroRNA-503 inhibits the G1/S transition by downregulating cyclin D3 and E2F3 in hepatocellular carcinoma. J Transl Med. 2013, 11 (1): 195-5876-11-195
Xu YY, Wu HJ, Ma HD, Xu LP, Huo Y, Yin LR: MicroRNA-503 suppresses proliferation and cell-cycle progression of endometrioid endometrial cancer by negatively regulating cyclin D1. FEBS J. 2013, 280 (16): 3768-3779.
Zhou B, Ma R, Si W, Li S, Xu Y, Tu X, Wang Q: MicroRNA-503 targets FGF2 and VEGFA and inhibits tumor angiogenesis and growth. Cancer Lett. 2013, 333 (2): 159-169.
Eichner LJ, Perry MC, Dufour CR, Bertos N, Park M, St-Pierre J, Giguere V: miR-378(*) mediates metabolic shift in breast cancer cells via the PGC-1beta/ERRgamma transcriptional pathway. Cell Metab. 2010, 12 (4): 352-361.
Lambertini E, Lolli A, Vezzali F, Penolazzi L, Gambari R, Piva R: Correlation between Slug transcription factor and miR-221 in MDA-MB-231 breast cancer cells. BMC Cancer. 2012, 12: 445-455.
Chen Y, Banda M, Speyer C, Smith J, Rabson A, Gorski DH: Regulation of the expression and activity of the antiangiogenic homeobox gene GAX/MEOX2 by ZEB2 and microRNA-221. Mol Cell Biol. 2010, 30 (15): 3902-3913.
Mansueto G, Forzati F, Ferraro A, Pallante P, Bianco M, Esposito F, Iaccarino A, Troncone G, Fusco A: Identification of a New Pathway for Tumor Progression: MicroRNA-181b Up-Regulation and CBX7 Down-Regulation by HMGA1 Protein. Genes Cancer. 2010, 1 (3): 210-224.
Jin L, Wessely O, Marcusson EG, Ivan C, Calin GA, Alahari SK: Prooncogenic factors miR-23b and miR-27b are regulated by Her2/Neu, EGF, and TNF-alpha in breast cancer. Cancer Res. 2013, 73 (9): 2884-2896.
Wang Y, Rathinam R, Walch A, Alahari SK: ST14 (suppression of tumorigenicity 14) gene is a target for miR-27b, and the inhibitory effect of ST14 on cell growth is independent of miR-27b regulation. J Biol Chem. 2009, 284 (34): 23094-23106.
Bhushan L, Kandpal RP: EphB6 receptor modulates micro RNA profile of breast carcinoma cells. PLoS One. 2011, 6 (7): e22484-
Liu S, Patel SH, Ginestier C, Ibarra I, Martin-Trevino R, Bai S, McDermott SP, Shang L, Ke J, Ou SJ, Heath A, Zhang KJ, Korkaya H, Clouthier SG, Charafe-Jauffret E, Birnbaum D, Hannon GJ, Wicha MS: MicroRNA93 regulates proliferation and differentiation of normal and malignant breast stem cells. PLoS Genet. 2012, 8 (6): e1002751-
Chen HY, Lin YM, Chung HC, Lang YD, Lin CJ, Huang J, Wang WC, Lin FM, Chen Z, Huang HD, Shyy JY, Liang JT, Chen RH: miR-103/107 promote metastasis of colorectal cancer by targeting the metastasis suppressors DAPK and KLF4. Cancer Res. 2012, 72 (14): 3631-3641.
Farazi TA, Horlings HM, Ten Hoeve JJ, Mihailovic A, Halfwerk H, Morozov P, Brown M, Hafner M, Reyal F, van Kouwenhove M, Kreike B, Sie D, Hovestadt V, Wessels LF, van de Vijver MJ, Tuschl T: MicroRNA sequence and expression analysis in breast tumors by deep sequencing. Cancer Res. 2011, 71 (13): 4443-4453.
Heneghan HM, Miller N, Kelly R, Newell J, Kerin MJ: Systemic miRNA-195 differentiates breast cancer from other malignancies and is a potential biomarker for detecting noninvasive and early stage disease. Oncologist. 2010, 15 (7): 673-682.
Jedlinski DJ, Gabrovska PN, Weinstein SR, Smith RA, Griffiths LR: Single nucleotide polymorphism in hsa-mir-196a-2 and breast cancer risk: a case control study. Twin Res Hum Genet. 2011, 14 (5): 417-421.
Savad S, Mehdipour P, Miryounesi M, Shirkoohi R, Fereidooni F, Mansouri F, Modarressi MH: Expression analysis of MiR-21, MiR-205, and MiR-342 in breast cancer in Iran. Asian Pac J Cancer Prev. 2012, 13 (3): 873-877.
Chen Y, Zhang J, Wang H, Zhao J, Xu C, Du Y, Luo X, Zheng F, Liu R, Zhang H, Ma D: miRNA-135a promotes breast cancer cell migration and invasion by targeting HOXA10. BMC Cancer. 2012, 12: 111-2407-12-11.
Shen L, Li J, Xu L, Ma J, Li H, Xiao X, Zhao J, Fang L: miR-497 induces apoptosis of breast cancer cells by targeting Bcl-w. Exp Ther Med. 2012, 3 (3): 475-480.
Jiang D, Wang X, Liu X, Li F: Gene delivery of cyclin-dependent kinase inhibitors p21 (Waf1) and p27 (Kip1) suppresses proliferation of MCF-7 breast cancer cells in vitro. Breast Cancer. 2013, Epub ahead of print
Yang J, Yoshizawa K, Nandi S, Tsubura A: Protective effects of pregnancy and lactation against N-methyl-N-nitrosourea-induced mammary carcinomas in female Lewis rats. Carcinogenesis. 1999, 20 (4): 623-628.
Chen D, Sun Y, Wei Y, Zhang P, Rezaeian A, Teruya-Feldstein J, Gupta S, Liang H, Lin H, Hung M, Ma L: LIFR is a breast cancer metastasis suppressor upstream of the Hippo-YAP pathway and a prognostic marker. Nat Med. 2012, 18 (10): 1511-1517.
Patel L, Pass I, Coxon P, Downes C, Smith S, Macphee C: Tumor suppressor and anti-inflammatory actions of PPARgamma agonists are mediated via upregulation of PTEN. Curr Biol. 2001, 11 (10): 764-768.
Pignatelli M, Cocca C, Santos A, Perez-Castillo A: Enhancement of BRCA1 gene expression by the peroxisome proliferator-activated receptor gamma in the MCF-7 breast cancer cell line. Oncogene. 2003, 22 (35): 5446-5450.
Yin Y, Yuan H, Zeng X, Kopelovich L, Glazer R: Inhibition of peroxisome proliferator-activated receptor gamma increases estrogen receptor-dependent tumor specification. Cancer Res. 2009, 69 (2): 687-694.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/644/prepub
We thank Texas Tech University Health Sciences Center Paul L. Foster School of Medicine, El Paso, TX for supporting this research.
The authors declare that they have no competing interests.
SN designed and carried out majority of the experiments, RS performed the immunoblotting analysis, VR performed the in silico analysis, RL helped in the performance of cell sorting and editing the manuscript, AA and TB helped in the revising the manuscript critically and LR conceived the idea, designed and supervised the project and helped to draft the manuscript. All authors read and approved the final manuscript.