Identification of gene-based responses in human blood cells exposed to alpha particle radiation
BMC Medical Genomics volume 7, Article number: 43 (2014)
The threat of a terrorist-precipitated nuclear event places humans at danger for radiological exposures. Isotopes which emit alpha (α)-particle radiation pose the highest risk. Currently, gene expression signatures are being developed for radiation biodosimetry and triage with respect to ionizing photon radiation. This study was designed to determine if similar gene expression profiles are obtained after exposures involving α-particles.
Peripheral blood mononuclear cells (PBMCs) were used to identify sensitive and robust gene-based biomarkers of α-particle radiation exposure. Cells were isolated from healthy individuals and were irradiated at doses ranging from 0-1.5 Gy. Microarray technology was employed to identify transcripts that were differentially expressed relative to unirradiated cells 24 hours post-exposure. Statistical analysis identified modulated genes at each of the individual doses.
Twenty-nine genes were common to all doses with expression levels ranging from 2-10 fold relative to control treatment group. This subset of genes was further assessed in independent complete white blood cell (WBC) populations exposed to either α-particles or X-rays using quantitative real-time PCR. This 29 gene panel was responsive in the α-particle exposed WBCs and was shown to exhibit differential fold-changes compared to X-irradiated cells, though no α-particle specific transcripts were identified.
Current gene panels for photon radiation may also be applicable for use in α-particle radiation biodosimetry.
Nuclear terrorism is a global concern with illicit trafficking events involving nuclear material on the rise . Between 1993 and 2011, the International Atomic Energy Association (IAEA) documented 2164 nuclear material incidents or malicious acts and 588 involved the theft or loss of nuclear or radioactive materials. A further 18 of these involved plutonium or highly enriched uranium. Such events highlight the potential for radioactive material to fall into the wrong hands and potentially being used for the fabrication of a radiological dispersal device (RDD) . Los Alamos National Laboratory has conducted a thorough review of RDD source material and has postulated that four of the nine isotopes most likely to be employed are α-particle emitters, primarily due to their minimal shielding requirements and ease of concealment . There is also the concern of long-term contamination due to these isotopes’ long half-lives and the severe biological damage that can occur from a minimal dose of exposure.
In recent years, much work has gone towards developing strategies for radiation biodosimetry with a specific focus on photon radiation [4–8]. However, the strategies currently employed for photon radiation may not provide adequate dose estimates for α-particle exposures. Unlike photon radiation, α-particles travel a short distance (40-70 μm) and create very dense ionizing tracks as they traverse a medium. They typically cause an energy deposition of 160 keV•μm-1 for 2.5 MeV α-particles in comparison to 2.0 keV•μm-1 for low linear energy transfer (LET) X-rays [9, 10]. Therefore, α-particles produce more significant biological effects when compared to equal absorbed doses from low-LET photon radiation, which are more sparsely ionizing [11–13]. This difference in ionization density may provide a means of distinguishing radiation type based on the magnitude of the biological response.
Development of biomarker-based biodosimetry has been put forth as one of the key priority development areas for nuclear threat countermeasures , and microarray data/gene based profiling has served as timely and minimally invasive means to address this priority area . There have been several studies examining the gene expression profiles of human cells using functional genomics platforms for photon radiation [16–18]. However, the availability of similar gene tools for high LET radiation types, such as α-particles, remains limited. To date, the majority of α-particle transcriptional studies have been performed in vitro using transformed or normal cell types [19–21]. There has also been a selected few studies that have profiled genomic changes and compared the responses following exposure of cells to different radiation types [22, 23]. Microarray studies in our own lab using epidermal keratinocytes exposed to both α-particle and X-ray radiation have also shown transcriptional differences between these radiation types .
In the present study, genomic strategies were employed to identify biomarkers of α-particle radiation exposure. Circulating peripheral blood mononuclear cells (PBMCs) were isolated from normal, healthy volunteers and exposed to α-particle radiation. Twenty-four hours post-exposure, the expression of transcripts was assessed using Illumina bead array technology and these responses were compared to non-irradiated controls. Dose-responsive genes were then further assessed in independent white blood cell (WBC) populations exposed to either α-particles or X-rays.
All procedures were approved by Health Canada’s Research Ethics Committee and a flow chart delineating the experimental sequence is outlined in Figure 1. Briefly, peripheral blood from healthy, non-smoking volunteers was drawn via periphery venipuncture with informed consent from all subjects into either 5 × 10 ml EDTA (for gene analysis) or 2 × 4 ml lithium heparin (for plasma analysis) vacutainer tubes (Becton Dickinson and Company, Franklin Lakes, NJ). A total of 6 male and 6 female donors participated. Before any further processing, a 100 μl whole blood sample was used for a complete blood count (CBC) via automatic haemocytometer (Beckman Coulter, Mississauga, ON).
Peripheral blood mononuclear cell isolation
PBMCs were isolated from whole blood for an initial global screening of gene transcripts using microarray technology. A similar isolation procedure was employed as described by Boyum . Briefly, 15 ml of Histopaque-1077 sucrose gradient (Sigma-Aldrich, MO, USA) was pipetted into the upper chamber of an Accuspin Tube (Sigma-Aldrich). The tube was centrifuged (800 × g) for 30 sec to ensure that the Histopaque was below the frit layer. Freshly drawn whole blood was pipetted into the upper chamber of tube. The tube was then centrifuged 800 × g for 15 minutes. The band of mononuclear cells was transferred to an alternate centrifuge tube and washed with 10 ml of isotonic phosphate buffered saline (PBS) three times. Pelleted cells were then resuspended in RPMI-1640 media supplemented with 10% fetal bovine serum (FBS), 2 mM L-glutamine and 100 U penicillin & 100 μg streptomycin/ml (Sigma-Aldrich).
Total white blood cell isolation
Following initial experiments on microarray analysis of PBMCs, further studies using qPCR were conducted on complete WBC populations. WBCs were isolated from whole blood using Histopaque-1119 (Sigma-Aldrich). Twenty-five milliliters of whole blood was gently poured over 15 ml of Histopaque-1119 and spun at 1000 × g for 10 minutes. This resulted in erythrocyte sedimentation below the Histopaque gradient and a total white blood cell population above the gradient. This volume (~15 ml) was then transferred to a new 50 ml tube and diluted 1:2 with phosphate buffered saline (PBS). The resulting WBC pellet was then washed twice more with 10 ml PBS and resuspended in RPMI-1640 media supplemented with 10% FBS, 2 mM L-glutamine and 100U penicillin & 100 μg streptomycin/ml.
PBMCs and WBC irradiations
Either Isolated PBMCs or WBCs were seeded at total cell density of 8-10 × 106 cells in 2ml of media and were cultured on thin Mylar based plastic dishes (MD) (Chemplex Industries, Palm City, FL, USA), which allowed for penetration of α-particles. Cells were allowed to settle for 10 minutes before performing the irradiations. Irradiations were performed at doses of 0 (control), 0.5, 1.0 or 1.5 Gy using Americium (241Am) electroplated discs with an activity level of 66.0 kBq ± 3% (dose rate of 0.98 ± 0.01 Gy/h, LET of 127.4 ± 0.4 keV/μm). The absorbed dose of α-particle radiation to which the cells were exposed was calculated using the GEANT4 v.9.1 Monte Carlo tool-kit . Cells destined for X-ray radiation at doses of 0 (control) 2, 5 or 10 Gy were exposed using the X-RAD 320 X-ray irradiation system (Precision X-ray, Inc., North Branford, CT, USA) at a higher dose rate of 0.98 ± 0.05 Gy/min. Exposures were performed in duplicate and pooled. Twenty-four hours following irradiation, a 50 μl aliquot of cells was assessed for cellular viability using the Trypan Blue viability assay (Bio-Rad, Hercules, CA), and a 100 μl aliquot was used for a CBC via automatic haemocytometer. The remainder of the cells were spun down and used for RNA extraction.
H2AX phosphorylation assay
H2AX phosphorylation was assessed using flow cytometry following a modified protocol by MacPhail et al. . Thirty minutes after exposure, WBC suspensions (5 × 105 cells per sample) were fixed with 10% formaldehyde (Fisher Scientific, USA) and incubated for 10 min at room temperature. The cells were then washed and re-suspended in 1 ml cold (−40°C), 70% methanol (Fisher Scientific) in 1x PBS and stored at -40°C overnight or up to two weeks. One ml of cold TBS (tris-phosphate buffered saline, 0.0154 M Trizma Hydrochloride (Sigma–Aldrich Canada), 0.5 M NaCl (Fisher Scientific), pH 7.4) was then added to each sample, mixed well, centrifuged (8 min, 400 × g, 4°C) and re-suspended in 1 ml of cold TST (TBS serum triton, 96% TBS, 4% FBS (Sigma–Aldrich), 0.1% Triton X-100 (Sigma–Aldrich)). The samples were incubated on ice for 10 min, centrifuged (5 min, 400 × g, 4°C) and re-suspended in 200 μl of anti-γ-H2AX-fluorescein isothiocyanate (FITC) antibody (Millipore, USA) diluted 1:500 in TST. After 2 h incubation on ice in the dark, 1 mL of TBS with 2% FBS was added. The samples were then centrifuged (5 min, 400 × g, 4°C), re-suspended in 250 μl TBS with 2% FBS. Immediately prior to analysis by flow cytometry, 2 μl of 1 mg/mL propidium iodide (PI) was added to each sample. For flow cytometry analysis, data acquisition was set to analyze 2 × 104 cells from the whole cell population as identified by a forward scatter (FSC) vs. side scatter (SSC) dot plot. All debris under the FSC and SSC threshold were excluded from the analysis. The γ-H2AX response was measured by assessing the increased level of intracellular fluorescence characterized in the cells, as determined by the geometric mean of the intensity peak of the anti-γ-H2AX-FITC (channel number) of the γ-H2AX positive cells. All samples were analyzed on a BD FACSCalibur flow cytometer (BD Biosciences, San Jose, CA, USA).
Twenty-four hours post-radiation exposure or negative control conditions, RNA extractions were performed on either PBMCs or WBCs. In the case of PBMCs, the cells were resuspended in 350 μl of Buffer RLT containing 1% β-Mercaptoethanol (Qiagen’s RNeasy Mini kit; Qiagen Inc, Mississauga, ON) then frozen at -80°C until processed. Frozen lysates were thawed on ice and mixed well by pipetting. The lysate was transferred directly onto a QIAshredder spin column (Qiagen Inc), placed in a 2 ml collection tube and centrifuged for 2 min at ~12,000 g. A volume of 350 μL of 70% ethanol was added. Total RNA was then extracted using the RNeasy Mini kit according to the manufacturer’s instructions (Qiagen Inc), with the addition of Qiagen’s On-Column RNase-free DNase (Qiagen Inc) to eliminate any remaining DNA contamination. In the case of the WBCs, RNA extractions were performed using QIAzol reagent (Qiagen) and following manufacturer’s instructions. Briefly, 700 μl of QIAzol reagent was added to the cell pellet and then homogenized via up-and-down pipetting of the mixture 50 times. After room temperature incubation for 15 minutes, 140 μl of chloroform was added for phase separation. The aqueous layer containing RNA was then removed and precipitated with 100% ethanol. Total RNA was then isolated using the miRNeasy column purification kit. All total RNA sample concentrations and RNA quality were determined using both an Agilent 2100 Bioanalyzer and RNA Nanochips (Agilent Technologies Canada Inc., Mississauga, ON) and spectrophotometrically (OD ratio of A260:A280) using a Nanodrop (Fisher Scientific, Ottawa, ON). All extracted PBMC RNA samples were determined to be of good quality (RNA Integrity Number ≥ 8.0) with minimal degradation and stored at -80°C until further analysis. WBC RNA from α-particle exposed samples was determined to be of good quality (RNA Integrity Number ≥ 9.2) with three samples being excluded from sample analysis due to insufficient RNA yield.
An input of 200 ng of PBMC mRNA was used for whole genome analysis following the Illumina(r) Whole Genome Expression Profiling Assay Guide (11317302 Rev. A). Samples were hybridized on Illumina human-12 v2 RNA BeadChips. BeadChips were imaged and quantified with the Illumina iScan scanner and data was processed with Illumina GenomeStudio v2010.2.8.11.
An input of 200 ng of WBC miRNA expression was profiled using the nCounter system (NanoString Technologies, Seattle, WA) which profiles the expression levels of miRNAs. This was performed using the human miRNA expression assay (version 2) according to manufacturer’s instructions and read using the nCounter digital analyser.
Quantitative real time-polymerase chain reaction (qPCR) validation
Selected genes deemed statistically significant (as described below in the Statistical analysis section) by microarray analysis or nCounter system were further assessed by qPCR. Total RNA (400 ng mRNA and 200 ng miRNA) isolated from cells were reverse transcribed into complementary DNA (cDNA) using the RT2 First Strand Kit (Qiagen) or miScript Kit respectively. Gene profiling was performed according to the manufacturer’s instructions using custom RT2-profiler PCR arrays (Qiagen). Reactions were prepared in 96-well plates and performed using a spectrofluorometric thermal cycler (Biorad iCycler; Hercules, CA). The relative expression of each gene was determined by using the comparative threshold (Ct) method .
Customized gene array panel
A total of 84 unique identifiers were used for the development of a customized 384-well format gene array panel (Table 1). This panel was comprised of genes that were shown by microarray technology to be dose-responsive and also expressed at 1.0 and 1.5 Gy following exposure of PBMC to α-particle radiation. This panel also included a negative control gene (GNG7), housekeeping genes (ACTB, GAPDH, GUSB, B2M) and selected few genes derived from the work of Paul and Amundson  which were not identified in this study as statistically significant by microarray analysis (Table 1). SABiosciences (Qiagen) designed the primers and provided a 384 well-format platform that was compatible for use on the LightCycler 480 real-time PCR system (Roche, Mississauga, ON). A high-throughput PCR platform, comprising the Caliper Zephyr Compact Liquid Handling Station, the Caliper Twister II plate handler (PerkinElmer, Woodbridge, ON) and the Lightcycler 480 was employed with custom protocols developed in the Inhalation Toxicology Laboratory of Health Canada. This system allowed for the screening of 144 samples in a one-week time-span.
Microarrays statistical analysis was performed as follows. Data background correction was done within GenomeStudio (Illumina), and then exported to the lumi R package. Data was then normalized via quantile method, rendering the distribution of probe intensities of each array in a set of arrays equivalent. Normalized data was then log2 transformed for statistical comparisons. Linear models for microarray data (LIMMA) was employed to identify differentially expressed gene signatures between the different exposure conditions for the PBMC microarray, nCounter miRNA and the WBC qPCR datasets [28, 29]. In brief, this method involves fitting a linear model for each gene in the data and moderating the standard error via an empirical Bayes method. This is used to estimate the moderated t-statistics/F-statistics for each gene, shrinking the standard error towards a common value. This test is similar to an ANOVA for each gene with the exception that standard deviations are moderated across genes, allowing more stable inference for each gene. Moderated standard deviations are a compromise between individual genewise standard deviations and overall pooled standard deviations. Multiple comparison false discovery rate (FDR) was evaluated using the Benjamini-Hochberg (BH) method .
PBMC qPCR data was analysed for statistical significance without multiple correction comparison using gene-wise ANOVA as there was a-priori reason for gene analysis. Hierarchical clustering was performed using the WBC qPCR data using dChip (http://www.hsph.harvard.edu/cli/complab/dchip). This software was used to cluster the different exposure conditions by gene signature and group genes by similarity of expression patterns. The distance between genes is measured as 1- r (Pearson correlation coefficient).
To ensure that the WBC suspensions were undergoing irradiation and sustaining DNA damage, a biological assay indicative of DNA damage, the phosphorylation of H2AX, was employed. Thirty minutes post-exposure, cells were assessed for the expression of γ-H2AX, a marker indicative of DNA double strand breaks. A dose-dependent increase in the γ-H2AX signal was observed following exposure to α-particle radiation as seen by a plot of the geometric mean of this signal as a function of dose and the pronounced shift in the curve (Figure 2). Statistically significant responses were obtained at the medium (1.0 Gy) and high (1.5 Gy) doses tested (p < 0.01) relative to the non-irradiated control treatment group. A bi-modal shaped curve was observed at the lowest dose of α-particle radiation which with increasing doses transitioned to a mono-modal curve. At the highest dose (1.5 Gy) an approximate 3-fold increase in γ-H2AX signal was observed relative to the control sample. As a positive control, isolated leukocytes were irradiated with X-rays at a high dose rate (1 Gy/min) and a dose range of 2-10 Gy. A plot of this response indicated X-rays to induce more double strand breaks as seen by the marked increase in γ-H2AX signal with dose of radiation relative to α-particle treated cells.
Genomic profiling was performed on RNA extracted from isolated PBMCs 24 hr post-exposure. In order to identify biomarkers, statistical stringency was prioritized to mine for reliable genes using a Benjamini Hochberg (BH) false discovery rate (FDR) correction. All differentially expressed genes were filtered on flagged spots and a BH FDR corrected p-value <0.05. A numeric summary of the gene responses at each of the doses is provided in Table 2. Overall, there was a pronounced induction of transcriptional response, with the majority of genes being up-regulated in the presence of the radiation insult. Escalating doses induced an increasing number of transcripts with 30, 69 and 137 genes differentially modified at 0.5, 1.0 and 1.5 Gy respectively. A Venn diagram was constructed to provide a quantitative representation of the similarities and differences in expression profiles at each of the doses (Figure 3). Twenty-nine genes were shown to be differentially expressed at all three doses with expression levels ranging from 2-10 fold. Thirty-nine genes were common between differentially expressed gene sets at both the medium (1.0 Gy) and high (1.5 Gy) dose. The range in expression levels of these 68 genes is summarized as a heat map which delineates the genes by degree of fold change (Figure 4).
All 68 genes observed to be differentially expressed were further validated using qPCR. A comparison of the responses using the two technologies showed a similar trend in differential expression levels. As shown in Table 3, all genes that exhibited a significant response across the 3 doses using microarray analysis were also observed to exhibit a similar trend using qPCR. However, approximately 20% of the total validated targets were shown to be non-significant using qPCR. This is not unexpected as in comparison to microarrays, qPCR may report different statistical assignments . In contrast, there were also a sub-set of 9 genes which showed statistical significance at the 0.5 Gy dose via qPCR but not using microarray analysis.
Custom qPCR panel
A customized gene array panel was developed to more cost-effectively assess whether the differentially responding genes observed from the microarray dataset were expressed in the WBC population, remained differentially expressed following exposure to α-particle radiation and were responsive to X-ray exposures.
WBCs were isolated from 12 healthy individuals and exposed to α-particles and X-rays. Total white blood cell counts were typically in the range of 5-10 × 106 cells/mL. The viability of the cells was assessed using the Trypan Blue viability assay pre- and post-irradiation. The cells remained viable (above 98%) and no significant changes in blood cell counts or populations subsets were observed post-irradiation relative to unirradiated cells (Table 4). Twenty-four hours after irradiation, RNA was extracted and reverse transcribed to cDNA. A comparison of the differential gene responses obtained from the qPCR of WBCs and microarray analysis of PBMCs showed similar fold change and statistical significance for the majority of the 68 transcripts that were assessed (Table 5). These responding genes were compared to those obtained for X-ray exposed cells. Box-plots of the responding genes from the WBC qPCR dataset allowed for a visual comparison of the two radiation types and the range in inter-individual variability between transcripts (Figure 5). Overall, all genes responsive in α-particle treated cells were also observed to be expressed in X-irradiated cells. The data displayed minimal variability between control treatment groups under varied radiation exposure conditions. Furthermore, the majority of genes displayed dose-response trends for both α-particle and X-ray radiation. Hierarchical clustering (Figure 6) was further used to display groupings and make class distinctions. The two control groups clustered together as expected and showed a distinct trend relative to the other exposed groups. The lowest dose of radiation (0.5 Gy α-particle) also clustered with the unexposed groups. The next largest classification comprised the remaining exposure groups, in which the 2 Gy X-ray, 1.0 and 1.5 Gy α-particle and 5 and 10 Gy X-ray were classified further from the controls in order of descending similarity respectively. The subsequent clustering of the 1.0 and 1.5 Gy α-particle exposure and the 5 and 10 Gy X-ray exposure together suggests that it is possible to make distinctions between high X-ray radiation doses and α-particle doses using a clustering algorithm.
The screening of ~800 miRNA transcripts using NanoString profiling resulted in a minimal number of responding targets. Only one miRNA (miR-34a) was observed to be differentially modulated (p < 0.05) across all doses following α-particle radiation exposure in WBC. This target was dose-responsive and subsequently validated using qPCR (Table 6). miR-34a was up-regulated over 2.5 fold in all exposed samples and had a similar 3 fold induction in the 1.0 and 1.5 Gy α-particle doses. This target was not specific to α-particle exposure as expression was also observed with X-ray irradiation at all three doses tested.
The overarching goal of this research was to identify genes responsive to α-particle radiation exposure for the purposes of developing effective triage tools for use in a population radiation exposure scenario. To date, the majority of genomic-based biomarker radiation studies have focused on photon radiation. Although there is a large body of work concerning radiation exposure and cytogenetic end-points, it was postulated that α-particle radiation may elicit differential cellular response due to its characteristic physical properties, which differ from photon radiation. This may potentially provide more accurate dose estimates for exposures and allow for differentiation between radiation types.
Circulating blood cells represent a sensitive target for early radiation damage and are easily accessible. Isolated leukocytes from healthy individuals were ex-vivo irradiated at a dose range of 0-1.5 Gy at 0.98 Gy/h. These doses were selected based on their relevance to an actual radiological dispersal device scenario , where the dose deposition is approximately 0.5 Gy per α-particle track . Furthermore, previous studies from our laboratory have shown observable biological damage at this dose-range and a time-point of 24 hours post-exposure [34–36]. The X-ray exposure doses and dose rates were selected based on the work by Paul and Amundson . Although whole blood was employed by Paul and Amundson, the limitations of our exposure system only allowed for the use of isolated leukocytes. It was observed that the red blood cells in whole blood minimized the transversal of α-particles through other cell types as measured by DNA double strand breaks via γ-H2AX (data not shown). Leukocytes exposed to α-particles displayed lower γ-H2AX intensities than the X-ray exposed samples possibly due to several factors. Firstly, the samples irradiated with X-rays were exposed to higher overall doses which increase the probability of inducing DNA double strand breaks. Secondly, the dose rates of exposure were markedly different. The α-particle exposure system is limited to a dose rate of ~1 Gy/hr and the X-ray exposure was performed at a dose rate of ~1 Gy/minute. This means the α-particle doses were delivered over a protracted amount of time relative to the X-rays. It has been previously documented that γ-H2AX foci reach a peak 30 minutes after exposure and then are resolved as repair is induced . The time-scale of the higher doses (1.0, 1.5 Gy) of α-particle radiation exposure is between ~1-1.5 hours. Thus, there is a degree of repair occurring as the cells undergo irradiation and the measured signal is lowered due to resolved H2AX foci. There is an eventual equilibrium between induction and repair; it would be expected that the signal intensity would be higher for an acute exposure of an equivalent dose. All of the X-ray exposures were conducted on the order of minutes, with the highest (10 Gy) dose approximating 10 minutes.
A preliminary microarray screening following the exposure of isolated PBMC to α-particles identified 29 genes responsive at all doses tested (0.5-1.5 Gy) and 39 which responded only to the medium (1.0 Gy) and high (1.5 Gy) dose. These genes were validated using qPCR and shown through pathway analysis to be associated with signaling pathways centered around p53 and GADD45A, consistent with a DNA damage response (data not shown). To confirm the validity of this gene set in a more physiologically relevant population of cells, the complete white blood cell population was harvested. Furthermore, the ability of this gene panel to discriminate between radiation qualities was concurrently assessed using X-ray exposures. For this purpose, a customized gene panel was constructed using genes identified as significantly modulated by microarray analysis and included a selected few identified in the literature as being X-ray-responsive, primarily from the work conducted by Paul and Amundson .
The customized gene panel confirmed the validity of our initial microarray results. Primarily, genes differentially responsive in PBMC were also observed to be significantly expressed in the total WBC population using qPCR. However, a selected few genes were shown to lack statistical significance at 0.5 Gy, most likely due to the use of stringent statistical analysis accounting for multiple statistical comparisons and FDR correction using Benjamini-Hochberg testing . When no FDR correction was employed, the qPCR dataset were more comparable to the PBMC microarray results. This is not unexpected as in comparison to microarrays, qPCR may report different statistical assignments [31, 38].
Further visualization of the data using box plots allowed for an assessment of the gene responses with respect to individual variability. All control treatment groups displayed low inter-individual variability for the majority of differentially expressed transcripts, particularly between radiation types, highlighting the potential for these transcripts to be strong biomarkers. Furthermore, the majority of genes displayed dose-response trends for both α-particle and X-ray radiation. This is further highlighted in the hierarchical clustering of the dataset. Inputting all qPCR data resulted in the classification of treatment groups based on exposure doses. Clustering of the 1.0 and 1.5 Gy α-particle exposures together and the 5 and 10 Gy X-ray exposure suggests potential for distinctions to be made between high X-ray radiation doses and comparably lower α-particle doses.
Overall, the gene-expression changes induced by α-particle radiation were not distinct from the X-ray responses. Although all genes modulated in the α-particle exposed WBCs were also observed after X-ray exposure, there were differing degrees of induction. In a selected few genes, α-particle doses of 1.0 Gy at a dose rate of ~1 Gy/hr were able to elicit the same fold induction as a 10 Gy X-ray dose at ~ 1 Gy/min. Thus, despite there being an order of magnitude difference in dose and a greater than fifty fold difference in dose rate between the radiation exposures, a similar cellular response was observed. It has been previously established that the lesions caused by α-particle tracks display different repair kinetics and fidelity . Only a 24 hour time point was examined in this study and there may be pronounced temporal differences in gene expression resulting from the differing degrees of damage and repair between the radiation types.
To date, Turtoi et al.,  is the only other group to examine α-particle radiation induced genomic-wide transcriptional effects in isolated blood cells. This group employed a harvest time of 1hr post-irradiation using a dose range of 0.05 - 1.6 Gy of α-particle radiation. Three hundred and thirty nine genes were shown to be differentially modified with 54% up-regulated and 46% down-regulated. In comparison, the present study identified fewer genes, the majority of which were upregulated. These differences may be attributed to experimental conditions, as Turtoi et al., used varying dose-rates and a post-irradiation harvest time of 1 hour. As well, their gene responses were obtained from only one individual, whereas the present study used 12 different donors.
In summary, genomic strategies were employed for the identification of gene-based responses in PBMCS and WBCs exposed to α-particle radiation. Genomic screening of PBMCs exposed to α-particle radiation identified twenty-nine transcripts that responded at 0.5, 1.0 and 1.5 Gy and thirty-nine genes were shown to be differentially modulated at exposures of 1.0 and 1.5 Gy. Subsequent comparison using WBCs with high dose-rate X-ray radiation showed that both radiation types elicited similar gene responses with varying degree of fold induction. No α-particle exclusive gene modulations were identified. Therefore, current gene panels for photon radiation may also be applicable for use in α-particle radiation biodosimetry. Future work includes testing the gene panel in an in vivo environment, using radiotherapy patients undergoing either total body irradiations or α-particle radiation therapy.
Fetal bovine serum
Triphosphate buffered saline
TBS serum triton
Phosphate buffered saline
Royal Park Medical Institute
Analysis of variance
Quantitative real-time polymerase chain reaction
Radiological dispersal devices
Linear energy transfer
Ingenuity pathway analysis
Complete blood counts.
Bhattacharjee Y: An unending mission to contain the stuff of nuclear nightmares. Science. 2010, 328 (5983): 1222-1224.
Kazi R: Pakistan’s HEU-based nuclear weapons programme and nuclear terrorism: a reality check. Strateg Anal. 2009, 33 (6): 861-876.
Van Tuyle G, Strub T, O’Brien H, Mason C, Gitomer S: Reducing RDD concerns related to large radiological source applications. 2003, Los Alamos, NM: Los Alamos National Laboratory, vol. LA-UR-03-6664
Sullivan JM, Prasanna PGS, Grace MB, Wathen LK, Wallace RL, Koerner JF, Coleman CN: Assessment of biodosimetry methods for a mass-casualty radiological incident: medical response and management considerations. Health Phys. 2013, 105 (6): 540-554.
Amundson SA, Bittner M, Meltzer P, Trent J, Fomace AJ: Biological indicators for the identification of ionizing radiation exposure in humans. Expert Rev Mol Diagn. 2001, 1 (2): 211-219.
Ainsbury EA, Bakhanova E, Barquinero JF, Brai M, Chumak V, Correcher V, Darroudi F, Fattibene P, Gruel G, Guclu I, Horn S, Jaworska A, Kulka U, Lindholm C, Lloyd D, Longo A, Marrale M, Monteiro Gil O, Oestreicher U, Pajic J, Rakic B, Romm H, Trompier F, Veronese I, Voisin P, Vral A, Whitehouse CA, Wieser A, Woda C, Wojcik A, et al: Review of retrospective dosimetry techniques for external ionising radiation exposures. Radiat Prot Dosim. 2011, 147 (4): 573-592.
Lloyd DC, Edwards AA, Moquet JE, Guerrero-Carbajal YC: The role of cytogenetics in early triage of radiation casualties. Appl Radiat Isot. 2000, 52 (5): 1107-1112.
De Lemos Pinto MMP, Santos NFG, Amaral A: Current status of biodosimetry based on standard cytogenetic methods. Radiat Environ Biophys. 2010, 49 (4): 567-581.
Tabocchini MA, Campa A, Dini V: DNA and cellular effects of charged particles. Health Phys. 2012, 103 (5): 547-555.
Dingfelder M: Track-structure simulations for charged particles. Health Phys. 2012, 103 (5): 590-595.
Dingfelder M: Track-structure simulations for charged particles. Health Physics. 2012, 103 (5): 590-595.
Jostes RF: Genetic, cytogenetic, and carcinogenic effects of radon: a review. Mutat Res Rev Genet Toxicol. 1996, 340 (2–3): 125-139.
Hall EJ, Hei TK: Genomic instability and bystander effects induced by high-LET radiation. Oncogene. 2003, 22 (45 REV. ISS. 5): 7034-7042.
Pellmar TC, Rockwell S, and the Radiological/Nuclear Threat Countermeasures Working Group: Priority list of research areas for radiological nuclear threat countermeasures. Radiat Res. 2005, 163 (1): 115-123.
Chaudhry MA: Biomarkers for human radiation exposure. J Biomed Sci. 2008, 15 (5): 557-563.
Paul S, Amundson SA: Development of gene expression signatures for practical radiation biodosimetry. Int J Radiat Oncol Biol Phys. 2008, 71 (4): 1236-1244.
Tucker JD, Joiner MC, Thomas RA, Grever WE, Bakhmutsky MV, Chinkhota CN, Smolinski JM, Divine GW, Auner GW: Accurate gene expression-based biodosimetry using a minimal set of human gene transcripts. Int J Radiat Oncol Biol Phys. 2014 Mar 15, 88 (4): 933-939.
Omaruddin RA, Roland TA, Wallace HJ, Chaudhry MA: Gene expression as a biomarker for human radiation exposure. Hum Cell. 2013 Mar, 26 (1): 2-7.
Ghandhi SA, Yaghoubian B, Amundson SA: Global gene expression analyses of bystander and alpha particle irradiated normal human lung fibroblasts: synchronous and differential responses. BMC Med Genet. 2008, 1 (1): 63.
Kalanxhi E, Dahle J: Transcriptional responses in irradiated and bystander fibroblasts after low dose α-particle radiation. Int J Radiat Biol. 2012, 88 (10): 713-719.
Chauhan V, Howland M, Mendenhall A, O’Hara S, Stocki TJ, McNamee JP, Wilkins RC: Effects of alpha particle radiation on gene expression in human pulmonary epithelial cells. Int J Hyg Environ Health. 2012a, 215 (5): 522-535.
Danielsson A, Claesson K, Parris TZ, Helou K, Nemes S, Elmroth K, Elgqvist J, Jensen H, Hultborn R: Differential gene expression in human fibroblasts after alpha-particle emitter 211At compared with 60Co irradiation. Int J Radiat Biol. 2013, 89 (4): 250-258.
Chauhan V, Howland M, Greene HB, Wilkins RC: Transcriptional and secretomic profiling of epidermal cells exposed to alpha particle radiation. Open biochem J. 2012, 6: 103.
Böyum A: Isolation of leucocytes from human blood. Further observations. Methylcellulose, dextran, and ficoll as erythrocyteaggregating agents. Scand J Clin Lab Investig Suppl. 1968, 97: 31-50.
Beaton LA, Burn TA, Stocki TJ, Chauhan V, Wilkins RC: Development and characterization of an in vitro alpha radiation exposure system. Phys Med Biol. 2011, 56 (12): 3645-3658.
MacPhail SH, Banath JP, Yu TY, Chu EH, Lambur H, Olive PL: Expression of phosphorylated histone H2AX in cultured cell lines following exposure to X-rays. Int J Radiat Biol. 2003, 79 (5): 351-358.
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3 (6): 1101-1108.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (3): Epub 2004 Feb 12
Smyth GK: Bioinformatics and computational biology solutions using R and Bioconductor. Limma: linear models for microarray data. 2005, XIX:397-420.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statist Soc B. 1995, 57 (1): 289-300.
Wang Y, Barbacioru C, Hyland F, Xiao W, Hunkapiller KL, Blake J, Chan F, Gonzalez C, Zhang L, Samaha RR: Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC Genomics. 2006, 7: 59.
Durante M, Manti L: Estimates of radiological risk from a terrorist attack using plutonium. Radiat Environ Biophys. 2002, 41 (2): 125-130.
Lorimore S, Goodhead D, Wright E: Inactivation of haemopoietic stem cells by slow α-particles. Int J Radiat Biol. 1993, 63 (5): 655-660.
Chauhan V, Howland M: Genomic profiling of a human leukemic monocytic cell-line (THP-1) exposed to alpha particle radiation. Sci World J. 2012, 2012: 205038.
Chauhan V, Howland M, Chen J, Kutzner B, Wilkins RC: Differential effects of alpha-particle radiation and X-irradiation on genes associated with apoptosis. Radiol Res Pract. 2011, 2011: 679806.
Chauhan V, Howland M, Wilkins R: A comparative assessment of cytokine expression in human-derived cell lines exposed to alpha particles and X-rays. ScientificWorldJournal. 2012, 2012: 609295.
Paull TT, Rogakou EP, Yamazaki V, Kirchgessner CU, Gellert M, Bonner WM: A critical role for histone H2AX in recruitment of repair factors to nuclear foci after DNA damage. Curr Biol. 2000, 10 (15): 886-895.
Skrzypski M: Quantitative reverse transcriptase real-time polymerase chain reaction (qRT-PCR) in translational oncology: lung cancer perspective. Lung Cancer. 2008, 59 (2): 147-154.
Leatherbarrow EL, Harper JV, Cucinotta FA, O’Neill P: Induction and quantification of γ-H2AX foci following low and high LET-irradiation. Int J Radiat Biol. 2006, 82 (2): 111-118.
Turtoi A, Brown I, Schläger M, Schneeweiss FHA: Gene expression profile of human lymphocytes exposed to 211at α particles. Radiat Res. 2010, 174 (2): 125-136.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/43/prepub
This work was supported by the Canadian Safety and Security Programme. The authors would like to thank Sami Qutob and Lindsay Beaton for editorial comments and Matthew Rodrigues for reference formatting.
The authors declare they have no competing interests.
VC contributed to the conception and design of the study, acquisition of data and analysis and interpretation of data and revising of manuscript. VC drafted the manuscript and provided final approval for publication. MH contributed to the execution of experiments, data analysis and interpretation and revising of manuscript. RW contributed to H2AX data analysis and final manuscript revision. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Chauhan, V., Howland, M. & Wilkins, R. Identification of gene-based responses in human blood cells exposed to alpha particle radiation. BMC Med Genomics 7, 43 (2014). https://doi.org/10.1186/1755-8794-7-43
- Gene expression