Validating a Four-gene Set for H-ARS Severity Prediction in Peripheral Blood Samples of Irradiated Rhesus Macaques

Increased radiological and nuclear threats require preparedness. Our earlier work identified a set of four genes (DDB2, FDXR, POU2AF1 and WNT3), which predicts severity of the hematological acute radiation syndrome (H-ARS) within the first three days postirradiation In this study of 41 Rhesus macaques (Macaca mulatta, 27 males, 14 females) irradiated with 5.8-7.2 Gy (LD29-50/60), including some treated with gamma-tocotrienol (GT3, a radiation countermeasure) we independently validated these genes as predictors in both sexes and examined them after three days. At the Armed Forces Radiobiology Research Institute/Uniformed Services University of the Health Sciences, peripheral whole blood (1 ml) of Rhesus macaques was collected into PAXgene® Blood RNA tubes pre-irradiation after 1, 2, 3, 35 and 60 days postirradiation, stored at -80°C for internal experimental analyses. Leftover tubes from these already ongoing studies were kindly provided to Bundeswehr Institute of Radiobiology. RNA was isolated (QIAsymphony), converted into cDNA, and for further gene expression (GE) studies quantitative RT-PCR was performed. Differential gene expression (DGE) was measured relative to the pre-irradiation Rhesus macaques samples. Within the first three days postirradiation, we found similar results to human data: 1. FDXR and DDB2 were up-regulated, FDXR up to 3.5-fold, and DDB2 up to 13.5-fold in the median; 2. POU2AF1 appeared down regulated around tenfold in nearly all Rhesus macaques; 3. Contrary to human data, DDB2 was more up-regulated than FDXR, and the difference of the fold change (FC) ranged between 2.4 and 10, while the median fold changes of WNT3, except days 1 and 35, were close to 1. Nevertheless, 46% of the Rhesus macaques showed down-regulated WNT3 on day one postirradiation, which decreased to 12.2% on day 3 postirradiation. Considering the extended phase, there was a trend towards decreased fold changes at day 35, with median-fold changes ranging from 0.7 for DDB2 to 0.1 for POU2AF1, and on day 60 postirradiation, DGE in surviving animals was close to pre-exposure values for all four genes. In conclusion, the diagnostic significance for radiation-induced H-ARS severity prediction of FDXR, DDB2, and POU2AF1 was confirmed in this Rhesus macaques model. However, DDB2 showed higher GE values than FDXR. As shown in previous studies, the diagnostic significance of WNT3 could not be reproduced in Rhesus macaques; this could be due to the choice of animal model and methodological challenges.

Validating Radiosensitivity with Pre-Exposure Differential Gene Expression in Peripheral Blood Predicting Survival and Non-Survival in a Second Irradiated Rhesus Macaque Cohort

Radiosensitivity differs in humans and possibly in closely related nonhuman primates. The reasons for variation in radiosensitivity are not well known. In an earlier study, we examined gene expression (GE) pre-radiation in peripheral blood among male (n = 62) and female (n = 60) rhesus macaques (n = 122), which did or did not survive (up to 60 days) after whole-body exposure of 7.0 Gy (LD66/60). Eight genes (CHD5, CHI3L1, DYSF, EPX, IGF2BP1, LCN2, MBOAT4, SLC22A4) revealed significant associations with survival. Access to a second rhesus macaque cohort (males = 40, females = 23, total n = 63) irradiated with 5.8-7.2 Gy (LD29-50/60) and some treated with gamma-tocotrienol (GT3, a radiation countermeasure) allowed us to validate these gene expression changes independently. Total RNA was isolated from whole blood samples and examined by quantitative RT-PCR on a 96-well format. cycle threshold (Ct)-values normalized to 18S rRNA were analyzed for their association with survival. Regardless of the species-specific TaqMan assay, similar results were obtained. Two genes (CHD5 and CHI3L1) out of eight revealed a significant association with survival in the second cohort, while only CHD5 (involved in DNA damage response and proliferation control) showed mean gene expression changes in the same direction for both cohorts. No expected association of CHD5 GE with dose, treatment, or sex could be established. Instead, we observed significant associations for those comparisons comprising pre-exposure samples with CHD5 Ct values ≤ 11 (total n = 17). CHD5 Ct values ≤ 11 in these comparisons were mainly associated with increased frequencies (61-100%) of non-survivors, a trend which depending on the sample numbers, reached significance (P = 0.03) in males and, accordingly, in females. This was also reflected by a logistic regression model including all available samples from both cohorts comprising CHD5 measurements (n = 104, odds ratio 1.38, 95% CI 1.07-1.79, P = 0.01). However, this association was driven by males (odds ratio 1.62, 95% CI 1.10-2.38, P = 0.01) and CHD5 Ct values ≤ 11 since removing low CHD5 Ct values from this model, converted to insignificance (P = 0.19). A second male subcohort comprising high CHD5 Ct values ≥ 14.4 in both cohorts (n = 5) appeared associated with survival. Removing these high CHD5 Ct values converted the model borderline significant (P = 0.051). Based on the probability function of the receiver operating characteristics (ROC) curves, 8 (12.3%) and 5 (7.7%) from 65 pre-exposure RNA measurements in males, death and survival could be predicted with a negative and positive predictive value ranging between 85-100%. An associated odds ratio reflected a 62% elevated risk for dying or surviving per unit change (Ct-value) in gene expression, considering the before-mentioned CHD5 thresholds in RNA copy numbers. In conclusion, we identified two subsets of male animals characterized by increased (Ct values ≤ 11) and decreased (Ct values ≥ 14.4) CHD5 GE copy numbers before radiation exposure, which independently of the cohort, radiation exposure or treatment appeared to predict the death or survival in males.

Development of a Point-of-Care Microfluidic RNA Extraction Slide for Gene Expression Diagnosis after Irradiation

In times of war, radiological/nuclear emergency scenarios have become a reemphasized threat. However, there are challenges in transferring whole-blood samples to laboratories for specialized diagnostics using RNA. This project aims to miniaturize the process of unwieldy conventional RNA extraction with its stationed technical equipment using a microfluidic-based slide (MBS) for point-of-care diagnostics. The MBS is thought to be a preliminary step toward the development of a so-called lab-on-a-chip microfluidic device. A MBS would enable early and fast field care combined with gene expression (GE) analysis for the prediction of hematologic acute radiation syndrome (HARS) severity or identification of RNA microbes. Whole blood samples from ten healthy donors were irradiated with 0, 0.5 and 4 Gy, simulating different ARS severity degrees. RNA quality and quantity of a preliminary MBS was compared with a conventional column-based (CB) RNA extraction method. GE of four HARS severity-predicting radiation-induced genes (FDXR, DDB2, POU2AF1 and WNT3) was examined employing qRT-PCR. Compared to the CB method, twice as much total RNA from whole blood could be extracted using the MBS (6.6 ± 3.2 µg vs. 12.0 ± 5.8 µg) in half of the extraction time, and all MBS RNA extracts appeared DNA-free in contrast to the CB method (30% were contaminated with DNA). Using MBS, RNA quality [RNA integrity number equivalent (RINe)] values decreased about threefold (3.3 ± 0.8 vs. 9.0 ± 0.4), indicating severe RNA degradation, while expected high-quality RINe ≥ 8 were found using column-based method. However, normalized cycle threshold (Ct) values, as well as radiation-induced GE fold-changes appeared comparable for all genes utilizing both methods, indicating that no RNA degradation took place. In summary, the preliminary MBS showed promising features such as: 1. halving the RNA extraction time without the burden of heavy technical equipment (e.g., a centrifuge); 2. absence of DNA contamination in contrast to CB RNA extraction; 3. reduction in blood required, because of twice the biological output of RNA; and 4. equal GE performance compared to CB, thus, increasing its appeal for later semi-automatic parallel field applications.

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.

RENEB Inter-Laboratory Comparison 2021: Inter-Assay Comparison of Eight Dosimetry Assays

Tools for radiation exposure reconstruction are required to support the medical management of radiation victims in radiological or nuclear incidents. Different biological and physical dosimetry assays can be used for various exposure scenarios to estimate the dose of ionizing radiation a person has absorbed. Regular validation of the techniques through inter-laboratory comparisons (ILC) is essential to guarantee high quality results. In the current RENEB inter-laboratory comparison, the performance quality of established cytogenetic assays [dicentric chromosome assay (DCA), cytokinesis-block micronucleus assay (CBMN), stable chromosomal translocation assay (FISH) and premature chromosome condensation assay (PCC)] was tested in comparison to molecular biological assays [gamma-H2AX foci (gH2AX), gene expression (GE)] and physical dosimetry-based assays [electron paramagnetic resonance (EPR), optically or thermally stimulated luminescence (LUM)]. Three blinded coded samples (e.g., blood, enamel or mobiles) were exposed to 0, 1.2 or 3.5 Gy X-ray reference doses (240 kVp, 1 Gy/min). These doses roughly correspond to clinically relevant groups of unexposed to low exposed (0-1 Gy), moderately exposed (1-2 Gy, no severe acute health effects expected) and highly exposed individuals (>2 Gy, requiring early intensive medical care). In the frame of the current RENEB inter-laboratory comparison, samples were sent to 86 specialized teams in 46 organizations from 27 nations for dose estimation and identification of three clinically relevant groups. The time for sending early crude reports and more precise reports was documented for each laboratory and assay where possible. The quality of dose estimates was analyzed with three different levels of granularity, 1. by calculating the frequency of correctly reported clinically relevant dose categories, 2. by determining the number of dose estimates within the uncertainty intervals recommended for triage dosimetry (±0.5 Gy or ±1.0 Gy for doses <2.5 Gy or >2.5 Gy), and 3. by calculating the absolute difference (AD) of estimated doses relative to the reference doses. In total, 554 dose estimates were submitted within the 6-week period given before the exercise was closed. For samples processed with the highest priority, earliest dose estimates/categories were reported within 5-10 h of receipt for GE, gH2AX, LUM, EPR, 2-3 days for DCA, CBMN and within 6-7 days for the FISH assay. For the unirradiated control sample, the categorization in the correct clinically relevant group (0-1 Gy) as well as the allocation to the triage uncertainty interval was, with the exception of a few outliers, successfully performed for all assays. For the 3.5 Gy sample the percentage of correct classifications to the clinically relevant group (≥2 Gy) was between 89-100% for all assays, with the exception of gH2AX. For the 1.2 Gy sample, an exact allocation to the clinically relevant group was more difficult and 0-50% or 0-48% of the estimates were wrongly classified into the lowest or highest dose categories, respectively. For the irradiated samples, the correct allocation to the triage uncertainty intervals varied considerably between assays for the 1.2 Gy (29-76%) and 3.5 Gy (17-100%) samples. While a systematic shift towards higher doses was observed for the cytogenetic-based assays, extreme outliers exceeding the reference doses 2-6 fold were observed for EPR, FISH and GE assays. These outliers were related to a particular material examined (tooth enamel for EPR assay, reported as kerma in enamel, but when converted into the proper quantity, i.e. to kerma in air, expected dose estimates could be recalculated in most cases), the level of experience of the teams (FISH) and methodological uncertainties (GE). This was the first RENEB ILC where everything, from blood sampling to irradiation and shipment of the samples, was organized and realized at the same institution, for several biological and physical retrospective dosimetry assays. Almost all assays appeared comparably applicable for the identification of unexposed and highly exposed individuals and the allocation of medical relevant groups, with the latter requiring medical support for the acute radiation scenario simulated in this exercise. However, extreme outliers or a systematic shift of dose estimates have been observed for some assays. Possible reasons will be discussed in the assay specific papers of this special issue. In summary, this ILC clearly demonstrates the need to conduct regular exercises to identify research needs, but also to identify technical problems and to optimize the design of future ILCs.

RENEB Inter-Laboratory Comparison 2021: The Gene Expression Assay

Early and high-throughput individual dose estimates are essential following large-scale radiation exposure events. In the context of the Running the European Network for Biodosimetry and Physical Dosimetry (RENEB) 2021 exercise, gene expression assays were conducted and their corresponding performance for dose-assessment is presented in this publication. Three blinded, coded whole blood samples from healthy donors were exposed to 0, 1.2 and 3.5 Gy X-ray doses (240 kVp, 1 Gy/min) using the X-ray source Yxlon. These exposures correspond to clinically relevant groups of unexposed, low dose (no severe acute health effects expected) and high dose exposed individuals (requiring early intensive medical health care). Samples were sent to eight teams for dose estimation and identification of clinically relevant groups. For quantitative reverse transcription polymerase chain reaction (qRT-PCR) and microarray analyses, samples were lysed, stored at 20°C and shipped on wet ice. RNA isolations and assays were run in each laboratory according to locally established protocols. The time-to-result for both rough early and more precise later reports has been documented where possible. Accuracy of dose estimates was calculated as the difference between estimated and reference doses for all doses (summed absolute difference, SAD) and by determining the number of correctly reported dose estimates that were defined as ±0.5 Gy for reference doses <2.5 Gy and ±1.0 Gy for reference doses >3 Gy, as recommended for triage dosimetry. We also examined the allocation of dose estimates to clinically/diagnostically relevant exposure groups. Altogether, 105 dose estimates were reported by the eight teams, and the earliest report times on dose categories and estimates were 5 h and 9 h, respectively. The coefficient of variation for 85% of all 436 qRT-PCR measurements did not exceed 10%. One team reported dose estimates that systematically deviated several-fold from reported dose estimates, and these outliers were excluded from further analysis. Teams employing a combination of several genes generated about two-times lower median SADs (0.8 Gy) compared to dose estimates based on single genes only (1.7 Gy). When considering the uncertainty intervals for triage dosimetry, dose estimates of all teams together were correctly reported in 100% of the 0 Gy, 50% of the 1.2 Gy and 50% of the 3.5 Gy exposed samples. The order of dose estimates (from lowest to highest) corresponding to three dose categories (unexposed, low dose and highest exposure) were correctly reported by all teams and all chosen genes or gene combinations. Furthermore, if teams reported no exposure or an exposure >3.5 Gy, it was always correctly allocated to the unexposed and the highly exposed group, while low exposed (1.2 Gy) samples sometimes could not be discriminated from highly (3.5 Gy) exposed samples. All teams used FDXR and 78.1% of correct dose estimates used FDXR as one of the predictors. Still, the accuracy of reported dose estimates based on FDXR differed considerably among teams with one team's SAD (0.5 Gy) being comparable to the dose accuracy employing a combination of genes. Using the workflow of this reference team, we performed additional experiments after the exercise on residual RNA and cDNA sent by six teams to the reference team. All samples were processed similarly with the intention to improve the accuracy of dose estimates when employing the same workflow. Re-evaluated dose estimates improved for half of the samples and worsened for the others. In conclusion, this inter-laboratory comparison exercise enabled (1) identification of technical problems and corrections in preparations for future events, (2) confirmed the early and high-throughput capabilities of gene expression, (3) emphasized different biodosimetry approaches using either only FDXR or a gene combination, (4) indicated some improvements in dose estimation with FDXR when employing a similar methodology, which requires further research for the final conclusion and (5) underlined the applicability of gene expression for identification of unexposed and highly exposed samples, supporting medical management in radiological or nuclear scenarios.

What We Have Learned from RENEB Inter-Laboratory Comparisons Since 2012 With Focus on ILC 2021

Inter-laboratory exercises are important tools within the European network for biological dosimetry and physical retrospective dosimetry (RENEB) to validate and improve the performance of member laboratories and to ensure an operational network with high quality standards for dose estimations in case of a large-scale radiological or nuclear event. In addition to the RENEB inter-laboratory comparison 2021, several inter-laboratory comparisons have been performed in the frame of RENEB for a number of assays in recent years. This publication gives an overview of RENEB inter-laboratory comparisons for biological dosimetry assays in the past and a final summary of the challenges and lessons learnt from the RENEB inter-laboratory comparison 2021. In addition, the dose estimates of all RENEB inter-laboratory comparisons since 2013 that have been conducted for the dicentric chromosome assay, the most established and applied assay, are compared and discussed.

RENEB Inter-Laboratory Comparison 2021: The FISH-Based Translocation Assay

Translocation analysis using fluorescence in situ hybridization (FISH) is the method of choice for dose assessment in case of chronic or past exposures to ionizing radiation. Although it is a widespread technique, unlike dicentrics, the number of FISH-based inter-laboratory comparisons is small. For this reason, although the current Running the European Network of Biological and Physical retrospective Dosimetry (RENEB) inter-laboratory comparison 2021 was designed as a fast response to a real emergency scenario, it was considered a good opportunity to perform an inter-laboratory comparison using the FISH technique to gain further experience. The Bundeswehr Institute of Radiobiology provided peripheral blood samples from one healthy human volunteer. Three test samples were irradiated with blinded doses of 0, 1.2, and 3.5 Gy, respectively. Samples were then sent to the seven participating laboratories. The FISH technique was applied according to the standard procedure of each laboratory. Both, the frequency of translocations and the estimated dose for each sample were sent to the coordinator using a special scoring sheet for FISH. All participants sent their results in due time. However, although it was initially requested to send the results based on the full analysis, evaluating 500 equivalent cells, most laboratories only sent the results based on triage, with a smaller number of analyzed cells. In the triage analysis, there was great heterogeneity in the number of equivalent cells scored. On the contrary, for the full analysis, this number was more homogeneous. For all three samples, one laboratory showed outlier yields compared to the other laboratories. Excluding these results, in the triage analysis, the frequency of translocations in sample no. 1 ranged from 0 to 0.013 translocations per cell, and for samples no. 2 and no. 3 the genomic mean frequency were 0.27 ± 0.03 and 1.47 ± 0.14, with a coefficient of variation of 0.29 and 0.23 respectively. Considering only results obtained in the triage analysis for sample no. 1, all laboratories, except one, classified this sample as the non-irradiated one. For sample no. 2, excluding the outlier value, the mean reported dose was 1.74 ± 0.16 Gy indicating a mean deviation of about 0.5 Gy to the delivered dose of 1.2 Gy. For sample no. 3 the mean dose estimated was 4.21 ± 0.21 Gy indicating a mean deviation of about 0.7 Gy to the delivered dose of 3.5 Gy. In the frame of RENEB, this is the second FISH-based inter-laboratory comparison. The whole exercise was planned as a response to an emergency, therefore, a triage analysis was requested for all the biomarkers except for FISH. Although a full analysis was initially requested for FISH, most of the laboratories reported only a triage-based result. The main reason is that it was not clearly stated what was required before starting the exercise. Results show that most of the laboratories successfully discriminated unexposed and irradiated samples from each other without any overlap. A good agreement in the observed frequencies of translocations was observed but there was a tendency to overestimate the delivered doses. Efforts to improve the harmonization of this technique and subsequent exercises to elucidate the reason for this trend should be promoted.

RENEB Inter-Laboratory Comparison 2021: The Cytokinesis-Block Micronucleus Assay

The goal of the RENEB inter-laboratory comparison 2021 exercise was to simulate a large-scale radiation accident involving a network of biodosimetry labs. Labs were required to perform their analyses using different biodosimetric assays in triage mode scoring and to rapidly report estimated radiation doses to the organizing institution. This article reports the results obtained with the cytokinesis-block micronucleus assay. Three test samples were exposed to blinded doses of 0, 1.2 and 3.5 Gy X-ray doses (240 kVp, 13 mA, ∼75 keV, 1 Gy/min). These doses belong to 3 triage categories of clinical relevance: a low dose category, for no exposure or exposures inferior to 1 Gy, requiring no direct treatment of subjects; a medium dose category, with doses ranging from 1 to 2 Gy, and a high dose category, after exposure to doses higher than 2 Gy, with the two latter requiring increasing medical attention. After irradiation the test samples (no. 1, no. 2 and no. 3) were sent by the organizing laboratory to 14 centers participating in the micronucleus assay exercise. Laboratories were asked to setup micronucleus cultures and to perform the micronucleus assay in triage mode, scoring 500 binucleated cells manually, or 1,000 binucleated cells in automated/semi-automated mode. One laboratory received no blood samples, but scored pictures from another lab. Based on their calibration curves, laboratories had to provide estimates of the administered doses. The accuracy of the reported dose estimates was further analyzed by the micronucleus assay lead. The micronucleus assay allowed classification of samples in the corresponding clinical triage categories (low, medium, high dose category) in 88% of cases (manual scoring, 88%; semi-automated scoring, 100%; automated scoring, 73%). Agreement between scoring laboratories, assessed by calculating the Fleiss' kappa, was excellent (100%) for semi-automated scoring, good (83%) for manual scoring and poor (53%) for fully automated scoring. Correct classification into triage scoring dose intervals (reference dose ±0.5 Gy for doses ≤2.5 Gy, or reference dose ±1 Gy for doses >2.5 Gy), recommended for triage biodosimetry, was obtained in 79% of cases (manual scoring, 73%; semi-automated scoring, 100%; automated scoring, 67%). The percentage of dose estimates whose 95% confidence intervals included the reference dose was 58% (manual scoring, 48%; semi-automated scoring, 72%; automated scoring, 60%). For the irradiated samples no. 2 and no. 3, a systematic shift towards higher dose estimations was observed. This was also noticed with the other cytogenetic assays in this intercomparison exercise. Accuracy of the rapid triage modality could be maintained when the number of manually scored cells was scaled down to 200 binucleated cells. In conclusion, the micronucleus assay, preferably performed in a semi-automated or manual scoring mode, is a reliable technique to perform rapid biodosimetry analysis in large-scale radiation emergencies.

RENEB Inter-Laboratory Comparison 2021: The Dicentric Chromosome Assay

After large-scale radiation accidents where many individuals are suspected to be exposed to ionizing radiation, biological and physical retrospective dosimetry assays are important tools to aid clinical decision making by categorizing individuals into unexposed/minimally, moderately or highly exposed groups. Quality-controlled inter-laboratory comparisons of simulated accident scenarios are regularly performed in the frame of the European legal association RENEB (Running the European Network of Biological and Physical retrospective Dosimetry) to optimize international networking and emergency readiness in case of large-scale radiation events. In total 33 laboratories from 22 countries around the world participated in the current RENEB inter-laboratory comparison 2021 for the dicentric chromosome assay. Blood was irradiated in vitro with X rays (240 kVp, 13 mA, ∼75 keV, 1 Gy/min) to simulate an acute, homogeneous whole-body exposure. Three blood samples (no. 1: 0 Gy, no. 2: 1.2 Gy, no. 3: 3.5 Gy) were sent to each participant and the task was to culture samples, to prepare slides and to assess radiation doses based on the observed dicentric yields from 50 manually or 150 semi-automatically scored metaphases (triage mode scoring). Approximately two-thirds of the participants applied calibration curves from irradiations with γ rays and about 1/3 from irradiations with X rays with varying energies. The categorization of the samples in clinically relevant groups corresponding to individuals that were unexposed/minimally (0-1 Gy), moderately (1-2 Gy) or highly exposed (>2 Gy) was successfully performed by all participants for sample no. 1 and no. 3 and by ≥74% for sample no. 2. However, while most participants estimated a dose of exactly 0 Gy for the sham-irradiated sample, the precise dose estimates of the samples irradiated with doses >0 Gy were systematically higher than the corresponding reference doses and showed a median deviation of 0.5 Gy (sample no. 2) and 0.95 Gy (sample no. 3) for manual scoring. By converting doses estimated based on γ-ray calibration curves to X-ray doses of a comparable mean photon energy as used in this exercise, the median deviation decreased to 0.27 Gy (sample no. 2) and 0.6 Gy (sample no. 3). The main aim of biological dosimetry in the case of a large-scale event is the categorization of individuals into clinically relevant groups, to aid clinical decision making. This task was successfully performed by all participants for the 0 Gy and 3.5 Gy samples and by 74% (manual scoring) and 80% (semi-automatic scoring) for the 1.2 Gy sample. Due to the accuracy of the dicentric chromosome assay and the high number of participating laboratories, a systematic shift of the dose estimates could be revealed. Differences in radiation quality (X ray vs. γ ray) between the test samples and the applied dose effect curves can partly explain the systematic shift. There might be several additional reasons for the observed bias (e.g., donor effects, transport, experimental conditions or the irradiation setup) and the analysis of these reasons provides great opportunities for future research. The participation of laboratories from countries around the world gave the opportunity to compare the results on an international level.

Gene Expression Changes in a Prefinal Health Stage of Lethally Irradiated Male and Female Rhesus Macaques

Radiation-induced gene expression (GE) changes can be used for early and high-throughput biodosimetry within the first three days postirradiation. However, is the method applicable in situations such as the Alexander Litvinenko case or the Goiania accident, where diagnosis occurred in a prefinal health stage? We aimed to characterize gene expression changes in a prefinal health stage of lethally irradiated male and female rhesus macaques. Peripheral blood was drawn pre-exposure and at the prefinal stage of male and female animals, which did not survive whole-body exposure with 700 cGy (LD66/60). RNA samples originated from a blinded randomized Good Laboratory Practice study comprising altogether 142 irradiated rhesus macaques of whom 60 animals and blood samples (15 samples for both time points and sexes) were used for this analysis. We evaluated GE on 34 genes widely used in biodosimetry and prediction of the hematological acute radiation syndrome severity (H-ARS) employing quantitative real-time polymerase chain reaction (qRT-PCR). These genes were run in duplicate and triplicate and altogether 96 measurements per time point and sex could be performed. In addition, 18S ribosomal RNA (rRNA) was measured to depict the ribosome/transcriptome status as well as for normalization purposes and 16S rRNA was evaluated as a surrogate for bacteremia. Mean differential gene expression (DGE) was calculated for each gene and sex including all replicate measurements and using pre-exposure samples as the reference. From 34 genes, altogether 27 genes appeared expressed. Pre-exposure samples revealed no signs of bacteremia and 18S rRNA GE was in the normal range in all 30 samples. Regarding prefinal samples, 46.7% and 40% of animals appeared infected in females and males, respectively, and for almost all males this was associated with out of normal range 18S rRNA values. The total number of detectable GE measurements was sixfold (females) and 15-fold (males) reduced in prefinal relative to pre-exposure samples and about tenfold lower in 80% of prefinal compared to pre-exposure samples (P < 0.0001). An overall 11-fold (median) downregulation in prefinal compared to pre-exposure samples was identified for most of the 27 genes and even FDXR appeared 4-14-fold downregulated in contrast to a pronounced up-regulation according to cited work. This pattern of overall downregulation of almost all genes and the rapid reduction of detectable genes at a prefinal stage was found in uninfected animals with normal range 18S rRNA as well. In conclusion, in a prefinal stage after lethal radiation exposure, the ribosome/transcriptome status remains present (based on normal range 18S rRNA values) in 60-67% of animals, but the whole transcriptome activity in general appears silenced and cannot be used for biodosimetry purposes, but probably as an indicator for an emerging prefinal health stage.

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 (ED₅₀) 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.

Early molecular markers for retrospective biodosimetry and prediction of acute health effects

Radiation-induced biological changes occurring within hours and days after irradiation can be potentially used for either exposure reconstruction (retrospective dosimetry) or the prediction of consecutively occurring acute or chronic health effects. The advantage of molecular protein or gene expression (GE) (mRNA) marker lies in their capability for early (1-3 days after irradiation), high-throughput and point-of-care diagnosis, required for the prediction of the acute radiation syndrome (ARS) in radiological or nuclear scenarios. These molecular marker in most cases respond differently regarding exposure characteristics such as e.g. radiation quality, dose, dose rate and most importantly over time. Changes over time are in particular challenging and demand certain strategies to deal with. With this review, we provide an overview and will focus on already identified and used mRNA GE and protein markers of the peripheral blood related to the ARS. These molecules are examined in light of 'ideal' characteristics of a biomarkers (e.g. easy accessible, early response, signal persistency) and the validation degree. Finally, we present strategies on the use of these markers considering challenges as their variation over time and future developments regarding e.g. origin of samples, point of care and high-throughput diagnosis.

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.

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.

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.

Analysis of Gene Expression Changes in PHA-M Stimulated Lymphocytes - Unraveling PHA Activity as Prerequisite for Dicentric Chromosome Analysis

Dicentric chromosome analysis (DCA) is the gold standard for individual radiation dose assessment. However, DCA is limited by the time-consuming phytohemagglutinin (PHA)-mediated lymphocyte activation. In this study using human peripheral blood lymphocytes, we investigated PHA-associated whole genome gene expression changes to elucidate this process and sought to identify suitable gene targets as a means of meeting our long-term objective of accelerating cell cycle kinetics to reduce DCA culture time. Human peripheral whole blood from three healthy donors was separately cultured in RPMI/FCS/antibiotics with BrdU and PHA-M. Diluted whole blood samples were transferred into PAXgene tubes at 0, 12, 24 and 36 h culture time. RNA was isolated and aliquots were used for whole genome gene expression screening. Microarray results were validated using qRT-PCR and differentially expressed genes [significantly (FDR corrected) twofold different from the 0 h value reference] were analyzed using several bioinformatic tools. The cell cycle positions and DNA-synthetic activities of lymphocytes were determined by analyzing the correlated total DNA content and incorporated BrdU level with flow cytometry after continued BrdU incubation. From 42,545 transcripts of the whole genome microarray 47.6%, on average, appeared expressed. The number of differentially expressed genes increased linearly from 855 to 2,858 and 4,607 at 12, 24 and 36 h after PHA addition, respectively. Approximately 2-3 times more up- than downregulated genes were observed with several hundred genes differentially expressed at each time point. Earliest enrichment was observed for gene sets related to the nucleus (12 h) followed by genes assigned to intracellular structures such as organelles (24 h) and finally genes related to the membrane and the extracellular matrix were enriched (36 h). Early gene expression changes at 12 h, in particular, were associated with protein classes such as chemokines/cytokines (e.g., CXCL1, CXCL2) and chaperones. Genes coding for biological processes involved in cell cycle control (e.g., MYBL2, RBL1, CCNA, CCNE) and DNA replication (e.g., POLA, POLE, MCM) appeared enriched at 24 h and later, but many more biological processes (42 altogether) showed enrichment as well. Flow cytometry data fit together with gene expression and bioinformatic analyses as cell cycle transition into S phase was observed with interindividual differences from 12 h onward, whereas progression into G₂ as well as into the second G₁ occurred from 36 h onward after activation. Gene set enrichment analysis over time identifies, in particular, two molecular categories of PHA-responsive gene targets (cytokine and cell cycle control genes). Based on that analysis target genes for cell cycle acceleration in lymphocytes have been identified ( CDKN1A/B/C, RBL-1/RBL-2, E2F2, Deaf-1), and it remains undetermined whether the time expenditure for DCA can be reduced by influencing gene expression involved in the regulatory circuits controlling PHA-associated cell cycle entry and/or progression at a specific early cell cycle phase.

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.

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.

Comparing seven mitogens with PHA-M for improved lymphocyte stimulation in dicentric chromosome analysis for biodosimetry

Dicentric chromosome analysis (DCA) is the gold standard for individual radiation dose estimation. Two limiting factors of DCA are the time-consuming lymphocyte stimulation and proliferation using the lectin PHA-M and the upper dose limit of individual dose assessment of ∼4 Gy. By measuring the mitotic index (MI), the authors investigated systematically whether the stimulation of lymphocytes can be improved after administration of alternative (and combined) mitogens. The authors compared the lymphocyte stimulation effectiveness of the traditionally used PHA-M (from Phaseolus vulgaris) with seven cited mitogens by determination of MIs: five lectins namely CNA (concanavalin A), PW (pokeweed), LMA (Maackia amurensis), LTV (T. vulgaris), PHA-L (P. vulgaris) as well as LPS (lipopolysaccharide, Escherichia coli) and SLO (streptolysine O, Streptococcus pyogenes) were applied. The conventional protocol using PHA-M for lymphocyte stimulation proved to be superior over lower/higher PHA-M concentrations as well as seven other mitogens administered either alone or combined with SLO or LPS.

Detection of hypoxia markers in the cerebellum after a traumatic frontal cortex injury: a human postmortem gene expression analysis

PURPOSE

The response to traumatic brain injury (TBI) is complex and induces various biological pathways in all brain regions that contribute to bad outcomes. The cerebellar hypoxia after a frontal cortex injury may potentiate the pathophysiological impacts of TBI. Therefore, a gene expression analysis was conducted to determine the influence of hypoxia on TBIs.

MATERIAL AND METHODS

Total RNA, including microRNAs, was isolated from the cerebellum of individuals who had died from severe frontal cortex injuries or due to natural causes of death (reference group).

RESULTS

From a total of 19,596 genes, an average of 59.56% messenger RNAs (mRNAs) appeared expressed with 42 of them showing significant >2-fold differences of upregulated (n = 18) and downregulated (n = 24) genes. The validity of 14 candidate genes (with low p values and high fold differences or based on cited literature) was confirmed using qRT-PCR (Spearman correlation r(2) = 0.93). Only four genes appeared to be either upregulated (FOSB and IL6) or downregulated (HSD11B1 and HSPA12B). From a total of 667 microRNAs, altogether, 248 microRNAs appeared expressed with 13 of them showing significant differences in the mean gene expression. The combination of two mRNAs (HSPA12B/FOSB or IL6/HSD11B1) or two microRNAs (either miR-138/miR-744 or miR-195/miR-324-5p) completely discriminated both groups, a finding unaltered by potential confounders such as age at biosampling, survival time, and the postmortem interval.

CONCLUSIONS

Cerebellar hypoxia markers are important to understand the pathophysiology of TBIs and could be used for therapeutic strategies or forensic purposes, e.g., to assess the severity of a brain injury.

Discriminating metastasised from non-metastasised seminoma based on transcriptional changes in primary tumours using NGS

Background:

We aimed to better discriminate (occult) metastasised from non-metastasised seminoma based on transcriptional changes of small RNAs in the primary tumour.

Methods:

Total RNAs including small RNAs were isolated from five testicular tumours of each, lymphogenic, occult and non-metastasised patients. Next-generation sequencing (SOLID, Life Technologies) was used to examine transcriptional changes. Small RNAs showing ⩾50 reads and a significant ⩾2-fold difference using non-metastasised tumours as the reference group were examined in univariate logistic regression analysis and combinations of two small RNAs were further examined using support vector machines.

Results:

On average, 1.3 × 10⁷, 1.4 × 10⁷ and 1.7 × 10⁷ small RNA reads were detectable in non-metastasised, occult and lymphogenic metastasised seminoma, respectively, of which 30–32% remained after trimming. Between 59 and 68% represented annotated reads and between 8.6 and 11% were annotated small RNA tags. Of them, 137 small RNAs showed>50 reads and a two-fold difference to the reference. In univariate analysis, 32–38 small RNAs significantly discriminated lymphogenic/occult from non-metastasised seminoma, and among these different comparisons, it were the same small RNAs in 51–88%. Many combinations of two of these small RNAs allowed a complete discrimination of metastasised from non-metastasised seminoma irrespective of the metastasis subtype.

Conclusions:

Metastasised and non-metastasised seminoma can be completely discriminated with a combination of two small RNAs.

Iodine-131 dose-dependent gene expression: alterations in both normal and tumour thyroid tissues of post-Chernobyl thyroid cancers

BACKGROUND

A strong, consistent association between childhood irradiation and subsequent thyroid cancer provides an excellent model for studying radiation carcinogenesis.

METHODS

We evaluated gene expression in 63 paired RNA specimens from frozen normal and tumour thyroid tissues with individual iodine-131 (I-131) doses (0.008-8.6 Gy, no unirradiated controls) received from Chernobyl fallout during childhood (Ukrainian-American cohort). Approximately half of these randomly selected samples (32 tumour/normal tissue RNA specimens) were hybridised on 64 whole-genome microarrays (Agilent, 4 × 44 K). Associations between I-131 dose and gene expression were assessed separately in normal and tumour tissues using Kruskal-Wallis and linear trend tests. Of 155 genes significantly associated with I-131 after Bonferroni correction and with ≥2-fold increase per dose category, we selected 95 genes. On the remaining 31 RNA samples these genes were used for validation purposes using qRT-PCR.

RESULTS

Expression of eight genes (ABCC3, C1orf9, C6orf62, FGFR1OP2, HEY2, NDOR1, STAT3, and UCP3) in normal tissue and six genes (ANKRD46, CD47, HNRNPH1, NDOR1, SCEL, and SERPINA1) in tumour tissue was significantly associated with I-131. PANTHER/DAVID pathway analyses demonstrated significant over-representation of genes coding for nucleic acid binding in normal and tumour tissues, and for p53, EGF, and FGF signalling pathways in tumour tissue.

CONCLUSION

The multistep process of radiation carcinogenesis begins in histologically normal thyroid tissue and may involve dose-dependent gene expression changes.

Laboratory intercomparison of gene expression assays

The possibility of a large-scale acute radiation exposure necessitates the development of new methods that could provide rapid individual dose estimates with high sample throughput. The focus of the study was an intercomparison of laboratories' dose-assessment performances using gene expression assays. Lithium-heparinized whole blood from one healthy donor was irradiated (240 kVp, 1 Gy/min) immediately after venipuncture at approximately 37°C using single X-ray doses. Blood samples to establish calibration curves (0.25-4 Gy) as well as 10 blinded test samples (0.1-6.4 Gy) were incubated for 24 h at 37°C supplemented with an equal volume of medium and 10% fetal calf serum. For quantitative reverse transcription polymerase chain reaction (qRT-PCR), samples were lysed, stored at -20°C and shipped on ice. For the Chemical Ligation Dependent Probe Amplification methodology (CLPA), aliquots were incubated in 2 ml CLPA reaction buffer (DxTerity), mixed and shipped at room temperature. Assays were run in each laboratory according to locally established protocols. The mean absolute difference (MAD) of estimated doses relative to the true doses (in Gy) was calculated. We also merged doses into binary categories reflecting aspects of clinical/diagnostic relevance and examined accuracy, sensitivity and specificity. The earliest reported time on dose estimates was <8 h. The standard deviation of technical replicate measurements in 75% of all measurements was below 11%. MAD values of 0.3-0.5 Gy and 0.8-1.3 Gy divided the laboratories contributions into two groups. These fourfold differences in accuracy could be primarily explained by unexpected variances of the housekeeping gene (P = 0.0008) and performance differences in processing of calibration and blinded test samples by half of the contributing laboratories. Reported gene expression dose estimates aggregated into binary categories in general showed an accuracies and sensitivities of 93-100% and 76-100% for the groups, with low MAD and high MAD, respectively. In conclusion, gene expression-based dose estimates were reported quickly, and for laboratories with MAD between 0.3-0.5 Gy binary dose categories of clinical significance could be discriminated with an accuracy and sensitivity comparable to established cytogenetic assays.

Laboratory intercomparison on the γ-H2AX foci assay

The focus of the study is an intercomparison of laboratories' dose-assessment performances using the γ-H2AX foci assay as a diagnostic triage tool for rapid individual radiation dose assessment. Homogenously X-irradiated (240 kVp, 1 Gy/min) blood samples for establishing calibration data (0.25-4 Gy) as well as blinded test samples (0.1-6.4 Gy) were incubated at 37°C for 2 and 24 h (repair time) and sent to the participants. The foci assay was performed according to protocols individually established in participating laboratories and therefore varied. The time taken to report dose estimates was documented for each laboratory. Additional information concerning laboratory organization/characteristics as well as assay performance was collected. The mean absolute difference (MAD) of estimated doses relative to the actual doses was calculated and radiation doses were merged into four triage categories reflecting clinical relevance to calculate accuracy, sensitivity and specificity. First γ-H2AX based dose estimates were reported 7 h after sample receipt. Estimates were similarly accurate for 2 and 24 h repair times, providing scope for its use in the early phase of a radiation exposure incident. Equal accuracy was achieved by scoring 20, 30, 40 or 50 cells per sample. However, MAD values of 0.5-0.7 Gy and 1.3-1.7 Gy divided the data sets into two groups, driven mainly by the considerable differences in foci yields between calibration and blind samples. Foci yields also varied dramatically between laboratories, highlighting reproducibility issues as an important caveat of the foci assay. Nonetheless, foci counts could distinguish high- and low-dose samples in all data sets and binary dose categories of clinical significance could be discriminated with satisfactory accuracy (mean 84%, ±0.03 SEM). Overall, the results suggest that the γ-H2AX assay is a useful tool for rapidly screening individuals for significant exposures that occurred up to at least 24 h earlier, and may help to prioritize cytogenetic dosimetry follow-up.

Comparison of established and emerging biodosimetry assays

Rapid biodosimetry tools are required to assist with triage in the case of a large-scale radiation incident. Here, we aimed to determine the dose-assessment accuracy of the well-established dicentric chromosome assay (DCA) and cytokinesis-block micronucleus assay (CBMN) in comparison to the emerging γ-H2AX foci and gene expression assays for triage mode biodosimetry and radiation injury assessment. Coded blood samples exposed to 10 X-ray doses (240 kVp, 1 Gy/min) of up to 6.4 Gy were sent to participants for dose estimation. Report times were documented for each laboratory and assay. The mean absolute difference (MAD) of estimated doses relative to the true doses was calculated. We also merged doses into binary dose categories of clinical relevance and examined accuracy, sensitivity and specificity of the assays. Dose estimates were reported by the first laboratories within 0.3-0.4 days of receipt of samples for the γ-H2AX and gene expression assays compared to 2.4 and 4 days for the DCA and CBMN assays, respectively. Irrespective of the assay we found a 2.5-4-fold variation of interlaboratory accuracy per assay and lowest MAD values for the DCA assay (0.16 Gy) followed by CBMN (0.34 Gy), gene expression (0.34 Gy) and γ-H2AX (0.45 Gy) foci assay. Binary categories of dose estimates could be discriminated with equal efficiency for all assays, but at doses ≥1.5 Gy a 10% decrease in efficiency was observed for the foci assay, which was still comparable to the CBMN assay. In conclusion, the DCA has been confirmed as the gold standard biodosimetry method, but in situations where speed and throughput are more important than ultimate accuracy, the emerging rapid molecular assays have the potential to become useful triage tools.

Gene expression comparisons performed for biodosimetry purposes on in vitro peripheral blood cellular subsets and irradiated individuals

We examined the benefit of gene expression analysis on peripheral blood cellular subsets of different radiosensitivity to elucidate their utility as biodosimeters for estimation of dose in irradiated individuals. Peripheral mononucleated cells were isolated from 18 healthy volunteers employing density separation in a CPT-NH tube. Peripheral mononucleated cells were cultured in RPMI 1640 medium containing 10% autologous serum and were irradiated with 0.1-1 Gy (240 kV, 13 mA, X rays at 1 Gy/min). A low-dose study was performed with isolated peripheral mononucleated cells from one healthy donor in three independent experiments. Peripheral mononucleated cells were irradiated at 0 (sham), 1, 2.5 and 5 cGy (70 kV, 13 mA X rays at 1 cGy/min) and gene expression was measured 24 and 48 h after irradiation. After irradiation, CD4(+) or CD8(+) cells were isolated by magnetic beads in independent experiments. RNA from lymphocyte subsets and peripheral mononucleated cells was isolated after 24 and 48 h and converted into cDNA. Gene expression of GADD45, CDKN1A, DDB2, PCNA, BAX and ATF3 were determined using RTQ-PCR. Data were analyzed employing linear and logistic regression analysis. The same examinations were performed in 5 individuals either diagnosed using CT scans (up to 4.3 cGy) or by administering (F-18)-fluoro-2-deoxy-d-glucose (F-18 FDG, 0.6 cGy). Methodological, intra- and inter-individual variability in 90-95% of measurements did not exceed the introduced twofold change over sham-irradiated control values in peripheral mononucleated cells and CD4(+) cells, and therefore no false positive results were observed. Dose reconstruction in peripheral mononucleated cells in opposite to CD4(+) lymphocytes required fewer genes and appeared more efficient (R-square = 84.8% compared to 51.8%). In vitro samples exposed to 10 cGy could be completely discriminated from sham-irradiated samples without individual pre-exposure controls, which coincided with our preliminary in vivo results. However, in vitro differential gene expression was measured relative to control values and did not differ significantly at 24 and 48 h after irradiation in contrast to our preliminary in vivo data. In addition, below 5 cGy in vitro data did not show reproducible significant changes in gene expression, which was opposite to our preliminary in vivo data. Therefore a twofold change in gene expression over control sufficiently controls for different sources of variance, and measuring gene expression in peripheral mononucleated cell for biological dosimetry purposes appears superior over measurements in lymphocyte subsets. The increased gene expression measured after low absorbed doses in vivo and in vitro might indicate a particular applicability of this method for a low-level radiation scenario in the absence of individual pre-exposure controls. However, the constant gene expression values measured up to 48 h in our in vitro model at doses >10 cGy, and the absence of reproducible and statistically significant gene expression changes below 5 cGy contrast to the preliminary in vivo results performed at similar doses. Therefore, measurements with our in vitro models should be interpreted cautiously.