Gene Expression Changes in Irradiated Baboons: A Summary and Interpretation of a Decade of Findings

Radiation Research Society Journal-based Historical Review of the Use of Biomarkers for Radiation Dose and Injury Assessment: Acute Health Effects Predictions

A multiple-parameter based approach using radiation-induced clinical signs and symptoms, hematology changes, cytogenetic chromosomal aberrations, and molecular biomarkers changes after radiation exposure is used for biodosimetry-based dose assessment. In the current article, relevant milestones from Radiation Research are documented that forms the basis of the current consensus approach for diagnostics after radiation exposure. For example, in 1962 the use of cytogenetic chromosomal aberration using the lymphocyte metaphase spread dicentric assay for biodosimetry applications was first published in Radiation Research. This assay is now complimented using other cytogenetic chromosomal aberration assays (i.e., chromosomal translocations, cytokinesis-blocked micronuclei, premature chromosome condensation, γ-H2AX foci, etc.). Changes in blood cell counts represent an early-phase biomarker for radiation exposures. Molecular biomarker changes have evolved to include panels of organ-specific plasma proteomic and blood-based gene expression biomarkers for radiation dose assessment. Maturation of these assays are shown by efforts for automated processing and scoring, development of point-of-care diagnostics devices, service laboratories inter-comparison exercises, and applications for dose and injury assessments in radiation accidents. An alternative and complementary approach has been advocated with the focus to de-emphasize "dose" and instead focus on predicting acute or delayed health effects. The same biomarkers used for dose estimation (e.g., lymphocyte counts) can be used to directly predict the later developing severity degree of acute health effects without performing dose estimation as an additional or intermediate step. This review illustrates contributing steps toward these developments published in Radiation Research.

Neutrophil to lymphocyte ratio in captive olive baboons (Papio anubis): The effects of age, sex, rearing, stress, and pregnancy

In apes and humans, neutrophil to lymphocyte ratio (NLR) can be used as a predictive indicator of a variety of clinical conditions, longevity, and physiological stress. In chimpanzees specifically, NLR systematically varies with age, rearing, sex, and premature death, indicating that NLR may be a useful diagnostic tool in assessing primate health. To date, just one very recent study has investigated NLR in old world monkeys and found lower NLR in males and nursery-reared individuals, as well as a negative relationship between NLR and disease outcomes. Given that baboons are increasingly used as research models, we aimed to characterize NLR in baboons by providing descriptive data and examinations of baboon NLR heritability, and of the relationships between NLR, age, rearing, and sex in 387 olive baboons (Papio anubis) between 6 months and 19 years of age. We found that (1) mother-reared baboons had higher NLRs than nursery-reared baboons; (2) females had higher NLRs than males; and (3) there was a quadratic relationship between NLR and age, such that middle-aged individuals had the highest NLR values. We also examined NLR as a function of transport to a new facility using a subset of the data. Baboons exhibited significantly higher transport NLRs compared to routine exam NLRs. More specifically, adult baboons had higher transport NLRs than routine NLRs, whereas juveniles showed no such difference, suggesting that younger animals may experience transport stress differently than older animals. We also found that transport NLR was heritable, whereas routine NLR was not, possibly suggesting that stress responses (as indicated in NLR) have a strong genetic component. Consistent with research in humans and chimpanzees, these findings suggest that NLR varies with important biological and life history variables and that NLR may be a useful health biomarker in baboons.

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.

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.

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.

Effects of cyclophosphamide and mitomycin C on radiation-induced transcriptional biomarkers in human lymphoblastoid cells

PURPOSE

Ionizing radiation (IR)-induced transcriptional changes are considered a potential biodosimetry for dose evaluation and health risk monitoring of acute or chronic radiation exposure. It is crucial to understand the impact of confounding factors on the radiation-responsive gene expressions for accurate and reproducible dose assessment. This study aims to explore the potential influence of exposures to chemotherapeutic agents such as cyclophosphamide (CP) and mitomycin C (MMC) on IR-induced transcriptional biomarkers.

METHODS

The human B lymphoblastoid cells (AHH-1) were exposed to 0, 20, 50, 100, 200 and 500 μg/ml CP or 0, 0.025, 0.05, 0.1 and 1 μg/ml MMC, respectively. The appropriate concentrations of CP and MMC were added for 1 h before irradiation with 0, 2, 4 and 6 Gy of 60Co γ-rays at a dose rate of 1 Gy/min. Cell viability was evaluated by CCK-8 assay. The gene expression responses of 18 radiation-induced transcriptional biomarkers were examined at 24 h after exposures to CP and MMC, respectively. The expression levels of five crucial DNA interstrand crosslinks (ICLs) repair genes were also evaluated. The biodosimetry models were established based on the specific radiation-responsive gene combinations.

RESULTS

The baseline transcriptional levels of the 18 selected genes were slightly affected by CP treatment in the absence of IR, while the transcript responses to IR could be inhibited as the concentration of CP up to 50 μg/ml. MMC treatment up-regulated the background levels in most radiation-responsive gene expressions. Of 18 genes, only the relative mRNA expression levels of CDKN1A and BBC3 were repressed after treatment with IR and MMC in combination. The relative mRNA level of RAD51 was significantly up-regulated after exposure to CP, while the expression of FANCD2, RAD51 and BLM showed an overall increase in response to MMC treatment. After irradiation, the relative mRNA expression levels of FANCD2, BRCA2 and RAD51 exhibited dose-dependent increases in IR alone and MMC treatment groups. In addition, the biodosimetry models were established using 2-4 radiation-responsive genes based on different radiation exposure scenarios.

CONCLUSION

Our findings suggested that IR-induced gene expression changes were slightly affected after exposure to a relatively low concentration of CP and MMC. Gene expression combinations might improve the broad applicability of transcriptional biodosimetry across diverse radiation exposures.

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.

Transcriptome profile changes in the jejunum of nonhuman primates exposed to supralethal dose of total- or partial-body radiation

The risk of exposure of the general public or military personnel to high levels of ionizing radiation from nuclear weapons or radiological accidents is a dire national security matter. The development of advanced molecular biodosimetry methods, those that measure biological response, such as transcriptomics, to screen large populations of radiation-exposed victims is key to improving survival outcomes during radiological mass casualty scenarios. In this study, nonhuman primates were exposed to either 12.0 Gy cobalt-60 gamma (total-body irradiation, TBI) or X-ray (partial-body irradiation, PBI) 24 h after administration of a potential radiation medical countermeasure, gamma-tocotrienol (GT3). Changes in the jejunal transcriptomic profiles in GT3-treated and irradiated animals were compared to healthy controls to assess the extent of radiation damage. No major effect of GT3 on radiation-induced transcriptome at this radiation dose was identified. About 80% of the pathways with a known activation or repression state were commonly observed between both exposures. Several common pathways activated due to irradiation include FAK signaling, CREB signaling in the neurons, phagosome formation, and G-protein coupled signaling pathway. Sex-specific differences associated with excessive mortality among irradiated females were identified in this study, including Estrogen receptor signaling. Differential pathway activation was also identified across PBI and TBI, pointing towards altered molecular response for different degrees of bone marrow sparing and radiation doses. This study provides insight into radiation-induced changes in jejunal transcriptional profiles, supporting the investigation for the identification of biomarkers for radiation injury and countermeasure efficacy.

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.

Longitudinal multi-omic changes in the transcriptome and proteome of peripheral blood cells after a 4 Gy total body radiation dose to Rhesus macaques

BACKGROUND

Non-human primates, such as Rhesus macaques, are a powerful model for studies of the cellular and physiological effects of radiation, development of radiation biodosimetry, and for understanding the impact of radiation on human health. Here, we study the effects of 4 Gy total body irradiation (TBI) at the molecular level out to 28 days and at the cytogenetic level out to 56 days after exposure. We combine the global transcriptomic and proteomic responses in peripheral whole blood to assess the impact of acute TBI exposure at extended times post irradiation.

RESULTS

The overall mRNA response in the first week reflects a strong inflammatory reaction, infection response with neutrophil and platelet activation. At 1 week, cell cycle arrest and re-entry processes were enriched among mRNA changes, oncogene-induced senescence and MAPK signaling among the proteome changes. Influenza life cycle and infection pathways initiated earlier in mRNA and are reflected among the proteomic changes during the first week. Transcription factor proteins SRC, TGFβ and NFATC2 were immediately induced at 1 day after irradiation with increased transcriptional activity as predicted by mRNA changes persisting up to 1 week. Cell counts revealed a mild / moderate hematopoietic acute radiation syndrome (H-ARS) reaction to irradiation with expected lymphopenia, neutropenia and thrombocytopenia that resolved within 30 days. Measurements of micronuclei per binucleated cell levels in cytokinesis-blocked T-lymphocytes remained high in the range 0.27-0.33 up to 28 days and declined to 0.1 by day 56.

CONCLUSIONS

Overall, we show that the TBI 4 Gy dose in NHPs induces many cellular changes that persist up to 1 month after exposure, consistent with damage, death, and repopulation of blood cells.

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.

Lung transcriptome of nonhuman primates exposed to total- and partial-body irradiation

The focus of radiation biodosimetry has changed recently, and a paradigm shift for using molecular technologies of omic platforms in addition to cytogenetic techniques has been observed. In our study, we have used a nonhuman primate model to investigate the impact of a supralethal dose of 12 Gy radiation on alterations in the lung transcriptome. We used 6 healthy and 32 irradiated animal samples to delineate radiation-induced changes. We also used a medical countermeasure, γ-tocotrienol (GT3), to observe any changes. We demonstrate significant radiation-induced changes in the lung transcriptome for total-body irradiation (TBI) and partial-body irradiation (PBI). However, no major influence of GT3 on radiation was noted in either comparison. Several common signaling pathways, including PI3K/AKT, GADD45, and p53, were upregulated in both exposures. TBI activated DNA-damage-related pathways in the lungs, whereas PTEN signaling was activated after PBI. Our study highlights the various transcriptional profiles associated with γ- and X-ray exposures, and the associated pathways include LXR/RXR activation in TBI, whereas pulmonary wound-healing and pulmonary fibrosis signaling was repressed in PBI. Our study provides important insights into the molecular pathways associated with irradiation that can be further investigated for biomarker discovery.

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).

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.

Impact of radiation therapy on healthy tissues

Radiation therapy has a fundamental role in the management of cancers. However, despite a constant improvement in radiotherapy techniques, the issue of radiation-induced side effects remains clinically relevant. Mechanisms of acute toxicity and late fibrosis are therefore important topics for translational research to improve the quality of life of patients treated with ionizing radiations. Tissue changes observed after radiotherapy are consequences of complex pathophysiology, involving macrophage activation, cytokine cascade, fibrotic changes, vascularization disorders, hypoxia, tissue destruction and subsequent chronic wound healing. Moreover, numerous data show the impact of these changes in the irradiated stroma on the oncogenic process, with interplays between tumor radiation response and pathways involved in the fibrotic process. The mechanisms of radiation-induced normal tissue inflammation are reviewed, with a focus on the impact of the inflammatory process on the onset of treatment-related toxicities and the oncogenic process. Possible targets for pharmacomodulation are also discussed.

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.