BACKGROUND
a unique feature of the space radiation environment is the presence of galactic cosmic rays  and solar particle events .
the former involves protons and fully ionized atomic nuclei such as 56fe, while the latter includes predominantly low to medium energy protons with a small heavy ion component.
these radiation exposures may pose a significant hazard to deep space flight crews during and following the mission.
the hazards associated with the space environment will likely impact many organ systems.
in the central nervous system , radiation exposure significantly affects the hippocampus , a structure critical for memory function.
for example, object recognition memory  <cit> , which uses a  <dig> h interval between learning and memory assessment to test hippocampal function  <cit> , is sensitive to effects of irradiation  <cit> .
object recognition memory is impaired in 2-month-old mice two weeks following 56fe ion irradiation  <cit>  and twelve weeks following irradiation with protons  <cit> .
it is unclear whether object recognition is also sensitive to effects of charged particles in mice irradiated at  <dig> months of age, a biological age relevant to the typical ages of astronauts during space missions.
the mechanisms mediating the detrimental effects of a charged particle exposure such as 56fe ions on hippocampus-dependent cognitive function are not clear, but might be associated with changes in hippocampal networks involved in synaptic plasticity and memory.
age-related changes in immediate early gene activity-regulated cytoskeleton-associated protein  in the hippocampus have been reported  <cit> .
our work showed that head-only 56fe ion irradiation of two-month-old mice alters hippocampal expression of arc  <cit> .
arc expression has been used to provide important insight into the post-transcriptional infrastructure of gene expression involved in synaptic plasticity and memory  <cit>  .
when neurons are engaged in information processing, arc is rapidly transcribed and can be visualized and quantified after ~ <dig> min.
subsequently, the mrna is translocated to the cytoplasm where it remains detectable for ~20– <dig> min after the initial transcription.
thus, two different cellular compartments  can be clearly distinguished, providing means to identify which neurons were active during distinct behavioral experiences  <cit> .
arc mrna localization is a useful tool to test memory function.
during spatial exploration, when an animal enters a specific spatial location to which the cell responds , hippocampal pyramidal neurons in the ca areas display robust firing.
the place fields create an environment-specific "map", which is believed critical in establishing the context of a memory  <cit> .
place fields are highly reproducible and stable across multiple visits to the same environment.
however, this stability is not seen when different environments are visited.
analysis of the temporal kinetics of arc mrna allows study of neuronal activity and network stability by detecting the number of active cells forming environment-specific maps and capturing the environmental specificity of place fields  <cit> .
when animals visit the same environment twice, ~ <dig> % of cells expressing arc mrna in the nucleus also express arc mrna in the cytoplasm, which indicates that the same neurons are active during both explorations.
in contrast, when two different environments are visited, the proportion of cells showing double labeling is very low  <cit> .
in addition to network stability, other molecular mechanisms are likely involved in mediating cognitive functions.
one of these is dna methylation because it correlates closely with gene expression.
hippocampal changes in cytosine methylation, a major epigenetic modification involving the addition of a methyl group to cytosine , play a key role in regulating expression of genes required for spatial learning and memory  <cit> .
a second form of dna methylation is hydroxymethylcytosine , which is derived from 5mc by the action of three ten-eleven translocation  enzymes   <cit> .
tet <dig> is believed most important for brain function because its levels are higher than those of tet <dig> or tet <dig> in this tissue  <cit> .
the levels and localizations of 5mc and 5hmc are high and exceptionally dynamic in the brain during development and aging  <cit> .
the exact roles for 5mc and 5hmc in the brain remain to be elucidated, but their high and dynamic levels in neurons and correlation with gene expression suggest strongly that they play important roles throughout the lifespan.
we recently reported that a single exposure to proton irradiation caused long-term tissue specific changes in dna methylation patterns that were present  <dig> weeks following exposure and that mapped to tissue specific pathways   <cit> .
that study demonstrated that the brain epigenome is responsive to ionizing radiation.
in the current study, we assessed the effects of whole body 56fe ion irradiation at early  and late  time points on hippocampus-dependent memory and hippocampal network stability; and whether these effects are associated with changes in hippocampal dna methylation  and gene expression.
methods
animals and study design
six-month-old c57bl/ <dig> j male mice  were obtained from jackson laboratories, bar harbor maine.
the biological age of the mice was selected to be relevant to the biological age of astronauts during space missions.
the mice were shipped from jackson laboratories to brookhaven national laboratory , upton, new york, and allowed to accommodate to the housing facility there for one week.
subsequently, the mice were irradiated with  <dig> ,  <dig> , or  <dig>  gy of  <dig> mev 56fe or sham-irradiated .
for irradiation, mice were individually loaded into  <dig> ×  <dig> ×  <dig> cm plastic square enclosures with air holes and placed in a foam fixture in the beam line of the nasa space radiation laboratory .
they were exposed to a rectangular beam of approximately  <dig> x  <dig> cm.
the focused beam of high-energy was generated by the booster accelerator at bnl and transferred to the experimental beam line at the nsrl facility.
dose calibration was performed so that the desired dose could be delivered.
the exposure time for the  <dig>  gy dose was  <dig>  min, the exposure time for the  <dig>  gy dose was  <dig>  min, and the exposure time for the  <dig>  gy was  <dig>  min.
sham-irradiated mice were placed into the plastic enclosures for the same time as the irradiated mice.
mice were randomly assigned to the experimental groups.
the week after the irradiation or sham-irradiation, the mice were shipped to oregon health & science university  and were assigned to one of two time points .
cognitive testing started two or twenty weeks following irradiation.
following cognitive testing, the mice were killed by cervical dislocation followed by decapitation.
the hippocampus of one hemibrain was dissected for dna methylation and rna-seq analyses.
the other hemibrain was processed for arc mrna and tet <dig> immunohistochemical analyses.
all protocols were reviewed and approved by the institutional animal care and use committees  of ohsu and bnl and were in compliance with all federal regulations.
novel object recognition
the novel object recognition test was performed as described  <cit> .
the mice were habituated to an open field  for  <dig> times for  <dig> min each over three subsequent days.
on day  <dig>  the mice were placed in the open field containing two identical objects and they were allowed to freely explore for  <dig> min.
on day  <dig>  the mice were placed again in the open field, but one familiar object was replaced with a novel object.
the mice were allowed to explore for  <dig> min.
movement and time spent exploring each object was recorded using ethovision xt video tracking system  and hand scored by a researcher blinded to the treatment of the mice.
the percent time exploring the novel object, out of the total time exploring the novel and familiar objects on day  <dig>  was used to assess novel object recognition.
for each group, the preference for the novel versus the familiar object was assessed.
the open field arena and objects were cleaned with  <dig> % acetic acid between mice and trials.
water maze
spatial learning and memory were assessed in the water maze as described  <cit> .
a circular pool  was filled with water made opaque with nontoxic chalk  and mice were trained to locate a submerged platform.
to determine if irradiation affected the ability to swim or learn the water maze task, mice were first trained to locate a clearly marked platform .
the mice were subsequently trained to locate the platform when it was hidden beneath the surface of opaque water .
training during the hidden platform sessions  required the mice to learn the location of the hidden platform based on extra-maze cues.
for both visible and hidden sessions, there were two daily sessions, morning and afternoon, which were 2-h apart.
each session consisted of two trials .
a trial ended when the mice located the platform.
mice that failed to locate the platform within  <dig> s were led to the platform by placing a finger in front of their swim path.
the mice were taken out of the pool after they were physically on the platform for a minimum of  <dig> s. during visible platform training, the platform was moved to a different quadrant of the pool for each session.
for the hidden platform training, the platform location was kept constant.
the mice were placed into the water facing the edge of the pool in one of nine randomized locations.
the start location was changed for each trial.
the swimming patterns of the mice were recorded with noldus ethovision video tracking software  set at six samples/s.
the time to locate the platform  and cumulative distance to the target were used as measures of performance for the visible and hidden sessions.
because swim speeds can influence the time it takes to reach the platform, they were also analyzed.
to measure spatial memory retention, probe trials  were conducted  <dig> h after the last hidden trial of each mouse on the first and second day of hidden platform training and  <dig> h following the last hidden trial of each mouse on the third day of hidden platform training and the cumulative distance to the target location during the probe trials was analyzed.
hippocampal network stability
three days after the last water maze test day, exploration of identical or different environments was used to study the stability of hippocampal networks  <cit> .
the method called catfish  relies on the precise temporal kinetics of the ieg arc, which has been used to provide important insight into the post-transcriptional infrastructure of gene expression involved in synaptic plasticity and memory  <cit> , .
when neurons are engaged in information processing, arc is rapidly transcribed and can be visualized and quantified after ~ <dig> min.
subsequently, the mrna is translocated to the cytoplasm where it remains detectable for ~20- <dig> min after the initial transcription.
ultimately, the mrna is translocated to tagged synapses for protein synthesis.
thus, two different cellular compartments  can be clearly distinguished, providing means to identify which neurons were active during distinct behavioral experiences  <cit> .
eight mice from each experimental radiation condition were placed individually into a novel environment  and allowed to explore for five minutes, as described  <cit> .
environment a is a square open field .
after exploration, mice were returned to their cage for  <dig> min, returned to the same environment for an additional  <dig> min .
another  <dig> mice of each experimental radiation condition explored environment a for  <dig> min, and  <dig> min later they were placed in a different environment , a circular arena  <dig> cm in diameter, and allowed to explore for  <dig> min .
following the last environmental exposure, the mice were killed by cervical dislocation and the brains quickly removed, as described above.
using catfish, arc mrna appear as discrete nuclear foci , and/or as diffuse mrna in the cytoplasm .
nuclear and cytoplasmic arc can be distinguished using intronic and full-length probe respectively labeled with digoxigenin or fluorescein.
arc staining was classified as: a) none ; b) cytoplasmic arc staining only; c) nuclear arc staining only; or d) both arc-foci/arc-cyto .
we can determine if the neurons responding  to an initial experience  are the same neurons that respond to a second and identical experience as expected based on sham-irradiated animals , and if the hippocampal networks activated during exploration of the first environment are statistically independent from those activated by the exploration of a different environment .
for the ca <dig> and ca <dig> regions of the hippocampus, the percentage of neurons showing both foci and cytoplasmic staining was analyzed.
as the pattern of percentages of arc-positive neurons in both environments might reflect differences in the total number of arc-positive neurons  this number was also determined.
for the dentate gyrus, the percent of neurons expressing arc only in the cytoplasm or expressing both arc-foci/arc-cyto  was analyzed.
this measure was analyzed as it reflects the total response to the environment.
tet <dig> immunohistochemistry
following cervical dislocation and decapitation, hemibrains were rapidly removed and frozen using isopentane and dry ice, as described  <cit>  and stored at − <dig> °c.
the hemibrains were shipped to the university of california, san francisco, and cut at  <dig> μm using a cryostat.
slides were sent to ohsu for analysis of tet <dig> immunoreacttivity using a specific primary antibody from santa cruz biotechnology .
briefly, sections  were rinsed in phosphate buffered saline , and incubated in  <dig> % normal goat serum  in pbs with  <dig>  % triton x- <dig> .
next, sections were incubated overnight with primary antibodies against tet <dig>  in  <dig> % ngs in pbs-tx.
the next day, tissue sections were washed in pbs, incubated for  <dig>  h in donkey anti-rabbit alexa  <dig>  in  <dig> % ngs in pbs-tx and again rinsed in pbs.
analysis of tet <dig> immunoreactivity was performed using an olympus ix <dig> confocal microscope equipped with slidebook software.
images of hippocampal regions  and the cortex were captured within  <dig> sections  using a 20x objective .
tet <dig> immunoreactivity was quantified within fixed area frames; ca <dig> , ca <dig> , dentate gyrus , and cortex .
background threshold levels were set and applied to all images for comparison.
pixel intensities above this threshold were used for quantification measures .
the total intensity was also quantified as a measure of overall pixel intensity within a specific brain region.
dna methylation
dna was isolated from the hippocampus.
antibodies against 5mc and 5hmc were used to immunoprecipitate sonicated dna preparations for methyl-dna immunoprecipitation  and hydroxymethyl-dna immunoprecipitation .
these antibodies were used to precipitate genomic regions that are enriched for either 5mc or 5hmc.
following immunoprecipitation, high-throughput genomic sequencing and segmentation analyses were used to identify enriched genomic regions.
for dip-seq library preparation, rnase-treated dna was isolated using the qiagen allprep dna/rna protocol.
the dna was sonicated using a cole parmer cpx- <dig> sonicator  and polished using the dna terminator end repair kit .
dna fragments were a-tailed using klenow exo-  and ligated to un-methylated ht trueseq indexed adapters and purified.
the resulting purified dna was denatured at  <dig> °c, resuspended in  <dig> ul of dip ip buffer, and immunoprecipitated with  <dig> μg of the highly specific 5-methylcytosine antibody  or  <dig> ul of 5-hydroxymethylcytosine  antibody and dynal anti-mouse igg beads.
beads were rinsed  <dig> times with ip buffer, eluted with  <dig> % sds at room temperature and the eluted dna purified and subjected to limited amplification .
libraries were sequenced on the hiseq <dig> platform at the ohsu massively parallel sequencing shared resource or the oregon state university center for genome research.
dip-seq regions methylated above “background” were identified using an optimized sliding window parameter and enriched regions selected over background models via a montecarlo-permutation test  <cit> .
rna-seq
to facilitate direct comparison of dip-seq data with gene expression data, rna-seq was used to profile transcription from the same animals used for the dip-seq experiments.
rna was isolated using the nebnext poly a selection kit .
illumina high-throughput sequencing technology was used to profile rna levels in an unbiased manner.
differential methylation may occur at novel regions and the unbiased nature of rna-seq analyses enables analysis of associated un-annotated transcription.
for illumina rna-seq library preparation, we used the nebnext ultra kit according to the manufacturers specifications .
illumina data were mapped to the uc santa cruz assembly using bowtie  <cit> .
for rna-seq analyses, tags that overlap with known refseq  were counted using r scripts  <cit> .
significance was assessed using the deseq <dig> package  <cit> .
the storey q-test was used to adjust for multiple comparisons  <cit> .
libraries were sequenced on the hiseq <dig> platform at the ohsu massively parallel sequencing shared resource.
real time pcr
for real-time pcr, 20– <dig> bp primers were designed using mit's primer <dig> software  using standard parameters and all quantitation utilized standard curve real-time pcr.
primers will be provided upon request.
pcrs  contained  <dig> μl 10× pcr buffer ,  <dig>  mm mgcl <dig>   <dig> μm dntp ,  <dig> – <dig>  μm primer , 1× sybr green i , and  <dig> u platinum taq .
pcr was run on an opticon op <dig>  for one cycle at  <dig> °c,  <dig> s, and 30– <dig> cycles at  <dig> °c,  <dig> s;  <dig> °c– <dig> °c,  <dig> s. rt-pcr experiments were normalized to 18s rna levels  and were expressed as ng of gel-purified  amplicon.
bioinformatics and statistical analyses
all behavioral, cognitive, and arc data are shown as mean ± sem.
the statistical analyses of the data were performed using spss™  and graphpad prism™  software packages.
to compare effects of radiation across groups, anova was used with irradiation as the between factor, followed up by dunnett’s posthoc tests or t-tests when appropriate.
to analyze locomotor activity over three days, repeated measures anova was used.
to assess effects of radiation on arc measures across dose conditions in mice exposed to a specific environmental condition, one-way anovas were used.
to compare exploration of the objects and the percentage of arc-positive cells and total number of arc cells following exposure to the two different environmental conditions, 2-sided t-tests were used.
all figures were generated using graphpad prism software.
we considered p <  <dig>  as statistically significant.
single read sequence data was mapped to the mouse reference genome  using the bowtie algorithm using standard flags except for allowing  <dig> mismatches  <cit> .
sequences that map to a single location were selected and domains enriched for 5mc or 5hmc were selected using a parameter-optimized monte-carlo-based segmentation algorithm  <cit> .
a  <dig> bp sliding-window was selected based on iterative analyses that maximized the number of enriched regions.
a comparison of different high-throughput sequencing based methods to study dna methylation concluded that medip-seq covers ~  <dig> % of genomic cpgs  <cit> .
for statistical comparisons of biological samples, regions of methylation enrichment were merged and differences in methylation interrogated with fdr-adjusted chi-square or negative binomial statistics  <cit> .
statistical and visualization studies involved the r programming language and bioconductor packages  <cit> .
gene ontology analyses utilized the bioconductor goseq package, which adjusts for rna-seq length bias artifacts  <cit> .
for gene ontology analyses the top  <dig> dmrs  or dhrs   within  <dig> kb of a transcriptional start site were non-redundantly annotated.
unless otherwise stated, overlap between dmrs and rna-seq data was analyzed using a similar relational approach.
pathway analyses involved standard bioconductor packages .
dip sequence-tag heatmaps were generated in r by plotting median-normalized dip-seq tag density in gene bodies and indicated flanking regions with color-maps scaled to the  <dig> % quantile.
statistical analyses of pathway data were conducted via fdr-adjusted fisher exact or ks-tests.
we considered  p <  <dig>  as statistically significant.
RESULTS
effects of 56fe ion irradiation on cognitive performance
when activity levels in the open field over three days were analyzed at  <dig> and  <dig> weeks following sham-irradiation and irradiation, all groups habituated to the open field and showed higher activity levels in the open field on the first day than on subsequent days and there was a small, but significant effect of 56fe ion irradiation at some doses at the  <dig> week time point .
when object recognition was assessed at the 2-week time point, detrimental effects of 56fe ion irradiation were seen with a non-linear dose response .
mice irradiated with either  <dig>  gy or  <dig>  gy were impaired in object recognition and showed no preference for exploring the novel object.
in contrast, sham-irradiated mice and those irradiated with  <dig>  gy showed object recognition and explored the novel object significantly more than the familiar one.
interestingly, the detrimental effects of 56fe ion irradiation on object recognition were transient because all dose groups explored the novel object significantly more than the familiar one at the 20-week time point .
in contrast to object recognition, no group differences were observed for spatial learning and memory in the water maze at either time point.fig.
<dig>
a object recognition at the 2-week time point.
sham-irradiated mice and those irradiated with 56fe ions  at  <dig>  gy showed object recognition and explored the novel object more than the familiar one but mice irradiated with either  <dig>  gy or  <dig>  gy did not.
b object recognition at the 20-week time point.
all groups explored the novel object significantly more than the familiar one.
n =  <dig> mice/dose.
*p <  <dig> ; **p <  <dig>  versus familiar object
effects of 56fe ion irradiation on network stability in the ca <dig> region of the hippocampus
to determine whether the impairments in object recognition were associated with reduced ability of hippocampal neurons to recognize similar environments and to discriminate distinct environments, we assessed the effects of 56fe ion irradiation on network stability.
this ability was analyzed by assessing the percentage of arc-positive neurons expressing both arc-positive foci in the nucleus and arc staining in the cytoplasm.
this combination represents the population of neurons activated by both experiences .
at the 2-week time point, the percentage of arc-positive neurons expressing arc mrna in the nucleus and cytoplasm in the ca <dig> region of the hippocampus of sham-irradiated mice and mice irradiated with 56fe at  <dig>  gy was significantly higher following exposure twice to the same environment, as opposed to exposure to two different environments .
in contrast, the percentage of dual arc-positive neurons in the ca <dig> region of the hippocampus was not different in these two environmental conditions in mice irradiated with  <dig>   or  <dig>  gy .
thus, the arc mrna and the novel object recognition results were consistent because in both cases significant radiation effects were observed for the  <dig>  and  <dig>  gy dose groups, but not for the  <dig>  gy dose group.
there were also trends towards an effect of dose on the percentage of arc-positive cells following exposure twice to the same environments , a difference between  <dig> and  <dig>  gy  and  <dig> and  <dig>  gy , but these did not reach significance.fig.
<dig> representative images for arc catfish data in area ca <dig>  a fluorescence image showing arc expression following exploration of identical  environments .
scale bar: 50 μm.
example of neurons with foci only , cytoplasmic only  or cytoplasmic and foci arc expression  .
the arc mrna is illustrated in red and cell nuclei are indicated in blue
fig.
<dig>
a the percentage of arc-positive neurons in the ca <dig> region of the hippocampus at the 2-week time point.
the percentage of arc-positive neurons in the ca <dig> region of the hippocampus of sham-irradiated mice and mice irradiated with 56fe at  <dig>  gy was higher following exposure to twice the same environment than two different environments.
b. the total number of arc-positive neurons in the ca <dig> region of the hippocampus at the 2-week time point.
there was no significant difference in the total number of arc-positive neurons following exposure to twice the same environment than two different environments but there was a significant difference between  <dig> and  <dig>  gy.
*p <  <dig>  versus same environment; 0
p <  <dig>  versus 0 gy.
n = 5- <dig> mice/dose/environment
as the pattern of percentages of arc-positive neurons in both environments might reflect differences in the total number of arc-positive neurons , this number was also determined .
no significant difference in the total number of arc-positive neurons was observed following exposure twice to the same environment, as compared with the two different environments.
a trend towards a higher total number of arc-positive neurons was observed in sham-irradiated mice, but it did not reach significance .
however, similar to what was seen for the percentage of arc-positive neurons , there were also trends towards an effect of dose on the total number of arc-positive cells following twice exposure to the same environments , a difference between  <dig> and  <dig>  gy  and a trend towards a difference between  <dig> and  <dig>  gy .
these data indicate that group differences in the total number of arc-positive cells cannot account for the pattern of the percentage of arc-positive neurons in the ca <dig> region of the hippocampus.
the significant reduction in levels of dual cytoplasmic and nuclear arc staining observed in the ca <dig> region of sham irradiated mice exposed to different environments, as compared to the same environment twice , was not observed in the ca <dig> region of the hippocampus .fig.
<dig> the percentage of arc-positive neurons in the ca <dig> region of the hippocampus at the 2-week time point.
there were no effects of dose or trend towards group difference between dose conditions in mice following exposure to twice the same environment and following exposure to different environments.
b the total number of arc-positive neurons in the ca <dig> region of the hippocampus at the 2-week time point.
there were no effects of dose or trend towards group difference between dose conditions in mice following exposure to twice the same environment and following exposure to different environments.
n = 5- <dig> mice/dose/environment
expression of arc protein in the hippocampal dentate gyrus was altered in mice irradiated with 56fe ions at two months of age  <cit> .
to determine whether 56fe ion irradiation also affect arc mrna levels in the dentate gyrus, we analyzed the total number of arc-positive neurons in the dentate gyrus activated by both behavioral exposures.
there was no effect of radiation on the percent of neurons expressing arc only in the cytoplasm or expressing both arc-foci and showing cytoplasmic staining in the dentate gyrus following exposure to twice the same environment or two different environments .
these data suggest that the changes reported before may only occur at the translational level in de dentate gyrus.
at  <dig> weeks following 56fe irradiation, there were no effects of radiation on the percentage of arc-positive neurons or the total number of arc-positive neurons in the ca <dig> region of the hippocampus or a difference in the percentage of arc-positive neurons following exposure to twice the same environment than two different environments.
dna methylation in hippocampus following sham and 56fe irradiation
we next assessed whether 56fe ion irradiation is associated with alterations in cytosine methylation  and/or cytosine hydroxymethylation  in the hippocampus.
as there were profound differences in the effects of  <dig>  gy and  <dig>  gy 56fe ion irradiation on object recognition and network stability in the hippocampus, we compared dna methylation in these two radiation conditions with sham-irradiation.
figure  <dig> illustrates regions with significantly increased  and decreased  5hmc and  <dig> mc at the 2-  and 20-week time points .
at both time points, there were more regions with up- and downregulated 5mc at the  <dig>  than the  <dig>  gy dose.
we validated our replicate 5hmc dip-seq analysis and found that the  <dig> regions tested showed the expected directional trend and that  <dig> of  <dig> showed statistically significant increases .
dip density histograms illustrating 5mc and 5hmc signal at refseq genes sorted by rna-seq gene expression levels at the 2- and 20-week time points are shown in fig.
5b and d, respectively.
data shown are from sham-irradiated mice.
data from irradiated mice showed a similar pattern.
as expected, there was no evidence of global changes in 5mc or 5hmc signal near refseq genes.
the antibodies used to pull down 5mc and 5hmc regions do not cross react .fig.
<dig> number of genetic regions with significantly increased  and decreased  5hmc and  <dig> mc at the 2-  and 20-week  time points.
at both time points, there were more regions with up- and downregulated 5mc at the  <dig>  than  <dig>  gy.
dip density histograms illustrating 5mc and 5hmc signal at refseq genes sorted by rna-seq gene expression levels at the 2-week  and 20-week  time points.
e real-time pcr validation of significantly upregulated dip-seq regions at the  <dig>  gy 2wk time point .
real-time pcr was conducted with genomic input standard curves and  <dig> dip biological replicates.
*denotes p <  <dig>  and #denotes p <  <dig> .
f the antibodies used to pull down 5mc and 5hmc regions do not cross react
figure 6a shows the gene ontology analysis for the  <dig>  gy 56fe dose at the 2-week time point.
in general, regions with decreased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation .
blue highlighted categories likely represent a generalized, as opposed to tissue specific, response to radiation exposure  <cit> .
at the  <dig>  gy  and  <dig>  gy  doses, regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc.
at  <dig>  gy, regions with increased 5hmc had significant overlap with pathways linked to decreased 5mc, while regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc .
regions with decreased 5mc also showed less enrichment for categories linked to neuronal function or maturation.
interestingly, ontology categories involved in transcriptional regulation, cell adhesion, and cell morphogenesis were more significantly represented at  <dig>  gy, while gene categories linked to neuronal process growth and maturation were overrepresented at  <dig>  gy.
in total, the results from the 2-week time point suggest that the dynamic regulation of 5hmc near genes correlates with decreases in 5mc.fig.
<dig>
a gene ontology analysis for the  <dig>  gy 56fe dose at the 2-week time point.
regions with decreased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation.
there was an overlap in non-neuronal categories 56fe ion irradiation with that seen 22 weeks following proton irradiation.
regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc.
green highlights neuronal categories and blue highlights non-neuronal categories also observed 22 weeks following proton irradiation.
b gene ontology analysis for the  <dig>  gy 56fe dose at the 2-week time point.
regions with decreased and increased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation.
regions with decreased 5hmc were associated with neuron projection, synapse, neuron differentiation, axon guidance, and axonogenesis.
regions with increased 5hmc were associated with neuron projection development, neuron projection and differentiation, presynaptic membrane, synaptic transmission and synapse, dendrite, axon.
c gene ontology analysis at the 20-week time point.
comparing  <dig>  and  <dig>  gy, there is an overlap in neuronal categories between regions with upregulated 5hmc .
blue highlighted categories are non-neuronal categories that were also observed for the 20-week time point following proton irradiation.
the significance of pathways was lower at the 20-week than the 2-week time point.
there were no significant gene categories at an fdr <  <dig>  for 5mc.
d venn diagrams depict the overlap between significant 5hmc and 5mc dmrs  at the indicated doses and indicated time points.
a at the 2-week time point and  <dig>  gy, regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc.
in addition, regions with increase 5hmc had highly significant overlap with pathways linked to increased  <dig> mc.
b at the 2-week time point and  <dig>  gy, regions with increased 5hmc had highly significant overlap with pathways linked to decreased 5mc, while regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc.
c at the 20-week time point at  <dig>  gy, regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc and regions with increased 5hmc and significant overlap with pathways linked to regions to increased 5mc.
d at the 20-week time point at  <dig>  gy, regions with decreased 5hmc only had highly significant overlap with pathways linked to increased 5mc
the gene ontology analysis for 5hmc at the 20-week time point is shown in fig.
6c.
comparing  <dig>  and  <dig>  gy, overlap is observed in neuronal categories between regions with upregulated 5hmc .
blue highlighted categories are non-neuronal categories enriched across conditions that likely represent a generalized response to radiation exposure  <cit> .
in general, the fdr-adjusted significance of shared neuronal and non-neuronal pathways was several orders of magnitude lower at the 20-week than the 2-week time point .
nonetheless, a striking observation for 5hmc changes, particularly for the  <dig>  gy dose, is that the neuronal gene ontology pathways affected at  <dig> weeks remain altered at  <dig> weeks.
this observation demonstrates stability over time for 5hmc changes induced by the radiation exposure.
in contrast, there were no significant gene categories at an fdr-adjusted p < <dig>  for 5mc.
in total, the dna methylation data from the two doses and two time points suggest that changes in 5mc  seen at the 2-week time point were mostly transient, while 5hmc changes were mostly stable.
we next examined the relationship between individual genes associated with significant changes in 5mc or 5hmc .
genes associated with up or down-regulated 5mc regions largely showed highly significant overlap with up or down-regulated 5hmc regions.
interestingly, at the  <dig>  gy dose genes associated with upregulated 5hmc and 5mc showed the most overlap while this overlap was not significant at the  <dig>  gy dose.
it is conceivable that this relationship reflects increased tet-mediated demethylation at  <dig>  gy because regions with decreased 5mc and 5hmc had the most significant overlap at this dose.
gene expression in hippocampus following sham and 56fe irradiation
figure  <dig> illustrates heatmaps depicting significantly-regulated genes  at  <dig>   and  <dig>  gy  at the two week time point.
real-time pcr validation of significantly upregulated genes at the  <dig>  gy dose confirmed that all genes showed the expected increase with  <dig> of  <dig> significantly increased .
we examined the overlap between the top  <dig> significantly regulated rna-seq genes and top  <dig> annotated dip-seq regions.
interestingly, significant overlap at  <dig>  and  <dig>  gy doses  was only seen for regions where both 5hmc signal and rna transcription were positively correlated .
for the increased rna transcription condition, there was more enrichment in pathways for the  <dig>  gy  than the  <dig>  gy  dose.
for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases such as alzheimer’s disease, parkinson’s disease, and huntington’s disease at the  <dig>  gy  and  <dig>  gy  doses.
interestingly, the enrichment for neurodegenerative disease categories was markedly more significant at  <dig>  gy than  <dig>  gy .
the rna-seq data for the 20-week time point did not show enrichment for pathways related to neurodegenerative disease .
significantly-regulated kegg pathway data for the decreased rna transcription condition for the  <dig>  gy dose are illustrated for alzheimer’s disease , parkinson’s disease , and oxidative phosphorylation .fig.
<dig> rna-seq analysis of the 56fe ion-irradiated hippocampus at the 2-week time point.
heatmaps depict significantly-regulated genes  at the  <dig>   and  <dig>   gy dose.
c real-time pcr validation of significantly upregulated rna-seq regions at the  <dig>  gy 2wk time point .
real-time pcr was conducted with absolute standard curves and  <dig> cdna biological replicates.
* denotes p <  <dig>  and # denotes p <  <dig> .
d and e overlap of top  <dig> significantly-upregulated hmc regions and significantly-upregulated genes at the  <dig>   and  <dig>   gy dose .
at both doses, there was a highly significant overlap between regions with increased 5hmc and increased rna transcription.
f–i gene ontology for the indicated significantly differentially-regulated regions .
for the increased rna transcription condition, there was more enrichment in pathways for the  <dig>   than the  <dig>   gy dose.
for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases, such as alzheimer’s diseases, parkinson’s disease, and huntington’s disease, at the  <dig>   and  <dig>   gy doses
fig.
<dig> correlation between alterations in dna methylation and gene expression.
a and c heatmaps shows sequence density of significantly-regulated 5hmc regions sorted by rna-seq gene expression levels for the  <dig>   and  <dig>   gy dose.
note the cluster of upregulated dip-seq data in intragenic regions at the top of the plot.
the black bar denotes the position of scaled refseq intragenic regions.
b and d gene ontology analyses of the top  <dig> most up-regulated dip-seq regions  in the corresponding heatmaps of the  <dig>   and  <dig>   gy doses.
e diagram depicts a subset of genes present in the synapse gene ontology category in d. note the orthologs and preo/post-synaptic components.
f and g ucsc genome browser track depicts dip-seq sequence density at two differentially hydroxymethylated gene selected from panels b and d hatch marks may represent multiple sequences
to identify dna methylation and gene expression correlations, we sorted the density of 5mc and 5hmc regions that showed significant changes by gene expression.
interestingly, regions that showed significant increases in 5hmc at  <dig>   and  <dig>   gy clustered with highly expressed genes.
the heatmap cluster of increased 5hmc showed remarkable spatial specificity and no similar clustering was seen for other 5hmc comparisons or for 5mc .
to test whether the gene expression sorted 5hmc cluster might reflect biologically relevant pathways we performed gene ontology analyses.
we selected the top  <dig> 5hmc dip-seq regions from the heatmap cluster .
the clustered 5hmc dmrs were significantly enriched for synapse-associated go categories, with a very striking enrichment at the  <dig>  gy  as compared to  <dig>  gy  dose.
this result suggests increased synapse remodeling/repair at the  <dig>  gy dose.
the diagram in fig.
8e illustrates a subset of genes in the synapse gene ontology category at the  <dig>  gy dose.
ucsc genome browser track depicts 5hmc dip-seq data density at two differentially hydroxymethylated genes selected from the synapse gene ontology category in panels b an d .
we further explored the positive relationship between 5hmc dmrs and gene expression by analyzing the density of 5hmc and 5mc signals sorted by significantly regulated refseq genes .
this analysis utilizes the converse relationship of our analysis above and it also detected specific enrichment for increases in 5hmc in the gene-bodies of genes that were most significantly up-regulated by 56fe ion irradiation.
this was not seen for 5mc .
taken together, our results suggest that 56fe ion irradiation triggers the coordinated upregulation of intragenic 5hmc and transcription for a subset of genes that play a key role in synapse function and development.fig.
<dig> heatmaps depict the density of 5hmc and 5mc signals sorted by significantly regulated refseq genes.
note the enrichment for increases in 5hmc in the gene-bodies of genes that were most significantly up-regulated by radiation.
this was not seen for 5mc.
5’ and 3’ denote the boundaries of the scaled refseq annotation
tet <dig> immunohistochemistry
finally, we analyzed immunoreactivity of tet <dig>  one of the three tet enzymes responsible for converting 5mc to 5hmc  <cit>  and likely the most important for brain function  <cit> , in the hippocampus and cortex two weeks following 56fe ion irradiation.
in the ca <dig> region of the hippocampus, tet <dig> immunoreactivity was slightly but significantly lower following irradiation at  <dig>  gy than sham irradiation .
in addition, there was a trend towards less tet <dig> immunoreactivity in the ca <dig> of mice irradiated at  <dig>  gy than sham-irradiated mice .
although the general pattern of less tet <dig> immunoreactivity in irradiated than sham-irradiated mice was the same in the ca <dig> region of the hippocampus , dentate gyrus , and cortex , it did not reach significance.
there was also a trend towards less tet <dig> immunoreactivity in the cortex of mice irradiated at  <dig>  gy than sham-irradiated mice .
in contrast to the 2-week time point, there were no effects of 56fe irradiation on tet <dig> immunoreactivity at the 20-week time point.fig.
<dig> tet <dig> immunoreactivity two weeks following sham or 56fe-irradiation.
a representative images of tet <dig> immunoreactivity in the ca <dig> region of the hippocampus of a sham-irradiated mouse and a mouse irradiated with  <dig>  gy.
b tet <dig> immunoreactivity levels in the ca <dig> region of the hippocampus.
*p <  <dig>  versus sham-irradiation.
c tet <dig> immunoreactivity levels in the ca <dig> region of the hippocampus.
d tet <dig> immunoreactivity levels in the dentate gyrus.
e tet <dig> immunoreactivity levels in the ca <dig> region of the cortex.
n =  <dig> mice/dose.
scale bar: 100 μm
discussion
in this study, we examined the short- and long-term effects of whole body 56fe ion irradiation on hippocampus-dependent cognitive function, hippocampal dna methylation, and hippocampal rna expression.
at the two-week time point, object recognition and network stability  were impaired following irradiation at the  <dig>  and  <dig>  gy doses, but not following irradiation at the  <dig>  gy dose.
in contrast, no impairments in object recognition or network stability were seen at the 20-week time point for any of the tested doses.
consistent with this time-dependent pattern and suggesting substantial recovery, the significance of affected gene ontology for 5hmc was lower at the 20-week compared to the 2-week time point and no significant gene ontology pathways were observed for 5mc at the 20-week time point.
moreover, the 5hmc changes at the 2-week time point, but not the 5mc changes, tracked with gene expression changes detected with rna-seq.
thus, while 5hmc changes were stable and/or represented new transcriptional activity, the changes in 5mc appeared transient and not linked to transcription changes.
because 5hmc is often found at promoters and transcribed regions of active genes, unlike 5mc, which is found at the promoters of inactive genes  <cit> , these data support an important role for 5hmc in brain function following irradiation.
the results from the 2-week time point suggest that the dynamic regulation of 5hmc and 5mc are coordinated, consistent with coupled dna methylation-oxidation  and active demethylation.
both 5mc and 5hmc in the brain are relatively persistent epigenetic marks  <cit> , but our recent work showed specific responses to radiation exposure for both dna modifications  <cit> .
that observation suggests an important role for dna methylation  in brain function.
the data presented here suggest further that 5hmc plays the more critical role for the long-term response to radiation exposure.
consistent with an important role for 5hmc in hippocampal function, we observed highly significant overlap for regions with increased 5hmc and increased rna transcription at both the  <dig>  and  <dig>  gy doses, with more enrichment in pathways for the  <dig>  than the  <dig>  gy dose.
this observation is consistent with the dose-dependent pattern in cognitive injury and suggests that the greater genomic and trancriptional response to the radiation damage at  <dig>  gy prevented the cognitive impairment observed at  <dig>  gy.
for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases, such as alzheimer’s disease, parkinson’s disease , and huntington’s disease, and the decrease was greater at the  <dig>  gy dose; again consistent with increased cognitive impairment at this dose.
this effect was transient, however, because there was no significant enrichment for categories linked to neurodegenerative disease for the 20-week time point nor were cognitive impairments observed.
these data highlight the risk of transient radiation exposures leading to mission-critical cognitive impairments during space missions.
the gene ontology analyses of the top  <dig> most up-regulated dip-seq regions revealed that the clustered 5hmc dmrs were significantly enriched for synapse-associated go categories, with a more profound enrichment at the  <dig>  gy than the  <dig>  gy dose.
as noted above, this dose-dependent pattern suggests that cognitive performance in the mice was affected by the extent of the tissue response to the radiation injury.
the mice irradiated with  <dig>  gy were cognitively impaired with regard to object recognition and network stability, while those irradiated with  <dig>  gy were not.
combined, these data suggest that damage caused at the  <dig>  gy dose was sufficient to cause cognitive injury, but not sufficient to induce synapse remodeling/repair that would have mitigated the cognitive injury.
in contrast, we speculate that the higher level of damage at the  <dig>  gy dose induced a level of synapse remodeling/repair sufficient to mitigate the cognitive injury.
according to this model, cognitive changes are observed when the synapse remodeling/repair response is insufficient to repair radiation damage.
consistent with this model, down-regulated rna-seq gene pathways linked to neuronal degeneration were more significantly enriched at  <dig>  gy than at the  <dig>  gy dose.
thus, decreases in these pathways may be associated with cognitive injury and impairment.
unusual dose-dependent patterns, in which a lower dose causes more changes than a higher dose, are not restricted to 56fe ion irradiation.
when comparing the effects of 28si ion irradiation ,  <dig>  gy affected hippocampus-dependent contextual fear memory and this effect was associated with changes in synaptic plasticity not seen following irradiation at a dose of  <dig> gy  <cit> .
these unusual dose response patterns highlight the complexity of risk assessment regarding safe space radiation doses for astronauts.
the fact that the  <dig>  gy dose of 56fe ions is more detrimental for hippocampal function is especially concerning because this is a relevant dose for a deep space mission  <cit> .
the pattern of the dose response for the effects of 56fe ion irradiation on object recognition appears modulated by age at irradiation.
in the current study, in which c57bl/ <dig> j mice were irradiated with 56fe ions at  <dig> months of age to better reflect the typical age of astronauts during space missions, those irradiated with 56fe ions at  <dig>  gy were impaired in object recognition, while those irradiated at  <dig>  gy were not.
in contrast to these data, when c57bl/ <dig> j mice were irradiated at  <dig> months of age, both dose groups were impaired in object recognition  <cit> .
these data indicate that for some cognitive measures and dose exposures, older adult mice might be less sensitive to effects of 56fe ion irradiation than younger adult mice.
future efforts are warranted to assess why there is this age-dependency in sensitivity to develop cognitive injury following space irradiation.
the short- and long-term effects of 56fe ion irradiation on dna methylation could not be accounted for by radiation-induced changes in tet <dig> protein levels.
the tet <dig> protein levels were only reduced significantly in the ca <dig> region of the hippocampus following exposure to the  <dig>  gy dose.
comparing tet <dig> immunoreactivity across brain regions, the levels were particularly high in the ca <dig> region of the hippocampus and high in the dentate gyrus, with much lower levels in the ca <dig> region of the hippocampus and cortex.
the high levels of tet <dig> in the ca <dig> region of the hippocampus is interesting, as the effects of 56fe ion irradiation on network stability likely related to hippocampal place cells was also seen in this hippocampal area.
the high levels of tet <dig> in the hippocampus are consistent with other studies showing that brain tet  <dig> protein levels are higher than those of tet <dig> or tet <dig>  <cit> .
we did not see alterations in dna methylation or rna transcription of the arc gene following 56fe ion irradiation, suggesting that arc is not a direct target of epigenetic regulation under these exposure conditions.
this is not surprising because the observed changes in arc expression in the ca <dig> region of the hippocampus dealt with nuclear vs. cytoplasmic rna localization.
the consistency between the impairments in network stability, as assessed by arc mrna localization data in the ca <dig> region of the hippocampus and impairments in novel object recognition at  <dig>  but not  <dig>  gy, and the differences in dna methylation at these two doses is striking and underlines the value of dose dependency of hippocampal impairments following 56fe ion irradiation.
some of the alterations in dna methylation at the 2- and 20-week time point were not specific for 56fe ion irradiation and also seen at the 20-week time point following proton irradiation in another study  <cit> .
for example, the non-neuronal category of cell adhesion was seen as a significant pathway for  <dig>  gy 5hmc up and 5hmc down at the 2- and 20-week time points and for 5mc down at the 2-week time point.
these space irradiation-induced long-term changes in dna methylation might interact with age-dependent changes in dna methylation that are thought to be associated with the risk to develop various neurological conditions  and thereby increase the susceptibility of the brain to detrimental effects of additional radiation exposure during long duration missions.
the distinct hippocampal outcome measures used in this study showed differential susceptibility to detrimental effects of 56fe ion irradiation.
in contrast to impairments in object recognition and network stability, no cognitive changes were seen in spatial learning and memory in the water maze or in habituation to the open field.
there was a 24-h interval between training and assessment of memory retention in the object recognition test, while in the water maze training and assessment of memory retention occurred on the same day.
these data suggest that a longer time interval and potentially protein synthesis, which is required for long-term memory  <cit> , might be required to detect detrimental effects of 56fe ion irradiation on memory retention in 6-month-old mice.
however, spatial habituation learning also involved 24-h intervals and did not show profound radiation effects, while the assessments of network stability which involved two tests separated by only  <dig> min did reveal detrimental effects of 56fe ion irradiation.
in total, our cognitive results highlight the importance of considering multiple hippocampal outcome measures in assessing effects of space irradiation on hippocampal function.
CONCLUSIONS
in summary, 56fe ion irradiation at a dose as low as  <dig>  gy impairs object recognition and network stability, and is associated with short- and long-term alterations in 5hmc and short-term changes in 5mc in the hippocampus.
object recognition and network stability were not significantly affected at  <dig>  gy, a dose at which 5hmc dmrs were significantly enriched for synapse-associated gene ontology categories.
taken together, the  <dig>  and  <dig>  gy data suggest that synapse remodeling coupled with epigenomic remodeling explains how the hippocampus responds to radiation exposure, and that a potentially dangerous window of exposure exists that is sufficient to induce cognitive change without triggering a sufficient transcriptional response that would prevent this change.
our data also suggest that 5hmc remodeling is both more stable than 5mc remodeling and is highly correlated with significant changes in gene expression.
thus, our study suggests that 5hmc is a sensitive biomarker of beneficial and adverse outcome pathways.
future efforts are warranted to determine the effects of different kinds of irradiation astronauts will encounter in the space environment on hippocampus-dependent cognitive function and hippocampal dna methylation, how these changes interact, and finally how they occur.
additional files
additional file 1: figure s <dig>
a. open field activity at the 2-week time point.
there was an effect of day .
while all groups habituated to the open field and showed higher activity levels in the open field on the first day than subsequent days, the mice irradiated with 56fe ions  at  <dig>  gy moved less than sham-irradiated mice .
b. open field activity at the 20-week time point.
there was an effect of day , with higher activity levels on the first day than subsequent days.
on day  <dig>  activity levels were higher in mice irradiated with 56fe ions  at  <dig>  gy  and that on day  <dig>  mice irradiated with 56fe ions  at  <dig>  gy  or  <dig>  gy  were higher than those in sham-irradiated mice.
n =  <dig> mice/dose.
*p <  <dig>  versus sham-irradiation.
additional file 2: table s <dig>  the percentage of arc-positive neurons in the dentate gyrus at the 2-week time point <dig>
additional file 3: figure s <dig>
a. significantly regulated kegg pathway data for the decreased rna transcription condition for the  <dig>  gy dose are illustrated for alzheimer’s disease .
key molecules in ad identified included amyloid precursor protein , β-secretase , presenilin , insulin-degrading enzyme , and apolipoprotein e .
additional file 4: figure s <dig>  significantly regulated kegg pathway data for the decreased rna transcription condition for the  <dig>  gy dose are illustrated for parkinson’s disease .
key molecules in pd identified included the dopamine transporter , parkin, tyrosine hydroxylase , and molecules playing a role in mitochondrial pathways.
additional file 5: figure s <dig>  significantly regulated kegg pathway data for the decreased rna transcription condition for the  <dig>  gy dose are illustrated for oxidative phosphorylation.
key molecules identified included nadh dehydrogenase, cytochrome c oxidase, and f-type atpase.

 a unique feature of the space radiation environment is the presence of galactic cosmic rays and solar particle events. the former involves protons and fully ionized atomic nuclei such as 56fe, while the latter includes predominantly low to medium energy protons with a small heavy ion component. these radiation exposures may pose a significant hazard to deep space flight crews during and following the mission. the hazards associated with the space environment will likely impact many organ systems. in the central nervous system, radiation exposure significantly affects the hippocampus, a structure critical for memory function. for example, object recognition memory, which uses a h interval between learning and memory assessment to test hippocampal function, is sensitive to effects of irradiation. object recognition memory is impaired in 2-month-old mice two weeks following 56fe ion irradiation and twelve weeks following irradiation with protons. it is unclear whether object recognition is also sensitive to effects of charged particles in mice irradiated at months of age, a biological age relevant to the typical ages of astronauts during space missions. the mechanisms mediating the detrimental effects of a charged particle exposure such as 56fe ions on hippocampus-dependent cognitive function are not clear, but might be associated with changes in hippocampal networks involved in synaptic plasticity and memory. age-related changes in immediate early gene activity-regulated cytoskeleton-associated protein in the hippocampus have been reported. our work showed that head-only 56fe ion irradiation of two-month-old mice alters hippocampal expression of arc. arc expression has been used to provide important insight into the post-transcriptional infrastructure of gene expression involved in synaptic plasticity and memory. when neurons are engaged in information processing, arc is rapidly transcribed and can be visualized and quantified after ~ min. subsequently, the mrna is translocated to the cytoplasm where it remains detectable for ~20– min after the initial transcription. thus, two different cellular compartments can be clearly distinguished, providing means to identify which neurons were active during distinct behavioral experiences. arc mrna localization is a useful tool to test memory function. during spatial exploration, when an animal enters a specific spatial location to which the cell responds, hippocampal pyramidal neurons in the ca areas display robust firing. the place fields create an environment-specific "map", which is believed critical in establishing the context of a memory. place fields are highly reproducible and stable across multiple visits to the same environment. however, this stability is not seen when different environments are visited. analysis of the temporal kinetics of arc mrna allows study of neuronal activity and network stability by detecting the number of active cells forming environment-specific maps and capturing the environmental specificity of place fields. when animals visit the same environment twice, ~ % of cells expressing arc mrna in the nucleus also express arc mrna in the cytoplasm, which indicates that the same neurons are active during both explorations. in contrast, when two different environments are visited, the proportion of cells showing double labeling is very low. in addition to network stability, other molecular mechanisms are likely involved in mediating cognitive functions. one of these is dna methylation because it correlates closely with gene expression. hippocampal changes in cytosine methylation, a major epigenetic modification involving the addition of a methyl group to cytosine, play a key role in regulating expression of genes required for spatial learning and memory. a second form of dna methylation is hydroxymethylcytosine, which is derived from 5mc by the action of three ten-eleven translocation enzymes. tet is believed most important for brain function because its levels are higher than those of tet or tet in this tissue. the levels and localizations of 5mc and 5hmc are high and exceptionally dynamic in the brain during development and aging. the exact roles for 5mc and 5hmc in the brain remain to be elucidated, but their high and dynamic levels in neurons and correlation with gene expression suggest strongly that they play important roles throughout the lifespan. we recently reported that a single exposure to proton irradiation caused long-term tissue specific changes in dna methylation patterns that were present weeks following exposure and that mapped to tissue specific pathways. that study demonstrated that the brain epigenome is responsive to ionizing radiation. in the current study, we assessed the effects of whole body 56fe ion irradiation at early and late time points on hippocampus-dependent memory and hippocampal network stability; and whether these effects are associated with changes in hippocampal dna methylation and gene expression. methods animals and study design six-month-old c57bl/ j male mice were obtained from jackson laboratories, bar harbor maine. the biological age of the mice was selected to be relevant to the biological age of astronauts during space missions. the mice were shipped from jackson laboratories to brookhaven national laboratory, upton, new york, and allowed to accommodate to the housing facility there for one week. subsequently, the mice were irradiated with,, or gy of mev 56fe or sham-irradiated. for irradiation, mice were individually loaded into × × cm plastic square enclosures with air holes and placed in a foam fixture in the beam line of the nasa space radiation laboratory. they were exposed to a rectangular beam of approximately x cm. the focused beam of high-energy was generated by the booster accelerator at bnl and transferred to the experimental beam line at the nsrl facility. dose calibration was performed so that the desired dose could be delivered. the exposure time for the gy dose was min, the exposure time for the gy dose was min, and the exposure time for the gy was min. sham-irradiated mice were placed into the plastic enclosures for the same time as the irradiated mice. mice were randomly assigned to the experimental groups. the week after the irradiation or sham-irradiation, the mice were shipped to oregon health & science university and were assigned to one of two time points. cognitive testing started two or twenty weeks following irradiation. following cognitive testing, the mice were killed by cervical dislocation followed by decapitation. the hippocampus of one hemibrain was dissected for dna methylation and rna-seq analyses. the other hemibrain was processed for arc mrna and tet immunohistochemical analyses. all protocols were reviewed and approved by the institutional animal care and use committees of ohsu and bnl and were in compliance with all federal regulations. novel object recognition the novel object recognition test was performed as described. the mice were habituated to an open field for times for min each over three subsequent days. on day the mice were placed in the open field containing two identical objects and they were allowed to freely explore for min. on day the mice were placed again in the open field, but one familiar object was replaced with a novel object. the mice were allowed to explore for min. movement and time spent exploring each object was recorded using ethovision xt video tracking system and hand scored by a researcher blinded to the treatment of the mice. the percent time exploring the novel object, out of the total time exploring the novel and familiar objects on day was used to assess novel object recognition. for each group, the preference for the novel versus the familiar object was assessed. the open field arena and objects were cleaned with % acetic acid between mice and trials. water maze spatial learning and memory were assessed in the water maze as described. a circular pool was filled with water made opaque with nontoxic chalk and mice were trained to locate a submerged platform. to determine if irradiation affected the ability to swim or learn the water maze task, mice were first trained to locate a clearly marked platform. the mice were subsequently trained to locate the platform when it was hidden beneath the surface of opaque water. training during the hidden platform sessions required the mice to learn the location of the hidden platform based on extra-maze cues. for both visible and hidden sessions, there were two daily sessions, morning and afternoon, which were 2-h apart. each session consisted of two trials. a trial ended when the mice located the platform. mice that failed to locate the platform within s were led to the platform by placing a finger in front of their swim path. the mice were taken out of the pool after they were physically on the platform for a minimum of s. during visible platform training, the platform was moved to a different quadrant of the pool for each session. for the hidden platform training, the platform location was kept constant. the mice were placed into the water facing the edge of the pool in one of nine randomized locations. the start location was changed for each trial. the swimming patterns of the mice were recorded with noldus ethovision video tracking software set at six samples/s. the time to locate the platform and cumulative distance to the target were used as measures of performance for the visible and hidden sessions. because swim speeds can influence the time it takes to reach the platform, they were also analyzed. to measure spatial memory retention, probe trials were conducted h after the last hidden trial of each mouse on the first and second day of hidden platform training and h following the last hidden trial of each mouse on the third day of hidden platform training and the cumulative distance to the target location during the probe trials was analyzed. hippocampal network stability three days after the last water maze test day, exploration of identical or different environments was used to study the stability of hippocampal networks. the method called catfish relies on the precise temporal kinetics of the ieg arc, which has been used to provide important insight into the post-transcriptional infrastructure of gene expression involved in synaptic plasticity and memory,. when neurons are engaged in information processing, arc is rapidly transcribed and can be visualized and quantified after ~ min. subsequently, the mrna is translocated to the cytoplasm where it remains detectable for ~20- min after the initial transcription. ultimately, the mrna is translocated to tagged synapses for protein synthesis. thus, two different cellular compartments can be clearly distinguished, providing means to identify which neurons were active during distinct behavioral experiences. eight mice from each experimental radiation condition were placed individually into a novel environment and allowed to explore for five minutes, as described. environment a is a square open field. after exploration, mice were returned to their cage for min, returned to the same environment for an additional min. another mice of each experimental radiation condition explored environment a for min, and min later they were placed in a different environment, a circular arena cm in diameter, and allowed to explore for min. following the last environmental exposure, the mice were killed by cervical dislocation and the brains quickly removed, as described above. using catfish, arc mrna appear as discrete nuclear foci, and/or as diffuse mrna in the cytoplasm. nuclear and cytoplasmic arc can be distinguished using intronic and full-length probe respectively labeled with digoxigenin or fluorescein. arc staining was classified as: a) none; b) cytoplasmic arc staining only; c) nuclear arc staining only; or d) both arc-foci/arc-cyto. we can determine if the neurons responding to an initial experience are the same neurons that respond to a second and identical experience as expected based on sham-irradiated animals, and if the hippocampal networks activated during exploration of the first environment are statistically independent from those activated by the exploration of a different environment. for the ca and ca regions of the hippocampus, the percentage of neurons showing both foci and cytoplasmic staining was analyzed. as the pattern of percentages of arc-positive neurons in both environments might reflect differences in the total number of arc-positive neurons this number was also determined. for the dentate gyrus, the percent of neurons expressing arc only in the cytoplasm or expressing both arc-foci/arc-cyto was analyzed. this measure was analyzed as it reflects the total response to the environment. tet immunohistochemistry following cervical dislocation and decapitation, hemibrains were rapidly removed and frozen using isopentane and dry ice, as described and stored at − °c. the hemibrains were shipped to the university of california, san francisco, and cut at μm using a cryostat. slides were sent to ohsu for analysis of tet immunoreacttivity using a specific primary antibody from santa cruz biotechnology. briefly, sections were rinsed in phosphate buffered saline, and incubated in % normal goat serum in pbs with % triton x-. next, sections were incubated overnight with primary antibodies against tet in % ngs in pbs-tx. the next day, tissue sections were washed in pbs, incubated for h in donkey anti-rabbit alexa in % ngs in pbs-tx and again rinsed in pbs. analysis of tet immunoreactivity was performed using an olympus ix confocal microscope equipped with slidebook software. images of hippocampal regions and the cortex were captured within sections using a 20x objective. tet immunoreactivity was quantified within fixed area frames; ca, ca, dentate gyrus, and cortex. background threshold levels were set and applied to all images for comparison. pixel intensities above this threshold were used for quantification measures. the total intensity was also quantified as a measure of overall pixel intensity within a specific brain region. dna methylation dna was isolated from the hippocampus. antibodies against 5mc and 5hmc were used to immunoprecipitate sonicated dna preparations for methyl-dna immunoprecipitation and hydroxymethyl-dna immunoprecipitation. these antibodies were used to precipitate genomic regions that are enriched for either 5mc or 5hmc. following immunoprecipitation, high-throughput genomic sequencing and segmentation analyses were used to identify enriched genomic regions. for dip-seq library preparation, rnase-treated dna was isolated using the qiagen allprep dna/rna protocol. the dna was sonicated using a cole parmer cpx- sonicator and polished using the dna terminator end repair kit. dna fragments were a-tailed using klenow exo- and ligated to un-methylated ht trueseq indexed adapters and purified. the resulting purified dna was denatured at °c, resuspended in ul of dip ip buffer, and immunoprecipitated with μg of the highly specific 5-methylcytosine antibody or ul of 5-hydroxymethylcytosine antibody and dynal anti-mouse igg beads. beads were rinsed times with ip buffer, eluted with % sds at room temperature and the eluted dna purified and subjected to limited amplification. libraries were sequenced on the hiseq platform at the ohsu massively parallel sequencing shared resource or the oregon state university center for genome research. dip-seq regions methylated above “background” were identified using an optimized sliding window parameter and enriched regions selected over background models via a montecarlo-permutation test. rna-seq to facilitate direct comparison of dip-seq data with gene expression data, rna-seq was used to profile transcription from the same animals used for the dip-seq experiments. rna was isolated using the nebnext poly a selection kit. illumina high-throughput sequencing technology was used to profile rna levels in an unbiased manner. differential methylation may occur at novel regions and the unbiased nature of rna-seq analyses enables analysis of associated un-annotated transcription. for illumina rna-seq library preparation, we used the nebnext ultra kit according to the manufacturers specifications. illumina data were mapped to the uc santa cruz assembly using bowtie. for rna-seq analyses, tags that overlap with known refseq were counted using r scripts. significance was assessed using the deseq package. the storey q-test was used to adjust for multiple comparisons. libraries were sequenced on the hiseq platform at the ohsu massively parallel sequencing shared resource. real time pcr for real-time pcr, 20– bp primers were designed using mit's primer software using standard parameters and all quantitation utilized standard curve real-time pcr. primers will be provided upon request. pcrs contained μl 10× pcr buffer, mm mgcl μm dntp, – μm primer, 1× sybr green i, and u platinum taq. pcr was run on an opticon op for one cycle at °c, s, and 30– cycles at °c, s; °c– °c, s. rt-pcr experiments were normalized to 18s rna levels and were expressed as ng of gel-purified amplicon. bioinformatics and statistical analyses all behavioral, cognitive, and arc data are shown as mean ± sem. the statistical analyses of the data were performed using spss™ and graphpad prism™ software packages. to compare effects of radiation across groups, anova was used with irradiation as the between factor, followed up by dunnett’s posthoc tests or t-tests when appropriate. to analyze locomotor activity over three days, repeated measures anova was used. to assess effects of radiation on arc measures across dose conditions in mice exposed to a specific environmental condition, one-way anovas were used. to compare exploration of the objects and the percentage of arc-positive cells and total number of arc cells following exposure to the two different environmental conditions, 2-sided t-tests were used. all figures were generated using graphpad prism software. we considered p < as statistically significant. single read sequence data was mapped to the mouse reference genome using the bowtie algorithm using standard flags except for allowing mismatches. sequences that map to a single location were selected and domains enriched for 5mc or 5hmc were selected using a parameter-optimized monte-carlo-based segmentation algorithm. a bp sliding-window was selected based on iterative analyses that maximized the number of enriched regions. a comparison of different high-throughput sequencing based methods to study dna methylation concluded that medip-seq covers ~ % of genomic cpgs. for statistical comparisons of biological samples, regions of methylation enrichment were merged and differences in methylation interrogated with fdr-adjusted chi-square or negative binomial statistics. statistical and visualization studies involved the r programming language and bioconductor packages. gene ontology analyses utilized the bioconductor goseq package, which adjusts for rna-seq length bias artifacts. for gene ontology analyses the top dmrs or dhrs within kb of a transcriptional start site were non-redundantly annotated. unless otherwise stated, overlap between dmrs and rna-seq data was analyzed using a similar relational approach. pathway analyses involved standard bioconductor packages. dip sequence-tag heatmaps were generated in r by plotting median-normalized dip-seq tag density in gene bodies and indicated flanking regions with color-maps scaled to the % quantile. statistical analyses of pathway data were conducted via fdr-adjusted fisher exact or ks-tests. we considered p < as statistically significant.  effects of 56fe ion irradiation on cognitive performance when activity levels in the open field over three days were analyzed at and weeks following sham-irradiation and irradiation, all groups habituated to the open field and showed higher activity levels in the open field on the first day than on subsequent days and there was a small, but significant effect of 56fe ion irradiation at some doses at the week time point. when object recognition was assessed at the 2-week time point, detrimental effects of 56fe ion irradiation were seen with a non-linear dose response. mice irradiated with either gy or gy were impaired in object recognition and showed no preference for exploring the novel object. in contrast, sham-irradiated mice and those irradiated with gy showed object recognition and explored the novel object significantly more than the familiar one. interestingly, the detrimental effects of 56fe ion irradiation on object recognition were transient because all dose groups explored the novel object significantly more than the familiar one at the 20-week time point. in contrast to object recognition, no group differences were observed for spatial learning and memory in the water maze at either time point.fig. a object recognition at the 2-week time point. sham-irradiated mice and those irradiated with 56fe ions at  gy showed object recognition and explored the novel object more than the familiar one but mice irradiated with either  gy or  gy did not. b object recognition at the 20-week time point. all groups explored the novel object significantly more than the familiar one. n =  mice/dose. *p < ; **p <  versus familiar object effects of 56fe ion irradiation on network stability in the ca region of the hippocampus to determine whether the impairments in object recognition were associated with reduced ability of hippocampal neurons to recognize similar environments and to discriminate distinct environments, we assessed the effects of 56fe ion irradiation on network stability. this ability was analyzed by assessing the percentage of arc-positive neurons expressing both arc-positive foci in the nucleus and arc staining in the cytoplasm. this combination represents the population of neurons activated by both experiences. at the 2-week time point, the percentage of arc-positive neurons expressing arc mrna in the nucleus and cytoplasm in the ca region of the hippocampus of sham-irradiated mice and mice irradiated with 56fe at gy was significantly higher following exposure twice to the same environment, as opposed to exposure to two different environments. in contrast, the percentage of dual arc-positive neurons in the ca region of the hippocampus was not different in these two environmental conditions in mice irradiated with or gy. thus, the arc mrna and the novel object recognition results were consistent because in both cases significant radiation effects were observed for the and gy dose groups, but not for the gy dose group. there were also trends towards an effect of dose on the percentage of arc-positive cells following exposure twice to the same environments, a difference between and gy and and gy, but these did not reach significance.fig. representative images for arc catfish data in area ca a fluorescence image showing arc expression following exploration of identical environments. scale bar: 50 μm. example of neurons with foci only, cytoplasmic only or cytoplasmic and foci arc expression. the arc mrna is illustrated in red and cell nuclei are indicated in blue fig. a the percentage of arc-positive neurons in the ca region of the hippocampus at the 2-week time point. the percentage of arc-positive neurons in the ca region of the hippocampus of sham-irradiated mice and mice irradiated with 56fe at  gy was higher following exposure to twice the same environment than two different environments. b. the total number of arc-positive neurons in the ca region of the hippocampus at the 2-week time point. there was no significant difference in the total number of arc-positive neurons following exposure to twice the same environment than two different environments but there was a significant difference between and  gy. *p <  versus same environment; 0 p <  versus 0 gy. n = 5- mice/dose/environment as the pattern of percentages of arc-positive neurons in both environments might reflect differences in the total number of arc-positive neurons, this number was also determined. no significant difference in the total number of arc-positive neurons was observed following exposure twice to the same environment, as compared with the two different environments. a trend towards a higher total number of arc-positive neurons was observed in sham-irradiated mice, but it did not reach significance. however, similar to what was seen for the percentage of arc-positive neurons, there were also trends towards an effect of dose on the total number of arc-positive cells following twice exposure to the same environments, a difference between and gy and a trend towards a difference between and gy. these data indicate that group differences in the total number of arc-positive cells cannot account for the pattern of the percentage of arc-positive neurons in the ca region of the hippocampus. the significant reduction in levels of dual cytoplasmic and nuclear arc staining observed in the ca region of sham irradiated mice exposed to different environments, as compared to the same environment twice, was not observed in the ca region of the hippocampus.fig. the percentage of arc-positive neurons in the ca region of the hippocampus at the 2-week time point. there were no effects of dose or trend towards group difference between dose conditions in mice following exposure to twice the same environment and following exposure to different environments. b the total number of arc-positive neurons in the ca region of the hippocampus at the 2-week time point. there were no effects of dose or trend towards group difference between dose conditions in mice following exposure to twice the same environment and following exposure to different environments. n = 5- mice/dose/environment expression of arc protein in the hippocampal dentate gyrus was altered in mice irradiated with 56fe ions at two months of age. to determine whether 56fe ion irradiation also affect arc mrna levels in the dentate gyrus, we analyzed the total number of arc-positive neurons in the dentate gyrus activated by both behavioral exposures. there was no effect of radiation on the percent of neurons expressing arc only in the cytoplasm or expressing both arc-foci and showing cytoplasmic staining in the dentate gyrus following exposure to twice the same environment or two different environments. these data suggest that the changes reported before may only occur at the translational level in de dentate gyrus. at weeks following 56fe irradiation, there were no effects of radiation on the percentage of arc-positive neurons or the total number of arc-positive neurons in the ca region of the hippocampus or a difference in the percentage of arc-positive neurons following exposure to twice the same environment than two different environments. dna methylation in hippocampus following sham and 56fe irradiation we next assessed whether 56fe ion irradiation is associated with alterations in cytosine methylation and/or cytosine hydroxymethylation in the hippocampus. as there were profound differences in the effects of gy and gy 56fe ion irradiation on object recognition and network stability in the hippocampus, we compared dna methylation in these two radiation conditions with sham-irradiation. figure  illustrates regions with significantly increased and decreased 5hmc and mc at the 2- and 20-week time points. at both time points, there were more regions with up- and downregulated 5mc at the than the gy dose. we validated our replicate 5hmc dip-seq analysis and found that the regions tested showed the expected directional trend and that of showed statistically significant increases. dip density histograms illustrating 5mc and 5hmc signal at refseq genes sorted by rna-seq gene expression levels at the 2- and 20-week time points are shown in fig. 5b and d, respectively. data shown are from sham-irradiated mice. data from irradiated mice showed a similar pattern. as expected, there was no evidence of global changes in 5mc or 5hmc signal near refseq genes. the antibodies used to pull down 5mc and 5hmc regions do not cross react.fig. number of genetic regions with significantly increased and decreased 5hmc and mc at the 2- and 20-week time points. at both time points, there were more regions with up- and downregulated 5mc at the than  gy. dip density histograms illustrating 5mc and 5hmc signal at refseq genes sorted by rna-seq gene expression levels at the 2-week and 20-week time points. e real-time pcr validation of significantly upregulated dip-seq regions at the  gy 2wk time point. real-time pcr was conducted with genomic input standard curves and dip biological replicates. *denotes p <  and #denotes p < . f the antibodies used to pull down 5mc and 5hmc regions do not cross react figure 6a shows the gene ontology analysis for the gy 56fe dose at the 2-week time point. in general, regions with decreased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation. blue highlighted categories likely represent a generalized, as opposed to tissue specific, response to radiation exposure. at the gy and gy doses, regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc. at gy, regions with increased 5hmc had significant overlap with pathways linked to decreased 5mc, while regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc. regions with decreased 5mc also showed less enrichment for categories linked to neuronal function or maturation. interestingly, ontology categories involved in transcriptional regulation, cell adhesion, and cell morphogenesis were more significantly represented at gy, while gene categories linked to neuronal process growth and maturation were overrepresented at gy. in total, the results from the 2-week time point suggest that the dynamic regulation of 5hmc near genes correlates with decreases in 5mc.fig. a gene ontology analysis for the  gy 56fe dose at the 2-week time point. regions with decreased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation. there was an overlap in non-neuronal categories 56fe ion irradiation with that seen 22 weeks following proton irradiation. regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc. green highlights neuronal categories and blue highlights non-neuronal categories also observed 22 weeks following proton irradiation. b gene ontology analysis for the  gy 56fe dose at the 2-week time point. regions with decreased and increased 5mc were associated with pathways in neuronal categories linked to axon guidance, axogenesis, neuronal development and differentiation. regions with decreased 5hmc were associated with neuron projection, synapse, neuron differentiation, axon guidance, and axonogenesis. regions with increased 5hmc were associated with neuron projection development, neuron projection and differentiation, presynaptic membrane, synaptic transmission and synapse, dendrite, axon. c gene ontology analysis at the 20-week time point. comparing and  gy, there is an overlap in neuronal categories between regions with upregulated 5hmc. blue highlighted categories are non-neuronal categories that were also observed for the 20-week time point following proton irradiation. the significance of pathways was lower at the 20-week than the 2-week time point. there were no significant gene categories at an fdr <  for 5mc. d venn diagrams depict the overlap between significant 5hmc and 5mc dmrs at the indicated doses and indicated time points. a at the 2-week time point and  gy, regions with both decreased and increased 5hmc had highly significant overlap with pathways linked to decreased 5mc. in addition, regions with increase 5hmc had highly significant overlap with pathways linked to increased mc. b at the 2-week time point and  gy, regions with increased 5hmc had highly significant overlap with pathways linked to decreased 5mc, while regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc. c at the 20-week time point at  gy, regions with decreased 5hmc had highly significant overlap with pathways linked to increased 5mc and regions with increased 5hmc and significant overlap with pathways linked to regions to increased 5mc. d at the 20-week time point at  gy, regions with decreased 5hmc only had highly significant overlap with pathways linked to increased 5mc the gene ontology analysis for 5hmc at the 20-week time point is shown in fig. 6c. comparing and gy, overlap is observed in neuronal categories between regions with upregulated 5hmc. blue highlighted categories are non-neuronal categories enriched across conditions that likely represent a generalized response to radiation exposure. in general, the fdr-adjusted significance of shared neuronal and non-neuronal pathways was several orders of magnitude lower at the 20-week than the 2-week time point. nonetheless, a striking observation for 5hmc changes, particularly for the gy dose, is that the neuronal gene ontology pathways affected at weeks remain altered at weeks. this observation demonstrates stability over time for 5hmc changes induced by the radiation exposure. in contrast, there were no significant gene categories at an fdr-adjusted p < for 5mc. in total, the dna methylation data from the two doses and two time points suggest that changes in 5mc seen at the 2-week time point were mostly transient, while 5hmc changes were mostly stable. we next examined the relationship between individual genes associated with significant changes in 5mc or 5hmc. genes associated with up or down-regulated 5mc regions largely showed highly significant overlap with up or down-regulated 5hmc regions. interestingly, at the gy dose genes associated with upregulated 5hmc and 5mc showed the most overlap while this overlap was not significant at the gy dose. it is conceivable that this relationship reflects increased tet-mediated demethylation at gy because regions with decreased 5mc and 5hmc had the most significant overlap at this dose. gene expression in hippocampus following sham and 56fe irradiation figure  illustrates heatmaps depicting significantly-regulated genes at and gy at the two week time point. real-time pcr validation of significantly upregulated genes at the gy dose confirmed that all genes showed the expected increase with of significantly increased. we examined the overlap between the top significantly regulated rna-seq genes and top annotated dip-seq regions. interestingly, significant overlap at and gy doses was only seen for regions where both 5hmc signal and rna transcription were positively correlated. for the increased rna transcription condition, there was more enrichment in pathways for the gy than the gy dose. for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases such as alzheimer’s disease, parkinson’s disease, and huntington’s disease at the gy and gy doses. interestingly, the enrichment for neurodegenerative disease categories was markedly more significant at gy than gy. the rna-seq data for the 20-week time point did not show enrichment for pathways related to neurodegenerative disease. significantly-regulated kegg pathway data for the decreased rna transcription condition for the gy dose are illustrated for alzheimer’s disease, parkinson’s disease, and oxidative phosphorylation.fig. rna-seq analysis of the 56fe ion-irradiated hippocampus at the 2-week time point. heatmaps depict significantly-regulated genes at the and gy dose. c real-time pcr validation of significantly upregulated rna-seq regions at the  gy 2wk time point. real-time pcr was conducted with absolute standard curves and cdna biological replicates. * denotes p <  and # denotes p < . d and e overlap of top significantly-upregulated hmc regions and significantly-upregulated genes at the and gy dose. at both doses, there was a highly significant overlap between regions with increased 5hmc and increased rna transcription. f–i gene ontology for the indicated significantly differentially-regulated regions. for the increased rna transcription condition, there was more enrichment in pathways for the than the gy dose. for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases, such as alzheimer’s diseases, parkinson’s disease, and huntington’s disease, at the and gy doses fig. correlation between alterations in dna methylation and gene expression. a and c heatmaps shows sequence density of significantly-regulated 5hmc regions sorted by rna-seq gene expression levels for the and gy dose. note the cluster of upregulated dip-seq data in intragenic regions at the top of the plot. the black bar denotes the position of scaled refseq intragenic regions. b and d gene ontology analyses of the top most up-regulated dip-seq regions in the corresponding heatmaps of the and gy doses. e diagram depicts a subset of genes present in the synapse gene ontology category in d. note the orthologs and preo/post-synaptic components. f and g ucsc genome browser track depicts dip-seq sequence density at two differentially hydroxymethylated gene selected from panels b and d hatch marks may represent multiple sequences to identify dna methylation and gene expression correlations, we sorted the density of 5mc and 5hmc regions that showed significant changes by gene expression. interestingly, regions that showed significant increases in 5hmc at and gy clustered with highly expressed genes. the heatmap cluster of increased 5hmc showed remarkable spatial specificity and no similar clustering was seen for other 5hmc comparisons or for 5mc. to test whether the gene expression sorted 5hmc cluster might reflect biologically relevant pathways we performed gene ontology analyses. we selected the top 5hmc dip-seq regions from the heatmap cluster. the clustered 5hmc dmrs were significantly enriched for synapse-associated go categories, with a very striking enrichment at the gy as compared to gy dose. this result suggests increased synapse remodeling/repair at the gy dose. the diagram in fig. 8e illustrates a subset of genes in the synapse gene ontology category at the gy dose. ucsc genome browser track depicts 5hmc dip-seq data density at two differentially hydroxymethylated genes selected from the synapse gene ontology category in panels b an d. we further explored the positive relationship between 5hmc dmrs and gene expression by analyzing the density of 5hmc and 5mc signals sorted by significantly regulated refseq genes. this analysis utilizes the converse relationship of our analysis above and it also detected specific enrichment for increases in 5hmc in the gene-bodies of genes that were most significantly up-regulated by 56fe ion irradiation. this was not seen for 5mc. taken together, our results suggest that 56fe ion irradiation triggers the coordinated upregulation of intragenic 5hmc and transcription for a subset of genes that play a key role in synapse function and development.fig. heatmaps depict the density of 5hmc and 5mc signals sorted by significantly regulated refseq genes. note the enrichment for increases in 5hmc in the gene-bodies of genes that were most significantly up-regulated by radiation. this was not seen for 5mc. 5’ and 3’ denote the boundaries of the scaled refseq annotation tet immunohistochemistry finally, we analyzed immunoreactivity of tet one of the three tet enzymes responsible for converting 5mc to 5hmc and likely the most important for brain function, in the hippocampus and cortex two weeks following 56fe ion irradiation. in the ca region of the hippocampus, tet immunoreactivity was slightly but significantly lower following irradiation at gy than sham irradiation. in addition, there was a trend towards less tet immunoreactivity in the ca of mice irradiated at gy than sham-irradiated mice. although the general pattern of less tet immunoreactivity in irradiated than sham-irradiated mice was the same in the ca region of the hippocampus, dentate gyrus, and cortex, it did not reach significance. there was also a trend towards less tet immunoreactivity in the cortex of mice irradiated at gy than sham-irradiated mice. in contrast to the 2-week time point, there were no effects of 56fe irradiation on tet immunoreactivity at the 20-week time point.fig. tet immunoreactivity two weeks following sham or 56fe-irradiation. a representative images of tet immunoreactivity in the ca region of the hippocampus of a sham-irradiated mouse and a mouse irradiated with  gy. b tet immunoreactivity levels in the ca region of the hippocampus. *p <  versus sham-irradiation. c tet immunoreactivity levels in the ca region of the hippocampus. d tet immunoreactivity levels in the dentate gyrus. e tet immunoreactivity levels in the ca region of the cortex. n =  mice/dose. scale bar: 100 μm discussion in this study, we examined the short- and long-term effects of whole body 56fe ion irradiation on hippocampus-dependent cognitive function, hippocampal dna methylation, and hippocampal rna expression. at the two-week time point, object recognition and network stability were impaired following irradiation at the and gy doses, but not following irradiation at the gy dose. in contrast, no impairments in object recognition or network stability were seen at the 20-week time point for any of the tested doses. consistent with this time-dependent pattern and suggesting substantial recovery, the significance of affected gene ontology for 5hmc was lower at the 20-week compared to the 2-week time point and no significant gene ontology pathways were observed for 5mc at the 20-week time point. moreover, the 5hmc changes at the 2-week time point, but not the 5mc changes, tracked with gene expression changes detected with rna-seq. thus, while 5hmc changes were stable and/or represented new transcriptional activity, the changes in 5mc appeared transient and not linked to transcription changes. because 5hmc is often found at promoters and transcribed regions of active genes, unlike 5mc, which is found at the promoters of inactive genes, these data support an important role for 5hmc in brain function following irradiation. the results from the 2-week time point suggest that the dynamic regulation of 5hmc and 5mc are coordinated, consistent with coupled dna methylation-oxidation and active demethylation. both 5mc and 5hmc in the brain are relatively persistent epigenetic marks, but our recent work showed specific responses to radiation exposure for both dna modifications. that observation suggests an important role for dna methylation in brain function. the data presented here suggest further that 5hmc plays the more critical role for the long-term response to radiation exposure. consistent with an important role for 5hmc in hippocampal function, we observed highly significant overlap for regions with increased 5hmc and increased rna transcription at both the and gy doses, with more enrichment in pathways for the than the gy dose. this observation is consistent with the dose-dependent pattern in cognitive injury and suggests that the greater genomic and trancriptional response to the radiation damage at gy prevented the cognitive impairment observed at gy. for the decreased rna transcription condition, there was a highly significant enrichment for categories linked to neurodegenerative diseases, such as alzheimer’s disease, parkinson’s disease, and huntington’s disease, and the decrease was greater at the gy dose; again consistent with increased cognitive impairment at this dose. this effect was transient, however, because there was no significant enrichment for categories linked to neurodegenerative disease for the 20-week time point nor were cognitive impairments observed. these data highlight the risk of transient radiation exposures leading to mission-critical cognitive impairments during space missions. the gene ontology analyses of the top most up-regulated dip-seq regions revealed that the clustered 5hmc dmrs were significantly enriched for synapse-associated go categories, with a more profound enrichment at the gy than the gy dose. as noted above, this dose-dependent pattern suggests that cognitive performance in the mice was affected by the extent of the tissue response to the radiation injury. the mice irradiated with gy were cognitively impaired with regard to object recognition and network stability, while those irradiated with gy were not. combined, these data suggest that damage caused at the gy dose was sufficient to cause cognitive injury, but not sufficient to induce synapse remodeling/repair that would have mitigated the cognitive injury. in contrast, we speculate that the higher level of damage at the gy dose induced a level of synapse remodeling/repair sufficient to mitigate the cognitive injury. according to this model, cognitive changes are observed when the synapse remodeling/repair response is insufficient to repair radiation damage. consistent with this model, down-regulated rna-seq gene pathways linked to neuronal degeneration were more significantly enriched at gy than at the gy dose. thus, decreases in these pathways may be associated with cognitive injury and impairment. unusual dose-dependent patterns, in which a lower dose causes more changes than a higher dose, are not restricted to 56fe ion irradiation. when comparing the effects of 28si ion irradiation, gy affected hippocampus-dependent contextual fear memory and this effect was associated with changes in synaptic plasticity not seen following irradiation at a dose of gy. these unusual dose response patterns highlight the complexity of risk assessment regarding safe space radiation doses for astronauts. the fact that the gy dose of 56fe ions is more detrimental for hippocampal function is especially concerning because this is a relevant dose for a deep space mission. the pattern of the dose response for the effects of 56fe ion irradiation on object recognition appears modulated by age at irradiation. in the current study, in which c57bl/ j mice were irradiated with 56fe ions at months of age to better reflect the typical age of astronauts during space missions, those irradiated with 56fe ions at gy were impaired in object recognition, while those irradiated at gy were not. in contrast to these data, when c57bl/ j mice were irradiated at months of age, both dose groups were impaired in object recognition. these data indicate that for some cognitive measures and dose exposures, older adult mice might be less sensitive to effects of 56fe ion irradiation than younger adult mice. future efforts are warranted to assess why there is this age-dependency in sensitivity to develop cognitive injury following space irradiation. the short- and long-term effects of 56fe ion irradiation on dna methylation could not be accounted for by radiation-induced changes in tet protein levels. the tet protein levels were only reduced significantly in the ca region of the hippocampus following exposure to the gy dose. comparing tet immunoreactivity across brain regions, the levels were particularly high in the ca region of the hippocampus and high in the dentate gyrus, with much lower levels in the ca region of the hippocampus and cortex. the high levels of tet in the ca region of the hippocampus is interesting, as the effects of 56fe ion irradiation on network stability likely related to hippocampal place cells was also seen in this hippocampal area. the high levels of tet in the hippocampus are consistent with other studies showing that brain tet protein levels are higher than those of tet or tet. we did not see alterations in dna methylation or rna transcription of the arc gene following 56fe ion irradiation, suggesting that arc is not a direct target of epigenetic regulation under these exposure conditions. this is not surprising because the observed changes in arc expression in the ca region of the hippocampus dealt with nuclear vs. cytoplasmic rna localization. the consistency between the impairments in network stability, as assessed by arc mrna localization data in the ca region of the hippocampus and impairments in novel object recognition at but not gy, and the differences in dna methylation at these two doses is striking and underlines the value of dose dependency of hippocampal impairments following 56fe ion irradiation. some of the alterations in dna methylation at the 2- and 20-week time point were not specific for 56fe ion irradiation and also seen at the 20-week time point following proton irradiation in another study. for example, the non-neuronal category of cell adhesion was seen as a significant pathway for gy 5hmc up and 5hmc down at the 2- and 20-week time points and for 5mc down at the 2-week time point. these space irradiation-induced long-term changes in dna methylation might interact with age-dependent changes in dna methylation that are thought to be associated with the risk to develop various neurological conditions and thereby increase the susceptibility of the brain to detrimental effects of additional radiation exposure during long duration missions. the distinct hippocampal outcome measures used in this study showed differential susceptibility to detrimental effects of 56fe ion irradiation. in contrast to impairments in object recognition and network stability, no cognitive changes were seen in spatial learning and memory in the water maze or in habituation to the open field. there was a 24-h interval between training and assessment of memory retention in the object recognition test, while in the water maze training and assessment of memory retention occurred on the same day. these data suggest that a longer time interval and potentially protein synthesis, which is required for long-term memory, might be required to detect detrimental effects of 56fe ion irradiation on memory retention in 6-month-old mice. however, spatial habituation learning also involved 24-h intervals and did not show profound radiation effects, while the assessments of network stability which involved two tests separated by only min did reveal detrimental effects of 56fe ion irradiation. in total, our cognitive results highlight the importance of considering multiple hippocampal outcome measures in assessing effects of space irradiation on hippocampal function.  in summary, 56fe ion irradiation at a dose as low as gy impairs object recognition and network stability, and is associated with short- and long-term alterations in 5hmc and short-term changes in 5mc in the hippocampus. object recognition and network stability were not significantly affected at gy, a dose at which 5hmc dmrs were significantly enriched for synapse-associated gene ontology categories. taken together, the and gy data suggest that synapse remodeling coupled with epigenomic remodeling explains how the hippocampus responds to radiation exposure, and that a potentially dangerous window of exposure exists that is sufficient to induce cognitive change without triggering a sufficient transcriptional response that would prevent this change. our data also suggest that 5hmc remodeling is both more stable than 5mc remodeling and is highly correlated with significant changes in gene expression. thus, our study suggests that 5hmc is a sensitive biomarker of beneficial and adverse outcome pathways. future efforts are warranted to determine the effects of different kinds of irradiation astronauts will encounter in the space environment on hippocampus-dependent cognitive function and hippocampal dna methylation, how these changes interact, and finally how they occur. additional files additional file 1: figure s a. open field activity at the 2-week time point. there was an effect of day. while all groups habituated to the open field and showed higher activity levels in the open field on the first day than subsequent days, the mice irradiated with 56fe ions at gy moved less than sham-irradiated mice. b. open field activity at the 20-week time point. there was an effect of day, with higher activity levels on the first day than subsequent days. on day activity levels were higher in mice irradiated with 56fe ions at gy and that on day mice irradiated with 56fe ions at gy or gy were higher than those in sham-irradiated mice. n = mice/dose. *p < versus sham-irradiation. additional file 2: table s the percentage of arc-positive neurons in the dentate gyrus at the 2-week time point additional file 3: figure s a. significantly regulated kegg pathway data for the decreased rna transcription condition for the gy dose are illustrated for alzheimer’s disease. key molecules in ad identified included amyloid precursor protein, β-secretase, presenilin, insulin-degrading enzyme, and apolipoprotein e. additional file 4: figure s significantly regulated kegg pathway data for the decreased rna transcription condition for the gy dose are illustrated for parkinson’s disease. key molecules in pd identified included the dopamine transporter, parkin, tyrosine hydroxylase, and molecules playing a role in mitochondrial pathways. additional file 5: figure s significantly regulated kegg pathway data for the decreased rna transcription condition for the gy dose are illustrated for oxidative phosphorylation. key molecules identified included nadh dehydrogenase, cytochrome c oxidase, and f-type atpase. 
moreover, because the mechanisms mediating the detrimental effects of 56fe ions on hippocampal function are unclear, we examined changes in hippocampal networks involved in synaptic plasticity and memory, gene expression, and epigenetic changes in cytosine methylation and hydroxymethylation that could accompany changes in gene expression. consistent with this pattern, the significance of pathways for gene categories for 5hmc was lower, though not eliminated, at the 20-week time point compared to the 2-week time point. similarly, significant changes were observed for 5mc gene pathways at the 2-week time point, but no significant gene categories were observed at the 20-week time point. only the 5hmc changes tracked with gene expression changes. we assessed the effects of whole body 56fe ion irradiation at early and late time points on hippocampus-dependent memory and hippocampal network stability, and whether these effects are associated with epigenetic changes in hippocampal dna methylation and gene expression. at the two-week time point, object recognition and network stability were impaired following irradiation at the and gy dose, but not following irradiation at the gy dose. no impairments in object recognition or network stability were seen at the 20-week time point at any irradiation dose used. hippocampusdna methylationcognitionnetwork stabilitygene expressionhttp://dx.doi.org/ /100000104national aeronautics and space administrationnnj12zsa001nraber jacob issue-copyright-statement© the author 2016 here we examined object recognition in 6-month-old mice irradiated with 56fe ions, a biological age more relevant to the typical ages of astronauts. astronauts are exposed to 56fe ions that may pose a significant health hazard during and following prolonged missions in deep space. we showed previously that object recognition requiring the hippocampus, a structure critical for cognitive function, is affected in 2-month-old mice irradiated with 56fe ions.
 astronauts are exposed to 56fe ions that may pose a significant health hazard during and following prolonged missions in deep space. we showed previously that object recognition requiring the hippocampus, a structure critical for cognitive function, is affected in 2-month-old mice irradiated with 56fe ions. here we examined object recognition in 6-month-old mice irradiated with 56fe ions, a biological age more relevant to the typical ages of astronauts. moreover, because the mechanisms mediating the detrimental effects of 56fe ions on hippocampal function are unclear, we examined changes in hippocampal networks involved in synaptic plasticity and memory, gene expression, and epigenetic changes in cytosine methylation and hydroxymethylation that could accompany changes in gene expression. we assessed the effects of whole body 56fe ion irradiation at early and late time points on hippocampus-dependent memory and hippocampal network stability, and whether these effects are associated with epigenetic changes in hippocampal dna methylation and gene expression.  at the two-week time point, object recognition and network stability were impaired following irradiation at the and gy dose, but not following irradiation at the gy dose. no impairments in object recognition or network stability were seen at the 20-week time point at any irradiation dose used. consistent with this pattern, the significance of pathways for gene categories for 5hmc was lower, though not eliminated, at the 20-week time point compared to the 2-week time point. similarly, significant changes were observed for 5mc gene pathways at the 2-week time point, but no significant gene categories were observed at the 20-week time point. only the 5hmc changes tracked with gene expression changes.  dose- and time-dependent epigenomic remodeling in the hippocampus following 56fe ion exposure correlates with behavioral changes. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords hippocampusdna methylationcognitionnetwork stabilitygene expressionhttp://dx.doi.org/ /100000104national aeronautics and space administrationnnj12zsa001nraber jacob issue-copyright-statement© the author 2016 