Dynamic stochastic deep learning approaches for predicting geometric changes in head and neck cancer.

Objective.Modern radiotherapy stands to benefit from the ability to efficiently adapt plans during treatment in response to setup and geometric variations such as those caused by internal organ deformation or tumor shrinkage. A promising strategy is to develop a framework, which given an initial state defined by patient-attributes, can predict future states based on pre-learned patterns from a well-defined patient population.Approach.Here, we investigate the feasibility of predicting patient anatomical changes, defined as a joint state of volume and daily setup changes, across a fractionated treatment schedule using two approaches. The first is based on a new joint framework employing quantum mechanics in combination with deep recurrent neural networks, denoted QRNN. The second approach is developed based on a classical framework, which models patient changes as a Markov process, denoted MRNN. We evaluated the performance characteristics of these two approaches on a dataset of 125 head and neck cancer patients, which was supplemented by synthetic data generated using a generative adversarial network. Model performance was evaluated using area under the receiver operating characteristic curve (AUC) scores.Main results.The MRNN framework had slightly better performance than the QRNN framework, with MRNN (QRNN) validation AUC scores of 0.742±0.021 (0.675±0.036), 0.709±0.026 (0.656±0.021), 0.724±0.036 (0.652±0.044), and 0.698±0.016 (0.605±0.035) for system state vector sizes of 4, 6, 8, and 10, respectively. Of these, only the results from the two higher order states had statistically significant differences(p<0.05).A similar trend was also observed when the models were applied to an external testing dataset of 20 patients, yielding MRNN (QRNN) AUC scores of 0.707 (0.623), 0.687 (0.608), 0.723 (0.669), and 0.697 (0.609) for states vectors sizes of 4, 6, 8, and 10, respectively.Significance.These results suggest that both stochastic models have potential value in predicting patient changes during the course of adaptive radiotherapy.

A situational awareness Bayesian network approach for accurate and credible personalized adaptive radiotherapy outcomes prediction in lung cancer patients.

PURPOSE: A situational awareness Bayesian network (SA-BN) approach is developed to improve physicians' trust in the prediction of radiation outcomes and evaluate its performance for personalized adaptive radiotherapy (pART). METHODS: 118 non-small-cell lung cancer patients with their biophysical features were employed for discovery (n = 68) and validation (n = 50) of radiation outcomes prediction modeling. Patients' important characteristics identified by radiation experts to predict individual's tumor local control (LC) or radiation pneumonitis with grade ≥ 2 (RP2) were incorporated as expert knowledge (EK). Besides generating an EK-based naïve BN (EK-NBN), an SA-BN was developed by incorporating the EK features into pure data-driven BN (PD-BN) methods to improve the credibility of LC or / and RP2 prediction. After using area under the free-response receiver operating characteristics curve (AU-FROC) to assess the joint prediction of these outcomes, their prediction performances were compared with a regression approach based on the expert yielded estimates (EYE) penalty and its variants. RESULTS: In addition to improving the credibility of radiation outcomes prediction, the SA-BN approach outperformed the EYE penalty and its variants in terms of the joint prediction of LC and RP2. The value of AU-FROC improves from 0.70 (95% CI: 0.54-0.76) using EK-NBN, to 0.75 (0.65-0.82) using a variant of EYE penalty, to 0.83 (0.75-0.93) using PD-BN and 0.83 (0.77-0.90) using SA-BN; with similar trends in the validation cohort. CONCLUSIONS: The SA-BN approach can provide an accurate and credible human-machine interface to gain physicians' trust in clinical decision-making, which has the potential to be an important component of pART.

Quantum-inspired algorithm for radiotherapy planning optimization.

PURPOSE: Modern inverse radiotherapy treatment planning requires nonconvex, large-scale optimizations that must be solved within a clinically feasible timeframe. We have developed and tested a quantum-inspired, stochastic algorithm for intensity-modulated radiotherapy (IMRT): quantum tunnel annealing (QTA). By modeling the likelihood probability of accepting a higher energy solution after a particle tunneling through a potential energy barrier, QTA features an additional degree of freedom (the barrier width, w) not shared by traditional stochastic optimization methods such as Simulated Annealing (SA). This additional degree of freedom can improve convergence rates and achieve a more efficient and, potentially, effective treatment planning process. METHODS: To analyze the character of the proposed QTA algorithm, we chose two stereotactic body radiation therapy (SBRT) liver cases of variable complexity. The "easy" first case was used to confirm functionality, while the second case, with a more challenging geometry, was used to characterize and evaluate the QTA algorithm performance. Plan quality was assessed using dose-volume histogram-based objectives and dose distributions. Due to the stochastic nature of the solution search space, extensive tests were also conducted to determine the optimal smoothing technique, ensuring balance between plan deliverability and the resulting plan quality. QTA convergence rates were investigated in relation to the chosen barrier width function, and QTA and SA performances were compared regarding sensitivity to the choice of solution initializations, annealing schedules, and complexity of the dose-volume constraints. Finally, we investigated the extension from beamlet intensity optimization to direct aperture optimization (DAO). Influence matrices were calculated using the Eclipse scripting application program interface (API), and the optimizations were run on the University of Michigan's high-performance computing cluster, Flux. RESULTS: Our results indicate that QTA's barrier-width function can be tuned to achieve faster convergence rates. The QTA algorithm reached convergence up to 46.6% faster than SA for beamlet intensity optimization and up to 26.8% faster for DAO. QTA and SA were ultimately found to be equally insensitive to the initialization process, but the convergence rate of QTA was found to be more sensitive to the complexity of the dose-volume constraints. The optimal smoothing technique was found to be a combination of a Laplace-of-Gaussian (LOG) edge-finding filter implemented as a penalty within the objective function and a two-dimensional Savitzky-Golay filter applied to the final iteration; this achieved total monitor units more than 20% smaller than plans optimized by commercial treatment planning software. CONCLUSIONS: We have characterized the performance of a stochastic, quantum-inspired optimization algorithm, QTA, for radiotherapy treatment planning. This proof of concept study suggests that QTA can be tuned to achieve faster convergence than SA; therefore, QTA may be a good candidate for future knowledge-based or adaptive radiation therapy applications.

Development of a Fully Cross-Validated Bayesian Network Approach for Local Control Prediction in Lung Cancer.

The purpose of this study is to demonstrate that a Bayesian network (BN) approach can explore hierarchical biophysical relationships that influence tumor response and predict tumor local control (LC) in non-small-cell lung cancer (NSCLC) patients before and during radiotherapy from a large-scale dataset. Our BN building approach has two steps. First, relevant biophysical predictors influencing LC before and during the treatment are selected through an extended Markov blanket (eMB) method. From this eMB process, the most robust BN structure for LC prediction was found via a wrapper-based approach. Sixty-eight patients with complete feature information were used to identify a full BN model for LC prediction before and during the treatment. Fifty more recent patients with some missing information were reserved for independent testing of the developed pre- and during-therapy BNs. A nested cross-validation (N-CV) was developed to evaluate the performance of the two-step BN approach. An ensemble BN model is generated from the N-CV sampling process to assess its similarity with the corresponding full BN model, and thus evaluate the sensitivity of our BN approach. Our results show that the proposed BN development approach is a stable and robust approach to identify hierarchical relationships among biophysical features for LC prediction. Furthermore, BN predictions can be improved by incorporating during treatment information.

A multiobjective Bayesian networks approach for joint prediction of tumor local control and radiation pneumonitis in nonsmall-cell lung cancer (NSCLC) for response-adapted radiotherapy.

PURPOSE: Individualization of therapeutic outcomes in NSCLC radiotherapy is likely to be compromised by the lack of proper balance of biophysical factors affecting both tumor local control (LC) and side effects such as radiation pneumonitis (RP), which are likely to be intertwined. Here, we compare the performance of separate and joint outcomes predictions for response-adapted personalized treatment planning. METHODS: A total of 118 NSCLC patients treated on prospective protocols with 32 cases of local progression and 20 cases of RP grade 2 or higher (RP2) were studied. Sixty-eight patients with 297 features before and during radiotherapy were used for discovery and 50 patients were reserved for independent testing. A multiobjective Bayesian network (MO-BN) approach was developed to identify important features for joint LC/RP2 prediction using extended Markov blankets as inputs to develop a BN predictive structure. Cross-validation (CV) was used to guide the MO-BN structure learning. Area under the free-response receiver operating characteristic (AU-FROC) curve was used to evaluate joint prediction performance. RESULTS: Important features including single nucleotide polymorphisms (SNPs), micro RNAs, pretreatment cytokines, pretreatment PET radiomics together with lung and tumor gEUDs were selected and their biophysical inter-relationships with radiation outcomes (LC and RP2) were identified in a pretreatment MO-BN. The joint LC/RP2 prediction yielded an AU-FROC of 0.80 (95% CI: 0.70-0.86) upon internal CV. This improved to 0.85 (0.75-0.91) with additional two SNPs, changes in one cytokine and two radiomics PET image features through the course of radiotherapy in a during-treatment MO-BN. This MO-BN model outperformed combined single-objective Bayesian networks (SO-BNs) during-treatment [0.78 (0.67-0.84)]. AU-FROC values in the evaluation of the MO-BN and individual SO-BNs on the testing dataset were 0.77 and 0.68 for pretreatment, and 0.79 and 0.71 for during-treatment, respectively. CONCLUSIONS: MO-BNs can reveal possible biophysical cross-talks between competing radiotherapy clinical endpoints. The prediction is improved by providing additional during-treatment information. The developed MO-BNs can be an important component of decision support systems for personalized response-adapted radiotherapy.

Unraveling biophysical interactions of radiation pneumonitis in non-small-cell lung cancer via Bayesian network analysis.

BACKGROUND: In non-small-cell lung cancer radiotherapy, radiation pneumonitis≥grade 2 (RP2) depends on patients' dosimetric, clinical, biological and genomic characteristics. METHODS: We developed a Bayesian network (BN) approach to explore its potential for interpreting biophysical signaling pathways influencing RP2 from a heterogeneous dataset including single nucleotide polymorphisms, micro RNAs, cytokines, clinical data, and radiation treatment plans before and during the course of radiotherapy. Model building utilized 79 patients (21 with RP2) with complete data, and model testing used 50 additional patients with incomplete data. A developed large-scale Markov blanket approach selected relevant predictors. Resampling by k-fold cross-validation determined the optimal BN structure. Area under the receiver-operating characteristics curve (AUC) measured performance. RESULTS: Pre- and during-treatment BNs identified biophysical signaling pathways from the patients' relevant variables to RP2 risk. Internal cross-validation for the pre-BN yielded an AUC=0.82 which improved to 0.87 by incorporating during treatment changes. In the testing dataset, the pre- and during AUCs were 0.78 and 0.82, respectively. CONCLUSIONS: Our developed BN approach successfully handled a high number of heterogeneous variables in a small dataset, demonstrating potential for unraveling relevant biophysical features that could enhance prediction of RP2, although the current observations would require further independent validation.

Decreased Lung Perfusion After Breast/Chest Wall Irradiation: Quantitative Results From a Prospective Clinical Trial.

PURPOSE: To quantify lung perfusion changes after breast/chest wall radiation therapy (RT) using pre- and post-RT single photon emission computed tomography/computed tomography (SPECT/CT) attenuation-corrected perfusion scans; and correlate decreased perfusion with adjuvant RT dose for breast cancer in a prospective clinical trial. METHODS AND MATERIALS: As part of an institutional review board-approved trial studying the impact of RT technique on lung function in node-positive breast cancer, patients received breast/chest wall and regional nodal irradiation including superior internal mammary node RT to 50 to 52.2 Gy with a boost to the tumor bed/mastectomy scar. All patients underwent quantitative SPECT/CT lung perfusion scanning before RT and 1 year after RT. The SPECT/CT scans were co-registered, and the ratio of decreased perfusion after RT relative to the pre-RT perfusion scan was calculated to allow for direct comparison of SPECT/CT perfusion changes with delivered RT dose. The average ratio of decreased perfusion was calculated in 10-Gy dose increments from 0 to 60 Gy. RESULTS: Fifty patients had complete lung SPECT/CT perfusion data available. No patient developed symptoms consistent with pulmonary toxicity. Nearly all patients demonstrated decreased perfusion in the left lung according to voxel-based analyses. The average ratio of lung perfusion deficits increased for each 10-Gy increment in radiation dose to the lung, with the largest changes in regions of lung that received 50 to 60 Gy (ratio 0.72 [95% confidence interval 0.64-0.79], P<.001) compared with the 0- to 10-Gy region. For each increase in 10 Gy to the left lung, the lung perfusion ratio decreased by 0.06 (P<.001). CONCLUSIONS: In the assessment of 50 patients with node-positive breast cancer treated with RT in a prospective clinical trial, decreased lung perfusion by SPECT/CT was demonstrated. Our study allowed for quantification of lung perfusion defects in a prospective cohort of breast cancer patients for whom attenuation-corrected SPECT/CT scans could be registered directly to RT treatment fields for precise dose estimates.

The big data effort in radiation oncology: Data mining or data farming?

Although large volumes of information are entered into our electronic health care records, radiation oncology information systems and treatment planning systems on a daily basis, the goal of extracting and using this big data has been slow to emerge. Development of strategies to meet this goal is aided by examining issues with a data farming instead of a data mining conceptualization. Using this model, a vision of key data elements, clinical process changes, technology issues and solutions, and role for professional societies is presented. With a better view of technology, process and standardization factors, definition and prioritization of efforts can be more effectively directed.

Priority-driven plan optimization in locally advanced lung patients based on perfusion SPECT imaging.

PURPOSE: Limits on mean lung dose (MLD) allow for individualization of radiation doses at safe levels for patients with lung tumors. However, MLD does not account for individual differences in the extent or spatial distribution of pulmonary dysfunction among patients, which leads to toxicity variability at the same MLD. We investigated dose rearrangement to minimize the radiation dose to the functional lung as assessed by perfusion single photon emission computed tomography (SPECT) and maximize the target coverage to maintain conventional normal tissue limits. METHODS AND MATERIALS: Retrospective plans were optimized for 15 patients with locally advanced non-small cell lung cancer who were enrolled in a prospective imaging trial. A staged, priority-based optimization system was used. The baseline priorities were to meet physical MLD and other dose constraints for organs at risk, and to maximize the target generalized equivalent uniform dose (gEUD). To determine the benefit of dose rearrangement with perfusion SPECT, plans were reoptimized to minimize the generalized equivalent uniform functional dose (gEUfD) to the lung as the subsequent priority. RESULTS: When only physical MLD is minimized, lung gEUfD was 12.6 ± 4.9 Gy (6.3-21.7 Gy). When the dose is rearranged to minimize gEUfD directly in the optimization objective function, 10 of 15 cases showed a decrease in lung gEUfD of >20% (lung gEUfD mean 9.9 ± 4.3 Gy, range 2.1-16.2 Gy) while maintaining equivalent planning target volume coverage. Although all dose-limiting constraints remained unviolated, the dose rearrangement resulted in slight gEUD increases to the cord (5.4 ± 3.9 Gy), esophagus (3.0 ± 3.7 Gy), and heart (2.3 ± 2.6 Gy). CONCLUSIONS: Priority-driven optimization in conjunction with perfusion SPECT permits image guided spatial dose redistribution within the lung and allows for a reduced dose to the functional lung without compromising target coverage or exceeding conventional limits for organs at risk.

Parotid glands dose-effect relationships based on their actually delivered doses: implications for adaptive replanning in radiation therapy of head-and-neck cancer.

PURPOSE: Doses actually delivered to the parotid glands during radiation therapy often exceed planned doses. We hypothesized that the delivered doses correlate better with parotid salivary output than the planned doses, used in all previous studies, and that determining these correlations will help make decisions regarding adaptive radiation therapy (ART) aimed at reducing the delivered doses. METHODS AND MATERIALS: In this prospective study, oropharyngeal cancer patients treated definitively with chemoirradiation underwent daily cone-beam computed tomography (CBCT) with clinical setup alignment based on the C2 posterior edge. Parotid glands in the CBCTs were aligned by deformable registration to calculate cumulative delivered doses. Stimulated salivary flow rates were measured separately from each parotid gland pretherapy and periodically posttherapy. RESULTS: Thirty-six parotid glands of 18 patients were analyzed. Average mean planned doses was 32 Gy, and differences from planned to delivered mean gland doses were -4.9 to +8.4 Gy, median difference +2.2 Gy in glands in which delivered doses increased relative to planned. Both planned and delivered mean doses were significantly correlated with posttreatment salivary outputs at almost all posttherapy time points, without statistically significant differences in the correlations. Large dispersions (on average, SD 3.6 Gy) characterized the dose-effect relationships for both. The differences between the cumulative delivered doses and planned doses were evident at first fraction (r=.92, P<.0001) because of complex setup deviations (eg, rotations and neck articulations), uncorrected by the translational clinical alignments. CONCLUSIONS: After daily translational setup corrections, differences between planned and delivered doses in most glands were small relative to the SDs of the dose-saliva data, suggesting that ART is not likely to gain measurable salivary output improvement in most cases. These differences were observed at first treatment, indicating potential benefit for more complex setup corrections or adaptive interventions in the minority of patients with large deviations detected early by CBCT.

FusionArc optimization: a hybrid volumetric modulated arc therapy (VMAT) and intensity modulated radiation therapy (IMRT) planning strategy.

PURPOSE: To introduce a hybrid volumetric modulated arc therapy/intensity modulated radiation therapy (VMAT/IMRT) optimization strategy called FusionArc that combines the delivery efficiency of single-arc VMAT with the potentially desirable intensity modulation possible with IMRT. METHODS: A beamlet-based inverse planning system was enhanced to combine the advantages of VMAT and IMRT into one comprehensive technique. In the hybrid strategy, baseline single-arc VMAT plans are optimized and then the current cost function gradients with respect to the beamlets are used to define a metric for predicting which beam angles would benefit from further intensity modulation. Beams with the highest metric values (called the gradient factor) are converted from VMAT apertures to IMRT fluence, and the optimization proceeds with the mixed variable set until convergence or until additional beams are selected for conversion. One phantom and two clinical cases were used to validate the gradient factor and characterize the FusionArc strategy. Comparisons were made between standard IMRT, single-arc VMAT, and FusionArc plans with one to five IMRT∕hybrid beams. RESULTS: The gradient factor was found to be highly predictive of the VMAT angles that would benefit plan quality the most from beam modulation. Over the three cases studied, a FusionArc plan with three converted beams achieved superior dosimetric quality with reductions in final cost ranging from 26.4% to 48.1% compared to single-arc VMAT. Additionally, the three beam FusionArc plans required 22.4%-43.7% fewer MU∕Gy than a seven beam IMRT plan. While the FusionArc plans with five converted beams offer larger reductions in final cost--32.9%-55.2% compared to single-arc VMAT--the decrease in MU∕Gy compared to IMRT was noticeably smaller at 12.2%-18.5%, when compared to IMRT. CONCLUSIONS: A hybrid VMAT∕IMRT strategy was implemented to find a high quality compromise between gantry-angle and intensity-based degrees of freedom. This optimization method will allow patients to be simultaneously planned for dosimetric quality and delivery efficiency without switching between delivery techniques. Example phantom and clinical cases suggest that the conversion of only three VMAT segments to modulated beams may result in a good combination of quality and efficiency.

The dosimetric impact of prostate rotations during electromagnetically guided external-beam radiation therapy.

PURPOSE: To study the impact of daily rotations and translations of the prostate on dosimetric coverage during radiation therapy (RT). METHODS AND MATERIALS: Real-time tracking data for 26 patients were obtained during RT. Intensity modulated radiation therapy plans meeting RTOG 0126 dosimetric criteria were created with 0-, 2-, 3-, and 5-mm planning target volume (PTV) margins. Daily translations and rotations were used to reconstruct prostate delivered dose from the planned dose. D95 and V79 were computed from the delivered dose to evaluate target coverage and the adequacy of PTV margins. Prostate equivalent rotation is a new metric introduced in this study to quantify prostate rotations by accounting for prostate shape and length of rotational lever arm. RESULTS: Large variations in prostate delivered dose were seen among patients. Adequate target coverage was met in 39%, 65%, and 84% of the patients for plans with 2-, 3-, and 5-mm PTV margins, respectively. Although no correlations between prostate delivered dose and daily rotations were seen, the data showed a clear correlation with prostate equivalent rotation. CONCLUSIONS: Prostate rotations during RT could cause significant underdosing even if daily translations were managed. These rotations should be managed with rotational tolerances based on prostate equivalent rotations.

Penalization of aperture complexity in inversely planned volumetric modulated arc therapy.

PURPOSE: Apertures obtained during volumetric modulated arc therapy (VMAT) planning can be small and irregular, resulting in dosimetric inaccuracies during delivery. Our purpose is to develop and integrate an aperture-regularization objective function into the optimization process for VMAT, and to quantify the impact of using this objective function on dose delivery accuracy and optimized dose distributions. METHODS: An aperture-based metric ("edge penalty") was developed that penalizes complex aperture shapes based on the ratio of MLC side edge length and aperture area. To assess the utility of the metric, VMAT plans were created for example paraspinal, brain, and liver SBRT cases with and without incorporating the edge penalty in the cost function. To investigate the dose calculation accuracy, Gafchromic EBT2 film was used to measure the 15 highest weighted apertures individually and as a composite from each of two paraspinal plans: one with and one without the edge penalty applied. Films were analyzed using a triple-channel nonuniformity correction and measurements were compared directly to calculations. RESULTS: Apertures generated with the edge penalty were larger, more regularly shaped and required up to 30% fewer monitor units than those created without the edge penalty. Dose volume histogram analysis showed that the changes in doses to targets, organs at risk, and normal tissues were negligible. Edge penalty apertures that were measured with film for the paraspinal plan showed a notable decrease in the number of pixels disagreeing with calculation by more than 10%. For a 5% dose passing criterion, the number of pixels passing in the composite dose distributions for the non-edge penalty and edge penalty plans were 52% and 96%, respectively. Employing gamma with 3% dose/1 mm distance criteria resulted in a 79.5% (without penalty)/95.4% (with penalty) pass rate for the two plans. Gradient compensation of 3%/1 mm resulted in 83.3%/96.2% pass rates. CONCLUSIONS: The use of the edge penalty during optimization has the potential to markedly improve dose delivery accuracy for VMAT plans while still maintaining high quality optimized dose distributions. The penalty regularizes aperture shape and improves delivery efficiency.

Inverse-optimized 3D conformal planning: minimizing complexity while achieving equivalence with beamlet IMRT in multiple clinical sites.

PURPOSE: Inverse planned intensity modulated radiation therapy (IMRT) has helped many centers implement highly conformal treatment planning with beamlet-based techniques. The many comparisons between IMRT and 3D conformal (3DCRT) plans, however, have been limited because most 3DCRT plans are forward-planned while IMRT plans utilize inverse planning, meaning both optimization and delivery techniques are different. This work avoids that problem by comparing 3D plans generated with a unique inverse planning method for 3DCRT called inverse-optimized 3D (IO-3D) conformal planning. Since IO-3D and the beamlet IMRT to which it is compared use the same optimization techniques, cost functions, and plan evaluation tools, direct comparisons between IMRT and simple, optimized IO-3D plans are possible. Though IO-3D has some similarity to direct aperture optimization (DAO), since it directly optimizes the apertures used, IO-3D is specifically designed for 3DCRT fields (i.e., 1-2 apertures per beam) rather than starting with IMRT-like modulation and then optimizing aperture shapes. The two algorithms are very different in design, implementation, and use. The goals of this work include using IO-3D to evaluate how close simple but optimized IO-3D plans come to nonconstrained beamlet IMRT, showing that optimization, rather than modulation, may be the most important aspect of IMRT (for some sites). METHODS: The IO-3D dose calculation and optimization functionality is integrated in the in-house 3D planning/optimization system. New features include random point dose calculation distributions, costlet and cost function capabilities, fast dose volume histogram (DVH) and plan evaluation tools, optimization search strategies designed for IO-3D, and an improved, reimplemented edge/octree calculation algorithm. The IO-3D optimization, in distinction to DAO, is designed to optimize 3D conformal plans (one to two segments per beam) and optimizes MLC segment shapes and weights with various user-controllable search strategies which optimize plans without beamlet or pencil beam approximations. IO-3D allows comparisons of beamlet, multisegment, and conformal plans optimized using the same cost functions, dose points, and plan evaluation metrics, so quantitative comparisons are straightforward. Here, comparisons of IO-3D and beamlet IMRT techniques are presented for breast, brain, liver, and lung plans. RESULTS: IO-3D achieves high quality results comparable to beamlet IMRT, for many situations. Though the IO-3D plans have many fewer degrees of freedom for the optimization, this work finds that IO-3D plans with only one to two segments per beam are dosimetrically equivalent (or nearly so) to the beamlet IMRT plans, for several sites. IO-3D also reduces plan complexity significantly. Here, monitor units per fraction (MU/Fx) for IO-3D plans were 22%-68% less than that for the 1 cm × 1 cm beamlet IMRT plans and 72%-84% than the 0.5 cm × 0.5 cm beamlet IMRT plans. CONCLUSIONS: The unique IO-3D algorithm illustrates that inverse planning can achieve high quality 3D conformal plans equivalent (or nearly so) to unconstrained beamlet IMRT plans, for many sites. IO-3D thus provides the potential to optimize flat or few-segment 3DCRT plans, creating less complex optimized plans which are efficient and simple to deliver. The less complex IO-3D plans have operational advantages for scenarios including adaptive replanning, cases with interfraction and intrafraction motion, and pediatric patients.

Evaluation of multiple breathing states using a multiple instance geometry approximation (MIGA) in inverse-planned optimization for locoregional breast treatment.

PURPOSE: Although previous work demonstrated superior dose distributions for left-sided breast cancer patients planned for intensity-modulated radiation therapy (IMRT) at deep inspiration breath hold compared with conventional techniques with free-breathing, such techniques are not always feasible to limit the impact of respiration on treatment delivery. This study assessed whether optimization based on multiple instance geometry approximation (MIGA) could derive an IMRT plan that is less sensitive to known respiratory motions. METHODS AND MATERIALS: CT scans were acquired with an active breathing control device at multiple breath-hold states. Three inverse optimized plans were generated for eight left-sided breast cancer patients: one static IMRT plan optimized at end exhale, two (MIGA) plans based on a MIGA representation of normal breathing, and a MIGA representation of deep breathing, respectively. Breast and nodal targets were prescribed 52.2 Gy, and a simultaneous tumor bed boost was prescribed 60 Gy. RESULTS: With normal breathing, doses to the targets, heart, and left anterior descending (LAD) artery were equivalent whether optimizing with MIGA or on a static data set. When simulating motion due to deep breathing, optimization with MIGA appears to yield superior tumor-bed coverage, decreased LAD mean dose, and maximum heart and LAD dose compared with optimization on a static representation. CONCLUSIONS: For left-sided breast-cancer patients, inverse-based optimization accounting for motion due to normal breathing may be similar to optimization on a static data set. However, some patients may benefit from accounting for deep breathing with MIGA with improvements in tumor-bed coverage and dose to critical structures.

Adaptive diffusion smoothing: a diffusion-based method to reduce IMRT field complexity.

Inverse-planned intensity modulated radiation therapy (IMRT) is often able to achieve complex treatment planning goals that are unattainable with forward three-dimensional (3D) conformal planning. However, the common use of IMRT has introduced several new challenges. The potentially high degree of modulation in IMRT beams risks the loss of some advantages of 3D planning, such as excellent target coverage and high delivery efficiency. Previous attempts to reduce beam complexity by smoothing often result in plan degradation because the smoothing algorithm cannot distinguish between areas of desirable and undesirable modulation. The purpose of this work is to introduce and evaluate adaptive diffusion smoothing (ADS), a novel procedure designed to preferentially reduce IMRT beam complexity. In this method, a discrete diffusion equation is used to smooth IMRT beams using diffusion coefficients, automatically defined for each beamlet, that dictate the degree of smoothing allowed for each beamlet. This yields a method that can distinguish between areas of desirable and undesirable modulation. The ADS method has been incorporated into our optimization system as a weighted cost function penalty, with two diffusion coefficient definitions designed to promote: (1) uniform smoothing everywhere or (2) smoothing based on cost function gradients with respect to the plan beamlet intensities. The ADS method (with both coefficient types) has been tested in a phantom and in two clinical examples (prostate and head/neck). Both types of diffusion coefficients produce plans with reduced modulation and minimal dosimetric impact, but the cost function gradient-based coefficients show more potential for reducing beam modulation without affecting dosimetric plan quality. In summary, adaptive diffusion smoothing is a promising tool for ensuring that only the necessary amount of beam modulation is used, promoting more efficient and accurate IMRT planning, QA, and delivery.

Lexicographic ordering: intuitive multicriteria optimization for IMRT.

Optimization problems in IMRT inverse planning are inherently multicriterial since they involve multiple planning goals for targets and their neighbouring critical tissue structures. Clinical decisions are generally required, based on tradeoffs among these goals. Since the tradeoffs cannot be quantitatively determined prior to optimization, the decision-making process is usually indirect and iterative, requiring many repetitive optimizations. This situation becomes even more challenging for cases with a large number of planning goals. To address this challenge, a multicriteria optimization strategy called lexicographic ordering (LO) has been implemented and evaluated for IMRT planning. The LO approach is a hierarchical method in which the planning goals are categorized into different priority levels and a sequence of sub-optimization problems is solved in order of priority. This prioritization concept is demonstrated using two clinical cases (a simple prostate case and a relatively complex head and neck case). In addition, a unique feature of LO in a decision support role is discussed. We demonstrate that a comprehensive list of planning goals (e.g., approximately 23 for the head and neck case) can be optimized using only a few priority levels. Tradeoffs between different levels have been successfully prohibited using the LO method, making the large size problem representations simpler and more manageable. Optimization time needed for each level was practical, ranging from approximately 26 s to approximately 217 s. Using prioritization, the LO approach mimics the mental process often used by physicians as they make decisions handling the various conflicting planning goals. This method produces encouraging results for difficult IMRT planning cases in a highly intuitive manner.

How extensive of a 4D dataset is needed to estimate cumulative dose distribution plan evaluation metrics in conformal lung therapy?

The purpose of this study was to investigate the number of intermediate states required to adequately approximate the clinically relevant cumulative dose to deforming/moving thoracic anatomy in four-dimensional (4D) conformal radiotherapy that uses 6 MV photons to target tumors. Four patients were involved in this study. For the first three patients, computed tomography images acquired at exhale and inhale were available; they were registered using B-spline deformation model and the computed transformation was further used to simulate intermediate states between exhale and inhale. For the fourth patient, 4D-acquired, phase-sorted datasets were available and each dataset was registered with the exhale dataset. The exhale-inhale transformation was also used to simulate intermediate states in order to compare the cumulative doses computed using the actual and the simulated datasets. Doses to each state were calculated using the Dose Planning Method (DPM) Monte Carlo code and dose was accumulated for scoring on the exhale anatomy via the transformation matrices for each state and time weighting factors. Cumulative doses were estimated using increasing numbers of intermediate states and compared to simpler scenarios such as a "2-state" model which used only the exhale and inhale datasets or the dose received during the average phase of the breathing cycle. Dose distributions for each modeled state as well as the cumulative doses were assessed using dose volume histograms and several treatment evaluation metrics such as mean lung dose, normal tissue complication probability, and generalized uniform dose. Although significant "point dose" differences can exist between each breathing state, the differences decrease when cumulative doses are considered, and can become less significant yet in terms of evaluation metrics depending upon the clinical end point. This study suggests that for certain "clinical" end points of importance for lung cancer, satisfactory predictions of accumulated total dose to be received by the distorting anatomy can be achieved by calculating the dose to but a few (or even simply the average) phases of the breathing cycle.

Reporting and analyzing statistical uncertainties in Monte Carlo-based treatment planning.

PURPOSE: To investigate methods of reporting and analyzing statistical uncertainties in doses to targets and normal tissues in Monte Carlo (MC)-based treatment planning. METHODS AND MATERIALS: Methods for quantifying statistical uncertainties in dose, such as uncertainty specification to specific dose points, or to volume-based regions, were analyzed in MC-based treatment planning for 5 lung cancer patients. The effect of statistical uncertainties on target and normal tissue dose indices was evaluated. The concept of uncertainty volume histograms for targets and organs at risk was examined, along with its utility, in conjunction with dose volume histograms, in assessing the acceptability of the statistical precision in dose distributions. The uncertainty evaluation tools were extended to four-dimensional planning for application on multiple instances of the patient geometry. All calculations were performed using the Dose Planning Method MC code. RESULTS: For targets, generalized equivalent uniform doses and mean target doses converged at 150 million simulated histories, corresponding to relative uncertainties of less than 2% in the mean target doses. For the normal lung tissue (a volume-effect organ), mean lung dose and normal tissue complication probability converged at 150 million histories despite the large range in the relative organ uncertainty volume histograms. For "serial" normal tissues such as the spinal cord, large fluctuations exist in point dose relative uncertainties. CONCLUSIONS: The tools presented here provide useful means for evaluating statistical precision in MC-based dose distributions. Tradeoffs between uncertainties in doses to targets, volume-effect organs, and "serial" normal tissues must be considered carefully in determining acceptable levels of statistical precision in MC-computed dose distributions.