Regulation of early signaling and gene expression in the α-particle and bystander response of IMR-90 human fibroblasts
© Ghandhi et al. 2010
Received: 22 March 2010
Accepted: 29 July 2010
Published: 29 July 2010
Skip to main content
© Ghandhi et al. 2010
Received: 22 March 2010
Accepted: 29 July 2010
Published: 29 July 2010
The existence of a radiation bystander effect, in which non-irradiated cells respond to signals from irradiated cells, is well established. To understand early signaling and gene regulation in bystander cells, we used a bio-informatics approach, measuring global gene expression at 30 minutes and signaling pathways between 30 minutes and 4 hours after exposure to α-particles in IMR-90 fibroblasts.
We used whole human genome microarrays and real time quantitative PCR to measure and validate gene expression. Microarray analysis was done using BRB-Array Tools; pathway and ontology analyses were done using Ingenuity Pathway Analysis and PANTHER, respectively. We studied signaling in irradiated and bystander cells using immunoblotting and semi-quantitative image analysis.
Gene ontology suggested signal transduction and transcriptional regulation responding 30 minutes after treatment affected cell structure, motility and adhesion, and interleukin synthesis. We measured time-dependent expression of genes controlled by the NF-κB pathway; matrix metalloproteinases 1 and 3; chemokine ligands 2, 3 and 5 and interleukins 1β, 6 and 33. There was an increased response of this set of genes 30 minutes after treatment and another wave of induction at 4 hours. We investigated AKT-GSK3β signaling and found both AKT and GSK3β are hyper-phosphorylated 30 minutes after irradiation and this effect is maintained through 4 hours. In bystander cells, a similar response was seen with a delay of 30 minutes. We proposed a network model where the observed decrease in phosphorylation of β-catenin protein after GSK3β dependent inactivation can trigger target gene expression at later times after radiation exposure
These results are the first to show that the radiation induced bystander signal induces a widespread gene expression response at 30 minutes after treatment and these changes are accompanied by modification of signaling proteins in the PI3K-AKT-GSK3β pathway.
Non-targeted effects could significantly enhance risks associated with exposure to low doses of ionizing radiation, which occurs in clinical and environmental contexts. It has been established that signals from irradiated cells travel through medium and cellular junctions to produce changes in gene expression [1, 2], ROS production  and moderate damage to DNA in bystander cells as measured by micronucleus formation . Although there is no direct epidemiological evidence for these risks in humans, the potential importance of bystander effects is highlighted by the recent demonstration of radiation bystander carcinogenesis in a mouse model . In primary fibroblasts the major players transmitting and maintaining signals between cells after irradiation appear to be soluble growth factors, cytokines, reactive oxygen species and extracellular matrix proteins [6, 7]. A wealth of information exists on cellular events that occur 4 hours and later, including studies on gene expression [1, 2, 8–10] cytokine production , γ-H2AX measurement of DNA damage  and chromosomal end-points [13, 14] in directly irradiated cells and bystander cells. However, the events that precede these and other well-studied bystander effects on chromosomes  and DNA damage  are yet to be elucidated. In studies using γ-H2AX as a marker of radiation induced DNA double strand breaks, the response in bystander cells was observed within 20-30 minutes after treatment . A recent study confirmed a burst in cellular ROS levels 30 minutes after irradiation, followed by a 1-2 hr window during which double strand break repair foci were induced in α-particle irradiated and bystander cells . Other studies on signaling in bystanders have proposed that an early increase in production of reactive radicals  and TNFα  after irradiation can induce a cytokine cascade, which is consistent with the large number of signaling and stress response genes induced in this study.
In the present study, we focused on early responses to understand primary events that are more proximal to the bystander signal. 30 minutes after exposure, both irradiated and bystander cells showed a burst of gene expression changes. Gene ontology and pathway analyses of differentially expressed genes at 30 minutes after treatment suggested responses that affect cell structure and motility, signal transduction, transcriptional regulation and cell-to-cell communication. We validated the microarray results by quantitative real-time PCR and found that there was good concordance between these two methods. We were also interested in time-dependent patterns of gene expression and focused our studies on genes that showed induction at 30 minutes in both irradiated and bystander populations. The selected genes encode proteins that are transcriptional targets of NF-κB, and time course analysis of mRNA levels further supported our previous suggestion [2, 3] that this signaling pathway is activated in bystanders in a synchronized manner from 30 minutes onward. From our earlier study that focused on the 4-hour transcriptional response , we predicted the involvement of β-catenin activation in gene expression in irradiated cells. We have now investigated protein modifications in the AKT-GSK3β signaling pathway upstream of β-catenin transcriptional activation. Our results show that the radiation signal can trigger a cascade of changes in AKT-GSK3β-βcatenin pathway almost concomitantly with a widespread gene expression response as early as 30 minutes after exposure.
Early passage (population doubling <35) IMR-90 human lung fibroblasts (Coriell repository, NJ) were sub-cultured in Dulbecco's modified Eagle's medium (Gibco) and Ham's F10 medium in a 1:1 mixture plus 15% fetal bovine serum. Mylar-bottomed culture dishes were prepared as described previously . An inner dish with a base of 38-micron-thick Mylar strips was inserted into a larger dish with a 6-micron Mylar base. The 38-micron Mylar completely shields the alpha particles so that only cells on the thinner Mylar areas of the dish were directly irradiated. Cells seeded in these dishes formed a contiguous layer. Cells were exposed to 0 (sham irradiated) or 50 cGy 4He ions (125 keV per micron) as simulated alpha particles using the track segment mode of the 5.5-MV Singletron accelerator at the Radiological Research Accelerator Facility of Columbia University. This dose corresponds to an average fluence of approximately 6-12 alpha particles per cell in the irradiated sections of the dish and the probability of a cell not receiving a particle is less than 0.25%. Four independent experiments were conducted.
Directly irradiated (outer dish) and bystander (inner dish) cells were separated at specified times after irradiation and RNA was isolated using Ribopure (Ambion, Life Technologies). RNA concentrations were measured using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific) and RNA quality was monitored with the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). All RNA samples had RNA integrity numbers >9.0  and 260 nm/280 nm absorbance ratios >2.
Directly irradiated (outer dish) and bystander (inner dish) cells were separated at specified times (30 minutes, 1 hour and 4 hours) after irradiation and trypsinized. For whole cell lysates, cells were collected, washed and lysed in 25% glycerol, 40 mM HEPES at pH 7.5, 1 mM DTT, 0.35 M NaCl, 0.5% NP-40 and Protease inhibitor mixture (HALT, Thermo Scientific). Separation of nuclear and cytoplasmic fractions was performed as recommended in the manufacturer's protocol for NucBuster™ from EMD Biosciences (Darmstadt, Germany). Protein concentrations were determined using the bicinchoninic acid method (Pierce) and measured using the Nanodrop-1000 spectrophotometer (Thermo Scientific). 50 micrograms of protein was used for western analysis and separated on 10% polyacrylamide gels. Primary antibodies were from Cell Signaling Technology, Boston: anti-AKT (cat# 9272), anti-phospho-AKT (S473) (cat# 9271), anti-GSK3β (cat# 9315), anti-phospho-GSK3β (S9) (cat #9336), anti-β-catenin (cat# 9562) and anti-phospho-β-catenin (S33/37/T41) (cat# 9561). Other antibodies were purchased from Millipore: anti-actin (cat# MAB1501) and Sigma: anti-TBP (cat# T1827). Secondary antibodies were conjugated to horseradish peroxidase and signals were detected using enhanced chemi-luminescence (Amersham, GE). Relevant bands were quantified by densitometry using Image J, background corrected and normalized to actin levels, then compared to time matched controls.
RNA isolations were performed in parallel across irradiated, bystander and sham-irradiated samples, so that all samples were collected from one sub-cultivated pool of IMR-90 cells that were seeded from a single cryo-vial. After treatment, cells were lysed in pools from Mylar dishes at 30-minute, 1-hour, 2-hour, 4-hour, 6-hour and 24-hour time points. We repeated the experiment four times to provide four biological replicates. We analyzed the 30-minute RNA pools by microarray hybridization. Cyanine-3 (Cy3) labeled cRNA was prepared from 0.3 μg RNA using the One-Color Low RNA Input Linear Amplification PLUS kit (Agilent). Dye incorporation and cRNA yield were monitored with the NanoDrop ND-1000 Spectrophotometer (Thermo Scientific). 1.5 μg of cRNA (>9 pmol Cy3 per μg cRNA) was fragmented, hybridized to Agilent Whole Human Genome Oligo Microarrays (G4112F) using the Gene Expression Hybridization Kit, and washed following Agilent's recommendations. Slides were scanned with the Agilent DNA Microarray Scanner (G2505B). Default parameters of Feature Extraction Software 9.1 (Agilent) and grid version 014850_D_F_20090416 were used for image analysis, data extraction, background correction, and flagging of non-uniform features.
Background corrected intensities were log2 transformed and median-normalized in BRB-Array Tools, Version 3.8.0 . Non-uniform outliers or features not significantly above background intensity in 40% or more of the hybridizations were filtered out, leaving 27,576 features. A further filter requiring a minimum 1.5-fold change in at least 20% of the hybridizations was then applied yielding a final set of 6911 features that were used for subsequent analyses. The microarray data are available through the Gene Expression Omnibus database using accession number GSE18760.
BRB-Array Tools was used to identify genes that were differentially expressed between controls and directly or bystander irradiated cells using a random-variance paired t-test, an improvement over the standard t-test that permits sharing information among genes about within-class variation without assuming that all genes have the same variance . The test compares the differences in mean log-intensities between classes relative to the expected variation in mean differences computed from the independent samples. Genes with p-values less than 0.005 were considered statistically significant. The false discovery rate (FDR) was also estimated for each gene using the method of Benjamini and Hochberg , to control for false positives.
The High-Capacity cDNA Archive Kit (Life Technologies, Foster City, CA) was used to prepare cDNA from total RNA. A custom low-density TaqMan array (Life Technologies, Foster City, CA) was designed using validated assays (additional file 1). Genes for inclusion on the low-density array (LDA) were selected on the basis of differential expression and low FDR, and seven previously selected endogenous control genes  were also included. For gene validation studies, 100 ng cDNA was used as input for LDAs. Quantitative real time PCR reactions were performed with the ABI 7900 Real Time PCR System using Universal PCR Master Mix (Life Technologies), with initial activation at 50°C for 120 seconds and 94.5°C for 10 minutes followed by 40 cycles of 97°C for 30 seconds and 59.7°C for 60 seconds.
Relative fold-inductions were calculated by the ΔΔCT method as previously used  using SDS version 2.3 software (Life Technologies). We applied geNorm  to the seven endogenous control genes on the LDAs to determine the most appropriate genes for normalizing the results. The LDA data was normalized to the geometric mean of peptidylprolyl isomerase A (PPIA) and ubiquitin C (UBC) gene expression levels.
The genes responding significantly (p < 0.005) to either direct alpha particle or bystander irradiation were imported into DAVID, the database for annotation, visualization and integrated discovery http://david.abcc.ncifcrf.gov/. The genes/proteins in our list were mapped to DAVID identifiers, and then functionally annotated using the DAVID biological processes and molecular function categories. The number of genes in each functional classification category was compared against the number of genes from the NCBI human genome in that category. The one-tailed Fisher exact t-test probability value was used to statistically determine over- or under- representation of classification categories, Bonferroni corrected p-values less than 0.05 were considered significant.
The sets of genes significantly responding to direct or bystander irradiation (p < 0.005) were also imported into Ingenuity Pathways Analysis (IPA) (Ingenuity® Systems, http://www.ingenuity.com) to analyze network interactions between the genes. The imported genes were mapped onto a global molecular network developed from information contained in the Ingenuity Pathways Knowledge Base. Networks of these significantly responding genes were then algorithmically generated based on their connectivity. The biological functions that were most significant to these networks were determined, and Fischer's exact test was used to calculate p-values determining the probability that each biological function assigned to a network was due to chance alone. We also identified the IPA canonical pathways that were most significant within the differentially expressed gene sets.
In four independent experiments, RNA was extracted from control, directly irradiated, and bystander IMR-90 fibroblasts 30 minutes after exposure and was hybridized to Agilent Human whole genome arrays. Using the class comparison tool of BRB-Array Tools , we identified genes with significantly different expression compared with controls. In directly irradiated cells, 709 genes were differentially expressed (p < 0.005 and false discovery rate (FDR) < 5%; additional file 2). In bystander cells 407 genes were differentially expressed (p < 0.005 and FDR < 8.5%; additional file 3). 293 genes responded significantly to both direct and bystander irradiation.
Gene ontology analysis on biological processes
PANTHER Biological Processes
7.63 × 10-4
2.24 × 10-2
BP00103:Cell surface receptor mediated signal transduction
9.22 × 10-3
4.38 × 10-3
BP00285:Cell structure and motility
1.05 × 10-2
2.98 × 10-4
2.00 × 10-2
6.16 × 10-8
2.00 × 10-2
1.51 × 10-2
BP00044:mRNA transcription regulation
5.66 × 10-7
6.18 × 10-4
BP00104:G-protein mediated signaling
8.51 × 10-4
BP00020:Fatty acid metabolism
2.32 × 10-3
1.92 × 10-2
BP00292:Other carbon metabolism
2.66 × 10-2
2.84 × 10-2
2.88 × 10-2
3.00 × 10-2
4.60 × 10-2
Gene ontology analysis on molecular functions.
PANTHER Molecular Functions
MF00213:Non-receptor serine/threonine protein kinase
8.20 × 10-3
2.39 × 10-2
MF00261:Actin binding cytoskeletal protein
4.62 × 10-2
2.05 × 10-4
MF00231:Microtubule binding motor protein
3.23 × 10-3
3.62 × 10-3
4.49 × 10-3
MF00264:Microtubule family cytoskeletal protein
2.70 × 10-2
MF00224:KRAB box transcription factor
4.08 × 10-2
MF00262:Non-motor actin binding protein
4.10 × 10-2
In this study we focused on early events to better understand the early signaling that shapes the radiation bystander effect. Time dependent analyses of specific genes such as PTGS2, IL8 and BCL2A1 had previously shown dramatic biphasic changes with a large response as early as 30 minutes after exposure . This suggested that the re-programming of gene transcription is initiated prior to 4 hours after exposure. We also wanted to investigate signaling pathways in order to identify candidates for the upstream initiating signals that may lead to the gene expression changes observed at later time points. Therefore, we investigated the early burst of gene expression, looking for potential signaling nodes using a bioinformatics approach in conjunction with more conventional analyses of protein modifications of transduction pathways.
We measured global transcriptional changes at 30 minutes after treatment and found even more differentially expressed genes (709 in irradiated cells and 407 in bystanders) than we had previously reported at 4 hours, when 197 genes were changed significantly in irradiated cells and 137 genes were changed significantly in bystander cells . The general patterns of gene expression were confirmed by qRT-PCR. Many genes that showed increases in expression after 4 hours were also induced at the 30-minute time in both irradiated cells (127 genes) and bystanders (120 genes). In contrast to the prominent role of p53 previously observed in the 4-hour response, however, there was no such pronounced p53 response at 30 minutes (Figure 1a) in the irradiated cells.
Of the 709 genes affected in directly irradiated cells, 416 genes were not significantly altered in bystanders. Gene ontology analysis using DAVID indicated that the top functional categories of the direct-exposure-only genes were in mRNA transcription regulation and activation of gene expression. These functions were not significant in bystanders at early times indicating that irradiation rapidly induced processes that altered the transcriptional landscape within hit cells, in addition to the overall stress response characterized by the release of cytokines and inflammatory regulators. The early bystander response was associated with cytokine signaling, cell structure changes and extra-cellular communication, suggesting that bystander fibroblasts were responding by altering their own cellular environment (Table 1). We had already reported that the gene expression responses of cyclooxygenase2 (PTGS2), interleukin-8 (IL8) and BCL2 related protein A1 (BCL2A1) were comparable at 30 minutes and 4 hours after exposure . These genes are known targets for regulation by the NF-κB transcription factor, and we predicted that the biphasic response in mRNA levels could be due to modulation by inhibitors of NF-κB. In this study, we have verified that this was a frequent pattern of gene expression. Many genes that responded with an increase in mRNA levels at 4 hours showed a strong response at 30 minutes as well (Figures 1 and 2), suggesting co-ordinate regulation of genes with this biphasic pattern. These genes did not include canonical p53-response genes, which generally show peak induction at approximately 4 hours after irradiation, depending on the model system used. We further verified time dependent patterns of mRNA changes for 11 additional genes across the 24 hour time interval and found that they all showed peaks of induction after 30 minutes and 4 or 6 hours with a sharp decrease to control levels at the time between 1 and 2 hours, Figure 2a. This pattern of modular gene expression could be attributed to co-regulation of genes under the control of a transcription factor, such as NF-κB, acting as an early signal transducer [3, 29].
We compared the two gene sets, 709 genes in irradiated and 407 genes from bystander, for significant biological functions in IPA, which categorizes genes by mechanism- and disease-related associations. Gene expression changes in both irradiated and bystander cells appeared to be relatively similar in top-scoring biological functions categories such as cancer, cellular growth, proliferation, and cell death, which are very broad descriptions for functions associated with regulation of cell numbers and tissue homeostasis (additional file 4). 301 cancer-related genes in irradiated cells gave a p-value range from 1.2 × 10-17 to 5.08 × 10-4 for the cancer functional category. For 199 differentially expressed genes from bystanders; the p-value range for the same category was 2.45 × 10-16 to 2.36 × 10-4. This indicated similar enrichment of this biological response in both gene sets. However, a closer look at the individual genes in the two conditions revealed differences, which were expected because direct irradiation would cause cells to have a more pronounced DNA damage and cell cycle response. IPA grouped 70 genes from 407 bystander responding genes under the functional category "cell-to-cell signaling" In IPA these genes were sub-grouped further by specific biological activities related to cell-to-cell signaling. We ranked these sub-categories by decreasing order of number of genes affected and the highest ranked functions in bystander cells correlated with increases in activation, adhesion, communication, signaling, binding and stimulation of normal cells. From IPA and ontology analyses taken together, the mRNA changes suggest stimulation of cellular response in bystanders by signaling molecules as early as 30 minutes after treatment (additional file 4).
Pathway analysis of genes affected in irradiated cells after 4 hours suggested that in addition to p53 and NF-κB transcription factors, β-catenin (most probably in concert with LEF1) could be an important trans-activator of gene expression. We chose to focus on the AKT-GSK3β-β-catenin axis of signaling in normal fibroblasts because of its central role in response to growth factors, changes in the medium, cytokine signaling, cell structure changes and stress due to production of reactive radicals, all of which are components of the bystander response . Previous studies on bystander cells after carbon ion or alpha particle irradiation have also implicated activation of Akt signaling based on up-regulation of target gene [9, 31]. The activation of β-catenin as a nuclear activator of transcription is known to follow GSK3β phosphorylation by Wnt signaling or alternatively to be regulated through PKB/AKT activation [30, 32]. Cells responded to irradiation by rapid activation of AKT via phosphorylation and concomitant inactivation by phosphorylation of GSK3β at 1 to 4 hours. In bystanders, there was a similar pattern of GSK3β inactivation, but with a lag of 30 minutes or more. Cytoplasmic levels of phosphorylated β-catenin decreased gradually, which is considered a good indicator of trans-location to the nucleus  and activation of gene expression. This effect was also observed in bystanders, so β-catenin activity as a transcriptional co-activator of gene expression may occur in both populations. We have summarized these findings in Figure 3, which is a proposed pathway model of AKT-GSK3β-β-catenin signaling pathway after irradiation. External signals such as growth factors, interleukins 6 and 8, which were induced early and secreted by irradiated cells could trigger activation of AKT in all cells  initiating the signaling cascade outlined in Figure 3, which results in gene expression activation of β-catenin target genes such as matrix metalloproteinase 3 (MMP3), Fas (FAS) and metallothionein 1L (MT1L). In irradiated cells at 4 hours, this transcriptional regulator potentially regulated expression of MT1L (+3.9 fold change), DKK1 (+3.0 fold change), MMP3 (+5.0 fold change), LAMC2 (+2.4 fold change) and FAS (+2.3 fold change) (Figure 3b). In bystander cells, only MMP3 (+4.0 fold change) and MT1L (+3.8 fold change) were detected as differentially expressed at 4 hours. Of the genes predicted to be β-catenin targets, only MMP3 showed the predicted NF-κB biphasic response that was similar in both irradiated and bystander cells. Although there is overlap and redundancy in the roles of both transcriptional regulators and target genes, our results suggest that this mechanism is active in both irradiated and bystander cells. The apparent lag between activation of this pathway in directly irradiated and bystander cells suggests that although it is clearly not the initiating signal, this pathway may play a role in the maintenance and later development of the bystander state. There is evidence for a lag in signaling between irradiated and bystander cells as observed by apoptosis induction in human fibroblasts , γH2AX foci induction after media transfer  and mutation induction in bystanders after treatment with conditioned media . The mechanism driving the burst of cytokine and other gene expression seen in bystander cells at 30 minutes post-exposure has yet to be clearly elucidated.
We also observed some changes in AKT expression in our study. Although there was an overall increase in phosphorylation in both irradiated and bystander cells, total AKT protein levels were slightly decreased (Figure 4a). This could be a consequence of the significant down-regulation of AKT2 mRNA observed at 30 minutes (0.4 ± 0.1 fold change) (additional files 2 and 3) and 4 hours (0.5 ± 0.1 fold change) in both irradiated cells and bystanders. The AKT antibody we used detected all three isoforms of AKT, and the phospho-Ser473-AKT antibody recognized phosphorylation of Ser473 on AKT1 as well as the corresponding residues in AKT2 and AKT3. However, GSK3β is known to be a substrate for all AKT isoforms , and in our system inactivation of GSK3β could result from activation of all AKT isoforms. Further studies will be required to distinguish the role of decreased AKT2 mRNA specifically, and its contribution to phosphorylation and inactivation of GSK3β.
The importance of our study approach is that whole genome bio-informatics can help elucidate novel signaling networks that contribute to any phenotype. Exploring gene expression and signaling after irradiation led us to propose activation of an important signal transduction module, AKT signaling via β-catenin, as part of the cellular response in bystanders. Further verification of the role of β-catenin in bystanders will help us understand the complex nature of this response. In the field of radiation biology, there are few studies that use this approach, but a recent study on keratinocytes exposed to 1 cGy X-rays that identified the GATA3 transcriptional factor as a critical regulator of cellular radiation response at the genome level  is one example. Importantly, our study is the first to investigate the global transcriptional response within the first half hour after exposure with the goal of identifying signals more proximal to the generation of the bystander signal. This could have implications for the understanding of how cells respond in vivo and lead to understanding the rapid nature of the response.
A rapid and widespread transcriptional response occurs following irradiation in human fibroblast cells, and is communicated rapidly to bystanders. We found that cell structure and cell-cell communication processes are triggered very quickly in irradiated cells and these changes precede cell cycle gene responses that occur at a later time. Timing is important because release of free radicals and cytokines occurs quickly after irradiation and identifying early signaling events is critical to a full understanding of the range and extent of changes in bystander cells. The wave-like response of many genes following radiation suggests coordinate control of a large number of genes especially those that are controlled by NFκB. Activation of AKT-GSK3β-β-catenin signaling early after irradiation and activation of target genes later in irradiated and bystander cells is a novel finding implicating a new pathway in bystander response. The connection between this signal transduction module, with the potential to converge extra-cellular communication with intra-cellular alterations of proteins and gene expression adds one more piece to the puzzle of stress response in a cell population where not all cells experience a direct ionization event.
reactive oxygen species
tumor necrosis factor alpha
nuclear factor kappa-light-chain-enhancer of activated B cells
v-akt murine thymoma viral oncogene homolog
glycogen synthase kinase 3 beta
Biometric research branch
false discovery rate
quantitative real time Reverse Transcription-Polymerase Chain Reaction
The database for annotation, visualization and integrated discovery.
Analyses were performed using BRB-Array Tools developed by Dr. Richard Simon and Amy Peng Lam. We would like to thank Mr. Stephen Marino and the RARAF staff for assistance with irradiations. The Radiological Research Accelerator Facility (RARAF) is an NIH supported Resource Center through NIBIB grant EB-002033. This work was supported by NIH grant CA-49062.
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.