Impact of inflammatory signaling on radiation biodosimetry: mouse model of inflammatory bowel disease

Background Ionizing Radiation (IR) is a known pro-inflammatory agent and in the process of development of biomarkers for radiation biodosimetry, a chronic inflammatory disease condition could act as a confounding factor. Hence, it is important to develop radiation signatures that can distinguish between IR-induced inflammatory responses and pre-existing disease. In this study, we compared the gene expression response of a genetically modified mouse model of inflammatory bowel disease (Il10−/−) with that of a normal wild-type mouse to potentially develop transcriptomics-based biodosimetry markers that can predict radiation exposure in individuals regardless of pre-existing inflammatory condition. Results Wild-type (WT) and Il10−/− mice were exposed to whole body irradiation of 7 Gy X-rays. Gene expression responses were studied using high throughput whole genome microarrays in peripheral blood 24 h post-irradiation. Analysis resulted in identification of 1962 and 1844 genes differentially expressed (p < 0.001, FDR < 10%) after radiation exposure in Il10−/− and WT mice respectively. A set of 155 genes was also identified as differentially expressed between WT and Il10−/− mice at the baseline pre-irradiation level. Gene ontology analysis revealed that the 155 baseline differentially expressed genes were mainly involved in inflammatory response, glutathione metabolism and collagen deposition. Analysis of radiation responsive genes revealed that innate immune response and p53 signaling processes were strongly associated with up-regulated genes, whereas B-cell development process was found to be significant amongst downregulated genes in the two genotypes. However, specific immune response pathways like MHC based antigen presentation, interferon signaling and hepatic fibrosis were associated with radiation responsive genes in Il10−/− mice but not WT mice. Further analysis using the IPA prediction tool revealed significant differences in the predicted activation status of T-cell mediated signaling as well as regulators of inflammation between WT and Il10−/− after irradiation. Conclusions Using a mouse model we established that an inflammatory disease condition could affect the expression of many radiation responsive genes. Nevertheless, we identified a panel of genes that, regardless of disease condition, could predict radiation exposure. Our results highlight the need for consideration of pre-existing conditions in the population in the process of development of reliable biodosimetry markers. Electronic supplementary material The online version of this article (10.1186/s12864-019-5689-y) contains supplementary material, which is available to authorized users.


Background
Chronic inflammation is a part of many human diseases [1]. Studies reveal that approximately 5-7% of the western population suffers from some sort of chronic inflammatory condition [1,2]. Ionizing radiation is a known pro-inflammatory agent [3]. Thus, pre-existing inflammatory disease could potentially confound radiation biodosimetry approaches used for medical triage during a radiological event.
The inflammatory aspect of radiation has been studied extensively in recent years [4][5][6][7]. Radiation induced inflammation has previously been reported in the atomic bomb survivors [8]. Lately, there has been renewed interest in the role of inflammation in radiotherapy [5,9]. On one hand, radiation-mediated immunotoxicity is one of the limiting factors for radiotherapy while on the other hand, radiation facilitated increased immunogenicity is being used to enhance tumor treatment [9][10][11]. Various factors govern the degree of inflammation after radiation exposure. Genetic background of an individual and preexisting disease conditions can greatly affect the outcome of a radiation exposure [12][13][14].
Any accidental or intentional large-scale radiological event will result in a large number of individuals requiring medical attention, which will require an effective biodosimetry tool for efficient distribution of medical resources [15,16]. An individual with an underlying inflammatory condition could show heightened immune responses to radiation exposure and thus obscure some biodosimetry methods. Thus, it is necessary to consider the immune modulating effects of radiation and the potential impact of underlying inflammatory conditions in the development of effective biodosimetry tools for triage purposes.
Various biodosimetric methods have been proposed and evaluated for large-scale radiation exposure scenarios [16]. Gene expression represents an emerging approach to biodosimetry, and could provide estimates of absorbed dose and severity of radiation-induced injury [17][18][19]. It has been shown that genes involved in inflammation and immune response are among the strongest early responders following ionizing radiation exposure [18,19]. In addition, radiation exposure triggers a long-lived inflammatory response, with up regulation of inflammatory markers and immune dysfunction that can persist for decades [3,8].
In this study, we used Il10 −/− , a mouse model of human inflammatory bowel disease (IBD), and compared its transcriptomic responses to radiation with those of normal wild-type (WT) mice. Interleukin (IL) 10 is a potent anti-inflammatory cytokine and plays a significant role in regulating the expression of Th1 cytokines, MHC mediated antigen presentation and antibody production from B-cells [20]. After radiation exposure, this cytokine plays an important role in regulating inflammation as well as free radical mediated damage in the neighboring non-exposed cells and tissues [21,22]. In the absence of IL-10 cytokine, the Il10 −/− mice become highly susceptible to enteric bacterial pathogens, and show a heightened inflammatory response to pathogen infection [23,24]. Similar to IBD patients, the Il10 −/− mice have been shown to have higher levels of IL23-producing macrophages, which through their interaction with circulating T-cells cause production of inflammatory cytokines [25]. In terms of radiation responses, Il10 −/− mice have been shown to readily develop colitis following a smaller dose of radiation compared to normal WT-mice [Unpublished data]. Thus, we have used this genotype to represent those individuals in the population who have an underlying chronic inflammatory condition that could potentially skew the results from gene expression based biodosimetry during mass triage after a radiation emergency. This study was undertaken as an initial investigation into the potential impact of chronic inflammation on gene expression based radiation biodosimetry markers.

Animals and irradiation
Wild-type (WT) C57BL/6 and Il10 −/− mice (B6.129P2-Il10 tm1Cgn /J) were obtained from Jackson Laboratory and bred at Georgetown University as het x het. The study was conducted according to the Georgetown University Institutional Animal Care and Use Committee (IACUC). All the mice (including the WT) were housed in the same room with similar bedding, food, and water under specific pathogen-free (SPF) housing conditions. In addition, the experiment was conducted in the same time for all mice, minimizing the potential of microbiome alterations in the facility. Furthermore, no other mice besides the ones in this study were housed in that room during the time of the experiment. Male littermates between 8 and 10 weeks of age were whole body irradiated with 7 Gy X-rays using an X-Rad 320 X-ray machine (Precision X-ray Inc., Brandford, CT). Sham-irradiated control mice (0 Gy) were treated the same as the exposed animals except the source was not switched on. The dose was selected based on the LD 50/30 value of WT-C57BL/6 mice (background strain for Il10 −/− mice), of 8Gy for X-ray irradiation [26]. It was reduced slightly (7 Gy) as the Il10 −/− mice appeared in early experiments to be more sensitive to radiation. Mice were euthanized 24 h after irradiation using carbon dioxide inhalation followed by cervical dislocation, per IACUC standard procedures. Approximately~0.4 ml blood was collected from each mouse using cardiac puncture method, into a 15 ml centrifuge tube containing 1.6 ml PAXgene Blood RNA stabilization and lysis solution (PreAnalytix GmBH, catalog # 762165), and mixed thoroughly. The tubes were then frozen at − 80°C , and shipped to Columbia University for RNA isolation. At the time of sacrifice, a drop of blood was also collected in EDTA Microtainer® tubes (BD medical, catalog #365974) for complete blood counts (CBC) using a Genesis hematology system (Oxford Science). Blood samples from 8 animals were used for CBC analysis.

RNA isolation
Blood samples in PAXgene solution were allowed to reach room temperature for 2 h before proceeding to RNA isolation. RNA was purified following the PAXgene RNA kit recommendations with on-column DNase I treatment. Globin RNA was reduced using the Ambion GLOBINclear-mouse/rat kit (Thermo Fisher Scientific). RNA yields were quantified using the NanoDrop ND1000 spectrophotometer (Thermo Fisher Scientific) and RNA quality was checked by the 2100 Bioanalyzer (Agilent). High quality RNA with an RNA integrity number of at least 7.0 was used for microarray hybridization.

Microarray hybridization
For both genotypes, RNA isolated from control and irradiated mouse blood samples (n = 5) was labeled with Cyanine-3 (Cy3) using the One-Color Low-Input Quick Amp Labeling Kit (Agilent Technologies, Santa Clara, CA) following manufacturer's instructions. Study size (n = 5) was selected to provide 90% power to detect two fold changes, based on an analysis of our previous mouse microarray data using the methods of Lee and Whitmore [27]. A total of 100 ng of input RNA was used for each labeling process. The labeled RNA was purified using the RNAeasy mini kit (Qiagen). Purified RNA having a Cy3 specific activity more than 8 was used for microarray hybridization. A total of 1.65 μg of labeled RNA for each sample was fragmented and hybridized to mouse GE4x44 v2 microarrays (part no. G4846A; Agilent Technologies). The hybridization was carried out at 65°C for 17 h in a hybridization oven followed by washing as recommended by the supplier (Agilent). Slides were scanned with the Agilent DNA microarray scanner (G2505B) and the images were analyzed with Feature Extraction software [v10.7] (Agilent) using default parameters for background correction and flagging non-uniform features.

Microarray data analysis
Background corrected hybridization intensities were imported into BRB-ArrayTools, v. 4.3.2 (NCI, Biometric Research Branch, Bethesda, MD) [28] log2-transformed and median normalized. Non-uniform outliers or features not significantly above background intensity in 25% or more of the samples, as well as features not changing at least 1.5-fold in 20% or more of the samples were filtered out, giving 16,000 features that were used in subsequent analyses. The microarray data is available through the NCBI Gene Expression Omnibus (series no. GSE114142; http://www.ncbi.nlm.nih.gov/geo/query/acc. cgi?acc=GSE114142).
BRB-Array Tools was used to identify genes that were differentially expressed after radiation exposure in WT and Il10 −/− mice. Genes with P < 0.001 and having a false-discovery rate of less than 10% as calculated using the Benjamini-Hochberg method were considered significantly differentially expressed.

Class prediction
We used class prediction methods in BRB-Array Tools to select genes and build predictors of radiation exposure status. The greedy pairs method [29] for feature selection was used to identify the 12 top-performing pairs of genes for discrimination between irradiated and control samples using the wild-type data as the training set. Seven classification methods (compound covariate predictor, linear discriminant analysis, 1-and 3-nearest neighbors, nearest centroid, support vector machines, and Bayesian compound covariate predictor) were used with the selected feature sets to predict the irradiation status of the remaining test samples (Il10 −/− ). The percentage correct classification was calculated for each approach. The same approach was then repeated using the Il10 −/− data as the training set and the wild-type data as the test set.

Gene ontology and network analysis
The lists of differentially expressed genes were analyzed using the functional annotation tool of the Database for Annotation Visualization and Integrated Discovery (DA-VID; v 6.7) [30]. Gene ontology terms and biological functions with a Benjamini-corrected P value < 0.05 were considered significantly over-represented within a gene list. Genes that were significantly differentially expressed in the two mouse strains were also imported into Ingenuity Pathway Analysis (IPA from QIAGEN Inc., www.qiagen.com/ingenuity) software and analyzed using IPA Core analysis and Comparative Analysis Tools. IPA uses curated information on the published relationships between gene products to predict networks and association between genes in a list. The upstream regulator analysis specifically uses information about the relationship between the activity of potential upstream regulatory factors and the expression changes of the measured genes to make predictions on the regulatory status of the upstream molecule. IPA generates a z-score for each factor in the upstream regulator analysis and for prediction of activation or inhibition state of biological functions. The IPA default cutoff of z ≥ 2 was used to predict activation and z ≤ − 2 to predict inhibition.

Quantitative real-time RT-PCR
Real-time quantitative RT-PCR (qRT-PCR) was performed for selected genes using Taqman chemistry and the ABI 7900 Real Time PCR System. The Globin cleared purified RNA from 5 control and irradiated animals was used for cDNA synthesis using the High-Capacity cDNA Archive Kit (Life Technologies). Gene expression assays (primer/ probe sets) were purchased from Thermo Fisher for the following genes: The ΔΔCT method was used to calculate expression relative to controls, using normalization to Actb expression. Each reaction was run in triplicate with 5 control and 5 irradiated samples, and means were compared using an unpaired t-test.

Results
Effect of radiation on blood cell counts of WT and Il10 −/− mice In the present study, we compared gene expression in the peripheral blood of WT and Il10 −/− mice with and without exposure to ionizing radiation. As different blood cell types have different susceptibility to radiation, we first analyzed the effect of radiation on blood cell counts in the two genotypes 24 h post irradiation.
The total number of white blood cells (WBC) and relative percentages of major blood cell population (i.e. lymphocytes, neutrophils and monocytes) is shown in Fig. 1. Overall, we found a decrease in WBC counts after radiation exposure, which was significant in WT, but not in Il10 −/− , mice (Fig. 1a). Although post-irradiation WBC numbers were similar in both genotypes, the Il10 −/− controls showed a higher degree of variability. Analysis of different blood cellular components revealed significant decreases in the lymphocyte percentage in the blood of both WT and Il10 −/− mice after radiation exposure (Fig. 1b). There was an increase in neutrophil percentage 24 h after radiation exposure, which was statistically significant in WT mice (P < 0.01) but not in Il10 −/− mice (Fig. 1c). There was no significant difference in the blood monocyte percentages after radiation exposure in either genotype (Fig. 1d). Interestingly, we observed significant differences in the baseline percentage of neutrophils between WT and Il10 −/− mice pre-irradiation (Fig. 1c). In a separate experiment, we found that at a lower dose (2Gy) of radiation exposure, most of the decrease in the WBC counts had been recovered by Day 30 post irradiation in both genotypes (Additional file 1). We did not find any significant differences in red blood cell (RBC) and platelet counts after radiation exposure in WT or Il10 −/− mice.
Analysis using the IPA network analysis tool [QIAGEN Inc., www.qiagen.com/ingenuity] revealed that LPS (lipopolysaccharide)/IL-1 signaling was the most significantly enriched (P = 7.23 × 10 − 8 ) canonical pathway in this 155 gene set ( Fig. 2b; Additional file 2C). This represents a broad inflammatory response, as well as a central pathway linking other responses, including inflammation and glutathione metabolism, that were most affected in the gene ontology analysis of these genes (Fig. 2c). Some of these baseline genes were validated using real-time PCR and are represented in Additional file 3.
Gene Ontology (GO) analysis with DAVID also revealed a level of similarity in terms of significant enrichment of radiation responsive gene ontology terms in WT and Il10 −/− mice (Fig. 4). We found innate immune response and p53 signaling were significantly enriched among radiation responsive genes in both WT and Il10 −/ − mice (Fig. 4a). However, additional immune-related GO terms, such as MHC-I antigen presentation, interferon signaling, and hepatic fibrosis, were also enriched among up-regulated genes of Il10 −/− mice but not in WT after radiation exposure. Also, down-regulated genes in WT, but not Il10 −/− , mice were significantly enriched with MHC-II protein complex terms, while DNA damage response and DNA repair terms were significantly enriched among down-regulated genes of Il10 −/− mice after radiation exposure (Fig. 4b).
We then used the IPA comparison analysis tool to assess and compare categories of biological processes activated or inhibited in the two genotypes after radiation exposure, as predicted from gene expression changes. We used default z-score ≥ 2 for activation or z-score ≤ − 2 for inhibition as the significance cutoff for our analyses (Fig. 5). We found canonical pathways involved in immune response to infection and B-cell signaling showing similar trends of activity in the two genotypes after radiation exposure. However, the predicted activation status of T-cell signaling differed between the two genotypes after radiation exposure. The T-cell signaling pathway was predicted to be activated in Il10 −/− mice but inhibited in WT mice in response to radiation (Additional file 5). In addition, the Th-1 response, which is associated with inflammation, appeared to be activated in Il10 −/− mice but inhibited in WT mice. Hormone based signaling pathways also seemed to have different activity status after radiation exposure in the two genotypes (Fig. 5a).
We also performed a similar predictive analysis for upstream regulators of the gene expression response to radiation in the two strains and found specific regulators of radiation response and immune function showing similar activity trends in the two genotypes. Significant activation of TP53, LPS, STAT1, IRF3 and IKBKB was indicated in the gene expression response of both mouse types. However, BNIP3L, IL27 and BCL2L1 showed opposite trends of activity in the two genotypes (Fig. 5b). We also identified a set of pathways and regulators that appeared significant in one genotype but not the other. These may suggest specific differences in the radiation response of WT and Il10 −/− mice.

Impact of inflammatory response on prediction of radiation exposure
Since one of our broad objectives is the development of gene expression signatures predictive of radiation exposure, we used microarray data to build and test classifiers using a greedy pairs method with 7 different algorithms.  To test the ability of a gene set selected from WT data to identify radiation exposure in the presence of chronic inflammation, a set of 24 genes was identified using a greedy pairs algorithm and the WT samples as a training data set ( Table 2). The data from the Il10 −/− mice were then used as a test set, and all 7 algorithms predicted the exposure status correctly with 100% accuracy (P-value 0.05). Many of the genes selected in the classifier were well-established radiation responsive genes with roles in the Trp53 pathway or in inflammation and immune response. In fact, "Immunity" was the only GO category significantly overrepresented (Benjamini corrected p-value = 0.025 by DAVID analysis) among the selected signature genes. As genes with immune functions formed such a significant part of the signature, we also used the Il10 −/− data to train a classifier (Table 3) and tested it on the WT data to see what impact this would have. The new classifier included a higher proportion of immunity related genes ("Immunity" overrepresented with a Benjamini corrected p-value = 0.0019 by DAVID analysis), and performed with 100% accuracy on the Il10 −/− samples during cross validation. Despite the increased contribution from immunity related genes, all 7 algorithms also predicted the exposure status of the WT test samples correctly with 100% accuracy. Finally, we validated expression levels of six of the classifier genes using quantitative real-time PCR. We chose four Trp53-responsive genes, Cdkn1a, Aen, Ccng1 and Phlda3 and two inflammatory genes, Ifit1 and Ifit3, for this analysis. As observed in the microarray results, we found a greater relative fold change of Tp53-responsive genes in WT compared to Il10 −/− mice after radiation exposure, and a greater relative fold change of inflammatory cytokines in Il10 −/− mice (Fig. 6). The inclusion of a 2 Gy dose measurement of Cdkn1a, Phlda3, Ifitm1 and Ifit3 (Additional file 6) suggests that the relative responses of the Trp53 regulated genes in the two strains are similar at low and high doses. In contrast, the heightened response of inflammatory genes in the Il10 −/− mice, although still significant for Ifitm1 at 2 Gy, was much greater at the higher dose. A more detailed delineation of the interaction between radiation dose and inflammatory gene expression in these two genotypes would be of interest.

Discussion
Exposure to ionizing radiation has long been known to cause production of various immune molecules like cytokines, DAMPs (Damage-associated molecular patterns) and cellular factors in the irradiated tissue, and to result in inflammation [5,6]. The degree of radiation mediated inflammation varies in individuals and depends on various factors including genetic, life-style and pre-existing disease conditions [14]. Inflammation has been indicated as a potential confounding factor in the development of IR biomarkers in epidemiological studies [31]. Previous transcriptomic studies have shown that immune and inflammation-related genes dominate the gene expression responses seen after in vivo radiation exposure [17,31] and therefore, inflammation related genes have been proposed as biomarkers of radiation exposure or to be used in biodosimetric panels [32]. This raises the concern that inflammation may affect radiation responses, and may impact the use of some genes in biodosimetric panels for use in triage after an accidental or deliberate radiological event. Several investigators have used treatment with LPS, a general component of the bacterial cell wall that is a potent inducer of inflammatory response [30,31], as a model for acute inflammation or infection. LPS has been shown to induce expression of some well-established radiation responsive genes and pathways [33][34][35], furthering concern about the impact of inflammatory states. As an initial approach to examining the potential impact of a pre-existing chronic inflammatory condition on radiation biodosimetry, we have used a well-established mouse model (Il10 −/− ) of human chronic inflammatory bowel disease (IBD) and compared its gene expression with and without radiation exposure with that of wild type (WT) littermates. Upstream regulators responding to radiation in WT and Il10 −/− mice. IPA based statistical z-score was used to predict activation (z ≥ 2, orange) or suppression (z ≤ − 2, blue) of a particular pathway or regulator. Grey color indicates no significant activity First, we compared blood counts in the two genotypes before and after irradiation. Of note was the significant difference in the baseline level of neutrophils in Il10 −/− compared to WT mice. Neutrophils are generally associated with inflammatory response. In addition, a higher level of neutrophils has been reported in IBD patients, which correlates with the severity of their clinical manifestations [36,37]. The higher level of neutrophils in Il10 −/− mice is similar to that in individuals suffering from IBD. The increased proportions of neutrophils contribute to a pro-inflammatory blood immune environment, which creates a concern of possible false positives or interference with radiation response for dosimetric biomarkers.
We next compared the baseline gene expression levels in un-irradiated Il10 −/− mice with those in WT mice (Fig. 2). The genes that showed significant differences in expression were mainly involved in inflammatory immune response, as might be expected, as well as glutathione metabolism and collagen deposition. In a recent publication, Mak et al. [38] showed that glutathione metabolism plays an important role in the activation process of inflammatory T-cells in response to immune stimuli. Excessive deposition of collagen, a major constituent of the extracellular matrix, has also been associated with inflammation and the development of fibrosis [39], functions also associated with radiation exposure. We found LPS / IL-1 signaling, a major regulator of gene expression responses in IBD patients [40], to be a significant pathway connecting the processes of inflammation, glutathione metabolism and fibrosis (Fig. 2c). These results suggest that chronic inflammatory conditions can indeed trigger some of the same gene expression patterns indicative of radiation response in the absence of exposure. Transcriptomic analysis of radiation responses in the two genotypes revealed broad similarity among the gene ontology terms most significantly affected by radiation exposure. Unsurprisingly, these included p53 signaling and some immune related functions (Fig. 4). Although significantly enriched among radiation responsive genes in both genotypes, the canonical radiation response of p53 signaling did appear to be slightly muted in the Il10 −/− mice, as implied by a slightly higher P value compared to WT-mice. This suggests that fewer p53-regulated genes responded significantly to radiation in the mutant mice. In confirming the radiation induction of individual genes, we found that both Cdkn1a and Aen were up regulated by radiation exposure to a lesser degree in the Il10 −/− mice, compared to WT mice (Fig. 6). There is some evidence in the literature of a regulatory link between IL10 and p53 in inflammatory processes. IL10-induced apoptosis in mouse bone marrow-derived mast cells has been shown to depend on activation of p53 [41]. Levels of IL10 and p53 were also found to be correlated in patients with autoimmune disease, with high levels of IL10 apparently leading to overexpression of p53 [42]. Taken together, these observations suggest that IL10 may contribute to the regulation of p53, leading to a blunted radiation response of genes downstream of p53 in the absence of IL10. The observed enrichment of genes with functions related to B-cell development among the down-regulated genes of both genotypes could be explained by the observed overall decrease in lymphocyte levels after radiation exposure (Fig. 1), particularly as B cells are among the most radiosensitive hematopoietic cells [43].
We observed that specific immune functions like MHC based antigen presentation and interferon signaling were particularly associated with radiation responsive genes in Il10 −/− mice but not WT mice. Antigen presentation mainly occurs through major histocompatibility complex (MHC) receptors on the cell surface, and is a very important process in activation of T-cells, Fig. 6 Gene expression measured by qRT-PCR. Expression levels of 6 genes (Cdkn1a, Aen, Phlda3, Ccng1, Ifitm1, Ifit3) showing significant difference in expression in Il10 −/− mice compared to WT after radiation exposure are presented here. The expression level for each gene was calculated relative to the expression levels in unexposed WT animals considered as 1. The expression was normalized against Beta-Actin gene expression. The data is represented as mean +/− SEM (n = 5). * Significant at P < 0.05, ** Significant at P < 0.01, *** Significant at P < 0.001, using unpaired t-test which mount an immune response by releasing cytokines or protein mediators of apoptosis [44]. The MHC-I is present on the cell surface of all the nucleated cells while MHC-II is present only in specific antigen presenting cells (APC), like macrophages and dendritic cells. Antigen presentation through MHC-II receptors results in activation of T-helper (Th) cells, which modulate the immune environment by releasing specific cytokines [44,45]. We found genes involved in MHC-II based antigen presentation to be specifically downregulated by radiation exposure in WT mice but upregulated in Il10 −/− mice. This indicates a differential modulation of the antigen presentation mechanism and its downstream immune response in the two genotypes after radiation exposure.
It has been reported that similar to IBD patients, Il10 −/− mice have higher levels IL23 producing macrophages, which through their interaction with T-cells via MHC-II receptors cause production of various inflammatory cytokines [25]. Intriguingly, IPA predicted a significant response to radiation of T-cell signaling in both genotypes, but in opposite directions. In Il10 −/− mice the gene expression pattern suggested a significant radiation activation of T-cell response, while in WT mice this pathway was predicted to be inhibited. Similarly, the T-helper type-1 (Th-1) immune response pathway was predicted to follow the same response pattern, with the pathway appearing to be activated by radiation exposure in Il10 −/− mice and inhibited in WT mice. Th1 mediated signaling has been reported to be the major causes of inflammation in people suffering from IBD as well as in animal models of IBD [24,40,45]. Production of interferons like IFN-gamma is a hallmark of Th1 immune response, and is also required for driving T-helper cell populations towards the Th1 type of immune response [46]. We found significant upregulation of interferon and its associated genes in Il10 −/− mice but not in WT after radiation exposure. This type of immune response also plays an important role against viral and bacterial infection. We found a significant association of these specific gene ontology terms specifically with the radiation response in Il10 −/− mice ( Fig. 4; Additional file 5). Overall, we found that radiation exposure results in higher inflammatory responses in Il10 −/− mice compared to WT, and this appears to be mediated by the T-cell population.

Impact on biodosimetry markers
Since our laboratory is interested in the development of biodosimetry markers for radiation exposure, we also used our data to select and test panels of genes that could distinguish radiation exposure in individual samples regardless of pre-existing inflammatory disease. The panels included p53 responsive genes like Cdkn1a and Aen as well as immune responsive genes like Ifit1 and Ifit3. Cdkn1a expression is known to respond widely to a number of stresses, and has been previously shown to respond to the inflammation inducing agent LPS [32]. In our current study, we found a significant difference between WT and Il10 −/− expression levels of Cdkn1a and other p53 regulated genes after radiation exposure, with the Cdkn1a radiation response being significantly muted in Il10 −/− mice compared to WT. The reverse was true of several immune related genes, which showed a larger magnitude of response to radiation exposure in the Il10 −/− mice compared to WT mice. Although in our current study these differences in response did not limit our ability to detect a relatively damaging dose of radiation, these results do suggest that excessive reliance on genes in these pathways may have an impact on the ultimate goal of dose discrimination. Care must also be taken to ensure that the elevated expression of inflammatory genes noted in the unirradiated Il10 −/− mice does not lead to misidentification of such inflammation-prone individuals as having received a small dose of radiation. Future studies across a range of doses will be needed to fully understand the impact of preexisting inflammation states on radiation biodosimetry using gene expression.
In a previous study, we had shown that genes selected from WT mice performed poorly in detecting exposure to an LD 50/30 radiation dose in two mouse models of DNA repair deficiency (Atm −/− and Prkdc scid ), but that inclusion of both WT and DNA repair deficient mice in the training set improved performance on both mutant and WT test sets to 100% [47]. We tested the ability of this earlier signature to predict radiation exposure in Il10 −/− mice. Although only one gene, Cdkn1a, was common between the earlier signature and either of the signatures reported in this study, we tested the WT plus repair mutant gene set on the dataset from this study and found that it was again able to predict the radiation exposure in Il10 −/− mice, and the new set of WT mice, with 100% efficiency (p < 0.005). This supports inclusion of individuals with a range of radiation sensitivities and underlying disease states in gene selection to produce more robust signatures for radiation biodosimetry.

Conclusions
Using a mouse model of inflammatory bowel disease, we have found that a chronic inflammatory state can impact both the basal expression levels and the relative radiation response of genes involved in immune and inflammatory pathways, as well as other key radiation response pathways, such as p53. While this did not impact the ability of gene expression signatures derived from either WT or pro-inflammatory Il10 −/− mice to correctly distinguish between unexposed and irradiated mice, it does raise concerns for the accurate determination of doses