
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Human Cancer Biology |
Authors' Affiliations: 1 Department of Radiation Oncology, University of California Davis Cancer Center, Sacramento, 2 Division of Biostatistics, 3 Department of Statistics, University of California Davis, Davis, and 4 Lawrence Livermore National Laboratory, Livermore, California
Requests for reprints: Zelanna Goldberg, Department of Radiation Oncology, University of California Davis Cancer Center, 4501 X Street, Suite G-149, Sacramento, CA 95817. Phone: 916-734-8172; Fax: 916-454-4614; E-mail: zgoldberg{at}ucdavis.edu.
| Abstract |
|---|
|
|
|---|
Experimental Design: In this article, we present the first data on the response of human tissue exposed in vivo to LDIR with precisely controlled and calibrated doses. We evaluated transcriptomic responses to a single exposure of LDIR in the normal skin of men undergoing therapeutic radiation for prostate cancer (research protocol, Health Insurance Portability and Accountability Actcompliant, Institutional Review Boardapproved). Using newly developed biostatistical tools that account for individual splice variants and the expected variability of temporal response between humans even when the outcome is measured at a single time, we show a dose-response pattern in gene expression in a number of pathways and gene groups that are biologically plausible responses to LDIR.
Results: Examining genes and pathways identified as radiation-responsive in cell culture models, we found seven gene groups and five pathways that were altered in men in this experiment. These included the Akt/phosphoinositide-3-kinase pathway, the growth factor pathway, the stress/apoptosis pathway, and the pathway initiated by transforming growth factor-ß signaling, whereas gene groups with altered expression included the keratins, the zinc finger proteins and signaling molecules in the mitogen-activated protein kinase gene group. We show that there is considerable individual variability in radiation response that makes the detection of effects difficult, but still feasible when analyzed according to gene group and pathway.
Conclusions: These results show for the first time that low doses of radiation have an identifiable biosignature in human tissue, irradiated in vivo with normal intact three-dimensional architecture, vascular supply, and innervation. The genes and pathways show that the tissue (a) does detect the injury, (b) initiates a stress/inflammatory response, (c) undergoes DNA remodeling, as suggested by the significant increase in zinc finger protein gene expression, and (d) initiates a "pro-survival" response. The ability to detect a distinct radiation response pattern following LDIR exposure has important implications for risk assessment in both therapeutic and national defense contexts.
Epidemiologic studies of exposed populations have evaluated global population exposure and health outcomes of the group, usually, although not exclusively, malignancy (15). A careful review of these studies reveals that, at the lowest exposure levels, no increase in detectable risk is observed. Whether this is a true threshold or simply an inherent limitation of the epidemiologic tools is the subject of vigorous debate.
The other approach to determining any potential biological activity in response to LDIR has been through classic in vitro cell biology, and more recently, molecular biology using genomics to characterize the response pattern. Classic end points such as cell lethality have not been sensitive enough to show any biological effects from such low dose exposures. Using genome-wide techniques on cell culture models, evidence of transcriptional changes following LDIR exposure has been shown. However, such responses are often subtle and inconsistent across models (6, 7).
Cell culture data of primary human keratinocytes and fibroblasts exposed to low or high doses of ionizing radiation have suggested a different response profile under each of these sets of circumstances (6, 7). Another group has examined brain samples from whole animal exposures and showed a transcriptomic response profile to LDIR that was qualitatively different than that seen at higher doses (8). Each of these studies has identified three distinct groups of transcripts: those that were altered by high- and low-dose exposures, those that were unique to high-dose, and those unique to the low-dose exposures. This suggests the possibility of a low-dose radiation response that is biologically dissimilar to higher doses of radiation.
While these results are intriguing, it remains unknown how applicable they are to the actual human response because the experiments were conducted on single cell layer, single cell type in vitro models, which fail to accommodate the known importance of cell-cell communication in radiation response. Furthermore, data suggesting that transcriptional changes of key response genes, e.g., protein kinase C, were markedly different in sensitivity and isoform when ionizing radiation was delivered in vitro or in vivo in a mouse model. This highlights the necessity for actual human data to begin assessment of the human responses to LDIR (9).
The studies described herein begin to address this information gap. We have developed a scientifically valid and ethically compliant model for direct evaluation of LDIR effects in normal human tissue, as well as a statistically validated tool to evaluate the response. Using genomic evaluation, with a priori selection of transcripts potentially sensitive to radiation-induced changes, we have (a) established the first data set of human responses to well-defined doses of LDIR, (b) developed a new methodology to address the natural heterogeneity of response in an unselected patient population, (c) developed a benchmark data set of transcript responses to LDIR in human skin, the first line of defense to radiation exposure from natural, medical, or terrorist sources, and (d) showed that although the in vitro evaluation of LDIR effects has identified genes and pathways activated by such exposure, it does not characterize the response in tissue. Although in vitro assays do serve an important role in the support of translational radiation response models, they cannot evaluate true radiation response in humans.
| Materials and Methods |
|---|
|
|
|---|
Tissue disruption. Each biopsy sample was loaded into a Lysing Matrix D tube (Qbiogene, Irvine, CA) containing 1 mL of a guanidine thiocyanate solution (5.1 mol/L guanidine thiocyanate, 50 mmol/L sodium citrate, 50 mmol/L EDTA, and 0.5% ß-mercaptoethanol). Each sample underwent three rounds of bead beating in the Fastprep beadbeater instrument at a setting of 6 m/s (model 120A, Qbiogene). Samples were kept on ice between each pulverization step. After tissue disruption, the samples were centrifuged at 12,000 x g for 5 minutes at 4°C. The supernatant was removed and transferred to a new 1.5 mL Eppendorf tube.
RNA extraction and quantitation. An equal volume (1 mL) of phenol chloroform (5:1, pH 4.7) containing 0.5% n-lauroylsarcosine was added to each sample and incubated for 5 minutes at room temperature (Sigma, St. Louis, MO). The samples were centrifuged at 12,000 x g for 10 minutes at 2°C. The aqueous phase was separated from the organic phase containing DNA and protein. An additional phenol chloroform (5:1, pH 4.7) extraction was done for further purification. One milliliter of isopropanol and 200 µL of 3 mol/L sodium acetate was added to each sample for precipitation and stored overnight at 20°C. The precipitated RNA samples were cleaned up using RNeasy columns (Qiagen, CA) by following the manufacturer's cleanup procedure. Modifications to the procedure included passing the RNA through the binding column twice to increase the binding efficiency and elution of the RNA with RNA storage solution (Ambion). The Turbo DNA-Free kit was used to remove any residual DNA contamination (Ambion). RNA integrity was verified with the Molecular Devices SpectraMax plus. The RNA samples averaged 1.6 µg of total RNA and the 260:280 nm ratios were between 1.7 and 1.9. The RNA was stored at 80°C.
Amplification and labeling. The RNA was processed for the Affymetrix arrays using their one-cycle in vitro transcription labeling kit (Affymetrix, Santa Clara, CA). Approximately 500 µg of mRNA was reverse-transcribed using a T7-Oligo(dT) promoter primer in the first-strand synthesis. Second-strand cDNA synthesis was done and the resulting double-stranded cDNA was purified. This cDNA served as a template in the in vitro transcription reaction. The in vitro transcription reaction was carried out in the presence of T7 RNA polymerase and a biotinylated nucleotide analogue/ribonucleotide mix. The biotinylated cRNA targets were then cleaned up and fragmented using standard Affymetrix materials and protocols. The cRNA was hybridized to the Human Genome U133 Plus 2.0 arrays.
Data set selection. To address the complexity of the data which would be generated by the Affymetrix GeneChips; a priori, we selected genes and pathways to be examined. These genes or gene pathways were chosen based on the published literature in radiation response, mostly, but not exclusively, at higher doses. To address the variable temporal response between individuals, any gene in which there was evidence in the literature to justify it as a radiation-responsive gene was evaluated, as was that entire gene pathway. This was to capture the data from all of the sample cohorts who were evaluated (biopsies obtained) at a fixed point in time after ionizing radiation exposure, recognizing that the rate of the transcriptional response would vary across the population. Below, we review the chosen genes/pathways and their rationales.
Ionizing radiation has been shown to activate nuclear factor
B through the degradation of inhibitor-
B. Therefore, the signaling pathways upstream of this (extracellular signal-regulated kinase and mitogen-activated protein kinase) were also considered in the pathways analysis (11). Cell cyclemodulating genes such as cyclins, cdc25 phosphatases, GADD45, ATM, and ATR as well as the chk-1 and chk-2 transcripts were also examined. The TP53 pathway, with upstream and downstream signaling and effector transcripts, was examined, as was transforming growth factor-ß and the ubiquitons (12, 13). Growth factors, such as basic fibroblast growth factor, and its downstream targets (RAS/MEK/mitogen-activated protein kinase/RSK) were assessed, as well as other apoptosis-related paths and survival signaling pathways (Akt/phosphoinositide-3-kinase, growth factors such as tumor necrosis factor, epidermal growth factor receptor, vascular endothelial growth factor, and their downstream effectors Ras, Raf-ERK) were evaluated (1418). DNA repair genes, in particular, the Rad 52 epistasis group, have been reported by some to be radiation-responsive, although others have questioned this relationship given the constituent expression of DNA repair proteins (6, 19). Inflammatory mediators such as cyclooxygenase-2, prostaglandin E2, and interleukins have also been reported (20). We also selected some tissue-specific genes with structural functions (i.e., keratins) given that we anticipated some tissue reorganization in response to ionizing radiation. This type of response was reported by Ding et al. with increased expression of cytoskeleton components, ANLN and KRT15, from normal human fibroblasts (6). DNA damage/remodeling via topoisomerase and zinc finger proteins has also been reported and hence these transcripts were also examined (21) along with genes involved in cell-cell signaling, such as GRAP2 and GPR51, and cytoskeleton, such as ANLN and KRT15.
Therefore, there were seven preidentified gene pathways analyzed. They include the Akt/phosphoinositide-3-kinase pathway (14, 18, 22, 23), chemokine pathway (9, 24), fibronectin pathway (14, 24, 25), growth factor/insulin pathway (9), inflammation pathway (9, 14, 20, 22, 24, 26), stress/apoptosis pathway (9, 13, 14, 20, 2224, 27), and the transforming growth factor-ß/cyclin/ubiquitin pathway (11, 13, 14, 18, 22, 24, 25, 28). There were 19 gene groups analyzed (see Table 1 ).
|
The analysis looked for dose-response patterns that were linear in dose either on the original cGy scale, or using a modified log dose (MLD) in which the positive exposures are coded as 0, 1, and 2, and the 0 exposure is coded as 1, as if it were 0.1 cGy. Because we expected (and found) variability between individuals as to timing of response on pathways and on the type of transcript within classes of genes, we did an analysis in which we aggregated results across individuals and transcripts in order to amplify possibly weak signals, without inducing false-positives.
Full details are given in Rocke (31) and Rocke et al. (32), but in summary, for each probe set and each individual, we conducted a linear regression and computed the t statistic testing for whether there is up-regulation, down-regulation, or neither. The t score will be positive if there is a trend towards up-regulation, and negative if there is an opposite trend. We examined the entire collection of t scores for a given gene group or pathway to determine the possibility of a weak signal in the direction of up-regulation or in the direction of down-regulation. We did this by testing the hypothesis that the collection of t scores had a mean of zero using the one-sample t test or a median of zero using the Wilcoxon rank-sum test. If this was rejected, it indicated that there was a trend towards up-regulation with dose of the gene group or pathway if the t scores are biased in a positive direction or towards down-regulation if the t scores are biased in a negative direction. Because of possible correlations in these test statistics, the most reliable P values are provided by resampling by repeatedly sampling random groups of transcripts of the same size. The empirical P value is the fraction of cases in the resampled gene sets in which the test statistic was more significant than the actual test statistic from the gene group or pathway. These were always two-sided P values.
| Results |
|---|
|
|
|---|
Tables 1 and 2 show the results for the preidentified gene groups and pathways. Of the 19 gene groups, 7 have empirical P < 0.0125 on either the regression on dose or the regression on MLD, using either the t test or the Wilcoxon test, a number far exceeding chance (the expected number out of 19 is <1, and the probability of getting 7 significant by chance is <0.0003). The largest effects are for the BCL-6 group, keratins, protein disulfide isomerases, S100, and zinc finger proteins.
|
0.00006. The effect is strongly exemplified by transforming growth factor-ß/cyclin/ubiquitin pathway activated in response to stress/apoptosis, inflammation, growth factor/insulin, and Akt/phosphoinositide-3-kinase pathways. Of the seven preidentified gene pathways analyzed, five showed a significant dose-response up-regulation to LDIR exposure. The two that did not were the chemokine and fibronectin pathways. There were 19 gene groups analyzed, of which seven showed radiation-responsive differential gene expression. The groups that were up-regulated include BCL6 (14), cytokines (9, 23, 25, 28), mitogen-activated protein kinases (14, 20, 23, 24, 35), and zinc finger proteins (27). Those that were down-regulated include the keratins (25), the protein disulfide isomerases (15), and S100 (36). Those that did not change significantly included BAG (37), BCL2 (38), BMPs (39, 40), BMP receptors (40), cyclins (41, 42), GADD45 (43), heat shock (41, 44), interleukins (45), RAD51 (19), serine/threonine kinases (46, 47), tumor necrosis factor-ß (48), and topoisomerases (21, 49). The radiation-responsive gene groups and pathways identified in our human biopsy studies have key regulatory roles in cell cycle arrest, DNA repair, inflammation, and apoptosis.
In most cases, gene groups and pathways that were significant over all patients were also significant for individual patients, although this did not always occur, which is consistent with the assumption of interindividual variability in response. In many cases, the individual tests were not significant, consistent with the lower statistical power of a test based on one individual, and in only a few cases was the effect significant in the opposite direction. This overall consistency gives additional confidence in the realism of the results. Details are given in Table 3 .
|
The microarray data from this study are available in Minimum Information About a Microarray Gene Experimentcompliant format from the National Center for Biotechnology Information Gene Expression Omnibus repository and include the 32. CEL files and the probe set summary data for all 32 arrays and 54,675 probe sets in a Microsoft Excel file, as well as the experimental metadata. The data sets defining the gene groups and pathways, and the programs in the R language that were used to process the data are available as supplementary files from the Journal web site. All of these files are also available on the second author's web site, http://www.idav.ucdavis.edu/~dmrocke.
| Discussion |
|---|
|
|
|---|
By devising a dosimetrically sound sampling method on individuals who are otherwise healthy but receiving therapeutic localized radiation for early stage prostate cancer, we have developed a new model system whereby direct evaluation of radiation effects can be determined. This model is ethically sound, Institutional Review Boardapproved, and allows real-time sampling of human tissue after in vivo radiation exposures of precisely calibrated radiation doses. This model is therefore unique in allowing the evaluation of whole tissue effects when exposure is under normal physiologic conditions.
Biopsies were obtained from patients 3 hours after an acute exposure to LDIR. This time point was chosen for practical considerations. More recently, emerging data have shown that low dose-response may be better defined at a later time point (6, 7, 35). Although none of the published data have examined an in vivo human model, it is reasonable to question whether a greater response may be seen after a longer interval. If so, this would suggest that secondary tissue effects, i.e., cellular responses to signaling from neighboring cells may be as important in tissue response to LDIR as the initial radiation deposition within the cells. Studies to better characterize the temporal response to LDIR in tissue, in a separate cohort of men, are ongoing and will be reported separately.
Using genome-wide analysis, preselected targeted genes, and pathways from in vitro studies, we have shown that human skin initiates a transcript program to enhance survival whereas simultaneously undergoing an inflammatory and stress response. Whether these pathways are all activated in the same cells within the tissue, or represent responses from different tissues (i.e., epithelial versus stromal versus vasculature) cannot be determined from this data set. Laser capture dissection of the biopsied tissue for separate analysis of each cell type can be undertaken now that benchmark data have been compiled showing that low-dose ionizing radiation is biologically active in a definable way in the human. These data were analyzed on the assumption of a linear dose-response. Nonlinearity at the lowest doses cannot be excluded.
The tissues sampled and examined were evaluated as whole tissue. Therefore, the genomic signature is a summative evaluation of the responses of thousands of cells of several lineages, containing both epithelial and stromal cells. Some of the cells would represent lineages that are classically "radiation-sensitive" whereas others were radiation-resistant. Although this model does not allow precise mechanistic intervention, it does offer a comparative advantage in giving us tissue level response which will be the defining level of human responsei.e., the tissue level response is more than just the summative whole of individual cell responses. The doses examined in this study, and those of interest to the radiation risk community are well below those that are frankly cytotoxic. Therefore, the tissue response will inevitably be an amalgamation of the primary response of the cell and the modulating effects of the cells surrounding it, both of the same lineage and others. This amalgamated response therefore incorporates the bystander response and represents the organismal response to the radiation insult. Aside from better determining the extent of true human response to low doses of ionizing radiation, it also provides a potential biosignature for evaluating skin exposure to low doses of ionizing radiation, well below doses expected to produce symptomatic radiation effects.
We have devised a methodology whereby an individual biopsy from these same sample types can be processed for both RNA and protein with both collected for analysis. This will allow a matched genomic and proteomic evaluation.5 These evaluations are ongoing, but will further extend our ability to define the tissue level response of human skin to low doses of ionizing radiation.
Defining the effects of low-dose radiation in humans requires a model with diverse genetic background with the subsequent statistical tool development to identify the effect of the agent above the background diversity of the study population. Therefore, we did these benchmark studies on an otherwise unselected population and devised statistically rigorous tools to address human complexity. Limitations of a sample set consisting of males with early stage prostate cancer are a clear lack of gender diversity and a limited age range. Radiation sensitivities in other populations, specifically the young and fetus' in utero, cannot be extrapolated from these data. This data set included samples from six Caucasians and two African-Americans.
The skin forms the first barrier of the body to radiation exposure and is an easy tissue to biopsy for scientific evaluation, clinical studies and potential field application for population screening for radiation exposure. It is a dynamic, complex system with multiple cell types. Because only healthy skin was biopsied (not previous scars or cancerous lesions), this radiation profile may represent a generalizable biosignature. The age distribution within the patient population sampled does not allow an analysis of age effects, and, as noted above, gender effects cannot be excluded. Relating this transcriptomic signature to a functional tissue end point for human health is part of an ongoing project.
These benchmark data show for the first time that low doses of radiation have an identifiable biosignature in human tissue, irradiated in vivo with normal intact three-dimensional architecture, vascular supply, and innervation. The genes and pathways show that the tissue (a) detects the injury, (b) initiates a stress/inflammatory response, (c) undergoes DNA remodeling, suggested by the significant increase in genes for zinc finger proteins being expressed, and (d) initiates a "pro-survival" response.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
5 Berglund S.R., Schwietert C.W., Jones A., Stern R.L., Lehmann J., and Goldberg Z. Optimized methodology for dual extraction of RNA and protein from small human skin biopsy, submitted for publication. ![]()
Received 11/30/05; revised 2/18/06; accepted 3/ 6/06.
| References |
|---|
|
|
|---|
rays induces specific gene regulations in normal human keratinocytes. Radiat Res 2005;163:62335.[CrossRef][Medline]
irradiation in vivo and ex vivo. Radiat Res 2003;159:4537.[CrossRef][Medline]
B through two distinct mechanisms. Proc Natl Acad Sci U S A 1998;95:130127.
radiation. Radiat Res 2001;156:3247.[CrossRef][Medline]
-irradiation-induced endothelial cell apoptosis, and vascular endothelial growth factor protects endothelial cells through the phosphoinositide 3-kinase-Akt-Bcl-2 pathway. J Biol Chem 2004;279:4335260.
B in radioresistance of TP53-inactive human keratinocytes. Cancer Res 2002;62:121321.
Irradiation up-regulates expression of B cell differentiation molecule CD23 by NF-
B activation. J Biochem Mol Biol 2004;37:50714.[Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |