David P. Gierga

COMPARISON OF TARGET REGISTRATION ERRORS FOR MULTIPLE IMAGE-GUIDED TECHNIQUES IN ACCELERATED PARTIAL BREAST IRRADIATION

PURPOSE External beam accelerated partial breast irradiation requires accurate localization of the target volume for each treatment fraction. Using the concept of target registration error (TRE), the performance of several methods of target localization was compared. METHODS AND MATERIALS Twelve patients who underwent external beam accelerated partial breast irradiation were included in this study. TRE was quantified for four methods of image guidance: standard laser-based setup, kilovoltage imaging of the chest wall, kilovoltage imaging of surgically implanted clips, and three-dimensional surface imaging of the breast. The use of a reference surface created from a free-breathing computed tomography scan and a reference surface directly captured with three-dimensional video imaging were compared. The effects of respiratory motion were also considered, and gating was used for 8 of 12 patients. RESULTS The median value of the TRE for the laser, chest wall, and clip alignment was 7.1 mm (n=94), 5.4 mm (n=81), and 2.4 mm (n=93), respectively. The median TRE for gated surface imaging based on the first fraction reference surface was 3.2 mm (n=49), and the TRE for gated surface imaging using the computed tomography-based reference surface was 4.9 mm (n=56). The TRE for nongated surface imaging using the first fraction reference surface was 6.2 mm (n=25). CONCLUSIONS The TRE of surface imaging using a reference surface directly captured with three-dimensional video and the TRE for clip-based setup were within 1 mm. Gated capture is important for surface imaging to reduce the effects of respiratory motion in accelerated partial breast irradiation.

Effects of organ motion on IMRT treatments with segments of few monitor units

Interplay between organ (breathing) motion and leaf motion has been shown in the literature to have a small dosimetric impact for clinical conditions (over a 30 fraction treatment). However, previous studies did not consider the case of treatment beams made up of many few-monitor-unit (MU) segments, where the segment delivery time (1-2 s) is of the order of the breathing period (3-5 s). In this study we assess if breathing compromises the radiotherapy treatment with IMRT segments of low number of MUs. We assess (i) how delivered dose varies, from patient to patient, with the number of MU per segment, (ii) if this delivered dose is identical to the average dose calculated without motion over the path of the motion, and (iii) the impact of the daily variation of the delivered dose as a function of MU per segment. The organ motion was studied along two orthogonal directions, representing the left-right and cranial-caudal directions of organ movement for a patient setup in the supine position. Breathing motion was modeled as sin(x), sin4(x), and sin6(x), based on functions used in the literature to represent organ motion. Measurements were performed with an ionization chamber and films. For a systematic study of motion effects, a MATLAB simulation was written to model organ movement and dose delivery. In the case of a single beam made up of one single segment, the dose delivered to point in a moving target over 30 fractions can vary up to 20% and 10% for segments of 10 MU and 20 MU, respectively. This dose error occurs because the tumor spends most of the time near the edges of the radiation beam. In the case of a single beam made of multiple segments with low MU, we observed 2.4%, 3.3%, and 4.3% differences, respectively, for sin(x), sin4(x), and sin6(x) motion, between delivered dose and motion-averaged dose for points in the penumbra region of the beam and over 30 fractions. In approximately 5-10% of the cases, differences between the motion-averaged dose and the delivered 30-fraction dose could reach 6%, 8% and 10-12%, respectively for sin(x), sin4(x), and sin6(x) motion. To analyze a clinical IMRT beam, two patient plans were randomly selected. For one of the patients, the beams showed a likelihood of up to 25.6% that the delivered dose would deviate from the motion-averaged dose by more than 1%. For the second patient, there was a likelihood of up to 62.8% of delivering a dose that differs by more than 1% from the motion-averaged dose and a likelihood of up to approximately 30% for a 2% dose error. For the entire five-beam IMRT plan, statistical averaging over the beams reduces the overall dose error between the delivered dose and the motion-averaged dose. For both patients there was a likelihood of up to 7.0% and 33.9% that the dose error was greater than 1%, respectively. For one of the patients, there was a 12.6% likelihood of a 2% dose error. Daily intrafraction variation of the delivered dose of more than 10% is non-negligible and can potentially lead to biological effects. We observed [for sin(x), sin4(x), and sin6(x)] that below 10-15 MU leads to large daily variations of the order of 15-35%. Therefore, for small MU segments, non-negligible biological effects can be incurred. We conclude that for most clinical cases the effects may be small because of the use of many beams, it is desirable to avoid low-MU segments when treating moving targets. In addition, dose averaging may not work well for hypo-fractionation, where fewer fractions are used. For hypo-fractionation, PDF modeling of the tumor motion in IMRT optimization may not be adequate.

Moving targets: detection and tracking of internal organ motion for treatment planning and patient set-up.

BACKGROUND AND PURPOSE Clinical target volumes of the thorax and abdomen are typically expanded to account for inter- and intrafractional organ motion. Usually, such expansions are based on clinical experience and planar observations of target motion during simulation. More precise, 4-dimensional motion margins for a specific patient may improve dose coverage of mobile targets and yet limit unnecessarily large field expansions. We are studying approaches to targeting moving tumors throughout the entire treatment process, from treatment planning to beam delivery. MATERIAL AND METHODS Radio-opaque markers were implanted under CT guidance in the liver at the gross tumor periphery. Organ motion during light respiration was volumetrically imaged by 4D Computed Tomography. Marker motion was also acquired by fluoroscopy and compared with 4DCT data. During treatment, daily diagnostic x-ray images were captured at end-exhale and -inhale for patient setup and target localization. RESULTS Based on the time-resolved CT data, target volumes can be designed to account for respiratory motion during treatment. Motion of the tumor as derived from 4DCT was consistent with fluoroscopic motion analysis. Radiographs acquired in the treatment room enabled millimeter-level patient set-up and assessment of target position relative to bony anatomy. Daily positional variations between bony anatomy and implanted markers were observed. CONCLUSIONS Image guided therapy, based on 4DCT imaging, fluoroscopic imaging studies, and daily gated diagonstic energy set-up radiographs is being developed to improve beam delivery precision. Monitoring internal target motion throughout the entire treatment process will ensure adequate dose coverage of the target while sparing the maximum healthy tissue.

Dosimetric impact of motion in free-breathing and gated lung radiotherapy: A 4D Monte Carlo study of intrafraction and interfraction effects

The purpose of this study was to investigate if interfraction and intrafraction motion in free-breathing and gated lung IMRT can lead to systematic dose differences between 3DCT and 4DCT. Dosimetric effects were studied considering the breathing pattern of three patients monitored during the course of their treatment and an in-house developed 4D Monte Carlo framework. Imaging data were taken in free-breathing and in cine mode for both 3D and 4D acquisition. Treatment planning for IMRT delivery was done based on the free-breathing data with the CORVUS (North American Scientific, Chatsworth, CA) planning system. The dose distributions as a function of phase in the breathing cycle were combined using deformable image registration. The study focused on (a) assessing the accuracy of the CORVUS pencil beam algorithm with Monte Carlo dose calculation in the lung, (b) evaluating the dosimetric effect of motion on the individual breathing phases of the respiratory cycle, and (c) assessing intrafraction and interfraction motion effects during free-breathing or gated radiotherapy. The comparison between (a) the planning system and the Monte Carlo system shows that the pencil beam algorithm underestimates the dose in low-density regions, such as lung tissue, and overestimates the dose in high-density regions, such as bone, by 5% or more of the prescribed dose (corresponding to approximately 3-5 Gy for the cases considered). For the patients studied this could have a significant impact on the dose volume histograms for the target structures depending on the margin added to the clinical target volume (CTV) to produce either the planning target (PTV) or internal target volume (ITV). The dose differences between (b) phases in the breathing cycle and the free-breathing case were shown to be negligible for all phases except for the inhale phase, where an underdosage of the tumor by as much as 9.3 Gy relative to the free-breathing was observed. The large difference was due to breathing-induced motion/deformation affecting the soft/lung tissue density and motion of the bone structures (such as the rib cage) in and out of the beam. Intrafraction and interfraction dosimetric differences between (c) free-breathing and gated delivery were found to be small. However, more significant dosimetric differences, of the order of 3%-5%, were observed between the dose calculations based on static CT (3DCT) and the ones based on time-resolved CT (4DCT). These differences are a consequence of the larger contribution of the inhale phase in the 3DCT data than in the 4DCT.

A Voluntary Breath-Hold Treatment Technique for the Left Breast With Unfavorable Cardiac Anatomy Using Surface Imaging

PURPOSE Breath-hold (BH) treatments can be used to reduce cardiac dose for patients with left-sided breast cancer and unfavorable cardiac anatomy. A surface imaging technique was developed for accurate patient setup and reproducible real-time BH positioning. METHODS AND MATERIALS Three-dimensional surface images were obtained for 20 patients. Surface imaging was used to correct the daily setup for each patient. Initial setup data were recorded for 443 fractions and were analyzed to assess random and systematic errors. Real time monitoring was used to verify surface placement during BH. The radiation beam was not turned on if the BH position difference was greater than 5 mm. Real-time surface data were analyzed for 2398 BHs and 363 treatment fractions. The mean and maximum differences were calculated. The percentage of BHs greater than tolerance was calculated. RESULTS The mean shifts for initial patient setup were 2.0 mm, 1.2 mm, and 0.3 mm in the vertical, longitudinal, and lateral directions, respectively. The mean 3-dimensional vector shift was 7.8 mm. Random and systematic errors were less than 4 mm. Real-time surface monitoring data indicated that 22% of the BHs were outside the 5-mm tolerance (range, 7%-41%), and there was a correlation with breast volume. The mean difference between the treated and reference BH positions was 2 mm in each direction. For out-of-tolerance BHs, the average difference in the BH position was 6.3 mm, and the average maximum difference was 8.8 mm. CONCLUSIONS Daily real-time surface imaging ensures accurate and reproducible positioning for BH treatment of left-sided breast cancer patients with unfavorable cardiac anatomy.

Addressing the growing cancer burden in the wake of the AIDS epidemic in Botswana: The BOTSOGO collaborative partnership.

Botswana has experienced a dramatic increase in HIV-related malignancies over the past decade. The BOTSOGO collaboration sought to establish a sustainable partnership with the Botswana oncology community to improve cancer care. This collaboration is anchored by regular tumor boards and on-site visits that have resulted in the introduction of new approaches to treatment and perceived improvements in care, providing a model for partnership between academic oncology centers and high-burden countries with limited resources.

An investigation of the feasibility of gadolinium for neutron capture synovectomy.

Neutron capture synovectomy (NCS) has been proposed as a possible treatment modality for rheumatoid arthritis. Neutron capture synovectomy is a two-part modality, in which a compound containing an isotope with an appreciable thermal neutron capture cross section is injected directly into the joint, followed by irradiation with a neutron beam. Investigations to date for NCS have focused on boron neutron capture synovectomy (BNCS), which utilizes the 10B(n,alpha)7Li nuclear reaction to deliver a highly localized dose to the synovium. This paper examines the feasibility of gadolinium, specifically 157Gd, as an alternative to boron as a neutron capture agent for NCS. This alternative modality is termed Gadolinium Neutron Capture Synovectomy, or GNCS. Monte Carlo simulations have been used to compare 10B and 157Gd as isotopes for accelerator-based NCS. The neutron source used in these calculations was a moderated spectrum from the 9Be(p,n) reaction at a proton energy of 4 MeV. The therapy time to deliver the NCS therapeutic dose of 10000 RBE-cGy, is 27 times longer when 157Gd is used instead of 10B. The skin dose to the treated joint is 33 times larger when 157Gd is used instead of 10B. Furthermore, the impact of using 157Gd instead of 10B was examined in terms of shielded whole-body dose to the patient. The effective dose is 202 mSv for GNCS, compared to 7.6 mSv for BNCS. This is shown to be a result of the longer treatment times required for GNCS; the contribution of the high-energy photons emitted from neutron capture in gadolinium is minimal. Possible explanations as to the relative performance of 157Gd and 10B are discussed, including differences in the RBE and range of boron and gadolinium neutron capture reaction products, and the relative values of the 10B and 157Gd thermal neutron capture cross section as a function of neutron energy.

Establishing and Delivering Quality Radiation Therapy in Resource-Constrained Settings: The Story of Botswana

There is a global cancer crisis, and it is disproportionately affecting resource-constrained settings, especially in low- and middle-income countries (LMICs). Radiotherapy is a critical and cost-effective component of a comprehensive cancer control plan that offers the potential for cure, control, and palliation of disease in greater than 50% of patients with cancer. Globally, LMICs do not have adequate access to quality radiation therapy and this gap is particularly pronounced in sub-Saharan Africa. Although there are numerous challenges in implementing a radiation therapy program in a low-resource setting, providing more equitable global access to radiotherapy is a responsibility and investment worth prioritizing. We outline a systems approach and a series of key questions to direct strategy toward establishing quality radiation services in LMICs, and highlight the story of private-public investment in Botswana from the late 1990s to the present. After assessing the need and defining the value of radiation, we explore core investments required, barriers that need to be overcome, and assets that can be leveraged to establish a radiation program. Considerations addressed include infrastructure; machine choice; quality assurance and patient safety; acquisition, development, and retention of human capital; governmental engagement; public-private partnerships; international collaborations; and the need to critically evaluate the program to foster further growth and sustainability.

Accuracy in breast shape alignment with 3D surface fitting algorithms

Surface imaging is in use in radiotherapy clinical practice for patient setup optimization and monitoring. Breast alignment is accomplished by searching for a tentative spatial correspondence between the reference and daily surface shape models. In this study, the authors quantify whole breast shape alignment by relying on texture features digitized on 3D surface models. Texture feature localization was validated through repeated measurements in a silicone breast phantom, mounted on a high precision mechanical stage. Clinical investigations on breast shape alignment included 133 fractions in 18 patients treated with accelerated partial breast irradiation. The breast shape was detected with a 3D video based surface imaging system so that breathing was compensated. An in-house algorithm for breast alignment, based on surface fitting constrained by nipple matching (constrained surface fitting), was applied. Results were compared with a commercial software where no constraints are utilized (unconstrained surface fitting). Texture feature localization was validated within 2 mm in each anatomical direction. Clinical data show that unconstrained surface fitting achieves adequate accuracy in most cases, though nipple mismatch is considerably higher than residual surface distances (3.9 mm vs 0.6 mm on average). Outliers beyond 1 cm can be experienced as the result of a degenerate surface fit, where unconstrained surface fitting is not sufficient to establish spatial correspondence. In the constrained surface fitting algorithm, average surface mismatch within 1 mm was obtained when nipple position was forced to match in the [1.5; 5] mm range. In conclusion, optimal results can be obtained by trading off the desired overall surface congruence vs matching of selected landmarks (constraint). Constrained surface fitting is put forward to represent an improvement in setup accuracy for those applications where whole breast positional reproducibility is an issue.

Development and construction of a neutron beam line for accelerator-based boron neutron capture synovectomy

A potential application of the 10B(n, alpha)7Li nuclear reaction for the treatment of rheumatoid arthritis, termed Boron Neutron Capture Synovectomy (BNCS), is under investigation. In an arthritic joint, the synovial lining becomes inflamed and is a source of great pain and discomfort for the afflicted patient. The goal of BNCS is to ablate the synovium, thereby eliminating the symptoms of the arthritis. A BNCS treatment would consist of an intra-articular injection of boron followed by neutron irradiation of the joint. Monte Carlo radiation transport calculations have been used to develop an accelerator-based epithermal neutron beam line for BNCS treatments. The model includes a moderator/reflector assembly, neutron producing target, target cooling system, and arthritic joint phantom. Single and parallel opposed beam irradiations have been modeled for the human knee, human finger, and rabbit knee joints. Additional reflectors, placed to the side and back of the joint, have been added to the model and have been shown to improve treatment times and skin doses by about a factor of 2. Several neutron-producing charged particle reactions have been examined for BNCS, including the 9Be(p,n) reaction at proton energies of 4 and 3.7 MeV, the 9Be(d,n) reaction at deuteron energies of 1.5 and 2.6 MeV, and the 7Li(p,n) reaction at a proton energy of 2.5 MeV. For an accelerator beam current of 1 mA and synovial boron uptake of 1000 ppm, the time to deliver a therapy dose of 10,000 RBEcGy ranges from 3 to 48 min, depending on the treated joint and the neutron producing charged particle reaction. The whole-body effective dose that a human would incur during a knee treatment has been estimated to be 3.6 rem or 0.75 rem, for 1000 ppm or 19,000 ppm synovial boron uptake, respectively, although the shielding configuration has not yet been optimized. The Monte Carlo design process culminated in the construction, installation, and testing of a dedicated BNCS beam line on the high-current tandem electrostatic accelerator at the Laboratory for Accelerator Beam Applications at the Massachusetts Institute of Technology.