|
|
|
|
Published online before print
August 12, 2003, 10.1101/gr.1240103 Genome Res. 13:2092-2100, 2003 ©2003 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/03 $5.00
Letter Transcriptional Response of Lymphoblastoid Cells to Ionizing RadiationDepartments of Pediatrics and Genetics, University of Pennsylvania, The Children's Hospital of Philadelphia, Philadelphia, Pennsylvania 19104, USA
The effects of ionizing radiation (IR) on the temporal transcriptional response of lymphoblastoid cells were investigated in this study. We used oligonucleotide microarrays to assess mRNA levels of genes in lymphoblastoid cells at various time points within 24 h following -irradiation. We
identified 319 and 816 IR-responsive genes following 3 Gy and 10 Gy of IR
exposure, respectively, with 126 genes in common between the two doses. A high
percentage of IR-responsive genes are involved in the control of cell cycle,
cell death, DNA repair, DNA metabolism, and RNA processing. We determined the
temporal expression profiles of the IR-responsive genes and assessed effects
of IR dose on this temporal pattern of expression. By combining dose-response
data with temporal profiles of expression, we have identified sets of
coordinately responding genes. Through a genomic approach, we characterized a
set of genes that are implicated in cellular adaptation to IR stress. These
findings will allow a better understanding of complex processes such as
radiation-induced carcinogenesis and the development of biomarkers for
radiation exposure.
Ionizing radiation (IR) is used extensively in medical diagnostic and treatment protocols. It is also present at low levels throughout the environment. IR poses a major threat to cells by compromising genomic integrity and cellular viability. The increased risk for developing malignancies associated with IR exposure has been well documented (for review, see Bast and Gansler 2000 The manner by which IR damages cells is dependent on the type of radiation. Electromagnetic radiation ionizes cellular components indirectly through the generation of highly reactive free radicals, whereas protons and other heavy particles are direct ionizing agents. The major consequence of IR exposure is the generation of single or double-stranded breaks in DNA, which result in a cascade of events involving a complex network of signal transduction and transcriptional regulation. Damage to DNA elicits a cellular stress response that includes DNA damage recognition and cell cycle arrest, followed by DNA repair or apoptosis. If any of these processes fail, mutations can accumulate in the genome, resulting in malignant transformation of somatic cells or heritable mutations in germ cells.
Previous studies have utilized microarrays to describe gene expression
changes associated with IR stress (Amundson et al.
1999
Transcriptional Profile Analysis of Lymphoblastoid Cells Exposed to 3 GY and 10 GY of Ionizing Radiation We irradiated lymphoblastoid cells from 10 unrelated individuals at 2 different doses, 3 Gy and 10 Gy. Cells were harvested prior to irradiation (0 h) and at 1, 2, 6, 12, and 24 h after IR exposure. For each time point and dose, RNA was extracted from each cell line and the samples from all 10 individuals were pooled. A separate reference sample consisting of unirradiated lymphoblastoid cells from a different set of 10 individuals was also prepared. Pooled RNA samples were labeled and hybridized onto Affymetrix Human GeneChips (U95A).
Results from the 3 Gy and 10 Gy experiments were analyzed separately. For
each gene, changes in expression levels were calculated by comparing the
signal intensity of the irradiated sample with that of the unirradiated
reference sample. Genes were considered IR responsive if their expression
levels postirradiation were significantly different from the reference
baseline as determined by the criteria modeled by d-Chip
(Li and Wong 2001 To assess the reproducibility of the microarray-based expression analysis, we compared the results from three independent hybridizations of the same RNA sample onto U95A GeneChip arrays. Labeling reactions and hybridizations were performed separately. The results from the three hybridizations were highly reproducible. By the same criteria used for analyzing the experimental data, only 11 genes of the total 12,559 genes displayed expression level differences in at least 1 of the 3 possible pair-wise comparisons. Of these genes, only 1 gene displayed greater than twofold expression difference between replicates.
IR-Responsive Genes We validated a group of genes by quantitative RTPCR using a different set of cell lines than those used in the microarray study (Fig. 1). In all cases, the quantitative RTPCR results confirmed the U95A GeneChip data.
Functional Characteristics of IR-Responsive Genes
IR-Responsive Genes in Common Between 3 GY and 10 GY IR Exposure
IR-Responsive Genes Specific for 3 GY Or 10 GY IR Exposure A large number of genes appear to respond specifically to only one of the IR doses. The higher radiation dose elicited transcriptional changes in a larger number of genes. The greater degree of insult caused by a higher radiation dose triggers a more complex response.
Specifically, 10 Gy of IR induced a number of DNA repair genes, including
BLM, ERCC4, NBS1, RAD51C, and XPC, which are not affected at
the 3 Gy dose. In addition, the p53-regulated genes, MDM2 and
PCNA, displayed increased expression levels. A higher radiation dose
also activated many cell death-related genes, including a large group of
anti-apoptotic genes (BAG2, BCL2, BCL2A1, BCL2L2, and
BNIP3). Although interferon transcripts did not show detectable
changes in levels, several transcripts for interferon-inducible proteins along
with those for interferon In contrast, a large number of RNA processing/modification genes show gene expression changes specifically at the 3 Gy dose. These include hnRNPs (HNRPC, HNRPD, HNRPH3, HNRPM, HNRPR, and HNRPU) and splicing factors (SFRS1, SFRS3, SFRS6, and SFRS10), which display a gradual decline in transcript levels through the 24-h period following IR exposure.
Irradiation Dose-Dependent Gene Expression Patterns Figure 3A shows that the temporal pattern of expression for certain genes remains quite similar between the two IR doses. The 10 Gy IR dose caused a larger magnitude of expression change for some of the genes, whereas other genes showed identical transcriptional profiles following either IR dose. Most of these genes with similar profiles at both doses are late IR-responsive genes, which exhibit gradual induction or repression kinetics.
Conversely, there are genes that display extremely different temporal
expression patterns at varied IR doses
(Fig. 3B). The majority of
these genes are early responders to IR damage. At the 3 Gy dose, these
IR-responsive genes exhibit transient expression changes that peak at 2 h
postirradiation. In contrast, at the 10 Gy dose, these genes exhibit more
rapid gene expression changes (occurring at 1 h postirradiation) that last
longer in duration. Experiments using quantitative RTPCR showed that
return to basal expression levels of some genes is reached after
Coordinated Expression Profiles of IR-Responsive Genes
To identify sets of genes that are highly correlated in their expression
patterns at both doses, we grouped the 126 IR-responsive genes that are in
common between the doses by similarity in their temporal expression patterns
using GeneCluster 2.0 (Golub et al.
1999
Several groups of known coregulated genes were found
(Fig. 5). For example, the
known p53-regulated genes, CDNK1A and GADD45A, clustered
together in both the 3 Gy and 10 Gy experiments
(Fig. 5A). Both of these genes
are important downstream effectors of the p53 pathway and are involved in
similar pathways. The gene for a hypothetical protein HSPC132 also clusters
with CDNK1A and GADD45A, which suggests that
HSPC132 may be regulated by p53, and may play an important role in
the p53 pathway. Similarly, the two tumor necrosis factor superfamily genes
TNFRSF6 (Fas/Apo-1) and TNFRSF10B (DR5) along with two genes
with zinc finger domains (ZNF148 and ALTE) were tightly
clustered together (Fig. 5B).
The tumor necrosis factor superfamily genes are involved in p53-dependent
induction of apoptosis, and ZNF148 is a transcription factor shown
recently to stabilize p53 (Bai and Merchant
2001
We recognize that the likelihood for small groups containing three or four genes to cluster together by chance is high. Therefore, we examined our data to identify larger clusters that are less likely to occur by chance. Within our SOM data, there are four groups of five genes, one group of six genes, and one group of nine genes that exhibit coordinated transcriptional IR response. Examples are shown in Figure 5, C and D.
DISCUSSION
The main effect of IR on cells is manifested as genotoxic stress resulting
from damaged DNA. Part of the cellular response involves stabilization of p53
protein. This increase in p53 protein levels then causes the induction of many
genes including ACTA2, CDKN1A, DDB2, FDXR, GADD45A, PIG3, TNFRSF6,
and TNFSF10B (Amundson et al.
1999 Comparison of gene expression profiles for the 3 Gy and 10 Gy IR-responsive genes revealed pertinent characteristics of the pathways involved in response to IR stress. The IR-responsive genes displayed either an early response to IR damage within the first 2 h after IR exposure or a late response that does not become apparent until after 6 h or more following IR treatment. A number of early IR-responsive genes exhibited marked differences in temporal expression pattern as a result of increased IR dose. Most of these genes displayed transient expression changes that peaked 2 h following 3 Gy of IR exposure and returned to basal levels by 6-h postirradiation. However, at 10 Gy, the same genes responded even faster by showing a rapid increase in gene expression at 1 h postirradiation, and the expression levels remained high throughout the 24-h time course. In addition, a higher radiation dose caused a larger magnitude of change in gene expression, although the increase in magnitude was not proportional to the increase in IR dose. The higher radiation dose generated more DNA damage, which likely caused a more rapid response by the early IR-responsive genes. Functionally, several early IR-responsive genes act as immediate effectors of check-point control and activators of repair and apoptotic pathways. Many of these genes are primary targets of p53, and others may be direct targets of signaling pathways involved in IR damage detection. These genes are crucial in the first steps of properly dealing with IR stress. Therefore, early IR-responsive genes must be quick to react to IR exposure in order to prevent propagation of the detrimental effects of IR. On the contrary, most late IR-responsive genes exhibited similar temporal gene expression patterns even with varied IR dose. These genes are most likely downstream targets of early IR-responsive effectors. A large number of these late responders were down-regulated gradually during the 24-h time course, which may be due to the shutdown of various pathways as a result of cell death. Other late IR-responsive genes are cell cycle-related or cell cycle-regulated genes. The apparent increase or decrease in expression levels of some of these genes is due to cell cycle arrest and synchronization of the cell population following DNA damage. We also found sets of genes whose expression profiles are highly correlated at the 3 Gy and 10 Gy doses. Although the pathways regulating the expression profiles of these genes are unknown, their coordinated responses suggest coregulation through common regulatory elements.
Specifically, following 10 Gy of IR treatment, several MAP kinase and MAP
kinase-related genes are transcriptionally induced. MAP kinases have been
implicated in a variety of biological responses, one of which is
stress-induced apoptosis. The MAP kinase pathway is comprised of three
distinct components, ERK, JNK, and p38. JNK and p38 are stimulated by various
stress and pathogenic insults, whereas ERK responds to mitogenic and
differentiation signals (for review, see
Herr and Debatin 2001
Many downstream targets of interferon were induced gradually following IR
treatment, even though interferon transcripts did not show detectable changes
in transcript levels. It is known that in response to viral infections,
interferon operates through the JAK-STAT pathway to mediate transcriptional
changes in target genes. This results in antiproliferative effects, which help
to suppress viral replication. In the case of IR stress, interferon activity
may be promoting the same effects to prevent propagation of DNA damage. In
addition, some studies have shown that interferon and other cytokines can
sensitize tumor cells to IR damage (Sirota
et al. 1996
Specific to the 3 Gy experiments, RNA processing genes including many
hnRNPs, and splicing factors were gradually down-regulated during the 24-h
time course. There are previous reports of cleavage of some hnRNP proteins as
a result of apoptosis, including IR-induced apoptosis
(Waterhouse et al. 1996 The findings in this study elucidated parts of the intricate network of genes that are involved in the IR-response. A better understanding of the molecular components and pathways involved in cellular IR stress response will improve our understanding of complex processes such as carcinogenesis and radiation sensitivity. Medically, some of the IR-induced genes can be used as molecular markers for IR exposure.
Tissue Culture Ten lymphoblastoid cell lines (Coriell Cell Repositories) from members of the Centre d'Etude du Polymorphisme Humain (CEPH) families were grown at a density of 5 x 105 cells/mL in RMPI 1620 with 15% FBS. These individuals are not known to be related. Equal numbers of males and females were chosen. Cells were irradiated at 3 Gy and 10 Gy in a 137Cs irradiator. Cells were harvested prior to irradiation (0 h) and at 1, 2, 6, 12, and 24 h after IR exposure. Total RNA was extracted using RNeasy mini-kit (QIAGEN). A reference sample consisting of unirradiated cells from a different set of 10 CEPH individuals was also prepared in the same manner.
Probe Preparation and Hybridization
Data Analysis
Quantitative RTPCR
We thank Drs. Gerd A. Blobel, Greg H. Enders, Warren J. Ewens, Alan M. Gewirtz, Haig H. Kazazian Jr., and Richard S. Spielman for their comments and suggestions. K.Y.J. is supported by NIH training grant GM08216. This work was supported by NIH grant DC00154 and the W.W. Smith Endowed Chair (V.G.C.). The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
[Supplemental material available online at www.genome.org.] Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.1240103. Article published online before print in August 2003.
1 Corresponding author.
Amundson, S.A., Bittner, M., Chen, Y., Trent, J., Meltzer, P., and Fornace Jr., A.J. 1999. Fluorescent cDNA microarray hybridization reveals complexity and heterogeneity of cellular genotoxic stress responses. Oncogene 18: 36663672.[CrossRef][Medline] Amundson, S.A., Do, K.T., Shahab, S., Bittner, M., Meltzer, P., Trent, J., and Fornace Jr., A.J. 2000. Identification of potential mRNA biomarkers in peripheral blood lymphocytes for human exposure to ionizing radiation. Radiat. Res. 154: 342346.[CrossRef][Medline]
Amundson, S.A., Lee, R.A., Kock-Paiz, C.A., Bittner, M.L., Meltzer,
P., Trent, J.M., and Fornace Jr., A.J. 2003. Differential
responses of stress genes to low dose-rate Ashburner, M., Ball, C.A., Blake, J.A., Botstein, D., Butler, H., Cherry, J.M., Davis, A.P., Dolinski, K., Dwight, S.S., Eppig, J.T., et al. 2000. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 25: 2529.[CrossRef][Medline]
Asher, G., Lotem, J., Cohen, B., Sachs, L., and Shaul, Y.
2001. Regulation of p53 stability and p53-dependent apoptosis by
NADH quinone oxidoreductase 1. Proc. Natl. Acad. Sci.
98:
11881193.
Bai, L. and Merchant, J.L. 2001. ZBP-89 promotes
growth arrest through stabilization of p53. Mol. Cell.
Biol. 21:
46704683. Bast, R.C. and Gansler, T.S. 2000. Cancer medicine e.5. Decker, Hamilton, Ont.; Lewiston, NY. Bulavin, D.V., Demidov, O.N., Saito, S., Kauraniemi, P., Phillips, C., Amundson, S.A., Ambrosino, C., Sauter, G., Nebreda, A.R., Anderson, C.W., et al. 2002. Amplification of PPM1D in human tumors abrogates p53 tumor-suppressor activity. Nat. Genet. 31: 210215.[CrossRef][Medline]
Cai, Y., Zhang, C., Nawa, T., Aso, T., Tanaka, M., Oshiro, S.,
Ichijo, H., and Kitajima, S. 2000. Homocysteine-responsive ATF3
gene expression in human vascular endothelial cells: Activation of c-Jun
NH(2)-terminal kinase and promoter response element.
Blood 96:
21402148. Chen, B.P., Wolfgang, C.D., and Hai, T. 1996. Analysis of ATF3, a transcription factor induced by physiological stresses and modulated by gadd153/Chop10. Mol. Cell. Biol. 16: 11571168.[Abstract]
Embree-Ku, M., Venturini, D., and Boekelheide, K.
2002. Fas is involved in the p53-dependent apoptotic response to
ionizing radiation in mouse testis. Biol. Reprod.
66:
14561461.
Fiscella, M., Zhang, H., Fan, S., Sakaguchi, K., Shen, S., Mercer,
W.E., Vande Woude, G.F., O'Connor, P.M., and Appella, E. 1997.
Wip1, a novel human protein phosphatase that is induced in response to
ionizing radiation in a p53-dependent manner. Proc. Natl. Acad.
Sci. 94:
60486053.
Golub, T.R., Slonim, D.K., Tamayo, P., Huard, C., Gaasenbeek, M.,
Mesirov, J.P., Coller, H., Loh, M.L., Downing, J.R., Caligiuri, M.A., et al.
1999. Molecular classification of cancer: Class discovery and
class prediction by gene expression monitoring.
Science 286:
531537.
Gong, B. and Almasan, A. 2000. Apo2 ligand/TNF-related
apoptosis-inducing ligand and death receptor 5 mediate the apoptotic signaling
induced by ionizing radiation in leukemic cells. Cancer
Res. 60:
57545760.
Gruninger, L., Cottin, E., Li, Y.X., Noel, A., Ozsahin, M., and
Coucke, P.A. 1999. Sensitizing human cervical cancer cells In
vitro to ionizing radiation with interferon
Herr, I. and Debatin, K.M. 2001. Cellular stress
response and apoptosis in cancer therapy. Blood
98:
26032614. Hwang, P.M., Bunz, F., Yu, J., Rago, C., Chan, T.A., Murphy, M.P., Kelso, G.F., Smith, R.A., Kinzler, K.W., and Vogelstein, B. 2001. Ferredoxin reductase affects p53-dependent, 5-fluorouracil-induced apoptosis in colorectal cancer cells. Nat. Med. 7: 11111117.[CrossRef][Medline] Kannan, K., Amariglio, N., Rechavi, G., Jakob-Hirsch, J., Kela, I., Kaminski, N., Getz, G., Domany, E., and Givol, D. 2001. DNA microarrays identification of primary and secondary target genes regulated by p53. Oncogene 20: 22252234.[CrossRef][Medline]
Khodarev, N.N., Park, J.O., Yu, J., Gupta, N., Nodzenski, E.,
Roizman, B., and Weichselbaum, R.R. 2001. Dose-dependent and
independent temporal patterns of gene responses to ionizing radiation in
normal and tumor cells and tumor xenografts. Proc. Natl. Acad.
Sci. 98:
1266512670.
Li, C. and Wong, W.H. 2001. Model-based analysis of
oligonucleotide arrays: Expression index computation and outlier detection.
Proc. Natl. Acad. Sci.
98:
3136. Li, J., Yang, Y., Peng, Y., Austin, R.J., van Eyndhoven, W.G., Nguyen, K.C., Gabriele, T., McCurrach, M.E., Marks, J.R., Hoey, T., et al. 2002. Oncogenic properties of PPM1D located within a breast cancer amplification epicenter at 17q23. Nat. Genet. 31: 133134.[CrossRef][Medline]
Liang, G., Wolfgang, C.D., Chen, B.P., Chen, T.H., and Hai, T.
1996. ATF3 gene. Genomic organization, promoter, and regulation.
J. Biol. Chem. 271:
16951701.
McKinney, L.C., Aquilla, E.M., Coffin, D., Wink, D.A., and
Vodovotz, Y. 2000. Ionizing radiation potentiates the induction
of nitric oxide synthase by interferon-
Nikitina, E.Y. and Gabrilovich, D.I. 2001. Combination
of Nishioka, A., Ogawa, Y., Kubonishi, I., Kataoka, S., Hamada, N., Terashima, M., Inomata, T., and Yoshida, S. 1999. An augmentation of Fas (CD95/APO-1) antigen induced by radiation: Flow cytometry analysis of lymphoma and leukemia cell lines. Int. J. Mol. Med. 3: 275278.[Medline]
Reap, E.A., Roof, K., Maynor, K., Borrero, M., Booker, J., and
Cohen, P.L. 1997. Radiation and stress-induced apoptosis: A role
for Fas/Fas ligand interactions. Proc. Natl. Acad.
Sci. 94:
57505755.
Schmidberger, H., Rave-Frank, M., Lehmann, J., Schweinfurth, S.,
Rehring, E., Henckel, K., and Hess, C.F. 1999. The combined
effect of interferon Sharpless, N.E. and DePinho, R.A. 2002. p53: Good cop/bad cop. Cell 110: 912.[CrossRef][Medline] Sheard, M.A. 2001. Ionizing radiation as a response-enhancing agent for CD95-mediated apoptosis. Int. J. Cancer 96: 213220.[CrossRef][Medline]
Sheikh, M.S., Carrier, F., Papathanasiou, M.A., Hollander, M.C.,
Zhan, Q., Yu, K., and Fornace Jr., A.J. 1997. Identification of
several human homologs of hamster DNA damage-inducible transcripts. Cloning
and characterization of a novel UV-inducible cDNA that codes for a putative
RNA-binding protein. J. Biol. Chem.
272:
2672026726.
Sheikh, M.S., Burns, T.F., Huang, Y., Wu, G.S., Amundson, S.,
Brooks, K.S., Fornace Jr., A.J., and el-Deiry, W.S. 1998.
p53-dependent and -independent regulation of the death receptor KILLER/DR5
gene expression in response to genotoxic stress and tumor necrosis factor
Sirota, N.P., Bezlepkin, V.G., Kuznetsova, E.A., Lomayeva, M.G.,
Milonova, I.N., Ravin, V.K., Gaziev, A.I., and Bradbury, R.J.
1996. Modifying effect in vivo of interferon Takekawa, M., Adachi, M., Nakahata, A., Nakayama, I., Itoh, F., Tsukuda, H., Taya, Y., and Imai, K. 2000. p53-inducible wip1 phosphatase mediates a negative feedback regulation of p38 MAPK-p53 signaling in response to UV radiation. EMBO J. 19: 65176526.[CrossRef][Medline]
Tamayo, P., Slonim, D., Mesirov, J., Zhu, Q., Kitareewan, S.,
Dmitrovsky, E., Lander, E.S., and Golub, T.R. 1999. Interpreting
patterns of gene expression with self-organizing maps: Methods and application
to hematopoietic differentiation. Proc. Natl. Acad.
Sci. 96:
29072912.
Tusher, V.G., Tibshirani, R., and Chu, G. 2001.
Significance analysis of microarrays applied to the ionizing radiation
response. Proc. Natl. Acad. Sci.
98:
51165121.
van der Houven van Oordt, W., Diaz-Meco, M.T., Lozano, J., Krainer,
A.R., Moscat, J., and Caceres, J.F. 2000. The
MKK(3/6)-p38-signaling cascade alters the subcellular distribution of hnRNP A1
and modulates alternative splicing regulation. J. Cell
Biol. 149:
307316.
Waterhouse, N., Kumar, S., Song, Q., Strike, P., Sparrow, L.,
Dreyfuss, G., Alnemri, E.S., Litwack, G., Lavin, M., and Watters, D.
1996. Heteronuclear ribonucleoproteins C1 and C2, components of
the spliceosome, are specific targets of interleukin 1
Zhao, R., Gish, K., Murphy, M., Yin, Y., Notterman, D., Hoffman,
W.H., Tom, E., Mack, D.H., and Levine, A.J. 2000. Analysis of
p53-regulated gene expression patterns using oligonucleotide arrays.
Genes & Dev. 14:
981993.
Received February 3, 2003;
accepted in revised format June 9, 2003.
This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||