RENEB Inter-Laboratory Comparison 2021: The Gene Expression Assay

Calibration curve for radiation dose estimation using FDXR gene expression biodosimetry - premises and pitfalls

PURPOSE

Radiation-induced alterations in gene expression show great promise for dose reconstruction and for severity prediction of acute health effects. Among several genes explored as potential biomarkers, FDXR is widely used due to high upregulation in white blood cells following radiation exposure. Nonetheless, the absence of a standardized protocols for gene expression-based biodosimetry is a notable gap that warrants attention to enhance the accuracy, reproducibility and reliability. The objective of this study was to evaluate the sensitivity of transcriptional biodosimetry to differences in protocols used by different laboratories and establish guidelines for the calculation of calibration curve using FDXR expression data.

MATERIAL AND METHODS

Two sets of irradiated blood samples generated during RENEB exercise were used. The first included samples irradiated with known doses including: 0, 0.25, 0.5, 1, 2, 3 and 4 Gy. The second set consisted of three 'blind' samples irradiated with 1.8 Gy, 0.4 Gy and a sham-irradiated sample. After irradiation, samples were incubated at 37 °C over 24 h and sent to participating laboratories, where RNA isolation and FDXR expression analysis by qPCR were performed using sets of primers/probes and reference genes specific for each laboratory. Calibration curves based on FDXR expression data were generated using non-linear and linear regression and used for dose estimation of 'blind' samples.

RESULTS

Dose estimates for sham-irradiated sample (0.020-0.024 Gy) and sample irradiated with 0.4 Gy (0.369-0.381 Gy) showed remarkable consistency across all laboratories, closely approximating the true doses regardless variation in primers/probes and reference genes used. For sample irradiated with 1.8 Gy the dose estimates were less precise (1.198-2.011 Gy) but remained within an acceptable margin for triage within the context of high dose range.

CONCLUSION

Methodological differences in reference genes and primers/probes used for FDXR expression measurement do not have a significant impact on the dose estimates generated, provided that all reference genes performed as expected and the primers/probes target a similar set of transcript variants. The preferred method for constructing a calibration curve based on FDXR expression data involves employing linear regression to establish a function that describes the relationship between the logarithm of absorbed dose and FDXR ΔCt values. However, one should be careful with using non-irradiated sample data as these cannot be accurately represented on a logarithmic scale. A standard curve generated using this approach can give reliable dose estimations in a dose range from 50 mGy to 4 Gy at least.

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.

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.

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.

Paper-Based Vertical Flow Immunoassay for the Point-of-Care Multiplex Detection of Radiation Dosimetry Genes

In a nuclear or radiological incident, first responders must quickly and accurately measure radiation exposure among civilians as medical countermeasures are radiation dose-dependent and time-sensitive. Although several approaches have been explored to measure absorbed radiation dose, there is an important need to develop point-of-care (POC) bioassay devices that can be used immediately to triage thousands of individuals potentially exposed to radiation. Here we present a proof-of-concept study showing the use of a paper-based vertical flow immunoassay (VFI) to detect radiation dosimetry genes. Using labeled primers during amplification and a multiplex membrane, our results showed that the nucleic acid VFI can simultaneously detect two biodosimetry genes, CDKN1A and DDB2, as well as one housekeeping gene MRPS5. The assay demonstrated good linearity and precision with an inter- and intra-assay coefficient of variance <20% and <10%, respectively. Moreover, the assay showed its ability to discriminate non-irradiated controls (0 Gy) from irradiated samples (1 + 2 Gy) with an overall sensitivity of 62.5% and specificity of 100% (AUC = 0.8672, 95% CI: 0.723-1.000; p = 0.004). Interestingly, the gene combination also showed a dose-dependent response for 0, 1, and 2 Gy, similar to data obtained by real-time PCR benchmark. These preliminary results suggest that a VFI platform can be used to detect simultaneously multiple genes that can be then quantified, thus offering a new approach for a POC biodosimetry assay that could be rapidly deployed on-site to test a large population and help triage and medical management after radiological event.

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

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.