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, two–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.
