BACKGROUND
the goal of radiation biodosimetry is to accurately predict radiation dose, as a surrogate of radiological injury, in a large-scale radiation emergency.
most approaches focus on the first week after an event, when radiological triage is most needed.
a wide range of different methods have been used to identify cellular responses to radiation, ranging from cytogenetic measurements, specifically dicentric assays, to high throughput proteomic, metabolomic, and transcriptomic profiling  <cit> .
several groups have reported gene expression-based analyses, with the development of predictive gene sets that correlate with dose .
radiation quality plays an important role in driving gene expression, and different signaling pathways may be triggered in response to the different types of irradiation.
in addition, various genes and biological functions are induced in cells at different time points after radiation exposure .
to date, most of the gene expression signatures for the purpose of radiation biodosimetry have been obtained in response to photon  exposure.
however, there have been a limited number of studies that have compared gene expression profiles in mice or human blood cells after exposures involving heavy ions and α-particles, the latter representing isotopes likely to be used in a radiological dispersal device  <cit> .
for example, using human peripheral blood mononuclear cells, chauhan et al.
<cit>  found that the gene expression profile induced by α-particle radiation was very similar to the x-ray responses, despite the fact that α-particles are characterized by a higher linear energy transfer coefficient compared with x-rays  <cit> .
in contrast, comparison of α-particle and x-ray irradiation impact on human tumor and endothelial cells  <cit> , as well as human epidermal keratinocytes  <cit>  produced distinct gene expression profiles.
an improvised nuclear device  detonation is capable of producing a mixture of γ-rays and fast neutrons.
furthermore, as neutrons generally have a higher relative biological effectiveness  for most physiological endpoints, it is important to understand the impact that neutrons would have on the biodosimetry methods that are being developed for medical triage purposes.
radiation affects multiple biological processes, including immune system and inflammatory responses, cell cycle progression, cell death, dna repair, as well as metabolism.
however, the impact of neutron radiation on any of these processes has not been studied in detail.
this study presents an initial characterization of gene expression responses following exposure to ind-spectrum neutrons or x-rays.
in the present study, gene expression signatures were analyzed after in vivo irradiation of mice with either neutrons or x-rays.
the doses used were selected based on the limited rbe information available for ind-spectrum neutrons.
using the in vitro cytokinesis-block micronucleus assay  <cit>  in peripheral human blood lymphocytes, xu et al.
<cit>  calculated the rbe of our neutron spectrum to be between  <dig> and  <dig> for micronucleus induction with  <dig> kvp x-rays used as the reference radiation, with similar results for micronucleus induction in mice .
we therefore exposed mice to  <dig>  gy or  <dig> gy of neutron radiation and, separately,  <dig> gy or  <dig> gy of x-ray radiation.
we describe biological processes significantly over-represented only in the neutron response, or only in the x-ray radiation response, as well as processes shared by both types of radiation.
the ultimate goal is to investigate whether we can separately estimate the photon and the neutron component after a mixed photon/neutron exposure, and to develop gene expression signatures capable of discriminating between neutron and photon exposures.
methods
animals and irradiation
all animal experiments were conducted in accordance with applicable federal and state guidelines and were approved by the columbia university institutional animal care and use committee  .
male c57bl/ <dig> mice were received from charles river laboratories  and quarantined for 2 weeks before irradiation at 8-10 weeks of age.
twelve mice were used for each treatment, with six being sacrificed  <dig> day after treatment, and the other six being sacrificed  <dig> days after treatment.
samples were lost from three mice  due to clotting during blood draw.
mice were either sham-irradiated or exposed to  <dig>  gy or  <dig> gy neutron beam or to  <dig> gy or  <dig> gy x-rays.
neutron irradiations were performed using the accelerator-based ind-spectrum neutron source  <cit>  at the radiological research accelerator facility .
briefly, a mixed beam of atomic and molecular ions of hydrogen and deuterium was accelerated to  <dig> mev and used to bombard a thick beryllium target.
the neutrons emitted at 60° to the ion beam axis have a spectrum that closely mimics the hiroshima spectrum at  <dig>  km from the epicenter  <cit> .
for irradiation, six to eighteen mice were placed in holders in adjacent positions on an eighteen position ferris wheel, rotating around the beryllium target at a distance of  <dig>  cm, at an angle of 60° to the particle beam.
the mouse holders, based on  <dig> ml conical tubes, are designed to maintain a horizontal orientation as the wheel rotates, providing an isotropic irradiation, while maintaining the mice in an upright orientation.
with a wheel rotation speed of about  <dig> min per revolution, the dose rate was adjusted so that the smallest dose was delivered in  <dig> rotations .
the mouse holders were reversed end-to-end midway through exposure, to correct for any angular variation in dose rate so that the front and back of the mouse receive equivalent doses.
when fewer than  <dig> mice were on the wheel, two  <dig> ml tubes containing lucite phantoms were placed at either end of the string of occupied mouse holders in order to ensure a uniform scatter dose.
irradiations were performed with a total beam current of  <dig> μa on the target, resulting in a dose rate of  <dig>  gy/h of neutrons and  <dig>  gy/h of γ rays.
in order to allow the mice to acclimate to the restraining conditions of the holders and reduce stress during treatment, the mice were placed in the holders in two separate sessions at  <dig> days and  <dig> day prior to irradiation.
on the second training session one day prior to the actual irradiation, the mice in the holders were placed on the irradiation wheel.
sham-irradiated control mice underwent the same holder acclimation protocol, and on the day of exposure they were placed in the holders on the wheel with the beam off for mock treatment.
controls were performed prior to neutron irradiation, to eliminate dose from activation of materials at the endstation.
for x-ray irradiation, mice were exposed to  <dig> gy or  <dig> gy of x-rays from a westinghouse coronado orthovoltage x-ray machine operating at  <dig> kvp and 15ma with a  <dig> mm cu + 1mm al filter.
the dose rate at the mouse location was  <dig>  gy/min, as determined using a victoreen model  <dig> condenser r meter with a 250r chamber.
after dose administration, mice were housed in micro-isolator cages until the time of sacrifice and blood draw.
no significant changes in cell death or constitution of white blood cells were observed after irradiation.
blood collection and rna isolation
blood was collected by cardiac puncture at the time of euthanasia at days  <dig> and  <dig> post-irradiation.
each sample  was added to a  <dig> ml centrifuge tube that contained  <dig>  ml of paxgene blood rna stabilization and lysis solution  and mixed thoroughly, while a small amount of blood was added to anti-coagulant containing tubes for blood count using a genesis hematology system .
after collection was complete, blood was mixed gently but thoroughly and the tubes were incubated at  <dig> °c for  <dig> h. blood samples were allowed to reach room temperature for  <dig> h before proceeding to rna isolation.
rna was purified following the paxgene rna kit recommendations with on-column dnase i treatment.
globin rna was reduced using the ambion globinclear-mouse/rat kit .
rna yields were quantified using the nanodrop nd <dig> spectrophotometer  and rna quality was checked by the  <dig> bioanalyzer .
high quality rna with an rna integrity number of at least  <dig>  was used for microarray hybridization.
microarray hybridization and data analysis
cyanine- <dig> labeled crna was prepared using the one-color low input quick amp labeling kit .
dye incorporation and crna yield were measured with a nanodrop nd <dig> spectrophotometer .
labeled crna was fragmented and hybridized to agilent mouse gene expression 4x44k v <dig> microarray kit .
slides were scanned with the agilent dna microarray scanner  and the images were analyzed with feature extraction software  using default parameters for background correction and flagging non-uniform features.
background-corrected hybridization intensities were imported into brb-arraytools, version  <dig> .
<dig>  log2-transformed and median normalized.
non-uniform outliers or features not significantly above background intensity in 25% or more of the hybridizations were filtered out.
in addition, a minimum  <dig> -fold change in at least 20% of the hybridizations was set as a requirement.
furthermore, probes were averaged to one probe per gene and duplicate features were reduced by selecting the one with maximum signal intensity.
a final set of  <dig>  features was used in subsequent analyses.
the microarray data is available through the gene expression omnibus with accession number gse <dig>
class comparison was conducted in brb-arraytools to identify genes differentially expressed  between radiation exposed samples and time-matched unirradiated controls using a random variance t-test  <cit> .
time-matched controls were also compared with each other but no significantly differentially expressed genes  were found.
genes with p-values less than  <dig>  were considered statistically significant.
the false discovery rate  was estimated for each gene by the method of benjamini and hochberg  <cit> , to control for false positives.
all genes used in this analysis had a fdr of less than  <dig> .
hierarchical clustering of microarray gene expression data was performed with the dynamic heatmap viewer of the brb-arraytools software using a one minus correlation metric and average linkage.
genes differentially expressed following exposure to  <dig> gy neutron, and  <dig> and  <dig> gy x-ray at either day  <dig> or day  <dig> post-irradiation  were used to construct the heatmap.
gene ontology analysis
lists of genes that were either significantly overexpressed or underexpressed compared with controls were imported separately into the database for annotation, visualization, and integrated discovery , version  <dig>   <cit> , to identify enriched biological processes and gene ontology  terms using the functional annotation tool.
benjamini corrected p values of <  <dig>  were considered significant.
to construct tables  <dig>   <dig>  and  <dig>  redundant go terms were grouped using the revigo software  <cit> .
quantitative rt-pcr
cdna was prepared from total globin-cleared rna using the high-capacity cdna archive kit .
quantitative real-time pcr  was performed for five genes  using taqman assays .
fzr <dig> , ccna <dig> , cdc25b , nusap <dig>  were pre-designed validated assays.
mouse ube2c primers were designed using the primerquest tool , and the sequences were as follows: forward primer, ctgctaggagaacccaacatc; reverse primer: gctggagacctgctttgaata; and probe: ctttgaacacacacgctgcggaac.
a β-actin assay  was also performed alongside as control.
the gene expression validation experiments were conducted with  <dig> ng cdna using universal pcr master mix  in an abi  <dig> real time pcr system.
after an initial activation at  <dig> °c for  <dig> min and  <dig> °c for  <dig> min, the pcr reaction was performed by  <dig> cycles of  <dig> °c for  <dig> s and  <dig> °c for  <dig> s. relative fold-induction was calculated by the 2-ΔΔct method  <cit> , using sds version  <dig>  .
data were normalized to β-actin gene expression levels.
RESULTS
microarray experiments
mice were either sham-irradiated or exposed to  <dig>  gy or  <dig> gy of neutron irradiation or  <dig> gy or  <dig> gy of x-ray irradiation.. all animals remained in apparent good health, with no adverse events noted during the course of the study.
total blood counts were within the normal range in controls and all animals prior to irradiation and at days  <dig> and  <dig> post irradiation.
global gene expression was measured in the blood of mice sacrificed  <dig> and  <dig> days post-irradiation using agilent whole mouse genome microarrays.
class comparison using brb-arraytools identified a total of  <dig> and  <dig> genes differentially expressed  < 5%) between unirradiated controls and neutron or x-ray exposed samples, respectively .
the number of differentially expressed genes varied with time, dose, and radiation quality .
following neutron exposures, genes were predominantly down-regulated, with nearly all differentially expressed genes being under expressed on day  <dig> post-exposure.
up- and down-regulated genes were more evenly divided after x-ray exposure.table  <dig> differentially expressed genes
significantly differentially expressed genes in mouse blood after neutron or x-ray treatment relative to unirradiated controls .
percent of upregulated and downregulated genes are shown.
n025d1:  <dig>  gy neutron, day 1; n025d7:  <dig>  gy neutron, day 7; n1d1:  <dig> gy neutron, day 1; n1d7:  <dig> gy neutron, day 7; x1d1:  <dig> gy x-rays, day 1; x1d7:  <dig> gy x-rays, day 7; x4d1:  <dig> gy x-rays, day 1; x4d7:  <dig> gy x-rays, day 7
we also sought to determine the extent to which neutrons and x-rays share common differentially expressed genes.
venn diagrams revealed the overlap of differentially expressed genes after  <dig> gy neutron or x-ray irradiation to be, for neutron,  <dig> % and  <dig> % at days  <dig> and  <dig>  respectively, and, for x-rays,  <dig> % and  <dig> % at days  <dig> and  <dig>  respectively.
comparing the overlap of differentially expressed genes between  <dig> gy neutron and  <dig> gy x-rays, we found the percentages to be  <dig> %,  <dig> %  and  <dig> %,  <dig> %  at days  <dig> and  <dig>  respectively .
in all,  <dig> genes were common to the response to  <dig> gy neutrons and both  <dig> and  <dig> gy x-rays at day  <dig> after exposure, and  <dig> were common to all three conditions at day  <dig> after exposure.
hierarchical clustering of these most consistently responsive genes was visualized as a heat map to compare the relative expression across samples .
patterns of differential expression were seen to vary by radiation quality and dose, and by time since exposure, showing general consistency between replicates.
control samples showed no obvious differences in expression between days  <dig> and  <dig> fig.
<dig> differentially expressed genes.
venn diagram showing overlap patterns of genes that are differentially expressed at days  <dig> and  <dig> in response to  <dig> gy neutron vs.  <dig> or  <dig> gy x-ray irradiation
fig.
<dig> heat map illustrating relative gene expression at days  <dig> and  <dig> post-irradiation.
hierarchical cluster analysis was performed on  <dig> probesets responding to both  <dig> gy neutrons and  <dig> and  <dig> gy x-rays relative to unexposed controls at either  <dig> or  <dig> days after exposure.
red indicates high expression, green indicates low expression as indicated in the color key.
each row represents one gene and each column represents an individual mouse, ordered by exposure dose and radiation type as labeled in the figure
finally, we searched for genes with sustained differential expression at both day  <dig> and day  <dig>  the higher doses produced the most sustained responses, with  <dig> genes differentially expressed at both times after  <dig> gy neutron and  <dig> after  <dig> gy x-ray exposure.
only  <dig> genes were regulated at both times after  <dig>  gy neutrons, and  <dig> after  <dig> gy x-ray exposure.
thirty-four genes were common to both times after  <dig> gy neutron and  <dig> and  <dig> gy x-ray exposure .
of these genes,  <dig> were uniformly downregulated.
interestingly, the remaining nine genes showed a bi-directional temporal response to irradiation.
these genes were downregulated at day  <dig> post-irradiation, whereas they were upregulated at day  <dig> .
one of the genes temporally regulated was e2f <dig> , a member of the e2f family of transcription factors that play important roles in the control of the cell cycle.
we also performed a search for other differentially expressed e2f factors and found that, unlike e2f <dig>  e2f <dig>  e2f <dig>  and e2f <dig> were downregulated only after  <dig> gy neutron exposure at day  <dig>  showing respectively a  <dig> -,  <dig> -, and  <dig> -fold change from control levels.
in contrast, e2f <dig> and e2f <dig> showed up regulation  after  <dig> gy x-rays at day  <dig> fig.
<dig> genes differentially expressed at both times.
a, venn diagram showing the number of genes differentially expressed at both time points after exposure to  <dig> gy neutrons or  <dig> or  <dig> gy x-rays.
b, list of bi-directionally regulated genes common to all conditions other than  <dig>  gy neutron  with fold change relative to controls for each dose and time.
n1d1:  <dig> gy neutron, day 1; n1d7:  <dig> gy neutron, day 7; x1d1:  <dig> gy x-rays, day 1; x1d7:  <dig> gy x-rays, day 7; x4d1:  <dig> gy x-rays, day 1; x4d7:  <dig> gy x-rays, day 7
gene ontology analysis
differentially expressed genes were functionally classified into gene ontology categories by the database for annotation, visualization, and integrated discovery   <cit>  analysis.
biological processes over-represented among differentially expressed genes with a benjamini-corrected p-value of less than  <dig>  were considered significant .
radiation exposure to the lower dose of neutrons  produced few biological processes, and these were primarily related to immune system response.
therefore, we focused our analysis on the response to the higher dose exposures .
the top biological processes among differentially expressed genes were related to immune system response, which was significantly enriched among genes downregulated in response to both  <dig> gy of x-rays and  <dig> gy neutron irradiation after  <dig> and  <dig> days, and after  <dig> gy x-rays at day  <dig>
a significant class of biological processes among downregulated genes responding to both x-rays  and, particularly, neutrons  was the one related to dna and rna metabolism and processing.
in the x-ray data, biological processes such as dna replication and regulation of transcription were significant only among genes downregulated on day  <dig> after  <dig> gy of x-rays.
the profile of dna/rna-related biological processes significant at day  <dig> after  <dig> gy neutron irradiation  was very similar to that of  <dig> gy x-rays.
however, at day  <dig>  the biological processes enriched at day  <dig> after neutron irradiation were not only retained, but additional dna/rna categories became significant, again only among down-regulated genes.
trna metabolism , rna localization  and transport , and noncoding rna metabolism  were among the new biological processes over-represented at day  <dig> table  <dig> dna and rna metabolism-related biological processes in response to x-ray exposure
rna processing
mrna processing
dna-dependent dna replication
regulation of gene expression, epigenetic
regulation of cell proliferation
significantly enriched biological processes among genes downregulated in response to  <dig> and  <dig> gy x-ray irradiation at days  <dig> and  <dig>  benjamini-corrected p values are shown.
processes unique to x-rays are in bold.
x4d1:  <dig> gy x-rays, day 1; x4d7:  <dig> gy x-rays, day  <dig>  up, upregulated; dn: downregulated
rna transport
nuclear-transcribed mrna catabolic process
trna metabolic process
rna splicing
ncrna metabolic process
rna localization
somatic cell dna recombination
rna catabolic process
gene silencing
regulation of dna metabolic process
significantly enriched biological processes among genes downregulated in response to  <dig> gy neutron irradiation at days  <dig> and  <dig>  benjamini-corrected p values are shown.
processes unique to neutrons are in bold.
n1d1:  <dig> gy neutron, day 1; n1d7:  <dig> gy neutron, day  <dig>  up, upregulated; dn: downregulated
arguably, the most pronounced differences in the transcriptomic profiles of blood cells after neutron and x-rays were in the area of cellular metabolism.
neutron irradiation resulted in the downregulation of genes involved in metabolic processes, including coenzyme, hexose, and lipid biosynthetic processes.
in contrast, exposure of mice to  <dig> gy x-rays resulted in an overrepresentation of metabolic processes among up-regulated genes, especially those involved in lipid, cofactor, and vitamin biological processes .
these processes were detected at day  <dig>  whereas they were no longer significant at day  <dig>  fatty acid metabolism , lipid biosynthesis , and sterol  and steroid  metabolism were among the top enriched biological processes.table  <dig> lipid-related biological processes
significant lipid biosynthesis processes among genes upregulated in response to  <dig> gy x-rays at day  <dig> with benjamini-corrected p values
dna repair pathways were also differentially regulated, showing significant enrichment among down-regulated genes, especially after neutron irradiation.
thus, mismatch repair  appeared downregulated  <dig> day after  <dig> gy neutrons, whereas at day  <dig>  besides mmr , base excision repair , nucleotide excision repair , and non-homologous end joining  were downregulated, as well.
in contrast, no dna repair pathways were significant after  <dig> gy x-rays, whereas mmr  and ber  were downregulated after  <dig> gy x-rays at days  <dig> and  <dig>  respectively.
so far, david analysis had suggested biological processes that appeared to be regulated in blood in response to radiation.
these processes were either significant at only one time point, or if they persisted at both  <dig> and  <dig> days post-irradiation, they were consistently either up- or down-regulated.
however, a careful examination of the biological processes revealed that the cell cycle regulation processes displayed a more complex temporal pattern.
specifically, cell cycle regulation processes were significant among down-regulated genes at day  <dig> after  <dig> gy x-rays exposure.
however the same processes were significant among up-regulated genes at day  <dig>  at day  <dig>   <dig> cell cycle genes were differentially underexpressed, whereas, at day  <dig>   <dig> cell cycle genes were overexpressed.
we searched for genes that appeared in both the downregulated and upregulated columns.
we discovered  <dig> genes, namely mns <dig> , fzr <dig> ), ube2c , ccna <dig> , anln , smc <dig>  and ncapd <dig> .
we examined the expression of these genes in the samples that had been exposed to  <dig> gy x-rays and  <dig> gy neutrons and found  <dig> genes  with the same trend of temporal regulation after  <dig> gy x-rays  and after  <dig> gy neutrons .
performing gene ontology analysis on the genes differentially expressed in samples exposed to  <dig> gy neutrons using a less stringent p value , we discovered that, similar to the x-ray samples, cell cycle biological processes were downregulated at day  <dig> and upregulated at day  <dig>  sixty-six genes with cell cycle gene ontology annotations were underexpressed at day  <dig> and  <dig> genes were upregulated at day  <dig>  besides ube2c, fzr <dig>  and ccna <dig> that were present in both the downregulated and upregulated lists, two more genes, cdc25b , and nusap <dig>  showed a  <dig> -fold  and  <dig> -fold  change at day  <dig>  and  <dig> -fold  and  <dig> -fold  at day  <dig> compared with controls.
when these genes were examined against the x-ray gene profiles, they were found to be bi-directionally temporally regulated, as well.
a heatmap of these five genes  is depicted in fig.
<dig>  the levels of gene expression among control  samples were not significantly different between days  <dig> and  <dig> fig.
<dig> heatmap illustrating relative expression of cell cycle genes bi-directionally regulated after neutron and x-ray irradiation.
the mean  fold change in gene expression relative to controls is color coded according to the scale bar at the bottom of the figure.
measurements were made by microarray analysis
quantitative real-time rt-pcr validation of cell cycle gene expression
we confirmed the expression pattern of the five cell cycle genes, shown by microarray analysis to be first under- then over-expressed after irradiation, by quantitative real-time pcr .
analysis of gene expression of ube2c, fzr <dig>  ccna <dig>  cdc25b, and nusap <dig> by qrt-pcr confirmed that these genes are temporally regulated by irradiation and reverse the direction of their change during the first week after exposure.
the fold-change of these genes in response to x-rays was in good agreement with the fold-change calculated by the microarray experiment, whereas qrt-pcr measurements indicated generally greater fold-changes than the microarrays in response to neutrons.fig.
<dig> gene expression measured by qrt-pcr.
expression of five genes  that were shown by microarray analysis to be bi-directionally regulated is depicted relative to controls and normalized to actb expression.
data represent the mean ± s.e.m.
.
n1d1:  <dig> gy neutron, day 1; n1d7:  <dig> gy neutron, day 7; x1d1:  <dig> gy x-rays, day 1; x1d7:  <dig> gy x-rays, day 7; x4d1:  <dig> gy x-rays, day 1; x4d7:  <dig> gy x-rays, day 7
discussion
the overall goal of this study was to identify differentially regulated genes in response to neutron or x-ray irradiation and perform a comparative analysis of biological processes between the two types of radiation at time points spanning the range of interest for biodosimetry.
our data suggest that the gene transcriptional response varied widely depending on radiation quality, dose, and time since exposure.
it should be noted that these characteristics of gene expression response may contribute to the apparently large number of “unique” genes responding to only one radiation quality.
previous studies have shown differences in the timing of gene expression responses at high and low doses, and following exposure to different radiation qualities.
it is likely that many of the genes seen to respond to only one radiation quality in this study would show a response to the other radiation quality in a different time-dose combination.
some of the observed differences may also be attributable to the different nature of x-rays vs. neutrons.
up to 2/ <dig> of damage from low-let  is due to indirect action, , whereas high-let neutrons cause direct damage to the dna , which is more complex and difficult to repair, and may result in different signaling responses.
a small number of genes showed similar changes after exposure to neutrons and x-rays, and displayed a bidirectional mode of regulation.
a similar temporal pattern of expression for some genes has been described previously in mice injected with 137cscl  <cit> .
in that study, genes were upregulated by day  <dig> or  <dig>  and then downregulated by day  <dig> or  <dig> after isotope administration.
many gene ontology categories, including actin and the cytoskeleton and integrin signaling pathways, showed the same temporal pattern.
in the present study, we identified a number of differentially expressed genes that were significantly different from controls at all radiation types and doses.
a number of these genes were downregulated at day  <dig>  and upregulated at day  <dig>  gene ontology analysis revealed that one biological function, cell cycle, was significant among down-regulated genes at day  <dig> and then significant among up-regulated genes at day  <dig>  all other biological functions were either significant at one dose and time of irradiation, or they were uniformly up- or down-regulated irrespective of time.
most of the genes differentially expressed in both neutron and x-ray exposures were related to immune response, and b and t cell physiology.
these genes were downregulated starting at day  <dig> and reduced expression persisted until the end of the experiments at day  <dig>  widespread decreased expression of immune function genes has been shown previously in both human blood irradiated ex vivo, as well as in vivo mouse peripheral blood following ionizing radiation exposure or 90sr as an internal emitter .
high-dose radiation  has been shown to disrupt immune cell functions, leading to increased cell death of blood cells in mice  <cit> .
moreover, in patients with acute radiation syndrome, hematopoietic cell proliferation is inhibited by radiation exposure  <cit> .
a preponderance of down-regulated genes has previously been associated with higher doses or later times after irradiation , perhaps reflecting greater damage or a failure of repair.
in our study, even  <dig>  gy of neutron radiation resulted in the downregulation of genes involved in immune cell function.
a significant number of genes that were downregulated by both neutron and x-ray exposure at day  <dig> were related to dna and mrna metabolism, gene transcription, rna processing and splicing.
however, at day  <dig>  these processes were no longer overrepresented following x-ray exposure, but they were present after neutron exposure.
moreover, in response to neutron exposure, these processes were enriched with additional related processes, such as trna metabolism/processing, rna transport/localization, and noncoding rna functions.
the contribution of dna and rna related functions in response to dna damage has only recently been appreciated.
a genome-wide sirna screen looking for modulators of dna damage signaling revealed that the largest number of hits were those targeting gene products responsible for nucleic acid metabolism, particularly those involved in mrna binding and processing  <cit> .
furthermore, a phosphoproteomic analysis showed a close link between genome stability and rna synthesis metabolism  <cit> .
likewise, it has been shown that the largest subset of atm/atr/dna-pk substrates identified in a phosphoproteomic screen were proteins linked to rna and dna metabolism, particularly those proteins involved in posttranscriptional mrna regulation  <cit> .
these observations, employing different experimental approaches, highlight the importance of regulatory circuits controlling rna metabolism and stability in dna repair and checkpoint function.
in addition to these findings, our study suggests that neutron  irradiation affects biological processes enriched in trna regulation and rna transportation and localization, as well as non-coding rna metabolism and processing.
trnas have been viewed as passive players involved in protein synthesis.
however, recent evidence suggests that they have more active roles and trna modulation represents a mechanism by which cells achieve altered expression of specific transcripts and proteins.
trna pools in cells can be divided into those that favor proliferation and those that support differentiation  <cit> .
as a result, modifications in trna and their corresponding enzymes are implicated in diseases, including diabetes and cancer.
for example, upregulation of certain trnas increases metastasis in breast cancer patients  <cit> .
furthermore, control of rna transport and localization would be expected to impact on the rate of protein translation  <cit> .
a difference observed between neutron and x-ray response was the enrichment in biological processes involved in lipid biosynthesis and metabolism that was seen only in response to x-ray exposure.
it has long been known that the cellular targets of ionizing radiation, such as x-rays, are not limited to nuclear dna, but that proteins and lipids in other cellular compartments, such as the plasma membrane  <cit> , are also affected.
the action of x-rays has been attributed to the generation of reactive oxygen species that oxidize dna, lipids and proteins  <cit> .
we can speculate that in response to x-rays, cells upregulate lipid, coenzyme, and vitamin biosynthetic and metabolic processes as a means of repairing the damage caused by x-ray irradiation to the cell membrane.
the latter two processes could also serve as anti-oxidant responses  <cit> .
in addition, fatty acid oxidation processes, which are also overrepresented in the x-ray irradiation, would be required by the cells to meet the energy demands of various metabolic processes.
in response to dna damage, cells activate the dna damage and repair signaling pathway.
dna damage that cannot be repaired efficiently leads to cell death or senescence.
although protein abundance and activity do not always follow gene expression changes, our gene ontology analysis suggests an apparent down regulation of several dna repair pathways after neutron but not after x-ray irradiation, especially at day  <dig>  perhaps reflecting a failure to repair the more complex damage resulting from high let radiation.
at day  <dig> after  <dig> gy neutron exposure, go analysis suggested suppression of mmr, whereas at day  <dig>  in addition to mmr, ner, ber, and nhej genes were all significantly over-represented among downregulated genes.
in contrast, exposure to  <dig> gy of x-rays transiently downregulated expression of mmr genes at day  <dig> and ber genes at day  <dig>  homologous recombination, which along with nhej constitutes the major dna double-strand break repair mechanism, was not significantly over-represented among differentially expressed genes.
it has been shown previously that high-let radiation induces complex dna damage that is not easily repaired and nhej is not involved .
more recently, it has been shown that high-let irradiation with protons or carbon ions causes a shift away from nhej toward hrr in the repair of double-strand breaks  <cit> .
consistent with these data, our observed downregulation of genes in the nhej and other dna repair pathways in response to neutron exposure may reflect the fact that these lesions are not repaired by these processes.
regulation at the gene expression level suggests a potential mechanism for favoring the homologous recombination pathway in the attempted repair of neutron damage, and is worthy of further investigation.
a major biological function that is affected by radiation is the cell cycle.
cell cycle regulating genes are important determinants of radiosensitivity and cell fate in response to dna damage.
we studied the effect of neutron radiation on mouse cell cycle-regulated genes and compared it with that of x-rays.
unlike genes in other biological processes, which were either up or downregulated after neutron or x-ray irradiation, many cell cycle genes showed a bidirectional expression based on time.
a group of  <dig> genes, namely ube2c, fzr <dig>  ccna <dig>  cdc25b, and nusap <dig>  were downregulated  <dig> day after irradiation, whereas the same genes were overexpressed  <dig> days post-irradiation.
these temporal changes were further confirmed by qrt-pcr.
a literature search revealed that these genes play important roles in the control of mitosis.
additionally, their protein products are related to the anaphase promoting complex/cyclosome  either as subunits of apc/c  or as substrates .
the apc/c is an e <dig> ubiquitin ligase, which is composed of at least  <dig> core subunits.
the apc/c is active during mitosis and g <dig> phase of the cell cycle.
because of its role in cell cycle regulation, apc/c is important for maintaining genomic integrity  <cit> .
furthermore, apc/c has been implicated in an array of diverse functions ranging from cell differentiation to apoptosis and senescence, as well as cellular metabolism, cell motility, and gene transcription through the degradation of specific substrates  <cit> .
apc/c targets a large repertoire of substrates and recruits them for ubiquitylation via one of two co-activators, cdc <dig> and cdh <dig>   <cit> .
the physiological role of cdh <dig> has been extensively studied in the context of human cancer, since downregulation of cdh <dig> has been reported in many cancers, including those of prostate, ovary, liver, brain, and during the malignant progression of a b-lymphoma cell line.
in mice, fzr <dig> heterozygosity results in the development of epithelial tumors, suggesting that fzr <dig> may be a haploinsufficient tumor suppressor  <cit> .
downregulation of cdh <dig> in post-mitotic neurons has been implicated in neurodegenerative diseases, such as alzheimer’s disease  <cit> .
in addition to its role in mitosis, cdh <dig> has important functions in mediating dna damage response to genotoxic stress  <cit>  that ensure genomic integrity  <cit> .
cdh1-null cells fail to maintain dna damage-induced g <dig> arrest and apc/ccdh <dig> is activated by x-irradiation-induced dna damage .
interestingly, the levels of mitotic cyclins in cdh1-/- cells after dna damage were similar to those of wild-type cells.
these data imply that cyclin a and cyclin b cannot be substrates for apc/ccdh <dig> when it is activated irregularly by dna damage at g <dig>  <cit> .
protein ubiquitination-mediated degradation involves two distinct steps: the covalent attachment of ubiquitin to proteins catalyzed by the sequential actions of the activating , conjugating , and ligating  enzymes, followed by the degradation of the poly-ubiquitylated protein by the 26s proteasome complex.
for apc/c the e <dig> enzymes are ube2c, which is one of the temporally controlled cell cycle genes in this study, and ube2s  <cit> , which is not differentially regulated in response to either neutrons or x-rays.
abundant experimental evidence has shown a role for ube2c in human tumor initiation and progression.
on the other hand, there are very few reports that implicate ube2c in dna damage response to radiation  <cit> .
moreover, the mechanistic details of ube2c response to radiation, as well as their pathophysiological significance remain unexplored.
apc/c regulates spindle formation by promoting the degradation of a number of spindle-binding proteins, including nusap <dig>  <cit> .
the nucleolar spindle-associated protein  <dig>  is a protein highly expressed in proliferating cells and interacts with microtubules  <cit> .
depletion of nusap <dig> caused faulty mitotic spindles, aberrant chromosome segregation, and defective cytokinesis.
overexpression of nusap <dig> caused microtubule bundling and cell cycle arrest at the g2/m checkpoint  <cit> .
cyclin a <dig> interacts with cyclin-dependent kinase  <dig> and controls essential functions in dna replication and cellular proliferation  <cit> .
cyclin a <dig> expression is associated with a poor prognosis in several types of cancer  <cit> .
cyclin a <dig> mrna as well as protein are cell cycle regulated  <cit>  with mrna and protein abundance increasing 4-fold and 20-fold, respectively, as cells progress from g <dig> to g <dig> phase.
apc/c degrades cyclin a <dig> at the end of mitosis, while mrna persists longer than the protein in cells.
to fully appreciate the significance of the temporal differential expression of apc-related genes in the radiation response, we will require knowledge about the status of protein levels and their posttranslational modifications .
however, we can speculate that these changes may be relevant to cell cycle progression, and especially mitosis, thus ensuring genomic stability after irradiation.
another cell cycle regulated gene that appears in the list of temporally bi-directionally expressed genes in the present study is e2f <dig>  while this gene was consistently downregulated at day  <dig> post-irradiation, it always appeared upregulated at day  <dig>  this was in sharp contrast with other e2f genes that showed no change at day  <dig> and downregulation  at day  <dig> post neutron irradiation, or significant upregulation  at day  <dig> post x-ray irradiation.
the e2f family of transcription factors has well known functions in the control of cell cycle, and e2f1- <dig>  especially, in promoting g1/s cell cycle transition and thus cell proliferation  <cit> .
target genes of e2f include several hundred genes that are involved not only in dna replication and cell cycle progression, but also in dna damage repair, apoptosis, differentiation and development  <cit> .
the role of e2fs in mitosis has been shown in cancer cells .
although the specific function of e2f <dig> in response to radiation has not been studied, e2f <dig> transcript and protein levels increase in response to genotoxic stress and maintain genomic stability in neuronal cells  <cit> .
why e2f <dig> shows bi-directional expression changes, whereas other e2f members  do not, is currently not known.
CONCLUSIONS
in the current work we identify genes that are differentially expressed following exposure to neutrons or x-rays, as well as genes that showed similar responses to the two radiation modalities.
in summary, our results show that genes involved in cell cycle regulation are differentially regulated by neutron and x-ray radiation.
however, a few cell cycle genes show a consistent temporal regulation that is common across the two radiation modalities.
these genes are functionally interconnected and play important roles during mitosis.
we have found differing patterns of gene expression response to x-rays and neutrons that vary with both dose and time since exposure.
these findings support the possibility of using gene expression to detect the neutron component of exposures resulting from an ind detonation, thus providing triage information more relevant to the actual radiation injury than an estimate of total dose alone.
further work will be needed to develop gene expression patterns specific to neutron exposure that would be useful for triage following an ind event.
the identification of biomarkers predictive of both dose and type of radiation would be an important advancement in biodosimetry to determine an individual’s exposure and allow more accurate triage for further medical treatment.
the ultimate goal of our neutron studies is to investigate whether we can separately estimate the photon and the neutron component after a mixed photon/neutron exposure, and to develop gene expression signatures capable of discriminating between photon/neutron and pure photon exposures.
following our initial investigation of the gene expression response to neutrons, reported here, our future studies will focus on mixed exposures, including the possibility of synergistic responses.
additional files
additional file 1: excel file with  <dig> tabs containing a summary of all the differentially expressed genes by assay time and radiation type and dose.
additional file 2: excel file with complete output of gene ontology  analyses summarized by dose.
results are in separate tabs for x-ray and neutron exposures.
abbreviations
berbase excision repair
daviddatabase for annotation, visualization, and integrated discovery
fdrfalse discovery rate
gogene ontology
indimprovised nuclear device
mmrmismatch repair
nernucleotide excision repair
nhejnon-homologous end joining
qrt-pcrquantitative real-time polymerase chain reaction
rberelative biological effectiveness

 the goal of radiation biodosimetry is to accurately predict radiation dose, as a surrogate of radiological injury, in a large-scale radiation emergency. most approaches focus on the first week after an event, when radiological triage is most needed. a wide range of different methods have been used to identify cellular responses to radiation, ranging from cytogenetic measurements, specifically dicentric assays, to high throughput proteomic, metabolomic, and transcriptomic profiling. several groups have reported gene expression-based analyses, with the development of predictive gene sets that correlate with dose. radiation quality plays an important role in driving gene expression, and different signaling pathways may be triggered in response to the different types of irradiation. in addition, various genes and biological functions are induced in cells at different time points after radiation exposure. to date, most of the gene expression signatures for the purpose of radiation biodosimetry have been obtained in response to photon exposure. however, there have been a limited number of studies that have compared gene expression profiles in mice or human blood cells after exposures involving heavy ions and α-particles, the latter representing isotopes likely to be used in a radiological dispersal device. for example, using human peripheral blood mononuclear cells, chauhan et al. found that the gene expression profile induced by α-particle radiation was very similar to the x-ray responses, despite the fact that α-particles are characterized by a higher linear energy transfer coefficient compared with x-rays. in contrast, comparison of α-particle and x-ray irradiation impact on human tumor and endothelial cells, as well as human epidermal keratinocytes produced distinct gene expression profiles. an improvised nuclear device detonation is capable of producing a mixture of γ-rays and fast neutrons. furthermore, as neutrons generally have a higher relative biological effectiveness for most physiological endpoints, it is important to understand the impact that neutrons would have on the biodosimetry methods that are being developed for medical triage purposes. radiation affects multiple biological processes, including immune system and inflammatory responses, cell cycle progression, cell death, dna repair, as well as metabolism. however, the impact of neutron radiation on any of these processes has not been studied in detail. this study presents an initial characterization of gene expression responses following exposure to ind-spectrum neutrons or x-rays. in the present study, gene expression signatures were analyzed after in vivo irradiation of mice with either neutrons or x-rays. the doses used were selected based on the limited rbe information available for ind-spectrum neutrons. using the in vitro cytokinesis-block micronucleus assay in peripheral human blood lymphocytes, xu et al. calculated the rbe of our neutron spectrum to be between and for micronucleus induction with kvp x-rays used as the reference radiation, with similar results for micronucleus induction in mice. we therefore exposed mice to gy or gy of neutron radiation and, separately, gy or gy of x-ray radiation. we describe biological processes significantly over-represented only in the neutron response, or only in the x-ray radiation response, as well as processes shared by both types of radiation. the ultimate goal is to investigate whether we can separately estimate the photon and the neutron component after a mixed photon/neutron exposure, and to develop gene expression signatures capable of discriminating between neutron and photon exposures. methods animals and irradiation all animal experiments were conducted in accordance with applicable federal and state guidelines and were approved by the columbia university institutional animal care and use committee. male c57bl/ mice were received from charles river laboratories and quarantined for 2 weeks before irradiation at 8-10 weeks of age. twelve mice were used for each treatment, with six being sacrificed day after treatment, and the other six being sacrificed days after treatment. samples were lost from three mice due to clotting during blood draw. mice were either sham-irradiated or exposed to gy or gy neutron beam or to gy or gy x-rays. neutron irradiations were performed using the accelerator-based ind-spectrum neutron source at the radiological research accelerator facility. briefly, a mixed beam of atomic and molecular ions of hydrogen and deuterium was accelerated to mev and used to bombard a thick beryllium target. the neutrons emitted at 60° to the ion beam axis have a spectrum that closely mimics the hiroshima spectrum at km from the epicenter. for irradiation, six to eighteen mice were placed in holders in adjacent positions on an eighteen position ferris wheel, rotating around the beryllium target at a distance of cm, at an angle of 60° to the particle beam. the mouse holders, based on ml conical tubes, are designed to maintain a horizontal orientation as the wheel rotates, providing an isotropic irradiation, while maintaining the mice in an upright orientation. with a wheel rotation speed of about min per revolution, the dose rate was adjusted so that the smallest dose was delivered in rotations. the mouse holders were reversed end-to-end midway through exposure, to correct for any angular variation in dose rate so that the front and back of the mouse receive equivalent doses. when fewer than mice were on the wheel, two ml tubes containing lucite phantoms were placed at either end of the string of occupied mouse holders in order to ensure a uniform scatter dose. irradiations were performed with a total beam current of μa on the target, resulting in a dose rate of gy/h of neutrons and gy/h of γ rays. in order to allow the mice to acclimate to the restraining conditions of the holders and reduce stress during treatment, the mice were placed in the holders in two separate sessions at days and day prior to irradiation. on the second training session one day prior to the actual irradiation, the mice in the holders were placed on the irradiation wheel. sham-irradiated control mice underwent the same holder acclimation protocol, and on the day of exposure they were placed in the holders on the wheel with the beam off for mock treatment. controls were performed prior to neutron irradiation, to eliminate dose from activation of materials at the endstation. for x-ray irradiation, mice were exposed to gy or gy of x-rays from a westinghouse coronado orthovoltage x-ray machine operating at kvp and 15ma with a mm cu + 1mm al filter. the dose rate at the mouse location was gy/min, as determined using a victoreen model condenser r meter with a 250r chamber. after dose administration, mice were housed in micro-isolator cages until the time of sacrifice and blood draw. no significant changes in cell death or constitution of white blood cells were observed after irradiation. blood collection and rna isolation blood was collected by cardiac puncture at the time of euthanasia at days and post-irradiation. each sample was added to a ml centrifuge tube that contained ml of paxgene blood rna stabilization and lysis solution and mixed thoroughly, while a small amount of blood was added to anti-coagulant containing tubes for blood count using a genesis hematology system. after collection was complete, blood was mixed gently but thoroughly and the tubes were incubated at °c for h. blood samples were allowed to reach room temperature for h before proceeding to rna isolation. rna was purified following the paxgene rna kit recommendations with on-column dnase i treatment. globin rna was reduced using the ambion globinclear-mouse/rat kit. rna yields were quantified using the nanodrop nd spectrophotometer and rna quality was checked by the bioanalyzer. high quality rna with an rna integrity number of at least was used for microarray hybridization. microarray hybridization and data analysis cyanine- labeled crna was prepared using the one-color low input quick amp labeling kit. dye incorporation and crna yield were measured with a nanodrop nd spectrophotometer. labeled crna was fragmented and hybridized to agilent mouse gene expression 4x44k v microarray kit. slides were scanned with the agilent dna microarray scanner and the images were analyzed with feature extraction software using default parameters for background correction and flagging non-uniform features. background-corrected hybridization intensities were imported into brb-arraytools, version. log2-transformed and median normalized. non-uniform outliers or features not significantly above background intensity in 25% or more of the hybridizations were filtered out. in addition, a minimum -fold change in at least 20% of the hybridizations was set as a requirement. furthermore, probes were averaged to one probe per gene and duplicate features were reduced by selecting the one with maximum signal intensity. a final set of features was used in subsequent analyses. the microarray data is available through the gene expression omnibus with accession number gse class comparison was conducted in brb-arraytools to identify genes differentially expressed between radiation exposed samples and time-matched unirradiated controls using a random variance t-test. time-matched controls were also compared with each other but no significantly differentially expressed genes were found. genes with p-values less than were considered statistically significant. the false discovery rate was estimated for each gene by the method of benjamini and hochberg, to control for false positives. all genes used in this analysis had a fdr of less than. hierarchical clustering of microarray gene expression data was performed with the dynamic heatmap viewer of the brb-arraytools software using a one minus correlation metric and average linkage. genes differentially expressed following exposure to gy neutron, and and gy x-ray at either day or day post-irradiation were used to construct the heatmap. gene ontology analysis lists of genes that were either significantly overexpressed or underexpressed compared with controls were imported separately into the database for annotation, visualization, and integrated discovery, version, to identify enriched biological processes and gene ontology terms using the functional annotation tool. benjamini corrected p values of <  were considered significant. to construct tables  and redundant go terms were grouped using the revigo software. quantitative rt-pcr cdna was prepared from total globin-cleared rna using the high-capacity cdna archive kit. quantitative real-time pcr was performed for five genes using taqman assays. fzr, ccna, cdc25b, nusap were pre-designed validated assays. mouse ube2c primers were designed using the primerquest tool, and the sequences were as follows: forward primer, ctgctaggagaacccaacatc; reverse primer: gctggagacctgctttgaata; and probe: ctttgaacacacacgctgcggaac. a β-actin assay was also performed alongside as control. the gene expression validation experiments were conducted with ng cdna using universal pcr master mix in an abi real time pcr system. after an initial activation at °c for min and °c for min, the pcr reaction was performed by cycles of °c for s and °c for s. relative fold-induction was calculated by the 2-ΔΔct method, using sds version. data were normalized to β-actin gene expression levels.  microarray experiments mice were either sham-irradiated or exposed to gy or gy of neutron irradiation or gy or gy of x-ray irradiation.. all animals remained in apparent good health, with no adverse events noted during the course of the study. total blood counts were within the normal range in controls and all animals prior to irradiation and at days and post irradiation. global gene expression was measured in the blood of mice sacrificed and days post-irradiation using agilent whole mouse genome microarrays. class comparison using brb-arraytools identified a total of and genes differentially expressed  < 5%) between unirradiated controls and neutron or x-ray exposed samples, respectively. the number of differentially expressed genes varied with time, dose, and radiation quality. following neutron exposures, genes were predominantly down-regulated, with nearly all differentially expressed genes being under expressed on day post-exposure. up- and down-regulated genes were more evenly divided after x-ray exposure.table differentially expressed genes significantly differentially expressed genes in mouse blood after neutron or x-ray treatment relative to unirradiated controls. percent of upregulated and downregulated genes are shown. n025d1: gy neutron, day 1; n025d7: gy neutron, day 7; n1d1: gy neutron, day 1; n1d7: gy neutron, day 7; x1d1: gy x-rays, day 1; x1d7: gy x-rays, day 7; x4d1: gy x-rays, day 1; x4d7: gy x-rays, day 7 we also sought to determine the extent to which neutrons and x-rays share common differentially expressed genes. venn diagrams revealed the overlap of differentially expressed genes after gy neutron or x-ray irradiation to be, for neutron, % and % at days and respectively, and, for x-rays, % and % at days and respectively. comparing the overlap of differentially expressed genes between gy neutron and gy x-rays, we found the percentages to be %, % and %, % at days and respectively. in all, genes were common to the response to gy neutrons and both and gy x-rays at day after exposure, and were common to all three conditions at day after exposure. hierarchical clustering of these most consistently responsive genes was visualized as a heat map to compare the relative expression across samples. patterns of differential expression were seen to vary by radiation quality and dose, and by time since exposure, showing general consistency between replicates. control samples showed no obvious differences in expression between days and fig. differentially expressed genes. venn diagram showing overlap patterns of genes that are differentially expressed at days and in response to gy neutron vs. or gy x-ray irradiation fig. heat map illustrating relative gene expression at days and post-irradiation. hierarchical cluster analysis was performed on probesets responding to both gy neutrons and and gy x-rays relative to unexposed controls at either or days after exposure. red indicates high expression, green indicates low expression as indicated in the color key. each row represents one gene and each column represents an individual mouse, ordered by exposure dose and radiation type as labeled in the figure finally, we searched for genes with sustained differential expression at both day and day the higher doses produced the most sustained responses, with genes differentially expressed at both times after gy neutron and after gy x-ray exposure. only genes were regulated at both times after gy neutrons, and after gy x-ray exposure. thirty-four genes were common to both times after gy neutron and and gy x-ray exposure. of these genes, were uniformly downregulated. interestingly, the remaining nine genes showed a bi-directional temporal response to irradiation. these genes were downregulated at day post-irradiation, whereas they were upregulated at day. one of the genes temporally regulated was e2f, a member of the e2f family of transcription factors that play important roles in the control of the cell cycle. we also performed a search for other differentially expressed e2f factors and found that, unlike e2f e2f e2f and e2f were downregulated only after gy neutron exposure at day showing respectively a -, -, and -fold change from control levels. in contrast, e2f and e2f showed up regulation after gy x-rays at day fig. genes differentially expressed at both times. a, venn diagram showing the number of genes differentially expressed at both time points after exposure to gy neutrons or or gy x-rays. b, list of bi-directionally regulated genes common to all conditions other than gy neutron with fold change relative to controls for each dose and time. n1d1: gy neutron, day 1; n1d7: gy neutron, day 7; x1d1: gy x-rays, day 1; x1d7: gy x-rays, day 7; x4d1: gy x-rays, day 1; x4d7: gy x-rays, day 7 gene ontology analysis differentially expressed genes were functionally classified into gene ontology categories by the database for annotation, visualization, and integrated discovery analysis. biological processes over-represented among differentially expressed genes with a benjamini-corrected p-value of less than were considered significant. radiation exposure to the lower dose of neutrons produced few biological processes, and these were primarily related to immune system response. therefore, we focused our analysis on the response to the higher dose exposures. the top biological processes among differentially expressed genes were related to immune system response, which was significantly enriched among genes downregulated in response to both gy of x-rays and gy neutron irradiation after and days, and after gy x-rays at day a significant class of biological processes among downregulated genes responding to both x-rays and, particularly, neutrons was the one related to dna and rna metabolism and processing. in the x-ray data, biological processes such as dna replication and regulation of transcription were significant only among genes downregulated on day after gy of x-rays. the profile of dna/rna-related biological processes significant at day after gy neutron irradiation was very similar to that of gy x-rays. however, at day the biological processes enriched at day after neutron irradiation were not only retained, but additional dna/rna categories became significant, again only among down-regulated genes. trna metabolism, rna localization and transport, and noncoding rna metabolism were among the new biological processes over-represented at day table dna and rna metabolism-related biological processes in response to x-ray exposure rna processing mrna processing dna-dependent dna replication regulation of gene expression, epigenetic regulation of cell proliferation significantly enriched biological processes among genes downregulated in response to and gy x-ray irradiation at days and benjamini-corrected p values are shown. processes unique to x-rays are in bold. x4d1: gy x-rays, day 1; x4d7: gy x-rays, day up, upregulated; dn: downregulated rna transport nuclear-transcribed mrna catabolic process trna metabolic process rna splicing ncrna metabolic process rna localization somatic cell dna recombination rna catabolic process gene silencing regulation of dna metabolic process significantly enriched biological processes among genes downregulated in response to gy neutron irradiation at days and benjamini-corrected p values are shown. processes unique to neutrons are in bold. n1d1: gy neutron, day 1; n1d7: gy neutron, day up, upregulated; dn: downregulated arguably, the most pronounced differences in the transcriptomic profiles of blood cells after neutron and x-rays were in the area of cellular metabolism. neutron irradiation resulted in the downregulation of genes involved in metabolic processes, including coenzyme, hexose, and lipid biosynthetic processes. in contrast, exposure of mice to gy x-rays resulted in an overrepresentation of metabolic processes among up-regulated genes, especially those involved in lipid, cofactor, and vitamin biological processes. these processes were detected at day whereas they were no longer significant at day fatty acid metabolism, lipid biosynthesis, and sterol and steroid metabolism were among the top enriched biological processes.table lipid-related biological processes significant lipid biosynthesis processes among genes upregulated in response to gy x-rays at day with benjamini-corrected p values dna repair pathways were also differentially regulated, showing significant enrichment among down-regulated genes, especially after neutron irradiation. thus, mismatch repair appeared downregulated day after gy neutrons, whereas at day besides mmr, base excision repair, nucleotide excision repair, and non-homologous end joining were downregulated, as well. in contrast, no dna repair pathways were significant after gy x-rays, whereas mmr and ber were downregulated after gy x-rays at days and respectively. so far, david analysis had suggested biological processes that appeared to be regulated in blood in response to radiation. these processes were either significant at only one time point, or if they persisted at both and days post-irradiation, they were consistently either up- or down-regulated. however, a careful examination of the biological processes revealed that the cell cycle regulation processes displayed a more complex temporal pattern. specifically, cell cycle regulation processes were significant among down-regulated genes at day after gy x-rays exposure. however the same processes were significant among up-regulated genes at day at day cell cycle genes were differentially underexpressed, whereas, at day cell cycle genes were overexpressed. we searched for genes that appeared in both the downregulated and upregulated columns. we discovered genes, namely mns, fzr ), ube2c, ccna, anln, smc and ncapd. we examined the expression of these genes in the samples that had been exposed to gy x-rays and gy neutrons and found genes with the same trend of temporal regulation after gy x-rays and after gy neutrons. performing gene ontology analysis on the genes differentially expressed in samples exposed to gy neutrons using a less stringent p value, we discovered that, similar to the x-ray samples, cell cycle biological processes were downregulated at day and upregulated at day sixty-six genes with cell cycle gene ontology annotations were underexpressed at day and genes were upregulated at day besides ube2c, fzr and ccna that were present in both the downregulated and upregulated lists, two more genes, cdc25b, and nusap showed a -fold and -fold change at day and -fold and -fold at day compared with controls. when these genes were examined against the x-ray gene profiles, they were found to be bi-directionally temporally regulated, as well. a heatmap of these five genes is depicted in fig.  the levels of gene expression among control samples were not significantly different between days and fig. heatmap illustrating relative expression of cell cycle genes bi-directionally regulated after neutron and x-ray irradiation. the mean fold change in gene expression relative to controls is color coded according to the scale bar at the bottom of the figure. measurements were made by microarray analysis quantitative real-time rt-pcr validation of cell cycle gene expression we confirmed the expression pattern of the five cell cycle genes, shown by microarray analysis to be first under- then over-expressed after irradiation, by quantitative real-time pcr. analysis of gene expression of ube2c, fzr ccna cdc25b, and nusap by qrt-pcr confirmed that these genes are temporally regulated by irradiation and reverse the direction of their change during the first week after exposure. the fold-change of these genes in response to x-rays was in good agreement with the fold-change calculated by the microarray experiment, whereas qrt-pcr measurements indicated generally greater fold-changes than the microarrays in response to neutrons.fig. gene expression measured by qrt-pcr. expression of five genes that were shown by microarray analysis to be bi-directionally regulated is depicted relative to controls and normalized to actb expression. data represent the mean ± s.e.m.. n1d1: gy neutron, day 1; n1d7: gy neutron, day 7; x1d1: gy x-rays, day 1; x1d7: gy x-rays, day 7; x4d1: gy x-rays, day 1; x4d7: gy x-rays, day 7 discussion the overall goal of this study was to identify differentially regulated genes in response to neutron or x-ray irradiation and perform a comparative analysis of biological processes between the two types of radiation at time points spanning the range of interest for biodosimetry. our data suggest that the gene transcriptional response varied widely depending on radiation quality, dose, and time since exposure. it should be noted that these characteristics of gene expression response may contribute to the apparently large number of “unique” genes responding to only one radiation quality. previous studies have shown differences in the timing of gene expression responses at high and low doses, and following exposure to different radiation qualities. it is likely that many of the genes seen to respond to only one radiation quality in this study would show a response to the other radiation quality in a different time-dose combination. some of the observed differences may also be attributable to the different nature of x-rays vs. neutrons. up to 2/ of damage from low-let is due to indirect action,, whereas high-let neutrons cause direct damage to the dna, which is more complex and difficult to repair, and may result in different signaling responses. a small number of genes showed similar changes after exposure to neutrons and x-rays, and displayed a bidirectional mode of regulation. a similar temporal pattern of expression for some genes has been described previously in mice injected with 137cscl. in that study, genes were upregulated by day or and then downregulated by day or after isotope administration. many gene ontology categories, including actin and the cytoskeleton and integrin signaling pathways, showed the same temporal pattern. in the present study, we identified a number of differentially expressed genes that were significantly different from controls at all radiation types and doses. a number of these genes were downregulated at day and upregulated at day gene ontology analysis revealed that one biological function, cell cycle, was significant among down-regulated genes at day and then significant among up-regulated genes at day all other biological functions were either significant at one dose and time of irradiation, or they were uniformly up- or down-regulated irrespective of time. most of the genes differentially expressed in both neutron and x-ray exposures were related to immune response, and b and t cell physiology. these genes were downregulated starting at day and reduced expression persisted until the end of the experiments at day widespread decreased expression of immune function genes has been shown previously in both human blood irradiated ex vivo, as well as in vivo mouse peripheral blood following ionizing radiation exposure or 90sr as an internal emitter. high-dose radiation has been shown to disrupt immune cell functions, leading to increased cell death of blood cells in mice. moreover, in patients with acute radiation syndrome, hematopoietic cell proliferation is inhibited by radiation exposure. a preponderance of down-regulated genes has previously been associated with higher doses or later times after irradiation, perhaps reflecting greater damage or a failure of repair. in our study, even gy of neutron radiation resulted in the downregulation of genes involved in immune cell function. a significant number of genes that were downregulated by both neutron and x-ray exposure at day were related to dna and mrna metabolism, gene transcription, rna processing and splicing. however, at day these processes were no longer overrepresented following x-ray exposure, but they were present after neutron exposure. moreover, in response to neutron exposure, these processes were enriched with additional related processes, such as trna metabolism/processing, rna transport/localization, and noncoding rna functions. the contribution of dna and rna related functions in response to dna damage has only recently been appreciated. a genome-wide sirna screen looking for modulators of dna damage signaling revealed that the largest number of hits were those targeting gene products responsible for nucleic acid metabolism, particularly those involved in mrna binding and processing. furthermore, a phosphoproteomic analysis showed a close link between genome stability and rna synthesis metabolism. likewise, it has been shown that the largest subset of atm/atr/dna-pk substrates identified in a phosphoproteomic screen were proteins linked to rna and dna metabolism, particularly those proteins involved in posttranscriptional mrna regulation. these observations, employing different experimental approaches, highlight the importance of regulatory circuits controlling rna metabolism and stability in dna repair and checkpoint function. in addition to these findings, our study suggests that neutron irradiation affects biological processes enriched in trna regulation and rna transportation and localization, as well as non-coding rna metabolism and processing. trnas have been viewed as passive players involved in protein synthesis. however, recent evidence suggests that they have more active roles and trna modulation represents a mechanism by which cells achieve altered expression of specific transcripts and proteins. trna pools in cells can be divided into those that favor proliferation and those that support differentiation. as a result, modifications in trna and their corresponding enzymes are implicated in diseases, including diabetes and cancer. for example, upregulation of certain trnas increases metastasis in breast cancer patients. furthermore, control of rna transport and localization would be expected to impact on the rate of protein translation. a difference observed between neutron and x-ray response was the enrichment in biological processes involved in lipid biosynthesis and metabolism that was seen only in response to x-ray exposure. it has long been known that the cellular targets of ionizing radiation, such as x-rays, are not limited to nuclear dna, but that proteins and lipids in other cellular compartments, such as the plasma membrane, are also affected. the action of x-rays has been attributed to the generation of reactive oxygen species that oxidize dna, lipids and proteins. we can speculate that in response to x-rays, cells upregulate lipid, coenzyme, and vitamin biosynthetic and metabolic processes as a means of repairing the damage caused by x-ray irradiation to the cell membrane. the latter two processes could also serve as anti-oxidant responses. in addition, fatty acid oxidation processes, which are also overrepresented in the x-ray irradiation, would be required by the cells to meet the energy demands of various metabolic processes. in response to dna damage, cells activate the dna damage and repair signaling pathway. dna damage that cannot be repaired efficiently leads to cell death or senescence. although protein abundance and activity do not always follow gene expression changes, our gene ontology analysis suggests an apparent down regulation of several dna repair pathways after neutron but not after x-ray irradiation, especially at day perhaps reflecting a failure to repair the more complex damage resulting from high let radiation. at day after gy neutron exposure, go analysis suggested suppression of mmr, whereas at day in addition to mmr, ner, ber, and nhej genes were all significantly over-represented among downregulated genes. in contrast, exposure to gy of x-rays transiently downregulated expression of mmr genes at day and ber genes at day homologous recombination, which along with nhej constitutes the major dna double-strand break repair mechanism, was not significantly over-represented among differentially expressed genes. it has been shown previously that high-let radiation induces complex dna damage that is not easily repaired and nhej is not involved. more recently, it has been shown that high-let irradiation with protons or carbon ions causes a shift away from nhej toward hrr in the repair of double-strand breaks. consistent with these data, our observed downregulation of genes in the nhej and other dna repair pathways in response to neutron exposure may reflect the fact that these lesions are not repaired by these processes. regulation at the gene expression level suggests a potential mechanism for favoring the homologous recombination pathway in the attempted repair of neutron damage, and is worthy of further investigation. a major biological function that is affected by radiation is the cell cycle. cell cycle regulating genes are important determinants of radiosensitivity and cell fate in response to dna damage. we studied the effect of neutron radiation on mouse cell cycle-regulated genes and compared it with that of x-rays. unlike genes in other biological processes, which were either up or downregulated after neutron or x-ray irradiation, many cell cycle genes showed a bidirectional expression based on time. a group of genes, namely ube2c, fzr ccna cdc25b, and nusap were downregulated day after irradiation, whereas the same genes were overexpressed days post-irradiation. these temporal changes were further confirmed by qrt-pcr. a literature search revealed that these genes play important roles in the control of mitosis. additionally, their protein products are related to the anaphase promoting complex/cyclosome either as subunits of apc/c or as substrates. the apc/c is an e ubiquitin ligase, which is composed of at least core subunits. the apc/c is active during mitosis and g phase of the cell cycle. because of its role in cell cycle regulation, apc/c is important for maintaining genomic integrity. furthermore, apc/c has been implicated in an array of diverse functions ranging from cell differentiation to apoptosis and senescence, as well as cellular metabolism, cell motility, and gene transcription through the degradation of specific substrates. apc/c targets a large repertoire of substrates and recruits them for ubiquitylation via one of two co-activators, cdc and cdh. the physiological role of cdh has been extensively studied in the context of human cancer, since downregulation of cdh has been reported in many cancers, including those of prostate, ovary, liver, brain, and during the malignant progression of a b-lymphoma cell line. in mice, fzr heterozygosity results in the development of epithelial tumors, suggesting that fzr may be a haploinsufficient tumor suppressor. downregulation of cdh in post-mitotic neurons has been implicated in neurodegenerative diseases, such as alzheimer’s disease. in addition to its role in mitosis, cdh has important functions in mediating dna damage response to genotoxic stress that ensure genomic integrity. cdh1-null cells fail to maintain dna damage-induced g arrest and apc/ccdh is activated by x-irradiation-induced dna damage. interestingly, the levels of mitotic cyclins in cdh1-/- cells after dna damage were similar to those of wild-type cells. these data imply that cyclin a and cyclin b cannot be substrates for apc/ccdh when it is activated irregularly by dna damage at g. protein ubiquitination-mediated degradation involves two distinct steps: the covalent attachment of ubiquitin to proteins catalyzed by the sequential actions of the activating, conjugating, and ligating enzymes, followed by the degradation of the poly-ubiquitylated protein by the 26s proteasome complex. for apc/c the e enzymes are ube2c, which is one of the temporally controlled cell cycle genes in this study, and ube2s, which is not differentially regulated in response to either neutrons or x-rays. abundant experimental evidence has shown a role for ube2c in human tumor initiation and progression. on the other hand, there are very few reports that implicate ube2c in dna damage response to radiation. moreover, the mechanistic details of ube2c response to radiation, as well as their pathophysiological significance remain unexplored. apc/c regulates spindle formation by promoting the degradation of a number of spindle-binding proteins, including nusap. the nucleolar spindle-associated protein is a protein highly expressed in proliferating cells and interacts with microtubules. depletion of nusap caused faulty mitotic spindles, aberrant chromosome segregation, and defective cytokinesis. overexpression of nusap caused microtubule bundling and cell cycle arrest at the g2/m checkpoint. cyclin a interacts with cyclin-dependent kinase and controls essential functions in dna replication and cellular proliferation. cyclin a expression is associated with a poor prognosis in several types of cancer. cyclin a mrna as well as protein are cell cycle regulated with mrna and protein abundance increasing 4-fold and 20-fold, respectively, as cells progress from g to g phase. apc/c degrades cyclin a at the end of mitosis, while mrna persists longer than the protein in cells. to fully appreciate the significance of the temporal differential expression of apc-related genes in the radiation response, we will require knowledge about the status of protein levels and their posttranslational modifications. however, we can speculate that these changes may be relevant to cell cycle progression, and especially mitosis, thus ensuring genomic stability after irradiation. another cell cycle regulated gene that appears in the list of temporally bi-directionally expressed genes in the present study is e2f while this gene was consistently downregulated at day post-irradiation, it always appeared upregulated at day this was in sharp contrast with other e2f genes that showed no change at day and downregulation at day post neutron irradiation, or significant upregulation at day post x-ray irradiation. the e2f family of transcription factors has well known functions in the control of cell cycle, and e2f1- especially, in promoting g1/s cell cycle transition and thus cell proliferation. target genes of e2f include several hundred genes that are involved not only in dna replication and cell cycle progression, but also in dna damage repair, apoptosis, differentiation and development. the role of e2fs in mitosis has been shown in cancer cells. although the specific function of e2f in response to radiation has not been studied, e2f transcript and protein levels increase in response to genotoxic stress and maintain genomic stability in neuronal cells. why e2f shows bi-directional expression changes, whereas other e2f members do not, is currently not known.  in the current work we identify genes that are differentially expressed following exposure to neutrons or x-rays, as well as genes that showed similar responses to the two radiation modalities. in summary, our results show that genes involved in cell cycle regulation are differentially regulated by neutron and x-ray radiation. however, a few cell cycle genes show a consistent temporal regulation that is common across the two radiation modalities. these genes are functionally interconnected and play important roles during mitosis. we have found differing patterns of gene expression response to x-rays and neutrons that vary with both dose and time since exposure. these findings support the possibility of using gene expression to detect the neutron component of exposures resulting from an ind detonation, thus providing triage information more relevant to the actual radiation injury than an estimate of total dose alone. further work will be needed to develop gene expression patterns specific to neutron exposure that would be useful for triage following an ind event. the identification of biomarkers predictive of both dose and type of radiation would be an important advancement in biodosimetry to determine an individual’s exposure and allow more accurate triage for further medical treatment. the ultimate goal of our neutron studies is to investigate whether we can separately estimate the photon and the neutron component after a mixed photon/neutron exposure, and to develop gene expression signatures capable of discriminating between photon/neutron and pure photon exposures. following our initial investigation of the gene expression response to neutrons, reported here, our future studies will focus on mixed exposures, including the possibility of synergistic responses. additional files additional file 1: excel file with tabs containing a summary of all the differentially expressed genes by assay time and radiation type and dose. additional file 2: excel file with complete output of gene ontology analyses summarized by dose. results are in separate tabs for x-ray and neutron exposures. abbreviations berbase excision repair daviddatabase for annotation, visualization, and integrated discovery fdrfalse discovery rate gogene ontology indimprovised nuclear device mmrmismatch repair nernucleotide excision repair nhejnon-homologous end joining qrt-pcrquantitative real-time polymerase chain reaction rberelative biological effectiveness 
neutron exposure resulted in mostly downregulated genes, whereas x-rays showed both down- and up-regulated genes. most biological processes significantly enriched at both timepoints were consistently represented by either under- or over-expressed genes. we describe, for the first time, the gene expression profile of mouse blood cells following exposure to neutrons. we have found that neutron radiation results in both distinct and common gene expression patterns compared with x-ray radiation. mice exposed to either or gy of neutron or or gy x-ray radiation were sacrificed at or days after exposure. whole genome microarray analysis identified and differentially expressed genes in the blood of mice exposed to neutron or x-ray radiation, respectively. cell cycle genes downregulated at day were mostly distinct from the cell cycle genes upregulated at day however, five cell cycle genes, fzr ube2c, ccna nusap and cdc25b, were both downregulated at day and upregulated at day in the event of an improvised nuclear device detonation, the prompt radiation exposure would consist of photons plus a neutron component that would contribute to the total dose.
 in the event of an improvised nuclear device detonation, the prompt radiation exposure would consist of photons plus a neutron component that would contribute to the total dose. as neutrons cause more complex and difficult to repair damage to cells that would result in a more severe health burden to affected individuals, it is paramount to be able to estimate the contribution of neutrons to an estimated dose, to provide information for those making treatment decisions.  mice exposed to either or gy of neutron or or gy x-ray radiation were sacrificed at or days after exposure. whole genome microarray analysis identified and differentially expressed genes in the blood of mice exposed to neutron or x-ray radiation, respectively. neutron exposure resulted in mostly downregulated genes, whereas x-rays showed both down- and up-regulated genes. a total of differentially expressed genes were regulated in response to all ≥ gy exposures at both times. of these, genes were consistently downregulated at days and whereas genes, including the transcription factor e2f showed bi-directional regulation; being downregulated at day while upregulated at day gene ontology analysis revealed that genes involved in nucleic acid metabolism processes were persistently downregulated in neutron irradiated mice, whereas genes involved in lipid metabolism were upregulated in x-ray irradiated animals. most biological processes significantly enriched at both timepoints were consistently represented by either under- or over-expressed genes. in contrast, cell cycle processes were significant among down-regulated genes at day but among up-regulated genes at day after exposure to either neutron or x-rays. cell cycle genes downregulated at day were mostly distinct from the cell cycle genes upregulated at day however, five cell cycle genes, fzr ube2c, ccna nusap and cdc25b, were both downregulated at day and upregulated at day  we describe, for the first time, the gene expression profile of mouse blood cells following exposure to neutrons. we have found that neutron radiation results in both distinct and common gene expression patterns compared with x-ray radiation. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords gene expressionmicroarraysneutronx-rayradiation biodosimetryhttp://dx.doi.org/http://dx.doi.org/ /100000060national institute of allergy and infectious diseasesu19ai067773http://dx.doi.org/http://dx.doi.org/ /100000070national institute of biomedical imaging and bioengineering5p41eb-002033issue-copyright-statement© the author 2017 