Skip to main content

Prostate cancer treated with brachytherapy; an exploratory study of dose-dependent biomarkers and quality of life



Low-dose-rate permanent prostate brachytherapy (PPB) is an attractive treatment option for patients with localised prostate cancer with excellent outcomes. As standard CT-based post-implant dosimetry often correlates poorly with late treatment-related toxicity, this exploratory (proof of concept) study was conducted to investigate correlations between radiation − induced DNA damage biomarker levels, and acute and late bowel, urinary, and sexual toxicity.


Twelve patients treated with 125I PPB monotherapy (145Gy) for prostate cancer were included in this prospective study. Post-implant CT based dosimetry assessed the minimum dose encompassing 90% (D90%) of the whole prostate volume (global), sub-regions of the prostate (12 sectors) and the near maximum doses (D0.1cc, D2cc) for the rectum and bladder. Six blood samples were collected from each patient; pre-treatment, 1 h (h), 4 h, 24 h post-implant, at 4 weeks (w) and at 3 months (m). DNA double strand breaks were investigated by staining the blood samples with immunofluorescence antibodies to γH2AX and 53BP1 proteins (γH2AX/53BP1). Patient self-scored quality of life from the Expanded Prostate Cancer Index Composite (EPIC) were obtained at baseline, 1 m, 3 m, 6 m, 9 m, 1 year (y), 2y and 3y post-treatment. Spearman’s correlation coefficients were used to evaluate correlations between temporal changes in γH2AX/53BP1, dose and toxicity.


The minimum follow up was 2 years. Population mean prostate D90% was 144.6 ± 12.1 Gy and rectal near maximum dose D0.1cc = 153.0 ± 30.8 Gy and D2cc = 62.7 ± 12.1 Gy and for the bladder D0.1cc = 123.1 ± 27.0 Gy and D2cc = 70.9 ± 11.9 Gy. Changes in EPIC scores from baseline showed high positive correlation between acute toxicity and late toxicity for both urinary and bowel symptoms. Increased production of γH2AX/53BP1 at 24 h relative to baseline positively correlated with late bowel symptoms. Overall, no correlations were observed between dose metrics (prostate global or sector doses) and γH2AX/53BP1 foci counts.


Our results show that a prompt increase in γH2AX/53BP1foci at 24 h post-implant relative to baseline may be a useful measure to assess elevated risk of late RT − related toxicities for PPB patients. A subsequent investigation recruiting a larger cohort of patients is warranted to verify our findings.


Prostate cancer is the most commonly diagnosed male cancer in the UK with an incidence of 103 cases per 100,000 (European age-standardised rates, 2010) and almost 47,300 patients were diagnosed in the UK alone in 2013 (Cancer Research UK). Initial treatment options include active surveillance [1], androgen deprivation therapy (ADT) [2], radical prostatectomy [3], external beam radiotherapy (EBRT) and brachytherapy with a patient’s Gleason score and tumour staging primarily dictating which treatment modality is optimal. For localised prostate cancer (T1-T2, Gleason score ≤ 7 and PSA < 20 ng/ml), low-dose-rate permanent prostate brachytherapy (PPB) is an attractive choice of treatment due to its convenience, favourable dose distribution to normal tissues, high success rate with regards to both progression free and overall survival [4, 5] and its potential higher likelihood of potency preservation when compared to prostatectomy [6].

Table 1 Patient basic baseline characteristics

Current post-implant dosimetric evaluation for PPB is calculated using CT derived models while patient outcomes and quality of life (QoL), in the form of function and bother, can be measured using patient − reported outcome measures such as the Expanded Prostate Cancer Index Composite (EPIC) questionnaire [7]. The use of blood − derived radiation biomarkers have provided in-vivo detection of irradiation for large volume acute radiation exposures [8,9,10], however their utility for small volume prolonged PPB exposures within brachytherapy has yet to be established. The histone H2AX is phosphorylated to γH2AX at the sites of radiation induced DNA double − strand breaks (DSBs) [11], where it co-localises with the DNA repair protein 53BP1. Loss of γH2AX/53 PB1 over time correlates well with the kinetics of DSBs repair, and has been extensively used to probe the kinetics of repair in a range of conditions [11,12,13,14,15]. γH2AX/53BP1 form distinct nuclear foci which can be readily detected and quantified using immunofluorescence microscopy. For acute doses, a minimum sensitivity down to 5 mGy was detected in patients undergoing a chest CT by comparing pre- and post- CT γH2AX foci numbers in peripheral blood lymphocytes [16]. With low background numbers of DSBs and a well characterised radiation response, peripheral blood lymphocytes are well established as an ideal ex-vivo cell model to use for biological dosimetry.

The level of detectable DNA damage in the lymphocyte population, which contains cells constantly circulating throughout the peripheral blood and also stationary in reservoirs throughout the body such as lymph nodes and spleen, may not correlate simply to the equivalent whole body dose delivered to the patient. There is also increasing evidence that doses received by the prostate and surrounding organs such as the bladder and rectum will display varying degrees of inter- and intra-patient variability [17, 18]. Whether radiation induced DNA damage to the lymphocyte population is dependent on where doses are targeted rather than simply the total dose delivered to all tissues is currently poorly understood.

The primary objective of this study was to investigate and quantify, for the first time, the induction and persistence of γH2AX/53BP1 foci in peripheral blood lymphocytes in patients undergoing low-dose-rate PPB for prostate cancer. Furthermore, we tested the hypothesis that levels of biological DNA damage biomarkers relative to baseline would correlate with patient reported QoL outcome measures determined using the EPIC questionnaire.


Patient cohort

Following ethical approval, 12 patients eligible for treatment with PPB for prostate cancer were identified and recruited for this study. Patients’ baseline characteristics are presented in Table 1. Written information was provided and potential subjects were given at least 7 days to decide upon participation. Inclusion criteria were as follows; patients with histologically confirmed adenocarcinoma of the prostate, who elected for low-dose-rate permanent prostate brachytherapy, WHO performance status 0 – 2, life expectancy of at least 12 months, age ≥ 18 years and a willingness to co-operate with follow-up. Exclusion criteria were; chronic bowel disorder or previous bowel surgery, rheumatoid arthritis, auto immune disease, an active bacterial, viral or fungal infection, renal impairment, diabetes or previous radiotherapy.


Prior to implantation, all patients underwent a detailed prostate volume study using a bi-planar trans-rectal probe 8848 (BK Medical Systems Inc. Washington, MA). All implants were planned using 125I stranded radioactive seeds with a reference air kerma rate of 0.458 – 0.496 U (U –μGyh−1m2) [19, 20]. The half-life of Iodine 125 is 59.43 days. The American Association of Physicists in Medicine Task Group 43 (TG-43 update) formalism was used in the planning and calculation of the final dosimetry [21]. All treatment plans were created to meet the objectives and constraints recommended by ESTRO/EAU/EORTC guidelines [22]. A Variseed 8.0 (Varian Medical Systems, Charlottesville, VA) treatment planning system (TPS) was used for both pre- and post- implant planning. The minimum prescription dose was 145 Gy for all patients, representing the planned dose that will encompass the entire target volume. The dose inside the target volume will be >145Gy, therefore planning constraints are used; no more that 55% of the target volume may receive ≥150% of dose (217.5Gy) and the volume receiving ≥200% (290Gy) should be limited to ≤ 20%. Post-implant CT scans (2.5 mm slice thickness) were obtained four weeks after implantation. CT scans were then transferred to the TPS and the structures including prostate, bladder and rectum were delineated. Although a pre-implantation plan is created to ensure 100% dose to the planning target volume, the practicalities of inserting needles into the prostate together with post-surgical oedema and seed migration can result in differences between the intended and delivered dose [23,24,25]. With this information, a dose distribution for each organ can be created for each patient; Fig. 1 is a typical example.

Fig. 1
figure 1

Typical dose distribution following 125I prostate brachytherapy. The rectum, prostate and bladder have been contoured on a CT scan 4 weeks after brachytherapy and reconstructed in 3D. Brachytherapy seeds within the prostate are highlighted green and a colourwash representing total dose has been applied to the surface of each organ

Post-implant dosimetry was evaluated for whole prostate (global) minimum dose to 90% of prostate (D90%) and the maximum dose, 0.1 cc (D0.1cc(Gy)) and 2 cc (D2cc (%)), received by the rectum and bladder. These point doses represent the highest dose delivered to the identified organ at risk and are described as by either absolute dose in Gy or as percentage of the prescribed dose and allow uniform comparison of acute and late toxicity to organ at risk dosimetry. In addition, sector analysis of 12 prostatic regions (sectors) was conducted. As described in more details in previous publications [25, 26], sector analysis involved dividing the prostate into three equal lengths along the cranio-caudal axis to create base, mid-gland and apex, then sub-dividing with a vertical plane and a horizontal plane creating right/left anterior and right/ left posterior sectors (also see Additional file 1: Figure S1). D90% was then determined for each sector of the prostate.

γH2AX and 53BP1 detection

In all the cases, peripheral venous blood was sampled in an EDTA tube and placed immediately on ice prior to transfer for analysis. For each patient, a sample was taken prior to brachytherapy to acquire a baseline and succeeding samples were gathered at 1 h, 4 h, 24 h, 4 weeks and 3 months post-seed implantation (a total of 6 blood samples for each patient). Blood samples at 4 weeks were always collected just before the acquisition of post-implant CT scans.

Separation of lymphocytes for γH2AX /53BP1 foci analysis

4 ml of whole blood was pipetted onto 3 ml of Ficoll-Plaque (GE Healthcare) and centrifuged at 4 oC, 400G, for 30 min. 1 ml of separated lymphocytes were aspirated and mixed with 10 ml of sterile phosphate buffered saline (PBS) to wash. The lymphocyte suspension was centrifuged at 4 oC, 100 g, for 20 min to form a cell pellet. The remaining PBS was decanted and the pellet re-suspended with 1 ml of fresh sterile PBS.

Fixation, staining of PBLs for γH2AX and 53BP1 foci

100 μl of re-suspended lymphocytes were loaded onto Superfrost Plus (VWR international) glass slides by cyto-spinning at 500 g for 10 min. The cells were then fixed with 4% formaldehyde at room temperature for 10 min. Each slide was then washed 3 times with PBS. The cells were permeabilized with a 0.5% (v/v) solution of Triton-X/PBS for 10 min. Blocking was achieved with 1% BSA in PBS for 1 hour at room temperature. Cells were then incubated for 1 h at room temperature with a solution of 1:5000 Mouse anti SER 139 H2AX (Millipore) and rabbit anti 53BP1 (Novus biologicals) in 0.05% Triton-X/PBS. Cells were then washed with 0.05% Triton-X/PBS. Secondary antibodies 1:2000 of goat anti-mouse Alexa-fluro 488 (Invitrogen) and goat anti rabbit 53BP1 in 0.05% Triton-X/PBS were applied for a further 1 h at room temperature. Cells were then washed with 0.05% Triton-X/PBS and air dried at room temperature. The slides were mounted with Prolong Antifade Gold with DAPI (Life Technologies) and coverslip.

Visualisation and counting of foci

Cells were visualised (by a single observer) using a Nikon fluorescent microscope at x63 magnification. Lymphocytes were identified by morphology after DAPI staining, other cell types were excluded from counting. γH2AX foci co-localising with 53BP1 foci were counted by eye in 200 cells for each time point.

Health related quality of life endpoints

Heath related QoL was assessed using the Expanded Prostate Cancer Index Composite (EPIC) instrument [7]. The EPIC questionnaire is divided into domains related to bowel, urinary, sexual and hormonal toxicities. Each domain is divided into functional and bother scores which can be combined to produce a summary score. For this study, bowel, urinary and sexual summary scores were assessed. A score of 100 indicates no problem; lower values indicate complications. For each patient, EPIC scores were obtained at baseline (pre brachytherapy), 1 month (m) post-brachytherapy, 3 m, 6 m, 9 m, 1 year (y), 2y and 3y. Primary endpoints were bowel-, urinary- and sexual- related toxicities after brachytherapy. Change-from baseline analysis was conducted using the EPIC scores. Acute toxicity was defined as those events (i.e. decrease of score from baseline) that presented within the first 3 months following treatment. Late toxicity was defined as events that developed at any time after treatment and not resolved at 2–3 years after treatment.

Statistical analysis

All correlation analysis was conducted in R (version 3.2.0) using Spearman’s correlation testing at a significance level p ≤0.05. Although multiple comparisons were conducted, no corrections were applied due to the exploratory/feasibility nature of this study (also see related section in the discussion). Moreover, to test if foci distributions corresponded to a Poisson distribution at each time-point, χ2 testing was used with a significance threshold of p = 0.05.


γH2AX and 53BP1 foci induction

Prior to implantation of 125I seeds into the prostate, mean (±SD) background foci (co-localising γH2AX/53BP1) was measured at 0.42 ± 0.20 foci per cell, ranging from 0.11 to 0.72 (Fig. 2b). A Shapiro-Wilk test confirmed the data was normally distributed (w = 0.943, p = 0.540). Following seed implantation, γH2AX/53BP1 foci numbers were significantly elevated as early as 1 h post-implantation (paired t-test p < 0.001), and remained significantly higher than background levels at 4 and 24 h, 4 weeks and 3 months post-implantation (Fig. 2b). Fig. 2c shows the mean number of foci per cell corrected for background.

Fig. 2
figure 2

a) Relative activity of 125I as a function of time post-seed-implantation. b) mean foci number per cell in patients' peripheralblood lymphocytes (the blue square on the vertical axis indicates mean background counts) and corrected for background (c). The error bars represent the standard error and the red dashed line represents the 125I half-life of 1426.3 h (59.43 days)

The highest foci numbers were observed 4 weeks (672 h) after irradiation, at 1.1 ± 0.1 foci per cell. This value is more than 2.5-fold above background foci levels (0.4 ± 0.1) and is significantly higher than 24 h values (paired t-test p = 0.02), suggesting that brachytherapy acts as a significant additional source of genetic stress in these lymphocytes. Based on ex vivo lymphocyte studies [11] showing an average of 15.5 foci induced per Gy with an average half-life of 2 h, an average of 0.7 excess foci per cell corresponds to an average dose-rate to the lymphocyte population of 4.5 cGy/h. Interestingly, this dose rate is comparable to the initial dose rate delivered to the prostate, 7.1 cGy/h (Additional file 2: Figure S2), despite the majority of lymphocytes not being within the exposed volume at any given time. This suggests either the impact of other stresses on lymphocyte DNA damage levels, or the formation of DSBs with a long half-life. This is further supported by the observation that γH2AX/53BP1 foci numbers continued to rise at 4 weeks (672 h) even after reduction in seeds activity due to natural decay (Fig. 2a), and showing only limited (statistically insignificant) decline at 3 months despite source activity falling by almost three-quarters.

As the increase in mean foci counts over time may have been due to an increased proportion of observed cells with foci or greater foci numbers in irradiated cells, the relative frequency distribution of foci per lymphocyte over time was plotted (Fig. 3). It can be seen that, at all time-points post-implantation, there is a gradual increase in both the number of hit cells, and the number of cells at each foci level, indicating some increase in damage in large number of cells rather than the presence of a small subpopulation of highly damaged cells. At no time point did more than 1% of all cells show more than 7 foci. However, while this damage may have been expected to follow a Poisson distribution, all foci distributions are significantly different to that predicted by a Poisson distribution (p < 0.0001), with an excess of undamaged cells, suggesting some degree of mixing between populations exposed to different stresses.

Fig. 3
figure 3

Foci distribution counts in patients’ peripheral blood lymphocytes pre- and post- seed implantation. The error bars represent the standard error

QoL scores over time

As demonstrated in Fig. 1, the main neighbouring organs at risk which receive high dose in prostate brachytherapy are the prostate, urethra, rectum and bladder and toxicities observed in patients undergoing brachytherapy are usually related to these organs. QoL in terms of patients’ bowel, urinary and sexual health was measured using the EPIC questionnaire before seed-implantation and up to three years post-treatment. Increased toxicities were observed post-implant peaking at 6 m for bowel symptoms (Fig. 4a), and at 1 month post-implant for both urinary (Fig. 4b) and sexual (Fig. 4c) symptoms.

Fig. 4
figure 4

Changes in quality-of-life scores over time for each domain studied on the Expanded Prostate Cancer Index Composite (EPIC) questionnaire. EPIC scores range from 0 to 100, with higher values representing a more favorable health-related quality of life. In the figure median with inter-quartile range IQR (boxes), and ± 1.5*IQR (whiskers) are presented. Dots are different data points and red dots are outliers

Global analysis

Population median prostate global D90% was (median (range)) 144.4 Gy (range: 125.4 Gy −160.3 Gy). The median of the highest point dose within the contoured rectum represented by D0.1cc was 145.2 Gy (range: 112.5 Gy – 212.5 Gy) and the highest dose delivered to a 2 cc of the rectum (D2cc) was 62.0% (range: 42.4% –84.8%). Bladder median D0.1cc was 123.0 Gy (81.3 Gy − 180.3 Gy) and D2cc =70.4% (54.0% − 91.7%).

Foci count, EPIC, sector analysis and dosimetric parameters: initial exploratory correlation analysis

Statistical correlation testing was conducted to explore potential signals of associations between the biological, patient − reported symptoms and dosimetric parameters using Spearman rank correlation coefficients. The results are presented in Figs. 5, 6 and 7. Background foci counts as well as the ratios at time point (t) relative to background (foci count prior to brachytherapy) are presented. Strong positive correlations were observed between mean foci counts at 24 h post-implantation and bowel toxicity at 1y and 2y, Fig. 5. This suggest that the larger the number of foci developed 24 h post-implant the more bowel symptoms the patients experienced. Although not statistically significant, a trend of positive correlations was also observed between the maximum dose to the rectum (rectum D0.1cc and D2cc) and late bowel toxicity. Several sector doses also correlated with early γH2AX/53BP1 foci counts (see Figs. 5, 6 and 7) but these correlations were not significant.

Fig. 5
figure 5

Bowel-related cross-correlation matrix of Spearman’s correlation coefficients (ρ). Upper triangle; ρ values with insignificant correlation (p ≤ 0.05) marked with an X. Lower triangle; colour coded coefficients. S(x) = prostatic sector number ‘x’, Rec = rectum, D0.1cc and D2cc = the doses received by the hottest 0.1 cc and 2 cc of the volume, respectively. Foci base = γH2AX/53BP1 foci count at baseline (pre-implant (t = 0 h)). Foci R t = ratio of foci count at time point t relative to baseline (0 h). EPIC(x) = change in EPIC score from baseline; positive values reflecting improvement, negative values deterioration and 0 no change. * For space consideration all ρ values are presented as percentages

Fig. 6
figure 6

Bladder-related cross-correlation matrix of Spearman’s correlation coefficients. (ρ). Upper triangle; ρ values with insignificant correlation (p ≤ 0.05) marked with an X. Lower triangle; colour coded coefficients. S(x) = prostatic sector number ‘x’, Bla = bladder, D0.1cc and D2cc = the doses received by the hottest 0.1 cc and 2 cc of the volume, respectively. Foci base = γH2AX/53BP1 foci count at baseline (pre-treatment (t = 0 h)). Foci R t = ratio of foci count at time point t relative to baseline (0 h). EPIC (x) = change in EPIC score from baseline; positive values reflecting improvement, negative values deterioration and 0 no change. * For space consideration all ρ values are presented as percentages

Fig. 7
figure 7

Sexual symptoms-related cross-correlation matrix of Spearman’s correlation coefficients (ρ). Upper triangle; ρ values with insignificant correlation (p ≤ 0.05) marked with an X. Lower triangle; colour coded coefficients. S (x) = prostatic sector number ‘x’. Foci base = γH2AX/53BP1 foci count at baseline (pre-treatment (t = 0 h)). Foci R t = ratio of foci count at time point t relative to baseline (0 h). EPIC (x) = change in EPIC score from baseline; positive values reflecting improvement, negative values deterioration and 0 no change. * For space consideration all ρ values are presented as percentages


By use of DNA DSB markers γH2AX and 53BP1, gamma-ray induced DSBs in peripheral blood lymphocytes of prostate cancer patients undergoing brachytherapy with 125I seeds can be readily visualised and quantified. X-ray induced DSBs have been detected numerous times in patients receiving EBRT [8, 9, 27, 28]. Moreover, in the field of nuclear medicine, Lassmann and colleagues have investigated the kinetics of γH2AX/53BP1 foci in peripheral leukocytes and lymphocytes in thyroid cancer patients being administered either 177Lu or 131I up to 144 hours post-infusion [12, 29, 30]. However, to our knowledge, this is the first study to demonstrate elevated DSBs in circulating lymphocytes in permanent prostate brachytherapy treatment and to study the kinetics of foci production in relation to late toxicity.

Lymphocytes within the peripheral blood circulation exist in equilibrium with large reservoirs of lymphocytes within the lungs and within extra vascular lymphatic tissue such as lymph nodes, the spleen and Peyer’s patches [31, 32]. Using ex-vivo FITC-labelled lymphocytes injected intravenously into sheep, Andrade et al. found these peripherally derived lymphocytes reached equilibrium with exchange to lymph tissue four hours post-injection [33]. Thus, at any given time, foci numbers measured in peripheral lymphocytes is a combination of both un-irradiated lymphocytes as well as those which have been irradiated by passing near the prostate at a range of earlier times. Redistribution of PBLs between intra and inter vascular spaces would be expected to reduce observed mean foci counts and increase the number of undamaged cells observed at early time points.

As expected, in the first hours post-implant an increase in foci count is observed, as levels of damage within the cells build up until the rate of repaired foci matches the induction of new damage by the radionuclide source. This process would be expected to reach equilibrium in a timeframe similar to that seen for DSB repair and PBL redistribution, both of which are typically on the order of a few hours. However, continual increases in DSB foci level are seen out to 24 hours and even four weeks post-implantation.

The continued increase to multiple weeks is particularly significant as by this time decay has reduced the activity of sources by more than a quarter. This persistent increase in damage is indicative of either extremely long-lived DNA damage within a proportion of PBL, or another effect increasing the rate of DSB induction even as the source decays. Foci kinetics at extremely long time points are poorly studied in lymphocytes, while the existence of fast and slow components of DNA repair are well established with reports of slow repair with half-lives of 38 hours [11]. It is possible to speculate that the current study is highlighting a proportion of DSBs repaired even slower or which totally failed to repair. Similarly, the effect of radiation on the distribution of lymphocytes between vascular and extra vascular − compartments is unknown. It is increasingly recognised that radiation derives an inflammatory response increasing lymphocytic infiltration in target tissues [34]. Any increase in PBL flow through the prostate would lead to a greater effective dose being seen by the lymphocytes, increasing the number of damaged cells and average foci numbers. Determining which of these effects dominates PBL foci kinetics requires further study.

One advantage DSB foci counting provides is that the dose received by the different structures (via DSBs) is over an integrated time post-implantation, whereas the post-treatment planning CT is a single snap shot of a patient’s physiology weeks later. The conventional method uses post-implant dose to predict radiation induced damage. An alternative method to assess radiation induced damage directly may lie in the ratio of foci observed between 0 h (and 4 h) and 24 h post − seed implantation.

In this study, positive correlations were found between increased early γH2AX/53BP1 foci production (24 h) and late bowel toxicity (EPIC 1y) (p = 0.035), EPIC 2 y (p = 0.001). Similar positive correlations were observed between increased foci count at 1 m relative to baseline and late bowel toxicity ((EPIC 1 y (p = 0.046), EPIC 2 y (p = 0.0189)).

Moreover, urinary toxicity positively correlated with the increase in doses received by prostatic base (sector 3, 4), D90% and the maximum dose received by the hottest 2 cc of the bladder. These observations are in agreement with results published by Hathout et al. that the dose received by the bladder neck had a strong prognostic power for both early and late toxicity [35]. It may be that the region of bladder irradiation plays an important role in predicating toxicity alongside both total dose and volume irradiated. Urethral toxicity is likely to dominate urinary toxicity, especially in brachytherapy patients [36]. For bladder and rectal doses, it’s expected that these organs will move and change position quite frequently during treatment, for example, due to bladder filling and emptying [37] and small bowel peristalsis, potentially significantly altering the prescribed dose received to either or both organs. Prostatic sector 6 dose positively correlated with increase in sexual symptoms at all times. Significant strong positive correlations were found between late and early urinary symptoms (Fig. 6). Similar correlations were also found for bowel toxicity (Fig. 5) but not for sexual symptoms (Fig. 7).

Several articles have demonstrated the association between late normal tissue toxicity after EBRT and insufficient repair of DNA double − strand breaks in breast, head and neck and prostate cancer [38, 39]. In a study investigating DNA damage response as a possible risk factor for late toxicity in prostate cancer (n = 61), van Oorschot et al. found no significant differences between patient’s γH2Ax immediately after RT [40]. However, significant differences in foci decay (at 24 h post-irradiation) were observed between radio-sensitive and radio-resistant prostate cancer patients [40]. This implies that less efficient repair of radiation − induced DSBs may contribute to the development of late toxicity. These observations have great potential in the field of personalized medicine as for many years different tests have been proposed to investigate individual patient sensitivity to RT.

Several studies have shown that γH2AX could be used for in-vivo dose estimation for EBRT [27, 28]. In contrast with EBRT, the dose delivered in brachytherapy is more confined to the prostate and it accumulates over a period of time. Moreover, in general, smaller patient’s volumes are exposed to radiation compared to EBRT. Therefore, unsurprisingly our results in this study for brachytherapy indicate a lack of correlation between foci counts at different time points and dosimetric measures as the radiobiological effect is different compared to EBRT.

As early identification of patients with high risk of toxicity is of clinical importance, we investigated initial DNA damage foci production in relation to late toxicity. Strong positive correlations were found between early foci production ratios (relative to baseline) and late rectal toxicity. If this finding is validated in future clinical studies, this assay could have clinical utility as methods to reduce rectal dose using, for example, peri-rectal hydrogel spacers for patients with high levels of circulating DSBs at early time-points after seed implant [41].

The authors acknowledge the limitations of the study with the small number of patients and the multiple comparisons conducted. There are many methods to adjust for multiple comparisons (e.g. Bonferroni correction, Ducan’s multiple range tests), however, it is overly conservative to apply these methods in the current exploratory study [42]. The results from this study can be used to help to generate hypotheses for testing in larger clinical trials. Furthermore, future investigations on the relationship between DNA damage markers γH2AX/53BP1 and other markers e.g. citrulline, ceramide, cytokines (CXCL1, CXCL6, CXCL8, CXCL10, TNFα)), at specific time points may be necessary to gain more insight into the mechanisms responsible for the observed PBL foci kinetics and the associated toxicities.


For the first time, we have demonstrated the detection and temporal changes of γH2AX/53BP1 DSB foci in peripheral blood lymphocytes for men treated with 125I permanent prostate brachytherapy. Radiation − induced γH2AX/53BP1 levels continue to increase for four weeks after seed implant and appear to reach equilibrium at less than 3 m where a clear repair signal was detected. In this small cohort of patients, levels of DSBs 24 h post-PPB seed implant positively correlated with late bowel toxicity. However, these findings require validation in larger patient cohorts to detect correlations with greater certainty.



Androgen deprivation therapy


Cubic centimetres


Double strand break


Expanded Prostate Cancer Index Composite


Permanent prostate brachytherapy


Prostate-specific antigen


Quality of life


  1. Klotz L, Vesprini D, Sethukavalan P, Jethava V, Zhang L, Jain S, et al. Long-term follow-up of a large active surveillance cohort of patients with prostate cancer. J Clin Oncol. 2015;33(3):272–7. doi:10.1200/JCO.2014.55.1192.

    Article  PubMed  Google Scholar 

  2. Bolla M, Collette L, Blank L, Warde P, Dubois JB, Mirimanoff R, et al. Long-term results with immediate androgen suppression and external irradiation in patients with locally advanced prostate cancer (an EORTC study): a phase III randomised trial. Lancet. 2002;360:103–6.

    Article  CAS  PubMed  Google Scholar 

  3. D’Amico AV, Whittington R, Malkowicz SB, Schultz D, Blank K, Broderick GA, et al. Biochemical outcome after radical prostatectomy, external beam radiation therapy, or interstitial radiation therapy for clinically localized prostate cancer. JAMA. 1998;280:969–74.

    Article  PubMed  Google Scholar 

  4. Miller DC, Sanda MG, Dunn RL, Montie JE, Pimentel H, Sandler HM, et al. Long-term outcomes among localized prostate cancer survivors: health-related quality-of-life changes after radical prostatectomy, external radiation, and brachytherapy. J Clin Oncol. 2005;23:2772–80.

    Article  PubMed  Google Scholar 

  5. Potters L, Morgenstern C, Calugaru E, Fearn P, Jassal A, Presser J, et al. 12-year outcomes following permanent prostate brachytherapy in patients with clinically localized prostate cancer. J Urol. 2005;173:1562–6.

    Article  PubMed  Google Scholar 

  6. Sanda MG, Dunn RL, Michalski J, Sandler HM, Northouse L, Hembroff L, et al. Quality of life and satisfaction with outcome among prostate-cancer survivors. N Engl J Med. 2008;358(12):1250–61. doi:10.1056/NEJMoa074311.

    Article  CAS  PubMed  Google Scholar 

  7. Wei JT, Dunn RL, Litwin MS, Sandler HM, Sanda MG. Development and validation of the expanded prostate cancer index composite (EPIC) for comprehensive assessment of health-related quality of life in men with prostate cancer. Urol. 2000;56:899–905.

    Article  CAS  PubMed  Google Scholar 

  8. Zahnreich S, Ebersberger A, Kaina B, Schmidberger H. Biodosimetry Based on γ-H2AX Quantification and Cytogenetics after Partial- and Total-Body Irradiation during Fractionated Radiotherapy. Radiat Res. 2015;183:432–46.

    Article  CAS  PubMed  Google Scholar 

  9. Woolf DK, Williams NR, Bakshi R, Madani SY, Eaton DJ, Fawcitt S, et al. Biological dosimetry for breast cancer radiotherapy: a comparison of external beam and intraoperative radiotherapy. Springerplus. 2014;3:329.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Rogakou EP, Pilch DR, Orr AH, Ivanova VS, Bonner WM. DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139. J Biol Chem. 1998;273:5858–68.

    Article  CAS  PubMed  Google Scholar 

  11. Horn S, Barnard S, Rothkamm K. (2011). Gamma-H2AX-Based Dose Estimation for Whole and Partial Body Radiation Exposure. Plos One. 6(9):e25113.

  12. Lassmann M, Hänscheid H, Gassen D, Biko J, Meineke V, Reiners C, et al. In-vivo formation of gamma-H2AX and 53BP1 DNA repair foci in blood cells after radioiodine therapy of differentiated thyroid cancer. J Nucl Med. 2010;51:1318–25.

    Article  PubMed  Google Scholar 

  13. Horn S, Barnard S, Brady D, Prise KM, Rothkamm K. Combined analysis of gamma-H2AX/53BP1 foci and caspase activation in lymphocyte subsets detects recent and more remote radiation exposures. Radiat Res. 2013;180(6):603–9. doi:10.1667/RR13342.1.

    Article  CAS  PubMed  Google Scholar 

  14. Grudzenski S, Raths A, Conrad S, Rübe CE, Löbrich M. Inducible response required for repair of low-dose radiation damage in human fibroblasts. Proc Natl Acad Sci U S A. 2010;10(32):14205–10. doi:10.1073/pnas.1002213107. PubMed PMID: 20660770; PubMed Central PMCID: PMC2922519.

    Article  Google Scholar 

  15. Martin OA, Ivashkevich A, Choo S, Woodbine L, Jeggo PA, Martin RF, Lobachevsky P. Statistical analysis of kinetics, distribution and co-localisation of DNA repair foci in irradiated cells: cell cycle effect and implications for prediction of radiosensitivity. DNA Repair (Amst). 2013;12(10):844–55. doi:10.1016/j.dnarep.2013.07.002.

    Article  CAS  Google Scholar 

  16. Rothkamm K, Balroop S, Shekhdar J, Fernie P, Goh V. Leukocyte DNA damage after multi-detector row CT: a quantitative biomarker of low-level radiation exposure. Radiology. 2007;242:244–51.

    Article  PubMed  Google Scholar 

  17. Melian E, Mageras GS, Fuks Z, Leibel SA, Niehaus A, Lorant H, et al. Variation in prostate position quantitation and implications for three-dimensional conformal treatment planning. Int J Radiat Oncol Biol Phys. 1997;38:73–81.

    Article  CAS  PubMed  Google Scholar 

  18. Pinkawa M, Asadpour B, Gagel B, Piroth MD, Holy R, Eble MJ. Prostate position variability and dose-volume histograms in radiotherapy for prostate cancer with full and empty bladder. Int J Radiat Oncol Biol Phys. 2006;64:856–61.

    Article  PubMed  Google Scholar 

  19. Nath R, Bice WS, Butler WM, Chen Z, Meigooni AS, Narayana V, et al. American Association of Physicists in Medicine. AAPM recommendations on dose prescription and reporting methods for permanent interstitial brachytherapy for prostate cancer: report of Task Group 137. Med Phys. 2009;36(11):5310–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. The Royal College of Radiologists. Quality Assurance Practice Guidelines for Transperineal LDR Permanent Seed Brachytherapy for Prostate Cancer. London: The Royal College of Radiologists; 2012.

  21. Rivard MJ, Coursey BM, DeWerd LA, Hanson WF, Huq MS, Ibbott GS, et al. Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys. 2004;31(3):633–74. Erratum in: Med Phys. 2004;31(12):3532–3.

    Article  PubMed  Google Scholar 

  22. Salembier C, Lavagnini P, Nickers P, Mangili P, Rijnders A, Polo A, et al. GEC ESTRO PROBATE Group. Tumour and target volumes in permanent prostate brachytherapy: a supplement to the ESTRO/EAU/EORTC recommendations on prostate brachytherapy. Radiother Oncol. 2007;83(1):3–10. Epub 2007 Feb 26.

    Article  PubMed  Google Scholar 

  23. Tejwani A, Bieniek E, Puckett L, Lavaf A, Guirguis A, Bennish A, et al. Case series analysis of post-brachytherapy prostate edema and its relevance to post-implant dosimetry. Post-implant prostate edema and dosimetry. J Contemp Brachytherapy. 2012;4:75–80.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Roberson PL, Narayana V, McShan DL, Winfield RJ, McLaughlin PW. Source placement error for permanent implant of the prostate. Med Phys. 1997;24(2):251–7.

    Article  CAS  PubMed  Google Scholar 

  25. Mohamed Yoosuf AB, Mitchell DM, Workman G, Jonnada S, Napier E, Jain S. Sector analysis provides additional spatial information on the permanent prostate brachytherapy learning curve. Brachytherapy. 2015;14(5):703–10. doi:10.1016/j.brachy.2015.05.009. Epub 2015 Jun 27.

    Article  PubMed  Google Scholar 

  26. Mohamed Yoosuf AB, Workman G, O’Toole MM, Straney M, Verghis R, Napier E, et al. Sector analysis of 125I permanent prostate brachytherapy provides a rapid and effective method of evaluating and comparing pre- and post-implant dosimetry. Brachytherapy. 2013;12:254–9.

    Article  PubMed  Google Scholar 

  27. Zwicker F, Swartman B, Sterzing F, Major G, Weber KJ, Huber PE, et al. Biological in-vivo measurement of dose distribution in patients’ lymphocytes by gamma-H2AX immunofluorescence staining: 3D conformal- vs. step-and-shoot IMRT of the prostate gland. Radiat Oncol. 2011;6:62. doi:10.1186/1748-717X-6-62.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zwicker F, Swartman B, Roeder F, Sterzing F, Hauswald H, Thieke C, et al. In vivo measurement of dose distribution in patients’ lymphocytes: helical tomotherapy versus step-and-shoot IMRT in prostate cancer. J Radiat Res. 2015;56(2):239–47. doi:10.1093/jrr/rru096. Epub 2014 Oct 31.

    Article  CAS  PubMed  Google Scholar 

  29. Eberlein U, Nowak C, Bluemel C, Buck AK, Werner RA, Scherthan H, et al. DNA damage in blood lymphocytes in patients after (177)Lu peptide receptor radionuclide therapy. Eur J Nucl Med Mol Imaging. 2015;42(11):1739–49. doi:10.1007/s00259-015-3083-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Eberlein U, Scherthan H, Bluemel C, Peper M, Lapa C, Buck AK, et al. DNA Damage in Peripheral Blood Lymphocytes of Thyroid Cancer Patients After Radioiodine Therapy. J Nucl Med. 2016;57(2):173–9. doi:10.2967/jnumed.115.164814.

    Article  CAS  PubMed  Google Scholar 

  31. Pabst R, Binns RM, Licence ST, Peter M. Evidence of a selective major vascular marginal pool of lymphocytes in the lung. Am Rev Respir Dis. 1987;136:1213–8.

    Article  CAS  PubMed  Google Scholar 

  32. Gowans JL. The lymphocyte—a disgraceful gap in medical knowledge. Immunol Today. 1996;17:288–91.

    Article  CAS  PubMed  Google Scholar 

  33. Andrade WN, Johnston MG, Hay JB. The relationship of blood lymphocytes to the recirculating lymphocyte pool. Blood. 1998;91:1653–61.

    CAS  PubMed  Google Scholar 

  34. Multhoff G, Radons J. Radiation, inflammation, and immune responses in cancer.Front. Oncol. 2012;2:58. doi:10.3389/fonc.2012.00058.

    Google Scholar 

  35. Hathout L, Folkert MR, Kollmeier MA, Yamada Y, Cohen GN, Zelefsky MJ. Dose to the bladder neck is the most important predictor for acute and late toxicity after low-dose-rate prostate brachytherapy: implications for establishing new dose constraints for treatment planning. Int J Radiat Oncol Biol Phys. 2014;90:312–9.

    Article  PubMed  Google Scholar 

  36. Thomas C, Keyes M, Liu M, Moravan V. Segmental urethral dosimetry and urinary toxicity in patients with no urinary symptoms before permanent prostate brachytherapy. Int J Radiat Oncol Biol Phys. 2008;72:447–55.

    Article  PubMed  Google Scholar 

  37. Stock RG, Stone NN. Importance of post-implant dosimetry in permanent prostate brachytherapy. Eur Urol. 2002;41:434–9.

    Article  PubMed  Google Scholar 

  38. Chua ML, Somaiah N, A’Hern R, Davies S, Gothard L, Yarnold J, et al. Residual DNA and chromosomal damage in ex vivo irradiated blood lymphocytes correlated with late normal tissue response to breast radiotherapy. Radiother Oncol. 2011;99(3):362–6. doi:10.1016/j.radonc.2011.05.071.

    Article  CAS  PubMed  Google Scholar 

  39. Goutham HV, Mumbrekar KD, Vadhiraja BM, Fernandes DJ, Sharan K, Kanive Parashiva G, et al. DNA double-strand break analysis by γ-H2AX foci: a useful method for determining the overreactors to radiation-induced acute reactions among head-and-neck cancer patients. Int J Radiat Oncol Biol Phys. 2012;84(5):e607–12. doi:10.1016/j.ijrobp.2012.06.041.

    Article  CAS  PubMed  Google Scholar 

  40. van Oorschot B, Hovingh SE, Moerland PD, Medema JP, Stalpers LJ, Vrieling H, et al. Reduced activity of double-strand break repair genes in prostate cancer patients with late normal tissue radiation toxicity. Int J Radiat Oncol Biol Phys. 2014;88(3):664–70. doi:10.1016/j.ijrobp.2013.11.219.

    Article  PubMed  Google Scholar 

  41. Vanneste BGL, Pijls-Johannesma M, van De Voorde L, van Lin EN, van de Beek K, van Loon J, et al. Spacers in radiotherapy treatment of prostate cancer: Is reduction of toxicity cost-effective? Radiother Oncol. 2015;114(2):276–81.

    Article  PubMed  Google Scholar 

  42. Bender R, Lange S. Adjusting for multiple testing--when and how? J Clin Epidemiol. 2001;54(4):343–9. Review.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors are grateful to the Friends of the Cancer Centre, the Belfast-Manchester Movember Centre of Excellence (CE013-2-004), funded in partnership with Prostate Cancer UK and the Northern Ireland Health and Social Care Trust R & D Office (COM/4965/14) for funding this work.


This work was supported by grants from Prostate Cancer UK, the Belfast-Manchester Movember Centre of Excellence (CE013-2-004), and the Northern Ireland Health and Social Care Trust R & D Office grant (COM/4965/14).

Availability of data and material

Not applicable.

Authors’ contributions

SO collated the data, conducted statistical analysis and drafted the final version of this manuscript. SH contributed in drafting the manuscript and in conducting γH2AX/53BP1 phosphorylation detection and quantification on blood samples. DB initiated blood samples collection, ran γH2AX/53BP1 phosphorylation detection and quantification and conducted preliminary analysis on the data. SMcM critically read the manuscript and provided valuable suggestions. AY was responsible for treatment planning and post-implant dosimetry. SJ and DM are the treating oncologists and were responsible for recruiting patients for this study. KC handled the EPIC forms and converted to digital format. CL designed and maintained the clinical database. AH and CMcG provided physics supervision. JO’S, KP and SJ designed this study and also supervised DB and SH. All authors read and approved the final version of this manuscript.

Competing interest

None to declare.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Approval for performing this study was obtained from ORECNI (Office for Research Ethics Committees Northern Ireland. Reference number: 11/NI/0149).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Sarah O. S. Osman.

Additional files

Additional file 1: Figure S1.

Twelve-sector analysis created by dividing the prostate into three equal lengths along the cranio-caudal axis to form the base, mid-gland, and apex and then subdivide with vertical and horizontal plans generating right/left posterior sectors (reused from (21) with permission from AB Mohamed Yoosuf). (PDF 61 kb)

Additional file 2: Figure S2.

Top; dose rate (cGy/h) as a function of time since implant for 125I monotherapy plotted using the equation, \( Dose\ rate\ (t)= m P D\left(\frac{ \ln 2}{t_{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}\right){.2}^{\left(\frac{t}{t_{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}\right)} \), where t is the elapsed time, mPD is the minimum peripheral dose (=145 Gy for 125I) and \( {t}_{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.} \) is the half-life (=59.43 days for 125I). Bottom; the time required to deliver relative fraction of the prescribed dose, \( Fractional\ dose(t)=1-{e}^{-\left(\frac{t. \ln 2}{t_{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}\right)} \). Reference: Dale RG. The applications of the linear-quadratic dose effect equation to fractionated and protracted therapy. Br J Radiol 1985; 58: 515–28. (PDF 178 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Osman, S.O.S., Horn, S., Brady, D. et al. Prostate cancer treated with brachytherapy; an exploratory study of dose-dependent biomarkers and quality of life. Radiat Oncol 12, 53 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: