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

Tin prefiltration in computed tomography does not significantly alter radiation-induced gene expression and DNA double-strand break formation

BACKGROUND

The tin (Sn) prefilter technique is a recently introduced dose-saving technique in computed tomography (CT). This study investigates whether there is an altered molecular biological response in blood cells using the tin prefiltering technique.

METHODS

Blood from 6 donors was X-irradiated ex-vivo with 20 mGy full dose (FD) protocols (Sn 150 kV, 150 kV, and 120 kV) and a tin prefiltered 16.5 mGy low dose (LD) protocol on a CT scanner. Biological changes were determined by quantification of γH2AX DNA double-strand break (DSB) foci, and differential gene expression (DGE) relative to unexposed samples were examined for seven known radiation-induced genes (FDXR, DDB2, BAX, CDKN1A, AEN, EDA2R, APOBEC3H) and 667 microRNAs (miRNA).

RESULTS

EDA2R and DDB2 gene expression (GE) increased 1.7-6-fold (p = 0.0004-0.02) and average DNA DSB foci value (0.31±0.02, p<0.0001) increased significantly relative to unexposed samples, but similarly for the applied radiation protocols. FDXR upregulation (2.2-fold) was significant for FD protocols (p = 0.01-0.02) relative to unexposed samples. miRNA GE changes were not significant (p = 0.15-1.00) and DGE were similar for the examined protocols (p = 0.10-1.00). An increased frequency of lower DGE values was seen in the Sn 150 kV LD protocol compared to the 120 kV FD and Sn 150 kV FD protocols (p = 0.001-0.008).

CONCLUSIONS

The current ex-vivo study indicates no changes regarding transcriptional and post-transcriptional DGE and DNA DSB induction when using the tin prefilter technique and even a significant tendency to lower radiation-induced DGE-changes due to the dose reduction of the tin prefilter with equal image quality compared to classical CT scan protocols was found.

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.

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

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

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

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

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

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

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

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

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.

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.

The DNA damage response to radiological imaging: from ROS and γH2AX foci induction to gene expression responses in vivo

Candidate ionising radiation exposure biomarkers must be validated in humans exposed in vivo. Blood from patients undergoing positron emission tomography-computed tomography scan (PET-CT) and skeletal scintigraphy (scintigraphy) was drawn before (0 h) and after (2 h) the procedure for correlation analyses of the response of selected biomarkers with radiation dose and other available patient information. FDXR, CDKN1A, BBC3, GADD45A, XPC, and MDM2 expression was determined by qRT-PCR, DNA damage (γH2AX) by flow cytometry, and reactive oxygen species (ROS) levels by flow cytometry using the 2', 7'-dichlorofluorescein diacetate test in peripheral blood mononuclear cells (PBMC). For ROS experiments, 0- and 2-h samples were additionally exposed to UVA to determine whether diagnostic irradiation conditioned the response to further oxidative insult. With some exceptions, radiological imaging induced weak γH2AX foci, ROS and gene expression fold changes, the latter with good coherence across genes within a patient. Diagnostic imaging did not influence oxidative stress in PBMC successively exposed to UVA. Correlation analyses with patient characteristics led to low correlation coefficient values. γH2AX fold change, which correlated positively with gene expression, presented a weak positive correlation with injected activity, indicating a radiation-induced subtle increase in DNA damage and subsequent activation of the DNA damage response pathway. The exposure discrimination potential of these biomarkers in the absence of control samples as frequently demanded in radiological emergencies, was assessed using raw data. These results suggest that the variability of the response in heterogeneous populations might complicate identifying individuals exposed to low radiation doses.

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.

The role of Evi/Wntless in exporting Wnt proteins

Intercellular communication by Wnt proteins governs many essential processes during development, tissue homeostasis and disease in all metazoans. Many context-dependent effects are initiated in the Wnt-producing cells and depend on the export of lipidated Wnt proteins. Although much focus has been on understanding intracellular Wnt signal transduction, the cellular machinery responsible for Wnt secretion became better understood only recently. After lipid modification by the acyl-transferase Porcupine, Wnt proteins bind their dedicated cargo protein Evi/Wntless for transport and secretion. Evi/Wntless and Porcupine are conserved transmembrane proteins, and their 3D structures were recently determined. In this Review, we summarise studies and structural data highlighting how Wnts are transported from the ER to the plasma membrane, and the role of SNX3-retromer during the recycling of its cargo receptor Evi/Wntless. We also describe the regulation of Wnt export through a post-translational mechanism and review the importance of Wnt secretion for organ development and cancer, and as a future biomarker.

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

Transcriptional Dynamics of DNA Damage Responsive Genes in Circulating Leukocytes during Radiotherapy

Simple Summary

In this study, the transcriptional response of a panel of radiation responsive genes was monitored over time in blood samples after radiation exposure in vivo. For this aim, cancer patients treated by radiotherapy were recruited after consent forms were obtained. Following the first fraction of radiotherapy, 2 mL blood samples were collected at different time points during the first 24h hours (before the second fraction was delivered) and at mid and end of treatment. Amongst the 9 genes studied, the gene FDXR stood out as the most sensitive and responsive to the low dose of radiation received from the localised radiation treatment by the circulating white blood cells. The activation of FDXR was found to depend on the volume of the body exposed with a peak of expression around 8–9 hours after irradiation was delivered. Finally results obtained ex vivo confirmed the results obtained in vivo.

Abstract

External beam radiation therapy leads to cellular activation of the DNA damage response (DDR). DNA double-strand breaks (DSBs) activate the ATM/CHEK2/p53 pathway, inducing the transcription of stress genes. The dynamic nature of this transcriptional response has not been directly observed in vivo in humans. In this study we monitored the messenger RNA transcript abundances of nine DNA damage-responsive genes (CDKN1A, GADD45, CCNG1, FDXR, DDB2, MDM2, PHPT1, SESN1, and PUMA), eight of them regulated by p53 in circulating blood leukocytes at different time points (2, 6–8, 16–18, and 24 h) in cancer patients (lung, neck, brain, and pelvis) undergoing radiotherapy. We discovered that, although the calculated mean physical dose to the blood was very low (0.038–0.169 Gy), an upregulation of Ferredoxin reductase (FDXR) gene transcription was detectable 2 h after exposure and was dose dependent from the lowest irradiated percentage of the body (3.5% whole brain) to the highest, (up to 19.4%, pelvic zone) reaching a peak at 6–8 h. The radiation response of the other genes was not strong enough after such low doses to provide meaningful information. Following multiple fractions, the expression level increased further and was still significantly up-regulated by the end of the treatment. Moreover, we compared FDXR transcriptional responses to ionizing radiation (IR) in vivo with healthy donors’ blood cells exposed ex vivo and found a good correlation in the kinetics of expression from the 8-hours time-point onward, suggesting that a molecular transcriptional regulation mechanism yet to be identified is involved. To conclude, we provided the first in vivo human report of IR-induced gene transcription temporal response of a panel of p53-dependant genes. FDXR was demonstrated to be the most responsive gene, able to reliably inform on the low doses following partial body irradiation of the patients, and providing an expression pattern corresponding to the % of body exposed. An extended study would provide individual biological dosimetry information and may reveal inter-individual variability to predict radiotherapy-associated adverse health outcomes.

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.

Development of Biomarkers for Radiation Biodosimetry and Medical Countermeasures Research: Current Status, Utility, and Regulatory Pathways

Biomarkers are important indicators of biological processes in health or disease. For this reason, they play a critical role in advanced development of radiation biodosimetry tools and medical countermeasures (MCMs). They can aid in the assessment of radiation exposure level, extent of radiation-induced injury, and/or efficacy of an MCM. This meeting report summarizes the presentations and discussions from the 2020 workshop titled, "Biomarkers in Radiation Biodosimetry and Medical Countermeasures," sponsored by the Radiation and Nuclear Countermeasures Program (RNCP) at the National Institute of Allergy and Infectious Diseases (NIAID). The main goals of this meeting were to: 1. Provide an overview on biomarkers and to focus on the state of science with regards to biomarkers specific to radiation biodosimetry and MCMs; 2. Understand developmental challenges unique to the role of biomarkers in the fields of radiation biodosimetry and MCM development; and 3. Identify existing gaps and needs for translational application.

Making the Case for Absorbed Radiation Response Biodosimetry - Utility of a High-Throughput Biodosimetry System

There is an unmet need to provide medical personnel with a Food and Drug Administration (FDA)-approved biodosimetry method for quantifying individualized absorbed dose response to inform treatment decisions for a very large patient population potentially exposed to ionizing radiation in the event of a nuclear incident. Validation of biodosimetry devices requires comparison of absorbed dose estimates to delivered dose as an indication of accuracy; however, comparison to delivered dose does not account for biological variability or an individual's radiosensitivity. As there is no FDA-cleared gene-expression-based biodosimetry method for determining biological response to radiation, results from accuracy comparisons to delivered dose yield relatively wide tolerance intervals or uncertainty. The Arizona State University Biodesign Institute is developing a high-throughput, automated real-time polymerase chain reaction (RT-PCR)-based biodosimetry system that provides absorbed dose estimates for patients exposed to 0-10 Gy from blood collected 1-7 days postirradiation. While the absorbed dose estimates result from a calibration against the actual exposed dose, the reported dose estimate is a measure of response to absorbed dose based on the exposure models used in developing the system. A central concern with biodosimetry test evaluation is how variability in the dose estimate results could affect medical decision-making, and if the biodosimetry test system performance is quantitatively sufficient to inform effective treatment. A risk:benefit analysis of the expected system performance in the proposed intended use environment was performed to address the potential medical utility of this biodosimetry system. Uncertainty analysis is based on biomarker variability in non-human primate (NHP) models. Monte Carlo simulation was employed to test multiple groups of biomarkers and their potential variation in response to determine uncertainty associated with dose estimate results. Dose estimate uncertainty ranges from ±1.2-1.7 Gy depending on the exposure dose over a range of 2-10 Gy. The risk:benefit of individualized absorbed dose estimates within the context of medical interventions after a nuclear incident is considered and the application of the biodosimetry system is evaluated in this framework. NHP dose-response relationships, as measured by clinical outcome end points, show expected biological and radiosensitivity responses in the primate populations tested and corroborate the biological variability observed in the reported absorbed dose estimate. Performance is examined in relationship to current clinical management and treatment recommendations, with evaluation of potential patient risk in over- and underestimating absorbed dose.

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.

Transcriptomics for radiation biodosimetry: progress and challenges

PURPOSE

Transcriptomic-based approaches are being developed to meet the needs for large-scale radiation dose and injury assessment and provide population triage following a radiological or nuclear event. This review provides background and definition of the need for new biodosimetry approaches, and summarizes the major advances in this field. It discusses some of the major model systems used in gene signature development, and highlights some of the remaining challenges, including individual variation in gene expression, potential confounding factors, and accounting for the complexity of realistic exposure scenarios.

CONCLUSIONS

Transcriptomic approaches show great promise for both dose reconstruction and for prediction of individual radiological injury. However, further work will be needed to ensure that gene expression signatures will be robust and appropriate for their intended use in radiological or nuclear emergencies.

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

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

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.

Role of blood derived cell fractions, temperature and sample transport on gene expression-based biological dosimetry

PURPOSE

For triage purposes following a nuclear accident or a terrorist event, gene expression biomarkers in blood have been demonstrated to be good bioindicators of ionizing radiation (IR) exposure and can be used to assess the dose received by exposed individuals. Many IR-sensitive genes are regulated by the DNA damage response pathway, and modulators of this pathway could potentially affect their expression level and therefore alter accurate dose estimations. In the present study, we addressed the potential influence of temperature, sample transport conditions and the blood cell fraction analyzed on the transcriptional response of the following radiation-responsive genes: FDXR, CCNG1, MDM2, PHPT1, APOBEC3H, DDB2, SESN1, P21, PUMA, and GADD45.

MATERIALS AND METHODS

Whole blood from healthy donors was exposed to a 2 Gy X-ray dose with a dose rate of 0.5 Gy/min (output 13 mA, 250 kV peak, 0.2 mA) and incubated for 24 h at either 37, 22, or 4 °C. For mimicking the effect of transport conditions at different temperatures, samples incubated at 37 °C for 24 h were kept at 37, 22 or 4 °C for another 24 h. Comparisons of biomarker responses to IR between white blood cells (WBCs), peripheral blood mononuclear cells (PBMCs) and whole blood were carried out after a 2 Gy X-ray exposure and incubation at 37 °C for 24 hours.

RESULTS

Hypothermic conditions (22 or 4 °C) following irradiation drastically inhibited transcriptional responses to IR exposure. However, sample shipment at different temperatures did not affect gene expression level except for SESN1. The transcriptional response to IR of specific genes depended on the cell fraction used, apart from FDXR, CCNG1, and SESN1.

CONCLUSION

In conclusion, temperature during the incubation period and cell fraction but not the storing conditions during transport can influence the transcriptional response of specific genes. However, FDXR and CCNG1 showed a consistent response under all the different conditions tested demonstrating their reliability as individual biological dosimetry biomarkers.

Transcriptome of rhesus macaque (Macaca mulatta) exposed to total-body irradiation

The field of biodosimetry has seen a paradigm shift towards an increased use of molecular phenotyping technologies including omics and miRNA, in addition to conventional cytogenetic techniques. Here, we have used a nonhuman primate (NHP) model to study the impact of gamma-irradiation on alterations in blood-based gene expression. With a goal to delineate radiation induced changes in gene expression, we followed eight NHPs for 60 days after exposure to 6.5 Gy gamma-radiation for survival outcomes. Analysis of differential gene expression in response to radiation exposure yielded 26,944 dysregulated genes that were not significantly impacted by sex. Further analysis showed an increased association of several pathways including IL-3 signaling, ephrin receptor signaling, ErbB signaling, nitric oxide signaling in the cardiovascular system, Wnt/β-catenin signaling, and inflammasome pathway, which were associated with positive survival outcomes in NHPs after acute exposure to radiation. This study provides novel insights into major pathways and networks involved in radiation-induced injuries that may identify biomarkers for radiation injury.

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.

High-throughput meta-analysis and validation of differentially expressed genes as potential biomarkers of ionizing radiation-response

BACKGROUND AND PURPOSE

The high-throughput analysis of gene expression in ionizing radiation (IR)-exposed human peripheral white blood cells (WBC) has emerged as a novel method for biodosimetry markers detection. We aimed to detect IR-exposure differential expressed genes (DEGs) as potential predictive biomarkers for biodosimetry and radioinduced-response.

MATERIALS AND METHODS

We performed a meta-analysis of raw data from public microarrays of ex vivo low linear energy transfer-irradiated human peripheral WBC. Functional enrichment and transcription factors (TF) detection from resulting DEGs were assessed. Six selected DEGs among studies were validated by qRT-PCR on mRNA from human peripheral blood samples from nine healthy human donors 24 h after ex vivo X-rays-irradiation.

RESULTS

We identified 275 DEGs after IR-exposure (parameters: |lfc| ≥ 0.7, q value <0.05), enriched in processes such as regulation after IR-exposure, DNA damage checkpoint, signal transduction by p53 and mitotic cell cycle checkpoint. Among these DEGs, DRAM1, NUDT15, PCNA, PLK2 and TIGAR were selected for qRT-PCR validation. Their expression levels significantly increased at 1-4 Gy respect to non-irradiated controls. Particularly, PCNA increased dose dependently. Curiously, TCF4 (Entrez Gene: 6925), detected as overrepresented TF in the radioinduced DEGs set, significantly decreased post-irradiation.

CONCLUSION

These six DEGs show potential to be proposed as candidates for IR-exposure biomarkers, considering their observed molecular radioinduced-response. Among them, TCF4, bioinformatically detected, was validated herein as an IR-responsive gene.

Two-miRNA-based finger-stick assay for estimation of absorbed ionizing radiation dose

Nuclear radiation and radioactive fallouts resulting from a nuclear weapon detonation or reactor accidents could result in injuries affecting multiple sensitive organs, defined as acute radiation syndrome (ARS). Rapid and early estimation of injuries to sensitive organs using markers of radiation response is critical for identifying individuals who could potentially exhibit ARS; however, there are currently no biodosimetry assays approved for human use. We developed a sensitive microRNA (miRNA)-based blood test for radiation dose reconstruction with ±0.5 Gy resolution at critical dose range. Radiation dose-dependent changes in miR-150-5p in blood were internally normalized by a miRNA, miR-23a-3p, that was nonresponsive to radiation. miR-23a-3p was not highly expressed in blood cells but was abundant in circulation and was released primarily from the lung. Our assay showed the capability for dose estimation within hours to 1 week after exposure using a drop of blood from mice. We tested this biodosimetry assay for estimation of absorbed ionizing radiation dose in mice of varying ages and after exposure to both improvised nuclear device (IND)-spectrum neutrons and gamma rays. Leukemia specimens from patients exposed to fractionated radiation showed depletion of miR-150-5p in blood. We bridged the exposure of these patients to fractionated radiation by comparing responses after fractionated versus single acute exposure in mice. Although validation in nonhuman primates is needed, this proof-of-concept study suggests the potential utility of this assay in radiation disaster management and clinical applications.

Overcoming challenges in human saliva gene expression measurements

Saliva, as a non-invasive and easily accessible biofluid, has been shown to contain RNA biomarkers for prediction and diagnosis of several diseases. However, systematic analysis done by our group identified two problematic issues not coherently described before: (1) most of the isolated RNA originates from the oral microbiome and (2) the amount of isolated human RNA is comparatively low. The degree of bacterial contamination showed ratios up to 1:900,000, so that only about one out of 900,000 RNA copies was of human origin, but the RNA quality (average RIN 6.7 + /− 0.8) allowed for qRT-PCR. Using 12 saliva samples from healthy donors, we modified the methodology to (1) select only human RNA during cDNA synthesis by aiming at the poly(A)+-tail and (2) introduced a pre-amplification of human RNA before qRT-PCR. Further, the manufacturer’s criteria for successful pre-amplification (Ct values ≤ 35 for unamplified cDNA) had to be replaced by (3) proofing linear pre-amplification for each gene, thus, increasing the number of evaluable samples up to 70.6%. When considering theses three modifications unbiased gene expression analysis on human salivary RNA can be performed.

CT Irradiation-induced Changes of Gene Expression within Peripheral Blood Cells

Computed tomography (CT) is a crucial element of medical imaging diagnostics. The widespread application of this technology has made CT one of the major contributors to medical radiation burden, despite the fact that doses per individual CT scan steadily decrease due to the advancement of technology. Epidemiological risk assessment of CT exposure is hampered by the fact that moderate adverse effects triggered by low doses of CT exposure are likely masked by statistical fluctuations. In light of these limitations, there is need of further insights into the biological processes induced by CT scans to complement the existing knowledge base of risk assessment. This prompted us to investigate the early transcriptomic response of ex vivo irradiated peripheral blood of three healthy individuals. Samples were irradiated employing a modern dual-source-CT-scanner with a tube voltage of 150 kV, resulting in an estimated effective dose of 9.6 mSv. RNA was isolated 1 h and 6 h after exposure, respectively, and subsequently analyzed by RNA deep sequencing. Differential gene expression analysis revealed shared upregulation of AEN, FDXR, and DDB2 6 h after exposure in all three probands. All three genes have previously been discussed as radiation responsive genes and have already been implicated in DNA damage response and cell cycle control after DNA damage. In summary, we substantiated the usefulness of AEN, FDXR, and DDB2 as RNA markers of low dose irradiation. Moreover, the upregulation of genes associated with DNA damage reminds one of the genotoxic nature of CT diagnostics even with the low doses currently applied.

New Approaches for Quantitative Reconstruction of Radiation Dose in Human Blood Cells

In the event of a nuclear attack or large-scale radiation event, there would be an urgent need for assessing the dose to which hundreds or thousands of individuals were exposed. Biodosimetry approaches are being developed to address this need, including transcriptomics. Studies have identified many genes with potential for biodosimetry, but, to date most have focused on classification of samples by exposure levels, rather than dose reconstruction. We report here a proof-of-principle study applying new methods to select radiation-responsive genes to generate quantitative, rather than categorical, radiation dose reconstructions based on a blood sample. We used a new normalization method to reduce effects of variability of signal intensity in unirradiated samples across studies; developed a quantitative dose-reconstruction method that is generally under-utilized compared to categorical methods; and combined these to determine a gene set as a reconstructor. Our dose-reconstruction biomarker was trained using two data sets and tested on two independent ones. It was able to reconstruct dose up to 4.5 Gy with root mean squared error (RMSE) of ± 0.35 Gy on a test dataset using the same platform, and up to 6.0 Gy with RMSE of ± 1.74 Gy on a test set using a different platform.

Individual Radiosensitivity in Oncological Patients: Linking Adverse Normal Tissue Reactions and Genetic Features

Introduction: Adverse effects of radiotherapy (RT) significantly affect patient's quality of life (QOL). The possibility to identify patient-related factors that are associated with individual radiosensitivity would optimize adjuvant RT treatment, limiting the severity of normal tissue reactions, and improving patient's QOL. In this study, we analyzed the relationships between genetic features and toxicity grading manifested by RT patients looking for possible biomarkers of individual radiosensitivity.

Methods: Early radiation toxicity was evaluated on 143 oncological patients according to the Common Terminology Criteria for Adverse Events (CTCAE). An individual radiosensitivity (IRS) index defining four classes of radiosensitivity (highly radiosensitive, radiosensitive, normal, and radioresistant) was determined by a G₂-chromosomal assay on ex vivo irradiated, patient-derived blood samples. The expression level of 15 radioresponsive genes has been measured by quantitative real-time PCR at 24 h after the first RT fraction, in blood samples of a subset of 57 patients, representing the four IRS classes.

Results: By applying univariate and multivariate statistical analyses, we found that fatigue was significantly associated with IRS index. Interestingly, associations were detected between clinical radiation toxicity and gene expression (ATM, CDKN1A, FDXR, SESN1, XPC, ZMAT3, and BCL2/BAX ratio) and between IRS index and gene expression (BBC3, FDXR, GADD45A, and BCL2/BAX).

Conclusions: In this prospective cohort study we found that associations exist between normal tissue reactions and genetic features in RT-treated patients. Overall, our findings can contribute to the identification of biological markers to predict RT toxicity in normal tissues.

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.

Measurement of γ-H2AX foci, miRNA-101, and gene expression as a means to quantify radiation-absorbed dose in cancer patients who had undergone radiotherapy

Radiological accidents and nuclear terrorism pose an increased threat to members of the public who, following such an event, would need to be assessed for medical care by fast triage. Assay methods such as chromosome aberrations (CA), cytokinesis-block micronucleus (CBMN) and fluorescence in situ hybridization (FISH) techniques have been well established for dose estimation and their potential for handling more samples has also been proved with automation. However, culturing of lymphocytes is an inevitable step, which limits the potential of these markers for triage. In vitro analysis of gamma-H2AX (γ-H2AX), gene and microRNA (miRNA) markers do not require culturing of lymphocytes, and as such have been suggested as attractive tools for triage. Despite studies reporting in vitro dose-response curves, limited evidence is available evaluating the suitability of these assays in real situations. In this study, we have measured the absorbed dose using γ-H2AX, gene (GADD45A, FDXR, and CDKN1A) and miRNA-101 expression in blood samples of cancer patients (n = 20) who had undergone partial-body radiotherapy and compared with the derived equivalent whole-body doses (EWBD). The obtained results from all patients showed a significant (p < 0.05) increase of γ-H2AX foci in post-irradiated as compared to pre-irradiated samples. Moreover, estimated doses using γ-H2AX foci showed a correlation with the derived EWBD (r² = 0.60, p = 0.0003) and was also shown to be dependent on the irradiated body volume. Consistent with γ-H2AX foci frequency, an increase in fold change expression of genes and miRNA-101 was observed. However, the estimated dose significantly varied among the subjects and showed poor correlation (r² = 0.09, 0.04, 0.01 and 0.03 for GADD45A, FDXR, CDKN1A and miRNA-101, respectively) with EWBD. The overall results suggest that the established in vitro γ-H2AX assay is suitable for the detection of radiation exposure and can also provide an estimate of the dose in in vivo irradiated samples. The genes and miRNA-101 markers showed increased expression; nevertheless, there is a need for further improvements to measure doses accurately using these markers.

Development of a biodosimeter for radiation triage using novel blood protein biomarker panels in humans and non-human primates

Purpose: In a significant nuclear event, hundreds of thousands of individuals will require rapid triage for absorbed radiation to ensure effective medical treatment and efficient use of medical resources. We are developing a rapid screening method to assess whether an individual received an absorbed dose of ≥2 Gy based on the analysis of a specific panel of blood proteins in a fingerstick blood sample.Materials and methods: We studied a data set of 1051 human blood samples obtained from radiotherapy patients, normal healthy individuals, and several special population groups. We compared the findings in humans with those from irradiation studies in non-human primates (NHPs).Results: We identified a panel of three protein biomarkers, salivary alpha amylase (AMY1), Flt3 ligand (FLT3L), and monocyte chemotactic protein 1 (MCP1), which are upregulated in human patients receiving fractionated doses of total body irradiation (TBI) therapy as a treatment for cancer. These proteins exhibited a similar radiation response in NHPs after single acute or fractionated doses of ionizing radiation.Conclusion: Our work provides confidence in this biomarker panel for biodosimetry triage using fingerstick blood samples and in the use of NHPs as a model for irradiated humans.

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.

Gene expression changes in human iPSC-derived cardiomyocytes after X-ray irradiation

Purpose: Radiation-induced heart disease caused by cardiac exposure to ionizing radiation comprises a variety of cardiovascular effects. Research in this field has been hampered by limited availability of clinical samples and appropriate test models. In this study, we wanted to elucidate the molecular mechanisms underlying electrophysiological changes, which we have observed in a previous study. Materials and methods: We employed RNA deep-sequencing of human-induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) 48 h after 5 Gy X-ray irradiation. By comparison to public data from hiPSC-CMs and human myocardium, we verified the expression of cardiac-specific genes in hiPSC-CMs. Results were validated by qRT-PCR. Results: Differentially gene expression analysis identified 39 and 481 significantly up- and down-regulated genes after irradiation, respectively. Besides, a large fraction of genes associated with cell cycle processes, we identified genes implicated in cardiac calcium homeostasis (PDE3B), oxidative stress response (FDXR and SPATA18) and the etiology of cardiomyopathy (SGCD, BBC3 and GDF15). Conclusions: Notably, observed gene expression characteristics specific to hiPSC-CMs might be relevant regarding further investigations of the response to external stressors like radiation. The genes and biological processes highlighted in our study present promising starting points for functional follow-up studies for which hiPSC-CMs could pose an appropriate cell model when cell type specific peculiarities are taken into account.

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

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