Promotion of variant human mammary epithelial cell outgrowth by ionizing radiation: an agent-based model supported by in vitro studies
- Rituparna Mukhopadhyay†1,
- Sylvain V Costes†1,
- Alexey V Bazarov2,
- William C Hines1,
- Mary Helen Barcellos-Hoff1, 3 and
- Paul Yaswen1Email author
© Mukhopadhyay et al.; licensee BioMed Central Ltd. 2010
Received: 26 August 2009
Accepted: 10 February 2010
Published: 10 February 2010
Most human mammary epithelial cells (HMEC) cultured from histologically normal breast tissues enter a senescent state termed stasis after 5 to 20 population doublings. These senescent cells display increased size, contain senescence associated β-galactosidase activity, and express cyclin-dependent kinase inhibitor, p16INK4A (CDKN2A; p16). However, HMEC grown in a serum-free medium, spontaneously yield, at low frequency, variant (v) HMEC that are capable of long-term growth and are susceptible to genomic instability. We investigated whether ionizing radiation, which increases breast cancer risk in women, affects the rate of vHMEC outgrowth.
Pre-stasis HMEC cultures were exposed to 5 to 200 cGy of sparsely (X- or γ-rays) or densely (1 GeV/amu 56Fe) ionizing radiation. Proliferation (bromodeoxyuridine incorporation), senescence (senescence-associated β-galactosidase activity), and p16 expression were assayed in subcultured irradiated or unirradiated populations four to six weeks following radiation exposure, when patches of vHMEC became apparent. Long-term growth potential and p16 promoter methylation in subsequent passages were also monitored. Agent-based modeling, incorporating a simple set of rules and underlying assumptions, was used to simulate vHMEC outgrowth and evaluate mechanistic hypotheses.
Cultures derived from irradiated cells contained significantly more vHMEC, lacking senescence associated β-galactosidase or p16 expression, than cultures derived from unirradiated cells. As expected, post-stasis vHMEC cultures derived from both unirradiated and irradiated cells exhibited more extensive methylation of the p16 gene than pre-stasis HMEC cultures. However, the extent of methylation of individual CpG sites in vHMEC samples did not correlate with passage number or treatment. Exposure to sparsely or densely ionizing radiation elicited similar increases in the numbers of vHMEC compared to unirradiated controls. Agent-based modeling indicated that radiation-induced premature senescence of normal HMEC most likely accelerated vHMEC outgrowth through alleviation of spatial constraints. Subsequent experiments using defined co-cultures of vHMEC and senescent cells supported this mechanism.
Our studies indicate that ionizing radiation can promote the outgrowth of epigenetically altered cells with pre-malignant potential.
Carcinogenic consequences of radiation exposure have historically been attributed to targeted effects - misrepaired DNA damage directly caused by dose-dependent ionization events in the cell of cancer origin. Radiation can also induce non-targeted effects - altered cytokines and signaling that affect the cellular composition and microenvironment of irradiated tissues , and non-mutational, but heritable changes that alter cell-cell interactions and induce persistent phenotypes associated with malignant progression [2–4]. The potential carcinogenic contribution of these non-targeted effects, which are typically not linearly proportional to radiation dose, has not been well studied, particularly in primary human epithelial cells.
In this study, we used primary cultures of human mammary epithelial cells (HMEC) as an experimental system to directly evaluate the potential of ionizing radiation to promote the outgrowth of cells bearing a pre-malignancy-associated epigenetic change. In serum-free growth medium, HMEC from histologically normal breast tissues arrest growth after 5 to 20 population doublings, exhibit senescent morphologies, and express p16INK4A (CDKN2A; p16) [5, 6]. This p16-dependent form of senescence, termed stasis, is distinguished by irreversible growth arrest with 2N DNA content (reviewed in ). Stasis is associated with poorly defined imbalances in signal transduction brought on by cell culture conditions or oncogene activation, but is not directly associated with DNA damage or dysfunctional telomeres [8, 9]. Stasis requires activation of another well-known tumor suppressor, pRB, which functions downstream of p16, and serves as a block to indefinite proliferation (immortality) - a prerequisite for malignant transformation. The heterogeneous p16 expression observed in human breast epithelial cells in situ , and frequent aberrations in the p16-pRB pathway in human tumors , suggest that conditions that influence its expression and silencing have physiological and pathological relevance.
HMEC cultured in a serum-free medium spontaneously yield rare variant (vHMEC) cells in which p16 genes are methylated and silenced at frequencies that differ among normal specimens from women . In previous studies, such vHMEC have been shown to be susceptible to genomic instability associated with telomere and centrosome dysfunction [6, 13, 14]. In some cases following carcinogen or oncogene exposure, these vHMEC give rise to immortalized clones bearing chromosomal aberrations commonly observed in primary human breast tumors . Here, we sought to determine whether low to moderate doses of radiation (5-200 cGy) would alter the rate of vHMEC outgrowth from non-malignant tissues. In replicate experiments using HMEC from four different specimens, we found that radiation caused shorter growth plateaus and significantly increased the rate at which p16(-) vHMEC grew out in long term cultures. Computer simulations using an agent-based model suggested that radiation accelerated selection of pre-existing vHMEC, a prediction that was confirmed experimentally. Thus non-targeted radiation effects can lower a critical cancer barrier by altering cell interactions to promote outgrowth of cells with pre-malignant phenotypes.
Materials and methods
Histologically normal breast tissues, obtained from surgically discarded reduction mammoplasty (B1400 and N17) and prophylactic mastectomy (B1389 & B1450) specimens with patients' informed consent and institutional review board approval, were provided by the UCSF Cancer Center and the Cooperative Human Tissue Network. Specimens B1400 and N17 were derived from disease-free women aged 50 and 17 years, respectively. Specimen B1389 was derived from a 53-year old woman with a family history of breast and ovarian cancer, diagnosed with non-proliferative fibrocystic disease. Specimen B1450 was derived from a woman of unspecified age, diagnosed with invasive ductal and lobular carcinoma in the contralateral breast. Tissue samples were minced and enzymatically dissociated using 0.1% w/v collagenase I in Dulbecco's Modified Eagle Medium at 37°C for 12 to 18 hrs. Small tissue fragments (organoids) remaining after digestion, were collected by centrifugation at 100 × g for two minutes. These organoids were stored frozen or seeded directly into Mammary Epithelial Cell Growth Medium (MEGM; Lonza, Walkersville, MD, USA). The resulting HMEC were cultured in serum-free MEGM medium as previously described  and verified to be mycoplasma-free (Bionique Testing Laboratories (Saranac Lake, NY, USA). HMEC were routinely subcultured when 80% confluent and reseeded at a densities of 5 × 103 cells/cm2. Total population doublings were estimated using the equation, Population Doublings (PD) = log (A/B)/log2, where A was the number of collected cells and B was the number of plated cells. The effects of plating efficiencies were not taken into account in these estimates, and the calculated values were therefore lower than actual values. In some experiments, 5-aza-2'-deoxycytidine was added at a concentration of 3.3 μM to the culture medium for 72 to 96 hours prior to harvest.
Passage 4 HMEC were routinely seeded at a density of 0.25 × 106 cells per T-25 flask. The cells were grown for seven to eight days to 60% confluence prior to radiation exposure. Control plates were sham-irradiated. After irradiation, the cells were allowed to recover for 48 hrs, then dissociated with trypsin and replated at densities of 0.25 × 106 cells per 100 mm dish. The cultures were then incubated for four to six weeks. When visual inspection indicated that the largest patches of vHMEC had attained diameters of 1.3 to 2.5 cm, all the cultures were harvested for analysis. To generate growth curves, both irradiated and unirradiated HMEC cultures were maintained in triplicate through Passage 9.
Low linear energy transfer (LET) exposures were conducted using either 160 kVp X-ray or 137Cs γ-irradiation. Cultures were exposed to single doses of 5, 25, 50, 100 or 200 cGy at dose rates of either 22.5 cGy/minute (5, 25, or 50 cGy) or 1 Gy/minute (100 or 200 cGy). High-LET radiation exposures were conducted using a 1 GeV/amu 56Fe ion source at the NASA Space Radiation Laboratory at Brookhaven National Laboratory (Upton, NY, USA).
Senescence associated β-galactosidase (SA-βgal) activity
For visualization, cells were washed, fixed, and incubated overnight at 37°C with X-gal chromogenic substrate at pH 6.0 as described . The cells were viewed and photographed using a phase contrast Nikon Eclipse TS100 microscope (Nikon Instruments Inc., Melville, NY, USA). For quantitation of SA-βgal (+) and (-) cells by flow cytometry, fluorescein digalactoside (FDG; Molecular Probes, Eugene, OR, USA) was used as a substrate as described [18–20]. Cells were dissociated and resuspended in phosphate buffer saline containing 5% fetal bovine serum (PBS/FBS). Diluted FDG solution was mixed with an equal volume of cell suspension, incubated at room temperature for three minutes and quenched by adding 10 volumes of PBS/FBS. Propidium iodide (1.25 μg/ml, Molecular Probes, Eugene, OR, USA) was added to exclude dead cells from analysis. Flow cytometry was performed using a FACScan flow cytometer and Cell Quest Pro software (Becton Dickinson, Franklin Lakes, NJ, USA). Determinations of the relative sensitivities of HMEC and vHMEC to radiation-induced senescence were estimated using the formula: , in which Ps(0) and Ps(D) were the respective percentages of SA-βgal(+) cells in a population before and after exposure to a radiation dose (D). The formula included the factor 1-Ps(0) to correct for SA-βgal(+) cells pre-existing in the population.
p16 immunohistochemistry and immunofluorescence
Cells were fixed with 4% paraformaldehyde for 30 minutes, washed and permeabilized with 0.1% Triton: 3% H2O2 for five minutes, incubated with a 1:300 dilution of primary antibodies (JC2, Neomarkers, Fremont, CA, USA), and washed three to four times. For visualization, the cells were incubated with VectaStain ABC reagent, followed by DAB substrate (Vector Laboratories, Burlingame, CA, USA). For quantitation, the cells were incubated with Alexa Fluor 488 conjugated goat anti-mouse IgG (1:200) secondary antibodies. Nuclei were counterstained with 0.5 ng/ml 4,6-diamidino-2-phenylindole (DAPI; Molecular Probes, Eugene, OR, USA). Images for individual channels were acquired and quantitated using a Cellomics Array Scan VTI (Thermo Scientific, Pittsburgh, PA, USA), and merged using Adobe Photoshop 6.0 software (Adobe Systems, San Jose, CA, USA).
Cells were incubated for the final 24 hr prior to harvest with 10 μM bromodeoxyuridine (BrdU; BD Biosciences, San Jose, CA, USA), then dissociated, washed with medium containing 2% FCS, and fixed with 95% ethanol. Nuclei were prepared by incubating the fixed cells with 0.8% pepsin for 20 minutes at 37°C, followed by an IFA buffer (10 mM HEPES, 150 mM NaCl, 4% FBS, and 0.1% sodium azide) wash containing 0.5% Tween. The nuclei were then incubated with a 1:5 dilution of fluoroscein-conjugated anti-BrdU antibodies in IFA (BD Biosciences, San Jose, CA, USA) on ice for 30 minutes, followed by treatment with RNase A (5 μg/ml) and propidium iodide (5 μg/ml). Flow cytometry was performed as described above.
Genomic DNA was isolated from cells using the Wizard Genomic DNA Isolation kit (Promega, Madison, WI, USA). A total of 500 to 1,000 ng of DNA was treated with bisulfite using the EZ DNA Methylation-Gold kit (Zymo Research, Orange, CA, USA). Polymerase chain reaction amplification of the p16 gene promoter region was performed using the primer set 5' TTT TTA GAG GAT TTG AGG GAT AGG 3' (-159 to -136) and 5' CTA CCT AAT TCC AAT TCC CCT ACA 3' (+209 to +233) [21, 22] and the following conditions: 95°C/two minutes × one cycle; 96°C/20 sec, 60°C/20 sec and 72°C/90 sec × 40 cycles; 72°C/five minutes × one cycle. PCR products were purified with the QIAquick gel extraction kit (Qiagen, Valencia, CA, USA) and cloned in the pGEMT vector (Promega, Madison, WI, USA). Individual clones were sequenced with M13 forward and/or reverse primers.
Total cell lysates were prepared in SDS/Tris buffer with protease inhibitors (Complete, Roche Applied Science, Indianapolis, IN, USA) and 20 to 30 μg/lane were separated on gradient SDS-PAGE gels (Invitrogen). After transfer, the separated proteins were incubated with a 1:1000 dilution of p16 antibodies (JC2, Neomarkers), followed by imaging using an Odyssey infrared imaging system (Licor Biosciences, Lincoln, Nebraska, USA).
Graphpad Prism 5 (GraphPad Software, Inc., La Jolla, CA, USA) and JMP IN 3 (JMP, Cary, NC, USA) software were used for all the statistical analyses.
Simulations were performed using Matlab software (The MathWorks, Natick, MA, USA). Agents were dispersed randomly in silico to simulate plating onto a 400 × 400 grid (1.5 mm/pixel) at initial densities of 1,200 agents/grid (20% confluence). The following rules were used for modeling the growth of HMEC: 1) agents could divide as long as there was open space surrounding them; 2) if there was no more open space, agents could still divide but would be compressed until they reached a minimum size; 3) when an agent had reached the maximum number of divisions allowed, it would stop dividing and its type changed permanently to a senescent type; 4) senescent agents could not reattach after detachment during subculture. The program was reinitiated using the resulting agents when the grid reached 80% saturation (approximately 6,000 agents/grid) to simulate re-plating. For each set of experimental conditions, five independent simulations were performed, leading to an average behavior with a relative standard error less than 10%. By representing the grid as a tensor whose vectors are properties of the included agents, the progression of the tensor could be visualized using the advanced imaging platform DIPimage (Image Processing Toolbox for Matlab, Delft University of Technology, Delft, The Netherlands). A growth plateau was defined as the period during which the rate of population doubling deviated significantly from an exponential rate. Plateau width was determined experimentally as the time to 80% confluence of HMEC cultures in which the majority of cells displayed morphological features of senescence immediately after plating.
The progeny of HMEC exposed to a moderate radiation dose form larger, more numerous patches of SA-βgal (-) and p16 (-) cells
Soluble factors induced by ionizing radiation, such as TGFβ, affect the growth of unirradiated cells . To determine whether the effect of radiation observed on the outgrowth of surviving cells was due to secreted factors, conditioned medium from irradiated or unirradiated cultures was added to the growth medium of unirradiated cultures. One set of irradiated and three sets of unirradiated B1400 or B1450 cultures were maintained separately. One unirradiated set served as a control and received only fresh medium at each feeding. A second unirradiated set received a 1:1 mixture of conditioned medium from irradiated cultures and fresh medium, while a third unirradiated set received a 1:1 mixture of conditioned medium from unirradiated cultures and fresh medium. Four replicates for each culture condition were maintained in this manner for four to six weeks. We did not observe any significant differences in the outgrowth of vHMEC in cultures supplemented with conditioned medium from irradiated or unirradiated cultures (data not shown). Thus soluble factors were not a likely cause of the differences observed.
We examined unirradiated and irradiated cultures 48 hrs following irradiation at 4p. The cultures exhibited similar cell densities and morphologies; areas of small proliferative cells were interspersed with large, flat non-mitotic cells in both (Figure 1B, upper panels). However, four to six weeks after irradiation, large uniform patches of small proliferative cells were evident earlier and were more numerous in the subcultured (5p) irradiated populations than in the unirradiated populations (Figure 1B, lower panels, and Figure 1C).
We next measured the expression of p16 protein by immunofluorescence in unirradiated and irradiated cultures harvested at identical times (Figure 2C). More than 20 random fields from unirradiated and irradiated B1450 and B1400 cultures were analyzed. Significantly more (70.5 ± 3.3% vs. 44.8 ± 3.5, P < 0.0001) cells were p16 (-) in the irradiated cultures than those in the unirradiated cultures (Figure 2D).
The p16 gene is silenced epigenetically in the progeny of both unirradiated and irradiated HMEC
To determine whether radiation resulted in qualitative and/or quantitative differences in the methylation of CpG sites in the p16 promoter, we used bisulfite sequencing to investigate the methylation status of 35 CpG sites spanning from -159 to +233, relative to the translation start site (+1), in the core CpG island region of the p16 gene in two specimens, B1450 and B1400, with and without irradiation (Figure 3C). The majority of the CpG sites were unmethylated in pre-stasis 3p HMEC cultures harvested prior to stasis, whereas p16 promoter methylation was significantly increased in both unirradiated and irradiated vHMEC cultures, as indicated by the % methylation of individual CpG sites (Figure 3D). The methylation patterns were qualitatively and quantitatively distinct in irradiated and unirradiated cultures, however there was no discernable quantitative or qualitative correlation among p16 promoter methylation and passage level or radiation exposure status.
Radiation of different doses and qualities promotes outgrowth of vHMEC
Densely ionizing 56Fe ions have been reported to be more potent inducers of genomic instability, cell transformation, and tumorigenesis than sparsely ionizing X-rays. The relative biological effectiveness (RBE) for cell killing of 1 GeV/amu 56Fe particles compared to X-rays was 1.80 at D10 (dose resulting in 10% survival by cell killing data) for HMEC . RBE values reported for other endpoints such as tumorigenesis in mice have been much higher [27, 28]. We investigated the effect of radiation quality on the outgrowth of vHMEC. The progeny of cells generated from three breast tissue specimens (B1450, B1400, N17) were irradiated with densely ionizing 1 GeV/amu 56Fe particles and compared to those irradiated with sparsely ionizing 137Cs γ-rays four to six weeks after exposure, using flow cytometry to measure SA-βgal activity and BrdU incorporation (Figure 4B and 4C) in 5 to 10 replicates from each specimen for each treatment. Although the response of individual specimens to radiation exposure varied, as expected, the cultures irradiated with 137Cs γ-rays contained significantly more SA-βgal (-) and BrdU (+) cells, as in the X-ray experiments, than unirradiated controls (P < 0.0001; Wilcoxon test). Notably, the responses to 1 GeV/amu 56Fe ions were comparable to those of γ-rays and were not dose dependent.
Agent-based modeling can be used to accurately simulate HMEC growth kinetics
We considered two possible explanations for our observations that radiation promoted vHMEC outgrowth; first that radiation induced additional vHMEC and second, that radiation accelerated the selection of pre-existing vHMEC, consistent with studies by Holst et al.  demonstrating their existence in human breast tissues. To evaluate the likelihood of these alternative mechanisms, we defined underlying assumptions to use an agent-based model (ABM) to simulate vHMEC outgrowth.
ABMs are computer simulations that represent systems as collections of autonomous decision-making entities called agents . Each agent is programmed to respond to situations using a set of contextual rules. These models are non-deterministic, are typically based on multiple simulations so that a range of behaviors can be established for a population, and have proven to be very useful in predicting emerging properties from complex systems [31–36]. Recently, ABMs have been used to predict the responses of cancer stem cell populations to ionizing radiation .
Based on experimental observations, we knew that the growth of HMEC is contact inhibited, that the area occupied by individual cells is reduced with increasing culture density, that all HMEC have limited replication capacity, that after a limited number of divisions most HMEC are capable of expressing p16 that initiates senescence, and that senescent cells do not readily replate following trypsinization. Thus we defined the following agent rules: 1) agents replicate as long as there is space surrounding them; 2) if there is no open space, agents can still replicate as smaller entities until reaching a minimum size; 3) when agents reach the maximum number of replications allowed, they stop dividing and their type changes permanently to a senescent type; 4) senescent agents do not reattach after detachment during subculture.
Agent-based modeling predicts the accelerated outgrowth of the progeny of irradiated HMEC
A major challenge is to understand how cellular responses to radiation are integrated in a multicellular context to affect human health. In this study, we show how radiation can promote the outgrowth of pre-malignant cells by accelerating senescence of normal cells. We used ABM to show that the rapidity with which radiation-resistant vHMEC selectively populate cultures can be predicted by the initial proportion of the vHMEC present in primary cultures, the differential rate at which sensitive HMEC succumb to radiation-induced loss of growth potential, and the physical space available for expansion.
Together with p53, the p16 protein functions as a sentinel that integrates various cellular signals and stresses to limit proliferation. Methylation of the p16 gene has been identified in situ in histologically normal mammary epithelial cells of disease-free women . A Luria-Delbrück fluctuation analysis suggests that such cells may be the source of vHMEC that continue to proliferate after most HMEC in long term cultures have undergone p16-associated stasis . Expression of p16 appears to be strongly selected against in many human solid tumors, including those of the breast, where the gene encoding p16 is deleted in approximately 20%  and inactivated epigenetically in an additional 20%  of cases. Interestingly, exposure of workers in nuclear weapons manufacturing facilities to plutonium, for example, has been strongly linked (P = 0.03) to methylation of the p16 gene in lung adenocarcinomas . Renal cell carcinomas from patients living in areas contaminated by the Chernobyl accident have also shown aberrant hypermethylation of the p16/p14 locus .
As we and others have observed, radiation itself does not appear to be a direct inducer of p16 (Additional file 1). Our new results indicate that radiation-induced stress can be integrated with p16-inducing factors that cause premature growth arrest and senescence. ABM simulations showed that radiation can advance the selection process, allowing vHMEC to overtake and fill the voids created by the prematurely senescing normal cells. A similar model, in which lesion growth is driven by opportunistic expansion of apoptosis-resistant p53 mutant cells, has recently been proposed for UVB-induced squamous cell carcinomas . The expansion of such a variant population in situ would be expected to expand the target size in which additional malignancy promoting aberrations could occur, especially since this population has been shown to be particularly susceptible to genomic instability caused by telomere dysfunction and centrosome irregularities [6, 13, 14].
While accurately modeling the accelerated outgrowth of vHMEC from irradiated cultures, the current ABM does not explain why, in some cases, the total number of population doublings achieved by vHMEC from irradiated cultures exceeded that of vHMEC from unirradiated controls. This outcome could be due to either a radiation-induced increase in the number of vHMEC or modification of the phenotype of the existing vHMEC. Since the phenotypic feature of the vHMEC that allows them to avoid or overcome stasis is the stable repression of p16 expression, one possibility is that radiation directly affects epigenetic process(es) involved in the initiation or maintenance of this repression. DNA methylation-associated silencing of mammalian genes has been proposed to be a process, during which instances of spontaneous gene reactivation are initially frequent, becoming progressively less frequent over the course of multiple cell divisions [43, 44]. In agreement with this hypothesis, recent experiments indicate that p16 gene methylation occurs progressively in clonal vHMEC populations after silencing has occurred . Irradiation may lead to acceleration of this process and thus reduce the frequency of spontaneous p16 reactivation. Indeed, while occasional p16(+) cells could be observed in 8p post-stasis cultures, they were more frequent in the progeny of unirradiated versus irradiated cultures (data not shown). The mechanism responsible for this difference remains to be discovered, but may involve radiation effects on sentinel proteins such as p53. Notably, radiation has been reported to relieve p53-mediated repression of DNA methyltransferase 1 (Dnmt1) in human colon carcinoma cells . Under certain circumstances, Dnmt1, which has been shown to be continually required for p16 repression , will catalyze de novo methylation of specific promoter CpG islands . Thus, radiation, acting through transient activation of p53, could cause de novo methylation and more stable silencing of the p16 gene.
We have found that radiation can indirectly promote the outgrowth of a putative pre-malignant breast cell population by accelerating senescence of normal breast cells. Differential sensitivity of these cell populations to radiation-induced senescence has not previously been demonstrated. While the mechanism responsible for this differential sensitivity is beyond the scope of our present investigation, the work advances a systems-based paradigm of carcinogenic activity by a prototypic genotoxic agent, ionizing radiation. The ABM simulation of well-described HMEC phenotypes in primary cultures helped to distinguish between induced vHMEC and altered population kinetics as plausible explanations for the observed response to radiation. ABM showed that differential sensitivity to radiation-induced senescence can at least partly explain the enhanced outgrowth of vHMEC from irradiated cultures. This simple in vitro system illustrates the concept that heterogeneous responses of individual cells within a population must be integrated to achieve a system-level understanding of radiation effects.
fetal bovine serum
human mammary epithelial cells
linear energy transfer
- p :
phosphate buffer saline
relative biological effectiveness
senescence associated β-galactosidase
This work was funded by a NASA Specialized Center of Research, a Breast Cancer and the Environment Research Center Grant U01 ES012801 from the National Institute of Environment Health Science (NIEHS) and the National Cancer Institute (NCI), and the Office of Health and Environmental Research, Health Effects Division, United States Department of Energy (contract no. DE-AC02-05CH11231). A.V. Bazarov was supported by Flight Attendant Medical Research Institute Young Clinical Investigator Award No. 032122 and Komen for the Cure Research Grant No. BCTR0707231 and W.C. Hines received support from Komen for the Cure postdoctoral fellowship PDF0707408.
- Barcellos-Hoff MH, Costes SV: A systems biology approach to multicellular and multi-generational radiation responses. Mutat Res. 2006, 597: 32-38.View ArticlePubMedGoogle Scholar
- Andarawewa KL, Erickson AC, Chou WS, Costes SV, Gascard P, Mott JD, Bissell MJ, Barcellos-Hoff MH: Ionizing radiation predisposes nonmalignant human mammary epithelial cells to undergo transforming growth factor beta induced epithelial to mesenchymal transition. Cancer Res. 2007, 67: 8662-8670. 10.1158/0008-5472.CAN-07-1294.View ArticlePubMedGoogle Scholar
- Park CC, Henshall-Powell RL, Erickson AC, Talhouk R, Parvin B, Bissell MJ, Barcellos-Hoff MH: Ionizing radiation induces heritable disruption of epithelial cell interactions. Proc Natl Acad Sci USA. 2003, 100: 10728-10733. 10.1073/pnas.1832185100.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsai KK, Chuang EY, Little JB, Yuan ZM: Cellular mechanisms for low-dose ionizing radiation-induced perturbation of the breast tissue microenvironment. Cancer Res. 2005, 65: 6734-6744. 10.1158/0008-5472.CAN-05-0703.View ArticlePubMedGoogle Scholar
- Hammond SL, Ham RG, Stampfer MR: Serum-free growth of human mammary epithelial cells: rapid clonal growth in defined medium and extended serial passage with pituitary extract. Proc Natl Acad Sci USA. 1984, 81: 5435-5439. 10.1073/pnas.81.17.5435.View ArticlePubMedPubMed CentralGoogle Scholar
- Romanov SR, Kozakiewicz BK, Holst CR, Stampfer MR, Haupt LM, Tlsty TD: Normal human mammary epithelial cells spontaneously escape senescence and acquire genomic changes. Nature. 2001, 409: 633-637. 10.1038/35054579.View ArticlePubMedGoogle Scholar
- Yaswen P, Stampfer MR: Molecular changes accompanying senescence and immortalization of cultured human mammary epithelial cells. Int J Biochem Cell Biol. 2002, 34: 1382-1394. 10.1016/S1357-2725(02)00047-X.View ArticlePubMedGoogle Scholar
- Ramirez RD, Morales CP, Herbert BS, Rohde JM, Passons C, Shay JW, Wright WE: Putative telomere-independent mechanisms of replicative aging reflect inadequate growth conditions. Genes Dev. 2001, 15: 398-403. 10.1101/gad.859201.View ArticlePubMedPubMed CentralGoogle Scholar
- Drayton S, Rowe J, Jones R, Vatcheva R, Cuthbert-Heavens D, Marshall J, Fried M, Peters G: Tumor suppressor p16INK4a determines sensitivity of human cells to transformation by cooperating cellular oncogenes. Cancer Cell. 2003, 4: 301-310. 10.1016/S1535-6108(03)00242-3.View ArticlePubMedGoogle Scholar
- Krishnamurthy J, Torrice C, Ramsey MR, Kovalev GI, Al-Regaiey K, Su L, Sharpless NE: Ink4a/Arf expression is a biomarker of aging. J Clin Invest. 2004, 114: 1299-1307.View ArticlePubMedPubMed CentralGoogle Scholar
- Bartek J, Bartkova J, Lukas J: The retinoblastoma protein pathway in cell cycle control and cancer. Exp Cell Res. 1997, 237: 1-6. 10.1006/excr.1997.3776.View ArticlePubMedGoogle Scholar
- Tlsty TD, Romanov SR, Kozakiewicz BK, Holst CR, Haupt LM, Crawford YG: Loss of chromosomal integrity in human mammary epithelial cells subsequent to escape from senescence. J Mammary Gland Biol Neoplasia. 2001, 6: 235-243. 10.1023/A:1011369026168.View ArticlePubMedGoogle Scholar
- Chin K, de Solorzano CO, Knowles D, Jones A, Chou W, Rodriguez EG, Kuo WL, Ljung BM, Chew K, Myambo K, Miranda M, Krig S, Garbe J, Stampfer M, Yaswen P, Gray JW, Lockett SJ: In situ analyses of genome instability in breast cancer. Nat Genet. 2004, 36: 984-988. 10.1038/ng1409.View ArticlePubMedGoogle Scholar
- McDermott KM, Zhang J, Holst CR, Kozakiewicz BK, Singla V, Tlsty TD: p16(INK4a) prevents centrosome dysfunction and genomic instability in primary cells. PLoS Biol. 2006, 4: e51-10.1371/journal.pbio.0040051.View ArticlePubMedPubMed CentralGoogle Scholar
- Nonet GH, Stampfer MR, Chin K, Gray JW, Collins CC, Yaswen P: The ZNF217 gene amplified in breast cancers promotes immortalization of human mammary epithelial cells. Cancer Res. 2001, 61: 1250-1254.PubMedGoogle Scholar
- Stampfer MR, Yaswen P: Culture models of human mammary epithelial cell transformation. J Mammary Gland Biol Neoplasia. 2000, 5: 365-378. 10.1023/A:1009525827514.View ArticlePubMedGoogle Scholar
- Dimri GP, Lee X, Basile G, Roskelley C, Medrano EE, Rubelji I, Pereira-Smith OM, Peacocke M, Campisi J: A novel biomarker identifies senescent human cells in culture and in aging skin in vivo. Proc Natl Acad Sci USA. 1995, 92: 9363-9367. 10.1073/pnas.92.20.9363.View ArticlePubMedPubMed CentralGoogle Scholar
- Fiering SN, Roederer M, Nolan GP, Micklem DR, Parks DR, Herzenberg LA: Improved FACS-Gal: flow cytometric analysis and sorting of viable eukaryotic cells expressing reporter gene constructs. Cytometry. 1991, 12: 291-301. 10.1002/cyto.990120402.View ArticlePubMedGoogle Scholar
- Nolan GP, Fiering S, Nicolas JF, Herzenberg LA: Fluorescence-activated cell analysis and sorting of viable mammalian cells based on beta-D-galactosidase activity after transduction of Escherichia coli lacZ. Proc Natl Acad Sci USA. 1988, 85: 2603-2607. 10.1073/pnas.85.8.2603.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang NC, Hu ML: A fluorimetric method using fluorescein di-beta-D-galactopyranoside for quantifying the senescence-associated beta-galactosidase activity in human foreskin fibroblast Hs68 cells. Anal Biochem. 2004, 325: 337-343. 10.1016/j.ab.2003.11.012.View ArticlePubMedGoogle Scholar
- Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB: Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci USA. 1996, 93: 9821-9826. 10.1073/pnas.93.18.9821.View ArticlePubMedPubMed CentralGoogle Scholar
- Wong DJ, Foster SA, Galloway DA, Reid BJ: Progressive region-specific de novo methylation of the p16 CpG island in primary human mammary epithelial cell strains during escape from M(0) growth arrest. Mol Cell Biol. 1999, 19: 5642-5651.View ArticlePubMedPubMed CentralGoogle Scholar
- Huschtscha LI, Noble JR, Neumann AA, Moy EL, Barry P, Melki JR, Clark SJ, Reddel RR: Loss of p16INK4 expression by methylation is associated with lifespan extension of human mammary epithelial cells. Cancer Res. 1998, 58: 3508-3512.PubMedGoogle Scholar
- Foster SA, Wong DJ, Barrett MT, Galloway DA: Inactivation of p16 in human mammary epithelial cells by cpG island methylation. Mol Cell Biol. 1998, 18: 1793-1801.View ArticlePubMedPubMed CentralGoogle Scholar
- Brenner AJ, Stampfer MR, Aldaz CM: Increased p16 expression with first senescence arrest in human mammary epithelial cells and extended growth capacity with p16 inactivation. Oncogene. 1998, 17: 199-205. 10.1038/sj.onc.1201919.View ArticlePubMedGoogle Scholar
- Sudo H, Garbe J, Stampfer MR, Barcellos-Hoff MH, Kronenberg A: Karyotypic instability and centrosome aberrations in the progeny of finite life-span human mammary epithelial cells exposed to sparsely or densely ionizing radiation. Radiat Res. 2008, 170: 23-32. 10.1667/RR1317.1.View ArticlePubMedGoogle Scholar
- Alpen EL, Powers-Risius P, Curtis SB, DeGuzman R: Tumorigenic potential of high-Z, high-LET charged-particle radiations. Radiat Res. 1993, 136: 382-391. 10.2307/3578551.View ArticlePubMedGoogle Scholar
- Dicello JF, Christian A, Cucinotta FA, Gridley DS, Kathirithamby R, Mann J, Markham AR, Moyers MF, Novak GR, Piantadosi S, Ricart-Arbona R, Simonson DM, Strandberg JD, Vazquez M, Williams JR, Zhang Y, Zhou H, Huso D: In vivo mammary tumourigenesis in the Sprague-Dawley rat and microdosimetric correlates. Phys Med Biol. 2004, 49: 3817-3830. 10.1088/0031-9155/49/16/024.View ArticlePubMedGoogle Scholar
- Holst CR, Nuovo GJ, Esteller M, Chew K, Baylin SB, Herman JG, Tlsty TD: Methylation of p16(INK4a) promoters occurs in vivo in histologically normal human mammary epithelia. Cancer Res. 2003, 63: 1596-1601.PubMedGoogle Scholar
- Bonabeau E: Agent-based modeling: methods and techniques for simulating human systems. Proc Natl Acad Sci USA. 2002, 99 (Suppl 3): 7280-7287. 10.1073/pnas.082080899.View ArticlePubMedPubMed CentralGoogle Scholar
- Markus M, Bohm D, Schmick M: Simulation of vessel morphogenesis using cellular automata. Math Biosci. 1999, 156: 191-206. 10.1016/S0025-5564(98)10066-4.View ArticlePubMedGoogle Scholar
- Kansal AR, Torquato S, Harsh GI, Chiocca EA, Deisboeck TS: Simulated brain tumor growth dynamics using a three-dimensional cellular automaton. J Theor Biol. 2000, 203: 367-382. 10.1006/jtbi.2000.2000.View ArticlePubMedGoogle Scholar
- Kansal AR, Torquato S, Harsh IG, Chiocca EA, Deisboeck TS: Cellular automaton of idealized brain tumor growth dynamics. Biosystems. 2000, 55: 119-127. 10.1016/S0303-2647(99)00089-1.View ArticlePubMedGoogle Scholar
- Dormann S, Deutsch A: Modeling of self-organized avascular tumor growth with a hybrid cellular automaton. In Silico Biol. 2002, 2: 393-406.PubMedGoogle Scholar
- Chen S, Ganguli S, Hunt CA: An agent-based computational approach for representing aspects of in vitro multi-cellular tumor spheroid growth. Proceedings of the 26th Annual International Conference of the IEEE EMBS: September 1-5, 2004; San Francisco, CA, USA. 2004, 691-694.Google Scholar
- Sun T, McMinn P, Holcombe M, Smallwood R, MacNeil S: Agent based modelling helps in understanding the rules by which fibroblasts support keratinocyte colony formation. PLoS ONE. 2008, 3: e2129-10.1371/journal.pone.0002129.View ArticlePubMedPubMed CentralGoogle Scholar
- Enderling H, Park D, Hlatky L, Hahnfeldt P: The importance of spatial distribution of stemness and proliferation state in determining tumor radioresponse. Math Model Nat Phenom. 2009, 4: 117-133. 10.1051/mmnp/20094305.View ArticleGoogle Scholar
- Sharpless NE: INK4a/ARF: a multifunctional tumor suppressor locus. Mutat Res. 2005, 576: 22-38.View ArticlePubMedGoogle Scholar
- Tlsty TD, Crawford YG, Holst CR, Fordyce CA, Zhang J, McDermott K, Kozakiewicz K, Gauthier ML: Genetic and epigenetic changes in mammary epithelial cells may mimic early events in carcinogenesis. J Mammary Gland Biol Neoplasia. 2004, 9: 263-274. 10.1023/B:JOMG.0000048773.95897.5f.View ArticlePubMedGoogle Scholar
- Belinsky SA, Klinge DM, Liechty KC, March TH, Kang T, Gilliland FD, Sotnic N, Adamova G, Rusinova G, Telnov V: Plutonium targets the p16 gene for inactivation by promoter hypermethylation in human lung adenocarcinoma. Carcinogenesis. 2004, 25: 1063-1067. 10.1093/carcin/bgh096.View ArticlePubMedGoogle Scholar
- Romanenko A, Morell-Quadreny L, Lopez-Guerrero JA, Pellin A, Nepomnyaschy V, Vozianov A, Llombart-Bosch A: The INK4a/ARF locus: role in cell cycle control for renal cell epithelial tumor growth after the Chernobyl accident. Virchows Arch. 2004, 445: 298-304. 10.1007/s00428-004-1056-7.View ArticlePubMedGoogle Scholar
- Chao DL, Eck JT, Brash DE, Maley CC, Luebeck EG: Preneoplastic lesion growth driven by the death of adjacent normal stem cells. Proc Natl Acad Sci USA. 2008, 105: 15034-15039. 10.1073/pnas.0802211105.View ArticlePubMedPubMed CentralGoogle Scholar
- Turker MS: Gene silencing in mammalian cells and the spread of DNA methylation. Oncogene. 2002, 21: 5388-5393. 10.1038/sj.onc.1205599.View ArticlePubMedGoogle Scholar
- Yates PA, Burman R, Simpson J, Ponomoreva ON, Thayer MJ, Turker MS: Silencing of mouse Aprt is a gradual process in differentiated cells. Mol Cell Biol. 2003, 23: 4461-4470. 10.1128/MCB.23.13.4461-4470.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Hinshelwood RA, Melki JR, Huschtscha LI, Paul C, Song JZ, Stirzaker C, Reddel RR, Clark SJ: Aberrant de novo methylation of the p16INK4A CpG island is initiated post gene silencing in association with chromatin remodeling and mimics nucleosome positioning. Hum Mol Genet. 2009, 18: 3098-3109.View ArticlePubMedGoogle Scholar
- Peterson EJ, Bogler O, Taylor SM: p53-mediated repression of DNA methyltransferase 1 expression by specific DNA binding. Cancer Res. 2003, 63: 6579-6582.PubMedGoogle Scholar
- Suzuki M, Sunaga N, Shames DS, Toyooka S, Gazdar AF, Minna JD: RNA interference-mediated knockdown of DNA methyltransferase 1 leads to promoter demethylation and gene re-expression in human lung and breast cancer cells. Cancer Res. 2004, 64: 3137-3143. 10.1158/0008-5472.CAN-03-3046.View ArticlePubMedGoogle Scholar
- Jair KW, Bachman KE, Suzuki H, Ting AH, Rhee I, Yen RW, Baylin SB, Schuebel KE: De novo CpG island methylation in human cancer cells. Cancer Res. 2006, 66: 682-692. 10.1158/0008-5472.CAN-05-1980.View ArticlePubMedGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.