Parallel beamlet dose calculation via beamlet contexts in a distributed multi-GPU framework

PURPOSE

Dose calculation is one of the most computationally intensive, yet essential tasks in the treatment planning process. With the recent interest in automatic beam orientation and arc trajectory optimization techniques, there is a great need for more efficient model-based dose calculation algorithms that can accommodate hundreds to thousands of beam candidates at once. Foundational work has shown the translation of dose calculation algorithms to graphical processing units (GPUs), lending to remarkable gains in processing efficiency. But these methods provide parallelization of dose for only a single beamlet, serializing the calculation of multiple beamlets and under-utilizing the potential of modern GPUs. In this paper, the authors propose a framework enabling parallel computation of many beamlet doses using a novel beamlet context transformation and further embed this approach in a scalable network of multi-GPU computational nodes.

METHODS

The proposed context-based transformation separates beamlet-local density and TERMA into distinct beamlet contexts that independently provide sufficient data for beamlet dose calculation. Beamlet contexts are arranged in a composite context array with dosimetric isolation, and the context array is subjected to a GPU collapsed-cone convolution superposition procedure, producing the set of beamlet-specific dose distributions in a single pass. Dose from each context is converted to a sparse representation for efficient storage and retrieval during treatment plan optimization. The context radius is a new parameter permitting flexibility between the speed and fidelity of the dose calculation process. A distributed manager-worker architecture is constructed around the context-based GPU dose calculation approach supporting an arbitrary number of worker nodes and resident GPUs. Phantom experiments were executed to verify the accuracy of the context-based approach compared to Monte Carlo and a reference CPU-CCCS implementation for single beamlets and broad beams composed by addition of beamlets. Dose for representative 4π beam sets was calculated in lung and prostate cases to compare its efficiency with that of an existing beamlet-sequential GPU-CCCS implementation. Code profiling was also performed to evaluate the scalability of the framework across many networked GPUs.

RESULTS

The dosimetric accuracy of the context-based method displays <1.35% and 2.35% average error from the existing serialized CPU-CCCS algorithm and Monte Carlo simulation for beamlet-specific PDDs in water and slab phantoms, respectively. The context-based method demonstrates substantial speedup of up to two orders of magnitude over the beamlet-sequential GPU-CCCS method in the tested configurations. The context-based framework demonstrates near linear scaling in the number of distributed compute nodes and GPUs employed, indicating that it is flexible enough to meet the performance requirements of most users by simply increasing the hardware utilization.

CONCLUSIONS

The context-based approach demonstrates a new expectation of performance for beamlet-based dose calculation methods. This approach has been successful in accelerating the dose calculation process for very large-scale treatment planning problems - such as automatic 4π IMRT beam orientation and VMAT arc trajectory selection, with hundreds of thousands of beamlets - in clinically feasible timeframes. The flexibility of this framework makes it as a strong candidate for use in a variety of other very large-scale treatment planning tasks and clinical workflows.

Robust beam orientation optimization for intensity-modulated proton therapy

PURPOSE

Dose conformality and robustness are equally important in intensity modulated proton therapy (IMPT). Despite the obvious implication of beam orientation on both dosimetry and robustness, an automated, robust beam orientation optimization algorithm has not been incorporated due to the problem complexity and paramount computational challenge. In this study, we developed a novel IMPT framework that integrates robust beam orientation optimization (BOO) and robust fluence map optimization (FMO) in a unified framework.

METHODS

The unified framework is formulated to include a dose fidelity term, a heterogeneity-weighted group sparsity term, and a sensitivity regularization term. The L2, 1/2-norm group sparsity is used to reduce the number of active beams from the initial 1162 evenly distributed noncoplanar candidate beams, to between two and four. A heterogeneity index, which evaluates the lateral tissue heterogeneity of a beam, is used to weigh the group sparsity term. With this index, beams more resilient to setup uncertainties are encouraged. There is a symbiotic relationship between the heterogeneity index and the sensitivity regularization; the integrated optimization framework further improves beam robustness against both range and setup uncertainties. This Sensitivity regularization and Heterogeneity weighting based BOO and FMO framework (SHBOO-FMO) was tested on two skull-base tumor (SBT) patients and two bilateral head-and-neck (H&N) patients. The conventional CTV-based optimized plans (Conv) with SHBOO-FMO beams (SHBOO-Conv) and manual beams (MAN-Conv) were compared to investigate the beam robustness of the proposed method. The dosimetry and robustness of SHBOO-FMO plan were compared against the manual beam plan with CTV-based voxel-wise worst-case scenario approach (MAN-WC).

RESULTS

With SHBOO-FMO method, the beams with superior range robustness over manual beams were selected while the setup robustness was maintained or improved. On average, the lowest [D95%, V95%, V100%] of CTV were increased from [93.85%, 91.06%, 70.64%] in MAN-Conv plans, to [98.62%, 98.61%, 96.17%] in SHBOO-Conv plans with range uncertainties. With setup uncertainties, the average lowest [D98%, D95%, V95%, V100%] of CTV were increased from [92.06%, 94.83%, 94.31%, 78.93%] in MAN-Conv plans, to [93.54%, 96.61%, 97.01%, 91.98%] in SHBOO-Conv plans. Compared with the MAN-WC plans, the final SHBOO-FMO plans achieved comparable plan robustness and better OAR sparing, with an average reduction of [Dmean, Dmax] of [6.31, 6.55] GyRBE for the SBT cases and [1.89, 5.08] GyRBE for the H&N cases from the MAN-WC plans.

CONCLUSION

We developed a novel method to integrate robust BOO and robust FMO into IMPT optimization for a unified solution of both BOO and FMO, generating plans with superior dosimetry and good robustness.

Shape constrained fully convolutional DenseNet with adversarial training for multiorgan segmentation on head and neck CT and low-field MR images

PURPOSE

Image-guided radiotherapy provides images not only for patient positioning but also for online adaptive radiotherapy. Accurate delineation of organs-at-risk (OARs) on Head and Neck (H&N) CT and MR images is valuable to both initial treatment planning and adaptive planning, but manual contouring is laborious and inconsistent. A novel method based on the generative adversarial network (GAN) with shape constraint (SC-GAN) is developed for fully automated H&N OARs segmentation on CT and low-field MRI.

METHODS AND MATERIAL

A deep supervised fully convolutional DenseNet is employed as the segmentation network for voxel-wise prediction. A convolutional neural network (CNN)-based discriminator network is then utilized to correct predicted errors and image-level inconsistency between the prediction and ground truth. An additional shape representation loss between the prediction and ground truth in the latent shape space is integrated into the segmentation and adversarial loss functions to reduce false positivity and constrain the predicted shapes. The proposed segmentation method was first benchmarked on a public H&N CT database including 32 patients, and then on 25 0.35T MR images obtained from an MR-guided radiotherapy system. The OARs include brainstem, optical chiasm, larynx (MR only), mandible, pharynx (MR only), parotid glands (both left and right), optical nerves (both left and right), and submandibular glands (both left and right, CT only). The performance of the proposed SC-GAN was compared with GAN alone and GAN with the shape constraint (SC) but without the DenseNet (SC-GAN-ResNet) to quantify the contributions of shape constraint and DenseNet in the deep neural network segmentation.

RESULTS

The proposed SC-GAN slightly but consistently improve the segmentation accuracy on the benchmark H&N CT images compared with our previous deep segmentation network, which outperformed other published methods on the same or similar CT H&N dataset. On the low-field MR dataset, the following average Dice's indices were obtained using improved SC-GAN: 0.916 (brainstem), 0.589 (optical chiasm), 0.816 (mandible), 0.703 (optical nerves), 0.799 (larynx), 0.706 (pharynx), and 0.845 (parotid glands). The average surface distances ranged from 0.68 mm (brainstem) to 1.70 mm (larynx). The 95% surface distance ranged from 1.48 mm (left optical nerve) to 3.92 mm (larynx). Compared with CT, using 95% surface distance evaluation, the automated segmentation accuracy is higher on MR for the brainstem, optical chiasm, optical nerves and parotids, and lower for the mandible. The SC-GAN performance is superior to SC-GAN-ResNet, which is more accurate than GAN alone on both the CT and MR datasets. The segmentation time for one patient is 14 seconds using a single GPU.

CONCLUSION

The performance of our previous shape constrained fully CNNs for H&N segmentation is further improved by incorporating GAN and DenseNet. With the novel segmentation method, we showed that the low-field MR images acquired on a MR-guided radiation radiotherapy system can support accurate and fully automated segmentation of both bony and soft tissue OARs for adaptive radiotherapy.

Performance Comparison of Knowledge-Based Dose Prediction Techniques Based on Limited Patient Data

Purpose:

The accuracy of dose prediction is essential for knowledge-based planning and automated planning techniques. We compare the dose prediction accuracy of 3 prediction methods including statistical voxel dose learning, spectral regression, and support vector regression based on limited patient training data.

Methods:

Statistical voxel dose learning, spectral regression, and support vector regression were used to predict the dose of noncoplanar intensity-modulated radiation therapy (4π) and volumetric-modulated arc therapy head and neck, 4π lung, and volumetric-modulated arc therapy prostate plans. Twenty cases of each site were used for k-fold cross-validation, with k = 4. Statistical voxel dose learning bins voxels according to their Euclidean distance to the planning target volume and uses the median to predict the dose of new voxels. Distance to the planning target volume, polynomial combinations of the distance components, planning target volume, and organ at risk volume were used as features for spectral regression and support vector regression. A total of 28 features were included. Principal component analysis was performed on the input features to test the effect of dimension reduction. For the coplanar volumetric-modulated arc therapy plans, separate models were trained for voxels within the same axial slice as planning target volume voxels and voxels outside the primary beam. The effect of training separate models for each organ at risk compared to all voxels collectively was also tested. The mean squared error was calculated to evaluate the voxel dose prediction accuracy.

Results:

Statistical voxel dose learning using separate models for each organ at risk had the lowest root mean squared error for all sites and modalities: 3.91 Gy (head and neck 4π), 3.21 Gy (head and neck volumetric-modulated arc therapy), 2.49 Gy (lung 4π), and 2.35 Gy (prostate volumetric-modulated arc therapy). Compared to using the original features, principal component analysis reduced the 4π prediction error for head and neck spectral regression (−43.9%) and support vector regression (−42.8%) and lung support vector regression (−24.4%) predictions. Principal component analysis was more effective in using all/most of the possible principal components. Separate organ at risk models were more accurate than training on all organ at risk voxels in all cases.

Conclusion:

Compared with more sophisticated parametric machine learning methods with dimension reduction, statistical voxel dose learning is more robust to patient variability and provides the most accurate dose prediction method.

Robust optimization for intensity-modulated proton therapy with soft spot sensitivity regularization

PURPOSE

Proton dose distribution is sensitive to uncertainties in range estimation and patient positioning. Currently, the proton robustness is managed by worst-case scenario optimization methods, which are computationally inefficient. To overcome these challenges, we develop a novel intensity-modulated proton therapy (IMPT) optimization method that integrates dose fidelity with a sensitivity term that describes dose perturbation as the result of range and positioning uncertainties.

METHODS

In the integrated optimization framework, the optimization cost function is formulated to include two terms: a dose fidelity term and a robustness term penalizing the inner product of the scanning spot sensitivity and intensity. The sensitivity of an IMPT scanning spot to perturbations is defined as the dose distribution variation induced by range and positioning errors. To evaluate the sensitivity, the spatial gradient of the dose distribution of a specific spot is first calculated. The spot sensitivity is then determined by the total absolute value of the directional gradients of all affected voxels. The fast iterative shrinkage-thresholding algorithm is used to solve the optimization problem. This method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The proposed sensitivity-regularized method (SenR) was implemented on both clinic target volume (CTV) and planning target volume (PTV). They were compared with conventional PTV-based optimization method (Conv) and CTV-based voxel-wise worst-case scenario optimization approach (WC).

RESULTS

Under the nominal condition without uncertainties, the three methods achieved similar CTV dose coverage, while the CTV-based SenR approach better spared organs at risks (OARs) compared with the WC approach, with an average reduction of [Dmean, Dmax] of [4.72, 3.38]  GyRBE for the SBT cases and [2.54, 3.33] GyRBE for the H&N cases. The OAR sparing of the PTV-based SenR method was comparable with the WC method. The WC method, and SenR approaches all improved the plan robustness from the conventional PTV-based method. On average, under range uncertainties, the lowest [D95%, V95%, V100%] of CTV were increased from [93.75%, 88.47%, 47.37%] in the Conv method, to [99.28%, 99.51%, 86.64%] in the WC method, [97.71%, 97.85%, 81.65%] in the SenR-CTV method and [98.77%, 99.30%, 85.12%] in the SenR-PTV method, respectively. Under setup uncertainties, the average lowest [D95%, V95%, V100%] of CTV were increased from [95.35%, 94.92%, 65.12%] in the Conv method, to [99.43%, 99.63%, 87.12%] in the WC method, [96.97%, 97.13%, 77.86%] in the SenR-CTV method, and [98.21%, 98.34%, 83.88%] in the SenR-PTV method, respectively. The runtime of the SenR optimization is eight times shorter than that of the voxel-wise worst-case method.

CONCLUSION

We developed a novel computationally efficient robust optimization method for IMPT. The robustness is calculated as the spot sensitivity to both range and shift perturbations. The dose fidelity term is then regularized by the sensitivity term for the flexibility and trade-off between the dosimetry and the robustness. In the stress test, SenR is more resilient to unexpected uncertainties. These advantages in combination with its fast computation time make it a viable candidate for clinical IMPT planning.

Fully automatic multi-organ segmentation for head and neck cancer radiotherapy using shape representation model constrained fully convolutional neural networks

PURPOSE

Intensity modulated radiation therapy (IMRT) is commonly employed for treating head and neck (H&N) cancer with uniform tumor dose and conformal critical organ sparing. Accurate delineation of organs-at-risk (OARs) on H&N CT images is thus essential to treatment quality. Manual contouring used in current clinical practice is tedious, time-consuming, and can produce inconsistent results. Existing automated segmentation methods are challenged by the substantial inter-patient anatomical variation and low CT soft tissue contrast. To overcome the challenges, we developed a novel automated H&N OARs segmentation method that combines a fully convolutional neural network (FCNN) with a shape representation model (SRM).

METHODS

Based on manually segmented H&N CT, the SRM and FCNN were trained in two steps: (a) SRM learned the latent shape representation of H&N OARs from the training dataset; (b) the pre-trained SRM with fixed parameters were used to constrain the FCNN training. The combined segmentation network was then used to delineate nine OARs including the brainstem, optic chiasm, mandible, optical nerves, parotids, and submandibular glands on unseen H&N CT images. Twenty-two and 10 H&N CT scans provided by the Public Domain Database for Computational Anatomy (PDDCA) were utilized for training and validation, respectively. Dice similarity coefficient (DSC), positive predictive value (PPV), sensitivity (SEN), average surface distance (ASD), and 95% maximum surface distance (95%SD) were calculated to quantitatively evaluate the segmentation accuracy of the proposed method. The proposed method was compared with an active appearance model that won the 2015 MICCAI H&N Segmentation Grand Challenge based on the same dataset, an atlas method and a deep learning method based on different patient datasets.

RESULTS

An average DSC = 0.870 (brainstem), DSC = 0.583 (optic chiasm), DSC = 0.937 (mandible), DSC = 0.653 (left optic nerve), DSC = 0.689 (right optic nerve), DSC = 0.835 (left parotid), DSC = 0.832 (right parotid), DSC = 0.755 (left submandibular), and DSC = 0.813 (right submandibular) were achieved. The segmentation results are consistently superior to the results of atlas and statistical shape based methods as well as a patch-wise convolutional neural network method. Once the networks are trained off-line, the average time to segment all 9 OARs for an unseen CT scan is 9.5 s.

CONCLUSION

Experiments on clinical datasets of H&N patients demonstrated the effectiveness of the proposed deep neural network segmentation method for multi-organ segmentation on volumetric CT scans. The accuracy and robustness of the segmentation were further increased by incorporating shape priors using SMR. The proposed method showed competitive performance and took shorter time to segment multiple organs in comparison to state of the art methods.

A novel optimization framework for VMAT with dynamic gantry couch rotation

Existing volumetric modulated arc therapy (VMAT) optimization using coplanar arcs is highly efficient but usually dosimetrically inferior to intensity modulated radiation therapy (IMRT) with optimized non-coplanar beams. To achieve both dosimetric quality and delivery efficiency, we proposed in this study, a novel integrated optimization method for non-coplanar VMAT (4πVMAT). 4πVMAT with direct aperture optimization (DAO) was achieved by utilizing a least square dose fidelity objective, along with an anisotropic total variation term for regularizing the fluence smoothness, a single segment term for imposing simple apertures, and a group sparsity term for selecting beam angles. Continuous gantry/couch angle trajectories were selected using the Dijkstra's algorithm, where the edge and node costs were determined based on the maximal gantry rotation speed and the estimated fluence map at the current iteration, respectively. The couch-gantry-patient collision space was calculated based on actual machine geometry and a human subject 3D surface. Beams leading to collision are excluded from the DAO and beam trajectory selection (BTS). An alternating optimization strategy was implemented to solve the integrated DAO and BTS problem. The feasibility of 4πVMAT using one full-arc or two full-arcs was tested on nine patients with brain, lung, or prostate cancer. The plan was compared against a coplanar VMAT (2πVMAT) plan using one additional arc and collimator rotation. Compared to 2πVMAT, 4πVMAT reduced the average maximum and mean organs-at-risk dose by 9.63% and 3.08% of the prescription dose with the same target coverage. R50 was reduced by 23.0%. Maximum doses to the dose limiting organs, such as the brainstem, the major vessels, and the proximal bronchus, were reduced by 8.1 Gy (64.8%), 16.3 Gy (41.5%), and 19.83 Gy (55.5%), respectively. The novel 4πVMAT approach affords efficient delivery of non-coplanar arc trajectories that lead to dosimetric improvements compared with coplanar VMAT using more arcs.

VMAT optimization with dynamic collimator rotation

PURPOSE

Although collimator rotation is an optimization variable that can be exploited for dosimetric advantages, existing Volumetric Modulated Arc Therapy (VMAT) optimization uses a fixed collimator angle in each arc and only rotates the collimator between arcs. In this study, we develop a novel integrated optimization method for VMAT, accounting for dynamic collimator angles during the arc motion.

METHODS

Direct Aperture Optimization (DAO) for Dynamic Collimator in VMAT (DC-VMAT) was achieved by adding to the existing dose fidelity objective an anisotropic total variation term for regulating the fluence smoothness, a binary variable for forming simple apertures, and a group sparsity term for controlling collimator rotation. The optimal collimator angle for each beam angle was selected using the Dijkstra's algorithm, where the node costs depend on the estimated fluence map at the current iteration and the edge costs account for the mechanical constraints of multi-leaf collimator (MLC). An alternating optimization strategy was implemented to solve the DAO and collimator angle selection (CAS). Feasibility of DC-VMAT using one full-arc with dynamic collimator rotation was tested on a phantom with two small spherical targets, a brain, a lung and a prostate cancer patient. The plan was compared against a static collimator VMAT (SC-VMAT) plan using three full arcs with 60 degrees of collimator angle separation in patient studies.

RESULTS

With the same target coverage, DC-VMAT achieved 20.3% reduction of R50 in the phantom study, and reduced the average max and mean OAR dose by 4.49% and 2.53% of the prescription dose in patient studies, as compared with SC-VMAT. The collimator rotation co-ordinated with the gantry rotation in DC-VMAT plans for deliverability. There were 13 beam angles in the single-arc DC-VMAT plan in patient studies that requires slower gantry rotation to accommodate multiple collimator angles.

CONCLUSIONS

The novel DC-VMAT approach utilizes the dynamic collimator rotation during arc delivery. In doing so, DC-VMAT affords more sophisticated intensity modulation, alleviating the limitation previously imposed by the square beamlet from the MLC leaf thickness and achieves higher effective modulation resolution. Consequently, DC-VMAT with a single arc manages to achieve superior dosimetry than SC-VMAT with three full arcs.

Technical Note: Iterative megavoltage CT (MVCT) reconstruction using block-matching 3D-transform (BM3D) regularization

PURPOSE

Megavoltage CT (MVCT) images are noisier than kilovoltage CT (KVCT) due to low detector efficiency to high-energy x rays. Conventional denoising methods compromise edge resolution and low-contrast object visibility. In this work, we incorporated block-matching 3D-transform shrinkage (BM3D) transformation into MVCT iterative reconstruction as nonlocal patch-wise regularization.

METHODS

The iterative reconstruction was achieved by adding to the existing least square data fidelity objective a regularization term, formulated as the L1 norm of the BM3D transformed image. A Fast Iterative Shrinkage-Thresholding Algorithm (FISTA) was adopted to accelerate CT reconstruction. The proposed method was compared against total variation (TV) regularization, BM3D postprocess method, and filtered back projection (FBP).

RESULTS

In the Catphan phantom study, BM3D regularization better enhances low-contrast objects compared with TV regularization and BM3D postprocess method at the same noise level. The spatial resolution using BM3D regularization is 2.79 and 2.55 times higher than that using the TV regularization at 50% of the modulation transfer function (MTF) magnitude, for the fully sampled reconstruction and down-sampled reconstruction, respectively. The BM3D regularization images show better bony details and low-contrast soft tissues, on the head and neck (H&N) and prostate patient images.

CONCLUSIONS

The proposed iterative BM3D regularization CT reconstruction method takes advantage of both the BM3D denoising capability and iterative reconstruction data fidelity consistency. This novel approach is superior to TV regularized iterative reconstruction or BM3D postprocess for improving noisy MVCT image quality.

Response to "in regard to "Tran A, Zhang J, Woods K, Yu V, Nguyen D, Gustafson G, Rosen L, Sheng K. Treatment planning comparison of IMPT, VMAT and 4π radiotherapy for prostate cases""

In regard to our recently published paper entitled “Treatment planning comparison of IMPT, VMAT and 4π radiotherapy for prostate cases”, a question was raised whether “4π” was used appropriately to describe the non-coplanar planning and delivery space. In this letter, the term use is explained from both theoretical and practical perspectives. It is concluded that the self-explanatory term provides a flexible description of non-coplanar radiotherapy with beam orientation optimization. Confusions with this term can be avoided by understanding the evolving and machine/patient specific nature of 4π planning,

4π plan optimization for cortical-sparing brain radiotherapy

BACKGROUND AND PURPOSE

Incidental irradiation of normal brain tissue during radiotherapy is linked to cognitive decline, and may be mediated by damage to healthy cortex. Non-coplanar techniques may be used for cortical sparing. We compared normal brain sparing and probability of cortical atrophy using 4π radiation therapy planning vs. standard fixed gantry intensity-modulated radiotherapy (IMRT).

MATERIAL AND METHODS

Plans from previously irradiated brain tumor patients ("original IMRT", n = 13) were re-planned to spare cortex using both 4π optimization ("4π") and IMRT optimization ("optimized IMRT"). Homogeneity index (HI), gradient measure, doses to cortex and white matter (excluding tumor), brainstem, optics, and hippocampus were compared with matching PTV coverage. Probability of three grades of post-treatment cortical atrophy was modeled based on previously established dose response curves.

RESULTS

With matching PTV coverage, 4π significantly improved HI by 27% (p = 0.005) and gradient measure by 8% (p = 0.001) compared with optimized IMRT. 4π optimization reduced mean and equivalent uniform doses (EUD) to all standard OARs, with 14-15% reduction in hippocampal EUD (p ≤ 0.003) compared with the other two plans. 4π significantly reduced dose to fractional cortical volumes (V50, V40 and V30) compared with the original IMRT plans, and reduced cortical V30 by 7% (p = 0.008) compared with optimized IMRT. White matter EUD, mean dose, and fractional volumes V50, V40 and V30 were also significantly lower with 4π (p ≤ 0.001). With 4π, probability of grade 1, 2 and 3 cortical atrophy decreased by 12%, 21% and 26% compared with original IMRT and by 8%, 14% and 3% compared with optimized IMRT, respectively (p ≤ 0.001).

CONCLUSIONS

4π radiotherapy significantly improved cortical sparing and reduced doses to standard brain OARs, white matter, and the hippocampus. This was achieved with superior PTV dose homogeneity. Such sparing could reduce the probability of cortical atrophy that may lead to cognitive decline.

Integrated beam orientation and scanning-spot optimization in intensity-modulated proton therapy for brain and unilateral head and neck tumors

PURPOSE

Intensity-Modulated Proton Therapy (IMPT) is the state-of-the-art method of delivering proton radiotherapy. Previous research has been mainly focused on optimization of scanning spots with manually selected beam angles. Due to the computational complexity, the potential benefit of simultaneously optimizing beam orientations and spot pattern could not be realized. In this study, we developed a novel integrated beam orientation optimization (BOO) and scanning-spot optimization algorithm for intensity-modulated proton therapy (IMPT).

METHODS

A brain chordoma and three unilateral head-and-neck patients with a maximal target size of 112.49 cm³ were included in this study. A total number of 1162 noncoplanar candidate beams evenly distributed across 4π steradians were included in the optimization. For each candidate beam, the pencil-beam doses of all scanning spots covering the PTV and a margin were calculated. The beam angle selection and spot intensity optimization problem was formulated to include three terms: a dose fidelity term to penalize the deviation of PTV and OAR doses from ideal dose distribution; an L1-norm sparsity term to reduce the number of active spots and improve delivery efficiency; a group sparsity term to control the number of active beams between 2 and 4. For the group sparsity term, convex L2,1-norm and nonconvex L2,1/2-norm were tested. For the dose fidelity term, both quadratic function and linearized equivalent uniform dose (LEUD) cost function were implemented. The optimization problem was solved using the Fast Iterative Shrinkage-Thresholding Algorithm (FISTA). The IMPT BOO method was tested on three head-and-neck patients and one skull base chordoma patient. The results were compared with IMPT plans created using column generation selected beams or manually selected beams.

RESULTS

The L2,1-norm plan selected spatially aggregated beams, indicating potential degeneracy using this norm. L2,1/2-norm was able to select spatially separated beams and achieve smaller deviation from the ideal dose. In the L2,1/2-norm plans, the [mean dose, maximum dose] of OAR were reduced by an average of [2.38%, 4.24%] and[2.32%, 3.76%] of the prescription dose for the quadratic and LEUD cost function, respectively, compared with the IMPT plan using manual beam selection while maintaining the same PTV coverage. The L2,1/2 group sparsity plans were dosimetrically superior to the column generation plans as well. Besides beam orientation selection, spot sparsification was observed. Generally, with the quadratic cost function, 30%~60% spots in the selected beams remained active. With the LEUD cost function, the percentages of active spots were in the range of 35%~85%.The BOO-IMPT run time was approximately 20 min.

CONCLUSION

This work shows the first IMPT approach integrating noncoplanar BOO and scanning-spot optimization in a single mathematical framework. This method is computationally efficient, dosimetrically superior and produces delivery-friendly IMPT plans.

Fraction-variant beam orientation optimization for non-coplanar IMRT

Conventional beam orientation optimization (BOO) algorithms for IMRT assume that the same set of beam angles is used for all treatment fractions. In this paper we present a BOO formulation based on group sparsity that simultaneously optimizes non-coplanar beam angles for all fractions, yielding a fraction-variant (FV) treatment plan. Beam angles are selected by solving a multi-fraction fluence map optimization problem involving 500-700 candidate beams per fraction, with an additional group sparsity term that encourages most candidate beams to be inactive. The optimization problem is solved using the fast iterative shrinkage-thresholding algorithm. Our FV BOO algorithm is used to create five-fraction treatment plans for digital phantom, prostate, and lung cases as well as a 30-fraction plan for a head and neck case. A homogeneous PTV dose coverage is maintained in all fractions. The treatment plans are compared with fraction-invariant plans that use a fixed set of beam angles for all fractions. The FV plans reduced OAR mean dose and D 2 values on average by 3.3% and 3.8% of the prescription dose, respectively. Notably, mean OAR dose was reduced by 14.3% of prescription dose (rectum), 11.6% (penile bulb), 10.7% (seminal vesicle), 5.5% (right femur), 3.5% (bladder), 4.0% (normal left lung), 15.5% (cochleas), and 5.2% (chiasm). D 2 was reduced by 14.9% of prescription dose (right femur), 8.2% (penile bulb), 12.7% (proximal bronchus), 4.1% (normal left lung), 15.2% (cochleas), 10.1% (orbits), 9.1% (chiasm), 8.7% (brainstem), and 7.1% (parotids). Meanwhile, PTV homogeneity defined as D 95/D 5 improved from .92 to .95 (digital phantom), from .95 to .98 (prostate case), and from .94 to .97 (lung case), and remained constant for the head and neck case. Moreover, the FV plans are dosimetrically similar to conventional plans that use twice as many beams per fraction. Thus, FV BOO offers the potential to reduce delivery time for non-coplanar IMRT.

Altered resting-state voxel-level whole-brain functional connectivity in depressed Parkinson's disease

BACKGROUND

Depression is one of the most common non-motor symptoms in Parkinson's disease (PD), but its pathogenesis is still not very clear. Recently, degree centrality, a voxel-level whole-brain functional connectivity (FC) analysis of resting-state functional magnetic resonance imaging, has provided the most promising way to explore the neural network mechanisms underlying depressed PD.

METHODS

Degree centrality, voxel-wise image and clinical symptoms correlation and secondary seed-based FC analyses were performed in twenty-seven drug-naïve, early stage depressed PD patients, 27 non-depressed PD patients and 27 healthy controls (HCs) to reveal voxel-level whole-brain FC changes in depressed PD.

RESULTS

Compared with the HCs, depressed PD and non-depressed PD patients shared similar brain degree centrality abnormalities mainly in the basal ganglia, insular cortex, motor cortices, default mode network, prefrontal gyrus and the cerebellum. However, compared with non-depressed PD, depressed PD showed degree centrality abnormalities in the right middle prefrontal gyrus, anterior cingulate cortices, supplementary motor cortices and cerebellum lobule VI. The right middle prefrontal gyrus degree centrality abnormalities were correlated with the clinical depression severity, and using it as a seed, a secondary seed-based FC analysis further revealed the FC changes in the anterior cingulate cortices and the cerebellum lobule VI.

CONCLUSIONS

Our findings revealed that dysfunction in extensive brain areas were involved in depressed PD, and among these regions, the right middle prefrontal gyrus, anterior cingulate cortices and the cerebellum may pose as pathogenesis hubs underlying depressed PD.

Cochlea-sparing acoustic neuroma treatment with 4π radiation therapy

Purpose

This study investigates whether 4π noncoplanar radiation therapy can spare the cochleae and consequently potentially improve hearing preservation in patients with acoustic neuroma who are treated with radiation therapy.

Methods and materials

Clinical radiation therapy plans for 30 patients with acoustic neuroma were included (14 stereotactic radiation surgery [SRS], 6 stereotactic radiation therapy [SRT], and 10 intensity modulated radiation therapy [IMRT]). The 4π plans were created for each patient with 20 optimal beams selected using a greedy column generation method and subsequently recalculated in Eclipse for comparison. Organ-at-risk (OAR) doses, homogeneity index, conformity, and tumor control probability (TCP) were compared. Normal tissue complication probability (NTCP) was calculated for sensorineural hearing loss (SNHL) at 3 and 5 years posttreatment. The dose for each plan was then escalated to achieve 99.5% TCP.

Results

4π significantly reduced the mean dose to both cochleae by 2.0 Gy (32%) for SRS, 3.2 Gy (29%) for SRT, and 10.0 Gy (32%) for IMRT. The maximum dose to both cochleae was also reduced with 4π by 1.6 Gy (20%), 2.2 Gy (15%), and 7.1 Gy (18%) for SRS, SRT, and IMRT plans, respectively. The reductions in mean/maximum brainstem dose with 4π were also statistically significant. Mean doses to other OARs were reduced by 19% to 56% on average. 4π plans had a similar CN and TCP, with a significantly higher average homogeneity index (0.93 vs 0.92) and significantly lower average NTCP for SNHL at both 3 years (30.8% vs 40.8%) and 5 years (43.3% vs 61.7%). An average dose escalation of approximately 116% of the prescription dose achieved 99.5% TCP, which resulted in 32.6% and 43.4% NTCP for SNHL at 3 years and 46.4% and 64.7% at 5 years for 4π and clinical plans, respectively.

Conclusions

Compared with clinical planning methods, optimized 4π radiation therapy enables statistically significant sparing of the cochleae in acoustic neuroma treatment as well as lowering of other OAR doses, potentially reducing the risk of hearing loss.