- 1Department of Radiotherapy, Veneto Institute of Oncology IOV–IRCCS, Padua, Italy
- 2Department of Biology, University of Padua, Padua, Italy
- 3National Laboratories of Legnaro, Italian Institute of Nuclear Physics (LNL-INFN), Padua, Italy
- 4Department of Molecular Medicine, University of Padua, Padua, Italy
Introduction: Adverse effects of radiotherapy (RT) significantly affect patient's quality of life (QOL). The possibility to identify patient-related factors that are associated with individual radiosensitivity would optimize adjuvant RT treatment, limiting the severity of normal tissue reactions, and improving patient's QOL. In this study, we analyzed the relationships between genetic features and toxicity grading manifested by RT patients looking for possible biomarkers of individual radiosensitivity.
Methods: Early radiation toxicity was evaluated on 143 oncological patients according to the Common Terminology Criteria for Adverse Events (CTCAE). An individual radiosensitivity (IRS) index defining four classes of radiosensitivity (highly radiosensitive, radiosensitive, normal, and radioresistant) was determined by a G2-chromosomal assay on ex vivo irradiated, patient-derived blood samples. The expression level of 15 radioresponsive genes has been measured by quantitative real-time PCR at 24 h after the first RT fraction, in blood samples of a subset of 57 patients, representing the four IRS classes.
Results: By applying univariate and multivariate statistical analyses, we found that fatigue was significantly associated with IRS index. Interestingly, associations were detected between clinical radiation toxicity and gene expression (ATM, CDKN1A, FDXR, SESN1, XPC, ZMAT3, and BCL2/BAX ratio) and between IRS index and gene expression (BBC3, FDXR, GADD45A, and BCL2/BAX).
Conclusions: In this prospective cohort study we found that associations exist between normal tissue reactions and genetic features in RT-treated patients. Overall, our findings can contribute to the identification of biological markers to predict RT toxicity in normal tissues.
Introduction
The development of radiation-induced complications following radiotherapy (RT) has a significant impact on treatment outcome and patient's quality of life (QOL). In the last decades the therapeutic ratio has improved due to advancements in RT technologies and use of radioprotectors, mitigators, and radiosensitizers (1). Nevertheless, radiation toxicity of normal tissues surrounding the tumor is a serious problem for ~5–10% of patients, who are affected by high intrinsic radiosensitivity (2–4). Evidence of radiosensitivity in vivo is given by burns and radiodermitis in the irradiated body parts, together with bystander effect in neighboring area (5). Several factors, including cellular composition, differentiation, cell renewal capacity, as well as cellular radiosensitivity, determine the severity of radiation toxicity (6). Patient-related factors are deeply linked to the risk of manifesting radiation toxicity, and reliable biological markers are still not available to predict the onset of severe side-effects after RT. Human response to ionizing radiation (IR) is individual and variable, being influenced by age, smoking, diabetes, collagen vascular disease and genotype (7). Moreover, multiple genetic pathways such as DNA damage repair, oxidative stress, radiation fibrogenesis and endothelial cell damage are implicated in adverse tissue reactions following radiotherapy (8). However, the molecular basis of individual radiosensitivity remains poorly understood, and the relationship between different indicators of radiation sensitivity is elusive.
RT causes cancer cell death mainly by IR-induced DNA Double Strand Breaks (DSBs). Formation of DSBs, the most severe damage for genome integrity, triggers a cascade of cellular events, collectively termed DNA-damage response (DDR), which involves sensing the damage, signal transduction to the effectors of DNA repair, cell cycle arrest and apoptosis induction (9–11). Radiation-induced DSBs are efficiently repaired to ensure the maintenance of genome integrity but when DNA repair is hampered, unrepaired DSBs can originate chromosome aberrations (12). Following irradiation, unrepaired DSBs can be quantified in metaphase spreads by the yield of chromatid breaks formed at G2-phase, which is inversely related to the efficiency of the G2-phase checkpoint (13). Thus, the individual level of radiosensitivity can be assessed in ex vivo irradiated human peripheral blood lymphocytes (PBLs) by applying a “G2-chromosomal assay” (14–17).
Increasing evidence supports the existence of individual response to IR-induced DNA damage, which can be related to mutations in key genes of DDR pathway or to individual capacity to modulate the expression of DDR genes after IR-exposure. In this regard, the expression of genes involved in DDR pathway may be variable between individuals and can impact on own radiation response.
Several studies attempted to find biomarkers able to predict the onset of radiation toxicity in normal tissues after RT. Individual radiosensitivity evaluated by using in vitro irradiated patient-derived blood lymphocytes has been found to correlate with normal tissue reactions (13, 18, 19), and single nucleotide polymorphisms (SNPs) have been associated with acute and late radiation-induced normal tissue injury in RT patients (19–22). Data concerning the association between gene expression changes and normal tissue radiation toxicity refer to in vitro irradiation studies (4, 23–25) or to single gene analysis in oncological patients treated with RT (26). To date, a relationship between radio-induced normal tissue adverse effects, in vitro chromosomal radiosensitivity and in vivo expression of a set of radioresponsive genes is not available in the same cohort of RT patients. Since future clinical protocols aim at ameliorating patient's QOL it is demanding to identify patient-related factors that are associated with individual radiosensitivity before patients undergo RT (27, 28).
In this explorative study, the clinical features of early radiation toxicity have been associated with an Individual Radiosensitivity (IRS) index, defining four classes of radiosensitivity (highly radiosensitive, radiosensitive, normal and radioresistant) based on a G2-chromosomal assay on patient-derived PBLs irradiated in vitro (15). The expression level of 15 selected radioresponsive genes belonging to DDR pathway has been measured in blood samples from a subgroup of patients, representing the different IRS classes, 24 h after the first RT fraction, as an additional variable of intrinsic radiosensitivity. Data of clinical and genetic features have been statistically analyzed to find possible genetic factors associated with individual radiation sensitivity.
Materials and Methods
Outline of the Study
In this prospective study, breast cancer (BC) and head and neck squamous cell carcinoma (HNSCC) patients were enrolled as representative of patients experiencing normal tissue reactions after RT. Data of toxicity grading in normal tissues, in vitro chromosomal radiosensitivity and in vivo RT-induced gene expression changes, have been integrated to identify possible biomarkers of radiosensitivity in patients undergoing RT (Figure 1). Overall, 143 oncological patients were enrolled: 124 (all females) affected by BC, and 19 (6 females and 13 males) affected by HNSCC.
Figure 1. Flow chart showing the experimental phases for the identification of biomarkers of individual radiosensitivity in patients undergoing radiotherapy (RT).
Patients
Patients with BC or HNSCC histological diagnosis undergoing RT were enrolled from 2015 to 2017 at the Department of Radiotherapy, Veneto Institute of Oncology IOV–IRCCS, Padua, Italy (IOV) upon evaluation and approval of the IOV-IRCCS Ethic Committee (CE IOV 2015/18; CE IOV 2016/04). Privacy rights of human subjects were observed; all the procedures were in accordance with relevant guidelines and regulations. All subjects gave written informed consent in accordance with the Declaration of Helsinki.
Patients were enrolled applying the following exclusion criteria: patients suffering from congenital syndromes predisposing to radiosensitivity (such as Ataxia-Telangiectasia, Bloom syndrome, Down's syndrome, Gorlin syndrome, Klinefelter syndrome, Retinoblastoma, Wilm's tumor, Xeroderma pigmentosum, Rothmund-Thomson syndrome, Li-Fraumeni syndrome, Dyskeratosis congenita, Familial dysplastic nevus syndrome, Common variable immune deficiency, Nijmegen Breakage Syndrome, Fanconi Anemia, albinism), previous RT and/or chemotherapy treatment or ongoing chemotherapy treatment, previous anticancer drug employment, significant comorbidities, diabetic patients affected by breast cancer, age ≤ 18 years.
BC patients received by 3DCRT (42.40–50 Gy/16–25 fractions) plus a boost dose of 10 Gy in 5 fractions to the tumor bed. HNSCC patients received up to 70 Gy, daily fraction 1.8–2.12 Gy/day, for 5 days/week on primitive tumor by VMAT or IMRT.
Adverse tissue reactions (dermatitis radiation, pain, pruritus, fatigue) have been recorded at the completion of RT treatment (t1) and 1 month later (t2), using the Common Terminology Criteria for Adverse Events (CTCAE) (version 4.03, http://ctep.cancer.gov/protocolDevelopment/adverse_effects.htm). Adverse effects were classified as: grade 0 (G0, no adverse effects), grade 1 (G1, mild), grade 2 (G2, moderate), grade 3 (G3, severe). At the Department of Radiotherapy of the IOV-IRCCS management of acute toxicity followed a standardized procedure. All patients were clinically evaluated before starting RT and no significant side-effects were complained by patients.
Chromosome-Based Radiosensitivity Assay
The G2-chromosomal assay was performed following a standardized protocol (15). Briefly, whole blood cultures were incubated for 72 h at 37°C, 5% CO2 before being irradiated with 1 Gy of gamma rays in a Gamma Beam A15 60Co panoramic source at the National Laboratories of Legnaro (I.N.F.N., Padua, Italy; dose rate: 0.5 Gy/min).
Immediately after irradiation, each culture was split in two and one was treated with 4 mM caffeine. After 20 min at 37°C, both cultures were incubated with Colcemid at concentration of 0.1 μg/mL for 60 min, then chromosome spreads were prepared according to standard cytogenetic procedures. With few exceptions, chromatid aberration yields were obtained by scoring for chromatid breaks and gaps 50 metaphases per culture, under a Zeiss AxioImager Z2 microscope coupled with MSearch-AutoCapt software (Metasystems, Altlussheim Germany). Following calculation of the in vitro individual radiosensitivity index (IRS = [1-(G2caf-G2)/G2caf] × 100%, simplified as IRS = (G2/G2caf) × 100%) patients were classified as: highly radiosensitive, HRS (IRS > 70), radiosensitive, RS (50 < IRS ≤ 70), normal, N (30 ≤ IRS ≤ 50), and radioresistant, RR (IRS < 30) (15, 29).
Gene Expression Analysis
Fifty-seven over 143 patients were randomly selected within the four IRS classes (HRS, RS, N, RR) in order to have comparable numbers of patients in each group. This sample size guarantees a high statistical power (power = 0.83) in identifying as significant (alpha < 0.1) genes with an effect equal to 1.1 among groups using either an ANOVA test or a Wilcoxon test. Two whole blood samples were collected from each patient: one immediately before the first fractionated RT dose and the second 24 h later. Samples were collected into PAXgene® Blood RNA tubes (PreAnalytiX GmbH, Qiagen, Venlo, The Netherlands) for immediate stabilization of intracellular RNA, and stored at −80°C. Total RNA was purified by using PAXgene® Blood RNA Kit 6 (PreAnalytiX GmbH, Qiagen, Venlo, The Netherlands) and quantified using the ND-1000 spectrophotometer (Nanodrop, Wilmington, DE, USA).
For mRNA detection, retrotranscription and quantitative real time-PCR (qRT-PCR) reactions were performed according to our established protocol (30, 31). The gene-specific primers for ATM, BAX, BBC3, BCL2, CCNG1, cMYC, DDB2, FDXR, GADD45A, MDM2, CDKN1A, PCNA, SESN1, XPC, and ZMAT3 genes and for GADPH as reference, can be found in Supplementary Table 1. Real-time PCR was performed using an Applied Biosystems 7500 Fast Real-Time PCR System according to the following amplification protocol: 95°C for 10 min, 95°C for 15 sec, 60°C for 60 s (40 cycles). qRT-PCR reactions were always performed in triplicates. The relative expression levels of mRNAs between irradiated (2 Gy) and non-irradiated (0 Gy) blood samples of the same patients were calculated using the comparative delta CT (threshold cycle number) method (2−ΔΔCT) implemented in the 7500 Real Time System software (32).
Statistical Analysis
ANOVA and ANOVA post-hoc with Bonferroni correction was used to assess the gene expression mean differences among groups of patients defined using clinical and IRS annotations. In case of group size lower than 10 patients Kruskall-Wallis/Wilcoxon test was used. Multivariate regression analyses were used to test the association of clinical annotations (explanatory variables) with IRS value (dependent variable). Unsupervised cluster analysis of gene expression data was performed using hierarchical cluster analysis with Euclidean distance and complete linkage. All the analyses were performed using the R programming language (version 3.4), and the Bioconductor software suite (version 3.6).
Results
Radiation-Induced Toxicity in RT-Treated Patients
All patients were evaluated for the onset of radiation toxicity at the completion of RT treatment (t1) and 1 month later (t2). The overall distribution of subjects suffering from dermatitis radiation, pain, pruritus and fatigue is reported in Figure 2; Supplementary Table 2. Moderate (G2) dermatitis radiation was recorded at t1 in 13.7 and 21% of BC and HNSCC patients, respectively, whereas severe (G3) dermatitis radiation was observed in 5.6 and 10.5% of BC and HNSCC, respectively. At t2, G2 dermatitis radiation was observed in 8% of BC patients and in 33% of HNSCC patients; G3 dermatitis radiation was manifested by 0.8% of BC patients. Pain of G2 grade was present in 5.6% of BC and in 15.8% of HNSCC patients at t1; 1 month later (t2), 2.4% of BC and 18% of HNSCC patients manifested G2 pain. At t1, G3 pain was rarely recorded in BC patients but affected 21% of the HNSCC patients; at t2 none of the BC and HNSCC patients suffered from pain of G3 grade. Pruritus of G2 and G3 grade was recorded at t1 in 9.7 and 3.2% of BC patients, respectively, while at t2, G2 pruritus was present in 8.9% of BC patients. Concerning HNSCC patients, 5.2 and 5.5% of them manifested G3 pruritus at t1 and t2, respectively. In summary, HNSCC patients manifested higher degrees of dermatitis radiation and pain at both t1 and t2, whereas pruritus appeared to be more pronounced in BC patients. Fatigue (G1 grade) was present in 50.8% (BC) and 57.8% (HNSCC) of patients, and at t2 in 33.0 and 44.5%, respectively (Figure 2).
Figure 2. Clinical variables of radiation toxicity in RT patients. RT-induced toxicity grading in breast cancer, BC (A) and in head and neck squamous cell carcinoma HNSSC (B) patients at completion of RT treatment (t1) and 1 month later (t2).
G2-Chromosomal Radiosensitivity in RT Patients
IRS values were determined in the complete patient cohort (143 subjects). In blood cultures exposed in vitro to 1 Gy [according to the standardized protocol developed by Pantelias and coworkers (15, 29)], the average yield of G2 chromatid breaks was 2.5 (with standard deviation SD = 0.084 and coefficient of variation CV = 3%); the average IRS was 38.1 (SD = 10.62; CV = 27.8%). Based on the observed distribution of individual IRS values, the four classes of individual radiosensitivity should be, respectively: RR < 27.48 (mean – SD); 27.48 ≤ N ≤ 48.72 (mean ± SD); 48.72 < RS ≤ 69.96 (mean + SD); HRS > 69.96 (mean + 3 × SD) (Figure 3A). As these values are in strict agreement with those proposed earlier (15), for further statistical analyses we used the published thresholds (see Materials and Methods). One patient resulted highly radiosensitive (HRS, 1%), 16 patients were classified as radiosensitive (RS, 11%), 95 patients as normal (N, 66%) and 31 as radioresistant (RR, 22%) (Figure 3B).
Figure 3. (A) Data distribution for IRS values. Vertical lines indicate the observed thresholds for the 4 classes of individual radiosensitivity, calculated as: RR, mean - SD (IRS < 27.48); N, mean ± SD (27.48 ≤ IRS ≤ 48.72); RS, mean + SD (48.72 < IRS ≤ 69.96); HRS, mean + 3 SD (IRS > 69.96). (B) Patient distribution according to the four IRS classes.
Gene Expression in Blood Samples of RT Patients
Fifteen radioresponsive genes belonging to DDR pathway (Table 1) were analyzed by qRT-PCR in blood samples from 57 patients, randomly selected from the whole cohort, in order to have comparable numbers of patients within the four IRS classes (HRS, RS, N, and RR). A summary of the clinical data of this group of patients is available in Supplementary Table 3.
Transcription of most genes was significantly induced after the first RT fraction (Figure 4A).
Figure 4. (A) Gene expression analysis by qRT-PCR in blood samples from RT patients. The relative mRNA quantification was performed by comparing irradiated vs. non-irradiated blood samples derived from the same patient. Values are mean ± SE and expressed in fold-change. The value “1” of non-irradiated control (light gray bars) is arbitrarily given when no change is observed (***p < 0.001; **p < 0.01; *p < 0.05). (B) Heatmap and unsupervised cluster analysis on the expression profiles of DDR genes analyzed in 57 RT patients. The key color bar indicates standardized gene expression levels (low levels are in red, high levels are in yellow). The annotation bars (upper part of the heatmap) indicate the four classes of IRS index and tumor types.
The unsupervised cluster analysis of gene expression profiles reported as a heatmap in Figure 4B did not reveal differences between the two types of cancer, although the response of DDR genes was variable across patients. Indeed, a group of patients is characterized by high expression values of MDM2, SESN1, BCL2, ATM, and ZMAT3, however this group did not show significant enrichment for any of the available clinical annotations. Finally, the heatmap did not show any association between IRS index and gene expression changes (Figure 4B).
Identification of Biomarkers of Radiosensitivity
Univariate and multivariate statistical analyses were performed looking for association between (i) clinical variables of radiation toxicity and IRS index; (ii) clinical variables of radiation toxicity and expression level of DDR genes; (iii) IRS index and expression level of DDR genes.
Relationships between clinical variables of radiation toxicity and IRS classes are shown in Figure 5. At the completion of RT treatment (t1), patients experiencing adverse effects were distributed within the four IRS classes (HRS, RS, N, and RR), without any significant relationship and without differences between tumor types. Instead, for fatigue the IRS mean values significantly differed between patients with and without such adverse effect. Specifically, the estimated IRS mean values were, respectively, 40.44 and 36.17, with a decrease of 4.27 in patients with G1 fatigue at t1 (p = 0.015, t-test). The significance was confirmed by multivariate linear regression model (adjusted for age and disease type) (Table 2). Neither the multivariate nor the univariate analyses showed significant association between IRS index and clinical variables of radiation toxicity at t2.
Figure 5. Association analysis between clinical variables of radiation toxicity and IRS index. Number of patients in the four IRS classes (HRS, RS, N, RR) per grades of radiation-induced adverse effects (G0, G1, G2, G3), at the completion of RT treatment (t1) and 1 month later (t2). NA, data Not Available.
Table 2. Multivariate linear regression model with IRS values and all clinical variables as covariates in 143 patients.
Significant and moderately significant associations between gene expression and clinical radiation toxicity are shown in Table 3. Dermatitis radiation at t1 was associated with a 1.88-fold change of FDXR expression in patients experiencing G3 toxicity vs. a 1.44-fold change in G2 patients. The presence of pain at t1 was associated with a decrease of SESN1 expression (0.92- vs. 1.36-fold change in the comparison presence-absence, and 0.97- vs. 1.36-fold change when comparing more specifically G1 vs. G0). Symptoms of pruritus resulted associated at t1 with a 0.79-fold change of XPC and with a 1.01-fold change of ZMAT3; at t2 pruritus was associated with a 0.87-fold change of ATM, and a lower BCL2/BAX ratio (respectively, 0.76- vs. 1.15-fold change). G1 fatigue resulted associated with a 1.19-fold change of CDKN1A (p21) at the second clinical evaluation.
By univariate analyses (Table 4) we found a moderate significant association between the RS class and BBC3 and FDXR expression (adjusted p = 0.069) and between RR class and GADD45A expression (adjusted p = 0.096). The BCL2/BAX ratio was also associated with the RS class (adjusted p = 0.017).
Discussion
Despite the advancements in understanding and preventing RT effects on normal tissue, injuries deriving from radiation therapy cannot be avoided (33–35). Inter-individual differences in radiosensitivity are due to different endogenous and exogenous factors (e.g., DNA repair capacity, age, diet, and life-style) as well as to the experimental endpoint (clinical radiation toxicity, chromosome aberrations, etc.) (15, 27, 36, 37). Assessing the intrinsic component of radiosensitivity before RT could predict toxicity risk and improve the QOL (27, 28). To this purpose, from a cohort of oncological patients we collected data concerning radiation toxicity in normal tissues, in vitro G2-chromosomal radiosensitivity and in vivo expression level of 15 selected radioresponsive genes of DDR pathway, to find possible associations between genetic features and clinical radiosensitivity (Figure 1). By univariate and multivariate statistical models we have looked at significant associations between clinical and molecular data, controlling for potential confounders and for the multiplicity of the tests. Remarkably, no statistical differences have found between tumor types, allowing us to discuss our data as a whole.
To the best of our knowledge, this is the first study assessing the relationship between the three experimental endpoints in a cohort of RT-treated oncological patients. It is noteworthy that all patients have been enrolled in the same Radiotherapy Unit (IOV-IRCSS). Specifically, in this explorative study we considered breast and head and neck cancer patients as representative of patients experiencing radiation toxicity. Indeed, symptoms of grade 2 acute skin toxicity are observed in 15–24% of breast cancer patients at the completion of RT treatment (35, 38) whereas dermatitis radiation continues to be one of the most common side effects of RT in head and neck cancers (33, 39). In the present study, we considered dermatitis radiation, pain, pruritus and fatigue that are adverse effects commonly manifested in BC and HNSCC patients after RT, while tumor-specific adverse effects were excluded. Overall, HNSCC patients manifested higher degrees of dermatitis radiation and pain both at t1 and t2, whereas pruritus was more pronounced in BC patients at both t1 and t2 (Figure 2). For patients experiencing the highest level of dermatitis radiation we verified ex post the lack of relation with the phototype (Fitzpatrick scale). Fatigue induced by RT is a common symptom experienced by patients that deeply affects their QOL (40). In our cohort, all patients manifesting fatigue were evaluated as G1 grade, with overlapping proportions irrespective of cancer type: at t1 50.8% (BC) and 57.8% (HNSCC), and at t2 33.0 and 44.5%, respectively (Figure 2).
Previous studies showed that clinical radiation toxicity is related to G2-chromosomal radiosensitivity of in vitro irradiated lymphocytes (13, 19, 41). Here, we followed the standardized G2-assay developed by Pantelias and Terzoudi (15) in which the G2-checkpoint efficiency is abrogated by caffeine (inhibitor of ATM kinase) to maximize the radio-induced chromosomal damage, i.e., simulating the condition of high radiosensitivity of AT (Ataxia Telangiectasia) patients. This leads to accurate estimations of the individual radiosensitivity (the IRS index) by calculating the percentage ratio between the yields of radio-induced chromatid breaks in presence or absence of the functional G2-checkpoint (15, 29). IRS values obtained in the present study were distributed in strict agreement with previously published data (15, 29), confirming the reproducibility of the standardized G2-assay for assessing individual radiosensitivity in vitro. Based on multivariate analyses, fatigue emerged as the only adverse effect strictly associated with IRS index. Interestingly, in patients displaying G1 vs. G0 fatigue but having same values of other predictors, the average IRS index differed for a value of 4.27. No other clinical reactions were found associated with IRS values in this statistical analysis (Table 2).
Clinical radiosensitivity can be associated with individual factors, such as abnormal transcriptional responses to DNA damage and with defects in DNA repair (42–44). In this regard, previous studies of gene expression profiling, carried out in patient-derived PBLs irradiated in vitro, succeeded to some extent in discriminating groups of patients with and without severe late radiotherapy toxicity (23). An association was also observed between early skin reaction and the transcriptional response of lymphoblastoid cells derived from patients with acute radiation toxicity (4). The candidate genes here analyzed, belonging to the DNA Damage Response (DDR) pathway, were chosen on the basis of our previous data showing significant changes in their expression level in human PBLs at 24 h after irradiation with 2 Gy of γ-rays (30, 31). Moreover, GADD45A, CDKN1A, DDB2 and XPC, together with FDXR gene are well-known radio-responsive genes (4, 26, 30, 31, 45, 46). Gene expression analyses have been carried out in 57 patients randomly selected within the cohort of 143 patients, as representative of the four IRS classes (HRS, RS, N, and RR). This sample size is adequate to identify significant differences of gene expression among the four IRS classes (statistical power 0.83), either by ANOVA or Wilcoxon test. The expression of DDR genes was on the whole significantly induced at 24 h after the first fractionated dose (Figure 4A), in accordance with their radioresponsiveness. Notably, FDXR expression showed a 1.6-fold increase, that is very similar to the ~1.7-fold increase reported at the same time point in four breast cancer and 8 HNSCC patients (26). In humans, FDXR expression is upregulated in a dose-dependent manner after irradiation, both ex vivo and in vivo, indicating that FDXR is a good biomarker for radiation exposure and for estimating in vivo dose (46–48). Moreover, FDXR belongs to a genetic signature for the early prediction of hematological acute radiation syndrome (47, 49). Unlike those authors who have found no association between FDXR expression and the hematological acute radiation syndrome in subjects undergoing RT, our work associates FDXR expression level with dermatitis radiation at the completion of RT treatment (G2 vs. G3 grades, adjusted p = 0.096) (Table 3). While the expression level of FDXR gene increased in patients experiencing a high grade of dermatitis radiation, those of SESN1, ATM, XPC, ZMAT3, CDKN1A genes decreased when radiation toxicity was manifested (Table 3). Of course, these findings may be explained by a more complex radiation response than that determined by the DDR genes here analyzed. Indeed, additional genes belonging to different pathways are expected to participate in the whole cellular response to radiation. Emerging evidence suggests that the response to radiation is differently regulated in normal vs. cancer cells/tissues, and even within organism, where maintaining the overall homeostasis is a priority (50, 51). Interestingly, the tight interplay between DDR and immune response seems a key feature shared by systems that differ for higher levels of complexity (51).
By integrating the data of chromosomal radiosensitivity and gene expression we found that IRS classes were associated with the expression level of three DDR genes. In particular, the increased expression of BBC3 and FDXR genes, involved in the apoptotic pathway, is associated with the RS class, whereas the increased expression of GADD45A, regulating cell growth and apoptosis, is associated with the RR class (Table 4). A role of apoptotic pathway in normal tissue radiation toxicity has been previously reported, indeed the T-lymphocyte apoptosis assay significantly predicted differences in late radiation toxicity (52). Also the balance between pro-apoptotic and anti-apoptotic members of the BCL-2 family has a clinical significance on chemotherapy sensitivity and survival (53, 54). In our cohort, the BCL2/BAX ratio resulted associated both with the presence of pruritus 1 month after the completion of RT treatment (Table 3), and with the class of patients classified as radiosensitive by means of IRS index (Table 4). This common function is reinforcing the predictive value of these genes, although further analyses are necessary to support that they are reliable biomarkers of radiation toxicity.
Given the exploratory and pilot nature of our work, we privileged a cohort of RT patients, enrolled, treated and clinically evaluated in the same clinical Institute, irrespective of the tumor site. However, despite our multivariate analysis shows no statistical differences between BC and HNSCC patients (Table 2; Figure 4) future studies would be desirable to validate our results in a new cohort, taking into account additional clinical variables such as breast size, body mass index, alcohol consumption, hypertension, smoking habit, which might be associated to acute skin toxicity (55–57). Moreover, dosimetric data, radiation treatment volumes and doses to specific organs at risk would be important information to be included in future studies. Clearly, since HNSCC patients are less frequent than BC patients and often must undergo chemotherapy either before or concomitant to RT, multicentric studies would be recommended to reach a large sample size for both diseases.
Conclusion
The possibility to identify patients that are sensitive to radiation and at risk of suffering adverse effects would help clinicians in tailoring the best RT protocol and improve patient's QOL. In this prospective cohort study, we found that symptoms of dermatitis radiation, pain, pruritus and fatigue were associated with the expression level of some genes of the DNA-damage response pathway (FDXR, SESN1, XPC, ZMAT3, ATM, BCL2/BAX, and CDKN1A). We also found that fatigue was significantly associated with IRS values; moreover, IRS classes resulted associated with the expression level of BBC3, FDXR, GADD45A, and BCL2/BAX genes.
Of course, radiation-induced side effects comprehend a complex cellular and tissue response that cannot be limited to the expression level of the DDR genes considered in this study, but it is rather regulated by a wide network of gene-interactions. The development of a reproducible and powerful assay to predict individual normal tissue radiosensitivity has been referred to as the “holy grail” of radiotherapy. Although several in vitro assays have been tested to identify reliable biomarkers able to predict normal tissue radiosensitivity, results obtained up to now are not informative enough. In this regard, our multidisciplinary approach can contribute to delineate the genetic features of patients manifesting different grades of radiation-induced toxicity.
Data Availability Statement
The datasets analyzed during the current study are available from the corresponding authors on reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by IOV-IRCCS Ethic Committee. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
LC, DZ, AR, and MM conceived and designed the study. EF, EG, FB, BE, and MR contributed in enrollment of patients, radiotherapy treatment, and clinical follow-up. EP and CP carried out the experiments of gene expression and participated in data evaluation. LB carried out the experiments of G2-chromosomal assay. EC and CR performed the statistical analyses and participated in data evaluation. MM and AR wrote the final manuscript. The definitive supervision of the paper was done by LC. All authors contributed to the initial draft, read, and approved the final manuscript.
Funding
This work was supported by the fund 5 × 1000 IOV- Demabio 1 Demabio 2 (Veneto Institute of Oncology IOV–IRCCS).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are grateful to Prof. Giuseppe Opocher (Veneto Institute of Oncology IOV–IRCCS) for scientific advice and support. We gratefully acknowledge Vito Barbieri (Department of Surgical, Oncological and Gastroenteric Sciences, University of Padova) for support with cell irradiation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00987/full#supplementary-material
Abbreviations
CTCAE, Common Terminology Criteria for Adverse Events; DDR, DNA-Damage Response; DSBs, Double-Strand Breaks; HNSCC, Head and Neck Squamous Cell Carcinoma; HRS, Highly Radiosensitive; IR, Ionizing Radiation; IRS, Individual Radiosensitivity; N, Normal; PBLs, Peripheral Blood Lymphocytes; QOL, Quality of Life; qRT-PCR, quantitative Real-Time PCR; RR, Radioresistant; RS, Radiosensitive; RT, Radiotherapy.
References
1. De Ruysscher D, Niedermann G, Burnet NG, Siva S, Lee AWM, Hegi-Johnson F. Radiotherapy toxicity. Nat Rev Dis Prim. (2019) 5:13. doi: 10.1038/s41572-019-0064-5
2. Burnet NG, Johansen J, Turesson I, Nyman J, Peacock JH. Describing patients' normal tissue reactions: concerning the possibility of individualising radiotherapy dose prescriptions based on potential predictive assays of normal tissue radiosensitivity. Int J Cancer. (1998) 79:606–13. doi: 10.1002/(SICI)1097-0215(19981218)79:6<606::AID-IJC9>3.0.CO;2-Y
3. Azria D, Lapierre A, Gourgou S, De Ruysscher D, Colinge J, Lambin P, et al. Data-based radiation oncology: design of clinical trials in the toxicity biomarkers era. Front Oncol. (2017) 7:1–11. doi: 10.3389/fonc.2017.00083
4. Rieger KE, Hong W-J, Tusher VG, Tang J, Tibshirani R, Chu G. Toxicity from radiation therapy associated with abnormal transcriptional responses to DNA damage. Proc Natl Acad Sci USA. (2004) 101:6635–40. doi: 10.1073/pnas.0307761101
5. Foray N, Colin C, Bourguignon M. 100 years of individual radiosensitivity: how we have forgotten the evidence. Radiology. (2012) 264:627–31. doi: 10.1148/radiol.12112560
6. The 2007 recommendations of the International Commission on Radiological Protection. ICRP publication 103. Ann ICRP. (2007) 37:1–332. doi: 10.1016/j.icrp.2007.10.003
7. AGIR. Human Radiosensitivity. Report of the Independent Advisory Group on Ionising Radiation. Chilton, Doc HPA, RCE-21 (2013).
8. Barnett GC, West CML, Dunning AM, Elliott RM, Coles CE, iPharoah PDP, et al. Normal tissue reactions to radiotherapy: towards tailoring treatment dose by genotype. Nat Rev Cancer. (2009) 9:134–42. doi: 10.1038/nrc2587
9. Aparicio T, Baer R, Gautier J. DNA double-strand break repair pathway choice and cancer. DNA Repair. (2014) 19:169–75. doi: 10.1016/j.dnarep.2014.03.014
10. Harper JW, Elledge SJ. The DNA damage response: ten years after. Mol Cell. (2007) 28:739–45. doi: 10.1016/j.molcel.2007.11.015
11. Goldstein M, Kastan MB. The DNA damage response: implications for tumor responses to radiation and chemotherapy. Annu Rev Med. (2015) 66:129–43. doi: 10.1146/annurev-med-081313-121208
12. Durante M, Bedford JS, Chen DJ, Conrad S, Cornforth MN, Natarajan AT, et al. From DNA damage to chromosome aberrations: joining the break. Mutat Res Toxicol Environ Mutagen. (2013) 756:5–13. doi: 10.1016/j.mrgentox.2013.05.014
13. Borgmann K, Hoeller U, Nowack S, Bernhard M, Röper B, Brackrock S, et al. Individual radiosensitivity measured with lymphocytes may predict the risk of acute reaction after radiotherapy. Int J Radiat Oncol. (2008) 71:256–64. doi: 10.1016/j.ijrobp.2008.01.007
14. Barber JB., Burrill W, Spreadborough AR, Levine E, Warren C, Kiltie AE, et al. Relationship between in vitro chromosomal radiosensitivity of peripheral blood lymphocytes and the expression of normal tissue damage following radiotherapy for breast cancer. Radiother Oncol. (2000) 55:179–86. doi: 10.1016/S0167-8140(99)00158-9
15. Pantelias GE, Terzoudi GI. A standardized G2-assay for the prediction of individual radiosensitivity. Radiother Oncol. (2011) 101:28–34. doi: 10.1016/j.radonc.2011.09.021
16. Bryant PE, Gray L, Riches AC, Steel CM, Finnon P, Howe O, et al. The G2 chromosomal radiosensitivity assay. Int J Radiat Biol. (2002) 78:863–6. doi: 10.1080/09553000210144484
17. Scott D, Spreadborough AR, Jones LA, Roberts SA, Moore CJ. Chromosomal radiosensitivity in G2-phase lymphocytes as an indicator of cancer predisposition. Radiat Res. (1996) 145:3–16. doi: 10.2307/3579189
18. Lobachevsky P, Leong T, Daly P, Smith J, Best N, Tomaszewski J, et al. Compromized DNA repair as a basis for identification of cancer radiotherapy patients with extreme radiosensitivity. Cancer Lett. (2016) 383:212–9. doi: 10.1016/j.canlet.2016.09.010
19. De Ruyck K, Van Eijkeren M, Claes K, Morthier R, De Paepe A, Vral A, et al. Radiation-induced damage to normal tissues after radiotherapy in patients treated for gynecologic tumors: association with single nucleotide polymorphisms in XRCC1, XRCC3, and OGG1 genes and in vitro chromosomal radiosensitivity in lymphocytes. Int J Radiat Oncol. (2005) 62:1140–9. doi: 10.1016/j.ijrobp.2004.12.027
20. Alam A, Mukhopadhyay ND, Ning Y, Reshko LB, Cardnell RJG, Alam O, et al. A preliminary study on racial differences in HMOX1, NFE2L2, and TGFβ1 gene polymorphisms and radiation-induced late normal tissue toxicity. Int J Radiat Oncol. (2015) 93:436–43. doi: 10.1016/j.ijrobp.2015.05.049
21. Oliva D, Nilsson M, Strandéus M, Andersson B-Å, Sharp L, Laytragoon-Lewin N, et al. Individual genetic variation might predict acute skin reactions in women undergoing adjuvant breast cancer radiotherapy. Anticancer Res. (2018) 38:6763–70. doi: 10.21873/anticanres.13047
22. Borghini A, Vecoli C, Mercuri A, Petruzzelli MF, D'Errico MP, Portaluri M, et al. Genetic risk score and acute skin toxicity after breast radiation therapy. Cancer Biother Radiopharm. (2014) 29:267–72. doi: 10.1089/cbr.2014.1620
23. Svensson JP, Stalpers LJ, Esveldt-van Lange RE, Franken NA, Haveman J, Klein B, et al. Analysis of gene expression using gene sets discriminates cancer patients with and without late radiation toxicity. PLoS Med. (2006) 3:e422. doi: 10.1371/journal.pmed.0030422
24. Finnon P, Kabacik S, MacKay A, Raffy C, A'Hern R, Owen R, et al. Correlation of in vitro lymphocyte radiosensitivity and gene expression with late normal tissue reactions following curative radiotherapy for breast cancer. Radiother Oncol. (2012) 105:329–36. doi: 10.1016/j.radonc.2012.10.007
25. Badie C, Dziwura S, Raffy C, Tsigani T, Alsbeih G, Moody J, et al. Aberrant CDKN1A transcriptional response associates with abnormal sensitivity to radiation treatment. Br J Cancer. (2008) 98:1845–51. doi: 10.1038/sj.bjc.6604381
26. O'Brien G, Cruz-Garcia L, Majewski M, Grepl J, Abend M, Port M, et al. FDXR is a biomarker of radiation exposure in vivo. Sci Rep. (2018) 8:684. doi: 10.1038/s41598-017-19043-w
27. Herskind C, Talbot CJ, Kerns SL, Veldwijk MR, Rosenstein BS, West CML. Radiogenomics: a systems biology approach to understanding genetic risk factors for radiotherapy toxicity? Cancer Lett. (2016) 382:95–109. doi: 10.1016/j.canlet.2016.02.035
28. Pavlopoulou A, Bagos PG, Koutsandrea V, Georgakilas AG. Molecular determinants of radiosensitivity in normal and tumor tissue: a bioinformatic approach. Cancer Lett. (2017) 403:37–47. doi: 10.1016/j.canlet.2017.05.023
29. Terzoudi GI, Hatzi VI, Barszczewska K, Manola KN, Stavropoulou C, Angelakis P, et al. G2-checkpoint abrogation in irradiated lymphocytes: a new cytogenetic approach to assess individual radiosensitivity and predisposition to cancer. Int J Oncol. (2009) 35:1223–30. doi: 10.3892/ijo_00000439
30. Mognato M, Celotti L. Modeled microgravity affects cell survival and HPRT mutant frequency, but not the expression of DNA repair genes in human lymphocytes irradiated with ionising radiation. Mutat Res. (2005) 578:417–29. doi: 10.1016/j.mrfmmm.2005.06.011
31. Girardi C, De Pittà C, Casara S, Sales G, Lanfranchi G, Celotti L, et al. Analysis of miRNA and mRNA expression profiles highlights alterations in ionizing radiation response of human lymphocytes under modeled microgravity. PLoS ONE. (2012) 7:e31293. doi: 10.1371/journal.pone.0031293
32. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. (2001) 29:e45. doi: 10.1093/nar/29.9.e45
33. Trotti A, Bellm LA, Epstein JB, Frame D, Fuchs HJ, Gwede CK, et al. Mucositis incidence, severity and associated outcomes in patients with head and neck cancer receiving radiotherapy with or without chemotherapy: a systematic literature review. Radiother Oncol. (2003) 66:253–62. doi: 10.1016/S0167-8140(02)00404-8
34. Peponi E, Glanzmann C, Willi B, Huber G, Studer G. Dysphagia in head and neck cancer patients following intensity modulated radiotherapy (IMRT). Radiat Oncol. (2011) 6:1. doi: 10.1186/1748-717X-6-1
35. Ciammella P, Podgornii A, Galeandro M, Micera R, Ramundo D, Palmieri T, et al. Toxicity and cosmetic outcome of hypofractionated whole-breast radiotherapy: predictive clinical and dosimetric factors. Radiat Oncol. (2014) 9:1–10. doi: 10.1186/1748-717X-9-97
36. Smart V, Curwen GB, Whitehouse CA, Edwards A, Tawn EJ. Chromosomal radiosensitivity: a study of the chromosomal G2assay in human blood lymphocytes indicating significant inter-individual variability. Mutat Res. (2003) 528:105–10. doi: 10.1016/S0027-5107(03)00076-9
37. Vral A, Thierens H, Baeyens A, De Ridder L. Chromosomal aberrations and in vitro radiosensitivity: intra-individual versus inter-individual variability. Toxicol Lett. (2004) 149:345–52. doi: 10.1016/j.toxlet.2003.12.044
38. Valli M, Cima S, Gaudino D, Cartolari R, Deantonio L, Frapolli M, et al. Skin and lung toxicity in synchronous bilateral breast cancer treated with volumetric-modulated arc radiotherapy: a mono-institutional experience. Clin Transl Oncol. (2019). doi: 10.1007/s12094-019-02077-z. [Epub ahead of print].
39. Hegedus F, Mathew LM, Schwartz RA. Radiation dermatitis: an overview. Int J Dermatol. (2017) 56:909–14. doi: 10.1111/ijd.13371
40. Hopwood P, Haviland JS, Sumo G, Mills J, Bliss JM, Yarnold JR. Comparison of patient-reported breast, arm, and shoulder symptoms and body image after radiotherapy for early breast cancer: 5-year follow-up in the randomised Standardisation of Breast Radiotherapy (START) trials. Lancet Oncol. (2010) 11:231–40. doi: 10.1016/S1470-2045(09)70382-1
41. Guogyte K, Plieskiene A, Ladygiene R, Vaisiunas Ž, Sevriukova O, Janušonis V, et al. Assessment of correlation between chromosomal radiosensitivity of peripheral blood lymphocytes after in vitro irradiation and normal tissue side effects for cancer patients undergoing radiotherapy. Genome Integr. (2017) 8:2. doi: 10.4103/2041-9414.198907
42. Andreassen CN, Alsner J, Overgaard M, Sørensen FB, Overgaard J. Risk of radiation-induced subcutaneous fibrosis in relation to single nucleotide polymorphisms in TGFB1, SOD2, XRCC1, XRCC3, APEX and ATM–a study based on DNA from formalin fixed paraffin embedded tissue samples. Int J Radiat Biol. (2006) 82:577–86. doi: 10.1080/09553000600876637
43. Andreassen CN, Alsner J. Genetic variants and normal tissue toxicity after radiotherapy: a systematic review. Radiother Oncol. (2009) 92:299–309. doi: 10.1016/j.radonc.2009.06.015
44. Goodarzi AA, Jeggo PA. Irradiation induced foci (IRIF) as a biomarker for radiosensitivity. Mutat Res Mol Mech Mutagen. (2012) 736:39–47. doi: 10.1016/j.mrfmmm.2011.05.017
45. Amundson SA, Do KT, Shahab S, Bittner M, Meltzer P, Trent J, et al. Identification of potential mRNA biomarkers in peripheral blood lymphocytes for human exposure to ionizing radiation. Radiat Res. (2000) 154:342–6. doi: 10.1667/0033-7587(2000)154[0342:iopmbi]2.0.co;2
46. Manning G, Kabacik S, Finnon P, Bouffler S, Badie C. High and low dose responses of transcriptional biomarkers in ex vivo X-irradiated human blood. Int J Radiat Biol. (2013) 89:512–22. doi: 10.3109/09553002.2013.769694
47. Port M, Majewski M, Herodin F, Valente M, Drouet M, Forcheron F, et al. Validating baboon ex vivo and in vivo radiation-related gene expression with corresponding human data. Radiat Res. (2018) 189:389–98. doi: 10.1667/RR14958.1
48. Abend M, Badie C, Quintens R, Kriehuber R, Manning G, Macaeva E, et al. Examining radiation-induced in vivo and in vitro gene expression changes of the peripheral blood in different laboratories for biodosimetry purposes: first RENEB gene expression study. Radiat Res. (2016) 185:109–23. doi: 10.1667/RR14221.1
49. Port M, Herodin F, Valente M, Drouet M, Lamkowski A, Majewski M, et al. First generation gene expression signature for early prediction of late occurring hematological acute radiation syndrome in baboons. Radiat Res. (2016) 186:39–54. doi: 10.1667/RR14318.1
50. Georgakilas AG, Pavlopoulou A, Louka M, Nikitaki Z, Vorgias CE, Bagos PG. Michalopoulos I. Emerging molecular networks common in ionizing radiation, immune and inflammatory responses by employing bioinformatics approaches. Cancer Lett. (2015) 368:164–72. doi: 10.1016/j.canlet.2015.03.021
51. Pateras IS, Havaki S, Nikitopoulou X, Vougas K, Townsend PA, Panayiotidis MI, et al. The DNA damage response and immune signaling alliance: is it good or bad? Nature decides when and where. Pharmacol Ther. (2015) 154:36–56. doi: 10.1016/j.pharmthera.2015.06.011
52. Ozsahin M, Crompton NEA, Gourgou S, Kramar A, Li L, Shi Y, et al. CD4 and CD8 T-lymphocyte apoptosis can predict radiation-induced late toxicity: a prospective study in 399 patients. Clin Cancer Res. (2005) 11:7426–33. doi: 10.1158/1078-0432.CCR-04-2634
53. Avgeris M, Stamati L, Kontos CK, Piatopoulou D, Marmarinos A, Xagorari M, et al. BCL2L12 improves risk stratification and prediction of BFM-chemotherapy response in childhood acute lymphoblastic leukemia. Clin Chem Lab Med. (2018) 56:2104–18. doi: 10.1515/cclm-2018-0507
54. Del Principe MI, Bo MD, Bittolo T, Buccisano F, Rossi FM, Zucchetto A, et al. Clinical significance of bax/bcl-2 ratio in chronic lymphocytic leukemia. Haematologica. (2016) 101:77–85. doi: 10.3324/haematol.2015.131854
55. Kole AJ, Kole L, Moran MS. Acute radiation dermatitis in breast cancer patients: challenges and solutions. Breast Cancer Targets Ther. (2017) 9:313–23. doi: 10.2147/BCTT.S109763
56. Pignol JP, Olivotto I, Rakovitch E, Gardner S, Sixel K, Beckham W, et al. A multicenter randomized trial of breast intensity-modulated radiation therapy to reduce acute radiation dermatitis. J Clin Oncol. (2008) 26:2085–92. doi: 10.1200/JCO.2007.15.2488
57. Twardella D, Popanda O, Helmbold I, Ebbeler R, Benner A, Von Fournier D, et al. Personal characteristics, therapy modalities and individual DNA repair capacity as predictive factors of acute skin toxicity in an unselected cohort of breast cancer patients receiving radiotherapy. Radiother Oncol. (2003) 69:145–53. doi: 10.1016/S0167-8140(03)00166-X
Keywords: radiotherapy, adverse effects, chromosomal radiosensitivity, gene expression, association analysis
Citation: Palumbo E, Piotto C, Calura E, Fasanaro E, Groff E, Busato F, El Khouzai B, Rigo M, Baggio L, Romualdi C, Zafiropoulos D, Russo A, Mognato M and Corti L (2019) Individual Radiosensitivity in Oncological Patients: Linking Adverse Normal Tissue Reactions and Genetic Features. Front. Oncol. 9:987. doi: 10.3389/fonc.2019.00987
Received: 24 June 2019; Accepted: 16 September 2019;
Published: 01 October 2019.
Edited by:
Tiziana Rancati, National Tumor Institute, ItalyReviewed by:
Simeng Suy, Georgetown University, United StatesMaria Carmen De Santis, National Tumor Institute, Italy
Copyright © 2019 Palumbo, Piotto, Calura, Fasanaro, Groff, Busato, El Khouzai, Rigo, Baggio, Romualdi, Zafiropoulos, Russo, Mognato and Corti. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Antonella Russo, antonella.russo@unipd.it; Maddalena Mognato, maddalena.mognato@unipd.it
†Present address: Elisa Palumbo, Department of Molecular Medicine, University of Padua, Padua, Italy
Celeste Piotto, European Molecular Biology Laboratory Australia, Australian Regenerative Medicine Institute, Monash University, Melbourne, VIC, Australia
Elena Fasanaro, Otolaryngology Unit, Cittadella Hospital, Cittadella, Italy
Michele Rigo, IRCCS Sacro Cuore Don Calabria Hospital, Cancer Care Center, Verona, Italy
‡These authors have contributed equally to this work