Cluster model incorporating heterogeneous dose distribution of partial parotid irradiation for radiotherapy induced xerostomia prediction with machine learning methods.

PURPOSE: A cluster model incorporating heterogeneous dose distribution within the parotid gland was developed and validated retrospectively for radiotherapy (RT) induced xerostomia prediction with machine learning (ML) techniques. METHODS: Sixty clusters were obtained at 1 Gy step size with threshold doses ranging from 1 to 60 Gy, for each of the enrolled 155 patients with HNC from three institutions. Feature clusters were selected with the neighborhood component analysis (NCA) and subsequently fed into four supervised ML models for xerostomia prediction comparison: support vector machines (SVM), k-nearest neighbor (kNN), naïve Bayes (NB), and random forest (RF). The predictive performance of each model was evaluated using cross validation resampling with the area-under-the-curves (AUC) of the receiver-operating-characteristic (ROC). The xerostomia predicting capacity using testing data was assessed with accuracy, sensitivity, and specificity for these models and three cluster connectivity choices. Mean dose based logistic regression served as the benchmark for evaluation. RESULTS: Feature clusters identified by NCA fell in three threshold dose ranges: 5-15Gy, 25-35Gy, and 45-50Gy. Mean dose predictive power was 15% lower than that of the cluster model using the logistic regression classifier. Model validation demonstrated that kNN model outperformed slightly other three models but no substantial difference was observed. Applying the fine-tuned models to testing data yielded that the mean accuracy from SVM, kNN and NB models were between 0.68 and 0.7 while that of RF was ∼0.6. SVM model yielded the best sensitivity (0.76) and kNN model delivered consistent sensitivity and specificity. This is consistent with cross validation. Clusters calculated with three connectivity choices exhibited minimally different predictions. CONCLUSION: Compared to mean dose, the proposed cluster model has shown its improvement as the xerostomia predictor. When combining with ML techniques, it could provide a clinically useful tool for xerostomia prediction and facilitate decision making during radiotherapy planning for patients with HNC.

Medical Physics Practice Guideline (MPPG) 11.a: Plan and chart review in external beam radiotherapy and brachytherapy.

A therapeutic medical physicist is responsible for reviewing radiation therapy treatment plans and patient charts, including initial treatment plans and new chart review, on treatment chart (weekly) review, and end of treatment chart review for both external beam radiation and brachytherapy. Task group report TG 275 examined this topic using a risk-based approach to provide a thorough analysis and guidance for best practice. Considering differences in resources and workflows of various clinical practice settings, the Professional Council of the American Association of Physicists in Medicine assembled this task group to develop a practice guideline on the same topic to provide a minimum standard that balances an appropriate level of safety and resource utilization. This medical physics practice guidelines (MPPG) thus provides a concise set of recommendations for medical physicists and other clinical staff regarding the review of treatment plans and patient charts while providing specific recommendations about who to be involved, and when/what to check in the chart review process. The recommendations, particularly those related to the initial plan review process, are critical for preventing errors and ensuring smooth clinical workflow. We believe that an effective review process for high-risk items should include multiple layers with collective efforts across the department. Therefore, in this report, we make specific recommendations for various roles beyond medical physicists. The recommendations of this MPPG have been reviewed and endorsed by the American Society of Radiologic Technologists and the American Association of Medical Dosimetrists.

Commissioning a 50-100 kV X-ray unit for skin cancer treatment.

This study provides the authors' experience along with dosimetric data from the commissioning of two Sensus SRT-100 50-100 kV X-ray units. Data collected during the commissioning process included: a) HVL, b) output (dose rate), c) applicator cone factors, and d) percentage depth dose. A Farmer-type chamber (PTW-N23333), and a thin-window parallel plate ion chamber (PTW-N23342) were used for dose rate measurements and dose profiles were measured with EBT3 GafChromic film. The average HVL values for 50, 70, and 100 kV of the two treatment units were found to be 0.52, 1.15, and 2.20 mm Al, respectively. The HVL's were 5%-9% lower when measured with the Farmer chamber, as compared to measurements with the parallel plate chamber, for energies of 70 and 100 kV. Dose rates were also measured to be 3%-4% lower with the Farmer chamber. The dose rate variation between the two units was found to be 2%-9% for 50, 70, and 100 kV. The dose uniformity over a circle of 2 cm diameter was within 4% in four cardinal directions; however, the dose profiles for the 5 cm applicator were nonuniform, especially in the cathode-anode direction. Measurements indicated as much as 15% lower dose for the 50 kV beam at field edge on the anode side, when normalized to the center. The crossline profile was relatively more symmetric, with a maximum deviation of 10% at the field edge. All ion chamber readings agreed with film measurements within 3%. The nonuniform profile produced by these units may introduce uncertainty in dose rate measurements, especially for larger applicators. Since there is no intrinsic tool (crosshair or field light) for alignment with the beam axis, the user should take care when positioning the chamber for output measurements. The data obtained with a Farmer-type chamber should be used cautiously and as a reference only for the SRT-100 X-ray treatment unit.

Dosimetric comparison of two afterloaders and treatment planning systems in vaginal high-dose-rate brachytherapy.

Purpose: In treatment planning for high-dose-rate (HDR) single-channel vaginal cylinder brachytherapy, dose distribution along the cylinder is influenced by the anisotropy of the source. Differences in anisotropy are due to differences in source dimensions and characteristics. In this study, we compared HDR vaginal cylinder brachytherapy treatment plans from two afterloader/treatment planning systems. Material and methods: Seventy-five plans with prescription to the surface were generated for cylinders in Varian BrachyVision and Elekta Oncentra. To understand the impact of source anisotropy on dose distribution to the surface of the cylinder, potential effect caused by differences in cylinder geometry between systems was eliminated by re-planning Varian cylinder using Elekta source model. Mean relative dose was calculated for each point as well as the dome and length of the cylinder. Related-samples Wilcoxon signed-rank tests were performed to compare the mean relative dose between systems. Results: Treatment plans with VariSource iX source and cylinder demonstrated 16.2% lower (p < 0.001) dose at the tip compared to Elekta v.3. Average dose to the points along the dome of cylinder was 128.4% ±17.9% prescription dose with VariSource iX source and cylinder, and 99.9% ±4.3% with Elekta v.3 source and cylinder. For the same cylinder geometry, the effect of source characteristics produced up to 36.8% difference in dose homogeneity. When cylinder types were planned with the same source, there was no significant difference in dose distribution. Conclusions: This study demonstrates that the effect of source characteristics produced up to 37% difference in dose homogeneity when comparing two afterloader/treatment planning systems, independent of cylinder geometry. This insight on variation in dose surrounding source system is imperative for dosimetry considerations. Depending on the choice of afterloader, the extent of EQD2 for tumor control versus normal tissue toxicity can vary.

RapidPlan for Knowledge-Based Planning of Malignant Pleural Mesothelioma.

PURPOSE: Treatment planning for malignant pleural mesothelioma is a challenging task due to the relatively large size of the target and the need to spare critical organs that overlap with or are within the target volume. We aimed to develop a knowledge-based model using RapidPlan (RP) for patients with 2 intact lungs. METHODS AND MATERIALS: Data from 57 patients treated with volumetric modulated arc therapy were chosen for training the dose estimation model at a single dose level. The prescription dose was 50.4 Gy in 1.8 Gy fractions. The model was validated on 23 new patients by comparing the clinical plan to the RP. Time taken to plan the RP was compared with that for the clinical plan. RESULTS: For similar target coverage and plan inhomogeneity, RP significantly improved the sparing of the contralateral lung, heart, stomach, esophagus, and ipsilateral kidney. On average, the contralateral lung V5 Gy and V10 Gy were reduced by 13.9% (P < .001) and 7.9% (P < .001), respectively. The mean heart dose was reduced by 5 Gy (P < .001) and V30 Gy by 9.1% (P < .001). Mean dose to the stomach and esophagus were both reduced by 5 Gy (P < .001), and the ipsilateral kidney V18 Gy by 4.1% (P < .001). Mean total lung dose was reduced by 0.8 Gy with RP, which enabled an increase in prescription dose by 1 fraction Absolute volume of ipsilateral lung was adequately spared by both techniques, while sparing of all other organs, namely the cord, liver, and bowel, was not compromised with RP. Time taken with RP was 20 minutes, 45 seconds versus at least 4 hours for an experienced treatment planner. CONCLUSIONS: The RP model for malignant pleural mesothelioma showed improved sparing of critical organs with a reduced treatment planning time and increased prescription dose.

Improving Dermoscopic Image Segmentation with Enhanced Convolutional-Deconvolutional Networks.

Automatic skin lesion segmentation on dermoscopic images is an essential step in computer-aided diagnosis of melanoma. However, this task is challenging due to significant variations of lesion appearances across different patients. This challenge is further exacerbated when dealing with a large amount of image data. In this paper, we extended our previous work by developing a deeper network architecture with smaller kernels to enhance its discriminant capacity. In addition, we explicitly included color information from multiple color spaces to facilitate network training and thus to further improve the segmentation performance. We participated and extensively evaluated our method on the ISBI 2017 skin lesion segmentation challenge. By training with the 2000 challenge training images, our method achieved an average Jaccard Index (JA) of 0:765 on the 600 challenge testing images, which ranked itself in the first place among 21 final submissions in the challenge.

Dose cluster model parameterization of the parotid gland in irradiation of head and neck cancer.

To explore the parotid normal tissue complication probability (NTCP) modeling with percolation-based dose clusters for head-and-neck patients receiving concomitant chemotherapy and radiation therapy. Cluster models incorporating the spatial dose distribution in the parotid gland were developed to evaluate the radiation induced complication. Cluster metrics including the mean cluster size (NMCS) and the largest cluster size both normalized by the gland volume (NSLC) were evaluated and scrutinized against the benchmark NTCP. Two fitting strategies to the Lyman-Kutcher-Burman (LKB) model using the maximum likelihood method were devised: the volume parameter n fixed at 1.0 (mean dose model) and unrestricted (full LKB model). The fitted parameters TD50 and m were assessed with the LKB NTCP models with the available xerostomia data. NSLC was a better metric than NMCS with reference to the LKB model and strong correlation (r ~ 0.95) was observed between NTCP and NSLC. The mean dose model returned the parameter TD50 (39.9 Gy) and m (0.4) from the NSLC of threshold dose at around 40 Gy. Drastically different TD50 and m values were obtained from the fittings via the full LKB model, where the threshold dose would be near 27 Gy. Bootstrapping analyses further confirmed the fitting outcomes. Strong correlation with the traditional NTCP models revealed that the cluster model could achieve what NTCP models attain and may offer additional information. Parameterization of the model indicated that the model could have different predictions from current clinical recommendations. Further investigation using toxicity data is under way to validate the cluster model.

Three-dimensional cluster formation and structure in heterogeneous dose distribution of intensity modulated radiation therapy.

PURPOSE: To investigate three-dimensional cluster structure and its correlation to clinical endpoint in heterogeneous dose distributions from intensity modulated radiation therapy. METHODS: Twenty-five clinical plans from twenty-one head and neck (HN) patients were used for a phenomenological study of the cluster structure formed from the dose distributions of organs at risks (OARs) close to the planning target volumes (PTVs). Initially, OAR clusters were searched to examine the pattern consistence among ten HN patients and five clinically similar plans from another HN patient. Second, clusters of the esophagus from another ten HN patients were scrutinized to correlate their sizes to radiobiological parameters. Finally, an extensive Monte Carlo (MC) procedure was implemented to gain deeper insights into the behavioral properties of the cluster formation. RESULTS: Clinical studies showed that OAR clusters had drastic differences despite similar PTV coverage among different patients, and the radiobiological parameters failed to positively correlate with the cluster sizes. MC study demonstrated the inverse relationship between the cluster size and the cluster connectivity, and the nonlinear changes in cluster size with dose thresholds. In addition, the clusters were insensitive to the shape of OARs. CONCLUSION: The results demonstrated that the cluster size could serve as an insightful index of normal tissue damage. The clinical outcome of the same dose-volume might be potentially different.

Automatic Skin Lesion Segmentation Using Deep Fully Convolutional Networks With Jaccard Distance.

Automatic skin lesion segmentation in dermoscopic images is a challenging task due to the low contrast between lesion and the surrounding skin, the irregular and fuzzy lesion borders, the existence of various artifacts, and various imaging acquisition conditions. In this paper, we present a fully automatic method for skin lesion segmentation by leveraging 19-layer deep convolutional neural networks that is trained end-to-end and does not rely on prior knowledge of the data. We propose a set of strategies to ensure effective and efficient learning with limited training data. Furthermore, we design a novel loss function based on Jaccard distance to eliminate the need of sample re-weighting, a typical procedure when using cross entropy as the loss function for image segmentation due to the strong imbalance between the number of foreground and background pixels. We evaluated the effectiveness, efficiency, as well as the generalization capability of the proposed framework on two publicly available databases. One is from ISBI 2016 skin lesion analysis towards melanoma detection challenge, and the other is the PH2 database. Experimental results showed that the proposed method outperformed other state-of-the-art algorithms on these two databases. Our method is general enough and only needs minimum pre- and post-processing, which allows its adoption in a variety of medical image segmentation tasks.

Improving Image Quality of On-Board Cone-Beam CT in Radiation Therapy Using Image Information Provided by Planning Multi-Detector CT: A Phantom Study.

PURPOSE: The aim of this study was to improve the image quality of cone-beam computed tomography (CBCT) mounted on the gantry of a linear accelerator used in radiation therapy based on the image information provided by planning multi-detector CT (MDCT). METHODS: MDCT-based shading correction for CBCT and virtual monochromatic CT (VMCT) synthesized using the dual-energy method were performed. In VMCT, the high-energy data were obtained from CBCT, while the low-energy data were obtained from MDCT. An electron density phantom was used to investigate the efficacy of shading correction and VMCT on improving the target detectability, Hounsfield unit (HU) accuracy and variation, which were quantified by calculating the contrast-to-noise ratio (CNR), the percent difference (%Diff) and the standard deviation of the CT numbers for tissue equivalent background material, respectively. Treatment plan studies for a chest phantom were conducted to investigate the effects of image quality improvement on dose planning. RESULTS: For the electron density phantom, the mean value of CNR was 17.84, 26.78 and 34.31 in CBCT, shading-corrected CBCT and VMCT, respectively. The mean value of %Diff was 152.67%, 11.93% and 7.66% in CBCT, shading-corrected CBCT and VMCT, respectively. The standard deviation within a uniform background of CBCT, shading-corrected CBCT and VMCT was 85, 23 and 15 HU, respectively. With regards to the chest phantom, the monitor unit (MU) difference between the treatment plan calculated using MDCT and those based on CBCT, shading corrected CBCT and VMCT was 6.32%, 1.05% and 0.94%, respectively. CONCLUSIONS: Enhancement of image quality in on-board CBCT can contribute to daily patient setup and adaptive dose delivery, thus enabling higher confidence in patient treatment accuracy in radiation therapy. Based on our results, VMCT has the highest image quality, followed by the shading corrected CBCT and the original CBCT. The research results presented in this study should be able to provide a route to reach a high level of image quality for CBCT imaging in radiation oncology.

A Feasibility Study of Tumor Motion Estimate With Regional Deformable Registration Method for 4-Dimensional Radiation Therapy of Lung Cancer.

This study aims to employ 4-dimensional computed tomography to quantify intrafractional tumor motion for patients with lung cancer to improve target localization in radiation therapy. A multistage regional deformable registration was implemented to calculate the excursion of gross tumor volume (GTV) during a breathing cycle. GTV was initially delineated on 0% phase of 4-dimensional computed tomography manually, and a subregion with 20 mm margin supplemented to GTV was generated with Eclipse treatment planning system (Varian Medical Systems, Palo Alto, California). The structures, together with the 4-dimensional computed tomography set, were exported into an in-house software, with which a 3-stage B-spline deformable registration was carried out to map the subregion and warp GTV contour to other breathing phases. The center of mass of the GTV was computed using the contours, and the tumor motion was appraised as the excursion of the center of mass between 0% phase and other phases. Application of the algorithm to the 10 patients showed that clinically satisfactory outcomes were achievable with a spatial accuracy around 2 mm for GTV contour propagation between adjacent phases and 3 mm between opposite phases. The tumor excursion was determined in the vast range of 1 mm through 1.6 cm, depending on the tumor location and tumor size. Compared to the traditional whole image-based registration, the regional method was found computationally a factor of 5 more efficient. The proposed technique has demonstrated its capability in extracting thoracic tumor motion and should find its application in 4-dimensional radiation therapy in the future to maximally utilize the available spatial-temporal information.

Robust breathing signal extraction from cone beam CT projections based on adaptive and global optimization techniques.

We present a study of extracting respiratory signals from cone beam computed tomography (CBCT) projections within the framework of the Amsterdam Shroud (AS) technique. Acquired prior to the radiotherapy treatment, CBCT projections were preprocessed for contrast enhancement by converting the original intensity images to attenuation images with which the AS image was created. An adaptive robust z-normalization filtering was applied to further augment the weak oscillating structures locally. From the enhanced AS image, the respiratory signal was extracted using a two-step optimization approach to effectively reveal the large-scale regularity of the breathing signals. CBCT projection images from five patients acquired with the Varian Onboard Imager on the Clinac iX System Linear Accelerator (Varian Medical Systems, Palo Alto, CA) were employed to assess the proposed technique. Stable breathing signals can be reliably extracted using the proposed algorithm. Reference waveforms obtained using an air bellows belt (Philips Medical Systems, Cleveland, OH) were exported and compared to those with the AS based signals. The average errors for the enrolled patients between the estimated breath per minute (bpm) and the reference waveform bpm can be as low as -0.07 with the standard deviation 1.58. The new algorithm outperformed the original AS technique for all patients by 8.5% to 30%. The impact of gantry rotation on the breathing signal was assessed with data acquired with a Quasar phantom (Modus Medical Devices Inc., London, Canada) and found to be minimal on the signal frequency. The new technique developed in this work will provide a practical solution to rendering markerless breathing signal using the CBCT projections for thoracic and abdominal patients.

Clinical outcome of vertebral compression fracture after single fraction spine radiosurgery for spinal metastases.

Vertebral compression fracture (VCF) occurs after stereotactic body radiation therapy (SBRT) for spine metastasis. Recently, single fraction radiosurgery (sfSRS) is used more frequently. The aim of this study is to determine the clinical outcome of VCF after sfSRS. Spinal instability neoplastic score (SINS) criteria were used to retrospectively score 143 consecutive vertebral segments in 79 patients treated with SRS. Follow-up MRI, pain, and neurologic assessments obtained every 3-6 months. Pain also scored at 7, 14, and 30 days after sfSRS. Follow up was 16 ± 18 months ±SD, range 3-78. Long-term radiographic control occurred in 94 % of cases. Pain improvement resulted within 7 days in 100 % of cases with severe pain and sustained long-term in 95 %. VCF occurred in 21 % of segments: 30 % were de novo VCF. The overall 1 year fracture free probability (1yFFP) was 76 %. Pre-existing VCF resulted in higher probability to progress: 1yFFP 90 versus 60 %. Symptoms presented in 6 % of cases with de novo VCF and 39 % with progressive. The former were treated with vertebral augmentation (VA), the latter with open surgery. Surgery/VA prior to SRS did not change risk of progressive VCF. Univariate but not multivariate analysis identified histology (colorectal), pre-existing VCF, and pain (severe) as significant predictors of VCF. In conclusion, sfSRS compares favourably to SBRT for radiographic and pain control with similar VCF risk. Patients with pre-existing VCF have a higher probability to progress, become symptomatic, and require surgery. These results may help discussing risk and benefits with patients undergoing sfSRS for spinal metastasis and developing new treatment algorithms.

The use of isodose levels to interpret radiation induced lung injury: a quantitative analysis of computed tomography changes.

BACKGROUND: Patients treated with stereotactic body radiation therapy (SBRT) for lung cancer are often found to have radiation-induced lung injury (RILI) surrounding the treated tumor. We investigated whether treatment isodose levels could predict RILI. METHODS: Thirty-seven lung lesions in 32 patients were treated with SBRT and received post-treatment follow up (FU) computed tomography (CT). Each CT was fused with the original simulation CT and treatment isodose levels were overlaid. The RILI surrounding the treated lesion was contoured. The RILI extension index [fibrosis extension index (FEI)] was defined as the volume of RILI extending outside a given isodose level relative to the total volume of RILI and was expressed as a percentage. RESULTS: Univariate analysis revealed that the planning target volume (PTV) was positively correlated with RILI volume at FU: correlation coefficient (CC) =0.628 and P<0.0001 at 1(st) FU; CE =0.401 and P=0.021 at 2(nd) FU; CE =0.265 and P=0.306 at 3(rd) FU. FEI -40 Gy at 1(st) FU was significantly positively correlated with FEI -40 Gy at subsequent FU's (CC =0.689 and P=6.5×10(-5) comparing 1(st) and 2(nd) FU; 0.901 and P=0.020 comparing 2(nd) and 3(rd) FU. Ninety-six percent of the RILI was found within the 20 Gy isodose line. Sixty-five percent of patients were found to have a decrease in RILI on the second 2(nd) CT. CONCLUSIONS: We have shown that RILI evolves over time and 1(st) CT correlates well with subsequent CTs. Ninety-six percent of the RILI can be found to occur within the 20 Gy isodose lines, which may prove beneficial to radiologists attempting to distinguish recurrence vs. RILI.

Tracking fuzzy borders using geodesic curves with application to liver segmentation on planning CT.

PURPOSE: This work aims to develop a robust and efficient method to track the fuzzy borders between liver and the abutted organs where automatic liver segmentation usually suffers, and to investigate its applications in automatic liver segmentation on noncontrast-enhanced planning computed tomography (CT) images. METHODS: In order to track the fuzzy liver-chestwall and liver-heart borders where oversegmentation is often found, a starting point and an ending point were first identified on the coronal view images; the fuzzy border was then determined as a geodesic curve constructed by minimizing the gradient-weighted path length between these two points near the fuzzy border. The minimization of path length was numerically solved by fast-marching method. The resultant fuzzy borders were incorporated into the authors' automatic segmentation scheme, in which the liver was initially estimated by a patient-specific adaptive thresholding and then refined by a geodesic active contour model. By using planning CT images of 15 liver patients treated with stereotactic body radiation therapy, the liver contours extracted by the proposed computerized scheme were compared with those manually delineated by a radiation oncologist. RESULTS: The proposed automatic liver segmentation method yielded an average Dice similarity coefficient of 0.930 ± 0.015, whereas it was 0.912 ± 0.020 if the fuzzy border tracking was not used. The application of fuzzy border tracking was found to significantly improve the segmentation performance. The mean liver volume obtained by the proposed method was 1727 cm(3), whereas it was 1719 cm(3) for manual-outlined volumes. The computer-generated liver volumes achieved excellent agreement with manual-outlined volumes with correlation coefficient of 0.98. CONCLUSIONS: The proposed method was shown to provide accurate segmentation for liver in the planning CT images where contrast agent is not applied. The authors' results also clearly demonstrated that the application of tracking the fuzzy borders could significantly reduce contour leakage during active contour evolution.

What is the optimal dose for 125I prostate implants? A dose-response analysis of biochemical control, posttreatment prostate biopsies, and long-term urinary symptoms.

PURPOSE: To define the optimal dose for 125I prostate implants by correlating post implant CT dosimetry findings with urinary symptoms, biochemical failure, and posttreatment biopsies. METHODS AND MATERIALS: Patients with T1-T2, Gleason score 2-6 prostate cancer treated with I-125 brachytherapy were analyzed. Group 1 (276 patients) was observed from 18 to 108 months (median, 34 months) and had urinary symptoms prospectively assessed using the International Prostate Symptom Score (IPSS) system. Group 2 (181 patients) observed from 24 to 108 months (median, 44 months) and did not receive hormonal therapy. Implant dose was defined as the D90 (dose delivered to 90% of the prostate on a dose-volume histogram). Patients were analyzed by dose categories: <140 Gy, 140 to <160 Gy, 160 to <180 Gy, and > or =180 Gy. In Group 1, the mean pre- to postimplant IPSS scores were compared in different dose categories by using a matched paired t test. In Group 2, the effect of dose on biochemical control was tested with actuarial methods by using the American Society for Therapeutic Radiology and Oncology definition and on local control with posttreatment biopsies (113 patients). RESULTS: A comparison of pre- with postimplant IPSS revealed no significant changes in scores in the dose groups <180 Gy except for small changes in urgency and bladder emptying in the dose group <140 Gy. In dose group >180 Gy, mean scores changed from 0.5 to 1.0 (p=0.002) for emptying, 0.76 to 1.29 (p=0.004) for weak stream, 0.24 to 0.51 (p=0.009) for straining, 1.55 to 1.82 (p=0.05) for nocturia, and 6.3 to 8.45 (p=0.0009) for the total score. Freedom from biochemical failure (FFBF) at 5 years was 68% for doses <140 Gy, 97% for 140 to <160 Gy, 98% for 160 to <180 Gy, and 95% for > or =180 Gy (p=0.0025). Overall, patients with doses <140 Gy (median follow-up, 66 months) had an FFBF of 68%, compared with 96% for patients with doses > or =140 Gy (median follow-up, 35 months; p=0.0002). Multivariate analysis found dose to be the most significant factor affecting FFBF. Positive biopsies were found in 23% for doses <140 Gy, 21% for 140 to <160 Gy, 10% for 160 to <180 Gy, and 8% for > or =180 Gy. Overall, biopsies were positive in 22% for doses <160 Gy vs. 9% for > or =160 Gy (p=0.05). CONCLUSIONS: Optimal 125I prostate implants should deliver a D90 of 140-180 Gy, on the basis of postimplant dosimetry. Doses of <140 Gy are associated with increased biochemical failure, and doses >180 Gy with a small increase in long-term urinary symptoms.

Prostate gland motion and deformation caused by needle placement during brachytherapy.

PURPOSE: To determine the extent of edge and gland position changes caused by needle insertion in patients undergoing prostate brachytherapy. METHODS AND MATERIALS: Nineteen patients with T1-T3 prostate cancer were implanted with the real-time method by using a two-phase peripheral loading technique. Serial contours of the prostate at 5-mm intervals were acquired by the dose-planning system. All of the peripheral needles were then placed and spaced 5-10 mm apart by using the largest transverse ultrasound image as the reference plane. The position of the probe was relocated at the zero plane, and the difference between the preneedle and postneedle zero plane was recorded as the difference in the z axis. Axial ultrasound images were again acquired. The second set of captured images, which matched in number the first set, was contoured over the previously contoured preneedle images. Prostate gland deformation and displacement were determined by comparing the preneedle contoured image with the images captured after needle placement. Deformation was determined by calculating the differences between the edges of the gland as measured at the major axis of the gland (x and y planes). Displacement was determined by measuring the differences between the center positions of the two contoured structures. Deformation and displacement were determined on each acquired 5-mm image. Differences were compared by student's t test. RESULTS: The mean preneedle prostate volume was 47 ml (range, 21.5-68.7 ml), compared with 48.1 ml (range, 19.4-80.3 ml; p = 0.228) after peripheral needle placement. A median of 16 (range, 12-19) peripheral needles were placed. The median change in the base position of the prostate was 1.5 cm (range of 0 to 3.0 cm; p = 0.0034). The mean x and y deformation was 6.8 mm (median, 7.9 mm; range, 4.3-8.1 mm) and 3.6 mm (median, 3.3 mm; range, 1.0-5.5 mm), respectively. The greatest deformation for any individual slice for x was 21.6 mm and for y was 15.3 mm. The mean number of slices that were found with a >2-, 5-, and 10-mm deformation in the x axis was 7 (range, 3-10), 4 (range, 1-3), and 1 (range, 0-4), respectively. Similar deformation in the y axis was found in 6 (range, 3-10), 2.5 (range, 0-6), and 0.3 (range, 0-2) slices. The mean x and y displacement was 1.9 mm (median, 1.8 mm; range, 0.3-6.6 mm) and 2.8 mm (median, 1.9 mm; range, 2-5.8 mm). The greatest displacement for any individual slice for x was 7 mm and for y was 10 mm. The mean number of slices with a displacement >2, 5, and 10 mm in the x axis was 5 (range, 1-10), 0.8 (range, 0-5), and 0, respectively. Similar displacement in the y axis was found in 5 (range, 0-9), 1.7 (range, 0-7), and 0 slices, respectively. CONCLUSIONS: Placing most needles in the periphery results in a minimal prostate volume increase, suggesting little need to overplan the implant when this method is used. However, significant edge and gland position changes caused by the needle insertion did occur. These changes may explain some of the difficulty in reproducing the preplan and should be taken into consideration for all types of prostate brachytherapy planning.

Comparison of intraoperative dosimetric implant representation with postimplant dosimetry in patients receiving prostate brachytherapy.

PURPOSE: To compare the results of intraoperative dosimetry with those of CT-based postimplant dosimetry in patients undergoing prostate seed implantation. METHODS AND MATERIALS: Seventy-seven patients with T1-T3 prostate cancer received an ultrasound-guided permanent seed implant (36 received (125)I, 7 (103)Pd, and 34 a partial (103)Pd implant plus external beam radiation therapy). The implantation was augmented with an intraoperative dosimetric planning system. After the peripheral needles were placed, 5-mm axial images were acquired into the treatment planning system. Soft tissue structures (prostate, urethra, and rectum) were contoured, and exact needle positions were registered. Seeds were placed with an applicator, and their positions were entered into the planning system. The dose distributions for the implant were calculated after interior needle and seed placement. Postimplant dosimetry was performed 1 month later on the basis of CT imaging. Prostate and urethral doses were compared, by using paired t tests, for the real-time dosimetry in the operating room (OR) and the postimplant dosimetry. RESULTS: The mean preimplant prostate volume was 39.8 cm(3), the postneedle planning volume was 41.5 cm(3) (p<0.001), and the 1-month CT volume was 43.6 cm(3) (p<0.001). The mean difference between the OR dose received by 90% of the prostate (D(90)) and the CT D(90) was 3.4% (95% confidence interval, 2.5-6.6%; p=0.034). The mean dose to 30% of the urethra was 120% of prescription in the OR and 138% on CT. The mean difference was 18% (95% confidence interval, 13-24%; p<0.001). CONCLUSIONS: Although small differences exist between the OR and CT dosimetry results, these data suggest that this intraoperative implant dosimetric representation system provides a close match to the actual delivered doses. These data support the use of this system to modify the implant during surgery to achieve more consistent dosimetry results.

Image-guided stereotactic body radiotherapy for lung tumors using BodyLoc with tomotherapy: clinical implementation and set-up accuracy.

We investigated the use of a BodyLoc immobilization and stereotactic localization device combined with TomoTherapy megavoltage CT (MVCT) in lung stereotactic body radiotherapy (SBRT) to reduce set-up uncertainty and treatment time. Eight patients treated with 3-5 fractions of SBRT were retrospectively analyzed. A BodyLoc localizer was used in both CT simulation for localization and the initial patient treatment set-up. Patients were immobilized with a vacuum cushion on the back and a thermoplastic body cast on the anterior body. Pretreatment MVCT from the TomoTherapy unit was fused with the planning kilovoltage CT (KVCT) before each fraction of treatment to determine interfractional set-up error. The comparison of two MVCTs during a fraction of treatment resulted in the intrafractional uncertainty of the treatment. A total of 224 target isocenter shifts were analyzed to assess these inter- and intrafractional tumor motions. We found that for interfractional shifts, the mean set-up errors and standard deviations were -1.1 +/- 2.8 mm, -2.5 +/- 8.7 mm, and 4.1 +/- 2.6 mm, for lateral, longitudinal, and vertical variation, respectively; the mean setup rotational variation was -0.3 +/- 0.7 degrees; and the maximum motion was 13.5 mm in the longitudinal direction. For intrafractional shifts, the mean set-up errors and standard deviations were -0.1 +/- 0.7 mm, -0.3 +/- 2.0 mm, and 0.5 +/- 1.1 mm for the lateral, longitudinal, and vertical shifts, respectively; the mean rotational variation was 0.1 +/- 0.2 degrees; and the maximum motion was 3.8 mm in the longitudinal direction. There was no correlation among patient characteristics, set-up uncertainties, and isocenter shifts, and the interfractional set-up uncertainties were larger than the intrafractional isocenter shift. The results of this study suggested that image-guided stereotactic body radiotherapy using the BodyLoc immobilization system with TomoTherapy can improve treatment accuracy.

Analysis of daily setup variation with tomotherapy megavoltage computed tomography.

The purpose of this study was to evaluate different setup uncertainties for various anatomic sites with TomoTherapy pretreatment megavoltage computed tomography (MVCT) and to provide optimal margin guidelines for these anatomic sites. Ninety-two patients with tumors in head and neck (HN), brain, lung, abdominal, or prostate regions were included in the study. MVCT was used to verify patient position and tumor target localization before each treatment. With the anatomy registration tool, MVCT provided real-time tumor shift coordinates relative to the positions where the simulation CT was performed. Thermoplastic facemasks were used for HN and brain treatments. Vac-Lok cushions were used to immobilize the lower extremities up to the thighs for prostate patients. No respiration suppression was administered for lung and abdomen patients. The interfractional setup variations were recorded and corrected before treatment. The mean interfractional setup error was the smallest for HN among the 5 sites analyzed. The average 3D displacement in lateral, longitudinal, and vertical directions for the 5 sites ranged from 2.2-7.7 mm for HN and lung, respectively. The largest movement in the lung was 2.0 cm in the longitudinal direction, with a mean error of 6.0 mm and standard deviation of 4.8 mm. The mean interfractional rotation variation was small and ranged from 0.2-0.5 degrees, with the standard deviation ranging from 0.7-0.9 degrees. Internal organ displacement was also investigated with a posttreatment MVCT scan for HN, lung, abdomen, and prostate patients. The maximum 3D intrafractional displacement across all sites was less than 4.5 mm. The interfractional systematic errors and random errors were analyzed and the suggested margins for HN, brain, prostate, abdomen, and lung in the lateral, longitudinal, and vertical directions were between 4.2 and 8.2 mm, 5.0 mm and 12.0 mm, and 1.5 mm and 6.8 mm, respectively. We suggest that TomoTherapy pretreatment MVCT can be used to improve the accuracy of patient positioning and reduce tumor margin.