Comparative analysis of three functional predictive assays in lymphocytes of patients with breast and gynaecological cancer treated by radiotherapy.

PURPOSE: There is a need to develop predictive tests that would allow identifying cancer patients with a high risk of developing side effects to radiotherapy. We compared the predictive value of three functional assays: the G(0) aberration assay, the G(2) aberration assay and the alkaline comet assay in lymphocytes of breast cancer and gynaecological cancer patients.
MATERIAL AND METHODS: Peripheral blood was collected from 35 patients with breast cancer and 34 patients with gynaecological cancer before the onset of therapy. Chromosomal aberrations were scored in lymphocytes irradiated in the G(0) or G(2) phase of the cell cycle. DNA repair kinetics was performed with the alkaline comet assay following irradiation of unstimulated lymphocytes. The results were compared with the severity of early and late side effects to radiotherapy.
RESULTS: No correlation was observed between the results of the assays and the severity of side effects. Moreover, each assay identified different patients as radiosensitive.
CONCLUSIONS: There is no simple correlation between the in vitro sensitivity of lymphocytes and the risk of developing early and late side effects.

Purpose

Radiotherapy patients show different degrees of adverse side effects and it is assumed that the risk of developing side effects is related to the genetically determined, intrinsic sensitivity of the patient [1-3]. From the perspective of the radiation protection of the patient it would be desirable to develop tests predicting the individual risk of developing side effects. In the 1990s numerous studies have been performed with the aim of finding a genetic fingerprint of intrinsic radiosensitivity (reviewed in [2]). More recently it was realised that the genetics behind the intrinsic radiosensitivity is very complex and its characterisation requires a systems biology approach rather than the search for genetic changes in a few loci [4].

An alternative approach in the search for a marker of individual radiosensitivity is a functional assay. Here, surrogate normal tissue – most often peripheral blood lymphocytes (PBL) – is drawn from a patient and the in vitro radiosensitivity is determined by assays that allow quantifying the level of DNA damage or the capacity of DNA repair. Notwithstanding some promising results [5], the bulk of data is controversial and the reasons for the inconsistency are not understood [6-8]. Also, it is clear that a single assay will not have a satisfactory predictive value [9]. But will combining different assays solve the problem? Do different functional assays complement each other?

The present study was undertaken in order to compare three functional assays: the G0 aberration assay, the G2 aberration assay and the alkaline comet assay. The assays were compared using PBL of patients treated by external beam radiotherapy for breast cancer and gynaecological cancers. PBL were collected before radiotherapy, exposed in vitro to a single dose of radiation and assessed for the in vitro sensitivity. The results were correlated with early and late side effects. Although the in vitro sensitivity of PBL showed marked inter-individual variabilities, no correlations were detected, neither between the assays nor with the extent of side effects.

Material and methods

Cancer patients

Breast cancer (BC) patients: the group consisted of 35 patients with breast cancer staged IIA to IIIB (median age 57). All patients underwent radical mastectomy and chemotherapy prior to radiotherapy. Conformal radiotherapy with a 3D planning system was performed using the modified inverse hockey stick technique [10].

Gynaecological cancer (GC) patients: the group consisted of 34 patients with cancer of endometrium (18 patients) and cervix (16 patients) of different stages (median age 63). All patients with endometrial cancer underwent surgery and were subsequently treated by external beam radiotherapy and high dose rate (HDR) brachytherapy without chemotherapy. Patients with cervix cancer received different forms of therapy: 5 patients underwent surgery, 12 received chemotherapy and all were treated by external beam radiotherapy plus HDR brachytherapy. The radiotherapy of gynaecological cancers was based on a 2D planning system. The patients were treated with the isocentric 4 field box technique.

The study was approved by the local ethical committee.

Lymphocyte collection and in vitro irradiation

All blood samples were collected in the Holy Cross Cancer Centre from patients before the onset of radiotherapy, immediately aliquoted in 1.5 ml Eppendorf tubes and cooled in ice water (0.8°C) in order to maintain reproducible irradiation conditions. All irradiations were carried out at the Holy Cross Cancer Centre using 6 MeV photons generated by a Siemens therapeutic linear accelerator. Blood samples were always irradiated at a distance of 80 cm from the source in a cold paraffin bolus to achieve electron equilibrium. A blood sample from one donor was always processed for all three assays.

Blood samples for the G0 assay and comet assay were irradiated with 2 Gy of photons shortly after collection and transported in an ice box back to the laboratory (located in a different part of the town) where they was immediately processed for the assays as described below. Blood for the G2 assay was not irradiated, but transported to the laboratory in an ice box and cultured as described below. After 69 hours cultures were cooled 0.8°C (melting ice), transported in an ice box to the Holy Cross Cancer Centre and irradiated with a dose of 1 Gy. Cooled cultures were transported back to the laboratory, placed back into the incubator and processed for the chromosomal aberration G2 assay as described below.

Chromosomal aberration assay

0.5 ml of whole blood was added to 4.5 ml of complete culture medium as described in [11]. 5-bromo-2-deoxyuridine (BrdU) was used only for the G0 aberration assay. Three test tubes were set up from each patient: one with blood irradiated for the G0 aberration assay, one with blood to be irradiated for the G2 aberration assay and one with unexposed blood. Test tubes were incubated at 37°C (5% CO2, 95% humidity). For the G0 aberration assay colcemid (Gibco 15212-046, 0.08 µg/ml) was added after 45 hours of culture time and cells were harvested 3 hours later. For the G2 aberration assay colcemid was added one hour after placing the culture back into the incubator (corresponding to about 70th hour of culture time) and cells were harvested 2 hours later.

Cells were harvested and fixed as described in [11]. Slides for the G0 aberration assay were stained with the fluorescence plus Giemsa (FPG) technique that allows to identify cells in the first postirradiation division [12]. Slides for the G2 assay were stained for 10 minutes with a 10% Giemsa solution. No cell proliferation control was needed because of the short time between exposure and harvest.

For metaphase analysis the slides were coded and the aberrations scored randomly in 100 mitoses. All aberrations (G0 assay: mainly dicentrics, rings and acentric fragments; G2 assay: mainly acentric fragments) were scored and expressed as sum of aberrations. Gaps (chromosome arm discontinuities that are shorter than the width of the chromosome arm) were not included in the analysis. Spontaneous aberrations were scored on slides harvested after 72 hours of culture.

Comet assay

Comet assay was performed according to the method of Singh et al. [13] and as described in detail in [14]. 50 cells per point were analysed for the Olive tail moment (OTM) using the CASP software [14]. Repair capacity was assessed by two ways. 1: the level of damage after 0 min repair was considered as 100% and the residual damage after specific times of repair were expressed as percent of the initial level. Each repair kinetic curve was used to calculate the area under the curve (AUC) by fitting the curve to the exponential decay function y = ae-bx with the help of SigmaPlot (Systat Software Inc., USA®). AUC was calculated with the help of GraphCalc (http://www.graphcalc.com). 2: For each repair time point an OTM value was classified as high (being in the top 33.3% of values), intermediate (middle 33.3% of values) and low (bottom 33.3% of values) with respect to pooled values from all patients (per patient group). This approach allowed constructing cluster heat maps, similarly as used for analysis of gene expression patterns [15].

Clinical classification of normal tissue reactions

Acute and late normal tissue reactions to radiotherapy were scored using the Radiation Therapy Oncology Group/European Organization for Research and Treatment of Cancer (EORTC/RTOG) grading system. Early reactions were scored during the treatment or immediately thereafter. Late reactions were scored periodically, starting six months post treatment. A lack of reaction was classified as grade 0, a mild reactions as grade 1, a moderate reaction as grade 2 and severe reaction as grade 3 or 4. For breast cancer patients reactions of the skin, larynx, lung and heart were recorded. For gynaecological patients reactions were scored in skin, colon, bladder and the peripheral blood (leucopoenia).

Statistical methods

The Mann-Whitney Rank Sum Test was used to compare the frequencies of aberrations and the AUC values between groups of patients. Relationships between the data were analysed by the Pearson product moment correlation. Relation between the aberration frequencies and the characteristics of the DNA repair kinetics were analysed by constructing cluster heat maps [15].

Results

Chromosomal aberrations

Spontaneous aberrations were all of chromatid type and the mean aberrations frequency was for each group was 1.5. The mean G0 aberration frequency in lymphocytes of BC patients was 40.7 and in lymphocytes of GC patients 32.4. This difference is statistically significant. The mean G2 aberrations frequencies were 37.5 for the BC patients and 30.6 for GC patients, the difference not being significant.

The results of the comet assay are shown as individual repair kinetic curves in Fig. 1. The median AUC values were 13.66 for BC patients and 14.80 for GC patients. This difference is not statistically significant. From the curves shown in Fig. 1 it is clear that the level of residual damage showed a stronger individual heterogeneity in lymphocytes of GC patients than in lymphocytes of BC patients.

Fig. 1

Individual repair kinetic curves obtained with the comet assay. A) Breast cancer patients, B) gynaecological cancer patients. AUC: area under curve. Thick line represents the median curve

Correlation between aberrations and comet assay

The correlation analyses between the different endpoints are shown in Fig. 2. No correlations were observed between G0 and G2 aberrations (panels A), G0 aberrations and AUC (panels B) and G2 aberrations and AUC (panels C). In order to analyse in more detail the relationship between aberrations and the kinetics of DNA repair we constructed cluster heat maps where the OTM values per repair time point and patient group were stratified into three groups: high, medium and low. A grey-scale hue was assigned to each level group. The alignments of grey scale patterns and sorted G0 aberration frequencies are shown in Fig. 3. If a given level of chromosomal aberrations correlated with a particular pattern of grey scale hues, this would be optically discernible. However, neither for the BC nor for the GC patients could a pattern be identified. The same situation was seen for G2 aberrations (not shown).

Fig. 2

Correlations between the frequencies of G0 aberrations, G2 aberrations and AUC values. Top panels: breast cancer patients, bottom panels: gynaecological cancer patients

Fig. 3

Relationship between the G0 aberration frequency in lymphocytes of breast cancer patients (left panel) and gynaecological cancer patients (right panel) and the kinetic of DNA repair measured with the comet assay. The grey hue is indicative of a high, medium or low value of OTM (Olive tail moment) with respect to pooled values per patient group and repair time point. Patients are sorted according to increasing aberration frequency

Correlation between aberrations, comet assay and side effects of therapy

The levels of chromosomal aberrations and the AUC values were plotted as function of the degree of severity of side effects. Side effects observed in different tissues were pooled and expressed as cumulative grades. The results are presented in Fig. 4 for early and late effects and both patient groups. Due to a larger number of tissues/organs at risk, the values of the cumulative grades were higher for GC patients as compared to BC patients. The severity of neither early not late effects correlated with the frequency of the aberrations or the AUC.

Fig. 4

Relationships between radiation-induced aberration frequencies plus AUC values and early (black symbols) and late (open symbols) reactions to radiotherapy, expressed as cumulative grades. Top panels: breast cancer patients, bottom panels: gynaecological cancer patients. Panels A: G0 aberrations, panels B: G2 aberrations, panels C: AUC

Discussion

The search for predictive tests for acute and late reactions to radiotherapy is based on the assumption that the major factor associated with a high risk of developing side effects is the genetically determined capacity to detect and repair DNA damage [1]. Although methods are now available to directly analyse genes and proteins involved in biological responses to radiation [16, 17], we chose to use three functional assays in PBL of patients because they integrate the influence of many genes that potentially play a role in the cellular radiation sensitivity.

We did not detect any relationship between the severity of acute and late effects and the in vitro sensitivity of PBL as assessed by the G0 aberration test, the G2 aberration test and the alkaline comet assay. This negative result is not totally unexpected in view of the conflicting results published by others. A number of groups have attempted to find a correlation between the chromosomal radiosensitivity of PBL and the reactions of healthy tissues to radiotherapy (reviewed in [6, 8]). A positive relationship was reported in the case of patients with cervix cancer [18-20], head and neck cancer [21], breast cancer [22] and prostate cancer [23]. In contrast, Slonina et al. [24] and Lisowska et al. [11] did not observe any correlation between the frequency of radiation-induced cytogenetic changes in PBL and of normal tissue reactions in patients with head and neck cancer. Similarly negative results were published for breast cancer [25] and prostate cancer [7]. The reasons for the conflicting results are not clear.

The application of the alkaline comet assay as a predictive test yielded similarly disagreeing results. Müller et al. [5] observed a delayed repair of DNA damage in PBL of cancer patients with various tumour types that developed severe side effects to radiotherapy. Similar results were published by Gabelova et al. [26] for cervical cancer and Alapetite et al. [27] for breast cancer. In contrast, Djuzenova et al. [28] did not find an impaired level of initial DNA damage, residual DNA damage and repair kinetics in PBL of breast cancer patients that developed severe side effects. This result is in accordance with our data. It is clear that the comet assay still requires validation as a test to detect individual radiosensitivity.

Our study not only showed a lack of correlation between the in vitro sensitivity of PBL and the severity of side effects, but also between the in vitro assays. That the results of the G0 and the G2 aberration assays do not correlate may not be surprising in view of the fact that different DNA double strand break repair pathways dominate during the G0 and G2 phases of the cell cycle [29]. A more intriguing result appears to be the lack of relationship between the G0 aberrations and the comet assay, since both assays were performed on unstimulated lymphocytes. The relationship was analysed by comparing the frequencies of aberrations with the AUC values. Because no correlation was observed, we performed a more detailed comparison by aligning the aberration frequencies against grey-scale hue diagrams that characterise the level of DNA damage following 0, 15, 30, 60 and 180 minutes of repair. This approach, called a cluster heat map was introduced by [15] to detect possible patterns in complex data sets. The rational for applying the heat map approach here was to test the possibility that an association exists between the aberrations frequency and a particular time point of DNA repair kinetics, measured by the comet assay. Such association would be visible by clustering of a hue in a heat map aligned according to increasing aberration frequency. However, no such association could be detected. An explanation of the lack of correlation between the frequency of aberrations and the comet assay results may be that, as suggested by Speit et al. [30], an enhanced frequency of cytogenetic damage in PBL results from an impaired fidelity of DNA repair rather than the kinetics of damage removal. Changes in the fidelity of DNA repair cannot be visualised by the comet assay.

In conclusion, our study failed to reveal a relationship between the risk of developing side effects to radiotherapy and the in vitro sensitivity of PBL as assessed by three functional assays. Each in vitro test identified different donors as radiosensitive. This could be either due to a different specificity of a test with respect to the DNA damage processing pathway or to other factors such as the physiological state of the blood donor or random experimental fluctuations. Thus, the usefulness of the functional tests for predicting healthy tissue reactions requires further investigation.