Quantifying variability in radiation dose due to respiratory-induced tumor motion

NuTracker: a coordinate-based neural network representation of lung motion for intrafraction tumor tracking with various surrogates in radiotherapy

Objective. Tracking tumors and surrounding tissues in real-time is critical for reducing errors and uncertainties during radiotherapy. Existing methods are either limited by the linear representation or scale poorly with the volume resolution. To address both issues, we propose a novel coordinate-based neural network representation of lung motion to predict the instantaneous 3D volume at arbitrary spatial resolution from various surrogates: patient surface, fiducial marker, and single kV projection.Approach. The proposed model, namely NuTracker, decomposes the 4DCT into a template volume and dense displacement fields (DDFs), and uses two coordinate neural networks to predict them from spatial coordinates and surrogate states. The predicted template is spatially warped with the predicted DDF to produce the deformed volume for a given surrogate state. The nonlinear coordinate networks enable representing complex motion at infinite resolution. The decomposition allows imposing different regularizations on the spatial and temporal domains. The meta-learning and multi-task learning are used to train NuTracker across patients and tasks, so that commonalities and differences can be exploited. NuTracker was evaluated on seven patients implanted with markers using a leave-one-phase-out procedure.Main results. The 3D marker localization error is 0.66 mm on average and <1 mm at 95th-percentile, which is about 26% and 32% improvement over the predominant linear methods. The tumor coverage and image quality are improved by 5.7% and 11% in terms of dice and PSNR. The difference in the localization error for different surrogates is small and is not statistically significant. Cross-population learning and multi-task learning contribute to performance. The model tolerates surrogate drift to a certain extent.Significance. NuTracker can provide accurate estimation for entire tumor volume based on various surrogates at infinite resolution. It is of great potential to apply the coordinate network to other imaging modalities, e.g. 4DCBCT and other tasks, e.g. 4D dose calculation.

The effect of local non-thermal plasma therapy on the cancer-immunity cycle in a melanoma mouse model

Melanoma remains a deadly cancer despite significant advances in immune checkpoint blockade and targeted therapies. The incidence of melanoma is also growing worldwide, which highlights the need for novel treatment options and strategic combination of therapies. Here, we investigate non-thermal plasma (NTP), an ionized gas, as a promising, therapeutic option. In a melanoma mouse model, direct treatment of tumors with NTP results in reduced tumor burden and prolonged survival. Physical characterization of NTP treatment in situ reveals the deposited NTP energy and temperature associated with therapy response, and whole transcriptome analysis of the tumor identified several modulated pathways. NTP treatment also enhances the cancer-immunity cycle, as immune cells in both the tumor and tumor-draining lymph nodes appear more stimulated to perform their anti-cancer functions. Thus, our data suggest that local NTP therapy stimulates systemic, anti-cancer immunity. We discuss, in detail, how these fundamental insights will help direct the translation of NTP technology into the clinic and inform rational combination strategies to address the challenges in melanoma therapy.

Correlation of Optical Surface Respiratory Motion Signal and Internal Lung and Liver Tumor Motion: A Retrospective Single-Center Observational Study

Purpose: Surface-guided radiation therapy (SGRT) application has limitations. This study aimed to explore the relationship between patient characteristics and their external/internal correlation to qualitatively assess the external/internal correlation in a particular patient. Methods: Liver and lung cancer patients treated with radiotherapy in our institution were retrospectively analyzed. The external/internal correlation were calculated with Spearman correlation coefficient (SCC) and SCC after support vector regression (SVR) fitting (SCC_svr). The relationship between the external/internal correlation and magnitudes of motion of the tumor and external marker (A_i, A_e), tumor volume V_t, patient age, gender, and tumor location were explored. Results: The external/internal motions of liver and lung cancer patients were strongly correlated in the S-I direction, with mean SCC_svr values of 0.913 and 0.813. The correlation coefficients between the external/internal correlations and the patients’ characteristics (A_i, A_e, V_t, and age) were all smaller than 0.5; A_i, A_e and liver tumor volumes were positively correlated with the strength of the external/internal correlation, while lung tumor volumes and patient age were negative. The external/internal correlations in males and females were roughly equal, and the external/internal correlations in patients with peripheral lung cancers were stronger than those in patients with central lung cancers. Conclusion: The external/internal correlation shows great individual differences. The effects of A_i, A_e, V_t, and age are weakly to moderately correlated. Our results suggest the necessity of individualized assessment of patient's external/internal motion correlation prior to the application of SGRT technique for breath motion monitoring.

A hierarchical model of abdominal configuration changes extracted from golden angle radial magnetic resonance imaging

Abdominal organs are subject to a variety of physiological forces that superimpose their effects to influence local motion and configuration. These forces not only include breathing, but can also arise from cyclic antral contractions and a range of slow configuration changes. To elucidate each individual motion pattern as well as their combined effects, a hierarchical motion model was built for characterization of these 3 motion modes (characterized as deformation maps between states) using golden angle radial MR signals. Breathing motions are characterized first. Antral contraction states are then reconstructed after breathing motion-induced deformation are corrected; slow configuration change states are further extracted from breathing motion-corrected image reconstructions. The hierarchical model is established based on these multimodal states, which can be either individually shown or combined to demonstrate any arbitrary composited motion patterns. The model was evaluated using 20 MR scans acquired from 9 subjects. Poor reproducibility of breathing motions both within as well as between scan sessions was observed, with an average intra-subject difference of 1.6 cycles min-1 for average breathing frequencies of 12.0 cycles min-1. Antral contraction frequency distributions were more stable than breathing, but also presented poor reproducibility between scans with an average difference of 0.3 cycles min-1 for average frequencies of 3.2 cycles min-1. The magnitudes of motions beyond breathing were found to be significant, with 14.4 and 33.8 mm maximal motions measured from antral contraction and slow configuration changes, respectively. Hierarchical motion models have potential in multiple applications in radiotherapy, including improving the accuracy of dose delivery estimation, providing guidance for margin creation, and supporting advanced decisions and strategies for immobilization, treatment monitoring and gating.

Real-time liver tracking algorithm based on LSTM and SVR networks for use in surface-guided radiation therapy

Background

Surface-guided radiation therapy can be used to continuously monitor a patient’s surface motions during radiotherapy by a non-irradiating, noninvasive optical surface imaging technique. In this study, machine learning methods were applied to predict external respiratory motion signals and predict internal liver motion in this therapeutic context.

Methods

Seven groups of interrelated external/internal respiratory liver motion samples lasting from 5 to 6 min collected simultaneously were used as a dataset, D_v. Long short-term memory (LSTM) and support vector regression (SVR) networks were then used to establish external respiratory signal prediction models (LSTMpred/SVRpred) and external/internal respiratory motion correlation models (LSTMcorr/SVRcorr). These external prediction and external/internal correlation models were then combined into an integrated model. Finally, the LSTMcorr model was used to perform five groups of model updating experiments to confirm the necessity of continuously updating the external/internal correlation model. The root-mean-square error (RMSE), mean absolute error (MAE), and maximum absolute error (MAX_AE) were used to evaluate the performance of each model.

Results

The models established using the LSTM neural network performed better than those established using the SVR network in the tasks of predicting external respiratory signals for latency-compensation (RMSE < 0.5 mm at a latency of 450 ms) and predicting internal liver motion using external signals (RMSE < 0.6 mm). The prediction errors of the integrated model (RMSE ≤ 1.0 mm) were slightly higher than those of the external prediction and external/internal correlation models. The RMSE/MAE of the fifth model update was approximately ten times smaller than that of the first model update.

Conclusions

The LSTM networks outperform SVR networks at predicting external respiratory signals and internal liver motion because of LSTM’s strong ability to deal with time-dependencies. The LSTM-based integrated model performs well at predicting liver motion from external respiratory signals with system latencies of up to 450 ms. It is necessary to update the external/internal correlation model continuously.

Using needle orientation sensing as surrogate signal for respiratory motion estimation in percutaneous interventions

PURPOSE

To develop and evaluate an approach to estimate the respiratory-induced motion of lesions in the chest and abdomen.

MATERIALS AND METHODS

The proposed approach uses the motion of an initial reference needle inserted into a moving organ to estimate the lesion (target) displacement that is caused by respiration. The needles position is measured using an inertial measurement unit (IMU) sensor externally attached to the hub of an initially placed reference needle. Data obtained from the IMU sensor and the target motion are used to train a learning-based approach to estimate the position of the moving target. An experimental platform was designed to mimic respiratory motion of the liver. Liver motion profiles of human subjects provided inputs to the experimental platform. Variables including the insertion angle, target depth, target motion velocity and target proximity to the reference needle were evaluated by measuring the error of the estimated target position and processing time.

RESULTS

The mean error of estimation of the target position ranged between 0.86 and 1.29 mm. The processing maximum training and testing time was 5 ms which is suitable for real-time target motion estimation using the needle position sensor.

CONCLUSION

The external motion of an initially placed reference needle inserted into a moving organ can be used as a surrogate, measurable and accessible signal to estimate in real-time the position of a moving target caused by respiration; this technique could then be used to guide the placement of subsequently inserted needles directly into the target.

Generalised polynomial chaos-based uncertainty quantification for planning MRgLITT procedures

PURPOSE

A generalised polynomial chaos (gPC) method is used to incorporate constitutive parameter uncertainties within the Pennes representation of bioheat transfer phenomena. The stochastic temperature predictions of the mathematical model are critically evaluated against MR thermometry data for planning MR-guided laser-induced thermal therapies (MRgLITT).

METHODS

The Pennes bioheat transfer model coupled with a diffusion theory approximation of laser tissue interaction was implemented as the underlying deterministic kernel. A probabilistic sensitivity study was used to identify parameters that provide the most variance in temperature output. Confidence intervals of the temperature predictions are compared to MR temperature imaging (MRTI) obtained during phantom and in vivo canine (n = 4) MRgLITT experiments. The gPC predictions were quantitatively compared to MRTI data using probabilistic linear and temporal profiles as well as 2-D 60 °C isotherms.

RESULTS

Optical parameters provided the highest variance in the model output (peak standard deviation: anisotropy 3.51 °C, absorption 2.94 °C, scattering 1.84 °C, conductivity 1.43 °C, and perfusion 0.94 °C). Further, within the statistical sense considered, a non-linear model of the temperature and damage-dependent perfusion, absorption, and scattering is captured within the confidence intervals of the linear gPC method. Multivariate stochastic model predictions using parameters with the dominant sensitivities show good agreement with experimental MRTI data.

CONCLUSIONS

Given parameter uncertainties and mathematical modelling approximations of the Pennes bioheat model, the statistical framework demonstrates conservative estimates of the therapeutic heating and has potential for use as a computational prediction tool for thermal therapy planning.