Proton Treatment Techniques for Posterior Fossa Tumors: Consequences for Linear Energy Transfer and Dose-Volume Parameters for the Brainstem and Organs at Risk.

PURPOSE: In proton therapy of posterior fossa tumors, at least partial inclusion of the brainstem in the target is necessary because of its proximity to the tumor and required margins. Additionally, the preferred beam geometry results in directing the field distal edge toward this critical structure, raising concerns for brainstem toxicity. Some treatment techniques place the beam's distal edge within the brainstem (dose-sparing techniques), and others avoid elevated linear energy transfer (LET) of the proton field by placing the distal edge beyond it (LET-sparing techniques). Hybrid approaches are also being used. We examine the dosimetric efficacy of these techniques, accounting for LET-dependent and dose-dependent variable relative biologic effectiveness (RBE) distributions. METHODS: Six techniques were applied in ependymoma cases: (a) 3-field dose-sparing; (b) 3-field LET-sparing; (c) 2-field dose-sparing, wide angles; (d) 2-field LET-sparing, wide angles; (e) 2-field LET-sparing, steep angles; and (f) 2-field LET-sparing with feathered distal end. Monte Carlo calculated dose, LET, and RBE-weighted dose distributions were compared. RESULTS: Decreased LET values in the brainstem by LET-sparing techniques were accompanied by higher, not statistically significant, median dose: 53.6 Gy(RBE), 53.4 Gy(RBE), and 54.3 Gy(RBE) for techniques (b), (d), and (e) versus 52.1 Gy(RBE) for technique (a). Accounting for variable RBE distributions, the brainstem volume receiving at least 55 Gy(RBE) increased from 72.5% for technique (a) to 80.3% for (b) (P<.01) and from 70.7% for technique (c) to 77.6% for (d) (P<.01). Less than 2%, but statistically significant, decrease in maximum variable RBE-weighted brainstem dose was observed for the LET-sparing techniques compared with the corresponding dose-sparing (P=.03 and .004). CONCLUSIONS: Extending the proton range beyond the brainstem to reduce LET results in clinically comparable maximum radiobiologic effective dose to this sensitive structure. However this method significantly increasing the brainstem volume receiving RBE-weighted dose higher than 55 Gy(RBE) with possible consequences based on known dose-volume parameters for increased toxicity.

Development of a Monte Carlo model for treatment planning dose verification of the Leksell Gamma Knife Perfexion radiosurgery system.

Detailed Monte Carlo (MC) modeling of the Leksell Gamma Knife (GK) Perfexion (PFX) collimator system is the only accurate ab initio approach appearing in the literature. As a different approach, in this work, we present a MC model based on film measurement. By adjusting the model parameters and fine-tuning the derived fluence map for each individual source to match the manufacturer's ring output factors, we created a reasonable virtual source model for MC simulations to verify treatment planning dose for the GK PFX radiosurgery system. The MC simulation model was commissioned by simple single shots. Dose profiles and both ring and collimator output factors were compared with the treatment planning system (TPS). Good agreement was achieved for dose profiles especially for the region of plateau (< 2%), while larger difference (< 5%) came from the penumbra region. The maximum difference of the calculated output factor was within 0.7%. The model was further validated by a clinical test case. Good agreement was obtained. The DVHs for brainstem and the skull were almost identical and, for the target, the volume covered by the prescription (12.5 Gy to 50% isodose line) was 95.6% from MC calculation versus 100% from the TPS.

Incidence of CNS Injury for a Cohort of 111 Patients Treated With Proton Therapy for Medulloblastoma: LET and RBE Associations for Areas of Injury.

BACKGROUND: Central nervous system (CNS) injury is a rare complication of radiation therapy for pediatric brain tumors, but its incidence with proton radiation therapy (PRT) is less well defined. Increased linear energy transfer (LET) and relative biological effectiveness (RBE) at the distal end of proton beams may influence this risk. We report the incidence of CNS injury in medulloblastoma patients treated with PRT and investigate correlations with LET and RBE values. METHODS AND MATERIALS: We reviewed 111 consecutive patients treated with PRT for medulloblastoma between 2002 and 2011 and selected patients with clinical symptoms of CNS injury. Magnetic resonance imaging (MRI) findings for all patients were contoured on original planning scans (treatment change areas [TCA]). Dose and LET distributions were calculated for the treated plans using Monte Carlo system. RBE values were estimated based on LET-based published models. RESULTS: At a median follow-up of 4.2 years, the 5-year cumulative incidence of CNS injury was 3.6% for any grade and 2.7% for grade 3+. Three of 4 symptomatic patients were treated with a whole posterior fossa boost. Eight of 10 defined TCAs had higher LET values than the target but statistically nonsignificant differences in RBE values (P=.12). CONCLUSIONS: Central nervous system and brainstem injury incidence for PRT in this series is similar to that reported for photon radiation therapy. The risk of CNS injury was higher for whole posterior fossa boost than for involved field. Although no clear correlation with RBE values was found, numbers were small and additional investigation is warranted to better determine the relationship between injury and LET.

Disregarding RBE variation in treatment plan comparison may lead to bias in favor of proton plans.

PURPOSE: Currently in proton radiation therapy, a constant relative biological effectiveness (RBE) equal to 1.1 is assumed. The purpose of this study is to evaluate the impact of disregarding variations in RBE on the comparison of proton and photon treatment plans. METHODS: Intensity modulated treatment plans using photons and protons were created for three brain tumor cases with the target situated close to organs at risk. The proton plans were optimized assuming a standard RBE equal to 1.1, and the resulting linear energy transfer (LET) distribution for the plans was calculated. In the plan evaluation, the effect of a variable RBE was studied. The RBE model used considers the RBE variation with dose, LET, and the tissue specific parameter α/ß of photons. The plan comparison was based on dose distributions, DVHs and normal tissue complication probabilities (NTCPs). RESULTS: Under the assumption of RBE=1.1, higher doses to the tumor and lower doses to the normal tissues were obtained for the proton plans compared to the photon plans. In contrast, when accounting for RBE variations, the comparison showed lower doses to the tumor and hot spots in organs at risk in the proton plans. These hot spots resulted in higher estimated NTCPs in the proton plans compared to the photon plans. CONCLUSIONS: Disregarding RBE variations might lead to suboptimal proton plans giving lower effect in the tumor and higher effect in normal tissues than expected. For cases where the target is situated close to structures sensitive to hot spot doses, this trend may lead to bias in favor of proton plans in treatment plan comparisons.

Statistical assessment of proton treatment plans under setup and range uncertainties.

PURPOSE: To evaluate a method for quantifying the effect of setup errors and range uncertainties on dose distribution and dose-volume histogram using statistical parameters; and to assess existing planning practice in selected treatment sites under setup and range uncertainties. METHODS AND MATERIALS: Twenty passively scattered proton lung cancer plans, 10 prostate, and 1 brain cancer scanning-beam proton plan(s) were analyzed. To account for the dose under uncertainties, we performed a comprehensive simulation in which the dose was recalculated 600 times per given plan under the influence of random and systematic setup errors and proton range errors. On the basis of simulation results, we determined the probability of dose variations and calculated the expected values and standard deviations of dose-volume histograms. The uncertainties in dose were spatially visualized on the planning CT as a probability map of failure to target coverage or overdose of critical structures. RESULTS: The expected value of target coverage under the uncertainties was consistently lower than that of the nominal value determined from the clinical target volume coverage without setup error or range uncertainty, with a mean difference of -1.1% (-0.9% for breath-hold), -0.3%, and -2.2% for lung, prostate, and a brain cases, respectively. The organs with most sensitive dose under uncertainties were esophagus and spinal cord for lung, rectum for prostate, and brain stem for brain cancer. CONCLUSIONS: A clinically feasible robustness plan analysis tool based on direct dose calculation and statistical simulation has been developed. Both the expectation value and standard deviation are useful to evaluate the impact of uncertainties. The existing proton beam planning method used in this institution seems to be adequate in terms of target coverage. However, structures that are small in volume or located near the target area showed greater sensitivity to uncertainties.

Patient-specific three-dimensional concomitant dose from cone beam computed tomography exposure in image-guided radiotherapy.

PURPOSE: The purpose of the present study was to quantify the concomitant dose received by patients undergoing cone beam computed tomography (CBCT) scanning in different clinical scenarios as a part of image-guided radiotherapy (IGRT) procedures. METHODS AND MATERIALS: We calculated the three-dimensional concomitant dose received as a result of CBCT scans in 6 patients representing different clinical scenarios: two pelvis, two head and neck, and two chest. We assessed the effect that a daily on-line IGRT strategy would have on the patient dose distribution, assuming 40 CBCT scans throughout the treatment course. The additional dose to the planning target volume margin region was also estimated. RESULTS: In the pelvis, a single CBCT scan delivered a mean dose to the femoral heads of 2-6 cGy and the rectum of 1-2 cGy. An additional dose to the planning target volume was within 1-3 cGy. In the chest, the mean dose to the planning target volume varied from 2.5 to 5 cGy. The lung and spinal cord planning organ at risk volume received ≤4 cGy and ≤5 cGy, respectively. In the head and neck, a single CBCT scan delivered a mean dose of 0.3 cGy, with bony structures receiving 0.5-0.8 cGy. The femoral heads received an additional dose of 1.5-2.5 Gy. A reduction of 20-30% in the mean dose to the organs at risk was achieved using bowtie filtration. In the head and neck, the dose to the eyes and brainstem was eliminated by decreasing the craniocaudal field size. CONCLUSIONS: The additional dose from on-line IGRT procedures can be clinically relevant. The organ dose can be significantly reduced with the use of appropriate patient-specific settings. The concomitant dose from CBCT should be accounted for and the acquisition settings optimized for optimal IGRT strategies on a patient basis.

A pre-clinical assessment of an atlas-based automatic segmentation tool for the head and neck.

BACKGROUND AND PURPOSE: Accurate conformal radiotherapy treatment requires manual delineation of target volumes and organs at risk (OAR) that is both time-consuming and subject to large inter-user variability. One solution is atlas-based automatic segmentation (ABAS) where a priori information is used to delineate various organs of interest. The aim of the present study is to establish the accuracy of one such tool for the head and neck (H&N) using two different evaluation methods. MATERIALS AND METHODS: Two radiotherapy centres were provided with an ABAS tool that was used to outline the brainstem, parotids and mandible on several patients. The results were compared to manual delineations for the first centre (EM1) and reviewed/edited for the second centre (EM2), both of which were deemed as equally valid gold standards. The contours were compared in terms of their volume, sensitivity and specificity with the results being interpreted using the Dice similarity coefficient and a receiver operator characteristic (ROC) curve. RESULTS: Automatic segmentation took typically approximately 7min for each patient on a standard PC. The results indicated that the atlas contour volume was generally within +/-1SD of each gold standard apart from the parotids for EM1 and brainstem for EM2 that were over- and under-estimated, respectively (within +/-2SD). The similarity of the atlas contours with their respective gold standard was satisfactory with an average Dice coefficient for all OAR of 0.68+/-0.25 for EM1 and 0.82+/-0.13 for EM2. All data had satisfactory sensitivity and specificity resulting in a favourable position in ROC space. CONCLUSIONS: These tests have shown that the ABAS tool exhibits satisfactory sensitivity and specificity for the OAR investigated. There is, however, a systematic over-segmentation of the parotids (EM1) and under-segmentation of the brainstem (EM2) that require careful review and editing in the majority of cases. Such issues have been discussed with the software manufacturer and a revised version is due for release.

Mapping radiation dose distribution on the fractional anisotropy map: applications in the assessment of treatment-induced white matter injury.

We describe a method to map whole brain radiation dose distribution on to diffusion tensor MR (DT-MR) fractional anisotropy (FA) images and illustrate its applications for studying dose-effect relationships and regional susceptibility in two childhood medulloblastoma survivors. To determine the FA changes voxel-by-voxel in white matter, the post-treatment follow-up FA maps were coregistered to baseline pre-treatment FA maps and automatic segmentation for white matter was carried out. DeltaFA maps representing relative FA change in white matter were hence generated for visual inspection and quantitative analysis. The radiation dose distribution, calculated from radiotherapy plan and exported as images, was coregistered to baseline FA images. DT-MR imaging and processing noise was small with root mean square value of 1.49% for mean DeltaFA. We evaluated the mean DeltaFA changes of regions-of-interest according to radiation dose regions to provide an estimate of the dose-response and found increasing reduction in mean DeltaFA with increasing radiation dose up to 45 Gy after which there was a reversal in the mean FA trend and mean FA approached baseline value. We also found more severe mean FA reduction in the frontal lobes compared to the parietal lobes despite the same radiation dose, suggesting regional susceptibility in the frontal lobe, and mean FA increase in the brainstem after radiation in both patients. We conclude that the method described may be useful in estimating dose-effect relationships and studying regional susceptibility of the brain to radiation in medulloblastoma survivors.

Monte Carlo evaluation of tissue inhomogeneity effects in the treatment of the head and neck.

PURPOSE: To use Monte Carlo dose calculation to assess the degree to which tissue inhomogeneities in the head and neck affect static field conformal, computed tomography (CT)-based 6-MV photon treatment plans. METHODS AND MATERIALS: We retrospectively studied the three-dimensional treatment plans that had been used for the treatment of 5 patients with tumors in the nasopharyngeal or paranasal sinus regions. Two patients had large surgical cavities. The plans were designed with a clinical treatment planning system that uses a measurement-based pencil-beam dose-calculation algorithm with an equivalent path-length inhomogeneity correction. Each plan employs conformally-shaped 6-MV photon beams. Patient anatomy and electron densities were obtained from the treatment planning CT images. For each plan, the dose distribution was recalculated with the Monte Carlo method, utilizing the same beam geometry and CT images. The Monte Carlo method accurately accounts for the perturbation effects of local tissue heterogeneities. The Monte Carlo calculated dose distributions were compared with those from the clinical treatment planning system. RESULTS: The degree to which tissue inhomogeneity affects the dose distributions of individual fields varies with the specific anatomic geometry, especially the size and location of air cavities in relation to the beam orientation and field size. Most of the beam apertures completely enclose the air cavities within or adjacent to the gross tumor volume (GTV). Equivalent squares (including blocking) ranged from approximately 5 to 9.5 cm. A common feature observed for individual fields is that the Monte Carlo calculated doses to tissue directly behind and within an air cavity are lower. However, after combining the fields employed in each treatment plan, the overall dose distribution shows only small differences between the two methods. For all 5 patients, the Monte Carlo calculated treatment plans showed a slightly lower dose received by the 95% of target volume (D(95)) than the plans calculated with the pencil-beam algorithm. The average difference in the target volume encompassed by the prescription isodose line was less than 2.2%. The difference between the dose-volume histograms (DVHs) of the GTV was generally small. For the brainstem and chiasm, the DVHs of the two plans were similar. For the spinal cord, differences in the details of the DHV and the dose to 1 cc (D(1cc)) of the structure were observed, with Monte Carlo calculation generally predicting increased dose indices to the spinal cord. However, these changes are not expected to be clinically significant. CONCLUSION: For 6-MV photons, the effects of both normal tissue inhomogeneities and surgical air cavities on the target coverage were adequately accounted for by conventional pencil beam methods for all of the cases studied. Although differences in details of the DVHs of the normal structures were observed, depending on whether Monte Carlo or pencil-beam algorithm was used for calculation, these differences are not expected to be clinically significant. In general, the pencil-beam calculation corrected for primary attenuation by the equivalent pathlength is a sufficiently accurate method for head-and-neck treatment planning using 6-MV photons.

[The otoneurological assessment of radiotherapy-induced brain stem involvement in patients with rhinopharyngeal neoplasms]. / Valutazione otoneurologica della sofferenza troncoencefalica radioindotta in pazienti con neoplasie della rinofaringe.

The authors report a study in which otoneurological tests were employed in order to determine the possibility of radio-induced alterations in the brainstem of patients with a rhinopharyngeal carcinoma, which extended to basicranial structures, who had undergone radiotherapy. The case report includes 16 patients; 10 males and 6 females, aged from 37 to 82, all with rhinopharyngeal tumors. All the subjects underwent Co 60 radiotherapy (44-68 Gy); the brainstem received from 40 to 100% of the total dose. Otoneurological evaluation prior to and following radiotherapy was performed employing pure tone audiometry, ABR, rotatory tests, saccadic eye movements, smooth pursuit. After RT treatment, ABR analysis revealed an abnormal wave I-V interpeak interval in 40% of the cases and pathologic in 37%. Smooth pursuit, saccades and sinusoidal rotation analysis showed important alterations respectively in 21%, 6% and 12% of the subjects. The most significant variations were in patients who received more than 60 Gy. The data gathered regarding abnormalities of otoneurological parameters indicate a probable close relationship between these modifications and precocious radio-induced brainstem damage.