Data Science Opportunities To Improve Radiotherapy Planning and Clinical Decision Making

Radiotherapy aims to achieve a high tumor control probability while minimizing damage to normal tissues. Personalizing radiotherapy treatments for individual patients, therefore, depends on integrating physical treatment planning with predictive models of tumor control and normal tissue complications. Predictive models could be improved using a wide range of rich data sources, including tumor and normal tissue genomics, radiomics, and dosiomics. Deep learning will drive improvements in classifying normal tissue tolerance, predicting intra-treatment tumor changes, tracking accumulated dose distributions, and quantifying the tumor response to radiotherapy based on imaging. Mechanistic patient-specific computer simulations ('digital twins') could also be used to guide adaptive radiotherapy. Overall, we are entering an era where improved modeling methods will allow the use of newly available data sources to better guide radiotherapy treatments.

Radiation Dose-Volume-Response Relationships for Adverse Events in Childhood Cancer Survivors: Introduction to the Scientific Issues in PENTEC

At its very core, radiation oncology involves a trade-off between the benefits and risks of exposing tumors and normal tissue to relatively high doses of ionizing radiation. This trade-off is particularly critical in childhood cancer survivors (CCS), in whom both benefits and risks can be hugely consequential due to the long life expectancy if the primary cancer is controlled. Estimating the normal tissue-related risks of a specific radiation therapy plan in an individual patient relies on predictive mathematical modeling of empirical data on adverse events. The Pediatric Normal-Tissue Effects in the Clinic (PENTEC) collaborative network was formed to summarize and, when possible, to synthesize dose-volume-response relationships for a range of adverse events incident in CCS based on the literature. Normal-tissue clinical radiation biology in children is particularly challenging for many reasons: (1) Childhood malignancies are relatively uncommon-constituting approximately 1% of new incident cancers in the United States-and biologically heterogeneous, leading to many small series in the literature and large variability within and between series. This creates challenges in synthesizing data across series. (2) CCS are at an elevated risk for a range of adverse health events that are not specific to radiation therapy. Thus, excess relative or absolute risk compared with a reference population becomes the appropriate metric. (3) Various study designs and quantities to express risk are found in the literature, and these are summarized. (4) Adverse effects in CCS often occur 30, 50, or more years after therapy. This limits the information content of series with even very extended follow-up, and lifetime risk estimates are typically extrapolations that become dependent on the mathematical model used. (5) The long latent period means that retrospective dosimetry is required, as individual computed tomography-based radiation therapy plans gradually became available after 1980. (6) Many individual patient-level factors affect outcomes, including age at exposure, attained age, lifestyle exposures, health behaviors, other treatment modalities, dose, fractionation, and dose distribution. (7) Prospective databases with individual patient-level data and radiation dosimetry are being built and will facilitate advances in dose-volume-response modeling. We discuss these challenges and attempts to overcome them in the setting of PENTEC.

Direct incorporation of patient-specific efficacy and toxicity estimates in radiation therapy plan optimization

PURPOSE

Current radiation therapy (RT) treatment planning relies mainly on pre-defined dose-based objectives and constraints to develop plans that aim to control disease while limiting damage to normal tissues during treatment. These objectives and constraints are generally population-based, in that they are developed from the aggregate response of a broad patient population to radiation. However, correlations of new biologic markers and patient-specific factors to treatment efficacy and toxicity provide the opportunity to further stratify patient populations and develop a more individualized approach to RT planning. We introduce a novel intensity-modulated radiation therapy (IMRT) optimization strategy that directly incorporates patient-specific dose response models into the planning process. In this strategy, we integrate the concept of utility-based planning where the optimization objective is to maximize the predicted value of overall treatment utility, defined by the probability of efficacy (e.g., local control) minus the weighted sum of toxicity probabilities. To demonstrate the feasibility of the approach, we apply the strategy to treatment planning for non-small cell lung cancer (NSCLC) patients.

METHODS AND MATERIALS

We developed a prioritized approach to patient-specific IMRT planning. Using a commercial treatment planning system (TPS), we calculate dose based on an influence matrix of beamlet-dose contributions to regions-of-interest. Then, outside of the TPS, we hierarchically solve two optimization problems to generate optimal beamlet weights that can then be imported back to the TPS. The first optimization problem maximizes a patient's overall plan utility subject to typical clinical dose constraints. In this process, we facilitate direct optimization of efficacy and toxicity trade-off based on individualized dose-response models. After optimal utility is determined, we solve a secondary optimization problem that minimizes a conventional dose-based objective subject to the same clinical dose constraints as the first stage but with the addition of a constraint to maintain the optimal utility from the first optimization solution. We tested this method by retrospectively generating plans for five previously treated NSCLC patients and comparing the prioritized utility plans to conventional plans optimized with only dose metric objectives. To define a plan utility function for each patient, we utilized previously published correlations of dose to local control and grade 3-5 toxicities that include patient age, stage, microRNA levels, and cytokine levels, among other clinical factors.

RESULTS

The proposed optimization approach successfully generated RT plans for five NSCLC patients that improve overall plan utility based on personalized efficacy and toxicity models while accounting for clinical dose constraints. Prioritized utility plans demonstrated the largest average improvement in local control (16.6%) when compared to plans generated with conventional planning objectives. However, for some patients, the utility-based plans resulted in similar local control estimates with decreased estimated toxicity.

CONCLUSION

The proposed optimization approach, where the maximization of a patient's RT plan utility is prioritized over the minimization of standardized dose metrics, has the potential to improve treatment outcomes by directly accounting for variability within a patient population. The implementation of the utility-based objective function offers an intuitive, humanized approach to biological optimization in which planning trade-offs are explicitly optimized.

Proton therapy delivery method affects dose-averaged linear energy transfer in patients

The dosimetric advantages of proton therapy have led to its rapid proliferation in recent decades. This has been accompanied by a shift in technology from older units that deliver protons by passive scattering (PS) to newer units that increasingly use pencil-beam scanning (PBS). The biologic effectiveness of proton physical dose purportedly rises with increasing dose-weighted average linear energy transfer (LET_D). The objective of this study was to determine the extent to which proton delivery methods affect LET_D. We calculated LET_Dfrom simple, dosimetrically matched, and clinical treatment plans with TOPAS Monte-Carlo transport code. Simple treatment plans comprised single fields of PS and PBS protons in a water phantom. We performed simulations of matched and clinical treatment plans by using the treatment and anatomic data obtained from a cohort of children with craniopharyngioma who previously received PS or PBS proton therapy. We compared the distributions of LET_Dfrom PS and PBS delivery methods in clinically relevant ROIs. Wilcoxon signed-rank tests comparing single fields in water revealed that the LET_Dvalues from PBS were significantly greater than those from PS inside and outside the targeted volume (p < 0.01). Statistical tests comparing LET_D-volume histograms from matched and clinical treatment plans showed that LET_Dwas generally greater for PBS treatment plans than for PS treatment plans (p < 0.05). In conclusion, the proton delivery method affects LET_Dboth inside and outside of the target volume. These findings suggest that PBS is more biologically effective than PS. Given the rapid expansion of PBS proton therapy, future studies are needed to confirm the applicability of treatment evaluation methods developed for PS proton therapy to those for modern PBS treatments to ensure their safety and effectiveness for the growing population of patients receiving proton therapy. This study uses data from two clinical trials: NCT01419067 and NCT02792582.

Generalized approach for radiotherapy treatment planning by optimizing projected health outcome: preliminary results for prostate radiotherapy patients

Research in cancer care increasingly focuses on survivorship issues, e.g. managing disease- and treatment-related morbidity and mortality occurring during and after treatment. This necessitates innovative approaches that consider treatment side effects in addition to tumor cure. Current treatment-planning methods rely on constrained iterative optimization of dose distributions as a surrogate for health outcomes. The goal of this study was to develop a generally applicable method to directly optimize projected health outcomes. We developed an outcome-based objective function to guide selection of the number, angle, and relative fluence weight of photon and proton radiotherapy beams in a sample of ten prostate-cancer patients by optimizing the projected health outcome. We tested whether outcome-optimized radiotherapy (OORT) improved the projected longitudinal outcome compared to dose-optimized radiotherapy (DORT) first for a statistically significant majority of patients, then for each individual patient. We assessed whether the results were influenced by the selection of treatment modality, late-risk model, or host factors. The results of this study revealed that OORT was superior to DORT. Namely, OORT maintained or improved the projected health outcome of photon- and proton-therapy treatment plans for all ten patients compared to DORT. Furthermore, the results were qualitatively similar across three treatment modalities, six late-risk models, and 10 patients. The major finding of this work was that it is feasible to directly optimize the longitudinal (i.e. long- and short-term) health outcomes associated with the total (i.e. therapeutic and stray) absorbed dose in all of the tissues (i.e. healthy and diseased) in individual patients. This approach enables consideration of arbitrary treatment factors, host factors, health endpoints, and times of relevance to cancer survivorship. It also provides a simpler, more direct approach to realizing the full beneficial potential of cancer radiotherapy.

Justification and optimization of radiation exposures: a new framework to aggregate arbitrary detriments and benefits

Myriad radiation effects, including benefits and detriments, complicate justifying and optimizing radiation exposures. The purpose of this study was to develop a comprehensive conceptual framework and corresponding quantitative methods to aggregate the detriments and benefits of radiation exposures to individuals, groups, and populations. In this study, concepts from the ICRP for low dose were integrated with clinical techniques focused on high dose to develop a comprehensive figure of merit (FOM) that takes into account arbitrary host- and exposure-related factors, endpoints, and time points. The study built on existing methods with three new capabilities: application to individuals, groups, and populations; extension to arbitrary numbers and types of endpoints; and inclusion of limitation, where relevant. The FOM was applied to three illustrative exposure situations: emergency response, diagnostic imaging, and cancer radiotherapy, to evaluate its utility in diverse settings. The example application to radiation protection revealed the FOM's utility in optimizing the benefits and risks to a population while keeping individual exposures below applicable regulatory limits. Examples in diagnostic imaging and cancer radiotherapy demonstrated the FOM's utility for guiding population- and patient-specific decisions in medical applications. The major finding of this work is that it is possible to quantitatively combine the benefits and detriments of any radiation exposure situation involving an individual or population to perform cost-effectiveness analyses using the ICRP key principles of radiation protection. This FOM fills a chronic gap in the application of radiation-protection theory, i.e., limitations of generalized frameworks to algorithmically justify and optimize radiation exposures. This new framework potentially enhances objective optimization and justification, especially in complex exposure situations.

Biological optimization for mediastinal lymphoma radiotherapy - a preliminary study

Purpose: In current radiotherapy (RT) planning and delivery, population-based dose-volume constraints are used to limit the risk of toxicity from incidental irradiation of organs at risks (OARs). However, weighing tradeoffs between target coverage and doses to OARs (or prioritizing different OARs) in a quantitative way for each patient is challenging. We introduce a novel RT planning approach for patients with mediastinal Hodgkin lymphoma (HL) that aims to maximize overall outcome for each patient by optimizing on tumor control and mortality from late effects simultaneously.Material and Methods: We retrospectively analyzed 34 HL patients treated with conformal RT (3DCRT). We used published data to model recurrence and radiation-induced mortality from coronary heart disease and secondary lung and breast cancers. Patient-specific doses to the heart, lung, breast, and target were incorporated in the models as well as age, sex, and cardiac risk factors (CRFs). A preliminary plan of candidate beams was created for each patient in a commercial treatment planning system. From these candidate beams, outcome-optimized (O-OPT) plans for each patient were created with an in-house optimization code that minimized the individual risk of recurrence and mortality from late effects. O-OPT plans were compared to VMAT plans and clinical 3DCRT plans.Results: O-OPT plans generally had the lowest risk, followed by the clinical 3DCRT plans, then the VMAT plans with the highest risk with median (maximum) total risk values of 4.9 (11.1), 5.1 (17.7), and 7.6 (20.3)%, respectively (no CRFs). Compared to clinical 3DCRT plans, O-OPT planning reduced the total risk by at least 1% for 9/34 cases assuming no CRFs and 11/34 cases assuming presence of CRFs.Conclusions: We developed an individualized, outcome-optimized planning technique for HL. Some of the resulting plans were substantially different from clinical plans. The results varied depending on how risk models were defined or prioritized.

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.