Introduction

Molecular Medicine Reports

1791-2997 1791-3004

D.A. Spandidos

10.3892/mmr.2017.6476

mmr-15-06-3599

Articles

GDF-15 gene expression alterations in human lymphoblastoid cells and peripheral blood lymphocytes following exposure to ionizing radiation

Shuang

1 * Zhang

Qing-Zhao

1 * Zhang

De-Qin

2 * Feng

Jiang-Bin

1 Luo

Qun

3 Lu

Xue

1 Wang

Xin-Ru

4 Li

Kun-Peng

5 Chen

De-Qing

1 Mu

Xiao-Feng

5 Gao

Ling

1 Liu

Qing-Jie

1China CDC Key Laboratory of Radiological Protection and Nuclear Emergency, National Institute for Radiological Protection, Chinese Center for Disease Control and Prevention, Beijing 100088, P.R. China 2Beijing Shijingshan Center for Disease Control and Prevention, Beijing 100043, P.R. China 3Department of Transfusion, Chinese PLA General Hospital, Beijing 100853, P.R. China 4Department of Clinical Laboratory, Second Artillery General Hospital PLA, Beijing 100088, P.R. China 5Department of Radiotherapy, General Hospital of Armed Police Forces, Beijing 100039, P.R. China

Correspondence to: Dr Qing-Jie Liu or Dr Ling Gao, China CDC Key Laboratory of Radiological Protection and Nuclear Emergency, National Institute for Radiological Protection, Chinese Center for Disease Control and Prevention, 2 Xinkang Street, Beijing 100088, P.R. China, E-mail: qjliu@nirp.cn, E-mail: gaolingxinxin@126.com *

Contributed equally

062017

12042017

15 6 3599 3606 04032016 20022017

2017

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

The identification of rapid, sensitive and high-throughput biomarkers is imperative in order to identify individuals harmed by radiation accidents, and accurately evaluate the absorbed doses of radiation. DNA microarrays have previously been used to evaluate the alterations in growth/differentiation factor 15 (GDF15) gene expression in AHH-1 human lymphoblastoid cells, following exposure to γ-rays. The present study aimed to characterize the relationship between the dose of ionizing radiation and the produced effects in GDF-15 gene expression in AHH-1 cells and human peripheral blood lymphocytes (HPBLs). GDF-15 mRNA and protein expression levels following exposure to γ-rays and neutron radiation were assessed by reverse transcription-quantitative polymerase chain reaction and western blot analysis in AHH-1 cells. In addition, alterations in GDF-15 gene expression in HPBLs following ex vivo irradiation were evaluated. The present results demonstrated that GDF-15 mRNA and protein expression levels in AHH-1 cells were significantly upregulated following exposure to γ-ray doses ranging between 1 and 10 Gy, regardless of the dose rate. A total of 48 h following exposure to neutron radiation, a dose-response relationship was identified in AHH-1 cells at γ-ray doses between 0.4 and 1.6 Gy. GDF-15 mRNA levels in HPBLs were significantly upregulated following exposure to γ-ray doses between 1 and 8 Gy, within 4–48 h following irradiation. These results suggested that significant time- and dose-dependent alterations in GDF-15 mRNA and protein expression occur in AHH-1 cells and HPBLs in the early phases following exposure to ionizing radiation. In conclusion, alterations in GDF-15 gene expression may have potential as a biomarker to evaluate radiation exposure.

ionizing radiation growth differentiation factor 15 gene expression dose-effect neutron γ-rays

Introduction

Gene expression is sensitive to environmental factors such as ionizing radiation. Radiation-induced alterations in gene expression profiles in humans and animal models have been analyzed in an effort to identify biomarkers suitable for the accurate estimation of the radiation dose in exposed individuals (1–10). DNA microarrays for gene expression profile analysis have been used as a novel screening method for the identification of novel biomarkers, as they are able to simultaneously assess alterations in the expression of thousands of genes (3,11–17). Quantitative polymerase chain reaction (qPCR) is often performed to validate microarray results. Alterations in the expression of numerous genes have been identified in irradiated cell lines and human peripheral blood lymphocytes (HPBLs) (18–23). Previous studies in mouse and human models have reported that the patterns of alteration in gene expression induced by radiation are specific, durable and can accurately reflect radiation exposure (24,25).

Previous studies have reported that alterations in the expression of certain genes following exposure to radiation are dose- or time-dependent within a specific dose range or duration of treatment (2,4,26–28). Various radiation-responsive genes, including growth arrest and DNA damage-inducible 45 (GADD45), cyclin-dependent kinase inhibitor 1A (CDKN1A), mouse double minute 2 homolog and proliferating cell nuclear antigen (PCNA), have been identified as downstream targets of tumor protein 53 (6,18,29). A previous study demonstrated that the γ-ray-induced GADD45 upregulation in HPBLs is a linear dose-response relationship after 24 and 48 h of 0–3 Gy irradiation (6). In addition, in human blood samples exposed to various γ-ray doses (0.5, 2, 5 and 8 Gy), alterations in the expression of five genes, including CDKN1A, were detected at 24 h post-irradiation (29). However, a specific pattern of radiation-induced alterations in gene expression, that is consistent across cell types, doses and duration of exposure, has yet to be characterized. The analysis of the alterations in the expression of radiation-responsive genes has potential as a strategy to detect radiation exposure; however, its limitations need to be addressed prior to its establishment as a method for the assessment of radiation exposure in clinical practice (30).

Identifying genes and dose-response models with the highest reproducibility is necessary to establish a gene expression signature database for rapid, reliable and high-throughput biodosimetry following radiation exposure (31–35). Gene expression studies in human primary fibroblasts, keratinocytes (36–38) and pulmonary epithelial cells (19) suggested that growth/differentiation factor 15 (GDF-15), which is also a target gene of p53, is a radiation-responsive gene. Radiation-induced alterations in GDF-15 expression are correlated with the time and dose of radiation exposure (19,34). Previous microarray analysis in our laboratory, demonstrated that GDF-15, among other genes, was consistently upregulated in AHH-1 human lymphoblastoid cells following exposure to 0.5, 3 and 8 Gy of ⁶⁰Co γ-rays (39).

In the present study, the relationship between alterations in GDF-15 expression and the dose of ionizing radiation was investigated in AHH-1 cells and HPBLs. Alterations in GDF-15 expression in AHH-1 cells were examined at various time points following exposure to a wide range of ⁶⁰Co γ-ray and neutron radiation doses. Since HPBLs are sensitive to radiation-induced damage and can be easily sampled, they were used to validate the results obtained in AHH-1 cells. Alterations in GDF-15 mRNA expression levels in HPBLs irradiated ex vivo with γ-rays were investigated at various time points following exposure, and using a wide range of radiation doses. In addition, baseline GDF-15 gene expression levels were quantified in HPBLs from healthy adult donors.

Materials and methods <sec> <title>Cell culture

The AHH-1 cell line, which is a human B lymphocyte cell line derived from a 33-year-old Caucasian male and immortalized by Epstein-Barr virus (40), was obtained from American Type Culture Collection (Manassas, VA, USA). AHH-1 cells were cultured in RPMI-1,640 medium (Thermo Fisher Scientific, Inc., Waltham, MA, USA), supplemented with 10% heat-inactivated fetal bovine serum (FBS; HyClone; GE Healthcare Life Sciences, Logan, UT, USA), 2 mM l-glutamine (Thermo Fisher Scientific, Inc.), 100 U/ml penicillin (Sigma-Aldrich; Merck KGaA, Darmstadt, Germany) and 100 mg/ml streptomycin (Sigma-Aldrich; Merck KGaA). The AHH-1 cell line is diploid and its population doubling time ranges between 16 and 19 h. AHH-1 cells were incubated at 37°C in a humidified 5% CO₂ atmosphere. AHH-1 cells in the exponential phase of growth were seeded in flasks at a density of 1×10⁷ cells/ml and prepared for irradiation.

Human peripheral blood samples

Blood samples from 73 healthy adult donors (39 males and 34 females; age, 20–60 years) were obtained for the quantification of baseline GDF-15 mRNA expression levels. The eligibility of the donors was evaluated using questionnaires and regular medical examination. Peripheral blood samples (4 ml) were collected from each subject. Radiation-induced alterations in GDF-15 gene expression were investigated in peripheral blood samples from 3 healthy adult male subjects (aged 26, 29 and 41 years). The subjects had no history of chronic disease, substance abuse, smoking, or toxic chemical exposure. In addition, they had not been exposed to radiation or had a history of viral infections during the months preceding the study.

The present study was approved by the Ethics Committee of the National Institute for Radiological Protection, Chinese Center for Disease Control and Prevention (Beijing, China). Written informed consent was obtained from all human subjects prior to enrollment in the present study. Experiments were conducted according to the principles outlined in the Declaration of Helsinki.

Sample irradiation

Irradiation with ⁶⁰Co γ-rays was performed in the Beijing Radiation Center (Beijing, China). The source radioactivity was 130 TBq. The exposure setup was calibrated by physical measurement using a tissue-equivalent ionizing chamber. The radiation dose rate was calculated using the source radioactivity and the distance between the source and sample: A dose rate of 1 Gy/min corresponded to a source-sample distance of 73.5 cm. The homogeneous irradiation field was 30×30 cm; the samples were placed within a 5 cm-radius circle and the uncertainty of the calibration was 1%. AHH-1 cells were seeded in flasks (1×10⁷ cells/ml) and irradiated. Radiation doses of 0, 1, 3, 5, 8 and 10 Gy, at a dose rate of 1 Gy/min, were used in the dose-response experiment; doses of 3, 5 and 8 Gy, at a dose rate of 0.5, 1 and 2 Gy/min were used to investigate the effect of various dose rates. Irradiated AHH-1 cells were incubated at 37°C for 8–168 h prior to further analysis.

Peripheral blood samples (20 ml) were collected from each subject via venipuncture into Vacutainers (BD Biosciences, Franklin Lakes, NJ, USA). A total of 5 different radiation doses were used (0, 1, 3, 5 and 8 Gy). Blood samples were divided into five plastic centrifuge tubes, one tube per dose group, and irradiated with a single dose of ⁶⁰Co γ-rays at a dose rate of 1 Gy/min. Following irradiation, blood samples were incubated at 37°C for 4–72 h in RPMI-1640 medium of equal volume, supplemented with 10% FBS.

Neutron irradiation was performed at the General Hospital of Armed Police Forces (Beijing, China), using a ²⁵²Cf neutron source. AHH-1 cells were irradiated with 0, 0.4, 0.8, 1.2 and 1.6 Gy, at dose-rate of 0.073 Gy/min. The dose-rate was calibrated using a tissue-equivalent ionizing chamber prior to irradiation. The homogeneous irradiation field was 10×10 cm; the samples were placed within a 5 cm-radius circle and the uncertainty of the calibration was 3%. Irradiated AHH-1 cells were incubated at 37°C for 48 h prior to further analysis.

Reverse transcription-qPCR (RT-qPCR)

Total RNA was extracted from AHH-1 cells using TRIzol^® reagent (Invitrogen; Thermo Fisher Scientific, Inc.) and reverse transcribed into cDNA using a Reverse Transcription system kit (Promega Corporation, Madison, WI, USA), according to the manufacturer's protocol. Individual PCR reactions were carried out in 20 µl samples in a 96-well plate (Applied Biosystems^™; Thermo Fisher Scientific, Inc.), containing 0.5 µl of each forward and reverse primer (10 µM each), 2 µl cDNA, 7 µl distilled water and 10 µl 2X SYBR-Green PCR Master Mix (Applied Biosystems^™; Thermo Fisher Scientific, Inc.). qPCR was performed using the 7500 Fast Real-Time PCR System (Applied Biosystems^™; Thermo Fisher Scientific, Inc.). Amplification was performed under the following conditions: Initial cycle at 95°C for 10 min, followed by 40 cycles at 95°C for 15 sec and at 60°C for 1 min. The following primers were used: GDF-15, forward 5′-GTTAGCCAAAGACTGCCACTG-3′, reverse 5′-CCTTGAGCCCATTCCACA-3′; and β-actin, forward 5′-ACTTAGTTGCGTTACACCCTTTCT-3′ and reverse 5′-GACTGCTGTCACCTTCACCGT-3′. The relative expression levels for each gene were normalized to β-actin and analyzed using the 2^−∆∆Cq method (41) with the 7500 Sequence Detection software. Each group, composed of three parallel samples, was analyzed after irradiation. All experiments were performed three times.

Lymphocytes were isolated from human blood samples using Ficoll-Paque PLUS (GE Healthcare Life Sciences, Chalfont, UK), according to the manufacturer's protocol. Total RNA was extracted from HPBLs using TRIzol^® reagent (Invitrogen; Thermo Fisher Scientific, Inc.) and reverse transcribed into cDNA using a Reverse Transcription system kit (Promega Corporation), according to the manufacturer's protocol. qPCR was performed on cDNA using a TaqMan^® MGB probe (Thermo Fisher Scientific, Inc.). PCR was performed using a standard curve established using 10-fold successively diluted known copy numbers of a GDF-15 standard. The standard was prepared from the pUC57 plasmid containing GDF-15 cDNA (Thermo Fisher Scientific, Inc.). The plasmid was quantified based on the cDNA concentration determined using UV spectroscopy; stock solutions were prepared in Tris-EDTA buffer that ranged in concentration between 10² and 10⁸ cDNA copies/µl. A standard cDNA curve for the range of 10²-10⁸ cDNA copies per reaction was generated by analyzing 2 µl of each dilution in triplicate using a TaqMan^® MGB probe. The oligonucleotide primers for GDF-15 and thermocycling conditions for qPCR were as aforementioned. Individual PCR reactions were carried out in 20 µl samples in a 96-well plate containing GDF-15 primers, 2 µl cDNA, 7 µl distilled water and 1 µl TaqMan^® Gene Expression Master Mix (Applied Biosystems^™; Thermo Fisher Scientific, Inc.). mRNA expression levels (mRNA copy numbers/µl) were calculated based on the standard curve. Each group, composed of three parallel samples, was analyzed after the irradiation. All experiments were performed three times.

Western blot analysis

AHH-1 cells were washed with PBS, lysed in lysis buffer containing 50 mM Tris-HCl (pH 7.4), 150 mM NaCl, 1 mM EDTA, 1% NP-40 and 10 µg/ml protein inhibitor (Roche Diagnostics, Basel, Switzerland), and centrifuged for 20 min at 5,000 × g and 4°C. Supernatants were collected and protein concentrations were determined using the bicinchoninic acid method. Equal amounts of extracted protein (50 µg) were separated by 12% SDS-PAGE and transferred onto nitrocellulose membranes (GE Healthcare Life Sciences, Shanghai, China). Membranes were blocked for 1 h with TBS containing 0.05% Tween-20 (TBST) with 5% non-fat dry milk at room temperature. Subsequently, membranes were incubated at 4°C overnight with rabbit polyclonal anti-GDF-15 (cat. no. BS3872; 1:500; Biogot Technology Co., Ltd., Nanjing, China) and mouse monoclonal anti-β-actin (cat. no. SC-47778; 1:1,000; Santa Cruz Biotechnology, Inc., Dallas, TX, USA) primary antibodies. Membranes were washed three times with TBST prior to incubation for 2 h with horseradish peroxidase-conjugated anti-rabbit IgG (cat. no. ZB-5301; 1:5,000; Beijing Zhongshan Golden Bridge Biotechnology Co., Ltd., Beijing, China) and horseradish peroxidase-conjugated anti-mouse IgG (cat. no. ZB-2305; 1:5,000; Beijing Zhongshan Golden Bridge Biotechnology Co., Ltd.) secondary antibodies at 37°C. Protein bands were visualized using an Enhanced Chemiluminescence kit (GE Healthcare Life Sciences), according to the manufacturer's protocol, and analyzed using the ChemiDoc^™ XRS+system (Bio-Rad Laboratories, Inc., Hercules, CA, USA). The scanned films were semi-quantified using the UVP EC3 Imaging system (UVP, Inc., Upland, CA, USA). GDF-15 protein levels were normalized to β-actin, which was used as the loading control.

Statistical analysis

Statistical analysis was performed using the SPSS software version 16.0 (SPSS Inc., Chicago, IL, USA). The statistical significance of the difference between groups was assessed by one-way analysis of variance (ANOVA). Least significant difference and Student-Newman-Keuls tests were used to compare individual groups following ANOVA. A χ² test was used to compare differences in gene expression levels between sexes and among age groups of healthy adults. A linear regression method was used to establish the dose-response equation. ANOVA was used to analyze the significance of the equation of the dose-response regression curve. Results are presented as the mean+standard error of the mean. P<0.05 was considered to indicate a statistically significant difference.

Results <sec> <title>Exposure to γ-rays alters GDF-15 mRNA expression levels in AHH-1 cells

Dose-dependent alterations in GDF-15 mRNA expression levels of irradiated AHH-1 cells were assessed using RT-qPCR at 8, 12, 24, 48, 72, 120 and 168 h after exposure to 0, 1, 3, 5, 8 and 10 Gy of γ radiation. The results demonstrated that GDF-15 mRNA expression levels were significantly upregulated in a dose-dependent manner at 8–120 h following γ-ray exposure (P<0.05; Fig. 1A). γ irradiation appeared to reach its maximum effect (~14-fold increase compared with the control group) at 10 Gy and 12 h post-irradiation. The dose-response curves describing the relationship between radiation dose and GDF-15 relative expression levels between 8 and 120 h post-exposure were fitted to a linear model (P<0.05; Table I). No dose-response relationship was apparent at 168 h post-irradiation.

To investigate the association between the radiation dose rate and GDF-15 mRNA expression levels, AHH-1 cells were irradiated with 3, 5 and 8 Gy of ⁶⁰Co γ-rays at dose rates of 0.5, 1 and 2 Gy/min. GDF-15 mRNA expression levels were assessed 48 h following γ-ray exposure. The present results suggested that GDF-15 expression depended on the radiation dose, regardless of the dose rate (Fig. 1B).

Exposure to γ-rays alters GDF-15 protein expression levels in AHH-1 cells

GDF-15 protein expression levels in AHH-1 cells were assessed using western blot analysis at 8, 12, 24 and 48 h after exposure to 0, 1, 3, 5, 8 and 10 Gy γ-rays. GDF-15 protein expression levels were significantly upregulated in a dose-dependent manner within a dose range of 1–8 Gy, 8–48 h following irradiation (P<0.05; Fig. 2). γ irradiation appeared to reach its maximum effect (~8-fold increase compared with the control group) at 8 Gy and 12 h post-irradiation. The dose-response curves describing the relationship between radiation dose and GDF-15 protein expression levels between 12 and 48 h following exposure were fitted to a linear model (P<0.05; Table II). No dose-response relationship was apparent at 8 and 72 h post-irradiation.

Exposure to neutron radiation alters GDF-15 mRNA and protein expression levels in AHH-1 cells

Alterations in GDF-15 gene expression in AHH-1 cells were assessed using RT-qPCR, 48 h following exposure to 0, 0.4, 0.8 and 1.6 Gy of ²⁵²Cf neutron radiation. GDF-15 mRNA expression levels were significantly upregulated in a dose-dependent manner 48 h following irradiation (P<0.05; Fig. 3A). Neutron irradiation appeared to reach its maximum effect (~12-fold increase compared with the control group) at 1.6 Gy. The fitted dose-response curve describing the relationship between radiation dose and GDF-15 mRNA expression levels followed a quadratic model, described by the following equation: y=4.85x²-1.80x+1.35 (R²=0.9769, P<0.05), where y represents the GDF-15 mRNA expression levels induced (fold change) by x Gy of neutron radiation.

GDF-15 protein expression levels were assessed in AHH-1 cells using western blot analysis 48 h following exposure to 0–1.6 Gy of neutron radiation. GDF-15 protein levels were significantly increased in a dose-dependent manner (P<0.05; Fig. 3B). Neutron irradiation appeared to reach its maximum effect (~4-fold increase compared with the control group) at 1.6 Gy. The dose-response curve describing the relationship between radiation dose and GDF-15 protein expression levels followed a linear model, described by the following equation: y=2.11x+1.15 (R²=0.954, P<0.05), where y represents the GDF-15 protein expression levels induced by x Gy of neutron radiation.

Baseline GDF-15 mRNA expression levels in peripheral blood samples from healthy adults

To investigate the baseline levels of GDF-15 gene expression in the general population, GDF-15 mRNA levels (mRNA copy number per µl) were assessed in peripheral blood samples, obtained from 73 healthy adult donors. The present results demonstrated that GDF-15 mRNA expression levels in males were not significantly different compared with in females (Fig. 4A). In addition, no significant differences in GDF-15 mRNA expression levels were revealed among different age groups; however, an increasing trend was observed within the 20–50 age range (Fig. 4B). The mean GDF-15 mRNA copy number was calculated and revealed to be 2573.3±2006.9 per µl.

Exposure to γ-rays alters GDF-15 mRNA expression levels in HPBLs

Radiation-induced alterations in GDF-15 gene expression were investigated in human peripheral blood samples obtained from three healthy donors. Samples were analyzed 4, 8, 12, 24, 48 and 72 h following ex vivo irradiation with 0, 1, 3, 5 and 8 Gy of γ-rays. HPBLs were isolated and GDF-15 mRNA copy numbers were assessed with no significant differences among the three donors at each radiation dose. GDF-15 mRNA copy numbers were significantly upregulated in a dose-dependent manner at 4–48 h following γ-ray exposure (P<0.05; Fig. 5). The dose-response curves describing the relationship between radiation dose and GDF-15 mRNA expression levels at different time-points were fitted to a linear model (P<0.05; Table III). At 72 h following irradiation, GDF-15 mRNA copy numbers were significantly upregulated; however, no dose-response relationship was apparent (P>0.05).

Discussion

Radiation-responsive genes, which exhibit alterations in their expression in response to radiation exposure, may have potential as biomarkers for the evaluation of the dose of absorbed radiation. Therefore, a reproducible dose-response relationship should be established between the absorbed dose and the effect on gene expression. Previous studies have suggested that alterations in the expression of GADD45, CDKN1A and PCNA induced by ionizing radiation are time- and dose-dependent in human cell lines, HPBLs (3,6,17,29) and mice (8); however, the time- and dose-response relationship is not consistent among different genes, forms of radiation and cell types. GDF-15 is a member of the tumor growth factor-β superfamily and a downstream target of p53 (42), and serves important roles in human diseases, such as cardiomyopathy (43). Previous studies have suggested that the GDF-15 gene may participate in γ-ray-induced DNA damage, as microarray analysis identified GDF-15 as a radiation-responsive gene (19,34,36,39).

The present study demonstrated that GDF-15 mRNA and protein expression levels in AHH-1 cells were upregulated in a dose-dependent manner by 1–10 Gy of γ-rays, within 8–120 and 8–48 h following irradiation, respectively. In addition, HPBLs exposed to 1–8 Gy of γ-rays also exhibited a dose-dependent increase in GDF-15 mRNA copy number at 4–48 h post-irradiation. To the best of our knowledge, this is the first study that has demonstrated that radiation-induced alterations in GDF-15 expression in AHH-1 cells and HPBLs consistently follow a dose-response relationship within a wide range of doses and post-exposure time points. The present findings are consistent with previous studies on radiation-induced alterations in the expression of other genes, including GADD45, CDKN1A and p53-inducible gene 3 (PIG3) (3,6,18,39). Furthermore, GDF-15 baseline mRNA copy number appeared to be lower, and the time-frame for detection of radiation-induced alterations in GDF-15 expression appeared to be longer compared with PIG3, which was investigated in our previous study (39). Therefore, alterations in GDF-15 expression may have potential as a novel radiation biomarker.

Alterations in GDF-15 mRNA and protein expression levels in AHH-1 cells were observed up to 120 and 48 h following γ-ray exposure, respectively. When GDF-15 mRNA expression was analyzed, the alterations in gene expression appeared to be lower at 120 h post-irradiation compared with other time points, which may be associated with longer cell culture time and cell growth status following exposure. The differences in the duration of the effects of radiation on GDF-15 mRNA and protein may be due to the gene expression regulation, which is characterized by extensive inter-correlations between the rates of transcription, translation and protein degradation. The changes in protein levels depend on the mRNA concentration, translation efficiency and degradation of the existing protein (44). Therefore, different dose-response curves should be used according to the time period after irradiation to estimate an appropriate dose of radiation.

It has previously been demonstrated that radiation-induced chromosomal abnormalities in human lymphocytes depend on the radiation dose rate (45). Notably, in the present study, radiation-induced alterations in GDF-15 gene expression in AHH-1 cells were not modulated by the dose rate (0.5–2 Gy/min) of γ-rays, within a dose range of 3–8 Gy. Further studies, using a wider range of dose rates, are required to elucidate the effect of the dose rate on GDF-15 gene expression.

It has previously been reported that alterations in GDF-15 gene expression in human pulmonary epithelial cells were dose- and time-dependent, following exposure to α-particle radiation (19). In accordance with a previous study (34), the present results demonstrated that GDF-15 expression was altered by γ-rays and neutron radiation. These results suggested that GDF-15 expression may have potential as a radiation biomarker, and may be used in the development of radiation-induced gene expression signatures for use in biodosimetry (45). In the present study, dose-response relationships were also established between GDF-15 mRNA and protein expression levels and the dose of neutron radiation, within 48 h after exposure. A previous study reported that chromosomal aberrations induced by high-linear energy transfer (LET) radiation conformed to a linear relationship (46). Conversely, in the present study, the dose-response curve describing neutron-induced alterations in GDF-15 expression conformed to a quadratic model.

In addition, the expression changes of GDF-15 mRNA and protein in AHH-1 cells induced by 1 Gy γ-rays at 48 h post-irradiated were increased by approximately three-fold and two-fold, respectively. For the neutron exposed AHH-1 cells, the expression levels of GDF-15 mRNA and protein following exposure to 1.2 Gy at 48 h were elevated to approximately four-fold and three-fold, respectively. This difference between γ-ray- and neutron-induced changes of GDF-15 mRNA and protein expression may be caused by the different radiation qualities. High-LET radiations such as neutron are highly toxic and exert a complex effect on the genomic integrity of cells, which may lead to more severe biological damage compared with that of γ-rays (17).

The present study also investigated the expression of GDF-15 in healthy adults. GDF-15 mRNA copy number did not differ significantly between sexes, or among age groups. These results indicated the consistency in GDF-15 expression, which is an essential characteristic for a potential biomarker of radiation exposure.

Alterations in GDF-15 expression induced by ⁶⁰Co γ-rays were observed for up to 120 h following exposure in AHH-1 cells, and 72 h in HPBLs. This difference may be attributed to a difference in radiation sensitivity between the two cell types. In a previous study, AHH-1 cells were irradiated during the exponential growth phase, whereas HPBLs were in non-cycling (G₀ cell cycle phase), and DNA repair capabilities differed between proliferating and non-proliferating cells (39). The present results suggested that dose-dependent alterations in GDF-15 expression could be observed within a limited time frame in HPBLs, and this may prove to be a limitation in future applications. Further studies, using larger sample sizes, are required to elucidate the relationship between the radiation dose and alterations in GDF-15 expression.

In conclusion, the present study demonstrated that time- and dose-dependent alterations were induced in GDF-15 mRNA and protein expression levels in AHH-1 cells by γ-rays and neutron radiation. In addition, GDF-15 mRNA expression levels in HPBLs isolated from healthy adults did not differ with age or sex. Furthermore, alterations in GDF-15 mRNA expression levels in HPBLs induced by γ-rays were revealed to be dose-dependent. These results suggested that the GDF-15 gene in HPBLs may have potential as a biomarker to evaluate radiation exposure. Further studies, using a larger number of human blood samples, are required to investigate the radiation-induced alterations in GDF-15 expression.

Acknowledgements

The authors would like to thank Dr Xiang-Yan Zhou, Dr Hui-Min Lv and Dr Wei Zhang from the National Institute for Radiological Protection, Chinese Center for Disease Control and Prevention (Beijing, China) for their help and assistance with the original study protocols. The present study was supported by the National Natural Scientific Foundation of China (grant nos. 30570551, 81172593 and 81573081).

References 1

Amundson

Bittner

Meltzer

Trent

Fornace

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

Radiat Res1566576612001

10.1667/0033-7587(2001)156[0657:IOGEAA]2.0.CO;2

11604088

Amundson

Shahab

Bittner

Meltzer

Trent

Fornace

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

Radiat Res1543423462000

10.1667/0033-7587(2000)154[0342:IOPMBI]2.0.CO;2

11012342

Amundson

Fornace

Monitoring human radiation exposure by gene expression profiling: Possibilities and pitfalls

Health phys8536422003

10.1097/00004032-200307000-00009

12852469

Brengues

Paap

Bittner

Amundson

Seligmann

Korn

Lenigk

Zenhausern

Biodosimetry on small blood volume using gene expression assay

Health Phy981791852010

10.1097/01.HP.0000346706.44253.5c

Fält

Holmberg

Lambert

Wennborg

Long-term global gene expression patterns in irradiated human lymphocytes

Carcinogenesis24183718452003

10.1093/carcin/bgg134

12951355

Grace

McLeland

Blakely

Real-time quantitative RT-PCR assay of GADD45 gene expression changes as a biomarker for radiation biodosimetry

Int J Radiat Biol78101110212002

10.1080/09553000210158056

12456288

Kabacik

Mackay

Tamber

Manning

Finnon

Paillier

Ashworth

Bouffler

Badie

Gene expression following ionising radiation: Identification of biomarkers for dose estimation and prediction of individual response

Int J Radiat Biol871151292011

10.3109/09553002.2010.519424

21067298

Wang

Chen

Shen

Min

Radiation dose effect of DNA repair-related gene expression in mouse white blood cells

Med Sci Monit17BR290BR2972011

10.12659/MSM.881976

21959603

3539470

Moulder

Post-irradiation approaches to treatment of radiation injuries in the context of radiological terrorism and radiation accidents: A review

Int J Radiat Biol803102004

10.1080/09553000310001642920

14761845

Roy

Gruel

Vaurijoux

Cell response to ionising radiation analysed by gene expression patterns

Ann Ist Super Sanita452722772009

19861732

Amundson

Fornace

Gene expression profiles for monitoring radiation exposure

Radiat Prot Dosimetry9711162001

10.1093/oxfordjournals.rpd.a006632

11763352

Fachin

Mello

Sandrin-Garcia

Junta

Donadi

Passos

Sakamoto-Hojo

Gene expression profiles in human lymphocytes irradiated in vitro with low doses of gamma rays

Radiat Res1686506652007

10.1667/RR0487.1

18088177

Fachin

Mello

Sandrin-Garcia

Junta

Ghilardi-Netto

Donadi

Passos

Sakamoto-Hojo

Gene expression profiles in radiation workers occupationally exposed to ionizing radiation

J Radiat Res5061712009

10.1269/jrr.08034

19218781

Filiano

Fathallah-Shaykh

Fiveash

Gage

Cantor

Kharbanda

Johnson

Gene expression analysis in radiotherapy patients and C57BL/6 mice as a measure of exposure to ionizing radiation

Radiat Res17649612011

10.1667/RR2419.1

21361780

Gruel

Lucchesi

Pawlik

Frouin

Alibert

Kortulewski

Zarour

Jacquelin

Gidrol

Tronik-Le Roux

Novel microarray-based method for estimating exposure to ionizing radiation

Radiat Res1667467562006

10.1667/RR0260.1

17067202

Kang

Park

Song

Jeoung

Cho

Kim

Bae

Lee

Possible biomarkers for ionizing radiation exposure in human peripheral blood lymphocytes

Radiat Res1593123192003

10.1667/0033-7587(2003)159[0312:PBFIRE]2.0.CO;2

12600233

Turtoi

Brown

Schläger

Schneeweiss

Gene expression profile of human lymphocytes exposed to (211)At alpha particles

Radiat Res1741251362010

10.1667/RR1659.1

20681779

Akerman

Rosenzweig

Domon

Tsai

Bishop

McGarrity

MacGregor

Sistare

Chen

Morris

Alterations in gene expression profiles and the DNA-damage response in ionizing radiation-exposed TK6 cells

Environ Mol Mutagen451882052005

10.1002/em.20091

15657912

Chauhan

Howland

Gene expression responses in human lung fibroblasts exposed to alpha particle radiation

Toxicol In Vitro28122212292014

10.1016/j.tiv.2014.06.001

24945610

Grace

McLeland

Gagliardi

Smith

Jackson

IIIBlakely

Development and assessment of a quantitative reverse transcription-PCR assay for simultaneous measurement of four amplicons

Clin Chem49146714752003

10.1373/49.9.1467

12928227

Jen

Cheung

Transcriptional response of lymphoblastoid cells to ionizing radiation

Genome Res13209221002003

10.1101/gr.1240103

12915489

403696

Paul

Amundson

Gene expression signatures of radiation exposure in peripheral white blood cells of smokers and non-smokers

Int J Radiat Biol877918012011

10.3109/09553002.2011.568574

21801107

3580190

Pogosova-Agadjanyan

Fan

Georges

Schwartz

Kepler

Lee

Suchanek

Cronk

Brumbaugh

Engel

Identification of radiation-induced expression changes in nonimmortalized human T cells

Radiat Res1751721842011

10.1667/RR1977.1

21268710

Dressman

Muramoto

Chao

Meadows

Marshall

Ginsburg

Nevins

Chute

Gene expression signatures that predict radiation exposure in mice and humans

PLoS Med4e1062007

10.1371/journal.pmed.0040106

17407386

1845155

Meadows

Dressman

Muramoto

Himburg

Salter

Wei

Ginsburg

Chao

Nevins

Chute

Gene expression signatures of radiation response are specific, durable and accurate in mice and humans

PLoS One3e19122008

10.1371/journal.pone.0001912

18382685

2271127

Amundson

Grace

McLeland

Epperly

Yeager

Zhan

Greenberger

Fornace

Human in vivo radiation-induced biomarkers: Gene expression changes in radiotherapy patients

Cancer Res64636863712004

10.1158/0008-5472.CAN-04-1883

15374940

Turtoi

Brown

Oskamp

Schneeweiss

Early gene expression in human lymphocytes after gamma-irradiation-a genetic pattern with potential for biodosimetry

Int J Radiat Biol843753872008

10.1080/09553000802029886

18464067

Turtoi

Schneeweiss

Effect of (211)At alpha-particle irradiation on expression of selected radiation responsive genes in human lymphocytes

Int J Radiat Biol854034122009

10.1080/09553000902838541

19382019

Paul

Amundson

Development of gene expression signatures for practical radiation biodosimetry

Int J Radiat Oncol Biol Phys71123612442008

10.1016/j.ijrobp.2008.03.043

18572087

2478639

Paul

Barker

Turner

McLane

Wolden

Amundson

Prediction of in vivo radiation dose status in radiotherapy patients using ex vivo and in vivo gene expression signatures

Radiat Res1752572652011

10.1667/RR2420.1

21388269

3134561

Badie

Kabacik

Balagurunathan

Bernard

Brengues

Faggioni

Greither

Lista

Peinnequin

Poyot

Laboratory intercomparison of gene expression assays

Radiat Res1801381482013

10.1667/RR3236.1

23886340

4117252

Boldt

Knops

Kriehuber

Wolkenhauer

A frequency-based gene selection method to identify robust biomarkers for radiation dose prediction

Int J Radiat Biol882672762012

10.3109/09553002.2012.638358

22233095

Brzóska

Kruszewski

Toward the development of transcriptional biodosimetry for the identification of irradiated individuals and assessment of absorbed radiation dose

Radiat Environ Biophys543533632015

10.1007/s00411-015-0603-8

25972268

4510913

Joiner

Thomas

Grever

Smolinski

Divine

Konski

Auner

Tucker

Developing point of care and high-throughput biological assays for determining absorbed radiation dose

Radiother Oncol1012332362011

10.1016/j.radonc.2011.05.068

21724286

Tucker

Joiner

Thomas

Grever

Bakhmutsky

Chinkhota

Smolinski

Divine

Auner

Accurate gene expression-based biodosimetry using a minimal set of human gene transcripts

Int J Radiat Oncol Biol Phys889339392014

10.1016/j.ijrobp.2013.11.248

24444760

Kis

Szatmári

Keszei

Farkas

Ésik

Lumniczky

Falus

Sáfrány

Microarray analysis of radiation response genes in primary human fibroblasts

Int J Radiat Oncol Biol Phys66150615142006

10.1016/j.ijrobp.2006.08.004

17069989

Rødningen

Overgaard

Alsner

Hastie

Børresen-Dale

Microarray analysis of the transcriptional response to single or multiple doses of ionizing radiation in human subcutaneous fibroblasts

Radiother Oncol772312402005

10.1016/j.radonc.2005.09.020

16297999

Warters

Packard

Kramer

Gaffney

Moos

Differential gene expression in primary human skin keratinocytes and fibroblasts in response to ionizing radiation

Radiat Res17282952009

10.1667/RR1677.1

19580510

2734195

Liu

Zhang

Feng

Wang

Chen

Gao

Dose-effect of ionizing radiation-induced PIG3 gene expression alteration in human lymphoblastoid AHH-1 cells and human peripheral blood lymphocytes

Int J Radiat Biol9171802015

10.3109/09553002.2014.938374

24991881

Crespi

Thilly

Assay for gene mutation in a human lymphoblast line, AHH-1, competent for xenobiotic metabolism

Mutat Res1282212301984

10.1016/0027-5107(84)90110-6

6088974

Livak

Schmittgen

Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method

Methods254024082001

10.1006/meth.2001.1262

11846609

Osada

Park

Liu

Trink

Sidransky

A p53-type response element in the GDF15 promoter confers high specificity for p53 activation

Biochem Biophys Res Commun3549139182007

10.1016/j.bbrc.2007.01.089

17276395

1975780

Montoro-García

Hernández-Romero

Jover

García-Honrubia

Vilchez

Casas

Martínez

Climent

Caballero

Valdés

Marín

Growth differentiation factor-15, a novel biomarker related with disease severity in patients with hypertrophic cardiomyopathy

Eur J Intern Med231691742012

10.1016/j.ejim.2011.08.022

22284249

de Sousa Abreu

Penalva

Marcotte

Vogel

Global signatures of protein and mRNA expression levels

Mol Biosyst5151215262009

20023718

4089977

Purrott

Reeder

The effect of changes in dose rate on the yield of chromosome aberrations in human lymphocytes exposed to gamma radiation

Mutat Res354374441976

10.1016/0027-5107(76)90206-2

934166

Upton

Environmental standards for ionizing radiation: Theoretical basis for dose-response curves

Environ Health Perspect5231391983

10.1289/ehp.835231

6653536

1569362

Figure 1.

Exposure to γ-rays altered GDF-15 mRNA expression levels in AHH-1 cells. (A) AHH-1 cells were irradiated with 0–10 Gy ⁶⁰Co γ-rays at a dose rate of 1 Gy/min. Cells were collected for analysis 8–168 h post-irradiation. Relative GDF-15 mRNA expression levels were assessed using RT-qPCR. GDF-15 mRNA expression levels were significantly upregulated 8–120 h following irradiation compared with the control group (0 Gy). (B) AHH-1 cells were irradiated with 3, 5 or 8 Gy of ⁶⁰Co γ-rays at a dose-rate of 0.5, 1 and 2 Gy/min. Cells were collected for analysis 48 h post-irradiation. Relative GDF-15 mRNA expression levels were assessed using RT-qPCR. Alterations in GDF-15 expression depended on radiation dose, regardless of the dose rate. Data are expressed as the mean ± standard error of the mean of three independent experiments. *P<0.05 vs. 0 Gy group. GDF, growth/differentiation factor; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

Figure 2.

Exposure to γ-rays altered GDF-15 protein expression levels in AHH-1 cells. AHH-1 cells were irradiated with 0–10 Gy of ⁶⁰Co γ-rays at a dose rate of 1 Gy/min. Cells were collected for western blot analysis 8–72 h post-irradiation. (A) Representative blots demonstrating alterations in GDF-15 protein expression levels. (B) Semi-quantification of western blot analysis results, normalized to β-actin. Relative GDF-15 protein expression levels were significantly upregulated following irradiation compared with the control group (0 Gy). Data are expressed as the mean ± standard error of the mean of three independent experiments. *P<0.05 vs. 0 Gy group. GDF, growth/differentiation factor.

Figure 3.

Exposure to neutron radiation altered GDF-15 mRNA and protein expression levels in AHH-1 cells. AHH-1 cells were irradiated with 0–1.6 Gy of ²⁵²Cf neutron radiation at a dose rate of 0.073 Gy/min. (A) GDF-15 mRNA expression levels were significantly upregulated 48 h following irradiation compared with the control group (0 Gy). (B) Representative blot demonstrating alterations in GDF-15 protein expression levels. GDF-15 protein expression levels were visibly increased 48 h following irradiation compared with the control group (0 Gy). GDF-15 protein expression levels were normalized to β-actin. Data are expressed as the mean ± standard error of the mean of three independent experiments. *P<0.05 vs. 0 Gy group. GDF, growth/differentiation factor.

Figure 4.

GDF-15 mRNA levels in peripheral blood samples of 73 healthy adult donors. (A) GDF-15 mRNA expression levels (mRNA copy numbers per µl) in males were not significantly different compared with in females. (B) No significant differences in GDF-15 mRNA expression levels were revealed among different age groups; however, an increasing trend was observed within the 20–50 age range. Data are expressed as the mean ± standard error of the mean of three independent experiments. GDF, growth/differentiation factor.

Figure 5.

Exposure to γ-rays altered GDF-15 mRNA expression levels in HPBLs from three healthy adult donors. HPBLs were irradiated with 1–8 Gy of ⁶⁰Co γ-rays at a dose rate of 1 Gy/min. Cells were collected for analysis 4–72 h post-irradiation. GDF-15 mRNA expression levels (mRNA copy number per µl) were significantly upregulated in a dose-dependent manner 4–48 h following irradiation compared with the control group (0 Gy). Data are expressed as the mean ± standard error of the mean of three independent experiments. *P<0.05 vs. 0 Gy group. GDF, growth/differentiation factor; HPBL, human peripheral blood lymphocyte.

Table I.

Equations describing the dose-response relationship between the absorbed dose (0–10 Gy) of ⁶⁰Co γ-rays (dose rate, 1 Gy/min) and GDF-15 mRNA expression levels in AHH-1 cells.

Time post-irradiation (h)	Dose-response equation	R²
8	y=0.62x+2.75	0.728
12	y=1.31x+0.54	0.963
24	y=0.63x+0.92	0.954
48	y=0.96x+1.77	0.978
72	y=1.09x+1.59	0.970
120	y=0.48x+0.65	0.914

y, relative GDF-15 mRNA expression level (fold change); x, absorbed dose (Gy). GDF, growth/differentiation factor.

Table II.

Equations describing the dose-response relationship between the absorbed dose (0–10 Gy) of ⁶⁰Co γ-rays (dose rate, 1 Gy/min) and GDF-15 protein expression levels in AHH-1 cells.

Time post-irradiation (h)	Dose-response equation	R²
8	–	–
12	y=1.81x-0.34	0.959
24	y=0.43x+0.53	0.962
48	y=0.33x+0.85	0.926
72	–	–
120	–	–

y, relative GDF-15 protein expression level (fold change); x, absorbed dose (Gy). GDF, growth/differentiation factor.

Table III.

Equations describing the dose-response relationship between the absorbed dose (0–8 Gy) of ⁶⁰Co γ-rays (dose rate, 1 Gy/min) and GDF-15 mRNA expression levels in human peripheral blood lymphocytes.

Time post-irradiation (h)	Dose-response equation	R²
4	y=7,510x+23,380	0.797
8	y=27,480x+47,200	0.970
12	y=14,786x+75,212	0.851
24	y=14,210x+79,197	0.955
48	y=28,707x+72,207	0.899

y, number of GDF-15 mRNA copies/µl; x, absorbed dose (Gy); GDF, growth/differentiation factor.