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.

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.