Generation of a Transcriptional Radiation Exposure Signature in Human Blood Using Long-Read Nanopore Sequencing

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.

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.

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

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

The Influence of Computed Tomography Contrast Agent on Radiation-Induced Gene Expression and Double-Strand Breaks

After nuclear scenarios, combined injuries of acute radiation syndrome (ARS) with, e.g., abdominal trauma, will occur and may require contrast-enhanced computed tomography (CT) scans for diagnostic purposes. Here, we investigated the effect of iodinated contrast agents on radiation-induced gene expression (GE) changes used for biodosimetry (AEN, BAX, CDKN1A, EDA2R, APOBEC3H) and for hematologic ARS severity prediction (FDXR, DDB2, WNT3, POU2AF1), and on the induction of double-strand breaks (DSBs) used for biodosimetry. Whole blood samples from 10 healthy donors (5 males, 5 females, mean age: 28 ± 2 years) were irradiated with X rays (0, 1 and 4 Gy) with and without the addition of iodinated contrast agent (0.016 ml contrast agent/ml blood) to the blood prior to the exposure. The amount of contrast agent was set to be equivalent to the blood concentration of an average patient (80 kg) during a contrast-enhanced CT scan. After irradiation, blood samples were incubated at 37°C for 20 min (DSB) and 8 h (GE, DSB). GE was measured employing quantitative real-time polymerase chain reaction. DSB foci were revealed by γH2AX + 53BP1 immunostaining and quantified automatically in >927 cells/sample. Radiation-induced differential gene expression (DGE) and DSB foci were calculated using the respective unexposed sample without supplementation of contrast agent as the reference. Neither the GE nor the number of DSB foci was significantly (P = 0.07-0.94) altered by the contrast agent application. However, for some GE and DSB comparisons with/without contrast agent, there were weakly significant differences (P = 0.03-0.04) without an inherent logic and thus are likely due to inter-individual variation. In nuclear events, the diagnostics of combined injuries can require the use of an iodinated contrast agent, which, according to our results, does not alter or influence radiation-induced GE changes and the quantity of DSB foci. Therefore, the gene expression and γH2AX focus assay can still be applied for biodosimetry and/or hematologic ARS severity prediction in such scenarios.

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

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

Integrating evolutionarily conserved mechanism of response to radiation for exploring novel Caenorhabditis elegans radiation-responsive genes for estimation of radiation dose associated with spaceflight

During space exploration, space radiation is widely recognized as an inescapable perilous stressor, owing to its capacity to induce genomic DNA damage and escalate the likelihood of detrimental health outcomes. Rapid and reliable estimation of space radiation dose holds paramount significance in accurately assessing the health risks associated with spaceflight. However, the identification of space radiation-responsive genes, with their potential to serve as early indicators for diagnosing radiation dose associated with spaceflight, continues to pose a significant challenge. In this study, based on the evolutionarily conserved mechanism of radiation response, an in silico analysis method of homologous comparison was performed to identify the Caenorhabditis elegans orthologues of human radiation-responsive genes with possible roles in the major processes of response to radiation, and thereby to explore the potential C. elegans radiation-responsive genes for evaluating the levels of space radiation exposure. The results showed that there were 60 known C. elegans radiation-responsive genes and 211 C. elegans orthologues of human radiation-responsive genes implicated in the major processes of response to radiation. Through an investigation of all available transcriptomic datasets obtained from space-flown C. elegans, it was observed that the expression levels of the majority of these putative C. elegans radiation-responsive genes identified in this study were notably changed across various spaceflight conditions. Furthermore, this study indicated that within the identified genes, 19 known C. elegans radiation-responsive genes and 40 newly identified C. elegans orthologues of human radiation-responsive genes exhibited a remarkable positive correlation with the duration of spaceflight. Moreover, a noteworthy presence of substantial multi-collinearity among the majority of these identified genes was observed. This observation lends support to the possibility of treating each identified gene as an independent indicator of radiation dose in space. Ultimately, a subset of 15 potential radiation-responsive genes was identified, presenting the most promising indicators for estimation of radiation dose associated with spaceflight in C. elegans.

Metabolomics and transcriptomics based multi-omics integration reveals radiation-induced altered pathway networking and underlying mechanism

Recent advancement in integrated multi-omics has significantly contributed to many areas of the biomedical field. Radiation research has also grasped uprising omics technologies in biomarker identification to aid in triage management. Herein, we have used a combinatorial multi-omics approach based on transcriptomics together with metabolomics and lipidomics of blood from murine exposed to 1 Gy (LD) and 7.5 Gy (HD) of total-body irradiation (TBI) for a comprehensive understanding of biological processes through integrated pathways and networking. Both omics displayed demarcation of HD group from controls using multivariate analysis. Dysregulated amino acids, various PC, PE and carnitine were observed along with many dysregulated genes (Nos2, Hmgcs2, Oxct2a, etc.). Joint-Pathway Analysis and STITCH interaction showed radiation exposure resulted in changes in amino acid, carbohydrate, lipid, nucleotide, and fatty acid metabolism. Elicited immune response was also observed by Gene Ontology. BioPAN has predicted Elovl5, Elovl6 and Fads2 for fatty acid pathways, only in HD group. Collectively, the combined omics approach facilitated a better understanding of processes uncovering metabolic pathways. Presumably, this is the first in radiation metabolomics that utilized an integrated omics approach following TBI in mice. Our work showed that omics integration could be a valuable tool for better comprehending the mechanism as well as molecular interactions.

Comparative Analysis of miRNA Expression after Whole-Body Irradiation Across Three Strains of Mice

Whole- or partial-body exposure to ionizing radiation damages major organ systems, leading to dysfunction on both acute and chronic timescales. Radiation medical countermeasures can mitigate acute damages and may delay chronic effects when delivered within days after exposure. However, in the event of widespread radiation exposure, there will inevitably be scarce resources with limited countermeasures to distribute among the affected population. Radiation biodosimetry is necessary to separate exposed from unexposed victims and determine those who requires the most urgent care. Blood-based, microRNA signatures have great potential for biodosimetry, but the affected population in such a situation will be genetically heterogeneous and have varying miRNA responses to radiation. Thus, there is a need to understand differences in radiation-induced miRNA expression across different genetic backgrounds to develop a robust signature. We used inbred mouse strains C3H/HeJ and BALB/c mice to determine how accurate miRNA in blood would be in developing markers for radiation vs. no radiation, low dose (1 Gy, 2 Gy) vs. high dose (4 Gy, 8 Gy), and high risk (8 Gy) vs. low risk (1 Gy, 2 Gy, 4 Gy). Mice were exposed to whole-body doses of 0 Gy, 1 Gy, 2 Gy, 4 Gy, or 8 Gy of X rays. MiRNA expression changes were identified using NanoString nCounter panels on blood RNA collected 1, 2, 3 or 7 days postirradiation. Overall, C3H/HeJ mice had more differentially expressed miRNAs across all doses and timepoints than BALB/c mice. The highest amount of differential expression occurred at days 2 and 3 postirradiation for both strains. Comparison of C3H/HeJ and BALB/c expression profiles to those previously identified in C57BL/6 mice revealed 12 miRNAs that were commonly expressed across all three strains, only one of which, miR-340-5p, displayed a consistent regulation pattern in all three miRNA data. Notably multiple Let-7 family members predicted high-dose and high-risk radiation exposure (Let-7a, Let-7f, Let-7e, Let-7g, and Let-7d). KEGG pathway analysis demonstrated involvement of these predicted miRNAs in pathways related to: Fatty acid metabolism, Lysine degradation and FoxO signaling. These findings indicate differences in the miRNA response to radiation across various genetic backgrounds, and highlights key similarities, which we exploited to discover miRNAs that predict radiation exposure. Our study demonstrates the need and the utility of including multiple animal strains in developing and validating biodosimetry diagnostic signatures. From this data, we developed highly accurate miRNA signatures capable of predicting exposed and unexposed subjects within a genetically heterogeneous population as quickly as 24 h of exposure to radiation.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

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.

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.

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.

Dose Variations Using an X-Ray Cabinet to Establish in vitro Dose-Response Curves for Biological Dosimetry Assays

In biological dosimetry, dose-response curves are essential for reliable retrospective dose estimation of individual exposure in case of a radiation accident. Therefore, blood samples are irradiated in vitro and evaluated based on the applied assay. Accurate physical dosimetry of the irradiation performance is a critical part of the experimental procedure and is influenced by the experimental setup, especially when X-ray cabinets are used. The aim of this study was to investigate variations and pitfalls associated with the experimental setups used to establish calibration curves in biological dosimetry with X-ray cabinets. In this study, irradiation was performed with an X-ray source (195 kV, 10 mA, 0.5 mm Cu filter, dose rate 0.52 Gy/min, 1st and 2nd half-value layer = 1.01 and 1.76 mm Cu, respectively, average energy 86.9 keV). Blood collection tubes were irradiated with a dose of 1 Gy in vertical or horizontal orientation in the center of the beam area with or without usage of an additional fan heater. To evaluate the influence of the setups, physical dose measurements using thermoluminescence dosimeters, electron paramagnetic resonance dosimetry and ionization chamber as well as biological effects, quantified by dicentric chromosomes and micronuclei, were compared. This study revealed that the orientation of the sample tubes (vertical vs. horizontal) had a significant effect on the radiation dose with a variation of -41% up to +49% and contributed to a dose gradient of up to 870 mGy inside the vertical tubes due to the size of the sample tubes and the associated differences in the distance to the focal point of the tube. The number of dicentric chromosomes and micronuclei differed by ~30% between both orientations. An additional fan heater had no consistent impact. Therefore, dosimetric monitoring of experimental irradiation setups is mandatory prior to the establishment of calibration curves in biological dosimetry. Careful consideration of the experimental setup in collaboration with physicists is required to ensure traceability and reproducibility of irradiation conditions, to correlate the radiation dose and the number of aberrations correctly and to avoid systematical bias influencing the dose estimation in the frame of biological dosimetry.

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.

The Potential of Omics in Biological Dosimetry

Biological dosimetry is an internationally recognized method for quantifying and estimating radiation dose following suspected or verified excessive exposure to ionising radiation. In severe radiation accidents where a large number of people are potentially affected, it is possible to distinguish irradiated from non-irradiated people in order to initiate appropriate medical care if necessary. In addition to severe incidents caused by technical failure, environmental disasters, military actions, or criminal abuse, there are also radiation accidents in which only one or a few individuals are affected in the frame of occupational or medical exposure. The requirements for biological dosimetry are fundamentally different for these two scenarios. In particular, for large-scale radiation accidents, pre-screening methods are necessary to increase the throughput of samples for a rough first-dose categorization. The rapid development and increasing use of omics methods in research as well as in individual applications provides new opportunities for biological dosimetry. In addition to the discovery and search for new biomarkers, dosimetry assays based on omics technologies are becoming increasingly interesting and hold great potential, especially for large-scale dosimetry. In the following review, the different areas of biological dosimetry, the problems in finding suitable biomarkers, the current status of biomarker research based on omics, the potential applications of assays using omics technologies, and also the limitations for the different areas of biological dosimetry are discussed.

The transcriptomic revolution and radiation biology

PURPOSE

This article will briefly review the origins and evolution of functional genomics, first describing the experimental technology, and then some of the approaches applied to data analysis and visualization. It will emphasize application of functional genomics to radiation biology, using examples from the author's work to illustrate several key types of analysis. It concludes with a look at non-coding RNA, alternative reading of the genome, and single-cell transcriptomics, some of the innovative areas that may help to shape future research in radiation biology and oncology.

CONCLUSIONS

Transcriptomic approaches have provided insight into many areas of radiation biology and medicine, and innovations in technology and data analysis approaches promise continued contributions to radiation science in the future.

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.

The future of biological dosimetry in mass casualty radiation emergency response, personalized radiation risk estimation and space radiation protection

PURPOSE

The aim of this brief personal, high level review is to consider the state of the art for biological dosimetry for radiation routine and emergency response, and the potential future progress in this fascinating and active field. Four areas in which biomarkers may contribute to scientific advancement through improved dose and exposure characterization, as well as potential contributions to personalized risk estimation, are considered: emergency dosimetry, molecular epidemiology, personalized medical dosimetry, and space travel.

CONCLUSION

Ionizing radiation biodosimetry is an exciting field which will continue to benefit from active networking and collaboration with the wider fields of radiation research and radiation emergency response to ensure effective, joined up approaches to triage; radiation epidemiology to assess long term, low dose, radiation risk; radiation protection of workers, optimization and justification of radiation for diagnosis or treatment of patients in clinical uses, and protection of individuals traveling to space.

A 4-Gene Signature of CDKN1, FDXR, SESN1 and PCNA Radiation Biomarkers for Prediction of Patient Radiosensitivity

The quest for the discovery and validation of radiosensitivity biomarkers is ongoing and while conventional bioassays are well established as biomarkers, molecular advances have unveiled new emerging biomarkers. Herein, we present the validation of a new 4-gene signature panel of CDKN1, FDXR, SESN1 and PCNA previously reported to be radiation-responsive genes, using the conventional G2 chromosomal radiosensitivity assay. Radiation-induced G2 chromosomal radiosensitivity at 0.05 Gy and 0.5 Gy IR is presented for a healthy control (n = 45) and a prostate cancer (n = 14) donor cohort. For the prostate cancer cohort, data from two sampling time points (baseline and Androgen Deprivation Therapy (ADT)) is provided, and a significant difference (p > 0.001) between 0.05 Gy and 0.5 Gy was evident for all donor cohorts. Selected donor samples from each cohort also exposed to 0.05 Gy and 0.5 Gy IR were analysed for relative gene expression of the 4-gene signature. In the healthy donor cohort, there was a significant difference in gene expression between IR dose for CDKN1, FXDR and SESN1 but not PCNA and no significant difference found between all prostate cancer donors, unless they were classified as radiation-induced G2 chromosomal radiosensitive. Interestingly, ADT had an effect on radiation response for some donors highlighting intra-individual heterogeneity of prostate cancer donors.

Statistical Modeling of High Dimensional Counts

Statistical modeling of count data from RNA sequencing (RNA-seq) experiments is important for proper interpretation of results. Here I will describe how count data can be modeled using count distributions, or alternatively analyzed using nonparametric methods. I will focus on basic routines for performing data input, scaling/normalization, visualization, and statistical testing to determine sets of features where the counts reflect differences in gene expression across samples. Finally, I discuss limitations and possible extensions to the models presented here.

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.

An ionising radiation-induced specific transcriptional signature of inflammation-associated genes in whole blood from radiotherapy patients: a pilot study

BACKGROUND

This communication reports the identification of a new panel of transcriptional changes in inflammation-associated genes observed in response to ionising radiation received by radiotherapy patients.

METHODS

Peripheral blood samples were taken with ethical approval and informed consent from a total of 20 patients undergoing external beam radiotherapy for breast, lung, gastrointestinal or genitourinary tumours. Nanostring nCounter analysis of transcriptional changes was carried out in samples prior and 24 h post-delivery of the 1st radiotherapy fraction, just prior to the 5th or 6th fraction, and just before the last fraction.

RESULTS

Statistical analysis with BRB-ArrayTools, GLM MANOVA and nSolver, revealed a radiation responsive panel of genes which varied by patient group (type of cancer) and with time since exposure (as an analogue for dose received), which may be useful as a biomarker of radiation response.

CONCLUSION

Further validation in a wider group of patients is ongoing, together with work towards a full understanding of patient specific responses in support of personalised approaches to radiation medicine.

The long and the short of it: unlocking nanopore long-read RNA sequencing data with short-read differential expression analysis tools

Application of Oxford Nanopore Technologies’ long-read sequencing platform to transcriptomic analysis is increasing in popularity. However, such analysis can be challenging due to the high sequence error and small library sizes, which decreases quantification accuracy and reduces power for statistical testing. Here, we report the analysis of two nanopore RNA-seq datasets with the goal of obtaining gene- and isoform-level differential expression information. A dataset of synthetic, spliced, spike-in RNAs (‘sequins’) as well as a mouse neural stem cell dataset from samples with a null mutation of the epigenetic regulator Smchd1 was analysed using a mix of long-read specific tools for preprocessing together with established short-read RNA-seq methods for downstream analysis. We used limma-voom to perform differential gene expression analysis, and the novel FLAMES pipeline to perform isoform identification and quantification, followed by DRIMSeq and limma-diffSplice (with stageR) to perform differential transcript usage analysis. We compared results from the sequins dataset to the ground truth, and results of the mouse dataset to a previous short-read study on equivalent samples. Overall, our work shows that transcriptomic analysis of long-read nanopore data using long-read specific preprocessing methods together with short-read differential expression methods and software that are already in wide use can yield meaningful results.

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.

In Vivo Validation of Alternative FDXR Transcripts in Human Blood in Response to Ionizing Radiation

Following cell stress such as ionising radiation (IR) exposure, multiple cellular pathways are activated. We recently demonstrated that ferredoxin reductase (FDXR) has a remarkable IR-induced transcriptional responsiveness in blood. Here, we provided a first comprehensive FDXR variant profile following DNA damage. First, specific quantitative real-time polymerase chain reaction (qPCR) primers were designed to establish dose-responses for eight curated FDXR variants, all up-regulated after IR in a dose-dependent manner. The potential role of gender on the expression of these variants was tested, and neither the variants response to IR nor the background level of expression was profoundly affected; moreover, in vitro induction of inflammation temporarily counteracted IR response early after exposure. Importantly, transcriptional up-regulation of these variants was further confirmed in vivo in blood of radiotherapy patients. Full-length nanopore sequencing was performed to identify other FDXR variants and revealed the high responsiveness of FDXR-201 and FDXR-208. Moreover, FDXR-218 and FDXR-219 showed no detectable endogenous expression, but a clear detection after IR. Overall, we characterised 14 FDXR transcript variants and identified for the first time their response to DNA damage in vivo. Future studies are required to unravel the function of these splicing variants, but they already represent a new class of radiation exposure biomarkers.

The 2019 Bo Lindell Laureate Lecture: On the use of interdisciplinary, stakeholder-driven, radiation protection research in support of medical uses of ionising radiation

Medical exposures form the largest manmade contributor to total ionising radiation exposure of the UK population. In recent years, new technologies have been developed to improve treatment and prognosis of individuals treated with radiation for diseases such as cancer. However, there is evidence of public, patient, and medical professional concern that radiation protection regulations and practices, as well as understanding of potential long-term adverse health effects of radiation exposure (in the context of other health risks), have not always 'kept pace' with technological developments in this field. This is a truly complex, multi-disciplinary problem for the modern world.The 'Radiation Theme' of the Public Health England and Newcastle University Health Protection Research Unit on 'Chemical and Radiation Threats and Hazards' is addressing this need, with a key focus on a genuinely interdisciplinary approach bringing together world-leading epidemiologists, radiation biologists, clinicians, statisticians, and artists. In addition, the project has a strong grounding in public, patient, and medical professional involvement in research. Similarly, the EU-CONCERT-funded LDLensRad project seeks to understand the mechanisms of action of low-dose ionising radiation in the lens of the eye, and the potential contribution to the development of cataract - in contemporary research, such projects will only be considered successful when they make use of expertise from a variety of fields and when they are able to demonstrate that the outputs are not only of benefit to society, but that society understands and welcomes the benefits. Finally, successful engagement, training, and retention of early career scientists within this field is crucial for sustainability of the research. Herein, the contribution of embedded interdisciplinary working, stakeholder involvement, and training of early career scientists to recent advancements in the field of medical (and wider) radiation protection research is discussed and considered.