Genome-wide gene expression profiling suggests distinct radiation susceptibilities in sporadic and post-Chernobyl papillary thyroid cancers

Papillary thyroid cancers (PTCs) incidence dramatically increased in the vicinity of Chernobyl. The cancer-initiating role of radiation elsewhere is debated. Therefore, we searched for a signature distinguishing radio-induced from sporadic cancers. Using microarrays, we compared the expression profiles of PTCs from the Chernobyl Tissue Bank (CTB, n=12) and from French patients with no history of exposure to ionising radiations (n=14). We also compared the transcriptional responses of human lymphocytes to the presumed aetiological agents initiating these tumours, γ-radiation and H2O2. On a global scale, the transcriptomes of CTB and French tumours are indistinguishable, and the transcriptional responses to γ-radiation and H2O2 are similar. On a finer scale, a 118 genes signature discriminated the γ-radiation and H2O2 responses. This signature could be used to classify the tumours as CTB or French with an error of 15–27%. Similar results were obtained with an independent signature of 13 genes involved in homologous recombination. Although sporadic and radio-induced PTCs represent the same disease, they are distinguishable with molecular signatures reflecting specific responses to γ-radiation and H2O2. These signatures in PTCs could reflect the susceptibility profiles of the patients, suggesting the feasibility of a radiation susceptibility test.

An increased incidence of thyroid carcinomas in children was first noticed in Belarus and Ukraine 4 years after the 1986 Chernobyl accident (Baverstock et al, 1992;Kazakov et al, 1992). Increased incidence has been observed since then in people exposed to fallout during childhood in these regions (United Nations Scientific Committee of the Effect of Atomic Radiation, 2000; Mahoney et al, 2004). The aggressiveness and morphology of these tumours (over 95% classified on the basis of their pathology as papillary thyroid carcinomas (PTCs)) appear to be related to the age of the patients at the time of the accident and to the lag between the accident and diagnosis, that is, the latency of the cancers .
Radiation is the only proven cause of PTC so far. Although the cause of PTC in patients not exposed to radiation remains uncertain, a growing body of evidence suggests that H 2 O 2 could play a role in the absence of radiation. Indeed, it is a potent DNAdamaging agent produced in large amounts during thyroid hormone synthesis (Corvilain et al, 2000). It causes DNA damage (guanine oxidation, single-and double-strand breaks) in human lymphocytes (Turner et al, 2003), hamster ovarian cells (Dahm-Daphi et al, 2000;Mondello et al, 2002), and in human, dog and sheep thyroid cells in primary culture (Chico Galdo et al, 2006). Hydrogen peroxide is believed to destroy follicular thyroid cells in myxoedematous endemic cretinism (Kohrle et al, 2005) and to cause cancers in the thyroid of Tg-a 1B AR mice (Ledent et al, 1997). Lack of protective systems, peroxiredoxin or glutathione peroxidases, in knockout mice lead to cancer (Neumann et al, 2003;Lee et al, 2006). Transfection of an H 2 O 2 -generating system transform epithelial cells (Chu et al, 1996). The spontaneous somatic mutation rate in normal mice and rat thyroid cells is substantially higher than in liver and lung cells (Corvilain et al, 1994). With a turnover of 8.5 years in adults (Coclet et al, 1989), thyrocytes have time to accumulate H 2 O 2 -induced DNA damages. Hydrogen peroxide has been found to play a role in several human cancers (Quinn et al, 2006). Thus, a number of arguments support a role of H 2 O 2 in the initiation of PTC, and in particular in patients not exposed to radiation.
The vast majority of PTCs harbour either a BRAF mutation (45%; Xing, 2005) or a RET/PTC rearrangement (35% in adults; Nikiforov, 2002), which are generally mutually exclusive (Soares et al, 2003). Both gene alterations result in the constitutive activation of the RAS -RAF -MAPK signalling pathway (Kimura et al, 2003;Soares et al, 2003). Gene-expression signatures separating BRAF from RET/PTC tumours have been reported, but the number of genes involved varies from a few dozens (Frattini et al, 2004) to several thousands (Giordano et al, 2005). Although early reports pointed at a lower BRAF mutation frequency in Chernobyl patients, recent evidence suggests that the BRAF mutation is associated with age and is more prevalent among older Chernobyl patients and/or among patients with longer latency tumours (Kumagai et al, 2004;Lima et al, 2004;Powell et al, 2005;Rosenbaum et al, 2005). Several research teams have reported higher frequencies of RET/PTC rearrangements in post-Chernobyl patients (Nikiforov et al, 1997). These higher frequencies could result from the fact that radiation induces double-strand breaks, and thus rearrangements rather than point mutations (Dahm-Daphi et al, 2000), or possibly to a differing molecular profile in childhood vs adult papillary carcinomas (Powell et al, 2005). The induction of RET/PTC rearrangements after in vitro irradiation of immortalised thyroid cells (Caudill et al, 2005) supports the former explanation. Whether the two best-characterised genetic alterations found in PTC are involved in a radiation signature remains an open question. In addition, radiation induces other unknown alterations.
In this paper, we have compared the gene-expression profiles of PTCs from adult French patients with no history of exposure to radiation and from adult Ukrainian patients exposed to Chernobyl fallout during childhood, and asked whether there is a geneexpression signature distinguishing radiation-induced from sporadic cancers. Our preliminary investigation suggested the absence of a large-scale radiation signature (Detours et al, 2005). We extend it here by using a more recent microarray technology, by covering more genes, by studying more patients and by establishing results with a wider range of statistical methods. We confirm that French and Chernobyl Tissue Bank (CTB) tumours have the same overall expression profiles and have indistinguishable BRAF and RET/PTC frequencies.
We also compared the transcriptional responses of human cells to the two most likely aetiological agents of CTB and French tumours; g-radiation and H 2 O 2 . The similarity of CTB and French tumours is mirrored by the similarity of the transcriptional responses to g-radiation and H 2 O 2 . However, subtle expression differences are exploitable to accurately classify these tumours according to their origin. Part of these expression differences includes genes involved in the differential response to H 2 O 2 and radiation, and genes involved in homologous recombination which suggests that different-and detectable-susceptibility profiles lead to sporadic and radiation-induced PTC.

Transcriptional and genetic data
Paired samples of tumoral and adjacent non-tumoral thyroid tissues were obtained from the CTB (www.chernobyltissuebank.com) and from patients undergoing surgery for thyroid disease at the Ambroise Pare Hospital (Boulogne, France). French tissues were immediately frozen in liquid nitrogen and stored at À801C until use. Diagnoses were made by the Department of Pathology at the Ambroise Pare Hospital or by the International Pathology Panel of the CTB. The protocol received approval from the Ethics Committees of the institutions. The detail of BRAF-RET/PTC status determination, RNA processing and microarray data preprocessing is available in Supplementary information file S3. Microarray data are available from the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo), accession number GSE3950. Jarzab et al (2005) data were downloaded from www.genomika.pl/ thyroidcancer/PTCCancerRes.html. We used the original MAS 5.0 normalised expression levels, took the log 2 of expression ratios and averaged over patients. The probes of the two platforms could be matched on the basis of their Entrez IDs for 4203 genes.

Unsupervised classification
Hierarchical clustering was computed with the R language function hclust with Ward linkage. Multidimensional scaling was computed with the R function isoMDS. Both methods were fed Pearson correlation distances as input.

Supervised classification
Support vector machine classification was run with linear kernel and cost ¼ 1 using the rfe 0.2 and e1071 1.5.9 packages for R. The generalised partial least-square (GPLS) implementation from package gpls 1.1.0 (Ding and Gentleman, 2004) for R was run with default parameters. Prediction analysis of microarray (Tibshirani et al, 2002) was run with threshold values in {1.0, 1.1, 1.2, y, 3.0} using pamr 1.25 for R. The random forest classification used default parameters from R package random-Forest 4.5.12 (Zhang et al, 2003). Generalised partial least-square and random forest (RF) were combined with an external genes selection procedure focusing on the n genes with the highest absolute t-statistics, with n in {1, 2 1 , 2 2 , y, 2 13 }. We adopted the inner/outer cross-validation scheme described in details in Ruschhaupt et al (2004) and implemented in the package MCRestimate 1.3.0 to prevent parameter and gene selection biases (Ambroise and McLachlan, 2002). Note that a simpler split-sample validation, in which samples are not recycled as in the current cross-validation protocol, would be suboptimal here because of the limited availability of CTB samples (Simon et al, 2003). A 13-fold cross-validation protocol with each round including parameter and gene selection, and classification was run. At each one of the 13 rounds, the best parameters (including signature size) were estimated by running a nested (inner) 12-fold cross-validation for each combination of parameters. Table 2 presents averages over 10 repetitions of the entire inner/outer cross-validation, each based on a different random 13-fold partitions of the data. The random error was computed by averaging the error of five runs of the complete classification procedure on data with CTB and French labels randomly assigned to samples. The same protocol was used for the classification on the basis of the 118 genes g-radiation vs H 2 O 2 signature, except that the number of genes, n, was chosen in {1, 5, 10, y, 118} and that the tested prediction analysis of microarray (PAM) thresholds were in {0.1, 0.2, y, 3.0}. Classifications on the basis of DNA repair signatures were run without gene selection, and therefore without inner crossvalidation. The PAM threshold was set to 0.5. All P-values were derived by running 1000 times the complete cross-validation with CTB and French labels assigned randomly to samples and counting how many runs produced classification error below the error obtained on the actual data.

Derivation of the c-radiation vs H 2 O 2 signature
We downloaded the Supplementary data set S2 of Amundson et al (2005) from the Oncogene web site (www.nature.com/onc/ index.html). Genes with expression values differing by 1.5-fold between the 2.5 Gy g-radiationand H 2 O 2 (200 mM)-treated TK6 cells were selected. To remove immune system-related genes, we downloaded the gcrma-processed version of the GNF human gene atlas (Su et al, 2004;symatlas.gnf.org), which contains expression profiles of normal tissues in most organs. We performed an unpaired two class Significance Analysis of Microarrays (SAM; Tusher et al, 2001) with class no. 1 including immune systemrelated tissues and white blood cells and class no. 2 including all other tissues. We selected the 20% top-ranking genes, which were all significant at qo0.05, and removed them from the g-radiation and H 2 O 2 signature.

Expression profiles and gene alteration status of PTCs from France and from the Chernobyl Tissue Bank
Expression profiles were determined for the tumours of 14 patients from France with no documented history of exposure to radiation, and 12 tumours from the Chernobyl Tissue Bank (see online Materials and Methods). CTB tumours are papillary cancers collected in young people who were exposed to the Chernobyl accident ( (Thomas et al, 2000), see patient information, Table 1). There are 9 tumours of classical subtype, 4 of follicular subtype and 1 of trabecular subtype among the 14 French PTC samples. There are 8 classical, 3 follicular and 1 solid subtypes among the 12 CTB PTC samples. Three French and four CTB mRNA samples (PTC6, PTC7, PTC11 and S405, S420, S422, S423) were reused from our earlier study (Detours et al, 2005).
The mRNA expression profiles of all tumours were determined with 12 000 EST (8000 genes) cDNA microarrays using patientmatched nontumoural adjacent tissues as controls. To assess the quality of the data, we compared our expression ratios averaged over samples with those of Jarzab et al (2005), who used the Affymetrix s platform. Pearson's correlation measured on the B4000 genes available and expressed in both platforms was 0.72 ( Figure 1).
Tumours were screened for the presence of a RET/PTC rearrangement and for BRAF V600E mutation (Table 1). A RET/ PTC rearrangement was found in 42% (5/12) of the CTB tumours and in 21% (3/14) of the French tumours. The difference between the two groups is not significant according to Fisher's exact test. The BRAF mutation is found in comparable proportions in French (36%, 5/14) and CTB tumours (41%, 5/12). None of these alterations was detected in 30% (8/26) of the tumours.

Chernobyl Tissue Bank and French PTCs have similar overall expression profiles
We first searched for global expression differences between CTB and French PTCs, that is, extensive differences detectable when all the genes present on our arrays are considered. Hierarchical    clustering based on all genes did not reveal a clear separation between these two classes of PTCs (Figure 2, upper panel).
Multidimensional scaling collapses the high-dimensional genes space into two dimensions while preserving the distance relationships between all pairs of samples (Figure 2, lower panel). Figure 2 confirms that French and CTB tumours have similar expression profiles when compared on a global scale although CTB tumours form a more compact group.
Four supervised classification algorithms find multigenes signatures separating CTB from French PTCs The absence of separation between CTB and French PTCs on the basis of all genes or at the level of individual genes, does not exclude that these tumours are distinguishable on the basis of a subset of genes. We investigated this possibility with a supervised classification approach (details in Materials and Methods). To strengthen the reliability of our conclusions, all the results were reproduced with four linear classification procedures: linear kernel support vector machines (LKSVM), GPLS, PAMs and RF. Each one included or was combined with a gene selection procedure, that is, a procedure to uncover multigenes signatures including as few genes as possible. All four approaches were tested using a rigorous inner/outer cross-validation procedure (Materials and Methods). It guaranteed that classification testing was performed on independent samples not used for classifier training. The cross-validation results are presented in Table 2A.
The best performer was GPLS. It misclassified 17% of CTB tumours as French PTC, 7% of French PTCs as CTB, resulting in an overall error rate of 12%. Running the same classification on data in which the 'CTB' and 'French' labels were randomly assigned to the 26 tumour samples led to high error rates of 45% (n ¼ 5, s.d. ¼ 12%, see Materials and Methods), as expected for random classification of slightly unbalanced classes (12 CTB and 14 French samples). Thus, the low error rates were unlikely to result from artefacts, including data overfitting. Figure 3 shows the 256 most classifying genes found by GPLS/t-statistics trained on all 26 samples (corresponding genes listed in Supplementary Table S1). The optimal signature size varied among the different crossvalidation runs from one gene to several thousands of genes, with a median of 256 genes. Such limited stability is widespread, including in large studies (Ein-Dor et al, 2005;Michiels et al, 2005). The three other classification procedures, LKSVM/RFE, RF/ t-test and PAM produced a global error of 15, 23 and 27%, respectively. Thus, undirected selections of classifying genes lead to separation CTB and French tumours.
Note that if the classification results were confined to a subtype, the accuracy would not be as low as 15%, it would be greater than 35%-the classical subtype is the largest, representing 65% of our tumours.

Hydrogen peroxide and c-radiation elicit similar transcriptional responses in lymphocytes
Because hydrogen peroxide, H 2 O 2 , is produced at high levels during thyroid hormone synthesis (Corvilain et al, 2000) and is a well-known DNA-damaging agent, we investigated the possibility that in the absence of an obvious external risk factor, for example radiation, French cancers must have occurred as a result of H 2 O 2 exposure. Amundson et al (2005) measured with microarrays the transcriptional responses of a B-lymphocyte cell line, TK6, to 13 stress agents. These included 10 DNA-damaging agents: H 2 O 2 , radiation (neutron and g-rays at 2.5 and 8 Gy), adriamycin, arsenite, campothecin, CdCl 2 , cisplatin, methyl methanesulphonate and UVB (280À320 nm). We downloaded the expression data published with the paper and produced the hierarchical clustering shown in Figure 4 (see online Materials and Methods). The responses to 200 mM H 2 O 2 and to 2.5 Gy g-radiation clustered together, that is, among 12 stress agents, including 10 DNAdamaging agents, g-radiation at 2.5 Gy elicited the transcriptional response that was the closest to that of H 2 O 2 . We concluded that these similar transcriptional responses reflect similar damages in the cells.

Chernobyl Tissue Bank and French tumours are accurately classified on the basis of genes regulated differently in c-radiation and H 2 O 2 responses
The transcriptional responses to g-radiation and H 2 O 2 are broadly similar; however, some genes are expressed differently between the two in vitro assays. We reasoned that these expression differences may mirror subtle underlying g-radiation and H 2 O 2 susceptibility differences between CTB and French tumours that could be used for classification. We found 293 genes in the 1451 published by Amundson et al (2005) with a fold change greater than 1.5 between the g-radiation (2.5 Gy) and the H 2 O 2 responses (200 mM). These responses were measured in B lymphocytes, whereas our goal was to classify thyroid tumours. Thus, we removed immune system-specific genes from the set of 293 genes (see online Materials and Methods). This filtering left 162 genes. Among them, 118 were spotted on our microarrays. They are listed in Supplementary Table S2 and will be referred to thereafter as the g-radiation vs H 2 O 2 signature. Note that it was derived independently of our PTC data.
Next, we applied the same four classification algorithms as above except that only the independently selected 118 genes were used. Error rates (Table 2B) were comparable to those obtained in Table 2A, where the classifying genes were selected from a list of 8000. Again, all four algorithms classified the tumours with an error rate p27%, GPLS/t-test and LKSVM/RFE being the most accurate with an error rate of 15%. This result shows a relation between the g-radiation vs H 2 O 2 signature and CTB and sporadic carcinomas distinction, which could reflect the underlying aetiology of CTB and French tumours.

Chernobyl Tissue Bank and French tumours are accurately classified on the basis of 13 genes involved in homologous recombination
To focus better on which elements of the DNA-damage response may differ between CTB and French tumours, we investigated if genes involved in the different DNA repair mechanisms led to accurate classification. We collected from the Human DNA Repair Genes database (Wood et al, 2001(Wood et al, , 2005, all the genes known to be involved in base-excision repair, mismatch-excision repair, nucleotide-excision repair, homologous recombination and nonhomologous end joining. The signature from each one of these five repair mechanisms was then used to classify the CTB and French tumours. These signatures contain few genes and were compiled from a source curated by DNA repair experts. Therefore, we skipped the gene selection step, which in turn alleviates the need for time-consuming internal cross-validation. The resulting computational gain made it tractable to run an additional statistical control: all five classification tasks were rerun 1000 times with the CTB and French labels randomly assigned to the tumours to estimate P-values, that is, the odds that the classification error was as low as the one observed with the actual data. Besides this, the classification proceeded exactly as above. The classification error rates were high for base-excision repair, mismatch-excision repair and nonhomologous end joining, regardless of the algorithm (not shown). The nucleotide-excision repair signature produced an error rate of 27% with RF, but B50% with GPLS, PAM and LKSVM. In contrast, the homologous recombination signature (Table 3) led to a classification below 31% for all four procedures, below 20% for two and equal to 15% for LKSVM (Table 2C).
The P-value for RF, 0.064, was slightly above the 0.05 significance standard. All the other P-values were highly significant and remained below 0.02 after adjusting for the fact that five classification tasks were being examined (using Bonferonni correction, i.e., multiplying the P-values by 5). This suggests that homologous recombination, which repairs double-strand breaks, operates differently in CTB and French tumours or in the associated normal tissues. None of the homologous recombination signature genes are part of the 118 genes of the g-radiation vs H 2 O 2 signature. Thus, the homologous recombination and g-radiation vs H 2 O 2 signatures are nonoverlapping. They are thus two different signatures supporting a link between radiation and the CTB/ French PTC expression differences.

DISCUSSION
We compared French and CTB tumours at the level of their global expression profiles, that is, of their overall phenotype. Hierarchical clustering and multidimensional scaling failed to uncover a largescale difference between them. Note that, would such difference exist, our preliminary study (Detours et al, 2005) would have revealed it. Thus, the conclusion of pathologists that sporadic and radiation-induced PTCs are the same type of lesions is supported by expression data.
The similarity of expression profiles on a global scale, as observed with hierarchical clustering performed on all genes, does  not preclude that small groups of genes differ between these profiles. Supervised classification is the tool of choice to evaluate whether a group of genes can be exploited to discriminate different classes of tumours (Allison et al, 2006). Four linear classification algorithms assigned the tumours to the French or CTB groups with an error ranging from 12 to 27%, and p15% for two algorithms. These figures are typical of properly designed microarray studies (Ntzani and Ioannidis, 2003), and compare very favourably with histopathological diagnosis accuracy in the field of thyroid tumours (Baloch et al, 2001;Hegedus, 2004;Clary et al, 2005). The stability of the gene lists uncovered through supervised classification is problematic, even in studies using hundreds of samples (Ein-Dor et al, 2005;Michiels et al, 2005). Clearly, much larger studies will be needed to list exactly and exhaustively the discriminating genes, and validate them over a larger group. Nevertheless, our results strongly suggest that such genes exist: accurate classification of CTB and French tumours is possible on the basis of their expression profiles. Initial reports of a low BRAF mutation frequencies in post-Chernobyl tumours (Nikiforova et al, 2004) and of a large impact of BRAF on gene expression (Giordano et al, 2005) raised the possibility of a radiation damage signature based on the mutational status of the tumours. The frequency of BRAF V600E mutation was similar, 38%, in our French and CTB tumours. Our analysis does not exclude the possibility of other damage signatures yet to be identified.
Radiation is a proven causing factor for PTC and a number of arguments support the view that H 2 O 2 -induced damage also contributes to initiate these tumours (see Introduction). Taking advantage of published data on the transcriptional responses of human lymphocytes to 13 stress agents (Amundson et al, 2005), we investigated how similar the responses to H 2 O 2 and g-radiation are. We found that among 10 genotoxic agents, H 2 O 2 at 200 mM elicits the response most similar to that of radiation at 2.5 Gy. This strengthens the argument for H 2 O 2 as a PTC-causing agent, as this similarity most probably mirrors a similarity of the damage inflicted by H 2 O 2 and radiation. This similarity is in line with the finding that French and CTB tumours have similar global profiles. Interestingly, Xiong et al (2005) demonstrated that the number of chromatid breaks per cell following g-irradiation was significantly higher in the lymphocytes of 57 PTC patients with no documented exposure to radiation than in the lymphocytes of healthy controls. This difference could be related to impaired homologous recombination as the 18067T allele variant of XRCC3 was more frequent in 134 thyroid cancer patients than in 166 healthy patients in another study .
Transcriptional responses to H 2 O 2 and g-radiation are similar relatively to other responses to genotoxic agents. However, 118 genes regulated differently in response to H 2 O 2 and radiation were uncovered and could be used to classify CTB and French tumours with an error as low as 15%. This is straightforward evidence that at least some of the genes associated with these tumours are also associated with the response to their presumed respective aetiological agent.
Next, we investigated whether French and CTB tumours could be classified on the basis of five signatures covering the genes involved in the five major DNA repair mechanisms: base-excision repair, mismatch-excision repair, nucleotide-excision repair, homologous recombination and nonhomologous end joining (Wood et al, 2001(Wood et al, , 2005. The homologous recombination signature, which shares no genes with the H 2 O 2 vs g-radiation signature, led to classification errors ranging from 15 to 31%. None of the other four signatures led to accurate classification. The specificity for the homologous recombination effect, and the good classification of CTB and French tumours using the 118 genes regulated differently in response to H 2 O 2 and radiation, make unlikely the confounding effect of age-or ethnicity-related factors. The fact that homologous recombination is involved in double-strand break repair fits the notion that radiation causes more double-strand breaks than H 2 O 2 . Nevertheless, although potential confounders are controlled for by the use of patient-matched adjacent tissues, they are not formally ruled out in our study. This will become possible in the future as tumours from younger Ukrainian patients born after 1987 become available. Thus, several independent gene-expression signatures separate our CTB and sporadic PTCs. These subtle expression differences between CTB and French tumours must be interpreted in light of the fact that the tumours investigated were removed 415 years after the Chernobyl accident. Thus, any discriminating gene-expression signature had to be sustained over this time interval. DNA damage resulting from radiation, however, is typically mostly repaired within a time scale of hours. Consequently, either the reported signatures are 'damage signatures', that is, they are late results, from radiation-induced DNA damage (e.g. non-or incorrectly repaired damage), and/or they are 'susceptibility signatures', that is, they mirror radiation susceptibility factors pre-existing to the accident. The fact that one of the signature relies on the relative response to the two postulated causing agents (g-rays and H 2 O 2 ), and that the other relies on double-strand break repair genes, suggests that these signatures are related to the tumour-initiating mechanisms. This and the longlasting presence of these signatures support the susceptibility signature model. The recent finding that different TP53 alleles are associated with radiation exposure in adult PTC from Russian-Ukrainian patients (Rogounovitch et al, 2006) also supports this view. The susceptibility model, and the corollary that radiation susceptibility varies among individuals, may partly explain why only a minority of the population most exposed to radiation in Ukraine and Belarus developed PTC.
Thus, we interpret our findings as evidence for different and detectable cancer susceptibility factors underlying CTB and French tumours, which leads to several testable predictions. Expression ratios of tumours with respect to patient-matched adjacent tissues were measured. Hence, we could uncover susceptibility signatures only to the extent that they manifest themselves differently in the cancers and their adjacent tissues. We anticipate that the direct comparison of expression levels instead of expression ratios could lead to a stronger signature, possibly involving more genes. In addition, a radiation susceptibility signature could be present in healthy cells of any type in post-Chernobyl cancer patients. This, then, suggests the possibility of developing an expression-based in vitro test for radiation susceptibility. Finally, large-scale studies could uncover the genetic or epigenetic variations underlying the phenotypic differences reported in this paper. These concepts and approaches may apply to other types of cancers.