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.

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.

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.

Examining potential confounding factors in gene expression analysis of human saliva and identifying potential housekeeping genes.

Isolation of RNA from whole saliva, a non-invasive and easily accessible biofluid that is an attractive alternative to blood for high-throughput biodosimetry of radiological/nuclear victims might be of clinical significance for prediction and diagnosis of disease. In a previous analysis of 12 human samples we identified two challenges to measuring gene expression from total RNA: (1) the fraction of human RNA in whole saliva was low and (2) the bacterial contamination was overwhelming. To overcome these challenges, we performed selective cDNA synthesis for human RNA species only by employing poly(A)+-tail primers followed by qRT-PCR. In the current study, this approach was independently validated on 91 samples from 61 healthy donors. Additionally, we used the ratio of human to bacterial RNA to adjust the input RNA to include equal amounts of human RNA across all samples before cDNA synthesis, which then ensured comparable analysis using the same base human input material. Furthermore, we examined relative levels of ten known housekeeping genes, and assessed inter- and intra-individual differences in 61 salivary RNA isolates, while considering effects of demographical factors (e.g. sex, age), epidemiological factors comprising social habits (e.g. alcohol, cigarette consumption), oral hygiene (e.g. flossing, mouthwash), previous radiological diagnostic procedures (e.g. number of CT-scans) and saliva collection time (circadian periodic). Total human RNA amounts appeared significantly associated with age only (P ≤ 0.02). None of the chosen housekeeping genes showed significant circadian periodicity and either did not associate or were weakly associated with the 24 confounders examined, with one exception, 60% of genes were altered by mouthwash. ATP6, ACTB and B2M represented genes with the highest mean baseline expression (Ct-values ≤ 30) and were detected in all samples. Combining these housekeeping genes for normalization purposes did not decrease inter-individual variance, but increased the robustness. In summary, our work addresses critical confounders and provides important information for the successful examination of gene expression in human whole saliva.

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.

Effect of dose rate on residual γ-H2AX levels and frequency of micronuclei in X-irradiated mouse lymphocytes.

The biological risks associated with low-dose-rate (LDR) radiation exposures are not yet well defined. To assess the risk related to DNA damage, we compared the yields of two established biodosimetry end points, γ-H2AX and micronuclei (MNi), in peripheral mouse blood lymphocytes after prolonged in vivo exposure to LDR X rays (0.31 cGy/min) vs. acute high-dose-rate (HDR) exposure (1.03 Gy/min). C57BL/6 mice were total-body irradiated with 320 kVP X rays with doses of 0, 1.1, 2.2 and 4.45 Gy. Residual levels of total γ-H2AX fluorescence in lymphocytes isolated 24 h after the start of irradiation were assessed using indirect immunofluorescence methods. The terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay was used to determine apoptotic cell frequency in lymphocytes sampled at 24 h. Curve fitting analysis suggested that the dose response for γ-H2AX yields after acute exposures could be described by a linear dependence. In contrast, a linear-quadratic dose-response shape was more appropriate for LDR exposure (perhaps reflecting differences in repair time after different LDR doses). Dose-rate sparing effects (P < 0.05) were observed at doses ≤2.2 Gy, such that the acute dose γ-H2AX and TUNEL-positive cell yields were significantly larger than the equivalent LDR yields. At the 4.45 Gy dose there was no difference in γ-H2AX expression between the two dose rates, whereas there was a two- to threefold increase in apoptosis in the LDR samples compared to the equivalent 4.45 Gy acute dose. Micronuclei yields were measured at 24 h and 7 days using the in vitro cytokinesis-blocked micronucleus (CBMN) assay. The results showed that MNi yields increased up to 2.2 Gy with no further increase at 4.45 Gy and with no detectable dose-rate effect across the dose range 24 h or 7 days post exposure. In conclusion, the γ-H2AX biomarker showed higher sensitivity to measure dose-rate effects after low-dose LDR X rays compared to MNi formation; however, confounding factors such as variable repair times post exposure, increased cell killing and cell cycle block likely contributed to the yields of MNi with accumulating doses of ionizing radiation.

Integrated interdisciplinary training in the radiological sciences.

The radiation sciences are increasingly interdisciplinary, both from the research and the clinical perspectives. Beyond clinical and research issues, there are very real issues of communication between scientists from different disciplines. It follows that there is an increasing need for interdisciplinary training courses in the radiological sciences. Training courses are common in biomedical academic and clinical environments, but are typically targeted to scientists in specific technical fields. In the era of multidisciplinary biomedical science, there is a need for highly integrated multidisciplinary training courses that are designed for, and are useful to, scientists who are from a mix of very different academic fields and backgrounds. We briefly describe our experiences running such an integrated training course for researchers in the field of biomedical radiation microbeams, and draw some conclusions about how such interdisciplinary training courses can best function. These conclusions should be applicable to many other areas of the radiological sciences. In summary, we found that it is highly beneficial to keep the scientists from the different disciplines together. In practice, this means not segregating the training course into sections specifically for biologists and sections specifically for physicists and engineers, but rather keeping the students together to attend the same lectures and hands-on studies throughout the course. This structure added value to the learning experience not only in terms of the cross fertilization of information and ideas between scientists from the different disciplines, but also in terms of reinforcing some basic concepts for scientists in their own discipline.

Radiation induced COX-2 expression and mutagenesis at non-targeted lung tissues of gpt delta transgenic mice.

BACKGROUND: Although radiation-induced bystander effects have been confirmed using a variety of endpoints, the mechanism(s) underlying these effects are not well understood, especially for in vivo study. METHODS: A 1-cm(2) area (1 cm × 1 cm) in the lower abdominal region of gpt delta transgenic mice was irradiated with 5 Gy of 300 keV X-rays, and changes in out-of-field lung and liver were observed. RESULTS: Compared with sham-treated controls, the Spi(-) mutation frequency increased 2.4-fold in non-targeted lung tissues at 24 h after partial body irradiation (PBIR). Consistent with dramatic Cyclooxygenase 2 (COX-2) induction in the non-targeted bronchial epithelial cells, increasing levels of prostaglandin, together with 8-hydroxydeoxyguanosine, in the out-of-field lung tissues were observed after PBIR. In addition, DNA double-strand breaks and apoptosis were induced in bystander lung tissues after PBIR. CONCLUSION: The PBIR induces DNA damage and mutagenesis in non-targeted lung tissues, especially in bronchial epithelial cells, and COX-2 has an essential role in bystander mutagenesis.

A nucleotide excision repair master-switch: p53 regulated coordinate induction of global genomic repair genes.

The tumor suppressor gene p53 is mutated in many human cancers. One of its major roles is as a transcription factor, and its many effector genes control key cellular processes including cell cycle checkpoints and apoptosis. An important role in DNA repair is also emerging for both p53 itself and some of its effector genes. The products of two p53-regulated genes, GADD45a and DDB2, are now known to participate in the global genomic repair (GGR) sub-pathway of nucleotide excision repair (NER). We recently reported the induction of a third GGR gene, XPC, following exposure of normal human peripheral blood lymphocytes to gamma-rays. We now show that XPC is induced in a variety of human cell lines in response to both ionizing and ultra-violet (UV) radiation and alkylating agents, and that this induction requires wild-type p53.

Induction of gene expression as a monitor of exposure to ionizing radiation.

The complex molecular responses to genotoxic stress are mediated by a variety of regulatory pathways. The transcription factor TP53 plays a central role in the cellular response to DNA-damaging agents such as ionizing radiation, but other pathways also play important roles. In addition, differences in radiation quality, such as the exposure to high-LET radiation that occurs during space travel, may influence the pattern of responses. The premise is developed that stress gene responses can be employed as molecular markers for radiation exposure using a combination of informatics and functional genomics approaches. Published studies from our laboratory have already demonstrated such transcriptional responses with doses of gamma rays as low as 2 cGy, and in peripheral blood lymphocytes (PBLs) irradiated ex vivo with doses as low as 20 cGy. We have also found several genes elevated in vivo 24 h after whole-body irradiation of mice with 20 cGy. Such studies should provide insight into the molecular responses to physiologically relevant doses, which cannot necessarily be extrapolated from high-dose studies. In addition, ongoing experiments are identifying large numbers of potential biomarkers using microarray hybridization and various irradiation protocols including expression at different times after exposure to low- and high-LET radiation. Computation-intensive informatics analysis methods are also being developed for management of the complex gene expression profiles resulting from these experiments. With further development of these approaches, it may be feasible to monitor changes in gene expression after low-dose radiation exposure and other physiological stresses that may be encountered during manned space flight, such as the planned mission to Mars.

Physiological function as regulation of large transcriptional programs: the cellular response to genotoxic stress.

The responses to ionizing radiation and other genotoxic environmental stresses are complex and are regulated by a number of overlapping molecular pathways. One such stress signaling pathway involves p53, which regulates the expression of over 100 genes already identified. It is also becoming increasingly apparent that the pattern of stress gene expression has some cell type specificity. It may be possible to exploit these differences in stress gene responsiveness as molecular markers through the use of a combined informatics and functional genomics approach. The techniques of microarray analysis potentially offer the opportunity to monitor changes in gene expression across the entire set of expressed genes in a cell or organism. As an initial step in the development of a functional genomics approach to stress gene analysis, we have recently demonstrated the utility of cDNA microarray hybridization to measure radiation-stress gene responses and identified a number of previously unknown radiation-regulated genes. The responses of some of these genes to DNA-damaging agents vary widely in cell lines from different tissues of origin and different genetic backgrounds. While this again highlights the importance of a cellular context to genotoxic stress responses, it also raises the prospect of expression-profiling of cell lines, tissues, and tumors. Such profiles may have a predictive value if they can define regions of 'expression space' that correlate with important endpoints, such as response to cancer therapy regimens, or identification of exposures to environmental toxins.

Gene expression profiles for monitoring radiation exposure.

Previous demonstrations that the dose, dose rate, radiation quality, and elapsed time since ionising radiation exposure result in variations in the response of stress genes suggest that gene expression signatures may be informative markers of radiation exposure. Defining sets of genes that ate specific for different outcomes of interest will be key to such an approach. A generalised post-exposure prolile may identify exposed individuals within a population, while more specific fingerprints may reveal details of a radiation exposure. Changes in gene expression in human cell lines occur after as little as 0.02 Gy rays, and in peripheral blood lymphocytes alter as little as 0.2 Gy. Diverse genes are also elevated in vivo in mice 24 h after 0.2 Gy irradiation. Ongoing microarray analyses meanwhile continue to identify large numbers of potential biomarkers from varied irradiation protocols. Development of computation-intensive informatics analysis methods will be needed for management of the complex gene expression profiles resulting from such experiments. Although the preliminary data are encouraging, significant work remains before meaningful correlations with risk or practical assessment of exposure can be made by gene expression profiling.

Biological indicators for the identification of ionizing radiation exposure in humans.

While the effects of acute high-dose irradiation are well-documented, less is known about the effects of low level chronic radiation exposure. Physical dosimetry cannot always be relied upon, so dose estimates and determination of past radiation exposure must often be based upon biological indicators. Some of the established methods used in the assessment of nuclear accidents are reviewed here, including cytogenetic analyses, mutation-based assays and electron spin resonance. As interest in research on low-level radiation exposures expands, there is an increasing need for new biomarkers that can identify exposed individuals in human populations. Developments in high-throughput gene expression profiling may enable future development of a rapid and noninvasive testing method for application to potentially exposed populations.

An informatics approach identifying markers of chemosensitivity in human cancer cell lines.

We have used a sensitive and reproducible method of measuring mRNA expression to compare basal levels of 10 transcripts in the 60 cell lines of the National Cancer Institute's in vitro anticancer drug screen (NCI-ACDS) under conditions of exponential growth. The strongest correlation among these target genes was between levels of CIP1/WAF1 and BAX. Levels of the three major growth arrest and DNA damage-inducible gene transcripts, (GADD34, GADD45, and GADD153), which are coordinately regulated in response to many stresses, were also correlated across the 60 cell lines. Although the stress induction of several of the transcripts studied here has been shown to be dependent on wild-type p53 status, basal levels of only CIP1/WAF1 and BAX were found to correlate with p53 status. As expected, basal expression of O6 alkyl guanine alkyl-transferase correlated well with resistance to O6-alkylating agents (r = -0.44) but not with resistance to alkylators with different mechanisms of action (r = -0.04). When basal expression levels of the 10 genes across the NCI-ACDS panel were compared with sensitivities to a panel of 122 standard chemotherapy agents, the most striking relationship was a strong negative correlation (r = -0.3) between basal BCL-X levels and sensitivity to drugs in all of the mechanistic classes except one class of antimetabolites. Sensitivities to a maximally diverse sample of 1200 from 70,000 compounds tested in the NCI-ACDS of agents were also negatively correlated with BCL-X levels. A novel application of factor analysis revealed that the newly discovered associations were independent of previously demonstrated sensitivity factors such as p53 mutation status and native population doubling time. A similar pattern of correlation was seen for Bcl-X(L) protein levels. Conversely, BAX and BCL2, two other genes associated with regulation of apoptosis, showed no overall correlation with drug sensitivities. This suggests that BCL-X may play a unique role in general resistance to cytotoxic agents, with the cell lines demonstrating relative resistance to 70,000 cytotoxic agents in the NCI-ACDS being characterized by high BCL-X expression.

Identification of potential mRNA biomarkers in peripheral blood lymphocytes for human exposure to ionizing radiation.

Since early in the Atomic Age, biological indicators of radiation exposure have been sought, but currently available methods are not entirely satisfactory. Using cDNA microarray hybridization to discover new potential biomarkers, we have identified genes expressed at increased levels in human peripheral blood lymphocytes after ex vivo irradiation. We recently used this technique to identify a large set of ionizing radiation-responsive genes in a human cell line (Oncogene 18, 3666-3672, 1999). The present set of radiation markers in peripheral blood lymphocytes was identified 24 h after treatment, and while the magnitude of mRNA induction generally decreased over time, many markers were still significantly elevated up to 72 h after irradiation. In all donors, the most highly responsive gene identified was DDB2, which codes for the p48 subunit of XPE, a protein known to play a crucial role in repair of ultraviolet (UV) radiation damage in DNA. Induction of DDB2, CDKN1A (also known at C1P1/WAF1) and XPC showed a linear dose-response relationship between 0.2 and 2 Gy at 24 and 48 h after irradiation, with less linearity at earlier or later times. These results suggest that relative levels of gene expressions in peripheral blood cells may provide estimated of environmental radiation exposures.

p53-mediated DNA repair responses to UV radiation: studies of mouse cells lacking p53, p21, and/or gadd45 genes.

Human cells lacking functional p53 exhibit a partial deficiency in nucleotide excision repair (NER), the pathway for repair of UV-induced DNA damage. The global genomic repair (GGR) subpathway of NER, but not transcription-coupled repair (TCR), is mainly affected by p53 loss or inactivation. We have utilized mouse embryo fibroblasts (MEFs) lacking p53 genes or downstream effector genes of the p53 pathway, gadd45 (Gadd45a) or p21 (Cdkn1a), as well as MEFs lacking both gadd45 and p21 genes to address the potential contribution of these downstream effectors to p53-associated DNA repair. Loss of p53 or gadd45 had a pronounced effect on GGR, while p21 loss had only a marginal effect, determined by measurements of repair synthesis (unscheduled DNA synthesis), by immunoassays to detect removal of UV photoproducts from genomic DNA, and by assays determining strand-specific removal of CPDs from the mouse dhfr gene. Taken together, the evidence suggests a role for Gadd45, but relatively little role for p21, in DNA repair responses to UV radiation. Recent evidence suggests that Gadd45 binds to UV-damaged chromatin and may affect lesion accessibility. MEFs lacking p53 or gadd45 genes exhibited decreased colony-forming ability after UV radiation and cisplatin compared to wild-type MEFs, indicating their sensitivity to DNA damage. We provide evidence that Gadd45 affects chromatin remodelling of templates concurrent with DNA repair, thus indicating that Gadd45 may participate in the coupling between chromatin assembly and DNA repair.

The complexity of radiation stress responses: analysis by informatics and functional genomics approaches.

Molecular responses to genotoxic stress are complex and are mediated by a variety of regulatory pathways. One key element in cellular response is the stress gene transcription factor p53, which can regulate nearly 100 genes that have already been identified. Although p53 plays a central role in the cellular response to DNA-damaging agents such as ionizing radiation (IR), other pathways can also have important roles. One example is the transcriptional responses associated with IR-induced apoptosis, where induction of some genes is limited to p53 wild-type (wt) cells that also have the ability to undergo rapid apoptosis after irradiation. In contrast, other genes are triggered after IR in lines undergoing rapid apoptosis regardless of p53 status. From this and other examples, it is apparent that the pattern of stress gene expression is cell type specific in both primary and transformed lines. The premise will be developed that such differences in stress gene responsiveness can be employed as molecular markers using a combination of informatics and functional genomics approaches. An example is given using the panel of lines of the NCI anticancer drug screen where both the p53 status and sensitivity to a large collection of cytotoxic agents have been determined. The utility of cDNA microarray hybridization to measure IR-stress gene responses has recently been demonstrated and a large number of additional IR-stress genes have been identified. The responses of some of these genes to IR and other DNA-damaging agents varied widely in cell lines from different tissues of origin and different genetic backgrounds, highlighting the importance of cellular context to genotoxic stress responses; this also highlights the need for informatics approaches to discover and prioritize hypotheses regarding the importance of particular cellular factors. The aim of this review is to demonstrate the utility of combining an informatics approach with functional genomics in the study of stress responses.

Induction of stress genes by low doses of gamma rays.

Using cells of a human myeloid tumor cell line (ML-1), we have detected induction of several stress-responsive genes by doses of gamma rays below 50 cGy. We found a linear dose-response relationship for induction of CDKN1A (formerly known as CIP1/WAF1) and GADD45 mRNA levels over the range of 2-50 cGy, with no evidence of a threshold for induction. Although exposures to 2 and 5 cGy did not result in any detectable reduction in cloning efficiency or increased apoptosis in ML-1 cells, these exposures did produce a transient delay of cells in the phases of the cell cycle in addition to the observed up-regulation of CDKN1A and GADD45. The relative induction of genes such as CDKN1A by radiation doses that produce little toxicity indicates that surviving cells do contribute significantly to the observed stress responses. These studies should provide insight into the molecular responses to physiologically relevant doses that cannot necessarily be extrapolated from high-dose studies.

Fluorescent cDNA microarray hybridization reveals complexity and heterogeneity of cellular genotoxic stress responses.

The fate of cells exposed to ionizing radiation (IR) may depend greatly on changes in gene expression, so that an improved view of gene induction profiles is important for understanding mechanisms of checkpoint control, repair and cell death following such exposures. We have used a quantitative fluorescent cDNA microarray hybridization approach to identify genes regulated in response to 7-irradiation in the p53 wild-type ML-1 human myeloid cell line. Hybridization of the array to fluorescently-labeled RNA from treated and untreated cells was followed by computer analysis to derive relative changes in expression levels of the genes present in the array, which agreed well with actual quantitative changes in expression. Forty-eight sequences, 30 not previously identified as IR-responsive, were significantly regulated by IR. Induction by IR and other stresses of a subset of these genes, including the previously characterized CIP1/ WAF1, MDM2 and BAX genes, as well as nine genes not previously reported to be IR-responsive, was examined in a panel of 12 human cell lines. Responses varied widely in cell lines with different tissues of origin and different genetic backgrounds, highlighting the importance of cellular context to genotoxic stress responses. Two of the newly identified IR-responsive genes, FRA-1 and ATF3, showed a p53-associated component to their IR-induction, and this was confirmed both in isogenic human cell lines and in mouse thymus. The majority of the IR-responsive genes, however, showed no indication of p53-dependent regulation, representing a potentially important class of stress-responsive genes in leukemic cells.