Radiotherapeutic treatment options for oligotopic malignant liver lesions

Background Several radiotherapeutic approaches for patients with oligotopic malignant liver lesions unfit for surgical resection exist. The most advanced competitive techniques are high-dose-rate (HDR) brachytherapy, Cyberknife, volume-modulated-arc therapy (VMAT) and Tomotherapy. We evaluated the optimal technique by a planning study for a single ablative dose with different lesion sizes. Methods We compared dose distributions of HDR-brachytherapy with stereotactic ablative radiotherapy using the Cyberknife, VMAT or Tomotherapy. Tumor-control-probabilities (TCP), normal-tissue-complication-probabilities (NTCP) were determined in a theoretical framework applying a single dose of 20 Gy (demanding 95% coverage) for intrahepatic lesions of 1–5 cm in size. We evaluated therapeutic ratios by TCP (mean dose in the lesion) relative to high-dose (conformality) or low-dose liver exposition in dependency on the lesion size for each technique. In addition, we considered treatment times and accuracy (clinical target volume vs planning target volume). Results HDR-brachtherapy has the highest therapeutic ratios with respect to high-dose as well as low-dose liver exposition even for extended lesions, and the Cyberknife being suited second best. However, for lesions ≥ 3 cm diameter the therapeutic ratios of all ablative techniques are increasingly converging, and better tolerance and shorter treatment times of noninvasive external techniques become more important. On the other hand, mean tumor doses of HDR-brachytherapy of near 60 Gy are unattainable by the other techniques gaining only 22–34 Gy, and the conformality of HDR-brachytherapy is still rather good for lesions ≥ 3 cm diameter. Conclusions HDR-brachytherapy is by far the most effective technique to treat intrahepatic lesions by a single fraction, but sparing of the surroundings declines with increasing lesion size and approaches the benchmarks of external beam radiosurgery techniques. External beam radiotherapy has the advantage to use suitable fractionation schedules.


Background
The liver is a common target for metastases from various solid cancers, especially from colorectal cancer and other gastrointestinal cancers. Surgical resection of metastases has been shown to improve survival in colorectal cancer [1], thus surgery is regarded as the preferred local treatment strategy for oligotopic colorectal cancer liver metastases. However, only less than one third of patients are appropriate candidates for surgical resection due to location and distribution of lesions within the liver, patient co-morbidities prohibiting surgical intervention or extrahepatic involvement [2,3]. Despite major advances in systemic treatments for colorectal cancer during the last decades still long-term control of liver disease is almost impossible without additional local treatment of metastatic disease [4]. Therefore, other local Open Access *Correspondence: pirus.ghadjar@charite.de 1 Charité -Universitätsmedizin Berlin, corporate member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Department of Radiation Oncology, Augustenburger Platz 1, 13353 Berlin, Germany Full list of author information is available at the end of the article treatment strategies for patients with oligotopic malignant liver lesions are warranted to complement systemic therapies in cases not qualifying for surgical resection [3,5]. There are some general restrictions to the use of radiofrequency ablation (RFA) and other thermoablative procedures, which include danger to thermosensitive structures, limitations to treat lesions in proximity of larger vessels due to cooling, ineffectivity to treat lesions with increased arterial tumor perfusion such as hepatocellular carcinoma or neuroendocrine tumors and limited local control of lesions > 3 cm size [5,6].
Image-guided high-dose-rate interstitial (HDR-) brachytherapy uses an iridium-192 source of < 1 mm diameter, which is moved in an implanted catheter array according to a treatment plan in order to generate an optimal dose distribution for a specified target volume and adjacent organs of risk. It overcomes the mentioned limitations of RFA [6] and was shown to be an effective treatment option [7][8][9][10]. Likewise, a variety of modern external beam radiation techniques are available using dedicated linear accelerators, competing with the interstitial rather accurate but invasive method. This includes intensity modulation radiotherapy (IMRT), rotational techniques such as volume modulated arc therapy (VMAT) or Tomotherapy and the Cyberknife, which are suitable for radiosurgery in a single fraction [3,11].
Other non-surgical treatment options do exist for advanced unresectable liver metastases including transarterial chemoembolization, selective internal radiation therapy and hepatic arterial infusion, which are mainly useful for disseminated liver disease. In this article, we focus on oligotopic malignant liver lesions that are treated in a single fraction as this is a logistical and practical advantage for the patient, however without evidence of improved oncological outcome compared to fractionated schedules.
It is the task of the radiation oncologist to ascertain a well-suited radiotherapy technique for a given target volume embedded in an individual anatomy with specified organs of risk or regions of interest (ROI). This article compares radiotherapeutic single dose approaches such as HDR-brachytherapy, rotational techniques (VMAT, Tomotherapy) and Cyberknife. For that reason we calculated and analyzed optimal dose distributions for a single dose of 20 Gy covering intrahepatic lesions between 1 and 5 cm diameter outlining the pros and cons for these different techniques.

Methods
We specified intrahepatic spherical liver lesions in a peripheral location of 1-5 cm diameter and segmented adjacent ROI (liver, right/left kidney, stomach, heart and right lung) to estimate toxicity. For all approaches treatment planning dose calculations were optimized demanding mean doses D mean as high as possible inside the lesions ensuring the prescribed standard dose of 20 Gy with target coverage of > 95% and minimal burden to the liver. Listed parameters include D mean in the lesions characterizing the TCP, and V 10Gy , V 5Gy , V 2Gy and V 1Gy of the liver, and maximum doses D max of other critical organs. The parameters for any other prescribed dose can easily be calculated by rule of three.
The following four radiation techniques were evaluated: HDR-brachytherapy We assumed one or five accurately implanted catheters in the lesion and used the planning system Brachyvision ® . The planning target volume (PTV) and clinical target volume (CTV) was set equal, because catheters do not move relative to the lesion. We note that catheter implantation requires an invasive image-guided procedure that must be accurate to the millimeter, and is demanding for patients and staff and involves some risks (bleeding, infection etc.). Cyberknife This radiosurgery system uses a compact linear accelerator on a robotic manipulator that can radiate from a large spherical angle using standard circular collimators of selected diameter. The principal advantage is the full three-dimensional performance and the integrated image guidance system that enables on-line tracking during radiation. An almost perfect tracking is implemented for Cyberknife based on a gold marker implantation that an image-guided invasive procedure similar to the catheter implantation for brachytherapy, but less critical in terms of accuracy. The PTV and CTV was set equal. An inverse planning algorithm is available. Volume Modulated Arc Therapy (VMAT) Using a standard linear accelerator (Siemens/Varian Medical Systems Inc., Palo Alto, 94304 CA, USA) with one or more arcs rotating around the iso-center typically located in the lesion with variable speed, intensity and multileaf collimator (MLC) setup. Leafs have thickness of 0.5 cm in the center (20 × 20 cm). To account for breathing-induced movement of the liver the diameter was enlarged by one cm to form the PTV. To restrict this safety margin to only 5 mm we assumed patient compliance to control respiratory motion either with shallow breathing or breathhold techniques in addition to cone-beam CT for image guidance. Optimization was performed with respect to dose specifications for the PTV and the ROI (limiting V X% , volume with doses ≥ X%) by use of the planning system Eclipse ® . Tomotherapy Using a binary (on/off ) collimator of 1 cm thickness that allows full intensity modula- Table 1 Parameters to calculate normal tissue complication probabilities (NTCP) for conventionally fractionated regimens Neuropathy/plexopathy tion slice by slice in a dedicated linear accelerator (Accuray Inc., Sunnyvale, CA 94089, USA). Tomotherapy is superior for complex target volumes and complicated constraints, but is theoretically nearly equivalent to VMAT for simple targets like spheres. Again the CTV was enlarged by one cm in diameter to form the PTV assuming adequate control of respiratory motion and megavoltage CT for image guidance. The planning system Tomoplan ® calculates optimized dose distributions with respect to the selected restrictions.
We used the Lyman/Kutcher/Brown (LKB) model to estimate normal tissue complication probabilities (NTCP) for the liver [12][13][14][15][16] in order to evaluate longterm toxicity. In short, the LKB-model exploits the tolerance dose TD 50/5 implying a specified complication in 50% in 5 years after irradiation of the entire organ. This dose is the mean of a Gaussian distribution of tolerance doses with standard deviation σ = m × TD 50/5 . Integration over this probability function yields the cumulative distribution function that is the sigmoid NTCP-curve. The parameter m depends on the critical organ (here liver) and determines the slope of the NTCP-curve. Another crucial parameter is n describing the volume dependencies of the tolerance doses. The partial volume is υ i = V i / V 0 with V 0 the volume of the entire organ. Then the tolerance dose TD 50/5 (υ i ) of the partial volume is larger than TD 50/5 (υ 0 = 1) and can be calculated utilizing the volume parameter n by the formula TD 50/5 (υ i ) = TD 50/5 (1)/υ i n .
For better comparison we calculated therapeutic ratios using the term D mean (lesion d)/(%liver > 10 Gy) comparing TCP with high-dose exposition (conformality) and the term D mean (lesion d)/(%liver > 1 Gy) comparing TCP with low-dose exposition (long-term risk). Table 1 is listing the volume depending parameters TD 5/5 , TD 50/5 , σ (for the entire volume 1 or a partial volume 1/3) as well as m, n and α/β for various organs or tissues in comparison to the liver. The tolerance doses TD and their volume dependencies assume a conventional fractionation [17]. By using the LKB-model and Table 1 we constructed the NTCP-curves for the liver in Fig. 1 with important reference points for illustration [18].

Results
We note that m = 0.25 for the liver is in the upper range, i.e. the tolerance doses are not as sharply defined. In consequence, for a rather moderate dose of 10 × 2 Gy → 20 Gy (TD 50/5 − 2σ) to the whole liver we expect a radio-induced liver disease (RILD) in more than 2%. On the other hand, n = 0.6 is relatively high, i.e. for partial volumes we achieve much higher tolerance doses. For a partial liver irradiation of one third of the liver the TD 5/5 (1/3) is shifting to 50 Gy and the TD 50/5 (1/3) even to 80 Gy.
We used the linear quadratic model to calculate biologically effective doses BEQ 2Gy (assuming conventional fractionation with daily doses d = 2 Gy) iso-effective to single doses D S for tissues with different α/β values [19] in order to compare the abovementioned techniques that apply high single doses in circumscribed lesions ( Table 2): In the last column for α/β = 10 Gy we added the isoeffective biologically effective doses BED (assuming a low-dose-rate regimen with d → 0), for easier comparison with published data (see "Discussion" section): The liver is particularly sensitive to large single doses because of the low α/β = 1.5 Gy. Therefore, only two fractions of 5 Gy to the entire liver (iso-effective to 20 Gy conventionally fractionated) might induce a RILD in as much as 2% of patients.
The dose distributions of the considered treatment techniques applied to a lesion with d = 2 cm are depicted in Fig. 2. Important DVH parameters for the lesions and the ROI liver, right kidney and stomach are listed in Table 3 for lesion sizes d = 1, 2, 3, 5 cm.
HDR-brachytherapy is by far most effective gaining intralesional mean doses D mean > 55 Gy, whereby increasing the number of catheters (5 vs 1 catheter) homogenizes the intralesional dose distribution and lowers D mean . We achieved much lower D mean of approximately 30 Gy by Cyberknife and Tomotherapy and only 22 Gy by VMAT. For lesion sizes < 3 cm the rather high TCP is combined with improved conformality in the high-dose range around 10 Gy, sparing of the liver in the low-dose range around 1-2 Gy and low exposition of the other surrounding organs (e.g. right kidney, stomach) in comparison to all external radiation techniques. For lesions > 3 cm D mean is still superior for HDR-brachytherapy, but dose sparing in liver and surroundings is continuously declining with increasing size. For lesions of 5 cm size the low dose exposition with 1-2 Gy is even higher for HDRbrachytherapy in comparison to all external techniques and equal for exposition with 5 Gy, but the strikingly high dose to the tumor remains. Interestingly, more catheters (5 vs 1 catheter) neither improve the effectiveness nor the sparing of surroundings.
HDR-brachytherapy achieves the highest conformality for small lesions that also the dose distribution reflects ( Table 3, Fig. 2). For example for the 1 cm lesion only 0.5 ml healthy liver is exposed to > 10 Gy versus 1.3 ml for Cyberknife and around 20 ml for VMAT/Tomotherapy. Also the sparing of the liver in the low dose range is excellent, e.g. for the 1 cm lesion only 70 ml healthy liver exposed to 1 Gy versus 200-400 ml with external techniques. We emphasize that these advantages require a nearly perfect positioning of the catheter(s). For all lesion sizes the conformality of the Cyberknife is between HDR-brachytherapy and VMAT/Tomotherapy, but low dose expositions (volume > 1-2 Gy) are in the same range for all external techniques and higher than for interstitial techniques.
The VMAT technique that is standard in modern linear accelerators performs with lower intralesional mean doses of ~ 22 Gy accompanied by the lowest liver exposition at 1 Gy for lesions ≥ 3 cm. Tomotherapy achieves higher intralesional mean doses ~ 30 Gy with a larger load to the surroundings accordingly. The isodoses for higher doses (5-10 Gy) are more compact for Tomotherapy than VMAT (Fig. 2).
For a conclusive comparison we plotted the therapeutic ratios of the radiotherapy techniques in dependency on lesion sizes for high intrahepatic doses > 10 Gy (Fig. 3, left) and low intrahepatic doses > 1 Gy (Fig. 3, right). The HDR-brachytherapy performs always superior, but these advantages decline with increasing lesion size. The Cyberknife comes second and the external standard techniques (VMAT, Tomotherapy) come third. While for a lesion of d = 1 cm the therapeutic ratios are far apart by two decades (high dose) or one decade (low dose), the ratios differ only by a factors 2-3 for d = 5 cm.
The performance of the (typically invasive) Cyberknife lies in between HDR-brachytherapy and VMAT/Tomotherapy, but again the advantage declines with increasing lesion size.
VMAT shows a slight advantage over Tomotherapy for small lesions, mainly due to finer leafs in the center (0.5 cm vs 1 cm), while Tomotherapy performes better for large and complex target volumes (as illustrated in Fig. 3, right). Table 3 shows also relevant differences in the irradiation times that might influence the treatment decision. Each Cyberknife treatment session consumes about 1 h in comparison to only 2 min for VMAT and 3-6 min for Tomotherapy. In case of HDR-brachytherapy treatment time is rather short for small lesions and increases considerably with lesion size, e.g. around 20 min for lesions of 5 cm assuming a recent exchange of the iridium-source and up to 40 min, if the source change is long ago (typically about two months).  In addition, we must consider for a treatment decision that feasibility, invasiveness and risks of the techniques greatly differ.

Discussion
We have compared four radiosurgery techniques for a standard setup with equal optimization conditions (see "Methods" section) and ascertained image-guided HDRbrachytherapy as by far the most effective radiotherapy technique to control locally oligotopic malignant liver lesions in a single session. This is an option for patients with metastases either unresectable and/or unfit for surgical resection and remains valid even for larger lesions in our planning study. However, sparing of the surroundings by HDR-brachytherapy declines with increasing lesion size and finally approaches the benchmarks of external beam radiation techniques. Notably, HDR-brachytherapy generates an inhomogeneous dose distribution with dose excesses in the target around the catheters achieving intralesional mean doses of > 55 Gy. This is important as rather high single doses are estimated to be required for reliable and long-term tumor control.
Ricke et al. [9] conducted a dose escalation study comparing prescribed target doses of 15 Gy, 20 Gy versus 25 Gy. They found a long-term local control (> 3 years) in 95% for 25 Gy (BED = 88 Gy) and only 65% for 20 Gy (BED = 60 Gy) according Table 2. Local control was also correlated with D min (> 20.5 Gy) indicating that the coverage is important because of the steep dose gradient at the tumor margin [9,20]. Herfarth et al. 2001 [21] achieved for liver tumors/metastases lower local control of 80% for 22-26 Gy single fraction in the reference point (80% of this dose enclosing the PTV).
Numerous clinical trials with various fractionation schemes were recently reviewed (3). Generally, a BED > 100 Gy was recommended to achieve a satisfactory long-term local control. In a careful evaluation a BED near 140 Gy was estimated for a long-term (> 1 year) local control of 90% for colorectal cancer that turned out as the most resistant cancer type [22]. Such high BEDs require single doses > 30 Gy ( Table 2) that is in our analysis most likely accomplished by HDR-brachytherapy.
The LQ-model confirms that rather high single doses might be necessary for local control. Table 4 illustrates that single doses of 20-30 Gy are only sufficient for Fig. 2 Dose distributions comparing the four radiotherapy techniques, if the standard single dose 20 Gy is given to a 2 cm lesion. The isodoses range from 150% (30 Gy, orange), 100% (20 Gy, white), 75% (15 Gy, yellow), 50% (green, 10 Gy), 25% (5 Gy, light blue), 10% (2 Gy, purple) to 5% (1 Gy, light green) according to the inserted scales. The ROIs are also shown such as liver (brown), stomach (blue) and kidneys (green, red). For such a small lesion HDR-brachytherapy is clearly superior in comparison to the external techniques radiation sensitive tumors (with α = 0.3 Gy −1 ) and/ or tumors with low α/β < 10 Gy. However, radioresistant tumors (with α < 0.2 Gy −1 and α/β > 5 Gy) might get missed with single doses below 30 Gy. Clinical data suggest that α/β of metastases are typically > 10 Gy [21] resulting in higher single doses > 30 Gy for reliable local control. Single equivalent dose of ≥ 45 Gy achieved a two-year local control of 100% compared with 54% for < 45 Gy [11]. Our planning study shows ( Table 3) that external methods have difficulty depositing such high single doses is the lesions.
HDR-brachytherapy performs with a steep dose gradient sparing the surrounding tissues and reducing their NTCP. Both features favor in particular HDR-brachytherapy for intrahepatic metastases or tumors, because the liver is sensitive to radiation exposure in a large volume with low or medium radiation doses. For the liver much care is required, if large volumes or even the entire organ are irradiated given the low TD 5/5 (entire liver) of 25 Gy (conventional fractionation).
CT-guided implantations of brachytherapy catheters and implantation of a gold marker for Cyberknife tracking are similar minimally-invasive procedures, sometimes comparable in the risk. In case of the Cyberknife the risk of gold marker implantation must be balanced against the benefit of enhanced precision by reducing the safety margin (taking CTV as PTV).
On the other hand, the quality of HDR-brachytherapy critically depends on the precise placing of the catheters, while the geometry of gold marker(s) near or in the target is not as critical. It must be considered that HDR-brachytherapy and Cyberknife are more invasive and includes some risks, if compared with non-invasive external radiation methods (VMAT, Tomotherapy). If we address single session treatment approaches using either HDR-brachytherapy or ablative radiosurgery, the non-invasive external radiation methods accomplish lower therapeutic ratios (Fig. 3), as they are not able to deposit sufficiently high doses in the lesion (> 30 Gy) and to spare simultaneously the surroundings. Advantages of HDR-brachytherapy and Cyberknife decline for lesions > 3 cm. It is to note as limitation that besides logistical and practical advantages of a single session treatment there is no data available showing improved oncological outcome compared to fractionated schedules. Using sophisticated treatment delivery methods, for instance a breath-hold technique, may allow the application of PTV margins < 10 mm which would improve the comparative dosimetry of the non-invasive external radiation methods. Further, we did not compare outcome data of HDR-brachytherapy with RFA and/or other thermoablative procedures.
Various fractionation schedules have been applied using external radiation techniques including the Cyberknife. For intralesional BED > 100 Gy the use of  Table 4 Tumor control with single doses D s = {20.7/β + (α/2β) 2 } 1/2 − (α/2β) depending on α and α/β fractionated stereotactic ablative radiotherapy SABR is common (3). An optimal SABR-schedule of 3 × 16 Gy was suggested [11] yielding BED = 125 Gy (for α/β = 10 Gy). Another analysis recommended 3 × 17 Gy for colorectal cancer metastases with BED = 138 Gy [22]. The possibility of SABR using fractionated schedules is an important asset for the non-invasive radiation techniques [23,24], which is particularly useful for larger lesions. On the other hand, treatment of a few circumscribed lesions in a single session appears attractive as well. The choice between the various options might finally depend on patient preferences.

Conclusions
HDR-brachytherapy is attractive to treat oligotopic liver metastases in a single session, because it is maximally effective and spares nevertheless the liver yielding the highest therapeutic ratios. Repeated treatments are therefore possible in case of recurrent metastases. VMAT and Tomotherapy are of interest for lesions ≥ 3 cm, because liver exposure becomes comparable to HDR-brachytherapy but tolerance and patient comfort are much better. In case of a single dose advantages of Cyberknife are limited, because it requires invasive procedures and long treatment times, but is inferior to HDR-brachytherapy. VMAT and Tomotherapy are comparable with slight advantages of VMAT for small lesions and improving performance of Tomotherapy for larger lesions. However, all external beam radiotherapy approaches, including Cyberknife, have the important advantage to enable any kind of fractionated schemes.