A Monte Carlo model for organ dose reconstruction of patients in pencil beam scanning (PBS) proton therapy for epidemiologic studies of late effects.

Significant efforts such as the Pediatric Proton/Photon Consortium Registry (PPCR) involving multiple proton therapy centers have been made to conduct collaborative studies evaluating outcomes following proton therapy. As a groundwork dosimetry effort for the late effect investigation, we developed a Monte Carlo (MC) model of proton pencil beam scanning (PBS) to estimate organ/tissue doses of pediatric patients at the Maryland Proton Treatment Center (MPTC), one of the proton centers involved in the PPCR. The MC beam modeling was performed using the TOPAS (TOol for PArticle Simulation) MC code and commissioned to match measurement data within 1% for range, and 0.3 mm for spot sizes. The established MC model was then tested by calculating organ/tissue doses for sample intracranial and craniospinal irradiations on whole-body pediatric computational human phantoms. The simulated dose distributions were compared with the treatment planning system dose distributions, showing the 3 mm/3% gamma index passing rates of 94%-99%, validating our simulations with the MC model. The calculated organ/tissue doses per prescribed doses for the craniospinal irradiations (1 mGy Gy-1 to 1 Gy Gy-1) were generally much higher than those for the intracranial irradiations (2.1 µGy Gy-1 to 0.1 Gy Gy-1), which is due to the larger field coverage of the craniospinal irradiations. The largest difference was observed at the adrenal dose, i.e. ∼3000 times. In addition, the calculated organ/tissue doses were compared with those calculated with a simplified MC model, showing that the beam properties (i.e. spot size, spot divergence, mean energy, and energy spread) do not significantly influence dose calculations despite the limited irradiation cases. This implies that the use of the MC model commissioned to the MPTC measurement data might be dosimetrically acceptable for patient dose reconstructions at other proton centers particularly when their measurement data are unavailable. The developed MC model will be used to reconstruct organ/tissue doses for MPTC pediatric patients collected in the PPCR.

Comparison of normal tissue dose calculation methods for epidemiological studies of radiotherapy patients.

Radiation dosimetry is an essential input for epidemiological studies of radiotherapy patients aimed at quantifying the dose-response relationship of late-term morbidity and mortality. Individualised organ dose must be estimated for all tissues of interest located in-field, near-field, or out-of-field. Whereas conventional measurement approaches are limited to points in water or anthropomorphic phantoms, computational approaches using patient images or human phantoms offer greater flexibility and can provide more detailed three-dimensional dose information. In the current study, we systematically compared four different dose calculation algorithms so that dosimetrists and epidemiologists can better understand the advantages and limitations of the various approaches at their disposal. The four dose calculations algorithms considered were as follows: the (1) Analytical Anisotropic Algorithm (AAA) and (2) Acuros XB algorithm (Acuros XB), as implemented in the Eclipse treatment planning system (TPS); (3) a Monte Carlo radiation transport code, EGSnrc; and (4) an accelerated Monte Carlo code, the x-ray Voxel Monte Carlo (XVMC). The four algorithms were compared in terms of their accuracy and appropriateness in the context of dose reconstruction for epidemiological investigations. Accuracy in peripheral dose was evaluated first by benchmarking the calculated dose profiles against measurements in a homogeneous water phantom. Additional simulations in a heterogeneous cylinder phantom evaluated the performance of the algorithms in the presence of tissue heterogeneity. In general, we found that the algorithms contained within the commercial TPS (AAA and Acuros XB) were fast and accurate in-field or near-field, but not acceptable out-of-field. Therefore, the TPS is best suited for epidemiological studies involving large cohorts and where the organs of interest are located in-field or partially in-field. The EGSnrc and XVMC codes showed excellent agreement with measurements both in-field and out-of-field. The EGSnrc code was the most accurate dosimetry approach, but was too slow to be used for large-scale epidemiological cohorts. The XVMC code showed similar accuracy to EGSnrc, but was significantly faster, and thus epidemiological applications seem feasible, especially when the organs of interest reside far away from the field edge.

Revisiting fetal dose during radiation therapy: evaluating treatment techniques and a custom shield.

To create a comprehensive dataset of peripheral dose (PD) measurements from a new generation of linear accelerators with and without the presence of a newly designed fetal shield, PD measurements were performed to evaluate the effects of depth, field size, distance from the field edge, collimator angle, and beam modi-fiers for common treatment protocols and modalities. A custom fetal lead shield was designed and made for our department that allows external beam treatments from multiple angles while minimizing the need to adjust the shield during patient treatments. PD measurements were acquired for a comprehensive series of static fields on a stack of Solid Water. Additionally, PDs from various clinically relevant treatment scenarios for pregnant patients were measured using an anthropomorphic phantom that was abutted to a stack of Solid Water. As expected, the PD decreased as the distance from the field edge increased and the field size decreased. On aver-age, a PD reduction was observed when a 90° collimator rotation was applied and/or when the tertiary MLCs and jaws defined the field aperture. However, the effect of the collimator rotation (90° versus 0°) in PD reduction was not found to be clini-cally significant when the tertiary MLCs were used to define the field aperture. In the presence of both the MLCs and the fetal shield, the PD was reduced by 58% at a distance of 10 cm from the field edge. The newly designed fetal shield may effectively reduce fetal dose and is relatively easy to setup. Due to its design, we are able to use a broad range of treatment techniques and beam angles. We believe the acquired comprehensive PD dataset collected with and without the fetal shield will be useful for treatment teams to estimate fetal dose and help guide decisions on treat-ment techniques without the need to perform pretreatment phantom measurements.

Improving treatment plan evaluation with automation.

The goal of this work is to evaluate the effectiveness of Plan-Checker Tool (PCT) which was created to improve first-time plan quality, reduce patient delays, increase the efficiency of our electronic workflow, and standardize and automate the phys-ics plan review in the treatment planning system (TPS). PCT uses an application programming interface to check and compare data from the TPS and treatment management system (TMS). PCT includes a comprehensive checklist of automated and manual checks that are documented when performed by the user as part of a plan readiness check for treatment. Prior to and during PCT development, errors identified during the physics review and causes of patient treatment start delays were tracked to prioritize which checks should be automated. Nineteen of 33checklist items were automated, with data extracted with PCT. There was a 60% reduction in the number of patient delays in the six months after PCT release. PCT was suc-cessfully implemented for use on all external beam treatment plans in our clinic. While the number of errors found during the physics check did not decrease, automation of checks increased visibility of errors during the physics check, which led to decreased patient delays. The methods used here can be applied to any TMS and TPS that allows queries of the database.

SafetyNet: Streamlining and Automating QA in radiotherapy.

Proper quality assurance (QA) of the radiotherapy process can be time-consuming and expensive. Many QA efforts, such as data export and import, are inefficient when done by humans. Additionally, humans can be unreliable, lose attention, and fail to complete critical steps that are required for smooth operations. In our group we have sought to break down the QA tasks into separate steps and to automate those steps that are better done by software running autonomously or at the instigation of a human. A team of medical physicists and software engineers worked together to identify opportunities to streamline and automate QA. Development efforts follow a formal cycle of writing software requirements, developing software, testing and commissioning. The clinical release process is separated into clinical evaluation testing, training, and finally clinical release. We have improved six processes related to QA and safety. Steps that were previously performed by humans have been automated or streamlined to increase first-time quality, reduce time spent by humans doing low-level tasks, and expedite QA tests. Much of the gains were had by automating data transfer, implementing computer-based checking and automation of systems with an event-driven framework. These coordinated efforts by software engineers and clinical physicists have resulted in speed improvements in expediting patient-sensitive QA tests.

A deep learning segmentation method to assess dose to organs at risk during breast radiotherapy.

Background and purpose: Contouring of organs at risk is important for studying health effects following breast radiotherapy. However, manual contouring is time-consuming and subject to variability. The purpose of this study was to develop a deep learning-based method to automatically segment multiple structures on breast radiotherapy planning computed tomography (CT) images. Materials and methods: We used data from 118 patients, including 90 diagnostic CT scans with expert structure delineations for training and 28 breast radiotherapy planning CT images for testing. The radiotherapy CT images also had expert delineations for evaluating performance. We targeted a total of eleven organs at risk including five heart substructures. Segmentation performance was evaluated using the metrics of Dice similarity coefficient (DSC), overlap fraction, volume similarity, Hausdorff distance, mean surface distance, and dose. Results: The average DSC achieved on the radiotherapy planning images was 0.94 ± 0.02 for the whole heart, 0.96 ± 0.02 and 0.97 ± 0.01 for the left and right lung, 0.61 ± 0.10 for the esophagus, 0.81 ± 0.04 and 0.86 ± 0.04 for left and right atrium, 0.91 ± 0.02 and 0.84 ± 0.04 for left and right ventricle, and 0.21 ± 0.11 for the left anterior descending artery (LAD), respectively. Except for the LAD, the median difference in mean dose to these structures was small with absolute (relative) differences < 0.1 Gy (6 %). Conclusions: Except for the LAD, our method demonstrated excellent performance and can be generalized to segment additional structures of interest.

Dose prescription and reporting in stereotactic body radiotherapy: A multi-institutional study.

BACKGROUND AND PURPOSE: Radiation dose prescriptions are foundational for optimizing treatment efficacy and limiting treatment-related toxicity. We sought to assess the lack of standardization of SBRT dose prescriptions across institutions. MATERIALS & METHODS: Dosimetric data from 1298 patients from 9 academic institutions treated with IMRT and VMAT were collected. Dose parameters D100, D98, D95, D50, and D2 were used to assess dosimetric variability. RESULTS: Disease sites included lung (48.3 %) followed by liver (29.7 %), prostate (7.5 %), spine (6.8 %), brain (4.1 %), and pancreas (2.5 %). The PTV volume in lung varied widely with bimodality into two main groups (22.0-28.7 cm3) and (48.0-67.1 cm3). A hot spot ranging from 120-150 % was noted in nearly half of the patients, with significant variation across institutions. A D50 ≥ 110 % was found in nearly half of the institutions. There was significant dosimetric variation across institutions. CONCLUSIONS: The SBRT prescriptions in the literature or in treatment guidelines currently lack nuance and hence there is significant variation in dose prescriptions across academic institutions. These findings add greater importance to the identification of dose parameters associated with improved clinical outcome comparisons as we move towards more hypofractionated treatments. There is a need for standardized reporting to help institutions in adapting treatment protocols based on the outcome of clinical trials. Dosimetric parameters are subsequently needed for uniformity and thereby standardizing planning guidelines to maximize efficacy, mitigate toxicity, and reduce treatment disparities are urgently needed.

A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy.

PURPOSE: To provide the first clinical test case for commissioning of 192 Ir brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow. ACQUISITION AND VALIDATION METHODS: A computational patient phantom model was generated from a clinical multi-catheter 192 Ir HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic 192 Ir HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution. DATA FORMAT AND USAGE NOTES: The dataset is available online at http://irochouston.mdanderson.org/rpc/BrachySeeds/BrachySeeds/index.html,https://doi.org/10.52519/00005. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations. POTENTIAL APPLICATIONS: The dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.

Development and Clinical Implementation of an Automated Virtual Integrative Planner for Radiation Therapy of Head and Neck Cancer.

Purpose: Head and neck (HN) radiation (RT) treatment planning is complex and resource intensive. Deviations and inconsistent plan quality significantly affect clinical outcomes. We sought to develop a novel automated virtual integrative (AVI) knowledge-based planning application to reduce planning time, increase consistency, and improve baseline quality. Methods and Materials: An in-house write-enabled script was developed from a library of 668 previously treated HN RT plans. Prospective hazard analysis was performed, and mitigation strategies were implemented before clinical release. The AVI-planner software was retrospectively validated in a cohort of 52 recent HN cases. A physician panel evaluated planning limitations during initial deployment, and feedback was enacted via software refinements. A final second set of plans was generated and evaluated. Kolmogorov-Smirnov test in addition to generalized evaluation metric and weighted experience score were used to compare normal tissue sparing between final AVI planner versus respective clinically treated and historically accepted plans. A t test was used to compare the interactive time, complexity, and monitor units for AVI planner versus manual optimization. Results: Initially, 86% of plans were acceptable to treat, with 10% minor and 4% major revisions or rejection recommended. Variability was noted in plan quality among HN subsites, with high initial quality for oropharynx and oral cavity plans. Plans needing revisions were comprised of sinonasal, nasopharynx, P-16 negative squamous cell carcinoma unknown primary, or cutaneous primary sites. Normal tissue sparing varied within subsites, but AVI planner significantly lowered mean larynx dose (median, 18.5 vs 19.7 Gy; P < .01) compared with clinical plans. AVI planner significantly reduced interactive optimization time (mean, 2 vs 85 minutes; P < .01). Conclusions: AVI planner reliably generated clinically acceptable RT plans for oral cavity, salivary, oropharynx, larynx, and hypopharynx cancers. Physician-driven iterative learning processes resulted in favorable evolution in HN RT plan quality with significant time savings and improved consistency using AVI planner.

Organ dose conversion coefficients in CT scans for Korean adult males and females.

Dose monitoring in CT patients requires accurate dose estimation but most of the CT dose calculation tools are based on Caucasian computational phantoms. We established a library of organ dose conversion coefficients for Korean adults by using four Korean adult male and two female voxel phantoms combined with Monte Carlo simulation techniques. We calculated organ dose conversion coefficients for head, chest, abdomen and pelvis, and chest-abdomen-pelvis scans, and compared the results with the existing data calculated from Caucasian phantoms. We derived representative organ doses for Korean adults using Korean CT dose surveys combined with the dose conversion coefficients. The organ dose conversion coefficients from the Korean adult phantoms were slightly greater than those of the ICRP reference phantoms: up to 13% for the brain doses in head scans and up to 10% for the dose to the small intestine wall in abdominal scans. We derived Korean representative doses to major organs in head, chest, and AP scans using mean CTDIvol values extracted from the Korean nationwide surveys conducted in 2008 and 2017. The Korean-specific organ dose conversion coefficients should be useful to readily estimate organ absorbed doses for Korean adult male and female patients undergoing CT scans.

Fetal dose from proton pencil beam scanning craniospinal irradiation during pregnancy: a Monte Carlo study.

Objective. We conducted a Monte Carlo study to comprehensively investigate the fetal dose resulting from proton pencil beam scanning (PBS) craniospinal irradiation (CSI) during pregnancy.Approach. The gestational-age dependent pregnant phantom series developed at the University of Florida (UF) were converted into DICOM-RT format (CT images and structures) and imported into a treatment planning system (TPS) (Eclipse v15.6) commissioned to a IBA PBS nozzle. A proton PBS CSI plan (prescribed dose: 36 Gy) was created on the phantoms. The TOPAS MC code was used to simulate the proton PBS CSI on the phantoms, for which MC beam properties at the nozzle exit (spot size, spot divergence, mean energy, and energy spread) were matched to IBA PBS nozzle beam measurement data. We calculated mean absorbed doses for 28 organs and tissues and whole body of the fetus at eight gestational ages (8, 10, 15, 20, 25, 30, 35, and 38 weeks). For contextual purposes, the fetal organ/tissue doses from the treatment planning CT scan of the mother's head and torso were estimated using the National Cancer Institute dosimetry system for CT (NCICT, Version 3) considering a low-dose CT protocol (CTDIvol: 8.97 mGy).Main results. The majority of the fetal organ/tissue doses from the proton PBS CSI treatment fell within a range of 3-6 mGy. The fetal organ/tissue doses for the 38 week phantom showed the largest variation with the doses ranging from 2.9 mGy (adrenals) to 8.2 mGy (eye lenses) while the smallest variation ranging from 3.2 mGy (oesophagus) to 4.4 mGy (brain) was observed for the doses for the 20 week phantom. The fetal whole-body dose ranged from 3.7 mGy (25 weeks) to 5.8 mGy (8 weeks). Most of the fetal doses from the planning CT scan fell within a range of 7-13 mGy, approximately 2-to-9 times lower than the fetal dose equivalents of the proton PBS CSI treatment (assuming a quality factor of 7).Significance. The fetal organ/tissue doses observed in the present work will be useful for one of the first clinically informative predictions on the magnitude of fetal dose during proton PBS CSI during pregnancy.

Comparison of out-of-field normal tissue dose estimates for pencil beam scanning proton therapy: MCNP6, PHITS, and TOPAS.

Monte Carlo (MC) methods are considered the gold-standard approach to dose estimation for normal tissues outside the treatment field (out-of-field) in proton therapy. However, the physics of secondary particle production from high-energy protons are uncertain, particularly for secondary neutrons, due to challenges in performing accurate measurements. Instead, various physics models have been developed over the years to reenact these high-energy interactions based on theory. It should thus be acknowledged that MC users must currently accept some unknown uncertainties in out-of-field dose estimates. In the present study, we compared three MC codes (MCNP6, PHITS, and TOPAS) and their available physics models to investigate the variation in out-of-field normal tissue dosimetry for pencil beam scanning proton therapy patients. Total yield and double-differential (energy and angle) production of two major secondary particles, neutrons and gammas, were determined through irradiation of a water phantom at six proton energies (80, 90, 100, 110, 150, and 200 MeV). Out-of-field normal tissue doses were estimated for intracranial irradiations of 1-, 5-, and 15-year-old patients using whole-body computational phantoms. Notably, the total dose estimates for each out-of-field organ varied by approximately 25% across the three codes, independent of its distance from the treatment volume. Dose discrepancies amongst the codes were linked to the utilized physics model, which impacts the characteristics of the secondary radiation field. Using developer-recommended physics, TOPAS produced both the highest neutron and gamma doses to all out-of-field organs from all examined conditions; this was linked to its highest yields of secondary particles and second hardest energy spectra. Subsequent results when using other physics models found reduced yields and energies, resulting in lower dose estimates. Neutron dose estimates were the most impacted by physics model choice, and thus the variation in out-of-field dose estimates may be even larger than 25% when considering biological effectiveness.

Validation and Comparison of Radiograph-Based Organ Dose Reconstruction Approaches for Wilms Tumor Radiation Treatment Plans.

Purpose: Our purpose was to validate and compare the performance of 4 organ dose reconstruction approaches for historical radiation treatment planning based on 2-dimensional radiographs. Methods and Materials: We considered 10 patients with Wilms tumor with planning computed tomography images for whom we developed typical historic Wilms tumor radiation treatment plans, using anteroposterior and posteroanterior parallel-opposed 6 MV flank fields, normalized to 14.4 Gy. Two plans were created for each patient, with and without corner blocking. Regions of interest (lungs, heart, nipples, liver, spleen, contralateral kidney, and spinal cord) were delineated, and dose-volume metrics including organ mean and minimum dose (Dmean and Dmin) were computed as the reference baseline for comparison. Dosimetry for the 20 plans was then independently reconstructed using 4 different approaches. Three approaches involved surrogate anatomy, among which 2 used demographic-matching criteria for phantom selection/building, and 1 used machine learning. The fourth approach was also machine learning-based, but used no surrogate anatomies. Absolute differences in organ dose-volume metrics between the reconstructed and the reference values were calculated. Results: For Dmean and Dmin (average and minimum point dose) all 4 dose reconstruction approaches performed within 10% of the prescribed dose (≤1.4 Gy). The machine learning-based approaches showed a slight advantage for several of the considered regions of interest. For Dmax (maximum point dose), the absolute differences were much higher, that is, exceeding 14% (2 Gy), with the poorest agreement observed for near-beam and out-of-beam organs for all approaches. Conclusions: The studied approaches give comparable dose reconstruction results, and the choice of approach for cohort dosimetry for late effects studies should still be largely driven by the available resources (data, time, expertise, and funding).

Early HPV ctDNA Kinetics and Imaging Biomarkers Predict Therapeutic Response in p16+ Oropharyngeal Squamous Cell Carcinoma.

PURPOSE: In locally advanced p16+ oropharyngeal squamous cell carcinoma (OPSCC), (i) to investigate kinetics of human papillomavirus (HPV) circulating tumor DNA (ctDNA) and association with tumor progression after chemoradiation, and (ii) to compare the predictive value of ctDNA to imaging biomarkers of MRI and FDG-PET. EXPERIMENTAL DESIGN: Serial blood samples were collected from patients with AJCC8 stage III OPSCC (n = 34) enrolled on a randomized trial: pretreatment; during chemoradiation at weeks 2, 4, and 7; and posttreatment. All patients also had dynamic-contrast-enhanced and diffusion-weighted MRI, as well as FDG-PET scans pre-chemoradiation and week 2 during chemoradiation. ctDNA values were analyzed for prediction of freedom from progression (FFP), and correlations with aggressive tumor subvolumes with low blood volume (TVLBV) and low apparent diffusion coefficient (TVLADC), and metabolic tumor volume (MTV) using Cox proportional hazards model and Spearman rank correlation. RESULTS: Low pretreatment ctDNA and an early increase in ctDNA at week 2 compared with baseline were significantly associated with superior FFP (P < 0.02 and P < 0.05, respectively). At week 4 or 7, neither ctDNA counts nor clearance were significantly predictive of progression (P = 0.8). Pretreatment ctDNA values were significantly correlated with nodal TVLBV, TVLADC, and MTV pre-chemoradiation (P < 0.03), while the ctDNA values at week 2 were correlated with these imaging metrics in primary tumor. Multivariate analysis showed that ctDNA and the imaging metrics performed comparably to predict FFP. CONCLUSIONS: Early ctDNA kinetics during definitive chemoradiation may predict therapy response in stage III OPSCC.

Randomized Phase II Study of Physiologic MRI-Directed Adaptive Radiation Boost in Poor Prognosis Head and Neck Cancer.

PURPOSE: We conducted a randomized phase II multicenter clinical trial to test the hypothesis that physiologic MRI-based radiotherapy (RT) dose escalation would improve the outcome of patients with poor prognosis head and neck cancer. PATIENTS AND METHODS: MRI was acquired at baseline and at RT fraction 10 to create low blood volume/apparent diffusion coefficient maps for RT boost subvolume definition in gross tumor volume. Patients were randomized to receive 70 Gy (standard RT) or 80 Gy to the boost subvolume (RT boost) with concurrent weekly platinum. The primary endpoint was disease-free survival (DFS) with significance defined at a one-sided 0.1 level, and secondary endpoints included locoregional failure (LRF), overall survival (OS), comparison of adverse events and patient reported outcomes (PRO). RESULTS: Among 81 randomized patients, neither the primary endpoint of DFS (HR = 0.849, P = 0.31) nor OS (HR = 1.19, P = 0.66) was significantly improved in the RT boost arm. However, the incidence of LRF was significantly improved with the addition of the RT boost (HR = 0.43, P = 0.047). Two-year estimates [90% confidence interval (CI)] of the cumulative incidence of LRF were 40% (27%-53%) in the standard RT arm and 18% (10%-31%) in the RT boost arm. Two-year estimates (90% CI) for DFS were 48% (34%-60%) in the standard RT arm and 57% (43%-69%) in the RT boost arm. There were no significant differences in toxicity or longitudinal differences seen in EORTC QLQ30/HN35 subscales between treatment arms in linear mixed-effects models. CONCLUSIONS: Physiologic MRI-based RT boost decreased LRF without a significant increase in grade 3+ toxicity or longitudinal PRO differences, but did not significantly improve DFS or OS. Additional improvements in systemic therapy are likely necessary to realize improvements in DFS and OS.

Application of an automatic segmentation method for evaluating cardiac structure doses received by breast radiotherapy patients.

BACKGROUND AND PURPOSE: Quantifying radiation dose to cardiac substructures is important for research on the etiology and prevention of complications following radiotherapy; however, segmentation of substructures is challenging. In this study we demonstrate the application of our atlas-based automatic segmentation method to breast cancer radiotherapy plans for generating radiation doses in support of late effects research. MATERIAL AND METHODS: We applied our segmentation method to contour heart substructures on the computed tomography (CT) images of 70 breast cancer patients who received external photon radiotherapy. Two cardiologists provided manual segmentation of the whole heart (WH), left/right atria, left/right ventricles, and left anterior descending artery (LAD). The automatically contours were compared with manual delineations to evaluate similarity in terms of geometry and dose. RESULTS: The mean Dice similarity coefficient between manual and automatic segmentations was 0.96 for the WH, 0.65 to 0.82 for the atria and ventricles, and 0.06 for the LAD. The mean average surface distance was 1.2 mm for the WH, 3.4 to 4.1 mm for the atria and ventricles, and 6.4 mm for the LAD. We found the dose to the cardiac substructures based on our automatic segmentation agrees with manual segmentation within expected observer variability. For left breast patients, the mean absolute difference in mean dose was 0.1 Gy for the WH, 0.2 to 0.7 Gy for the atria and ventricles, and 1.8 Gy for the LAD. For right breast patients, these values were 0.0 Gy, 0.1 to 0.4 Gy, and 0.4 Gy, respectively. CONCLUSION: Our automatic segmentation method will facilitate the development of radiotherapy prescriptive criteria for mitigating cardiovascular complications.

Early MRI Blood Volume Changes in Constrictor Muscles Correlate With Postradiation Dysphagia.

PURPOSE: Predicting individual patient sensitivity to radiation therapy (RT) for tumor control or normal tissue toxicity is necessary to individualize treatment planning. In head and neck cancer, radiation doses are limited by many nearby critical structures, including structures involved in swallowing. Previous efforts showed that imaging parameters correlate with RT dose; here, we investigate the role of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) blood volume (BV) changes in predicting dysphagia. METHODS AND MATERIALS: This study included 32 patients with locally advanced oropharyngeal squamous cell carcinoma treated with definitive chemoradiation on an institutional protocol incorporating baseline and early midtreatment DCE-MRI. BV maps of the pharyngeal constrictor muscles (PCM) were created, and BV increases midtreatment were correlated with the following parameters at 3 and 12 months post-RT: RT dose, Dynamic Imaging Grade of Swallowing Toxicity swallow score, aspiration frequency, European Organisation for Research and Treatment of Cancer HN35 patient-reported outcomes, physician-reported dysphagia, and feeding tube (FT) dependence. RESULTS: The mean BV to the PCMs increased from baseline to fraction 10, which was significant for the superior PCM (P = .006) and middle PCM (P < .001), with a trend in the inferior PCM where lower mean doses were seen (P = .077). The factors associated with FT dependence at 3 months included BV increases in the total PCM (correlation, 0.48; P = .006) and middle PCM (correlation, 0.50; P = .004). A post-RT increase in aspiration was associated with a BV increase in the superior PCM (correlation, 0.44; P = .013),and the increase in the total PCMs was marginally significant (correlation, 0.34; P = .06). The best-performing models of FT dependence (area under the receiver operating curve [AUC] = 0.84) and aspiration increases (AUC = 0.78) included BV increases as well as a mean RT dose to middle PCM. CONCLUSIONS: Our results suggest that midtreatment BV increases derived from DCE-MRI are an early predictor of dysphagia. Further investigation of these promising imaging markers to assess individual patient sensitivity to treatment and the patient's subsequent risk of toxicities is warranted to improve personalization of RT planning.

A dose voxel kernel method for rapid reconstruction of out-of-field neutron dose of patients in pencil beam scanning (PBS) proton therapy.

Monte Carlo (MC) radiation transport methods are used for dose calculation as 'gold standard.' However, the method is computationally time-consuming and thus impractical for normal tissue dose reconstructions for the large number of proton therapy patients required for epidemiologic investigations of late health effects. In the present study, we developed a new dose calculation method for the rapid reconstruction of out-of-field neutron dose to patients undergoing pencil beam scanning (PBS) proton therapy. The new dose calculation method is based on neutron dose voxel kernels (DVKs) generated by MC simulations of a proton pencil beam irradiating a water phantom (60 × 60 × 300 cm3), which was conducted using a MC proton therapy simulation code, TOPAS. The DVKs were generated for 19 beam energies (from 70 to 250 MeV with the 10 MeV interval) and three range shifter thicknesses (1, 3, and 5 cm). An in-house program was written in C++ to superimpose the DVKs onto a patient CT images according to proton beam characteristics (energy, position, and direction) available in treatment plans. The DVK dose calculation method was tested by calculating organ/tissue-specific neutron doses of 1- and 5-year-old whole-body computational phantoms where intracranial and craniospinal irradiations were simulated. The DVK-based doses generally showed reasonable agreement with those calculated by direct MC simulations with a detailed PBS model that were previously published, with differences mostly less than 30% and 10% for the intracranial and craniospinal irradiations, respectively. The computation time of the DVK method for one patient ranged from 1 to 30 min on a single CPU core of a personal computer, demonstrating significant improvement over the direct MC dose calculation requiring several days on high-performance computing servers. Our DVK-based dose calculation method will be useful when dosimetry is needed for the large number of patients such as for epidemiologic or clinical research.

Development and comprehensive commissioning of an automated brachytherapy plan checker.

PURPOSE: To present on the commissioning of an automated brachytherapy plan checker (BPC) for the evaluation of high-dose-rate brachytherapy treatment plans in support of standardized workflows and patient safety. METHODS AND MATERIALS: A BPC was developed using an applications programming interface in a commercial treatment planning system based on different inputs (e.g., regulations, professional society recommendations, and user feedback) and leveraged our experience with an in-house developed external beam plan checker. The BPC was commissioned using a comprehensive suite of test plans with known errors and anonymized clinical plans. RESULTS: During commissioning, the BPC was successfully executed on a total of 87 test plans. Commissioning tests spanned a range of treatment sites and evaluated that pass and fail states were correct. Administration settings were changed in a nonclinical database to evaluate tests involving the source and afterloader. Clinical testing of the BPC was then performed in parallel with a manual review process before clinical implementation. CONCLUSIONS: To commission the BPC for clinical use, a comprehensive suite of test plans was developed and used to ensure the BPC correctly detected and reported errors. A summary of the test plans is presented to help guide users developing similar automated tools. The BPC represents a process-improvement initiative designed to reduce errors and improve safety for brachytherapy patients. By using a comprehensive test suite for commissioning, tests are available for periodic quality assurance and after software upgrades.