Abstract
Background/Aim: The aim of this study was to investigate radiation-induced tumour vascular damage and its impact thereof on the outcome of stereotactic body radiotherapy (SBRT). Materials and Methods: Vessel densities in animal tumours before and after a single dose of 20 Gy were quantified and used as input for simulations of three-dimensional tumours with heterogeneous oxygenation. SBRT treatments of the modelled tumours in 1-8 fractions were simulated. The impact of vessel collapse on the outcome of SBRT was investigated by calculating tumour control probability (TCP) and the dose required to obtain a TCP of 50% (D50). Results: A radiation-induced increase of acute hypoxia in tumours during SBRT treatment could be simulated based on the experimental data. The D50 values for these tumours were higher than for the simulated tumours without vessel collapse. Conclusion: The vascular changes after high doses of radiation could compromise the outcome of SBRT by increasing tumour hypoxia.
- Hypoxia
- radiation-induced vascular damage
- stereotactic body radiotherapy
- SBRT
- tumour control probability
- TCP
- tumour oxygenation
The modern practice of fractionated radiotherapy is largely based on the so-called 5 R's of radiobiology: repair, repopulation, redistribution, reoxygenation and radiosensitivity (1). For stereotactic treatments such as stereotactic radiosurgery (SRS) and stereotactic body radiotherapy (SBRT) where the treatment is delivered in only one or few fractions, the validity and relevance of the 5 R's has been questioned (2, 3). For example, the steep dose gradients characteristic of stereotactic treatments greatly reduces the impact of normal tissue sparing classically expected from an increased number of fractions (4), while tumour cell repopulation during treatment can be neglected given the short overall treatment time in treatments in few fractions. Furthermore, the effect of redistribution can be disregarded given that high doses are likely to cause cell-cycle arrest (5). The limited treatment time and low number of fractions also imply that reoxygenation is limited, indicating that hypoxia has an even bigger impact on the outcome of SRS and SBRT compared to conventionally fractionated radiotherapy (3, 6, 7). While slow reoxygenation following chronic hypoxia is typically associated with tumour shrinkage, it is highly unlikely to happen during a hypofractionated treatment course where fast reoxygenation of acute hypoxia is expected to occur between individual fractions (8). Reoxygenation, and the closely related radiosensitivity of tumour cells, thus remain highly influential parameters in stereotactic radiotherapy. Recently, there has been an increasing interest in the study of additional effects occurring only at high doses per fraction, particularly the damage to the tumour vasculature. Based on experimental observations, it has been suggested that damage to the tumour vasculature could have a secondary effect resulting from the sudden loss of blood supply and/or from an increase in the radiosensitivity of tumour cells due to a cytotoxic effect derived from endothelial apoptosis (2, 9, 10). While the vascular effect has, thus, almost exclusively been discussed as a potential second cytotoxic factor for tumour cells, the observed vascular response after doses of 10-15 Gy appears to be highly dynamic (2). With an initial reduction in blood flow that persists for various lengths of time, the actual vascular effect with respect to the treatment efficacy could therefore be highly dependent on the time point at which treatment fractions are delivered with respect to each other. Therefore, instead of increased cell death, increased radioresistance could be expected during the limited time of e.g. an SBRT treatment if the percentage of acutely hypoxic cells is temporarily increased from one fraction to the next. Thus, there is a gap in understanding how dose-fractionation actually affects the tumor vascular microenvironment, and impacts on the tumor control. In this study, radiobiological experiments and modelling were combined to investigate the changes in the tumour vasculature resulting from high doses and the subsequent impact on the tumour control probability (TCP) in stereotactic radiotherapy treatments.
Materials and Methods
The study consisted of two parts, one experimental part in which data on tumour vasculature before and after a dose of 20 Gy was collected and analysed, and one modelling part in which the experimental data was used to simulate the effect from the observed vascular response at the level of the whole tumor.
Mouse irradiation, multiphoton imaging and analysis. All procedures for animal experiments were performed according to City of Hope guidelines and approved by the Institutional Animal Care and Use Committee. Four NOD/SCID 16-18 week old female mice were injected with 105 E0771 CEA+ breast cancer cells suspended 50 μl PBS and matrigel in equal parts. Cells were injected onto the mammary fat pad. Two of the mice were anaesthetised and treated with 20 Gy using a targeted image-guided dose-calculated irradiation system (Precision X-Ray, North Branford, CT, USA) (11). Treatment was performed 10 days post E0771 injection. One day post treatment, mice were anesthetized initially with 3% isoflurane at a flow rate of 2 l/min mice. Retro-orbital injections of 500 μg of 70 kDa, FITC-conjugated fluorescent Dextran (Sigma Aldrich, St. Louis, MO, USA) suspended in 100 μl PBS were given to mice for imaging tumourvascular lumen. Mice were placed on a custom built heated stage maintained at 37°C, and anaesthesia was reduced to 1.5% isoflurane at a flow rate of 1.5 l/min for surgery and imaging. A midline incision to the skin was made. The skin, tumour, and mammary fat pad were separated from the peritoneum preserving vasculature surrounding the tumour. A glass coverslip was placed over the upturned mammary fat pad and tumour to enable proper imaging. Phosphate buffered saline was applied under the coverslip to preserve tissue function. Multiphoton microscopy was performed using a Prairie Ultima multiphoton microscope (Bruker Corporation Billica, MA, USA). Fluorescent excitation was performed using a Chameleon Ultra II tunable Ti:Sapphire laser with 140 femtosecond pulses (Coherent, Santa Clara, CA, USA). A 10x objective (UPlanFL N 10X/0.30 ∞/−/FN26.5, Olympus, Tokyo, Japan) was used to acquire all images. Tiled three-dimensional images were acquired by collecting a grid series of 15% overlapping z-stack images with a slice spacing of 15 μm. Data acquisition using the Prairie Ultima microscope was handled by Prairieview 5.3 software. Tiled images were stitched together using a Fiji/ImageJ grid collection/stitching plugin (12, 13). Fiji/ImageJ was used to manually quantify vessels in the tumour. Vessel branches were counted from three dimensional images manually and normalised for image area. Six square subregions of approximately 0.85 mm2 were analysed for each tumour. The regions were then used for measurements of intra-tumour vessel density variation including the calculation of the median, mean and standard deviation in each tumour.
Radiobiological modelling. To simulate the effect of radiation given in few fractions from the experimentally observed vascular responses at the level of the whole tumor, a previously developed three-dimensional in silico tumour model of heterogeneous tumour oxygenation and radiation response was used (14). The input to this model is distributions of inter-vessel distances (IVDs) previously obtained through corrosion casting in animal solid tumours (15). In the model, the tumour is geometrically defined including potential sub-volumes to which different IVDs can be assigned. In this way, several IVDs distributions (with mean values of 60-160 μm) can be combined in the same tumour to simulate tumours with a broad range of oxygenation distributions. Blood vessel capillaries are randomly placed in the defined volumes according to the assigned IVDs distribution, and a fraction of the vessels is randomly removed in order to simulate acute hypoxia resulting from a temporary collapse or occlusion of capillaries. The three-dimensional tumour map of oxygen partial pressure (pO2) is then calculated by solving a differential equation describing the diffusion and consumption of oxygen from the vessels. Due to the random nature of the exact spatial distribution of active vessels, fast reoxygenation can be simulated by repeating the assignment of IVDs distributions and re-calculating the pO2 map (14).
In order to construct in silico tumours representative of the experimental tumours in the present study, different combinations of blood vessel distributions characterised by different IVDs resulting in the experimentally determined mean vessel densities and standard deviations were sought. While the simulations produce three-dimensional models of the whole tumour, the experimental measurements were obtained from a limited number of samples. Therefore, combinations of several IVDs were investigated for the purpose of finding a mean vessel density and standard deviation matching the experimental vessel density assuming equal weights between the IVDs.
In order to investigate the impact of radiation-induced vessel collapse on the treatment outcome, SBRT treatments delivered in 1-8 fractions were simulated by delivering a homogeneous dose to the modelled tumours representing both the untreated (i.e. without vessel collapse) and treated (with vessel collapse) experimental tumours. In order to construct dose-response curves, a range of doses was considered. The cell surviving fraction SFi in each voxel i was subsequently calculated using the linear-quadratic (LQ) model taking into account the pO2 on voxel level by modifying the radiosensitivity parameters α and β of the LQ model with oxygen modifying factors (OMFs) (14, 16):
(1)
In order to evaluate the simulated treatments, the tumour control probability (TCP) was calculated as the product of the individual TCPs in each voxel, assuming a homogeneous density N of clonogenic cells at the start of the treatment:
(2)
The parameter values used in the simulations were α=0.35 Gy−1, β=0.035 Gy−2 (α/β=10 Gy), OMFmax=3 and N=108 cells (17).
Results
Figure 1 and Table I show the experimentally determined vessel densities in four animal tumours. Both the highest and lowest average vessel density in the two untreated tumours were found to be 18.0 mm−2 and 138.0 mm−2, respectively. The average vessel densities in the treated tumours were between 71.7 mm−2 and 62.0 mm−2. In the untreated tumour with low average vessel density, the six samples were extremely similar in comparison with the considerable heterogeneity observed in the three other tumours (Figure 1). Given the well-known heterogeneous nature of tumour vasculature, this could indicate that the six samples are less representative of the tumour vasculature in the case of tumour 1, and the simulated tumour before treatment was therefore based on the data for tumour 2. Based on the comparison of the IVDs with the experimental results for the two treated tumours, the simulation of tumours after a high dose was based on the data for tumour 4. Two patterns of hypoxia representative of clinical findings were considered: one in which the tumour oxygenation decreases towards the centre of the tumour, and one with a more heterogeneous pattern of multiple islands not located in the centre of the tumour (18). Figure 2 shows the cross sections through the simulated tumours with and without radiation-induced vessel collapse.
Figure 2A and B should thus be considered to represent the same tumour with and without radiation-induced vessel collapse according to the experimental findings (and the same for the tumours illustrated in Figure 2C and D).
Figure 3 shows the dose-response curves obtained for the tumours with central hypoxia, and Figure 4 shows the dose-response curves obtained for the tumours in which islands ofhypoxia were present.
In Table II, the D50 values of the dose-response curves in Figures 3 and 4 are presented.
The values of D50 are in general agreement with the doses clinically prescribed in SBRT (19, 20). Furthermore, the D50 is higher for the tumours with radiation-induced vessel collapse as a result of the fact that they are more hypoxic.
Discussion
With an increasing use of stereotactic body radiotherapy and other advanced radiotherapy techniques employing fundamentally different dose distributions and fractionation schedules, the classical 5 R's describing the response of late- and early-reacting tissues to fractionated radiotherapy have been challenged in recent years. Several experimental and clinical studies have reported interesting observations of effects that seem to occur only at high doses per fraction. In melanoma patients treated with a combination of high-dose radiotherapy and the monoclonal antibody ipilimumab, an abscopal effect was observed as lesions outside of the radiation field regressed after the treatment (21). This has been suggested to be a result of enhanced anti-tumour immunity after irradiation, but the rather limited clinical findings so far warrant further investigation of whether this phenomenon is specifically attributed to SBRT-like treatments and if other tumour types besides melanoma exhibit the effect (3). In addition to an interesting potential interplay with the immune system, the administration of high doses has been associated with damage to the tumour vasculature which might have a secondary effect in stereotactic treatments. This has been hypothesised to happen for two reasons: secondary tumour cell death from the sudden loss of blood supply and increased radiosensitivity of tumour cells due to a cytotoxic effect derived from endothelial apoptosis (2, 9, 10). While the discussion has been focused on a contribution to cell death, it should be pointed out that vessel collapse also occurs spontaneously in tumours, causing perfusion-limited acute hypoxia (8). This leads to an increased resistance to radiation, as opposed to cell death from acute deprivation of nutrients. Acute hypoxia is a transient state that can be reversed by fast reoxygenation before the affected cells are permanently damaged. If radiation-induced vessel collapse can also be reversed in due time, an increase in the radioresistance of the tumour could thus equally well be expected following a high dose. In fact, radiation-induced damage to the tumour vasculature has been observed to be repaired and the vasculature is restored (2). This indicates that the actual vascular effect could imply a reduction rather than an increase in tumour cell death, and highlights the importance of the timing of e.g. SBRT fractions. A recent modelling study demonstrated this effect for an assumed vessel collapse of 35% caused by the first fraction and persisting through the remaining treatment (22). This led to a greater loss of tumour control probability than of double sized chronically hypoxic tumour compartment.
Vessel densities in the six subregions for each of the experimental tumours.
The mean value of vessel density in each tumour with the corresponding standard deviation (Stdv) and coefficient of variance (CV).
Illustration of (A) simulated tumour without vessel collapse, hypoxic fraction HF <5 mmHg=5%, (B) corresponding tumour with vessel collapse, HF=25%, (C) simulated tumour without vessel collapse, hypoxic fraction HF=5%, (D) corresponding tumour with vessel collapse, HF=15%.
In the present study, the response of the tumour vasculature to a high dose was studied experimentally, and the therapeutic implications thereof were investigated by radiobiological modelling based on the experimental data. Four model tumours representative of the experimental tumours were simulated, two rather well-oxygenated tumours with no presumed vessel collapse, and two more hypoxic tumours with radiation-induced vessel collapse. The D50 values of all simulated tumours were in agreement with the clinically prescribed doses in SBRT, with the exception of the single-fraction schedules in line with previous findings (6). For the tumours with vessel collapse, the D50 was consistently higher than that in the simulated tumours without vessel collapse.The biggest differences were observed for the most hypoxic tumour with a hypoxic fraction ≤5 mmHg of 25%. Although it is well-known that higher dose is required in the more hypoxic tumour (23), the fact that the simulated tumours in this study were based on experimental tumours with and without radiation-induced vessel collapse provides support to the the conclusion that the vascular effect from high doses could lead to hypoxia. The simulations in turn show that this could have an impact on the treatment efficacy of e.g. SBRT through an increase in radiation resistance of the tumour.
Dose-response curves for the simulated SBRT treatments of 1-8 fractions delivered to a tumour with centralised hypoxia with and without vessel collapse, (A) and (B) respectively.
Dose-response curves for the simulated SBRT treatments of 1-8 fractions delivered to a tumour with multi-island hypoxia with and without vessel collapse, (A) and (B) respectively.
Tumour control probability (TCP) for the simulated treatments.
In the analysis of the experimental data, one of the untreated tumours was excluded due to the very low vessel density and limited spread of sample values. As mentioned in the Results section, this could indicate that the samples acquired were not representative of the vasculature of the tumour, in which a higher degree of heterogeneity would be expected as has been observed in the other experimental tumours and in numerous previous studies (24, 25). More importantly, the inclusion of a homogeneously hypoxic tumour as representative of a case without radiation-induced vessel collapse would compromise the purposes of the present study, focusing on the potential detrimental effects resulting from increased tumour hypoxia as a result of significant treatment-induced vascular damage. Furthermore, in the clinical setting, this effect would not be of major (if any) concern for such highly hypoxic tumours, which would be considered high-risk and in potential need for hypoxia-targeted treatment adaptation before the start of the treatment. Nonetheless, the number of tumours included in the experimental part of this study is small, and the measurements of the tumor vasculature were perfomed only 48 h post-SBRT. In order to expand the analysis and better understand the effects derived from vascular damage following high doses radiation, further experiments and modelling studies are planned in the near future.
In this study, experimental data on the changes in the tumour vasculature following a high dose of radiation were used to simulate tumours with increased acute hypoxia resulting from radiation-induced vessel collapse. Depending on the extent of vessel collapse and the prescribed dose, vascular changes after high doses of radiation could compromise the outcome of SBRT by increasing tumour hypoxia.
Acknowledgements
Research reported in this publication is supported by the National Cancer Institute of the National Institutes of Health under award number P30CA033572 and partly supported by National Institutes of Health grant 1R01CA154491-01. One of the author acknowledges discussion and suggestions provided by Professor John Shively, City of Hope. Financial support from the Cancer Research Funds of Radiumhemmet (Solna, Sweden) is gratefully acknowledged.
Footnotes
Authors' Contributions
EKL performed the simulations and wrote the manuscript. JB and SH performed the experimental work and subsequent analysis. AD and ITD provided the model for tumour oxygenation and ITD produced input data for the simulations. All Authors reviewed the manuscript.
This article is freely accessible online.
Conflicts of Interest
The Authors have no conflicts of interests to declare.
- Received May 4, 2019.
- Revision received May 15, 2019.
- Accepted May 16, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved
