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.

The Influence of Computed Tomography Contrast Agent on Radiation-Induced Gene Expression and Double-Strand Breaks

After nuclear scenarios, combined injuries of acute radiation syndrome (ARS) with, e.g., abdominal trauma, will occur and may require contrast-enhanced computed tomography (CT) scans for diagnostic purposes. Here, we investigated the effect of iodinated contrast agents on radiation-induced gene expression (GE) changes used for biodosimetry (AEN, BAX, CDKN1A, EDA2R, APOBEC3H) and for hematologic ARS severity prediction (FDXR, DDB2, WNT3, POU2AF1), and on the induction of double-strand breaks (DSBs) used for biodosimetry. Whole blood samples from 10 healthy donors (5 males, 5 females, mean age: 28 ± 2 years) were irradiated with X rays (0, 1 and 4 Gy) with and without the addition of iodinated contrast agent (0.016 ml contrast agent/ml blood) to the blood prior to the exposure. The amount of contrast agent was set to be equivalent to the blood concentration of an average patient (80 kg) during a contrast-enhanced CT scan. After irradiation, blood samples were incubated at 37°C for 20 min (DSB) and 8 h (GE, DSB). GE was measured employing quantitative real-time polymerase chain reaction. DSB foci were revealed by γH2AX + 53BP1 immunostaining and quantified automatically in >927 cells/sample. Radiation-induced differential gene expression (DGE) and DSB foci were calculated using the respective unexposed sample without supplementation of contrast agent as the reference. Neither the GE nor the number of DSB foci was significantly (P = 0.07-0.94) altered by the contrast agent application. However, for some GE and DSB comparisons with/without contrast agent, there were weakly significant differences (P = 0.03-0.04) without an inherent logic and thus are likely due to inter-individual variation. In nuclear events, the diagnostics of combined injuries can require the use of an iodinated contrast agent, which, according to our results, does not alter or influence radiation-induced GE changes and the quantity of DSB foci. Therefore, the gene expression and γH2AX focus assay can still be applied for biodosimetry and/or hematologic ARS severity prediction in such scenarios.

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.

Lymphocyte Function at Baseline Could Be a New Predictor of Tumor Burden following Six Cycles of Radium-223 Therapy in Patients with Metastasized, Castration-Resistant Prostate Cancer

Previous data indicate that one cycle of treatment with radium-223 (223Ra) did not significantly impair lymphocyte function in patients with metastasized, castration-resistant prostate cancer. The aim of the current study was to assess in 21 patients whether six cycles of this therapy had an effect on lymphocyte proliferation and interferon-γ and interleukin (IL)-10 ELISpot results. Lymphocyte proliferation after stimulation with microbial antigens and the production of interferon-γ continuously decreased after six cycles of radionuclide therapy, reaching statistical significance (p < 0.05) at months 1, 2, 4, and/or 6 after therapy. One month after the last cycle of therapy, 67% of patients showed a decrease in tumor burden. The tumor burden correlated negatively with IL-10 secretion at baseline, e.g., after stimulation with tetanus antigen (p < 0.0001, r = -0.82). As determined by receiver operating characteristic (ROC) curve analysis, tetanus-specific IL-10 spots at baseline had the highest predictive value (p = 0.005) for tumor burden at month 6, with an area under the curve (AUC) of 0.90 (sensitivity 100%, specificity 78%). In conclusion, we observed an additive effect of treatment with 223Ra on immune function and found that IL-10 secretion at baseline predicted tumor burden at month 6 after treatment.

Magnetic Resonance Image Findings and Potential Anatomic Risk Factors for Chodromalacia in Children and Adolescents Suffering from Non-Overload Atraumatic Knee Pain in the Ambulant Setting

PURPOSE

To evaluate magnetic resonance image (MRI) findings in children and adolescents suffering from knee pain without traumatic or physical overload history and to identify potential anatomic risk factors.

MATERIAL AND METHODS

A total of 507 MRIs of 6- to 20-year-old patients (251 males; 256 females) were evaluated with regard to detectable pathologies of the knee. The results were compared to a control group without pain (n = 73; 34 males; 39 females). A binary logistic regression model and t-tests for paired and unpaired samples were used to identify possible risk factors and significant anatomic differences of the study population.

RESULTS

In 348 patients (68.6%), at least one pathology was detected. The most commonly detected finding was chondromalacia of the patellofemoral (PF) joint (n = 205; 40.4%). Chondral lesions of the PF joint occurred significantly more often in knee pain patients than in the control group (40% vs. 11.0%; p = 0.001), especially in cases of a patella tilt angle > 5° (p ≤ 0.001), a bony sulcus angle > 150° (p = 0.002), a cartilaginous sulcus angle > 150° (p = 0.012), a lateral trochlear inclination < 11° (p ≤ 0.001), a lateralised patella (p = 0.023) and a Wiberg type II or III patella shape (p = 0.019). Moreover, a larger patella tilt angle (p = 0.021), a greater bony sulcus angle (p = 0.042), a larger cartilaginous sulcus angle (p = 0.038) and a lower value of the lateral trochlear inclination (p = 0.014) were detected in knee pain patients compared to the reference group.

CONCLUSION

Chondromalacia of the PF joint is frequently observed in children and adolescents suffering from non-overload atraumatic knee pain, whereby a patella tilt angle > 5°, a bony sulcus angle > 150°, a cartilaginous sulcus angle > 150°, a lateral trochlear inclination < 11°, a lateralised patella and a Wiberg type II or III patella shape seem to represent anatomic risk factors.

Validation of genes for H-ARS severity prediction in leukemia patients - interspecies comparison, challenges, and promises

PURPOSE

In a previous baboon-study, a total of 29 genes were identified for clinical outcome prediction of the hematologic, acute, radiation, syndrome (H-ARS) severity. Among them, four genes (FDXR, DDB2, POU2AF1, WNT3) appeared promising and were validated in five leukemia patients. Within this study, we sought further in-vivo validation in a larger number of whole-body irradiated patients.

MATERIAL AND METHODS

Peripheral blood was drawn from 10 leukemia patients before and up to 3 days during a fractionated (2 Gy/day) total-body irradiation (TBI) with 2-12Gy. After RNA-isolation, gene expression (GE) was evaluated on 31 genes widely used in biodosimetry and H-ARS prediction employing qRT-PCR. A customized low-density-array (LDA) allowed simultanously analyzing all genes, the 96-well format further examined the four most promising genes. Fold-changes (FC) in GE relative to pre-irradiation were calculated.

RESULTS

Five patients suffering from acute-lymphoblastic-leukemia (ALL) respectively non-Hodgkin-lymphoma (NHL) revealed sufficient RNA-amounts and corresponding lymphocyte and neutrophile counts for running qRT-PCR, while acute-myeloid-leukemia (AML) and one myelofibrosis patient could not supply enough RNA. Generally, 1-2µg total RNA was isolated, whereas up to 10-fold differences in RNA-quantities (associated suppressed GE-changes) were identified among pre-exposure and exposure samples. From 31 genes, 23 were expressed in at least one of the pre-exposure samples. Relative to pre-exposure, the number of expressed genes could halve at 48 and 72h after irradiation. Using the LDA, 13 genes were validated in human samples. The four most promising genes (vid. sup.) were either undetermined or too close to pre-exposure. However, they were measured using the more sensitive 96-well format, except WNT3, which wasn´t detectable. As in previous studies, an opposite regulation in GE for FDXR in leukemia patients (up-regulated) relative to baboons (down-regulated) was reconfirmed. Radiation-induced GE-changes of DDB2 (up-regulated) and POU2AF1 (down-regulated) behaved similarly in both species. Hence, 16 out of 23 genes of two species showed GE-changes in the same direction, and up-regulated FDXR as in human studies were revalidated.

CONCLUSION

Identified genes for H-ARS severity prediction, previously detected in baboons, were validated in ALL but not in AML patients. Limitations related to leukemia type, associated reduced RNA amounts, suppressed GE changes, and methodological challenges must be considered as factors negatively affecting the total number of validated genes. Based on that, we propose additional controls including blood cell counts and preferably fluorescence-based RNA quantity measurements for selecting promising samples and using a more sensitive 96-well format for candidate genes with low baseline copy numbers.

Evaluating transport conditions of conventional, widely used EDTA blood tubes for gene expression analysis in comparison to expensive, specialized PAXgene tubes in preparedness for radiological and nuclear events

PURPOSE

Gene expression (GE) analysis of a radio-sensitive gene set (FDXR, DDB2, WNT3, POU2AF1) has been introduced in the last decade as an early and high-throughput prediction tool of later developing acute hematologic radiation syndrome (H-ARS) severity. The use of special tubes for RNA extraction from peripheral whole blood (PAXgene) represent an established standard in GE studies, although uncommonly used in clinics and not immediately available in the quantities needed in radiological/nuclear (R/N) incidents. On the other hand, EDTA blood tubes are widely utilized in clinical practice.

MATERIAL AND METHODS

Using blood samples from eleven healthy donors, we investigated GE changes associated with delayed processing of EDTA tubes up to 4 h at room temperature (RT) after venipuncture (simulating delays caused by daily clinical routine), followed by a subsequent transport time of 24 h at RT, 4 °C, and -20 °C. Differential gene expression (DGE) of the target genes was further examined after X-irradiation with 0 Gy and 4 Gy under optimal transport conditions.

RESULTS

No significant changes in DGE were observed when storing EDTA whole blood samples up to 4 h at RT and subsequently kept at 4 °C for 24 h which is in line with expected DGE. However, other storage conditions, such as -20 °C or RT, decreased RNA quality and/or (significantly) caused changes in DGE exceeding the known methodological variance of the qRT-PCR.

CONCLUSION

Our data indicate that the use of EDTA whole blood tubes for GE-based H-ARS severity prediction is comparable to the quality of PAXgene tubes, when processed ≤ 4 h after venipuncture and the sample is transported within 24 hours at 4 °C.

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.

Radiation-Induced Gene Expression Changes Used for Biodosimetry and Clinical Outcome Prediction: Challenges and Promises

As the war in Ukraine progresses, the radiological and nuclear threat has never been as real as now. The formation of life-threatening acute radiation syndrome (ARS), in particular after the deployment of a nuclear weapon or an attack on a nuclear power station, must be considered realistic. ARS is caused by massive cell death, leading to functional organ deficits and, via systemic inflammatory responses, finally aggravates into multiple organ failure. As a deterministic effect, the severity of the disease dictates the clinical outcome. Hence, predicting ARS severity via biodosimetry or alternative approaches appears straightforward. Because the disease occurs delayed, therapy starting as early as possible has the most significant benefit. A clinically relevant diagnosis should be carried out within the diagnostic time window of about 3 days after exposure. Biodosimetry assays providing retrospective dose estimations within this time frame will support medical management decision-making. However, how closely can dose estimates be associated with the later developing ARS severity degrees when considering dose as one among other determinants of radiation exposure and cell death? From a clinical/triage point of view, ARS severity degrees can be further aggregated into unexposed, weakly diseased (no acute health effects expected), and strongly diseased patient groups, with the latter requiring hospitalization as well as an early and intensive treatment. Radiation-induced gene expression (GE) changes occur early after exposure and can be quickly quantified. GE can be used for biodosimetry purposes. Can GE be used to predict later developing ARS severity degrees and allocate individuals to the three clinically relevant groups as well?

Ex-vivo dose response characterization of the recently identified EDA2R gene after low level radiation exposures and comparison with FDXR gene expression and the γH2AX focus assay

OBJECTIVE

Recently, promising radiation-induced EDA2R gene expression (GE) changes after low level radiation could be shown. Stimulated by that, in this study, we intended to independently validate these findings and to further characterize dose-response relationships in comparison to FDXR and the γH2AX-DNA double-strand break (DSB) focus assay, since both assays are already widely used for biodosimetry purposes.

MATERIALS AND METHODS

Peripheral blood samples from six healthy human donors were irradiated ex vivo (dose: ranging from 2.6 to 49.7 mGy). Subsequently, the fold-differences relative to the sham irradiated reference group were calculated. Radiation-induced changes in GE of FDXR and EDA2R were examined using the quantitative real-time polymerase-chain-reaction (qRT-PCR). DSB foci were quantified in 100 γH2AX + 53BP1 immunostained cells employing fluorescence microscopy. Examinations were performed at single time points enabling sufficient detection of both endpoints.

RESULTS

A significant increase in EDA2R GE relative to the unexposed control was observed in the range of 2.6 mGy (1.6-fold, p = .045) to 5.4 mGy (2.2-fold, p = .0002), whereas the copy numbers increased linearly up to 13.1-fold at 49.7 mGy. On the contrary, FDXR upregulation (2.2-fold) became significant after a 22.6 mGy exposure (p ≤ .02) and increased linearly up to 4-fold at 49.7 mGy. A significant increase in radiation-induced foci (relative to unexposed, RIF-fd) was observed after 11.3 mGy (RIF-fd: 1.5 ± 0.5, p ≤ .03), while the foci increased linearly up to 3-fold at 49.7 mGy. From this, the FDXR and RIF-fd slopes have shown comparability, while the EDA2R slope was five times higher. Nevertheless, the coefficient of variation (CV) of EDA2R was about 30% higher than for RIF-fd.

CONCLUSION

Higher radiation-induced EDA2R GE changes and a lower radiation detection level compared to RIF-fd and FDXR GE changes examined under optimal conditions ex vivo on human samples appear promising. Yet, our results represent just the beginning of further studies to be conducted in animal models for further time- and dose-dependent evaluation and additional examinations on radiologically examined patients to evaluate the impact of confounder, such as age, sex, social behavior, or diseases.

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.

Four Genes Predictive for the Severity of Hematological Damage Reveal a Similar Response after X Irradiation and Chemotherapy

Radiological and especially nuclear accidents and incidents pose a threat to populations. In such events, gene expression (GE) analysis of a set of 4 genes (FDXR, DDB2, POU2AF1, WNT3) is an emerging approach for early and high-throughput prediction of the later manifesting severity degrees of the hematological acute radiation syndrome (H-ARS). Validation of this gene set on radiation victims is difficult since these events are rare. However, chemotherapy (CTX) is widely used e.g., breast cancer patient treatment and pathomechanisms, as well as blood cell count changes are comparable among both exposure types. We wondered whether GE changes are similarly deregulated after CTX, which would be interpreted as a confirmation of our already identified gene set for H-ARS prediction after irradiation. We examined radiation-induced differential GE (DGE) of our gene set as a positive control using in vitro whole blood samples from ten healthy donors (6 females, 4 males, aged: 24-40 years). Blood was incubated in vitro for 8 h after X irradiation with 0 and 4 Gy (1 Gy/min). These data were compared with DGE measured in vivo in blood samples of 10 breast tumor CTX patients (10 females, aged: 39-71 years) before and 4 days after administration of cyclophosphamide and epirubicin. RNA was isolated, reverse transcribed and quantitative real-time polymerase-chain-reaction (qRT-PCR) was performed to assess DGE of FDXR, DDB2, POU2AF1 and WNT3 relative to the unexposed samples using TaqMan assays. After X irradiation, we found a significant upregulation (irrespective of sex) with mean fold changes of 21 (P < 0.001) and 7 (P < 0.001) for FDXR and DDB2 and a significant down-regulation with mean fold changes of 2.5 (P < 0.001) and 2 (P = 0.005) for POU2AF1 and WNT3, respectively. After CTX, a similar pattern was observed, although mean fold changes of up-regulated FDXR (6-fold, P < 0.001) and DDB2 (3-fold, P < 0.001) as well as down-regulated POU2AF1 (1.2-fold, P = 0.270) and WNT3 (1.3-fold, P = 0.069) appeared lower corresponding to less altered blood cell count changes observed after CTX compared to historic radiation exposure data. However, a subpopulation of CTX patients (n = 6) showed on average a significant downregulation of POU2AF1 (1.8-fold, P = 0.04) and WNT3 (2.1-fold, P = 0.008). In summary, the pattern of up-regulated GE changes observed in all CTX patients and down-regulated GE changes observed in a subgroup of CTX patients appeared comparable with an already identified gene set predictive for the radiation-induced H-ARS. This underlines the significance of in vivo GE measurements in CTX patients, employed as a surrogate model to further validate already identified radiation-induced GE changes predictive for the H-ARS.

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.

Identifying radiation responsive exon-regions of genes often used for biodosimetry and acute radiation syndrome prediction

Gene expression (GE) analysis of FDXR, DDB2, WNT3 and POU2AF1 is a promising approach for identification of clinically relevant groups (unexposed, low- and high exposed) after radiological/nuclear events. However, results from international biodosimetry exercises have shown differences in dose estimates based on radiation-induced GE of the four genes. Also, differences in GE using next-generation-sequening (NGS) and validation with quantitative real-time polymerase chain reaction (qRT-PCR) was reported. These discrepancies could be caused by radiation-responsive differences among exons of the same gene. We performed GE analysis with qRT-PCR using TaqMan-assays covering all exon-regions of FDXR, DDB2, WNT3 and POU2AF1. Peripheral whole blood from three healthy donors was X-irradiated with 0, 0.5 and 4 Gy. After 24 and 48 h a dose-dependent up-regulation across almost all exon-regions for FDXR and DDB2 (4–42-fold) was found. A down-regulation for POU2AF1 (two- to threefold) and WNT3 (< sevenfold) at the 3’-end was found at 4 Gy irradiation only. Hence, this confirms our hypothesis for radiation-responsive exon-regions for WNT3 and POU2AF1, but not for FDXR and DDB2. Finally, we identified the most promising TaqMan-assays for FDXR (e.g. AR7DTG3, Hs00244586_m1), DDB2 (AR47X6H, Hs03044951_m1), WNT3 (Hs00902258_m1, Hs00902257_m1) and POU2AF1 (Hs01573370_g1, Hs01573371_m1) for biodosimetry purposes and acute radiation syndrome prediction, considering several criteria (detection limit, dose dependency, time persistency, inter-individual variability).

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.

mRNA and small RNA gene expression changes in peripheral blood to detect internal Ra-223 exposure

PURPOSE

Excretion analysis is the established method for detection of incorporated alpha-emitting radionuclides, but it is laborious and time consuming. We sought a simplified method in which changes in gene expression might be measured in human peripheral blood to detect incorporated radionuclides. Such an approach could be used to quickly determine internal exposure in instances of a radiological dispersal device or a radiation accident.

MATERIALS AND METHODS

We evaluated whole blood samples from five patients with castration-resistant prostate cancer and multiple bone metastases (without visceral or nodal involvement), who underwent treatment with the alpha emitting isotope Radium-223 dichloride (Ra-223, Xofigo^®). Patients received about 4 MBq per cycle and, depending on survival and treatment tolerance, were followed for six months. We collected 24 blood samples approximately monthly corresponding to treatment cycle.

RESULTS

Firstly, we conducted whole genome screening of mRNAs (mRNA seq) and small RNAs (small RNA seq) using next generation sequencing in one patient at eight different time points during all six cycles of Ra-223-therapy. We identified 1900 mRNAs and 972 small RNAs (222 miRNAs) that were differentially up- or down-regulated during follow-up after the first treatment with Ra-223. Overall candidate RNA species inclusion criteria were a general (≥|2|-fold) change or with peaking profiles (≥|5|-fold) at specific points in time. Next we chose 72 candidate mRNAs and 101 small RNAs (comprising 29 miRNAs) for methodologic (n = 8 samples, one patient) and independent (n = 16 samples, four patients) validation by qRT-PCR. In total, 15 mRNAs (but no small RNAs) were validated by methodologic and independent testing. However, the deregulation occurred at different time points, showing a large inter-individual variability in response among patients.

CONCLUSIONS

This proof of concept provides support for the applicability of gene expression measurements to detect internalized alpha-emitting radionuclides, but further work is needed with a larger sample size. While our approach has merit for internal deposition monitoring, it was complicated by the severe clinical condition of the patients we studied.

Gene expression for biodosimetry and effect prediction purposes: promises, pitfalls and future directions - key session ConRad 2021

PURPOSE

In a nuclear or radiological event, an early diagnostic or prognostic tool is needed to distinguish unexposed from low- and highly exposed individuals with the latter requiring early and intensive medical care. Radiation-induced gene expression (GE) changes observed within hours and days after irradiation have shown potential to serve as biomarkers for either dose reconstruction (retrospective dosimetry) or the prediction of consecutively occurring acute or chronic health effects. The advantage of GE markers lies in their capability for early (1-3 days after irradiation), high-throughput, and point-of-care (POC) diagnosis required for the prediction of the acute radiation syndrome (ARS).

CONCLUSIONS

As a key session of the ConRad conference in 2021, experts from different institutions were invited to provide state-of-the-art information on a range of topics including: (1) Biodosimetry: What are the current efforts to enhance the applicability of this method to perform retrospective biodosimetry? (2) Effect prediction: Can we apply radiation-induced GE changes for prediction of acute health effects as an approach, complementary to and integrating retrospective dose estimation? (3) High-throughput and point-of-care diagnostics: What are the current developments to make the GE approach applicable as a high-throughput as well as a POC diagnostic platform? (4) Low level radiation: What is the lowest dose range where GE can be used for biodosimetry purposes? (5) Methodological considerations: Different aspects of radiation-induced GE related to more detailed analysis of exons, transcripts and next-generation sequencing (NGS) were reported.

Gene expression changes in male and female rhesus macaque 60 days after irradiation

PURPOSE

Transcriptome changes can be expected in survivors after lethal irradiation. We aimed to characterize these in males and females and after different cytokine treatments 60 days after irradiation.

MATERIAL AND METHODS

Male and female rhesus macaques (n = 142) received a whole-body exposure with 700 cGy, from which 60 animals survived. Peripheral whole blood was drawn pre-exposure and before sacrificing the surviving animals after 60 days.

RESULTS

We evaluated gene expression in a three-phase study design. Phase I was a whole-genome screening (NGS) for mRNAs using five pre- and post-exposure RNA samples from both sexes (n = 20). Differential gene expression (DGE) was calculated between samples of survivors and pre-exposure samples (reference), separately for males and females. 1,243 up- and down-regulated genes were identified with 30-50% more deregulated genes in females. 37 candidate mRNAs were chosen for qRT-PCR validation in phase II using the remaining samples (n = 117). Altogether 17 genes showed (borderline) significant (t-test) DGE in groups of untreated or treated animals. Nine genes (CD248, EDAR, FAM19A5, GAL3ST4, GCNT4, HBG2/1, LRRN1, NOG, SYT14) remained with significant changes and were detected in at least 50% of samples per group. Panther analysis revealed an overlap between both sexes, related to the WNT signaling pathway, cell adhesion and immunological functions. For phase III, we validated the nine genes with candidate genes (n = 32) from an earlier conducted study on male baboons. Altogether 14 out of 41 genes showed a concordantly DGE across both species in a bilateral comparison.

CONCLUSIONS

Sixty days after radiation exposure, we identified (1) sex and cytokine treatment independent transcriptional changes, (2) females with almost twice as much deregulated genes appeared more radio-responsive than males, (3) Panther analysis revealed an association with immunological processes and WNT pathway for both sexes.

Gene expression changes and DNA damage after ex vivo exposure of peripheral blood cells to various CT photon spectra

Dual-energy CT provides enhanced diagnostic power with similar or even reduced radiation dose as compared to single-energy CT. Its principle is based on the distinct physical properties of low and high energetic photons, which, however, may also affect the biological effectiveness and hence the extent of CT-induced cellular damage. Therefore, a comparative analysis of biological effectiveness of dual- and single-energy CT scans with focus on early gene regulation and frequency of radiation-induced DNA double strand breaks (DSBs) was performed. Blood samples from three healthy individuals were irradiated ex vivo with single-energy (80 kV and 150 kV) and dual-energy tube voltages (80 kV/Sn150kV) employing a modern dual source CT scanner resulting in Volume Computed Tomography Dose Index (CTDIvol) of 15.79-18.26 mGy and dose length product (DLP) of 606.7-613.8 mGy*cm. Non-irradiated samples served as a control. Differential gene expression in peripheral blood mononuclear cells was analyzed 6 h after irradiation using whole transcriptome sequencing. DSB frequency was studied by 53BP1 + γH2AX co-immunostaining and microscopic evaluation of their focal accumulation at DSBs. Neither the analysis of gene expression nor DSB frequency provided any evidence for significantly increased biological effectiveness of dual-energy CT in comparison to samples irradiated with particular single-energy CT spectra. Relative to control, irradiated samples were characterized by a significantly higher rate of DSBs (p < 0.001) and the shared upregulation of five genes, AEN, BAX, DDB2, FDXR and EDA2R, which have already been suggested as radiation-induced biomarkers in previous studies. Despite steadily decreasing doses, CT diagnostics remain a genotoxic stressor with impact on gene regulation and DNA integrity. However, no evidence was found that varying X-ray spectra of CT impact the extent of cellular damage.

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.

Identifying a Diagnostic Window for the Use of Gene Expression Profiling to Predict Acute Radiation Syndrome

In the event of a mass casualty radiological or nuclear scenario, it is important to distinguish between the unexposed (worried well), low-dose exposed individuals and those developing the hematological acute radiation syndrome (HARS) within the first three days postirradiation. In previous baboon studies, we identified altered gene expression changes after irradiation, which were predictive for the later developing HARS severity. Similar changes in the expression of four of these genes were observed using an in vitro human whole blood model. However, these studies have provided only limited information on the time frame of the changes after exposure in relationship to the development of HARS. In this study we analyzed the time-dependent changes in mRNA expression after in vitro irradiation of whole blood. Changes in the expression of informative mRNAs (FDXR, DBB2, POU2AF1 and WNT3) were determined in the blood of eight healthy donors (6 males, 2 females) after irradiation at 0 (control), 0.5, 2 and 4 Gy using qRT-PCR. FDXR expression was significantly upregulated (P < 0.001) 4 h after ≥0.5 Gy irradiation, with an 18-40-fold peak attained 4-12 h postirradiation which remained elevated (4-9-fold) at 72 h. DDB2 expression was upregulated after 4 h (fold change, 5-8, P < 0.001 at ≥ 0.5 Gy) and remained upregulated (3-4-fold) until 72 h (P < 0.001). The earliest time points showing a significant downregulation of POU2AF1 and WNT3 were 4 h (fold change = 0.4, P = 0.001, at 4 Gy) and 8 h (fold change = 0.3-0.5, P < 0.001, 2-4 Gy), respectively. These results indicate that the diagnostic window for detecting HARS-predictive changes in gene expression may be opened as early as 2 h for most (75%) and at 4 h postirradiation for all individuals examined. Depending on the RNA species studied this may continue for at least three days postirradiation.

Software Tools for the Evaluation of Clinical Signs and Symptoms in the Medical Management of Acute Radiation Syndrome-A Five-year Experience

A suite of software tools has been developed for dose estimation (BAT, WinFRAT) and prediction of acute health effects (WinFRAT, H-Module) using clinical symptoms and/or changes in blood cell counts. We constructed a database of 191 ARS cases using the METREPOL (n = 167) and the SEARCH-database (n = 24). The cases ranged from unexposed (RC0), to mild (RC1), moderate (RC2), severe (RC3), and lethal ARS (RC4). From 2015-2019, radiobiology students and participants of two NATO meetings predicted clinical outcomes (RC, H-ARS, and hospitalization) based on clinical symptoms. We evaluated the prediction outcomes using the same input datasets with a total of 32 teams and 94 participants. We found that: (1) unexposed (RC0) and mildly exposed individuals (RC1) could not be discriminated; (2) the severity of RC2 and RC3 were systematically overestimated, but almost all lethal cases (RC4) were correctly predicted; (3) introducing a prior education component for non-physicians significantly increased the correct predictions of RC, ARS, and hospitalization by around 10% (p<0.005) with a threefold reduction in variance and a halving of the evaluation time per case; (4) correct outcome prediction was independent of the software tools used; and (5) comparing the dose estimates generated by the teams with H-ARS severity reflected known limitations of dose alone as a surrogate for H-ARS severity. We found inexperienced personnel can use software tools to make accurate diagnostic and treatment recommendations with up to 98% accuracy. Educational training improved the quality of decision making and enabled participants lacking a medical background to perform comparably to experts.

Acute radiation syndrome-related gene expression in irradiated peripheral blood cell populations

PURPOSE

In a nuclear or radiological event, an early diagnostic tool is needed to distinguish the worried well from those individuals who may later develop life-threatenFing hematologic acute radiation syndrome. We examined the contribution of the peripheral blood's cell populations on radiation-induced gene expression (GE) changes.

MATERIALS AND METHODS

EDTA-whole-blood from six healthy donors was X-irradiated with 0 and 4Gy and T-lymphocytes, B-lymphocytes, NK-cells and granulocytes were separated using immunomagnetic methods. GE were examined in cell populations and whole blood.

RESULTS

The cell populations contributed to the total RNA amount with a ratio of 11.6 for T-lymphocytes, 1.2 for B-cells, 1.2 for NK-cells, 1.0 for granulocytes. To estimate the contribution of GE per cell population, the baseline (0Gy) and the radiation-induced fold-change in GE relative to unexposed was considered for each gene. The T-lymphocytes (74.8%/80.5%) contributed predominantly to the radiation-induced up-regulation observed for FDXR/DDB2 and the B-lymphocytes (97.1%/83.8%) for down-regulated POU2AF1/WNT3 with a similar effect on whole blood gene expression measurements reflecting a corresponding order of magnitude.

CONCLUSIONS

T- and B-lymphocytes contributed predominantly to the radiation-induced up-regulation of FDXR/DDB2 and down-regulation of POU2AF1/WNT3. This study underlines the use of FDXR/DDB2 for biodosimetry purposes and POU2AF1/WNT3 for effect prediction of acute health effects.

Preparing for a "dirty bomb" attack: the optimum mix of medical countermeasure resources

BACKGROUND

In radiological emergencies with radionuclide incorporation, decorporation treatment is particularly effective if started early. Treating all people potentially contaminated ("urgent treatment") may require large antidote stockpiles. An efficacious way to reduce antidote requirements is by using radioactivity screening equipment. We analyzed the suitability of such equipment for triage purposes and determined the most efficient mix of screening units and antidote daily doses.

METHODS

The committed effective doses corresponding to activities within the detection limits of monitoring portals and mobile whole-body counters were used to assess their usefulness as triage tools. To determine the optimal resource mix, we departed from a large-scale scenario (60,000 victims) and based on purchase prices of antidotes and screening equipment in Germany, we calculated efficiencies of different combinations of medical countermeasure resources by data envelopment analysis. Cost-effectiveness was expressed as the costs per life year saved and compared to risk reduction opportunities in other sectors of society as well as the values of a statistical life.

RESULTS

Monitoring portals are adequate instruments for a sensitive triage after cesium-137 exposure with a high screening throughput. For the detection of americium-241 whole-body counters with a lower daily screening capacity per unit are needed. Assuming that 1% of the potentially contaminated patients actually need decorporation treatment, an efficient resource mix includes 6 monitoring portals and 25 mobile whole-body counters. The optimum mix depends on price discounts and in particular the fraction of victims actually needing treatment. The cost-effectiveness of preparedness for a "dirty bomb" attack is less than for common health care, but costs for a life year saved are less than for many risk-reduction interventions in the environmental sector.

CONCLUSION

To achieve economic efficiency a high daily screening capacity is of major importance to substantially decrease the required amount of antidote doses. Among the determinants of the number of equipment units needed, the fraction of the potentially contaminated victims that actually needs treatment is the most difficult to assess. Judging cost-effectiveness of the preparedness for "dirty bomb" attacks is an issue of principle that must be dealt with by political leaders.

USING DICENTRIC DOSE ESTIMATES AND EARLY RADIATION-INDUCED BLOOD CELL COUNT CHANGES OF REAL CASE HISTORIES FOR VALIDATION OF THE HEMODOSE BIODOSIMETRY TOOL

HemoDose is a software tool, which estimates absorbed doses based on blood cell counts (BCC). The aim of our study was to validate HemoDose for early dose estimates. Dose estimates generated by HemoDose were compared with dose estimates stored in SEARCH from radiation victims. Moreover, BCC from unirradiated donors and corresponding HemoDose dose estimates were analysed. We employed linear or logistic regression analysis. There was a significant correlation between calculated doses by HemoDose based on single and multiple lymphocyte counts when omitting lowest and highest dose estimates. Furthermore, there was a significant correlation between calculated doses by HemoDose based on lymphocyte counts and the estimated doses based on DIC. The dose estimates by HemoDose based on lymphocyte counts and DIC showed a comparable correlation with HARS degrees 0 and 4. In conclusion, HemoDose dose estimation based on early lymphocyte counts appears to be a promising biodosimetry tool under certain considerations.

A New Smartphone Application to Predict Hematologic Acute Radiation Syndrome Based on Blood Cell Count Changes-The H-module App

Treatment regimens for acute radiation syndrome have been improved over the past years. The application of appropriate therapy relies on rapid and high-throughput tests ideally conducted in the first 3 d after a radiation exposure event. We have examined the utility of blood cell counts (BCCs) 3 d post irradiation to predict clinical outcome for hematologic acute radiation syndrome (HARS). The BCCs and HARS severity information originated from data available in the System-for-Evaluation-and-Archiving-of-Radiation Accidents-based-on-Case-Histories (SEARCH). We found an almost complete discrimination of unexposed (HARS score H0) vs. irradiated individuals during model development and validation (negative predictive value > 94%) when using BCC data for all 3 d. We also found that BCC data increased the correct prediction of exposed individuals from day 1 to day 3. We developed spreadsheets to calculate the likelihood of correct diagnoses of the worried-well, requirement of hospitalization (HARS 2-4), or development of severe hematopoietic syndrome (HARS 3-4). In two table-top exercises, we found the spreadsheets were confusing and cumbersome, so we converted the spreadsheets into a smartphone application, named the H-module App, designed for ease of use, wider dissemination, and accommodation of co-morbidities in the HARS severity prediction algorithm.

Characterization of Primary Human Dermal Fibroblasts to Ensure for Instance EMF Exposure Experiments under Comparable Cell Culture Condition

HDFa (human dermal fibroblasts) are used as cellular models for EMF exposure. To ensure reproducible in vitro experiments, comparable proliferation and differentiation cell conditions must exist, and different donors, passage numbers, culture time, and growth media must be considered. In this study, the authors cultured fibroblasts in DMEM or 106 medium. Growth curves, vitality, morphology, and gene expression of genes coding for proliferation (PCNA, CDKN2A, CDKN1A, SFN), differentiation (PDGFRA, TGM2, ACTA2, PDPN, NTN1, MGP, PPP1R14), and SFN target genes (TP63, MMP1, MMP3) were examined in both media and passage numbers 3-4, 5-6 and >6. At passages 3-4, proliferating cells can be observed in both media. While cells cultured in DMEM proliferate over the passages, from passage 5, cells in 106 medium persisted around the seeded number. TGM2 down-regulation over all passages in both media and cells morphology suggest papillary-type fibroblasts. Downregulation of SFN (negative regulator of mitotic translation and cell differentiation) coincided with proliferating fibroblasts over all examined conditions. Downstream SFN target genes in proliferating cells appeared upregulated (TP63) and downregulated (MMP1/MMP3), suggestive for a status characterized by increased stemnesses (upregulated TP63) and wound healing capacity (downregulated MMP1, MMP3). Resting cells (SFN control values) were associated with control values of TP63 and MMP1/MMP3 expression, suggesting a reduced stemness and wound healing capacity. In conclusion, a set of markers related to proliferation (SFN), differentiation (TGM2), stemnesses (TP63), and wound healing (MMP1/MMP3) allow a culture characterization so that cells under two different conditions can be exposed, thus enabling reproducible EMF experiments or experiments with other exposures.

Utility of gene expression studies in relation to radiation exposure and clinical outcomes: thyroid cancer in the Ukrainian-American cohort and late health effects in a MAYAK worker cohort

PURPOSE

We herein report on changes in gene expression after radiation exposure to iodine-131 from the Chernobyl accident in the Ukrainian-American thyroid cohort and to external gamma ray or internal plutonium exposure in the Mayak Production Association radiation workers.

MATERIALS AND METHODS

Taking advantage of access to tissue samples from the thyroid cancer cases in the Ukrainian-American cohort, our group tried to identify candidate genes to discriminate spontaneously occurring thyroid cancers from thyroid cancers caused by radiation exposure. We also examined gene expression changes in normal and cancerous thyroid tissue in relation to iodine-131 dose separately. Gene expression changes in the peripheral blood of radiation exposed Mayak workers were examined to elucidate the dose-to-gene and gene-to-health (e.g. cardiovascular disease) relationships.

CONCLUSIONS

Results of both projects are discussed under the aspect of dose-response relationships (dose-to-gene) and clinical outcome relationships (gene-to-effect) in light of how mechanistic data can be translated into actionable knowledge for radiation protection or clinical purposes.

Rapid High-Throughput Diagnostic Triage after a Mass Radiation Exposure Event Using Early Gene Expression Changes

Radiological exposure scenarios involving large numbers of people require a rapid and high-throughput method to identify the unexposed, and those exposed to low- and high-dose radiation. Those with high-dose exposure, e.g., >2 Gy and depending on host characteristics, may develop severe hematological acute radiation syndrome (HARS), requiring hospitalization and treatment. Previously, we identified a set of genes that discriminated these clinically relevant groups. In the current work, we examined the utility of gene expression changes to classify 1,000 split blood samples into HARS severity scores of H0, H1 and H2-4, with the latter indicating likely hospitalization. In several previous radiation dose experiments, we determined that these HARS categories corresponded, respectively, to doses of 0 Gy (unexposed), 0.5 Gy and 5 Gy. The main purpose of this work was to assess the rapidity of blood sample processing using targeted next-generation sequencing (NGS). Peripheral blood samples from two healthy donors were X-ray irradiated in vitro and incubated at 37°C for 24 h. A total of 1,000 samples were evaluated by laboratory personnel blinded to the radiation dose. Changes in gene expression of FDXR, DDB2, POU2AF1 and WNT3 were examined with qRT-PCR as positive controls. Targeted NGS (TREX) was used on all samples for the same four genes. Agreement using both methods was almost 78%. Using NGS, all 1,000 samples were processed within 30 h. Classification of the HARS severity categories corresponding to radiation dose had an overall agreement ranging between 90-97%. Depending on the end point, either a combination of all genes or FDXR alone (H0 HARS or unexposed) provided the best classification. Using this optimized automated methodology, we assessed 100× more samples approximately three times faster compared to standard cytogenetic studies. We showed that a small set of genes, rather than a complex constellation of genes, provided robust positive (97%) and negative (97%) predictive values for HARS categories and radiation doses of 0, 0.5 and 5 Gy. The findings of this study support the potential utility of early radiation-induced gene expression changes for high-throughput biodosimetry and rapid identification of irradiated persons in need of hospitalization.

Persistent mRNA and miRNA expression changes in irradiated baboons

We examined the transcriptome/post-transcriptome for persistent gene expression changes after radiation exposure in a baboon model. Eighteen baboons were irradiated with a whole body equivalent dose of 2.5 or 5 Gy. Blood samples were taken before, 7, 28 and 75–106 days after radiation exposure. Stage I was a whole genome screening for mRNA combined with a qRT-PCR platform for detection of 667 miRNAs. Candidate mRNAs and miRNAs differentially up- or down-regulated in stage I were chosen for validation in stage II using the remaining samples. Only 12 of 32 candidate genes provided analyzable results with two mRNAs showing significant 3–5-fold differences in gene expression over the reference (p < 0.0001). From 667 candidate miRNAs, 290 miRNA were eligible for analysis with 21 miRNAs independently validated using qRT-PCR. These miRNAs showed persistent expression changes on each day and over days 7–106 days after exposure (n = 7). In particular miR-212 involved in radiosensitivity and immune modulation appeared persistently and 48–77-fold up-regulated over the entire time period. We are finally trying to put our results into a context of clinical implications and provide possible hints on underlying molecular mechanisms to be examined in future studies.

Drug Coated Balloon Is Less Effective for Treatment of DES In-Stent Restenosis Both in Native Coronary Arteries and Saphenous Vein Grafts: Results From a Bicenter Registry

BACKGROUND

The paclitaxel drug coated balloon (DCB) is an established treatment for bare metal stent (BMS) in-stent restenosis (ISR) in native coronary arteries. The evidence of DCB-application for drug eluting stent (DES) ISR both in native coronaries and saphenous vein grafts (SVG) is limited. Aim of our study was to compare the differential efficacy of DCB for treatment of BMS- and DES-ISR in native coronary vessels and SVGs.

METHODS AND RESULTS

N = 135 DCB-treated patients with available follow up (FU) angiography were included in this retrospective study. Patients received treatment between April 2009 and March 2013 at 2 tertiary care hospitals in Germany. DCB was applied in BMS-ISR (n = 65; 48%) and DES-ISR (n = 70; 52%). DCB-treated lesions were located in native coronary arteries (n = 110; 81%; BMS-ISR: n = 58; 53%; DES-ISR: n = 52; 47%) and SVGs (n = 25; 19%; BMS-ISR: n = 7, 28%; DES-ISR: n = 18, 72%). Median FU was 12 months. Endpoints were binary restenosis and target lesion revascularization (TLR). Binary restenosis (29% vs. 57%; P < 0.01) and TLR (18% vs. 46%; P < 0.01) were significantly more frequent in DES-ISR versus BMS-ISR. In SVGs, TLR was required in 72% (DES-ISR) versus 14% (BMS-ISR); P = 0.02. In the Kaplan-Meier-analysis freedom from both endpoints was significantly decreased in the DES-lesions both in the total population (binary restenosis P < 0.01; TLR P < 0.01) and native coronaries (binary restenosis P = 0.02; TLR P = 0.04).

CONCLUSIONS

DCB treatment is less effective in DES-ISR than in BMS-ISR. The diminished efficacy of DCB treatment is even more pronounced in DES-ISR located within degenerated SVGs.