Applicability of Gene Expression in Saliva as an Alternative to Blood for Biodosimetry and Prediction of Radiation-induced Health Effects.

As the great majority of gene expression (GE) biodosimetry studies have been performed using blood as the preferred source of tissue, searching for simple and less-invasive sampling methods is important when considering biodosimetry approaches. Knowing that whole saliva contains an ultrafiltrate of blood and white blood cells, it is expected that the findings in blood can also be found in saliva. This human in vivo study aims to examine radiation-induced GE changes in saliva for biodosimetry purposes and to predict radiation-induced disease, which is yet poorly characterized. Furthermore, we examined whether transcriptional biomarkers in blood can also be found equivalently in saliva. Saliva and blood samples were collected in parallel from radiotherapy (RT) treated patients who suffered from head and neck cancer (n = 8) undergoing fractioned partial-body irradiations (1.8 Gy/fraction and 50-70 Gy total dose). Samples were taken 12-24 h before first irradiation and ideally 24 and 48 h, as well as 5 weeks after radiotherapy onset. Due to the low quality and quantity of isolated RNA samples from one patient, they had to be excluded from further analysis, leaving a total of 24 saliva and 24 blood samples from 7 patients eligible for analysis. Using qRT-PCR, 18S rRNA and 16S rRNA (the ratio being a surrogate for the relative human RNA/bacterial burden), four housekeeping genes and nine mRNAs previously identified as radiation responsive in blood-based studies were detected. Significant GE associations with absorbed dose were found for five genes and after the 2nd radiotherapy fraction, shown by, e.g., the increase of CDKN1A (2.0 fold, P = 0.017) and FDXR (1.9 fold increased, P = 0.002). After the 25th radiotherapy fraction, however, all four genes (FDXR, DDB2, POU2AF1, WNT3) predicting ARS (acute radiation syndrome) severity, as well as further genes (including CCNG1 [median-fold change (FC) = 0.3, P = 0.013], and GADD45A (median-FC = 0.3, P = 0.031)) appeared significantly downregulated (FC = 0.3, P = 0.01-0.03). A significant association of CCNG1, POU2AF1, HPRT1, and WNT3 (P = 0.006-0.04) with acute or late radiotoxicity could be shown before the onset of these clinical outcomes. In an established set of four genes predicting acute health effects in blood, the response in saliva samples was similar to the expected up- (FDXR, DDB2) or downregulation (POU2AF1, WNT3) in blood for up to 71% of the measurements. Comparing GE responses (PHPT1, CCNG1, CDKN1A, GADD45A, SESN1) in saliva and blood samples, there was a significant linear association between saliva and blood response of CDKN1A (R2 = 0.60, P = 0.0004). However, the GE pattern of other genes differed between saliva and blood. In summary, the current human in vivo study, (I) reveals significant radiation-induced GE associations of five transcriptional biomarkers in salivary samples, (II) suggests genes predicting diverse clinical outcomes such as acute and late radiotoxicity as well as ARS severity, and (III) supports the view that blood-based GE response can be reflected in saliva samples, indicating that saliva is a "mirror of the body" for certain but not all genes and, thus, studies for each gene of interest in blood are required for saliva.

PUM1 and PGK1 are Favorable Housekeeping Genes over Established Biodosimetry-related Housekeeping Genes such as HPRT1, ITFG1, DPM1, MRPS5, 18S rRNA and Others after Radiation Exposure.

In gene expression (GE) studies, housekeeping genes (HKGs) are required for normalization purposes. In large-scale inter-laboratory comparison studies, significant differences in dose estimates are reported and divergent HKGs are employed by the teams. Among them, the 18S rRNA HKG is known for its robustness. However, the high abundance of 18S rRNA copy numbers requires dilution, which is time-consuming and a possible source of errors. This study was conducted to identify the most promising HKGs showing the least radiation-induced GE variance after radiation exposure. In the screening stage of this study, 35 HKGs were analyzed. This included selected HKGs (ITFG1, MRPS5, and DPM1) used in large-scale biodosimetry studies which were not covered on an additionally employed pre-designed 96-well platform comprising another 32 HKGs used for different exposures. Altogether 41 samples were examined, including 27 ex vivo X-ray irradiated blood samples (0, 0.5, 4 Gy), six X-irradiated samples (0, 0.5, 5 Gy) from two cell lines (U118, A549), as well as eight non-irradiated tissue samples to encompass multiple biological entities. In the independent validation stage, the most suitable candidate genes were examined from another 257 blood samples, taking advantage of already stored material originating from three studies. These comprise 100 blood samples from ex vivo X-ray irradiated (0-4 Gy) healthy donors, 68 blood samples from 5.8 Gy irradiated (cobalt-60) Rhesus macaques (RM) (LD29/60) collected 0-60 days postirradiation, and 89 blood samples from chemotherapy-(CTx) treated breast tumor patients. CTx and radiation-induced GE changes in previous studies appeared comparable. RNA was isolated, converted into cDNA, and GE was quantified employing TaqMan assays and quantitative RT-PCR. We calculated the standard deviation (SD) and the interquartile range (IQR) as measures of GE variance using raw cycle threshold (Ct) values and ranked the HKGs accordingly. Dose, time, age, and sex-dependent GE changes were examined employing the parametrical t-test and non-parametrical Kruskal Wallis test, as well as linear regression analysis. Generally, similar ranking results evolved using either SD or IQR GE measures of variance, indicating a tight distribution of GE values. PUM1 and PGK1 showed the lowest variance among the first ten most suitable genes in the screening phase. MRPL19 revealed low variance among the first ten most suitable genes in the screening phase only for blood and cells, but certain comparisons indicated a weak association of MRPL19 with dose (P = 0.02-0.09). In the validation phase, these results could be confirmed. Here, IQR Ct values from, e.g., X-irradiated blood samples were 0.6 raw Ct values for PUM1 and PGK1, which is considered to represent GE differences as expected due to methodological variance. Overall, when compared, the GE variance of both genes was either comparable or lower compared to 18S rRNA. Compared with the IQR GE values of PUM1 and PGKI, twofold-fivefold increased values were calculated for the biodosimetry HKG HPRT1, and comparable values were calculated for biodosimetry HKGs ITFG1, MRPS5, and DPM1. Significant dose-dependent associations were found for ITFG1 and MRPS5 (P = 0.001-0.07) and widely absent or weak (P = 0.02-0.07) for HPRT1 and DPM1. In summary, PUM1 and PGK1 appeared most promising for radiation exposure studies among the 35 HKGs examined, considering GE variance and adverse associations of GE with dose.

Simulations of radioiodine exposure and protective thyroid blocking in a new biokinetic model of the mother-fetus unit at different pregnancy ages.

In the case of nuclear incidents, radioiodine may be released. After incorporation, it accumulates in the thyroid and enhances the risk of thyroidal dysfunctions and cancer occurrence by internal irradiation. Pregnant women and children are particularly vulnerable. Therefore, thyroidal protection by administering a large dose of stable (non-radioactive) iodine, blocking radioiodide uptake into the gland, is essential in these subpopulations. However, a quantitative estimation of the protection conferred to the maternal and fetal thyroids in the different stages of pregnancy is difficult. We departed from an established biokinetic model for radioiodine in pregnancy using first-order kinetics. As the uptake of iodide into the thyroid and several other tissues is mediated by a saturable active transport, we integrated an uptake mechanism described by a Michaelis-Menten kinetic. This permits simulating the competition between stable and radioactive iodide at the membrane carrier site, one of the protective mechanisms. The Wollf-Chaikoff effect, as the other protective mechanism, was simulated by adding a total net uptake block for iodide into the thyroid, becoming active when the gland is saturated with iodine. The model's validity was confirmed by comparing predicted values with results from other models and sparse empirical data. According to our model, in the case of radioiodine exposure without thyroid blocking, the thyroid equivalent dose in the maternal gland increases about 45% within the first weeks of pregnancy to remain in the same range until term. Beginning in the 12th pregnancy week, the equivalent dose in the fetal thyroid disproportionately increases over time and amounts to three times the dose of the maternal gland at term. The maternal and fetal glands' protection increases concomitantly with the amount of stable iodine administered to the mother simultaneously with acute radioiodine exposure. The dose-effect curves reflecting the combined thyroidal protection by the competition at the membrane carrier site and the Wolff-Chaikoff effect in the mother are characterized by a mean effective dose (ED50) of roughly 1.5 mg all over pregnancy. In the case of the fetal thyroid, the mean effective doses for thyroid blocking, taking into account only the competition at the carrier site are numerically lower than in the mother. Taking into account additionally the Wolff-Chaikoff effect, the dose-effect curves for thyroidal protection in the fetus show a shift to the left over time, with a mean effective dose of 12.9 mg in the 12th week of pregnancy decreasing to 0.5 mg at term. In any case, according to our model, the usually recommended dose of 100 mg stable iodine given at the time of acute radioiodine exposure confers a very high level of thyroidal protection to the maternal and fetal glands over pregnancy. For ethical reasons, the possibilities of experimental studies on thyroid blocking in pregnant women are extremely limited. Furthermore, results from animal studies are associated with the uncertainties related to the translation of the data to humans. Thus model-based simulations may be a valuable tool for better insight into the efficacy of thyroidal protection and improve preparedness planning for uncommon nuclear or radiological emergencies.

Examining potential confounding factors in gene expression analysis of human saliva and identifying potential housekeeping genes.

Isolation of RNA from whole saliva, a non-invasive and easily accessible biofluid that is an attractive alternative to blood for high-throughput biodosimetry of radiological/nuclear victims might be of clinical significance for prediction and diagnosis of disease. In a previous analysis of 12 human samples we identified two challenges to measuring gene expression from total RNA: (1) the fraction of human RNA in whole saliva was low and (2) the bacterial contamination was overwhelming. To overcome these challenges, we performed selective cDNA synthesis for human RNA species only by employing poly(A)+-tail primers followed by qRT-PCR. In the current study, this approach was independently validated on 91 samples from 61 healthy donors. Additionally, we used the ratio of human to bacterial RNA to adjust the input RNA to include equal amounts of human RNA across all samples before cDNA synthesis, which then ensured comparable analysis using the same base human input material. Furthermore, we examined relative levels of ten known housekeeping genes, and assessed inter- and intra-individual differences in 61 salivary RNA isolates, while considering effects of demographical factors (e.g. sex, age), epidemiological factors comprising social habits (e.g. alcohol, cigarette consumption), oral hygiene (e.g. flossing, mouthwash), previous radiological diagnostic procedures (e.g. number of CT-scans) and saliva collection time (circadian periodic). Total human RNA amounts appeared significantly associated with age only (P ≤ 0.02). None of the chosen housekeeping genes showed significant circadian periodicity and either did not associate or were weakly associated with the 24 confounders examined, with one exception, 60% of genes were altered by mouthwash. ATP6, ACTB and B2M represented genes with the highest mean baseline expression (Ct-values ≤ 30) and were detected in all samples. Combining these housekeeping genes for normalization purposes did not decrease inter-individual variance, but increased the robustness. In summary, our work addresses critical confounders and provides important information for the successful examination of gene expression in human whole saliva.

A comparison of thyroidal protection by iodine and perchlorate against radioiodine exposure in Caucasians and Japanese.

Radioactive iodine released in nuclear accidents may accumulate in the thyroid and by irradiation enhances the risk of cancer. Radioiodine uptake into the gland can be inhibited by large doses of stable iodine or perchlorate. Nutritional iodine daily intake may impact thyroid physiology, so that radiological doses absorbed by the thyroid as well as thyroid blocking efficacy may differ in Japanese with a very rich iodine diet compared to Caucasians. Based on established biokinetic-dosimetric models for the thyroid, we derived the parameters for Caucasians and Japanese to quantitatively compare the effects of radioiodine exposure and the protective efficacy of thyroid blocking by stable iodine at the officially recommended dosages (100 mg in Germany, 76 mg in Japan) or perchlorate. The maximum transport capacity for iodine uptake into the thyroid is lower in Japanese compared to Caucasians. For the same radioiodine exposure pattern, the radiological equivalent thyroid dose is substantially lower in Japanese in the absence of thyroid blocking treatments. In the case of acute radioiodine exposure, stable iodine is less potent in Japanese (ED50 = 41.6 mg) than in Caucasians (ED50 = 2.7 mg) and confers less thyroid protection at the recommended dosages because of a delayed responsiveness to iodine saturation of the gland (Wolff-Chaikoff effect). Perchlorate (ED50 = 10 mg in Caucasians) at a dose of 1000 mg has roughly the same thyroid blocking effect as 100 mg iodine in Caucasians, whereas it confers a much better protection than 76 mg iodine in Japanese. For prolonged exposures, a single dose of iodine offer substantially lower protection than after acute radioiodine exposure in both groups. Repetitive daily iodine administrations improve efficacy without reaching levels after acute radioiodine exposure and achieve only slightly better protection in Japanese than in Caucasians. However, in the case of continuous radioiodine exposure, daily doses of 1000 mg perchlorate achieve a high protective efficacy in Caucasians as well as Japanese (> 0.98). In Caucasians, iodine (100 mg) and perchlorate (1000 mg) at the recommended dosages seem alternatives in case of acute radioiodine exposure, whereas perchlorate has a higher protective efficacy in the case of longer lasting radioiodine exposures. In Japanese, considering protective efficacy, preference should be given to perchlorate in acute as well as prolonged radioiodine exposure scenarios.

Inter-laboratory comparison of gene expression biodosimetry for protracted radiation exposures as part of the RENEB and EURADOS WG10 2019 exercise.

Large-scale radiation emergency scenarios involving protracted low dose rate radiation exposure (e.g. a hidden radioactive source in a train) necessitate the development of high throughput methods for providing rapid individual dose estimates. During the RENEB (Running the European Network of Biodosimetry) 2019 exercise, four EDTA-blood samples were exposed to an Iridium-192 source (1.36 TBq, Tech-Ops 880 Sentinal) at varying distances and geometries. This resulted in protracted doses ranging between 0.2 and 2.4 Gy using dose rates of 1.5-40 mGy/min and exposure times of 1 or 2.5 h. Blood samples were exposed in thermo bottles that maintained temperatures between 39 and 27.7 °C. After exposure, EDTA-blood samples were transferred into PAXGene tubes to preserve RNA. RNA was isolated in one laboratory and aliquots of four blinded RNA were sent to another five teams for dose estimation based on gene expression changes. Using an X-ray machine, samples for two calibration curves (first: constant dose rate of 8.3 mGy/min and 0.5-8 h varying exposure times; second: varying dose rates of 0.5-8.3 mGy/min and 4 h exposure time) were generated for distribution. Assays were run in each laboratory according to locally established protocols using either a microarray platform (one team) or quantitative real-time PCR (qRT-PCR, five teams). The qRT-PCR measurements were highly reproducible with coefficient of variation below 15% in ≥ 75% of measurements resulting in reported dose estimates ranging between 0 and 0.5 Gy in all samples and in all laboratories. Up to twofold reductions in RNA copy numbers per degree Celsius relative to 37 °C were observed. However, when irradiating independent samples equivalent to the blinded samples but increasing the combined exposure and incubation time to 4 h at 37 °C, expected gene expression changes corresponding to the absorbed doses were observed. Clearly, time and an optimal temperature of 37 °C must be allowed for the biological response to manifest as gene expression changes prior to running the gene expression assay. In conclusion, dose reconstructions based on gene expression measurements are highly reproducible across different techniques, protocols and laboratories. Even a radiation dose of 0.25 Gy protracted over 4 h (1 mGy/min) can be identified. These results demonstrate the importance of the incubation conditions and time span between radiation exposure and measurements of gene expression changes when using this method in a field exercise or real emergency situation.

The 2019-2020 EURADOS WG10 and RENEB Field Test of Retrospective Dosimetry Methods in a Small-Scale Incident Involving Ionizing Radiation.

With the use of ionizing radiation comes the risk of accidents and malevolent misuse. When unplanned exposures occur, there are several methods which can be used to retrospectively reconstruct individual radiation exposures; biological methods include analysis of aberrations and damage of chromosomes and DNA, while physical methods rely on luminescence (TL/OSL) or EPR signals. To ensure the quality and dependability of these methods, they should be evaluated under realistic exposure conditions. In 2019, EURADOS Working Group 10 and RENEB organized a field test with the purpose of evaluating retrospective dosimetry methods as carried out in potential real-life exposure scenarios. A 1.36 TBq 192Ir source was used to irradiate anthropomorphic phantoms in different geometries at doses of several Gy in an outdoor open-air geometry. Materials intended for accident dosimetry (including mobile phones and blood) were placed on the phantoms together with reference dosimeters (LiF, NaCl, glass). The objective was to estimate radiation exposures received by individuals as measured using blood and fortuitous materials, and to evaluate these methods by comparing the estimated doses to reference measurements and Monte Carlo simulations. Herein we describe the overall planning, goals, execution and preliminary outcomes of the 2019 field test. Such field tests are essential for the development of new and existing methods. The outputs from this field test include useful experience in terms of planning and execution of future exercises, with respect to time management, radiation protection, and reference dosimetry to be considered to obtain relevant data for analysis.

Predicting the Radiation Sensitivity of Male and Female Rhesus Macaques Using Gene Expression.

Radiosensitivity differs in humans and likely among closely-related primates. Reasons for variation in radiosensitivity are not well known. We examined preirradiation gene expression in peripheral blood among male and female rhesus macaques which did or did not survive (up to 60 days) after whole-body irradiation with 700 cGy (LD66/60). RNA samples originated from a blinded randomized Good Laboratory Practice study in 142 irradiated rhesus macaques. Animals were untreated (placebo), or treated using recombinant human IL-12, G-CSF or combination of the two. We evaluated gene expression in a two-phase study design where phase I was a whole genome screen [next generation sequencing (NGS)] for mRNAs (RNA-seq) using five RNA samples from untreated male and female animals per group of survivor and non-survivor (total n = 20). Differential gene expression (DGE) was defined as a statistically significant and ≥2-fold up- or downregulation of mRNA species and was calculated between groups of survivors and non-survivors (reference) and by gender. Altogether 659 genes were identified, but the overlapping number of differentially expressed genes (DGE) observed in both genders was small (n = 36). Fifty-eight candidate mRNAs were chosen for independent validation in phase II using the remaining samples (n = 122) evaluated with qRT-PCR. Among the 58 candidates, 16 were of significance or borderline significance (t test) by DGE. Univariate and multivariate logistic regression analysis and receiver operating characteristic (ROC) curve analysis further refined and identified the most outstanding validated genes and gene combinations. For untreated male macaques, we identified EPX (P = 0.005, ROC=1.0), IGF2BP1 (P = 0.05, ROC=0.74) and the combination of EPX with SLC22A4 (P = 0.03, ROC=0.85) which appeared most predictive for the clinical outcome for treated and combined (untreated and treated) male macaque groups, respectively. For untreated, treated and both combined female macaque groups the same gene (MBOAT4, P = 0.0004, ROC = 0.81) was most predictive. Based on the probability function of the ROC curves, up to 74% of preirradiation RNA measurements predicted survival with a positive and negative predictive value ranging between 85-100% and associated odds ratios reflecting a 2-3-fold elevated risk for surviving per unit change (cycle threshold value) in gene expression. In conclusion, we identified gender-dependent genes and gene combinations in preirradiation blood samples for survival prediction after irradiation in rhesus macaques.

Differentiating Total- or Partial-Body Irradiation in Baboons Using mRNA Expression Patterns: A Proof of Concept.

To better predict clinical outcome after radiation exposure, it is very important to know the absorbed dose and body areas exposed. Previously we found that 22 miRNAs appeared to predict total- and partial-body irradiation (TBI and PBI, respectively) patterns and were suggestive of the percentage of the body exposed in a baboon model. Motivated by these results, we performed a similar analysis on the transcriptional level (mRNAs) using whole genome microarrays. From 17 irradiated baboons, blood samples were taken before, and at 1, 2, 7, 28 and 75-106 days postirradiation to an equivalent TBI dose of 2.5 or 5 Gy applied either to the total body or to different parts of the body such as the upper body (UBE) or left hemibody (LHB). We compared quantile normalized log2-transformed gene expression values with three exposure pattern comparisons, namely TBI vs. PBI, TBI vs. LHB and UBE vs. LHB using Kruskal-Wallis and logistic regression analysis for receiver-operator characteristic (ROC) calculation. We found several hundred significantly (P < 0.05) and ≥2-fold deregulated mRNAs per exposure pattern comparison with a peak of 163-860 mRNAs at day 28. Lower numbers on day 2 (60 mRNAs) and day 7 (91-162 mRNAs) were observed, with the lowest number of deregulated mRNAs at day 75-106 (22-58 mRNAs). The 14 most promising mRNAs (e.g., LTF, DEFA3, OLFM4) appeared 10.1-46.2-fold upregulated and the exposure groups were completely or almost completely discriminated (ROC between 0.8-1.0). Several of the mRNA gene expression changes were significantly associated with the percentage of the body exposed. The numbers of overlapping genes used for diagnosis on consecutive days postirradiation were mostly 0 or less than 10. Bioinformatic analysis confirmed that at each time point different biological processes predominated. Our results suggest mRNA changes over time may be used to retrospectively determine radiation exposure patterns as partial or total body. mRNA gene expression changes likely could be applied over a longer time frame (2-75 days postirradiation) than miRNA, but due to the transient gene expression changes a different set of candidate mRNAs appears to be required at each day after irradiation.

miRNA Expression Patterns Differ by Total- or Partial-Body Radiation Exposure in Baboons.

In a radiation exposure event, a likely scenario may include either total-body irradiation (TBI) or different partial-body irradiation (PBI) patterns. Knowledge of the exposure pattern is expected to improve prediction of clinical outcome. We examined miRNA species in 17 irradiated baboons receiving an upper-body, left hemibody or total-body irradiation of 2.5 or 5 Gy. Blood samples were taken before irradiation and at 1, 2, 7, 28 and 75-106 days after irradiation. Using a qRT-PCR platform for simultaneous detection of 667 miRNAs, we identified 55 miRNAs over all time points. Candidate miRNAs, such as miR-17, miR-128 or miR-15b, significantly discriminated TBI from different PBI exposure patterns, and 5-to-10-fold changes in gene expression were observed among the groups. A total of 22 miRNAs (including miR-17) revealed significant linear associations of gene expression changes with the percentage of the exposed body area (P < 0.0001). All these changes were primarily observed at day 7 postirradiation and almost no miRNAs were detected either before or after 7 days. A significant association in the reduction of lymphocyte counts in TBI compared to PBI animals corresponded with the number of miRNA candidates. This finding suggests that our target miRNAs predominantly originated from irradiated lymphocytes. In summary, gene expression changes in the peripheral blood provided indications of the exposure pattern and a suggestion of the percentage of the exposed body area.

A quantitative comparison of the chemo- and radiotoxicity of uranium at different enrichment grades.

OBJECTIVE: The radiotoxic effects of uranium are often in the focus of the public fears but the chemical toxic effects of uranium are reported to surpass radiation effects. As there is no uranium isotope that is not radioactive, it is not possible to study chemical effects fully independently from radiation effects. In order to quantitate and compare radio- and chemotoxicity, we determined the median lethal doses of uranium due to its chemical toxicity and calculated the absorbed radiological doses resulting from the ingestion or inhalation of corresponding amounts depending on the isotopic enrichment grade. Committed effective doses over 50 years are related to the stochastic health effects like cancer occurrence and can be converted to a loss of statistical life time (mean loss 0.4 day / mSv). The equivalent doses absorbed within a short time frame permits conclusion on the induction of deterministic effects (e.g. acute radiation sickness). METHOD: Simulations were based on the biokinetic models of the International Commission for Radioprotection and performed using Integrated Modules for Bioassay Analysis software. Results were compared with the doses given by the calculator of the WISE uranium project. The fractions of the total doses absorbed within different time periods were derived from the respective areas under the activity-time curves in the whole body. RESULTS: The distribution of the total dose on the organs and tissues depends on the invasion pathway and the solubility of the compound. In the case of inhalation, the absorption of the total dose is more protracted than after ingestion. The incorporation of depleted or natural uranium in lethal amounts due to nephrotoxicity does not lead to deterministic radiation effects and is associated with committed effective doses reaching at most about 200 mSv (proposed possible threshold for therapeutic interventions after accidental radionuclide incorporation). The inhalation of low enriched uranium leads to higher effective doses up to 690 mSv, but they are still insufficient to cause acute deterministic effects. Even highly enriched uranium seems not to induce radiation nephropathy, but deterministic effects on the hematopoetic system cannot be excluded in particularly sensitive patients. But the equivalent doses to the lungs associated with the inhalation of poorly soluble compounds of highly enriched uranium are in a range that may induce radiation pneumonitis. CONCLUSION: Our findings give clear evidence that for depleted and natural uranium chemical toxicity is much more marked than radiotoxicity. However, this conclusion must not be drawn for enriched and in particular highly enriched compounds that besides stochastic effects may even cause deterministic radiation effects.

Impact of Inter-Individual Variance in the Expression of a Radiation-Responsive Gene Panel Used for Triage.

In previous studies we determined a gene expression signature in baboons for predicting the severity of hematological acute radiation syndrome. We subsequently validated a set of eight of these genes in leukemia patients undergoing total-body irradiation. In the current study, we addressed the effect of intra-individual variability on the basal level of expression of those eight radiation-responsive genes identified previously, by examining baseline levels in 200 unexposed healthy human donors (122 males and 88 females with an average age of 46 years) using real-time PCR. In addition to the eight candidate genes ( DAGLA, WNT3, CD177, PLA2G16, WLS, POU2AF1, STAT4 and PRF1), we examined two more genes ( FDXR and DDB2) widely used in ex vivo whole blood experiments. Although significant sex- (seven genes) and age-dependent (two genes) differences in expression were found, the fold changes ranged only between 1.1-1.6. These were well within the twofold differences in gene expression generally considered to represent control values. Age and sex contributed less than 20-30% to the complete inter-individual variance, which is calculated as the fold change between the lowest (reference) and the highest Ct value minimum-maximum fold change (min-max FC). Min-max FCs ranging between 10-17 were observed for most genes; however, for three genes, min-max FCs of complete inter-individual variance were found to be 37.1 ( WNT3), 51.4 ( WLS) and 1,627.8 ( CD177). In addition, to determine whether discrimination between healthy and diseased baboons might be altered by replacing the published gene expression data of the 18 healthy baboons with that of the 200 healthy humans, we employed logistic regression analysis and calculated the area under the receiver operating characteristic (ROC) curve. The additional inter-individual variance of the human data set had either no impact or marginal impact on the ROC area, since up to 32-fold change gene expression differences between healthy and diseased baboons were observed.

Correlation of Radiation Dose Estimates by DIC with the METREPOL Hematological Classes of Disease Severity.

The degree of severity of hematologic acute radiation syndrome (HARS) may vary across the range of radiation doses, such that dose alone may be a less reliable predictor of clinical course. We sought to elucidate the relationship between absorbed dose and risk of clinically relevant HARS in humans. We used the database SEARCH (System for Evaluation and Archiving of Radiation Accidents based on Case Histories), which contains the histories of radiation accident victims. From 153 cases we extracted data on dose estimates using the dicentric assay to measure individual biological dosimetry. The data were analyzed according to the corresponding hematological response categories of clinical significance (H1-4). These categories are derived from the medical treatment protocols for radiation accident victims (METREPOL) and represent the clinical outcome of HARS based on severity categories ranging from 1-4. In addition, the category H0 represents a post-exposure hematological response that is within the normal range for nonexposed individuals. Age at exposure, gender and ethnicity were considered as potential confounders in unconditional cumulative logistic regression analysis. In most cases, victims were Caucasian (82.4%) and male (92.8%), who originated from either the Chernobyl (69.3%) or Goiânia (10.5%) accident, and nearly 60% were aged 20-40 years at time of exposure. All individuals were whole-body exposed (mean 3.8 Gy, stdev ±3.1), and single exposures were predominantly reported (79%). Seventy percent of victims in category H0 were exposed to ≤1 Gy, with rapidly decreasing proportions of H0 seen at doses up to 5 Gy. There were few HARS H4 cases reported at exposed dose of 1-2 Gy, while 82% of H4 cases received doses of >5 Gy. HARS H1-3 cases varied among dose ranges from 1-5 Gy. In summary, single whole-body radiation doses <1 Gy and >5 Gy corresponded in general with H0 and H3-4, respectively, and this was consistent with medical expectations. This underlines the usefulness of dose estimates for HARS prediction. However, whole-body doses between 1-5 Gy poorly corresponded to HARS H1-3. The dose range of 1-5 Gy was of limited value for medical decision-making regarding, e.g., hospitalization for H2-3, but not H1 and treatment decisions that differ between H1-3. Also, there were some H0 cases at high doses and H2-4 cases at low doses, thereby challenging an individual recommendation based solely on dose.

Validating Baboon Ex Vivo and In Vivo Radiation-Related Gene Expression with Corresponding Human Data.

The research for high-throughput diagnostic tests for victims of radio/nuclear incidents remains ongoing. In this context, we have previously identified candidate genes that predict risk of late-occurring hematologic acute radiation syndrome (HARS) in a baboon model. The goal of the current study was to validate these genes after radiation exposure in humans. We also examined ex vivo relative to in vivo measurements in both species and describe dose-response relationships. Eighteen baboons were irradiated in vivo to simulate different patterns of partial- or total-body irradiation (TBI), corresponding to an equivalent dose of 2.5 or 5 Sv. Human in vivo blood samples were obtained from patients exposed to different dose ranges: diagnostic computerized tomography (CT; 0.004-0.018 Sv); radiotherapy for prostate cancer (0.25-0.3 Sv); and TBI of leukemia patients (2 × 1.5 or 2 × 2 Sv, five patients each). Peripheral whole blood of another five baboons and human samples from five healthy donors were cultivated ex vivo and irradiated with 0-4 Sv. RNA was isolated pairwise before and 24 h after irradiation and converted into cDNA. Gene expression of six promising candidate genes found previously by us in a baboon model ( WNT3, POU2AF1, CCR7, ARG2, CD177, WLS), as well as three genes commonly used in ex vivo whole blood experiments ( FDXR, PCNA, DDB2) was measured using qRT-PCR. We confirmed the six baboon candidate genes in leukemia patients. However, expression for the candidate gene FDXR showed an inverse relationship, as it was downregulated in baboons and upregulated in human samples. Comparisons among the in vivo and ex vivo experiments revealed the same pattern in both species and indicated peripheral blood cells to represent the radiation-responsive targets causing WNT3 and POU2AF1 gene expression changes. CCR7, ARG2, CD177 and WLS appeared to be altered due to radiation-responsive targets other than the whole blood cells. Linear dose-response relationships of FDXR, WNT3 and POU2AF1 using human ex vivo samples corresponded with human in vivo samples, suggesting that ex vivo models for in vivo dose estimates can be used over a wide dose range (0.001-5 Sv for POU2AF1). In summary, we validated six baboon candidate genes in humans, but the FDXR measurements underscored the importance of independent assessments even when candidates from animal models have striking gene sequence homology to humans. Since whole blood cells represented the same radiation-responsive targets for FDXR, WNT3 and POU2AF1 gene expression changes, ex vivo cell culture models can be utilized for in vivo dose estimates over a dose range covering up to 3.5 log scales. These findings might be a step forward in the development of a gene expression-based high-throughput diagnostic test for populations involved in large-scale radio/nuclear incidents.

Rapid Prediction of Hematologic Acute Radiation Syndrome in Radiation Injury Patients Using Peripheral Blood Cell Counts.

Rapid clinical triage of radiation injury patients is essential for determining appropriate diagnostic and therapeutic interventions. We examined the utility of blood cell counts (BCCs) in the first three days postirradiation to predict clinical outcome, specifically for hematologic acute radiation syndrome (HARS). We analyzed BCC test samples from radiation accident victims (n = 135) along with their clinical outcome HARS severity scores (H1-4) using the System for Evaluation and Archiving of Radiation Accidents based on Case Histories (SEARCH) database. Data from nonirradiated individuals (H0, n = 132) were collected from an outpatient facility. We created binary categories for severity scores, i.e., 1 (H0 vs. H1-4), 2 (H0-1 vs. H2-4) and 3 (H0-2 vs. H3-4), to assess the discrimination ability of BCCs using unconditional logistic regression analysis. The test sample contained 454 BCCs from 267 individuals. We validated the discrimination ability on a second independent group comprised of 275 BCCs from 252 individuals originating from SEARCH (HARS 1-4), an outpatient facility (H0) and hospitals (e.g., leukemia patients, H4). Individuals with a score of H0 were easily separated from exposed individuals based on developing lymphopenia and granulocytosis. The separation of H0 and H1-4 became more prominent with increasing hematologic severity scores and time. On day 1, lymphocyte counts were most predictive for discriminating binary categories, followed by granulocytes and thrombocytes. For days 2 and 3, an almost complete separation was achieved when BCCs from different days were combined, supporting the measurement of sequential BCC. We found an almost complete discrimination of H0 vs. irradiated individuals during model validation (negative predictive value, NPV > 94%) for all three days, while the correct prediction of exposed individuals increased from day 1 (positive predictive value, PPV 78-89%) to day 3 (PPV > 90%). The models were unable to provide predictions for 10.9% of the test samples, because the PPVs or NPVs did not reach a 95% likelihood defined as the lower limit for a prediction. We developed a prediction model spreadsheet to provide early and prompt diagnostic predictions and therapeutic recommendations including identification of the worried well, requirement of hospitalization or development of severe hematopoietic syndrome. These results improve the provisional classification of HARS. For the final diagnosis, further procedures (sequential diagnosis, retrospective dosimetry, clinical follow-up, etc.) must be taken into account. Clinical outcome of radiation injury patients can be rapidly predicted within the first three days postirradiation using peripheral BCC.

Using Clinical Signs and Symptoms for Medical Management of Radiation Casualties - 2015 NATO Exercise.

The utility of early-phase (≤5 days) radiation-induced clinical signs and symptoms (e.g., vomiting, diarrhea, erythema and changes in blood cell counts) was examined for the prediction of later occurring acute radiation syndrome (ARS) severity and the development of medical management strategies. Medical treatment protocols for radiation accident victims (METREPOL) was used to grade ARS severities, which were assigned response categories (RCs). Data on individuals (n = 191) with mild (RC1, n = 45), moderate (RC2, n = 19), severe (RC3, n = 20) and fatal (RC4, n = 18) ARS, as well as nonexposed individuals (RC0, n = 89) were generated using either METREPOL (n = 167) or the system for evaluation and archiving of radiation accidents based on case histories (SEARCH) database (n = 24), the latter comprised of real-case descriptions. These data were converted into tables reflecting clinical signs and symptoms, and submitted to eight teams representing five participating countries. The teams were comprised of medical doctors, biologists and pharmacists with subject matter expertise. The tables comprised cumulated clinical data from day 1-3 and day 1-5 postirradiation. While it would have reflected a more realistic scenario to provide the data to the teams over the course of a 3- or 5-day period, the logistics of doing so proved too challenging. In addition, the team members participating in this exercise chose to receive the cumulated reports of day 1-3 and 1-5. The teams were tasked with predicting ARS incidence, ARS severity and the requirement for hospitalization for multiple cases, as well as providing the certainty of their diagnosis. Five of the teams also performed dose estimates. The teams did not employ harmonized methodologies, and the expertise among the members varied, as did the tools used and the means of analyzing the clinical data. The earliest report time was 3 h after the tables were sent to the team members. The majority of cases developing ARS (89.6% ± 3.3 SD) and requiring hospitalization (88.8% ± 4.6 SD) were correctly identified by all teams. Determination of ARS severity was particularly challenging for RC2-3, which was systematically overestimated. However, RC4 was correctly predicted at 94-100% by all teams. RC0 and RC1 ARS severities were more difficult to discriminate. When reported RCs (0-1 and 3-4) were merged, on average 89.6% (±3.3 SD) of all cases could be correctly classified. Comparisons on frequency distributions revealed no statistically significant differences among the following: 1. reported ARS from different teams (P > 0.2); 2. cases generated based on METREPOL or SEARCH (P > 0.5); or 3. results reported at day 3 and 5 postirradiation (P > 0.1). Dose estimates of all teams increased significantly along with ARS severity (P < 0.0001) as well as with dose estimates generated from dicentric chromosomal-aberration measurements available for SEARCH cases (P < 0.0001). In summary, early-phase radiation-induced clinical signs and symptoms proved to be useful for rapid and accurate assessment, with minor limitations, toward predicting life-threatening ARS severity and developing treatment management strategies.

The Impact of Time on Decorporation Efficacy After a "Dirty Bomb" Attack Studied by Simulation.

Background: In the case of a nuclear or radiological incident, there is a risk of external and internal contamination with radionuclides in addition to external irradiation. There is no consensus whether decorporation treatment should be initiated right away on spec or pending the results of internal dosimetry to determine the indication. Method: Based on biokinetic models for plutonium-239, americium-241 and cesium-137, the efficacy of a decorporation treatment using DTPA or Prussian blue was simulated depending on the initiation time and the duration of treatment for different invasion pathways and physicochemical properties of the inhaled compounds. Results: For the same level of radioactivity incorporated, the committed effective dose increases with the speed of the invasion process. The impact of the initiation time of a decorporation treatment is particularly important when the absorption of the radionuclide is fast. Even if started early after incorporation, the therapeutic efficacy is less for americium-241 or cesium-137 compared to plutonium-239. Therapeutic efficacy increases with treatment duration up to about 90 days for plutonium-239 and cesium-137, whereas a prolongation of the treatment over this limit may further enhance efficacy in the case of americium-241. Conclusion: In the case of a nuclear incident, several fractions with different but a priori unknown physicochemical properties may be inhaled. Thus, decorporation therapy should be started as soon as possible after the incorporation of the radionuclide(s), as a loss of efficacy caused by a delay of treatment initiation possibly cannot be compensated later on. Treatment should be pursued for several months.

Reconsidering Current Decorporation Strategies after Incorporation of Radionuclides.

In the case of a nuclear accident or a terrorist attack by a "dirty bomb," there is a risk of external and internal contamination with radionuclides in addition to external irradiation. Internal irradiation as a consequence of radionuclide incorporation is associated with a higher risk of stochastic radiation effects (e.g., tumors). Decorporation treatment will enhance the elimination of radionuclides and reduce the committed effective dose as a metric of stochastic health effects. Although treatment efficacy is better when started early, beginning the therapy without knowing the committed effective dose may unnecessarily expose the patient to the side effects of the medication. The question is: Delay the therapy to wait for the results of internal dosimetry or start the therapy promptly on spec? To prove insight into this question, a selective review of the literature was conducted. The importance of the initiation time of treatment in the efficacy of decorporation treatment can be explained with pharmacokinetic laws and first order processes determining the disposition of xenobiotics in the organism. Nevertheless, there is no internationally accepted standard on when to start a decorporation therapy (exception: iodide). The "precautionary approach," emphasizing the importance of the committed effective dose for the indication of treatment, is competing with the "urgent approach" advocating the administration of medication "a priori" within several hours. A review of the literature actually indicates that the most important drugs used for decorporation are well tolerated with few adverse effects. In consideration of the higher efficacy and the low side-effects of a short-term treatment, initiating decorporation therapy as soon as possible after internal contamination, even before the committed effective dose has been assessed, appears to be a reasonable approach. The decision of continuation or discontinuation of the therapy should be taken after internal dosimetry is completed on the basis of the committed effective dose.

First Generation Gene Expression Signature for Early Prediction of Late Occurring Hematological Acute Radiation Syndrome in Baboons.

We implemented a two-stage study to predict late occurring hematologic acute radiation syndrome (HARS) in a baboon model based on gene expression changes measured in peripheral blood within the first two days after irradiation. Eighteen baboons were irradiated to simulate different patterns of partial-body and total-body exposure, which corresponded to an equivalent dose of 2.5 or 5 Gy. According to changes in blood cell counts the surviving baboons (n = 17) exhibited mild (H1-2, n = 4) or more severe (H2-3, n = 13) HARS. Blood samples taken before irradiation served as unexposed control (H0, n = 17). For stage I of this study, a whole genome screen (mRNA microarrays) was performed using a portion of the samples (H0, n = 5; H1-2, n = 4; H2-3, n = 5). For stage II, using the remaining samples and the more sensitive methodology, qRT-PCR, validation was performed on candidate genes that were differentially up- or down-regulated during the first two days after irradiation. Differential gene expression was defined as significant (P < 0.05) and greater than or equal to a twofold difference above a H0 classification. From approximately 20,000 genes, on average 46% appeared to be expressed. On day 1 postirradiation for H2-3, approximately 2-3 times more genes appeared up-regulated (1,418 vs. 550) or down-regulated (1,603 vs. 735) compared to H1-2. This pattern became more pronounced at day 2 while the number of differentially expressed genes decreased. The specific genes showed an enrichment of biological processes coding for immune system processes, natural killer cell activation and immune response (P = 1 × E-06 up to 9 × E-14). Based on the P values, magnitude and sustained differential gene expression over time, we selected 89 candidate genes for validation using qRT-PCR. Ultimately, 22 genes were confirmed for identification of H1-3 classifications and seven genes for identification of H2-3 classifications using qRT-PCR. For H1-3 classifications, most genes were constantly three to fivefold down-regulated relative to H0 over both days, but some genes appeared 10.3-fold (VSIG4) or even 30.7-fold up-regulated (CD177) over H0. For H2-3, some genes appeared four to sevenfold up-regulated relative to H0 (RNASE3, DAGLA, ARG2), but other genes showed a strong 14- to 33-fold down-regulation relative to H0 (WNT3, POU2AF1, CCR7). All of these genes allowed an almost completely identifiable separation among each of the HARS categories. In summary, clinically relevant HARS can be independently predicted with all 29 irradiated genes examined in the peripheral blood of baboons within the first two days postirradiation. While further studies are needed to confirm these findings, this model shows potential relevance in the prediction of clinical outcomes in exposed humans and as an aid in the prioritizing of medical treatment.

Examining Radiation-Induced In Vivo and In Vitro Gene Expression Changes of the Peripheral Blood in Different Laboratories for Biodosimetry Purposes: First RENEB Gene Expression Study.

The risk of a large-scale event leading to acute radiation exposure necessitates the development of high-throughput methods for providing rapid individual dose estimates. Our work addresses three goals, which align with the directive of the European Union's Realizing the European Network of Biodosimetry project (EU-RENB): 1. To examine the suitability of different gene expression platforms for biodosimetry purposes; 2. To perform this examination using blood samples collected from prostate cancer patients (in vivo) and from healthy donors (in vitro); and 3. To compare radiation-induced gene expression changes of the in vivo with in vitro blood samples. For the in vitro part of this study, EDTA-treated whole blood was irradiated immediately after venipuncture using single X-ray doses (1 Gy/min(-1) dose rate, 100 keV). Blood samples used to generate calibration curves as well as 10 coded (blinded) samples (0-4 Gy dose range) were incubated for 24 h in vitro, lysed and shipped on wet ice. For the in vivo part of the study PAXgene tubes were used and peripheral blood (2.5 ml) was collected from prostate cancer patients before and 24 h after the first fractionated 2 Gy dose of localized radiotherapy to the pelvis [linear accelerator (LINAC), 580 MU/min, exposure 1-1.5 min]. Assays were run in each laboratory according to locally established protocols using either microarray platforms (2 laboratories) or qRT-PCR (2 laboratories). Report times on dose estimates were documented. The mean absolute difference of estimated doses relative to the true doses (Gy) were calculated. Doses were also merged into binary categories reflecting aspects of clinical/diagnostic relevance. For the in vitro part of the study, the earliest report time on dose estimates was 7 h for qRT-PCR and 35 h for microarrays. Methodological variance of gene expression measurements (CV ≤10% for technical replicates) and interindividual variance (≤twofold for all genes) were low. Dose estimates based on one gene, ferredoxin reductase (FDXR), using qRT-PCR were as precise as dose estimates based on multiple genes using microarrays, but the precision decreased at doses ≥2 Gy. Binary dose categories comprising, for example, unexposed compared with exposed samples, could be completely discriminated with most of our methods. Exposed prostate cancer blood samples (n = 4) could be completely discriminated from unexposed blood samples (n = 4, P < 0.03, two-sided Fisher's exact test) without individual controls. This could be performed by introducing an in vitro-to-in vivo correction factor of FDXR, which varied among the laboratories. After that the in vitro-constructed calibration curves could be used for dose estimation of the in vivo exposed prostate cancer blood samples within an accuracy window of ±0.5 Gy in both contributing qRT-PCR laboratories. In conclusion, early and precise dose estimates can be performed, in particular at doses ≤2 Gy in vitro. Blood samples of prostate cancer patients exposed to 0.09-0.017 Gy could be completely discriminated from pre-exposure blood samples with the doses successfully estimated using adjusted in vitro-constructed calibration curves.