| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Full paper |
1 Department of Radiation Medicine, The Ohio State University, Columbus, OH 43210, 2 Department of Radiation Oncology, Medical College of Wisconsin, Milwaukee, WI 53226, USA
Correspondence: Jian Z Wang, Department of Radiation Medicine, James Cancer Hospital and Solove Research Institute, The Ohio State University, 300 W 10th Ave, Rm 094, Columbus, OH 43210, USA. E-mail: wang.993{at}osu.edu
| Abstract |
|---|
|
|
|---|
= 0.15 Gy 1,
/
= 3.1 Gy and the number of clonogens K = 106
107 cells) determined in earlier studies were used to derive the COER for prostate cancer: COER = 1.4 with a standard confidence interval (CI) of (1.2, 1.8). The result is consistent with the in vitro OER measurements of human tumour cell lines under chronic hypoxia conditions. This implies that a higher dose is needed to overcome tumour hypoxia. For prostate tumours, the prescription dose required to overcome tumour hypoxia is 165 Gy (CI: 153
186 Gy) for permanent 125I implants and 88 Gy (CI: 74
118 Gy) in 2 Gy fractions for external-beam radiotherapy. The impact of LQ parameters on the calculations of COER and dose escalation was discussed. This study provides a preliminary estimate of the dose escalation needed to overcome tumour hypoxia based on clinical data. More clinical data with better statistics and longer follow-up time are required to further tune the radiobiological modelling of hypoxia for prostate cancer. | Introduction |
|---|
|
|
|---|
Recently, Movsas et al [13] and Parker et al [4] have demonstrated that hypoxic regions do exist in human prostate carcinoma and the treatment outcome of radiation therapy correlates with the degree of hypoxia. Oxygenation status (oxygen pressure, PO2) may thus be another underlying biological parameter beyond the classic prognostic factors (clinical stage, Gleason score and prostate specific antigen (PSA)) predicting a critical component of treatment failure in prostate cancer. Movsas et al [1] reported a significant difference in the biochemical control rate at 2 years (31% versus 92%) between two patient groups stratified by the prostate/muscle (P/M) PO2 ratio (<0.05 versus
0.05). Many other investigators have also found similar correlations between hypoxia and tumour progression [510]. It is generally believed that high-risk tumours tend to exhibit more severe hypoxia [3, 10]. Although it is well established that hypoxic cells/regions are less sensitive to radiation and that tumour hypoxia can result in clinical radiation therapy failure, it remains unclear how to overcome hypoxia in the clinical setting. Specific quantitative dose guidelines are needed that can be applied to the clinical setting, particularly to dose escalation trials aimed to combat hypoxia in prostate cancer.
The linear-quadratic (LQ) and tumour control probability (TCP) models extended to address tumour hypoxia were used to analyse the clinical data. Recently, the
/
ratio of the LQ model for prostate cancer has become a highly debated topic in the radiation therapy community [1122]. By taking into account the effect of tumour repopulation, Wang et al [14, 15] and Kal and Van Gellekom [19] took two different approaches, and obtained similar results of
/
ratio (around 3.1 Gy). Based on the clinical data of external-beam radiotherapy (EBRT) and low- and high-dose rate (LDR/HDR) brachytherapy [12, 13, 23], Wang et al [14, 15] have derived a new set of LQ parameters (
= 0.15 Gy1,
/
= 3.1 Gy and repair half-time Tr = 16 min) with estimated numbers of clonogens around 106
107, depending on the patient risk level. The new results provide reasonable estimates of radiosensitivity and the number of clonogens for human prostate tumours [14, 15, 2426]. The radiobiological modelling with these parameters provides a consistent interpretation for most currently available clinical data, including EBRT, permanent brachytherapy, HDR brachytherapy and their combinations [14, 15, 26].
Fowler et al [16, 21] argued that for slowly growing tumours (e.g. prostate cancer), the onset of tumour repopulation may be delayed significantly to over 200 days. Wang et al [17, 22] pointed out that there is no strong biological evidence for such arguments. The clinical results reported so far have been controversial. In a recently reported study, Perez et al [27] did observe an impact of elapsed treatment time on outcome of EBRT for T2 prostate tumours, implying the onset time of repopulation should be relatively short. Therefore in our study, the original assumption of the onset time of tumour repopulation will be followed for various sets of LQ parameters. Recently, the effect of hypoxia has been taken into account by Nahum et al [18] in their modelling study. However, their model parameters (
,
and oxygen enhancement ratio (OER)) were completely based on in vitro measurements, which may not be fully applicable to the clinical micro-environment of prostate cancer in vivo [28]. Tumour heterogeneity of radiosensitivity has been studied [20] and due to the large uncertainty existing in the current clinical data and the large number of parameters used, it is not possible to give a reasonable estimate of
/
ratio for any statistical significance using the heterogeneous model. For this reason, we chose the homogeneous model in this study.
The purpose of our study is to quantify the dose escalation to address tumour hypoxia in prostate cancer. The LQ and TCP models are used to interpret the clinical data recently reported from Fox Chase [1] based on the concept of OER. The specific aim of this analysis is to derive the clinical OER (COER) and to determine the dose escalation necessary to offset the hypoxia effect using either permanent brachytherapy or EBRT.
| Materials and methods |
|---|
|
|
|---|
|
|
and
characterize intrinsic radiosensitivity, G is the dose-rate factor accounting for sub-lethal damage repair,
is the effective tumour repopulation rate [
= ln(2)/Td, and Td is the effective clonogen doubling time]. A median potential tumour doubling time of 42 days measured in vivo from 7 prostate cancer patients was used [30]. This repopulation rate corresponds to a biological effective dose (BED) of 0.11 Gy per day (for
= 0.15 Gy 1), compatible with the results of slow-growing tumours presented in a repopulation-dose-equivalent study by Jones et al [31]. The repopulation time T can be replaced by (TTk), where Tk is the onset time of tumour repopulation. Because the actual value of this parameter is under debate [16, 17, 21, 22, 27] and no clinical data are available, we followed the original assumptions used by different investigators when they derived the LQ parameters from clinical data, i.e. Tk = 0 for the LQ parameters by Wang et al [14, 15] and by Kal and Van Gellekom [19], and Tk>200 days (i.e. ignoring the repopulation effect) for the LQ parameters by Brenner et al [11, 12] and by Fowler et al [13].
For EBRT, the total dose D equals nd and the dose-rate factor G equals 1/n, where n is the number of dose fractions and d is the dose per fraction. The treatment time T of EBRT can be approximated as the number of treatment fractions multiplied by 1.4 (5 fractions per week). For permanent brachytherapy, D, G and T have been derived in many papers (e.g. [14, 29]) as follows:
|
|
, D0 is the prescribed dose) and
is the decay constant for the implanted isotopes [
= ln(2)/Ts, Ts is the half-life time of the radioisotope]. Ts = 59.4 days for 125I implants.
The Poisson TCP model is often used to link cell killing of a radiotherapy scheme to treatment outcome [32],
|
|
Two representative sets of LQ parameters derived from clinical data [1115, 19] and one set of generic LQ parameters for tumour [10] were used in this study (Table 1
). The numbers shown in the parentheses of Table 1
indicate the standard confidence intervals (CI) [14] and they were used to estimate the uncertainties of the COER and dose escalations. The clonogen number K of 1.6x106 (CI: 5.6x104
8.8x107) derived from clinical data [14] was used in the calculations with the first parameter set (see the following subsection "Fox Chase clinical data"). Due to limited data for the other two parameter sets, the clonogen numbers were arbitrarily set to match the clinical data and the CIs were not considered. Our calculation was mainly based on the first parameter set. The later two parameter sets were used to test the model dependency of the final results presented in this study.
|
Given a certain amount of clonogens, to achieve the same TCP for two patient groups, the cell killing efficiency should be the same. That is:
|
|
h =
o/H,
h =
o/H2, and they are valid to apply to the LDR 125I brachytherapy and the EBRT with low dose fractions (d = 2 Gy). It appears that the hypoxia impacts more significantly on the quadratic term than on the linear term. Such relationships of the LQ parameters between aerated and hypoxic cells have been observed in in vitro experiments of low dose irradiation (0.53 Gy) for Chinese hamster V79-171 cells [37, 38]. Therefore, the LQ model extended to the hypoxia effect can be expressed as (Sh for hypoxic cells):
|
|
Because of the limited clinical data, we will focus on the effect of hypoxia on radiosensitivity only, and neglect the influence of hypoxia on repair time or repopulation rate. Hill et al presented data of animal and human carcinoma of the cervix showing no evidence for a difference in repopulation kinetics between aerated or hypoxic tumours [39, 40]. It is generally assumed that the transient hypoxia can be overcome by tumour reoxygenation by protracting radiotherapy (i.e. fractionated EBRT or LDR brachytherapy). The reoxygenation effect has been modelled in previous studies (e.g. [41]). Therefore, it is believed that the oxygen effect on the clinical outcome observed in the 125I data (see next sub-section) should be mainly due to the chronic hypoxia [36]. In this work, we use the above formulae to study the effect of chronic hypoxia.
For patients with tumours that are only partially hypoxic, i.e. containing oxygenated and hypoxic subvolumes, the two-compartment model was used to calculate the overall surviving fraction [18, 33, 35],
|
|
The uncertainty in estimating the COER originates from the uncertainties of the LQ/TCP parameters (as shown in Table 1
), the hypoxia fraction and the clinical data. The LQ/TCP parameters are strongly correlated to each other and constrained by clinical data; therefore, their uncertainties are not independent from each other [14] and have been considered in this study. The code autoEUD described in Wang and Li [26] was used to propagate the different uncertainties to determine the CIs of COER and dose escalation. For independent parameters, their contributions to the final CIs were calculated based on the root of the quadratic-summation.
For a given patient group of prostate cancer, based on Equations (1) to (6), we could calculate the TCP for a given radiation treatment or determine the required radiation dose D to achieve a given TCP. The software AutoEUD and the typical dosevolume histogram (DVH) of EBRT and brachytherapy described in a previous paper [26] were used to automate these calculations.
Fox Chase clinical data
Since 1999, Movsas et al have published a series of papers to address the hypoxia of prostate cancer [13]. For the first time, they demonstrated the existence of hypoxia in human prostate cancer based on the in vivo electrode measurements of oxygen levels, and that the degree of hypoxia in prostate cancer correlated to the treatment outcome of radiotherapy. 57 patients with localized prostate cancer were included in their study [1]. Before radiotherapy treatments, custom-made Eppendorf microelectrodes were used to measure the PO2 in both pathologically involved regions of the prostate and normal muscle (as an internal control). Real-time ultrasound imaging was used to guide the microelectrode measurement. For each patient, approximately 100 PO2 readings were obtained along 3 to 5 tracks in both prostate and muscle [3]. Following the Eppendorf PO2 measurement, 48 patients received permanent 125I brachytherapy. The dose prescribed to the prostate was 145 Gy. Nine patients were treated on a dose escalation protocol involving a 46 Gy (in 23 daily fractions) of EBRT to the pelvic region plus two HDR boost implants (8.75 Gy or 9.75 Gy per fraction). Hormonal therapy prior to radiation therapy was given to nine patients. In their study, the biological failure was defined as two consecutive rises in PSA without a return to baseline. With a median follow-up time of 19 months (ranging from 4 months to 31 months), the overall 2 year bNED (biologically no evidence of disease) was estimated as 81% for the entire patient group (as shown in Figure 1
of this paper and Figure 2
of [1]). In order to study the impact of hypoxia, the P/M ratio of PO2, which eliminates the potential interpatient and technical variations, was used to analyse the treatment outcome.
|
|
0.05 (aerated group): 31% versus 92% (p<0.0001) (see Figure 1
According to the patient characteristics of the three prognosis factors, patients of the Fox Chase data were stratified into the low- or intermediate-risk groups. Because most patients presented low-risk prognosis factors, the clonogen number (1.6x106) of low-risk patients was used for this study [14, 15]. Two radiotherapy modalities (125I brachytherapy vs combined EBRT and HDR) were used in the Fox Chase Study. Since the patient number of the combined modality (EBRT + HDR) is relatively small (9 patients) and the treatment outcome for low-risk patients is close for the two regimens, we ignored the difference of the two modalities and combined the data in the following study. For similar reasons, we also ignored the influence of hormone therapy applied to nine patients out of the total 57 patients. A representative DVH obtained from permanent 125I implant (Figure 1
of Ref. [26]) was used to account for the dose inhomogeneity of the Fox Chase data.
| Results |
|---|
|
|
|---|
1.8) (Table 2
|
= 0.04 Gy 1,
/
= 1.5 Gy, Tr = 1.9 h with repopulation ignored), only the hypoxic clonogens were considered [43]. The COER obtained with this parameter set is H = 1.5. Using the LQ parameters that are considered generic for tumours (i.e.
= 0.3 Gy 1,
/
= 10 Gy, Tr = 1 h and Td = 42 days), we obtained the COER: H = 1.2. Interestingly, the result derived with the parameters of Wang et al lies between the values obtained with these two sets of LQ parameters, and the differences among the COER values obtained with the three parameter sets are smaller than the CIs determined by various uncertainties. More clinical data with better statistics and better knowledge of the hypoxic fraction are required to further narrow down the CIs. To summarize the results, the COERs, as well as the required dose escalations (see the following subsection), derived with different LQ/TCP parameters are listed in Table 2
Dose escalation for hypoxic tumours
Based on the derived COER, the dose required to compensate for tumour hypoxia was calculated using the code autoEUD [26]. To bring the TCP for hypoxic tumours to the same level as that for the general population of prostate patients (TCP = 81%), the prescription dose required for permanent 125I implant should be 169 Gy. This dose may be used to target the hypoxia region if it can be identified. Otherwise, a reduced dose of 165 Gy (CI: 153
186 Gy) should be prescribed to the entire prostate, including the non-hypoxic regions, to achieve the same treatment outcome (a TCP of 81%) as shown in Figure 1
. This dose escalation is practically feasible, as demonstrated by a recent study [44].
The calculations performed so far provide a guide of dose escalation for a cohort of patients with hypoxic prostate tumours. If the hypoxic subvolume of individual tumours can be detected, for example by biological/functional imaging, the data presented may also be used to target these subvolumes and design the individualized radiation plan. Figure 2
shows the TCP and dose prescription as a function of the hypoxic fraction fh for individual prostate tumours. The TCP presented in Figure 2a
is for the dose prescription (145 Gy of 125I implant) used in the current clinical practice. Figure 2b
shows the dose escalation required to overcome the effect of hypoxia and to bring the TCP to 81% as that of the general population of prostate cancer patients. For individual patients with hypoxic fraction changing from 5% to 50%, the prescription dose is increased from 157 Gy to 181 Gy.
The above results obtained from the Fox Chase data (brachytherapy) data set may be used to derive the dose escalation of EBRT. Various clinical studies have shown that a biological equivalence exists between EBRT and permanent brachytherapy for prostate cancer [13, 26]. For patients with similar risk-levels, the 145 Gy brachytherapy yielded a similar treatment outcome as 71 Gy EBRT in 2 Gy fractions [13]. Based on the LQ and TCP models and the COER derived in this study, we calculated the EBRT dose required to compensate for the effect of hypoxia. This EBRT dose is found to be 88 Gy (CI: 74
118 Gy) in 2 Gy daily fractions.
Similarly, based on the COERs derived with the parameters of Brenner et al or with the generic LQ parameters for tumours (see Table 1
), we obtained the dose escalation necessary to compensate for prostate hypoxia. With Brenner et al's parameters, the prescription dose would be 193 Gy for 125I brachytherapy and 109 Gy for EBRT. With the generic tumour parameters, these doses would be 155 Gy for the 125I brachytherapy and 102 Gy for EBRT. Because the clinical data used to derive the COERs are from 125I brachytherapy, the dose escalations calculated for the 125I brachytherapy have smaller CIs and they match the sequence of the
/
ratio in the three parameter-sets. However, the extrapolated EBRT doses show large uncertainty and the results based on the parameters of Brenner et al and the generic tumour parameters are unusually high; therefore, the prediction of EBRT dose escalation has limited clinical value.
| Conclusion and discussion |
|---|
|
|
|---|
1.8). The result agrees with the in vitro OER measurements under chronic/continued hypoxia conditions [36]. To overcome tumour hypoxia, a dose escalation to 165 Gy (CI: 153
186 Gy) (instead of 145 Gy) is required for 125I permanent implants. A similar calculation was performed for EBRT. The EBRT dose needs to be escalated from 71 Gy to 88 Gy in 2 Gy fractions in order to achieve a TCP of 81% for the hypoxic tumour group. There is in vivo evidence showing that increasing levels of hypoxia in prostate cancer correlate significantly with the increasing clinical stage and patient age [2]. The P/M ratio of PO2 may be used as a prognostic factor for high-risk patients. Dose escalation in both EBRT and brachytherapy has been proven effective to treat high-risk patients of prostate cancer [12, 23, 26, 44, 45]. This study used similar approaches to address the unsatisfactory treatment outcome of hypoxic tumours. The required dose escalations (165 Gy for 125I implant and 88 Gy for EBRT) are practically deliverable [23, 44], and clinical trials with similar dose levels or even higher are being conducted in several institutions.
The results obtained in this work are based on the Fox Chase clinical data [1] that, for the first time, shows the existence of hypoxia in human prostate tumours and the correlation of the degree of hypoxia with treatment outcome. There are several limitations in the Fox Chase study that may compromise the accuracy in the current estimation of the COER and the corresponding target doses. First, the PO2 measurement with Eppendorf microeletrodes is a sampling method which measures only a sampled portion of the prostate. It may not reflect the overall PO2 distribution in the entire tumour. The method also fails to discriminate cell type and viability, resulting in PO2 readings from less significant tissues. Second, this method cannot distinguish the hypoxia types: chronic and transient. The transient hypoxia may be reduced or overcome by providing enough time for reoxygenation during the treatment course. Third, the study contained only a small number of patient samples with limited follow-up time, and furthermore mixed treatment modalities (two radiation methods and hormone therapy) were involved in the small data sample. Because of the limited statistical power, only two patient groups could be identified with an artificially determined threshold of the P/M ratio. More clinical data with larger numbers of patients, longer follow-up time and consistent treatment are certainly desired to further refine our conclusions.
Because of the limited clinical data, it was not possible to derive independent relationships between
,
and COER in this study. Instead, the relationships,
h =
o/H,
h =
o/H2, from the literature [5, 3335] were used. Although similar relationships have been observed in in vitro measurements [37, 38], the results obtained in this analysis would be affected if such relationships were not applicable to prostate cancer. The COER calculation was based on an average hypoxic fraction of 15% [5, 10, 42]. A change of this value would change the results of the COER as well as the dose escalation. There are also concerns about the OER changing with the irradiation dose fraction. In this study, we assume that the COER keeps constant for the LDR brachytherapy and the low-dose-fraction (2 Gy) EBRT. The use of alternative models, which include temporal variations in the hypoxic fraction, might also lead to quite different conclusions. Therefore, the results obtained in this study may be limited by the approximations used in the models and the uncertainties shown in the clinical data. Caution needs to be exercised in using the presented data for clinical decision-making purposes.
Received for publication January 4, 2006. Revision received June 2, 2006. Accepted for publication June 5, 2006.
| References |
|---|
|
|
|---|
/
ratio), similar to late-responding normal tissue. Int J Radiat Oncol Biol Phys 2002;52:613.[CrossRef][Medline]
/
for prostate cancer really low? Int J Radiat Oncol Biol Phys 2001;50:102131.[CrossRef][Medline]
/
ratio for prostate cancer? Int J Radiat Oncol Biol Phys 2003;55:194203.[CrossRef][Medline]
/
ratio for prostate cancer: what does the clinical outcome of HDR brachytherapy tell us? Int J Radiat Oncol Biol Phys 2003;57:11018.[CrossRef][Medline]
/
ratio for prostate cancer? In regard to Wang et al. [letter]. Int J Radiat Oncol Biol Phys 2003;57:5935.[Medline]
/
ratio for prostate cancer: in response to Dr. Fowler et al. [letter]. Int J Radiat Oncol Biol Phys 2003;57:5956.[CrossRef]
/
ratio. Int J Radiat Oncol Biol Phys 2003;57:391401.[CrossRef][Medline]
/
ratio for prostate cancer? Int J Radiat Oncol Biol Phys 2003;57:111621.[CrossRef][Medline]
/
estimates from homogeneous (individual) and heterogeneous (population) tumor control models for early stage prostate cancer. Med Phys 2003;30:283248.[CrossRef][Medline]
/
ratios for prostate cancer'. Phys Med Biol 2005;50:L14.[CrossRef][Medline]
/
ratios for prostate cancer'". Phys Med Biol 2005;50:L58.[CrossRef]
/
ratio of 1.5 from Brenner and Hall a modeling artifact? [letter]. Int J Radiat Oncol Biol Phys 2000;47:5368.[CrossRef][Medline]
/
values for prostate appear to be independent of modeling details [letter]. Int J Radiat Oncol Biol Phys 2000;47:275365.This article has been cited by other articles:
![]() |
BJR review of the year - 2006 Br. J. Radiol., March 1, 2007; 80(951): 147 - 151. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| BJR | DMFR | IMAGING | ALL BIR JOURNALS |