Kim, J-H., Moon, Y. R., Kim, J-S., Oh, M-H., Lee, J-W. and Chung, B. Y. Transcriptomic Profile of Arabidopsis Rosette Leaves during the Reproductive Stage after Exposure to Ionizing Radiation. Radiat. Res. 168, 267–280 (2007).
We attempted to obtain a transcriptomic profile of ionizing radiation-responsive genes in Arabidopsis plants using Affymetrix ATH1 whole-genome microarrays. The Arabidopsis plants were irradiated with 200 Gy γ rays at the early reproduction stage, 33 days after sowing. Rosette leaves were harvested during the postirradiation period from 36 to 49 days after sowing and used for the microarray analysis. The most remarkable changes in the genome-wide expression were observed at 42 days after sowing (9 days after the irradiation). We identified 2165 genes as γ-ray inducible and 1735 genes as γ-ray repressible. These numbers of affected genes were almost two to seven times higher than those at other times. In a comparison of the control and irradiated groups, we also identified 354 differentially expressed genes as significant by applying Welch's t test and fold change analysis. The gene ontology analysis showed that radiation up-regulated defense/ stress responses but down-regulated rhythm/growth responses. Specific expression patterns of 10 genes for antioxidant enzymes, photosynthesis or chlorophyll synthesis after irradiation were also obtained using real-time quantitative PCR analysis. We discuss physiological and genetic alterations in the antioxidative defense system, photosynthesis and chlorophyll metabolism after irradiation at the reproductive stage.
INTRODUCTION
The development of techniques for genetic evaluation of biochemical and physiological responses of plant cells to environmental stress factors has allowed enormous advances in our understanding of stress-inducible signaling mechanisms as well as offering promising genetic targets to generate transgenic plants that can tolerate multiple stresses (1–3). However, few advances have been reported for plant research in the field of radiation biology because of both experimental limitations and political/social restrictions. The lack of large-capacity facilities for irradiation and/or low-dose radiation fields (or greenhouses) for long-term research makes plants unattractive for radiation biology studies. Therefore, it has been unusual for ionizing radiation to be used as a mutagen to broaden the genetic spectrum of major crop plants (4). The mechanism of action by which radiation induces genetic alterations in cells has been elucidated by studies using in vitro systems (5–8). A recent report on radiation-induced DNA damages inducing mutations reinforced the concept that the sensitivity of genes to radiation depends on the base composition and/or the spatial localization of the gene on the chromosome (8). The scientific bases of this theory are universal enough to be applicable to both plant and animal systems. Genome-wide transcript profiling was performed using human cell lines exposed to radiation to elucidate changes in intracellular signaling (9). In that study, 19 genes were identified to be differentially expressed in cell cycle regulation, three in apoptosis, and four in nucleotide excision repair. This information increases our understanding of radiation responses at the gene expression level. In spite of some functional similarities of genes in plant and animal systems, however, differences in the genetic backgrounds have limited the use of expression profiles to identify cell signaling in both plants and animals (10). Accordingly, the need for genome-wide analyses of plant transcripts is increasing in the field of radiation biology. Transcriptional and post-transcriptional regulation to determine radiation-mediated signaling in plant cells has rarely been investigated on a genome-wide scale. In contrast, the physiological/ecological responses in irradiated plants have been investigated extensively (11– 15).
It has generally been presumed that radiation-induced damages to DNA occur through both direct energy deposition into the sugar-phosphate moiety and interactions with reactive oxygen species (ROS), generally hydroxyl radical (OH) produced by water radiolysis (8, 16, 17). In addition, cellular responses to radiation can occur through direct participation of ROS in cell signaling (18, 19) and/or indirect involvement of low-molecular-weight signaling factors released from the reactions of ROS and neighboring cellular components (20, 21). Therefore, ROS are critical mediators of intracellular signaling of radiation responses, and ROS-scavenging (antioxidant) enzyme activities are also crucial factors in propagation of radiation responses. Accordingly, studies on the physiology of stress using radiation may provide information about its ability to confer tolerance to multiple stresses to plants (11, 13, 15). Such studies may help improve the survival of plants under deteriorating environmental conditions, e.g., frequent droughts and floods, and increased exposure to UV radiation and may have applications in the detection and evaluation of radioactive materials near nuclear accident sites (15, 22). Profiling of radiation-induced/repressed transcripts can be used to increase the efficiencies of mutation and screening for desirable phenotypes in plants. To date, genome-wide changes in the expressions of genes in Arabidopsis have been documented after treatments with auxin and brassinosteroid (1), chitin (23), hypoxia (24), and oxidative stress-causing agents such as methyl viologen, Alternaria alternata toxin, 3-aminotriazole and ozone (2). Profiling of candidate genes involved in cell cycle regulation (25, 26) or wax biosynthesis (27) has been also reported using Arabidopsis microarrays. Arabidopsis is an attractive plant model for radiation biology because of its short generation time, genetic simplicity, and availability of whole genome information. Recently, the responses of hundreds of plant genes to radiation were characterized through a short-term global transcription analysis in Arabidopsis, and cis-acting elements in the promoter regions of the radiation-responsive genes were also analyzed (28–35). Most altered genes were involved in cell growth, cell division, development, translation, general metabolism, signal transduction, and stress/ defense response. However, a long-term genome-scale profiling of transcripts responding to radiation in plants at the reproductive stage has not been reported for any plant species, including Arabidopsis.
Massive doses of radiation have been shown to induce physical changes in plants, including enhancement of respiration, increase in ethylene production, induction of enzyme activities, and accumulation of sucrose and specific protein species (12). Photosynthesis and antioxidative defense systems have long been studied in relation to the radiation responses of various plant species (11, 13, 15, 30). These sensitive systems continue to be interesting subjects in plant radiation biology. However, the stress sensitivity of a plant can differ depending on the developmental stage (36, 37), and this is not limited to radiation. We first determined the developmental stage of Arabidopsis showing the most dramatic alterations in the phenotypic traits after exposure to γ rays at a dose of 200 Gy, which induced noticeable phenotypic changes in Arabidopsis seedlings in a preliminary experiment. Interestingly, the early reproduction stage within 4–5 weeks after seeds were sown was the most sensitive to radiation, with the irradiated plants showing inhibition of stem growth, retardation of leaf senescence, accumulation of chlorophyll and carotenoid, and elevation of photosynthetic saturation for light use. To gain more insight into how plants respond to radiation at the genome level, we carried out gene expression profiling using Arabidopsis leaves sampled at five times: 3, 6, 9, 13 and 16 days after the plants were irradiated with 200 Gy γ rays at the early reproduction stage. Affymetrix ATH1 oligonucleotide GeneChip arrays were used to cover the whole genome representing nearly all the nuclear, plastid and mitochondrial genes encoded in Arabidopsis. The goal was to generate an informative transcriptomic profile of plant genes responding to radiation at the reproductive stage through a long-term genome-scale analysis of Arabidopsis genes. The results allowed for the genome-wide identification of genes with significant positive or negative changes in expression (induced and repressed, respectively) in response to radiation in plants at the reproductive stage. This paper presents characteristics of radiation responses during the reproduction stage, identifying the radiation-responsive genes related to photosynthesis and antioxidant enzymes.
MATERIALS AND METHODS
Plant Materials and γ Irradiation
Arabidopsis thaliana ecotype Columbia (Col-0) was used in this study. Seedlings were grown for 33 days after sowing at 22/18°C (day/night) with a 16-h photoperiod in a compound soil mixture (vermiculite:peat moss:perlite = 1:1:1). They were exposed to 200 Gy γ rays (dose rate 50 Gy/h) generated by a 60Co γ irradiator (150 TBq of capacity; AECL, Canada) at the Korea Atomic Energy Research Institute. The irradiated seedlings were placed under the same growth conditions as described above. Then their rosette leaves were detached and immediately frozen in liquid nitrogen at 3, 6, 9, 13 and 16 days after irradiation, which corresponded to 36, 39, 42, 46 and 49 days after sowing, respectively. The leaves were stored at −80°C until RNA isolation.
RNA Preparation and Microarray Analysis
To minimize interplant variability, rosette leaves from a minimum of 10 plants were pooled for each RNA extraction. Total RNA was extracted from the leaves using Trizol (Invitrogen, Carlsbad, CA) and further purified using an RNeasy Mini Kit (Qiagen, Valencia, CA) following the manufacturer's specifications. The amount of RNA was determined by spectrophotometry at 260 nm, and its integrity was assessed by analyzing the ribosomal RNA bands after gel electrophoresis. Conversion of total RNA into double-stranded cDNA, generation of biotin-labeled cRNA from the double-stranded cDNA, and hybridization of the biotin-labeled cRNA with the GeneChip® Arabidopsis ATH1 Genome Array (Affymetrix, Santa Clara, CA) were performed and analyzed as recommended by Affymetrix.
Data Analysis
GenPlex™ v1.8 software (ISTECH Inc., Goyang city, South Korea) was used for data analysis. The MAS5 algorithm was used for expression summary and signal calculation. Global scaling normalization using a GCOS algorithm was performed, and then the normalized data were log2-transformed. Fold change and Welch's t test were applied to select the differentially expressed genes. The difference was more than twofold, and the significance level was 0.05. To better visualize and compare the two methods for differentially expressed genes, a volcano plot was used. The twofold differentially expressed genes were clustered using hierarchical clustering with Pearson correlation as a similarity measure and complete linkage as a linkage method (38–40). The top 200 differentially expressed genes with the greatest fold change were classified into functional subgroups by BINGO, which enables Gene Ontology significance analysis (41). Pathway analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, which is available online ( http://www.genome.jp/kegg/pathway.html).
Real-Time Quantitative PCR Analysis
One microgram of total RNA was reverse-transcribed into cDNA using AccuPower RT Premix (Bioneer, Daejeon, Korea) for 60 min at 42°C using 0.5 μg anchored oligo(dT)18V primers. Real-time PCR was carried out in a 25-μl reaction mixture containing 1 μl template cDNA, 12.5 μl SYBR®Premix Ex Taq™ (2×) (Takara Bio Inc., Otsu, Japan), 2 μl forward primer (10 μM), 2 μl reverse primer (10 μM), and 7.5 μl distilled water, using a Smart Cycler® II System (Cepheid, Sunnyvale, CA). Gene-specific primer sets used are listed in Table 1. The thermal profile consisted of 1 cycle at 95°C for 1 min followed by 45 cycles at 95°C for 5 s, 56°C for 20 s, and 72°C for 15 s. For each run, data acquisition and analysis were done using the Smart Cycler® software (version 2.0b, Cepheid. Then the relative transcript level was determined using the comparative (2−ΔΔCt) method (42, 43). Relative quantification of gene expression using the 2−ΔΔCt method was correlated with the absolute gene quantification obtained using calibration curves (44). Actin2 was used as an endogenous control gene to normalize for differences in the amounts of total RNA.
RESULTS AND DISCUSSION
Despite many recent analyses of radiation-regulated gene expression in various animal cells by cDNA or oligonucleotide microarrays (9, 45–47), there are no genome-wide profiles for the long-term expression of genes in plants exposed to radiation. We attempted to analyze radiation-regulated cellular responses at the transcriptional level in Arabidopsis plants, specifically leaves, during the reproduction stage.
Gene Expression of the Control and Irradiated Groups during the Reproduction Stage
We first aimed to determine which genes are expressed and/or regulated during the reproduction stage in the control and irradiated plants. Leaves were harvested at 3, 6, 9, 13 and 16 days after irradiation, which corresponded to 36, 39, 42, 46 and 49 days after sowing, respectively. Transcript levels were calculated relative to the respective 3-day samples and differentially expressed genes were selected with a twofold change and a significance level of 0.05 after log2 transformation of the normalized data. The 3-day samples were used as a control instead of the 0-day samples to exclude the effects of short-term fluctuations in the expression of genes after γ irradiation. Gene expression levels in leaves decrease in terms of the amounts of transcripts during the reproduction stage (36). In the present study, however, the total number of genes affected increased in the control leaves up to 13 days after irradiation and then decreased at 16 days after the irradiation, while the irradiated leaves showed the maximum expression change at 9 days after the irradiation (Fig. 1). The proportion of genes showing an increased transcript level was generally higher in the irradiated leaves than in the control leaves except for the 13-day samples.
The numbers of differentially expressed genes in the control and irradiated groups were 3347 and 3220 in combination during the postirradiation period, respectively, while only 152 and 207 differentially expressed genes existed in common throughout that period. These common differentially expressed genes were classified into two or three major groups using hierarchical clustering with Pearson correlation as a similarity measure and complete linkage as a linkage method (Figs. 2A and 3A). This result indicates that 6 days after the irradiation or 39 days after sowing can be a turning point in gene expression during the reproduction stage. Finally, the data set clusters suggest that although the majority of genes are either continuously up-regulated or down-regulated, there are also groups showing alternative patterns of expression (26). We identified two groups (Clusters 3 and 6) in the control leaves and three (Clusters 2, 4, and 6) in the irradiated ones (Figs. 2B and 3B).
The 200 differentially expressed genes with the greatest changes from the data used for Fig. 1 were classified into functional subgroups by BINGO, an algorithm for Gene Ontology (GO) significance analysis (41) (Table 2). The GO category (GOC) consisted of biological process (BP), cellular component (CC), and molecular function (MF). Genes with stress/stimulus-responsive GO terms were up-regulated in the control leaves, but cell wall/membrane- and development-related genes were down-regulated (Table 2A). This is probably due to the progress of leaf senescence. In contrast, the irradiated leaves showed up-regulation in a defense response-related GO term (Table 2B), which is in good agreement with previous reports of increased gene expression and activities of antioxidant enzymes in irradiated plants (13, 15, 31). Together with these results, the reduced expressions of genes with rhythmic process-, circadian rhythm-, and auxin stimulus-related GO terms explain the inhibition in stem growth and retardation in leaf senescence we observed in the irradiated leaves (supplementary Fig. S1). Moreover, the pronounced change in expression of the transcription regulator-related genes may imply that a broad spectrum of genes are affected in irradiated leaves. However, the pathway analysis of the differentially expressed genes demonstrated that the expressions of the photosynthesis-related genes are relatively constant in the irradiated leaves during the reproduction stage (Table 3).
Comparison of Gene Expression in Control and Irradiated Leaves during the Reproduction Stage
To investigate effects of γ radiation on genome-wide gene expression during the reproduction stage, transcript levels of genes in irradiated leaves were compared to those in the respective controls at 3, 6, 9, 13 and 16 days after irradiation. Differentially expressed genes with a twofold change were selected, and a significance level of 0.05 after log2 transformation of the normalized data was used. As discussed above for Table 2B, radiation caused a significant change in the expression of numerous genes (Fig. 4). The most striking change in gene expression was observed at 9 days after irradiation, when 2164 and 1735 genes were up- and down-regulated, respectively.
ROS are among the most decisive mediators in intracellular signaling of radiation responses. Cellular responses to radiation can occur through direct participation of ROS in cell signaling (18, 19) and/or indirect involvement of LMW signaling factors released from the reactions of ROS and neighboring cellular components (20, 21). Therefore, we attempted to analyze the expression of genes for ROS-scavenging antioxidant enzymes as crucial factors in the propagation of radiation responses (Table 4A). Some of them were confirmed by real-time quantitative PCR using the same RNA samples.
Transcription levels of cytosolic Cu/Zn-superoxide dismutase (cytCu/Zn-SOD, csd1, accession no. At1g08830), chloroplastic Cu/Zn-SOD (chlCu/Zn-SOD, csd2, accession no. At2g40100), and dehydroascorbate reductase (DHAR, dhar, accession no. At1g19570) genes were higher overall in the irradiated leaves than in the control ones during the postirradiation period, reaching maximum levels at 9 or 13 days after irradiation (Fig. 5A, B and E). In contrast, the expressions of the glutathione peroxidase (GP, gp, accession no. At4g31870) and catalase 2 (CAT2, cat2, accession no. At4g35090) genes remained lower in the irradiated leaves during the same period (Fig. 5F and G). Interestingly, the transcript levels of the chloroplastic Fe-SOD (chlFe-SOD, fsd1, accession no. At4g25100) and stromal ascorbate peroxidase (sAPX, sapx, accession no. At4g08390) genes fluctuated more dramatically, remaining lower in the irradiated leaves until 9 days after irradiation and becoming higher afterward (Fig. 5C and D). Although changes in the expression levels of the genes were not the same in the microarray and real-time PCR results due to the normalization for the former and the comparative expression for the latter using actin2 as an internal control gene, the overall expression patterns of the genes were similar in the two systems. These data are in good agreement with high SOD and low CAT activities in plants after an acute exposure to relatively high doses of radiation (15, 48). However, the distinct expression patterns of antioxidant enzyme genes during the postirradiation period may indicate the existence of a complicated and controlled gene regulation system for radiation-induced gene responses in plants. Similarly, the differential expression patterns for the genes encoding the SOD, APX and CAT isoenzymes have been reported in rice plants after γ irradiation under salt stress (22). The results shown in Table 5 suggest that the differentially expressed genes for ascorbate/glutathione metabolism may contribute to the radiation-induced stress and defense signaling (49). The high gene count for ribosome at 9 days after irradiation can be correlated with the maximum change in gene expression, as shown in Fig. 4.
Distinctive changes in photosynthesis and pigment metabolism of plants after radiation exposure have been documented (11, 15, 50, 51). In particular, a high dose of radiation could cause specific physical alterations in Arabidopsis seedlings at certain developmental stages, such as accumulation of chlorophyll and carotenoid and elevation in photosynthetic saturation for light use as well as inhibition of stem growth and retardation of leaf senescence. Although the expressions of the photosynthesis-related genes are relatively constant in the irradiated leaves during the reproduction stage compared to those at 3 days after irradiation (Table 3), their relative expressions compared to the respective controls showed significant radiation-induced changes in photosynthesis and pigment metabolism (Table 4B, Fig. 5H, I and J). Expression of glutamyl-tRNA reductase 1 (HEMA1, hemA1, accession no. At1g58290), NADPH-protochlorophyllide oxidoreductase A (PORA, porA, accession no. At5g54190), and NADPH-protochlorophyllide oxidoreductase B (PORB, porB, accession no. At4g27440) was increased during the postirradiation period, with maximum levels reached at different times. Differential but increased transcript levels of these genes, which encode two major enzymes for chlorophyll synthesis (52, 53), may contribute to the increased chlorophyll content of the irradiated leaves at 13 and 16 days after irradiation (supplementary Fig. S1). In contrast, genes of chlorophyll A-B binding proteins 2:4 (LHCB2:4, lhcb2:4, accession no. At3g27690) and 4.3 (LHCB4.3, lhcb4.3, accession no. At2g40100) showed opposite expression patterns; the former showed higher transcript levels in the irradiated leaves than in the controls during the postirradiation period while the latter showed the reverse. The alternative expression of lhcb2:4 can be correlated with retardation of leaf senescence in the irradiated leaves during the reproduction stage (supplementary Fig. S1).
Comparison of Gene Expression in the Control and Irradiated Groups throughout the Reproduction Stage
To compare transcript levels of genes in the control and irradiated groups during the postirradiation period, we carried out clustering analysis of 2856 genes below the 0.05 level of significance in Welch's t test (Fig. 6). The genes were divided into two subgroups (Clusters 1 and 2) based on the expression pattern. Expression of Cluster 1 genes (1570 genes) was decreased in the irradiated leaves, while that of Cluster 2 genes (1286 genes) was increased. Furthermore, filtered data for genes with more than a twofold change (533 genes) and changes with less than 0.05 significance in Welch's t test (2856 genes) included 354 differentially expressed genes, of which 225 were up-regulated and 129 down-regulated in the irradiated group (Fig. 7). However, only porB, lhcb4.3, cat2 and lhcb2:4 among the genes listed in Table 4 were selected in common with the 354 differentially expressed genes (supplementary Table S1). The other genes, especially the antioxidant genes confirmed in Fig. 5, were not selected as significant differentially expressed genes in the comparison of the control and irradiated groups. This can be attributed to the transience and fluctuation of their expressions, which changed throughout the postirradiation period. Instead, the involvement of two ribonucleotide reductase genes (accession nos. At2g21790 and At3g27060), which are critical for cell cycle progression, DNA damage repair, and plant development (54), implicated the most well-known genetic response for recovery in the irradiated leaves. This may be involved in the sudden increase of affected genes at 9 days after irradiation (Fig. 4).
In the present study, we showed changes in transcript levels of csd1, csd1, fsd1, sapx, dhar, cat2 and gp in leaves during the 16-day period after γ irradiation of plants at the early reproduction stage. Most of these enzymes belong to ascorbate-glutathione pathway, and they scavenge ROS or its derivatives (55). Interestingly, we found low transcript levels of cat2 and gp in the irradiated leaves. Although these results are partly predictable from our previous studies (15, 22), there are also opposing results with chronic radiation-induced responses (13). Similarly, the expression of cat2 was strongly induced by UV-light pulses (290 to 400 nm) in leaves of maize seedlings (56). Accordingly, the transcription of these genes appears to be inhibited after exposure to an acute dose of radiation. In contrast, the expression patterns of hemA1, porA, porB, lhcb2:4 and lhcb4.3 can be correlated with the physiological changes in the irradiated leaves, e.g., retardation in leaf senescence, accumulation of chlorophyll and carotenoid, and elevation of photosynthetic saturation for light use (Table 4 and supplementary Fig. S1). Together with cat2, however, porB and lhcb4.3 belong to the most distinguishable genes in terms of the expression change throughout the postirradiation period (data not presented).
Recently, Nagata et al. (31) reported the results of microarray analysis of genes that respond to ionizing radiation in Arabidopsis. Their study focused on the short-term effect of ionizing radiation on the expression of genes within 24 h after γ irradiation with a dose of 3 kGy, which makes the subsequent growth and survival of plants impossible. In contrast, we attempted to determine the long-term effects of ionizing radiation on the genome-wide expression of genes in Arabidopsis rosette leaves during the reproductive stage, which is responsible for the phenotypic and physiological changes in the irradiated leaves.
In conclusion, we generated genome-wide transcriptomic profiles of radiation-responsive genes in plants at the reproductive stage using the Arabidopsis genome array and revealed the changes in expression of the genes involving in antioxidative defense, photosynthesis and chlorophyll synthesis. Considering the importance of the ROS scavenging system and pigment metabolism in the radiation-induced responses of plants, our results could allow deeper insight into the long-term changes in the functional genes in plants after radiation exposure. Next we are going to use the expression profile to elucidate transcription factors orchestrating the expressions of radiation-responsive genes, especially those related to antioxidative defense and pigment metabolism.
SUPPLEMENTARY INFORMATION
Supplemental document containing additional information about methodology: http://dx.doi.org/10.1667/RR0963.1.s1 (10.1667_RR0963.1.s1.pdf).
FIG. S1. Phenotypic and physiological differences between the control and irradiated seedlings after γ irradiation: http://dx.doi.org/10.1667/RR0963.1.s2 (10.1667_RR0963.1.s2.pdf).
TABLE S1. Expression profile of 225 induced and 129 repressed genes in the irradiation groups from the 354 differentially expressed genes selected in Fig. 6: http://dx.doi.org/10.1667/RR0963.1.s3 (10.1667_RR0963.1.s3.pdf).
Acknowledgments
This project has been carried out under the Nuclear R&D Program by Ministry of Science and Technology, Republic of Korea.
REFERENCES
TABLE 1
Primer Sequences Used in Quantitative RT-PCR Experiment
TABLE 2
Gene Ontology Analysis of Differentially Expressed Genes in Control (A) and Irradiated (B) Samples during the Reproduction Stage
TABLE 3
Pathway Analysis of the Differentially Expressed Genes during the Reproduction Stage
TABLE 4
Radiation Responsive Genes during the Postirradiation Period
TABLE 4
(continued)
TABLE 5
Pathway Analysis of the Differentially Expressed Genes during the Postirradiation Period