Skip to main content
Journal of Radiation Research logoLink to Journal of Radiation Research
. 2014 Apr 11;55(5):902–911. doi: 10.1093/jrr/rru020

Clinical oxygen enhancement ratio of tumors in carbon ion radiotherapy: the influence of local oxygenation changes

Laura Antonovic 1,*, Emely Lindblom 1, Alexandru Dasu 2, Niels Bassler 3, Yoshiya Furusawa 4, Iuliana Toma-Dasu 1,5
PMCID: PMC4240637  PMID: 24728013

Abstract

The effect of carbon ion radiotherapy on hypoxic tumors has recently been questioned because of low linear energy transfer (LET) values in the spread-out Bragg peak (SOBP). The aim of this study was to investigate the role of hypoxia and local oxygenation changes (LOCs) in fractionated carbon ion radiotherapy. Three-dimensional tumors with hypoxic subvolumes were simulated assuming interfraction LOCs. Different fractionations were applied using a clinically relevant treatment plan with a known LET distribution. The surviving fraction was calculated, taking oxygen tension, dose and LET into account, using the repairable–conditionally repairable (RCR) damage model with parameters for human salivary gland tumor cells. The clinical oxygen enhancement ratio (OER) was defined as the ratio of doses required for a tumor control probability of 50% for hypoxic and well-oxygenated tumors. The resulting OER was well above unity for all fractionations. For the hypoxic tumor, the tumor control probability was considerably higher if LOCs were assumed, rather than static oxygenation. The beneficial effect of LOCs increased with the number of fractions. However, for very low fraction doses, the improvement related to LOCs did not compensate for the increase in total dose required for tumor control. In conclusion, our results suggest that hypoxia can influence the outcome of carbon ion radiotherapy because of the non-negligible oxygen effect at the low LETs in the SOBP. However, if LOCs occur, a relatively high level of tumor control probability is achievable with a large range of fractionation schedules for tumors with hypoxic subvolumes, but both hyperfractionation and hypofractionation should be pursued with caution.

Keywords: hypoxia, OER, TCP, RCR, carbon ion, fractionation, LOC

INTRODUCTION

Hypoxia is a common cause of treatment failure in photon radiation therapy (RT), as cells that lack oxygen are more resistant to radiation than well-oxygenated cells [1]. The oxygen effect is often quantified in terms of the oxygen enhancement ratio (OER), which is the ratio of doses required to achieve the same biological effect under hypoxic and oxic conditions [1]. Typically, a certain level of cell survival serves as an endpoint. For anoxic cells irradiated with photons, the OER is about 3 for a surviving fraction (SF) of 10% [1]. For low linear energy transfer (LET) carbon ions, the OER is about the same as for photons. With increasing LET, the OER decreases and reaches unity at high dose-averaged LET-values of ∼500 keV/µm [2]. This means that the sensitivity of cells to intermediate- and high-LET carbon ion irradiation is less dependent on the oxygen tension compared with the sensitivity to photon irradiation. Carbon ion RT is thus assumed to result in better control of hypoxic tumors.

In the treatment of tumors with carbon ions, the Bragg peak (BP) is spread out in order to achieve better target coverage. A common range of the dose-averaged LET in the target is 30–80 keV/µm. This range yields an OER of 2–3 [2]. Therefore, the oxygen effect might have an impact on treatment outcome not only in photon RT, but also in carbon ion RT.

Current clinical routine treatment planning does not differentiate between hypoxic and well-oxygenated cells either in carbon ion or in photon RT. In photon RT, however, several different methods for taking tumor oxygenation into account have been proposed, explored and evaluated throughout the years [1, 3, 4]. Regarding carbon ion RT, only very few studies in which the oxygenation of the tumor is considered are available. Bassler et al. have proposed LET-painting, where hypoxic tumor subvolumes are irradiated with higher LET than well-oxygenated regions [5, 6]. Recently, Scifoni et al. presented a possible way to include tumor oxygen status in TRiP (TReatment planning for Particles), the standard carbon ion treatment planning system at GSI (Gesellschaft für Schwerionenforschung) [7]. Uniform biological effect is achieved within the tumor by an optimization approach that compensates for the increase in survival level in the hypoxic regions by increasing the dose to these, or redistributing the LET of the beams.

However, there are so far no publications investigating the effects of fractionation on treatment outcome for carbon ion RT of hypoxic tumors. Since hypoxia possibly plays an important role in carbon ion RT, another benefit of fractionated treatment, apart from the repair of normal tissue between fractions, might be the reoxygenation of hypoxic cells [8].

The mechanism of reoxygenation that is usually referred to in the literature is the slow reoxygenation of chronically hypoxic cells caused by tumor shrinkage [1]. However, another type of reoxygenation, resulting from changes in acute hypoxia between fractions, should also be considered [911]. In a recent study by Lindblom [12], the influence of such local oxygenation changes (LOCs) on the outcome of stereotactic body radiotherapy (SBRT) using photons was investigated. The study showed that extreme hypofractionation (using one single fraction) should be pursued with caution, since it might result in reduced local control as a result of the reduced possibilities of interfraction LOCs. The interplay between total dose, number of fractions, and thus the frequency of LOCs and their impact on the expected local control, should be of interest for carbon ion therapy studies as well, especially when considering the increased tendency towards hypofractionation in this form of RT [1316].

Therefore, it is the aim of this study to theoretically investigate the influence of fractionation on the treatment outcome of tumors with heterogeneous oxygenations irradiated with carbon ions, assuming that LOCs take place between fractions.

MATERIALS AND METHODS

In silico tumor simulation and irradiation

Three-dimensional spherical tumors were simulated in silico with respect to oxygenation, number of clonogenic cells, and tumor size. Distributions of oxygen tension values within the tumor were based on a biologically relevant model of oxygen diffusion and consumption [17]. Distributions of intervessel distances, based on the experimental study by Konerding et al., with different average values were used for the simulation of hypoxic (160 µm) and well-oxygenated regions (100 µm) [18]. Acute hypoxia was simulated by randomly shutting down a fraction of the blood vessel positions. The details of the model can be found in previous publications [17, 19, 20]. A uniform clonogenic cell density was assumed in all calculations.

Two scenarios were considered for the tumor oxygenation: a well-oxygenated tumor, and a tumor containing a centrally positioned spherical hypoxic subvolume (Fig. 1). The diameters of the tumor and the hypoxic subvolume were 4.0 cm and 2.4 cm, respectively. The resulting hypoxic fractions (HFs) of the whole tumor and of the hypoxic subvolume were 16% and 63%, respectively, with a hypoxic threshold of 5 mmHg. The well-oxygenated tumor had a negligible hypoxic fraction (HF = 0.6%). Subsequent calculations were based on the distribution of oxygen tension values within the tumors.

Fig. 1.

Fig. 1.

Oxygen-tension distributions in the well-oxygenated tumor (left) and in the tumor with a centrally positioned hypoxic island (right).

A carbon ion treatment plan was made for the simulated tumor assuming a spherical planning target volume (PTV) with a diameter of 5 cm (0.5 cm PTV margin) centrally positioned in a large water phantom. One beam angle was used. The plan was optimized to give uniform biologically effective dose without taking the oxygenation of the tumor into account. The treatment planning system TRiP was used for the production of the plan [21]. Cross sections of the dose-averaged LET and absorbed dose distributions in the clinical target volume are presented in Fig. 2.

Fig. 2.

Fig. 2.

Cross sections of relative (physical) dose distribution (left) and dose-averaged LET distribution (right) in the clinical target volume. The circles indicate the clinical target volume position (full line) relative to the distributions as well as the center position of the hypoxic island (dashed line). The beam enters from the left.

In order to investigate the influence of the hypoxic island position along the beam axis, and thus relative to the LET and dose distribution, two additional simulation geometries were considered, one with the hypoxic island positioned in the proximal part of the tumor relative to the beam entrance and one in the distal part (± 8 mm off center) (Fig. 3).

Fig. 3.

Fig. 3.

Oxygen-tension distributions in the tumors with the hypoxic island positioned ±8 mm off center, proximally (left) and distally (right). The beam enters from the left.

Calculation of tumor control probability

The LET-parameterized repairable–conditionally repairable damage (RCR) cell-survival model was recently adapted to account for hypoxia by introducing a dose-modifying factor O~ [2224]. However, it can only be used to calculate the surviving fraction in either oxic or anoxic cell populations. In the present study, the model was further developed to take a continuous range of intermediate values of oxygen tension pO2 into account, by assuming the same dependence of O~ on pO2 for carbon ions as the dependence of the OER on pO2 for photons [24, 25, 7]. The dose-modifying factor O~(L,pO2) taking the LET L and pO2 into account, is thus given by the following equations:

O~(L,pO2)=O~(L)k+pO2k+O~(L)pO2 (1)

and

O~(L)=O~min+(O~maxO~min)eLLO~2, (2)

where O~max is the maximum dose-modifying factor in the absence of oxygen at minimal LET, O~min is the limit value at high LETs, LO~ is the LET where O~ is reduced to 37% of its maximum value, and k is the value of the oxygen tension corresponding to the half effect, which was here set to 2.5 mmHg [7, 24, 26]. The surviving fraction S is subsequently described as a function of fraction dose d, LET and oxygen tension, as in equation 3:

S(d,L,pO2)=ea(L)d/O~(L,pO2)+b(L)d/O~(L,pO2)ec(L)d/O~(L,pO2), (3)

where

a(L)=a0f(L)+a1eL/Ln, (4)
b(L)=b1eL/Ln, (5)
c(L)=c0f(L)+c1eL/Ln, (6)
f(L)=1eLLd(1+L/Ld)LdL (7)

and a0, a1, b1, c0, c1, Ln and Ld are free parameters depending on both the intrinsic radiosensitivity of the cells and the ion type [24].

In the present study, equations 2–7 were fitted to previously available and published cell-survival data from experiments of oxic and anoxic human salivary gland tumor (HSG) cells irradiated with carbon ions in the LET range of 22.5–501.5 keV/μm [2]. The dataset consisted of 19 and 20 survival curves, respectively, for the anoxic and well-oxygenated cells, with five to seven different doses for each LET and oxygen tension. The fitting was performed for the whole dataset simultaneously. Further details can be found in a previous publication [24]. The RCR parameters resulting from the fit were used in the subsequent calculations.

The choice of experimental data for this study is motivated by the use of the HSG cell line at the heavy-ion medical accelerator complex (HIMAC) at the National Institute of Radiological Sciences (NIRS) in Chiba and at the Gunma University, both in Japan, for planning biologic dose distributions in their clinical treatment planning systems [27].

The surviving fraction was calculated, voxel by voxel, depending on the oxygen tension, dose and LET of each voxel. The calculations were performed for different fraction sizes, the prescribed dose per fraction ranging from 0.7–10 Gy and the number of fractions, n, ranging from 3–36. LOCs were assumed to take place, meaning that the oxygen tension was assumed to change at the cellular level as a result of changes in acute hypoxia. This was achieved by randomly shutting down a fraction of vessels at each treatment fraction, while the average intervessel distance within each tumor region remained the same throughout the full course of treatment, leading to a redistribution of the pO2 values between voxels at each fraction [12, 17, 19, 20]. The overall geometry of the tumor with its hypoxic island, and the corresponding hypoxic fractions of the different regions, therefore, remained the same for all fractions (Fig. 1). For comparison, all simulations were also performed under the assumption that no local changes in oxygenation take place, i.e. that the oxygenation is static throughout the treatment.

The Poisson tumor control probability (TCP) dose–response model was used to predict the outcome from various combinations of fractionation schedules and oxygen distributions.

TCP=expi=1NvoxNij=1nSi,j(di,Li,pO2i,j), (8)

where Nvox is the total number of voxels, Ni, di and Li are the number of clonogenic cells, dose per fraction and LET in each voxel, and Si,j and pO2 i,j are the surviving fraction and oxygen tension in voxel i at fraction j, respectively. In the case of static oxygenation, the surviving fraction of each cell remains constant during the treatment, and the equation above reduces to:

TCP=exp{i=1NvoxNi[Si(di,Li,pO2i)]n} (9)

The voxel size was 200 µm, and the total number of clonogenic cells in the tumor was set to 106 in all simulations as the TCP values based on this number showed the best agreement with clinically reported values of local control [28]. The dose required to achieve a TCP of 50%, D50, was determined by fitting logit dose–response curves to the resulting TCP values [29].

In order to compare the effects of fractionation and LOCs on the treatment of well-oxygenated and hypoxic tumors, a clinical OER, OERclin, was defined as the ratio of D50 for well-oxygenated (D50,ox) and hypoxic (D50,hypox) tumors:

OERclin=D50,hypoxD50,ox. (10)

All computation steps required for determining the clinical OER are presented in a flow chart in Fig. 4.

Fig. 4.

Fig. 4.

Flow chart describing the computation steps required to determine the clinical OER. Input variables are in elliptic boxes and calculated variables in rectangular boxes.

RESULTS

Figure 5 shows the results of the fit to the HSG cell-survival data from Furusawa et al. for a selection of six different LETs ranging from 22.5–501.5 keV/μm [2]. The fitting parametersa0,a1,b1,c0,c1, Ln and Ld corresponding to the cell-survival curves shown in Fig. 5 are given in Table 1. The goodness-of-fit was assessed by calculating the root mean square deviation of the logarithm of the experimental cell survival from the analytical calculation, RMSD, which was 0.41, and the Pearson product-moment correlation coefficient, r2, which was 0.97. Figure 6 shows the resulting relationship betweenO~, L and pO2.

Fig. 5.

Fig. 5.

Oxic (red lines) and anoxic (blue lines) cell-survival curves resulting from the fit to experimental human salivary gland tumor oxic (dots) and hypoxic (crosses) cell-survival data using the LET-parameterized RCR cell-survival model modified to take oxygenation into account. The experimental data was obtained from from Furusawa et al [2].

Table 1:

Parameters a0, a1, b1, c0, c1, Ln, Ld, Õmin, Õmax and LÕ resulting from the simultaneous fit to the hypoxic and anoxic HSG cell-survival data shown in Fig. 5

a0 a1 b1 c0 c1 Ln Ld O~min O~max LO~
1Gy 1Gy 1Gy 1Gy 1Gy keVμm keVμm keVμm
6.5 2.4 3.4 6.5 0.7 120 76 1.2 3.1 114

Fig. 6.

Fig. 6.

The dose-modifying factor O~ as a function of LET and oxygen tension.

Figure 7a and b show dose–response curves in terms of TCP as a function of total prescribed dose D for the well-oxygenated and for the hypoxic tumors, respectively, assuming that the total dose is delivered in fractions of 0.7 Gy (= 2.1 GyE), 2 Gy (= 6 GyE) and 10 Gy (= 30 GyE). Curves resulting from simulations assuming static oxygenation and LOCs are shown in the same panels. As expected, there are negligible differences for the well-oxygenated tumor between the curves corresponding to the static oxygenation case and assuming LOCs (Fig. 7a).

Fig. 7.

Fig. 7.

TCP as a function of total prescribed dose for fraction doses of 0.7 Gy, 2 Gy and 10 Gy for: (a) a well-oxygenated tumor and (b) a hypoxic tumor.

For the hypoxic tumor, on the other hand, and for low fraction doses (d = 0.7 Gy and d = 2 Gy) there are large differences in the dose required to control the tumor between the case when LOCs take place and the static oxygenation case (Fig. 7b). The dose required to achieve a TCP of 50%, D50, is ∼10 Gy higher if the oxygenation is static compared with the case when LOCs are allowed to occur. For the highest fraction dose of 10 Gy, the difference in the resulting D50 between the two oxygen dynamics cases is much smaller.

For the hypoxic tumor, it should be noted that for d = 10 Gy and d = 2 Gy, the dose–response curves assuming LOCs almost coincide. For the lowest fraction dose (0.7 Gy), D50 is ∼5 Gy higher than for the two higher fraction doses. This means that if LOCs are assumed to take place, the total dose required to control the tumor is almost independent of the size of the fraction dose if the treatment is delivered with doses per fraction between 2 and 10 Gy. This observation is further illustrated by the comparison of the values of D50 that all lie in the interval between 15.5 Gy (d = 5 Gy) and 17.2 Gy (d = 10 Gy) for the fraction doses of 2 Gy, 3 Gy, 4 Gy, 5 Gy, 6 Gy, 8 Gy and 10 Gy. For the lower fraction doses of 1 Gy and 0.7 Gy, D50 is 19.2 Gy and 21.1 Gy, respectively, i.e. 4–5 Gy higher than for the higher fraction doses, as mentioned above.

Another aspect that could also be noted in Fig. 7b is the change in the slope of the dose–response curves when different scenarios are assumed with respect to the oxygenation of the tumor. The shallow curve obtained for d = 0.7 Gy and static oxygenation has a normalized slope of the curve of γ = 3.5. The curve is steeper when LOCs are assumed (γ = 7). The trend is the same for the other fraction doses, and this is explained by the fact that the LOCs lead to an overall increase in the homogeneity in the response of the cells, thus making the dose–response curve steeper, while the shallower curve reflects the larger heterogeneity in the response of the cells with different oxygenations kept constant during the course of the treatment [30].

In Fig. 8a, the relationship between D50 and the number of fractions n is shown. As expected, when the dose per fraction is lowered and the number of fractions is increased, D50 increases, owing to the repair of sublethal damages [31, 32]. However, the trend in the variation of D50 with the number of fractions has to be analyzed in conjunction with the difference in D50 between well-oxygenated and hypoxic tumors. Figure 8b shows the dependence of the clinical OER, defined as the ratio of D50 for hypoxic and well-oxygenated tumors, on the number of fractions. For the static oxygenation case, OERclin lies in the range 1.65–1.75 for all numbers of fractions investigated in this study (3–36). Assuming LOCs, OERclin is <1.2 for all numbers of fractions between nine and 36. In this range, the variation of OERclin with n is low, with the difference in OERclin between two adjacent points being ≤1.8%. If these values are compared with the corresponding ones for extremely hypofractionated schedules, the value of OERclin increases to 1.36 for n = 3, and the difference between two adjacent points becomes considerably larger, namely as large as 18%.

Fig. 8.

Fig. 8.

(a) Total dose required to achieve a TCP of 50%, D50, for different numbers of fractions for well-oxygenated and hypoxic tumors, with local oxygenation changes (LOCs) and with static oxygenation. (b) Clinical OER for different numbers of fractions assuming LOCs or static oxygenation.

In order to investigate the effect of the position of the hypoxic island in the beam direction (z-axis) on the TCP, dose–response curves were determined for a fraction dose of 1 Gy for the case when the hypoxic island is located close to the proximal edge of the tumor relative to the beam entrance or close to the distal edge, as illustrated in Fig. 3. The average values of LET L and the dose modifying factor O~ within the hypoxic island were also calculated for the three different positions and are presented in Table 2 together with the resulting absolute values of D50. When shifting the island position 8 mm off center closer to the proximal edge, where the average LET is lower than in the center of the tumor, D50 increased by 8.8% when assuming LOCs. When shifting the island 8 mm closer to the distal edge, where the LET is higher, D50, instead, dropped by 6.4% compared with the center position. For the static oxygenation case, D50 increased by 9.5% and decreased by 12.0%, respectively. In the x or y directions, on the other hand, the changes in D50 were ≤2.2% for any positional shifts of 8 mm relative to the center position. The explanation for the decrease in D50 when moving the hypoxic island closer to the distal edge is the increase in LET (and thus the decrease in O~) in this region (Fig. 2), as described by equations 1–9. The opposite holds when the hypoxic island is positioned closer to the proximal edge.

Table 2:

Values of the total dose required to achieve a TCP of 50%, D50, of the average LET L and of the average dose-modifying factor Õ of the hypoxic tumor with d = 1 Gy

Position of hypoxic island Average L (keV/µm) Average Õ D50 (Gy) Static D50 (Gy) LOCs
Proximal edge 50.5 1.46 34 20
Center 56.9 1.45 30 19
Distal edge 67.8 1.42 27 18

DISCUSSION

One of the most quoted advantages of carbon ion radiotherapy is the potential ability to control hypoxic tumors, as the oxygen effect is weaker for this type of radiation compared with conventional photon RT [33, 34]. However, this classic paradigm has recently been questioned, since the reported dose-averaged LET-values in the SOBP that covers a typical target are lower than the dose-averaged LET-values in the pristine Bragg peaks of the in vitro experiments that have been used to demonstrate the efficacy of carbon ions in killing anoxic cells [57, 33]. The increased awareness of the role that oxygen might play in carbon ion RT has led to propositions of LET-painting for counteracting hypoxia-related radioresistance and to the implementation of radiobiological models in a treatment planning system aiming to adapt the prescribed dose and beam components according to the oxygenation of the tumor [57].

Nonetheless, clinical experience in the use of carbon ion radiation therapy is growing, and it is demonstrating good outcomes for a large variety of tumor types and sites (including hypoxic tumors) [28, 35]. This indicates that, besides tumor oxygenation, dose and LET, further related aspects should be considered when analyzing clinical data or when predicting the outcome of various clinical schedules.

The present study adds to the previous efforts to include tumor oxygenation at the planning stage of treatment by analyzing the impact of hypoxia and reoxygenation between fractions on the predicted response of tumors with heterogeneous oxygenations. Our results are in agreement with the results of other theoretical studies indicating that a considerable increase in dose would be required to counteract the impact of hypoxic areas if oxygenation is static [57]. However, our results highlight the importance of also considering the dynamics of hypoxia.

The results of the present study clearly show that the dose needed to control hypoxic tumors assuming LOCs is considerably lower than if oxygenation is static, as is indicated by the difference in OERclin for the two cases. As expected, the magnitude of the improvement in TCP depends on the number of fractions (Figs 7 and 8). The improvement decreases for extremely hypofractionated schedules, simply because fewer interfraction oxygenation changes can take place. The OERclin drops rapidly between three and nine fractions, as shown in Fig. 8. If more than nine fractions are employed, the OERclin is rather independent of n. Lengthy schedules would, therefore, not bring a further improvement in OERclin, but they could increase the risk of accelerated proliferation taking place [36].

Another way to analyze the impact of fractionation is to investigate the dependence of the total dose required for tumor control on the size of the dose per fraction. It is well known that the total dose required for tumor control decreases with increasing size of the fraction dose. For well-oxygenated tumors, this is illustrated in Fig. 7a, which clearly shows the differences in dose–response for different fraction doses. This is also the case for hypoxic tumors with static oxygenation (Fig. 7b). However, for hypoxic tumors in which LOCs take place between fractions, the total dose required for tumor control is rather constant for a wide range of fraction doses (2–10 Gy). This means that the expected increase in total dose required for tumor control when going from high to low fraction doses is balanced by the beneficial effects of LOCs over a wide range of fraction doses. The expected increase in total dose due to repair between fractions is only found at fraction doses <2 Gy.

In a clinical study by Nakano et al., the carbon ion treatment outcome for hypoxic tumors was shown to be equally as good as the treatment outcome for well-oxygenated tumors [35]. The treatment was delivered in 24 fractions, which means that LOCs were allowed to occur. The high number of fractions might thus explain the good results for the hypoxic tumors, although the fraction doses employed in that study (0.7–1 Gy) are just below the threshold of what is efficient with respect to tumor control based on our study. The results of the study by Nakano et al. might, therefore, have been further improved, for both hypoxic and well-oxygenated tumors, if fraction doses were slightly increased.

Even though we have identified thresholds regarding both hypo- and hyperfractionation for the specific tumors of our study, no optimal fractionation schedule can be recommended in general for carbon ion radiotherapy of the wide variety of tumors encountered in the clinic. Our study shows that the position of the hypoxic island will have some impact on the dose required for tumor control, especially in the direction of the beam axis along which the average LET changes the most. It is well known that high-LET radiation is more efficient in killing both well-oxygenated and hypoxic cells compared with low-LET radiation [1]. As mentioned in the introduction, the protective effect of hypoxia is also weaker at high LETs. The increased effectiveness of the relatively high LET-values in the distal part of the tumor is therefore expected to be more pronounced for hypoxic than for well-oxygenated cells. Our results, showing that D50 decreases towards the distal part of the tumor, are therefore not surprising. The contribution of our results is, rather, that they give a quantitative estimate of the noted dependencies. Our results also show that the sensitivity of D50 to the position of the island is lower if LOCs are assumed. Other aspects that we did not fully explore, such as the size of the hypoxic island and the fraction of hypoxic cells, are also likely to influence the benefit of fractionating treatment with regard to LOCs.

The simulations in this study have been based on experimental observations regarding the spatial and temporal distribution of oxygen in the absence of tumor shrinkage [911]. These observations strongly support our assumption that LOCs do take place between fractions. Ljungkvist et al. demonstrated, on a microscopic scale, the heterogeneous distribution of well-oxygenated and hypoxic areas in tumor xenografts [9]. Furthermore, using a dual-staining technique, the same research group showed that the LOCs take place within several minutes, thereby demonstrating that the spatial distribution of oxygen in tumors is unlikely to be the same for two consecutive fractions separated by a few hours [10]. These observations were recently confirmed by Redler et al., who used electron paramagnetic resonance measurements to show that local oxygenation in tumors fluctuates within several tens of minutes [11]. With the background of these experimental observations, it appears probable that our results for the LOCs case are representative of clinical tumors.

Generally, it is thought that cells suffering from the two different types of hypoxia, chronic and acute, have the same induced radioresistance, but recently it has been suggested that the radiosensitivity may be modulated by the metabolic energy status of the cells [37]. Based on experimental data, it has also been suggested that chronically hypoxic cells might be more sensitive than acutely hypoxic cells [3840]. The present study, however, does not take into account the difference in sensitivity between the two types of hypoxic cells, since they cannot yet be distinguished in vivo.

Despite the possible improvements in TCP if LOCs are allowed for, it is important to note that hypoxia seems to be an important factor in modulating clinical outcome, as our results show that the OERclin is well above unity both for the static as well as for the dynamic oxygenation cases and for all fractionations. Unless slow reoxygenation of chronically hypoxic cells (caused by tumor shrinkage during the course of treatment) is assumed to take place on top of the local changes in acute hypoxia, it might therefore be beneficial to consider the oxygenation status of tumors at the planning stage of carbon ion RT.

CONCLUSION

In conclusion, our results suggest that hypoxic tumors could be controlled to almost the same degree as well-oxygenated tumors, if LOCs were allowed to take place between fractions. With regard to tumor control, extreme hypofractionation should thus be avoided in order to allow for interfraction LOCs, while extreme hyperfractionation should be employed only with great care, as it increases the total dose required for tumor control and might also lead to accelerated proliferation if overall treatment time is prolonged.

FUNDING

This work was supported by CIRRO—The Lundbeck Foundation Center for Interventional Research in Radiation Oncology, The Danish Council for Strategic Research, Novo Nordisk Fonden and by the Danish Cancer Society Project [ID#DP08023 to N.B.] and Cancer Research Funds of Radiumhemmet [124052 to E.L.]. Funding to pay the Open Access publication charges for this article was provided by Radiumhemmet.

ACKNOWLEDGEMENTS

Michael Krämer and Emanuele Scifoni of the Biophysics Group at GSI, Darmstadt, are thanked for providing access to TRiP98, and together with Walter Tinganelli for constructive comments.

REFERENCES

  • 1.Hall E, Giaccia A. Radiobiology for the Radiologist. 6th edn. Philadelphia: Lippincott Williams & Wilkins; 2006. [Google Scholar]
  • 2.Furusawa Y, Fukutsu K, Aoki M, et al. Inactivation of aerobic and hypoxic cells from three different cell lines by accelerated 3He-, 12C- and 20Ne-ion beams. Radiat Res. 2000;154:485–9. doi: 10.1667/0033-7587(2000)154[0485:ioaahc]2.0.co;2. [DOI] [PubMed] [Google Scholar]
  • 3.Horsman MR, Overgaard J. Hyperthermia: a potent enhancer of radiotherapy. Clin Oncol (R Coll Radiol) 2007;19:418–26. doi: 10.1016/j.clon.2007.03.015. [DOI] [PubMed] [Google Scholar]
  • 4.Toma-Dasu I, Uhrdin J, Antonovic L, et al. Dose prescription and treatment planning based on FMISO-PET hypoxia. Acta Oncol. 2012;51:222–30. doi: 10.3109/0284186X.2011.599815. [DOI] [PubMed] [Google Scholar]
  • 5.Bassler N, Jäkel O, Søndergaard CS, et al. Dose- and LET-painting with particle therapy. Acta Oncol. 2010;49:1170–76. doi: 10.3109/0284186X.2010.510640. [DOI] [PubMed] [Google Scholar]
  • 6.Bassler N, Toftegaard J, Lühr A, et al. LET-painting increases tumour control probability in hypoxic tumours. Acta Oncol. 2014;53:25–32. doi: 10.3109/0284186X.2013.832835. [DOI] [PubMed] [Google Scholar]
  • 7.Scifoni E, Tinganelli W, Weyrather W K, et al. Including oxygen enhancement ratio in ion beam treatment planning: model implementation and experimental verification. Phys Med Biol. 2013;58:3871–95. doi: 10.1088/0031-9155/58/11/3871. [DOI] [PubMed] [Google Scholar]
  • 8.Wang J, Li R, Guo C, et al. The influence of fractionation on cell survival and premature differentiation after carbon ion irradiation. J Radiat Res. 2008;49:391–8. doi: 10.1269/jrr.08012. [DOI] [PubMed] [Google Scholar]
  • 9.Ljungkvist A, Bussink J, Kaanders J, et al. Hypoxic cell turnover in different solid tumor lines. Int J Radiat Oncol Biol Phys. 2005;62:1157–68. doi: 10.1016/j.ijrobp.2005.03.049. [DOI] [PubMed] [Google Scholar]
  • 10.Ljungkvist A, Bussink J, Kaanders J, et al. Dynamics of hypoxia, proliferation and apoptosis after irradiation in a murine tumor model. Radiat Res. 2006;165:326–36. doi: 10.1667/rr3515.1. [DOI] [PubMed] [Google Scholar]
  • 11.Redler G, Epel B, Halpern H. Principal component analysis enhances SNR for dynamic electron paramagnetic resonance oxygen imaging of cycling hypoxia in vivo. Magn Reson Med. 2014;71:440–50. doi: 10.1002/mrm.24631. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 12.Lindblom E. The impact of hypoxia on tumour control probability in the high-dose range used in stereotactic body radiation therapy. Master's Thesis, Stockholm University, 2012.
  • 13.Miyamoto T, Baba M, Yamamoto N, et al. Curative treatment of stage I non-small-cell lung cancer with carbon ion beams using a hypofractionated regimen. Int J Radiat Oncol Biol Phys. 2007;67:750–58. doi: 10.1016/j.ijrobp.2006.10.006. [DOI] [PubMed] [Google Scholar]
  • 14.Miyamoto T, Baba M, Sugane T, et al. Carbon ion radiotherapy for stage I non-small cell lung cancer using a regimen of four fractions during 1 week. J Thorac Oncol. 2007;2:916–26. doi: 10.1097/JTO.0b013e3181560a68. [DOI] [PubMed] [Google Scholar]
  • 15.Okada T, Tsuji H, Kamada T, et al. Carbon ion radiotherapy in advanced hypofractionated regimens for prostate cancer: from 20 to 16 fractions. Int J Radiat Oncol Biol Phys. 2012;84:968–72. doi: 10.1016/j.ijrobp.2012.01.072. [DOI] [PubMed] [Google Scholar]
  • 16.Habermehl D, Debus J, Ganten T, et al. Hypofractionated carbon ion therapy delivered with scanned ion beams for patients with hepatocellular carcinoma – feasibility and clinical response. Radiat Oncol. 2013;8:59–66. doi: 10.1186/1748-717X-8-59. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 17.Dasu A, Toma-Dasu I, Karlsson M. Theoretical simulation of tumour oxygenation and results from acute and chronic hypoxia. Phys Med Biol. 2003;48:2829–42. doi: 10.1088/0031-9155/48/17/307. [DOI] [PubMed] [Google Scholar]
  • 18.Konerding MA, Malkusch W, Klapthor B. Evidence for characteristic vascular patterns in solid tumours: quantitative studies using corrosion casts. Br J Cancer. 1999;80:724–32. doi: 10.1038/sj.bjc.6690416. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 19.Dasu A, Toma-Dasu I, Karlsson M. The effects of hypoxia on the theoretical modelling of tumour control probability. Acta Oncol. 2005;44:563–71. doi: 10.1080/02841860500244435. [DOI] [PubMed] [Google Scholar]
  • 20.Toma-Dasu I, Dasu A, Brahme A. Dose prescription and optimisation based on tumour hypoxia. Acta Oncol. 2009;48:1181–92. doi: 10.3109/02841860903188643. [DOI] [PubMed] [Google Scholar]
  • 21.Krämer M, Scholz M. Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose. Phys Med Biol. 2000;45:3319–30. doi: 10.1088/0031-9155/45/11/314. [DOI] [PubMed] [Google Scholar]
  • 22.Lind BK, Persson LM, Edgren MR, et al. Repairable–conditionally repairable damage model based on dual Poisson processes. Radiat Res. 2003;160:366–75. doi: 10.1667/0033-7587(2003)160[0366:rrdmbo]2.0.co;2. [DOI] [PubMed] [Google Scholar]
  • 23.Wedenberg M, Lind BK, Toma-Daşu I, et al. Analytical description of the LET dependence of cell survival using the repairable–conditionally repairable damage model. Radiat Res. 2010;174:517–25. doi: 10.1667/RR2045.1. [DOI] [PubMed] [Google Scholar]
  • 24.Antonovic L, Brahme A, Furusawa Y, et al. Radiobiological description of the LET dependence of the cell survival of oxic and anoxic cells irradiated by carbon ions. J Radiat Res. 2013;54:18–26. doi: 10.1093/jrr/rrs070. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 25.Alper T, Howard-Flanders P. Role of oxygen in modifying the radiosensitivity of E. coli B. Nature. 1956;178:978–9. doi: 10.1038/178978a0. [DOI] [PubMed] [Google Scholar]
  • 26.Alper T. Cellular Radiobiology. Cambridge: Cambridge University Press; 1979. [Google Scholar]
  • 27.Kanai T, Matsufuji N, Miyamoto T, et al. Examination of GyE system for HIMAC carbon therapy. Int J Radiat Oncol Biol Phys. 2006;64:650–6. doi: 10.1016/j.ijrobp.2005.09.043. [DOI] [PubMed] [Google Scholar]
  • 28.Tsujii H, Mizoe J, Kamada T, et al. Clinical results of carbon ion RT at NIRS. J Radiat Res. 2007;48:A1–13. doi: 10.1269/jrr.48.a1. Suppl. [DOI] [PubMed] [Google Scholar]
  • 29.Dasu A, Toma-Dasu I. Prostate alpha/beta revisited – an analysis of clinical results from 14 168 patients. Acta Oncol. 2012;51:963–74. doi: 10.3109/0284186X.2012.719635. [DOI] [PubMed] [Google Scholar]
  • 30.Dasu A, Toma-Dasu I, Fowler JK. Should single or distributed parameters be used to explain the steepness of tumour control probability curves? Phys Med Biol. 2003;48:387–97. doi: 10.1088/0031-9155/48/3/308. [DOI] [PubMed] [Google Scholar]
  • 31.Elkind M, Whitmore G. Radiobiology of Cultured Mammalian Cells. New York: Gordon and Breach; 1967. [Google Scholar]
  • 32.Ellis F. Dose time and fractionation in radiotherapy. In: Elbert M, Howard A, editors. Current Topics in Radiation Research. Vol. 4. Amsterdam: North-Holland Publishing Co; 1968. pp. 359–97. [Google Scholar]
  • 33.Brahme A. Accurate description of the cell survival and biological effect at low and high doses and LET's. J Radiat Res. 2011;52:389–407. doi: 10.1269/jrr.10129. [DOI] [PubMed] [Google Scholar]
  • 34.Suit H, DeLaney T, Goldberg S, et al. Proton vs carbon ion beams in the definitive radiation treatment of cancer patients. Radiother Oncol. 2010;95:3–22. doi: 10.1016/j.radonc.2010.01.015. [DOI] [PubMed] [Google Scholar]
  • 35.Nakano T, Suzuki Y, Ohno T, et al. Carbon beam therapy overcomes the radiation resistance of uterine cervical cancer originating from hypoxia. Clin Cancer Res. 2006;12:2185–90. doi: 10.1158/1078-0432.CCR-05-1907. [DOI] [PubMed] [Google Scholar]
  • 36.Fowler JF, Tomé W, Fenwick J, et al. A challenge to traditional radiation oncology. Int J Radiat Oncol Biol Phys. 2004;60:1241–56. doi: 10.1016/j.ijrobp.2004.07.691. [DOI] [PubMed] [Google Scholar]
  • 37.Denekamp J, Dasu A. Inducible repair and the two forms of tumour hypoxia—time for a paradigm shift. Acta Oncol. 1999;38:903–18. doi: 10.1080/028418699432590. [DOI] [PubMed] [Google Scholar]
  • 38.Ling CC, Robinson E, Shrieve DC. Repair of radiation induced damage—dependence on oxygen and energy status. Int J Radiat Oncol Biol Phys. 1988;15:1179–86. doi: 10.1016/0360-3016(88)90201-5. [DOI] [PubMed] [Google Scholar]
  • 39.Pettersen EO, Wang H. Radiation-modifying effect of oxygen in synchronized cells pre-treated with acute or prolonged hypoxia. Int J Radiat Biol. 1996;70:319–26. doi: 10.1080/095530096145058. [DOI] [PubMed] [Google Scholar]
  • 40.Ma NY, Tinganelli W, Maier A, et al. Influence of chronic hypoxia and radiation quality on cell survival. J Radiat Res. 2013;54(Suppl 1):i13–22. doi: 10.1093/jrr/rrs135. [DOI] [PMC free article] [PubMed] [Google Scholar]

Articles from Journal of Radiation Research are provided here courtesy of Oxford University Press

RESOURCES