Inverse 4D conformal planning for lung SBRT using particle swarm optimization

Machine Learning Radiomics Model for External and Internal Respiratory Motion Correlation Prediction in Lung Tumor

Objectives: The complexity and specificity of lung tumor motion render it necessary to determine the external and internal correlation individually before applying indirect tumor tracking. However, the correlation cannot be determined from patient respiratory and tumor clinical characteristics before treatment. The purpose of this study is to present a machine learning model for an external/internal correlation prediction that is based on computed tomography (CT) radiomic features. Methods: 4-dimensional computed tomography (4DCT) images of 67 patients were collected retrospectively, and the external/internal correlation of lung tumor was calculated based on Spearman's rank correlation coefficient. Radiomic features were extracted from average intensity projection and the light gradient boosting machine (LightGBM)-based cross-validation (the recursive elimination method) was used for feature selection. The LightGBM framework forecasting models with classification thresholds 0.7, 0.8, and 0.9 are established using stratified 5-fold cross-validation. Model performance was assessed using receiver operating characteristics, sensitivity, and specificity. Results: There were 16, 18, and 13 features selected for models 0.7, 0.8, and 0.9, respectively. Texture features are of great importance in external/internal correlation prediction compared to other features in all models. The sensitivities of the predictions in models 0.7, 0.8, and 0.9 were 0.800 ± 0.126, 0.829 ± 0.140, and 0.864 ± 0.086, respectively. The specificities were 0.771 ± 0.114, 0.936 ± 0.0581, and 0.839 ± 0.101, whereas the area under the curve (AUC) was 0.837, 0.946, and 0.877, respectively. Conclusions: Our findings indicate that radiomics is an effective tool for respiratory motion correlation prediction, which can extract tumor motion characteristics. We proposed a machine learning framework for correlation prediction in the motion management strategy for lung tumor patients.

Combining Serial and Parallel Functionality in Functional Lung Avoidance Radiation Therapy

PURPOSE

Functional lung avoidance (FLA) radiation therapy (RT) aims to minimize post-RT pulmonary toxicity by preferentially avoiding dose to high-functioning lung (HFL) regions. A common limitation is that FLA approaches do not consider the conducting architecture for gas exchange. We previously proposed the functionally weighted airway sparing (FWAS) method to spare airways connected to HFL regions, showing that it is possible to substantially reduce risk of radiation-induced airway injury. Here, we compare the performance of FLA and FWAS and propose a novel method combining both approaches.

METHODS

We used breath-hold computed tomography (BHCT) and simulation 4-dimensional computed tomography (4DCT) from 12 lung stereotactic ablative radiation therapy patients. Four planning strategies were examined: (1) Conventional: no sparing other than clinical dose-volume constraints; (2) FLA: using a 4DCT-based ventilation map to delineate the HFL, plans were optimized to reduce mean dose and V_13.50 in HFL; (3) FWAS: we autosegemented 11 to 13 generations of individual airways from each patient's BHCT and assigned priorities based on the relative contribution of each airway to total ventilation. We used these priorities in the optimization along with airway dose constraints, estimated as a function of airway diameter and 5% probability of collapse; and (4) FLA + FWAS: we combined information from the 2 strategies. We prioritized clinical dose constraints for organs at risk and planning target volume in all plans. We performed the evaluation in terms of ventilation preservation accounting for radiation-induced damage to both lung parenchyma and airways.

RESULTS

We observed average ventilation preservation for FLA, FWAS, and FLA + FWAS as 3%, 8.5%, and 14.5% higher, respectively, than for Conventional plans for patients with ventilation preservation in Conventional plans <90%. Generalized estimated equations showed that all improvements were statistically significant (P ≤ .036). We observed no clinically relevant improvements in outcomes of the sparing techniques in patients with ventilation preservation in Conventional plans ≥90%.

CONCLUSIONS

These initial results suggest that it is crucial to consider the parallel and the serial nature of the lung to improve post-radiation therapy lung function and, consequently, quality of life for patients.

A dose volume constraint (DVC) projection-based algorithm for IMPT inverse planning optimization

PURPOSE

Provide a projection-based algorithm to solve the class of optimization problems encountered in intensity modulated proton therapy (IMPT). The algorithm can handle percentage dose-volume constraints that are usually found in such problems.

METHODS

To seek a feasible solution, the automatic relaxation method was used to project the spot weight vector onto the interval defined by lower and upper bound target dose constraints. The obtained solution was optimized separately based on the objective of each OAR in addition to maximizing the minimum target dose using the bisection search method using a stopping criterion of 10 cGy. The combined weight was used in the CQ algorithm to solve the split feasibility problem but with a special projection technique due to the non-convexity of dose volume constraints. The algorithm was applied to four clinical IMPT cases (meningioma, prostate, tongue, and oropharynx) and compared to the corresponding treatment plans optimized in Eclipse.

RESULTS

The treatment plans obtained, for the four cases, using the BCQ-ARM algorithm have dosimetric endpoints that are similar to their counterparts generated from Eclipse. The algorithm worked equally well with all cases, including the complex head and neck ones. The stopping criterion of 10 cGy results in making the generated plans slightly less optimal (ε-optimal) rather than optimal, but with the advantage of the possibility of generating a database of plans.

CONCLUSIONS

The application of the BCQ-ARM algorithm to different cases of IMPT plans with dose volume constraints was demonstrated. The algorithm is successful in generating plans that are dosimetrically equivalent to their corresponding Eclipse plans. Thus, it is suitable to generate optimized treatment plans in a clinically reasonable time frame. This article is protected by copyright. All rights reserved.

AAPM Task Group 264: The safe clinical implementation of MLC tracking in radiotherapy

The era of real-time radiotherapy is upon us. Robotic and gimbaled linac tracking are clinically established technologies with the clinical realization of couch tracking in development. Multileaf collimators (MLCs) are a standard equipment for most cancer radiotherapy systems, and therefore MLC tracking is a potentially widely available technology. MLC tracking has been the subject of theoretical and experimental research for decades and was first implemented for patient treatments in 2013. The AAPM Task Group 264 Safe Clinical Implementation of MLC Tracking in Radiotherapy Report was charged to proactively provide the broader radiation oncology community with (a) clinical implementation guidelines including hardware, software, and clinical indications for use, (b) commissioning and quality assurance recommendations based on early user experience, as well as guidelines on Failure Mode and Effects Analysis, and (c) a discussion of potential future developments. The deliverables from this report include: an explanation of MLC tracking and its historical development; terms and definitions relevant to MLC tracking; the clinical benefit of, clinical experience with and clinical implementation guidelines for MLC tracking; quality assurance guidelines, including example quality assurance worksheets; a clinical decision pathway, future outlook and overall recommendations.

Functionally weighted airway sparing (FWAS): a functional avoidance method for preserving post-treatment ventilation in lung radiotherapy

Recent changes to the guidelines for screening and early diagnosis of lung cancer have increased the interest in preserving post-radiotherapy lung function. Current investigational approaches are based on spatially mapping functional regions and generating regional avoidance plans that preferentially spare highly ventilated/perfused lung. A potentially critical, yet overlooked, aspect of functional avoidance is radiation injury to peripheral airways, which serve as gas conduits to and from functional lung regions. Dose redistribution based solely on regional function may cause irreparable damage to the 'supply chain'. To address this deficiency, we propose the functionally weighted airway sparing (FWAS) method. FWAS (i) maps the bronchial pathways to each functional sub-lobar lung volume; (ii) assigns a weighting factor to each airway based on the relative contribution of the sub-volume to overall lung function; and (iii) creates a treatment plan that aims to preserve these functional pathways. To evaluate it, we used four cases from a retrospective cohort of SAbR patients treated for lung cancer. Each patient's airways were auto-segmented from a diagnostic-quality breath-hold CT using a research virtual bronchoscopy software. A ventilation map was generated from the planning 4DCT to map regional lung function. For each terminal airway, as resolved by the segmentation software, the total ventilation within the sub-lobar volume supported by that airway was estimated and used as a function-based weighting factor. Upstream airways were weighted based on the cumulative volumetric ventilation supported by corresponding downstream airways. Using a previously developed model for airway radiosensitivity, dose constraints were determined for each airway corresponding to a <5% probability of airway collapse. Airway dose constraints, ventilation scores, and clinical dose constraints were input to a swarm optimization-based inverse planning engine to create a 3D conformal SAbR plan (CRT). The FWAS plans were compared to the patients' prescribed CRT clinical plans and the inverse-optimized clinical plans. Depending on the size and location of the tumour, the FWAS plan showed superior preservation of ventilation due to airflow preservation through open pathways (i.e. cumulative ventilation score from the sub-lobar volumes of open pathways). Improvements ranged between 3% and 23%, when comparing to the prescribed clinical plans, and between 3% and 35%, when comparing to the inverse-optimized clinical plans. The three plans satisfied clinical requirements for PTV coverage and OAR dose constraints. These initial results suggest that by sparing pathways to high-functioning lung subregions it is possible to reduce post-SAbR loss of respiratory function.

Exploiting tumor position differences between deep inspiration and expiration in lung stereotactic body radiation therapy planning

PURPOSE

We demonstrate proof of principle that normal tissue doses can be greatly reduced in lung stereotactic body radiation therapy (SBRT) for mobile tumors, if the delivered dose is split between opposite respiratory states.

METHODS

Patients that underwent 5 fraction lung SBRT at our institution and had deep inspiration breath hold (DIBH) and free breathing 4D computed tomography scans were included. Volumetric modulated arc therapy plans were generated on both respiratory phases and a third composite plan was generated delivering half the dose using the DIBH plan and the other half using the expiratory phase plan for each fraction. Computed tomography scans for the composite plan were fused based on ribs adjacent to the tumor to evaluate the dose volume histogram of critical structures.

RESULTS

Four patients with 4 total tumors had requisite planning scans available. Tumor size was between 0.7 to 2.9 cm and tumor movement 1.4 to 2.9 cm. Median reduction in the chest wall (CW) V30Gy for the composite plan was 74.6% (range 33.7 to 100%), 76.9% (range 32.9 to 100%), and 89.3% (range 69.5 to 100%) compared to the DIBH, expiration phase, and free breathing plans, respectively. Median reduction in CW maximum dose for the composite plan was 23.3% (range 0.27% to 46.4%), 23.5% (range 3.2 to 48.2%), and 23.4% (range 0.27% to 48.4%) compared to the DIBH, expiration phase, and free breathing plans, respectively. Greater reduction in CW maximum dose was observed when patients had no overlap in planning target volumes between DIBH and expiration phases (median reduction 43.9% for no overlap vs 2.7% with overlap). Between all plans, lung V20Gy absolute differences were within 1.3%. For 2 of 4 patients, the composite plan met constraints for 3 fraction SBRT, while standard plans did not.

CONCLUSIONS

We conclude that composite DIBH-expiration SBRT planning has the potential to improve organ at risk sparing.

Automatic VMAT planning for post-operative prostate cancer cases using particle swarm optimization: A proof of concept study

OBJECTIVE

To investigate the potential of Particle Swarm Optimization (PSO) for fully automatic VMAT radiotherapy (RT) treatment planning.

MATERIAL AND METHODS

In PSO a solution space of planning constraints is searched for the best possible RT plan in an iterative, statistical method, optimizing a population of candidate solutions. To identify the best candidate solution and for final evaluation a plan quality score (PQS), based on dose volume histogram (DVH) parameters, was introduced. Automatic PSO-based RT planning was used for N = 10 postoperative prostate cancer cases, retrospectively taken from our clinical database, with a prescribed dose of EUD = 66 Gy in addition to two constraints for rectum and one for bladder. Resulting PSO-based plans were compared dosimetrically to manually generated VMAT plans.

RESULTS

PSO successfully proposed treatment plans comparable to manually optimized ones in 9/10 cases. The median (range) PTV EUD was 65.4 Gy (64.7-66.0) for manual and 65.3 Gy (62.5-65.5) for PSO plans, respectively. However PSO plans achieved significantly lower doses in rectum D_2% 67.0 Gy (66.5-67.5) vs. 66.1 Gy (64.7-66.5, p = 0.016). All other evaluated parameters (PTV D_98% and D_2%, rectum V_40Gy and V_60Gy, bladder D_2% and V_60Gy) were comparable in both plans. Manual plans had lower PQS compared to PSO plans with -0.82 (-16.43-1.08) vs. 0.91 (-5.98-6.25).

CONCLUSION

PSO allows for fully automatic generation of VMAT plans with plan quality comparable to manually optimized plans. However, before clinical implementation further research is needed concerning further adaptation of PSO-specific parameters and the refinement of the PQS.

Inverse radiotherapy planning based on bioeffect modelling for locally advanced left-sided breast cancer

BACKGROUND AND PURPOSE

Treatment planning of radiotherapy (RT) for left-sided breast cancer is a challenging case. Several competing concerns are incorporated at present through protocol-defined dose-volume constraints, e.g. cardiac exposure and target coverage. Such constraints are limited by neglecting patient-specific risk factors (RFs). We propose an alternative RT planning method based solely on bioeffect models to minimize the estimated risks of breast cancer recurrence (BCR) and radiation-induced mortality endpoints considering patient-specific factors.

METHODS AND MATERIALS

Thirty-nine patients with left-sided breast cancer treated with comprehensive post-lumpectomy loco-regional conformal RT were included. An in-house particle swarm optimization (PSO) engine was used to choose fields from a large set of predefined fields and optimize monitor units to minimize the total risk of BCR and mortality caused by radiation-induced ischaemic heart disease (IHD), secondary lung cancer (SLC) and secondary breast cancer (SBC). Risk models included patient age, smoking status and cardiac risk and were developed using published multi-institutional data.

RESULTS

For the clinical plans the normal tissue complication probability, i.e. summed risk of IHD, SLC and SBC, was <3.7% and the risk of BCR was <6.1% for all patients. Median total decrease in mortality or recurrence achieved with individualized PSO plans was 0.4% (range, 0.06-2.0%)/0.5% (range, 0.11-2.2%) without/with risk factors.

CONCLUSIONS

Inverse RT plan optimization using bioeffect probability models allows individualization according to patient-specific risk factors. The modelled benefit when compared to clinical plans is, however, modest in most patients, demonstrating that current clinical plans are close to optimal. Larger gains may be achievable with morbidity endpoints rather than mortality.

Adaptive method for multicriteria optimization of intensity-modulated proton therapy

PURPOSE

Provide an adaptive multicriteria optimization (MCO) method for intensity-modulated proton therapy (IMPT) utilizing GPU technology. Previously described limitations of MCO such as Pareto approximation and limitation on the number of objectives were addressed.

METHODS

The treatment planning process for IMPT must account for multiple objectives, which requires extensive treatment planning resources. Often a large number of objectives (>10) are required. Hence the need for an MCO algorithm that can handle large number of objectives. The novelty of the MCO method presented here lies on the introduction of the adaptive weighting scheme that can generate a well-distributed and dense representation of the Pareto surface for a large number of objectives in an efficient manner. In our approach the generated Pareto surface is constructed for a set of DVH objectives. The MCO algorithm is based on the augmented weighted Chebychev metric (AWCM) method with an adaptive weighting scheme. This scheme uses the differential evolution (DE) method to generate a set of well-distributed Pareto points. The quality of the Pareto points' distribution in the objective space was assessed quantitatively using the Pareto sampling metric. The MCO algorithm was developed to perform multiple parallel searches to achieve a rapid mapping of the Pareto surface, produce clinically deliverable plans, and was implemented on a GPU cluster. The MCO algorithm was tested on two clinical cases with 10 and 18 objectives. For each case one of the MCO-generated plans was selected for comparison with the clinically generated plan. The MCO plan was randomly selected out of the set of MCO plans that had target coverage similar to the clinically generated plan and the same or better sparing of the organs at risk (OAR). Additionally, a validation study of the AWCM method vs the weighted sum method was performed.

RESULTS

The adaptive MCO algorithm generated Pareto points on the Pareto hypersurface in a fast (2-3 hr) and efficient manner for 2 cases with 10 and 18 objectives. The MCO algorithm generated a dense and well-distributed set of Pareto points on the objective space, and was able to achieve minimization of the Pareto sampling metric. The selected MCO plan showed an improvement of the DVH objectives in comparison to the clinically optimized plan in both cases. For case one, the MCO plan showed a 48% reduction of the 50% dose to OARs and a 16% reduction of the 1% dose to OARs. For case 2, the MCO plan showed a 72% reduction of the 50% dose to OARs and a 42% reduction of the 1% dose to OARs. The comparison of AWCM to WS showed that the AWCM method has a dosimetric advantage over WS for both patient cases.

CONCLUSION

We introduced an adaptive MCO algorithm for IMPT accelerated using GPUs. The algorithm is based on an adaptive method for generating Pareto plans in the objective space. We have shown that the algorithm can provide rapid and efficient mapping of the multicriteria Pareto surface with clinically deliverable plans.

Inverse-planned deliverable 4D-IMRT for lung SBRT

PURPOSE

We present a particle swarm optimization (PSO)-based technique to create deliverable four-dimensional (4D = 3D + time) intensity-modulated radiation therapy (IMRT) plans for lung stereotactic body radiotherapy (SBRT). The 4D planning concept uses respiratory motion as an additional degree of freedom to achieve further sparing of organs at risk (OARs). The 4D-IMRT plan involves the delivery of an order of magnitude more IMRT apertures (~15,000-20,000), with potentially large interaperture variations in the delivered fluence, compared to conventional (i.e., 3D) IMRT. In order to deliver the 4D plan in an efficient manner, we present an optimization-based aperture sequencing technique.

METHOD

A graphic processing unit (GPU)-enabled PSO-based inverse planning engine, developed and integrated with a research version of the Eclipse (Varian, Palo Alto, CA) treatment planning system (TPS), was employed to create 4D-IMRT plans as follows. Four-dimensional computed tomography scans (4DCTs) and beam configurations from clinical treatment plans of seven lung cancer patients were retrospectively collected, and in each case, the PSO engine iteratively adjusted aperture monitor unit (MU) weights for all beam apertures across all respiratory phases to optimize OAR dose sparing while maintaining planning target volume (PTV) coverage. We calculated the transition times from each aperture to all other apertures for each beam, taking into account the maximum leaf velocity of the multileaf collimator (MLC), and developed a mixed integer optimization technique for aperture sequencing. The goal of sequencing was to maximize delivery efficiency (i.e., minimize the time required to deliver the dose map) by accounting for leaf velocity, aperture MUs, and duration of each respiratory phase. The efficiency of the proposed delivery method was compared with that of a greedy algorithm which chose only from neighboring apertures for the subsequent steps in the sequence.

RESULTS

4D-IMRT-optimized plans achieved PTV coverage comparable to clinical plans while improving OAR sparing by an average of 39.7% for D max heart, 20.5% for D max esophagus, 25.6% for D max spinal cord, and 2.1% for V 13 lung (with D max standing for maximum dose and V 13 standing for volume receiving ≥ 13 Gy). Our mixed integer optimization-based aperture sequencing enabled the delivery to be performed in fewer cycles compared to the greedy method. This reduction was 89 ± 79 cycles corresponding to an improvement of 15.94 ± 8.01%, when considering respiratory cycle duration of 4 s, and 55 ± 33 cycles corresponding to an improvement of 15.14 ± 4.45%, when considering respiratory cycle duration of 6 s.

CONCLUSION

PSO-based 4D-IMRT represents an attractive technique to further improve OAR sparing in lung SBRT. Efficient delivery of a large number of sparse apertures (control points) introduces a challenge in 4D-IMRT treatment planning and delivery. Through judicious optimization of the aperture sequence across all phases, such delivery can be performed on a clinically feasible time scale.

Virtual Bronchoscopy-Guided Treatment Planning to Map and Mitigate Radiation-Induced Airway Injury in Lung SAbR

PURPOSE

Radiation injury to the bronchial tree is an important yet poorly understood potential side effect in lung stereotactic ablative radiation therapy (SAbR). We investigate the integration of virtual bronchoscopy in radiation therapy planning to quantify dosage to individual airways. We develop a risk model of airway collapse and develop treatment plans that reduce the risk of radiation-induced airway injury.

METHODS AND MATERIALS

Pre- and post-SAbR diagnostic-quality computerized tomography (CT) scans were retrospectively collected from 26 lung cancer patients. From each scan, the bronchial tree was segmented using a virtual bronchoscopy system and registered deformably to the planning CT. Univariate and stepwise multivariate Cox regressions were performed, examining factors such as age, comorbidities, smoking pack years, airway diameter, and maximum point dosage (D_max). Logistic regression was utilized to formulate a risk function of segmental collapse based on D_max and diameter. The risk function was incorporated into the objective function along with clinical dosage volume constraints for planning target volume (PTV) and organs at risk (OARs).

RESULTS

Univariate analysis showed that segmental diameter (P = .014) and D_max (P = .007) were significantly correlated with airway segment collapse. Multivariate stepwise Cox regression showed that diameter (P = .015), D_max (P < .0001), and pack/years of smoking (P = .02) were significant independent factors associated with collapse. Risk management-based plans enabled significant dosage reduction to individual airway segments while fulfilling clinical dosimetric objectives.

CONCLUSION

To our knowledge, this is the first systematic investigation of functional avoidance in lung SAbR based on mapping and minimizing doses to individual bronchial segments. Our early results show that it is possible to substantially lower airway dosage. Such dosage reduction may potentially reduce the risk of radiation-induced airway injury, while satisfying clinically prescribed dosimetric objectives.

Internal Motion Estimation by Internal-external Motion Modeling for Lung Cancer Radiotherapy

The aim of this study is to develop an internal-external correlation model for internal motion estimation for lung cancer radiotherapy. Deformation vector fields that characterize the internal-external motion are obtained by respectively registering the internal organ meshes and external surface meshes from the 4DCT images via a recently developed local topology preserved non-rigid point matching algorithm. A composite matrix is constructed by combing the estimated internal phasic DVFs with external phasic and directional DVFs. Principle component analysis is then applied to the composite matrix to extract principal motion characteristics, and generate model parameters to correlate the internal-external motion. The proposed model is evaluated on a 4D NURBS-based cardiac-torso (NCAT) synthetic phantom and 4DCT images from five lung cancer patients. For tumor tracking, the center of mass errors of the tracked tumor are 0.8(±0.5)mm/0.8(±0.4)mm for synthetic data, and 1.3(±1.0)mm/1.2(±1.2)mm for patient data in the intra-fraction/inter-fraction tracking, respectively. For lung tracking, the percent errors of the tracked contours are 0.06(±0.02)/0.07(±0.03) for synthetic data, and 0.06(±0.02)/0.06(±0.02) for patient data in the intra-fraction/inter-fraction tracking, respectively. The extensive validations have demonstrated the effectiveness and reliability of the proposed model in motion tracking for both the tumor and the lung in lung cancer radiotherapy.

Multi-GPU configuration of 4D intensity modulated radiation therapy inverse planning using global optimization

We report on the design, implementation and characterization of a multi-graphic processing unit (GPU) computational platform for higher-order optimization in radiotherapy treatment planning. In collaboration with a commercial vendor (Varian Medical Systems, Palo Alto, CA), a research prototype GPU-enabled Eclipse (V13.6) workstation was configured. The hardware consisted of dual 8-core Xeon processors, 256 GB RAM and four NVIDIA Tesla K80 general purpose GPUs. We demonstrate the utility of this platform for large radiotherapy optimization problems through the development and characterization of a parallelized particle swarm optimization (PSO) four dimensional (4D) intensity modulated radiation therapy (IMRT) technique. The PSO engine was coupled to the Eclipse treatment planning system via a vendor-provided scripting interface. Specific challenges addressed in this implementation were (i) data management and (ii) non-uniform memory access (NUMA). For the former, we alternated between parameters over which the computation process was parallelized. For the latter, we reduced the amount of data required to be transferred over the NUMA bridge. The datasets examined in this study were approximately 300 GB in size, including 4D computed tomography images, anatomical structure contours and dose deposition matrices. For evaluation, we created a 4D-IMRT treatment plan for one lung cancer patient and analyzed computation speed while varying several parameters (number of respiratory phases, GPUs, PSO particles, and data matrix sizes). The optimized 4D-IMRT plan enhanced sparing of organs at risk by an average reduction of [Formula: see text] in maximum dose, compared to the clinical optimized IMRT plan, where the internal target volume was used. We validated our computation time analyses in two additional cases. The computation speed in our implementation did not monotonically increase with the number of GPUs. The optimal number of GPUs (five, in our study) is directly related to the hardware specifications. The optimization process took 35 min using 50 PSO particles, 25 iterations and 5 GPUs.