Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations

Multivariable model for predicting acute oral mucositis during combined IMRT and chemotherapy for locally advanced nasopharyngeal cancer patients

INTRODUCTION/OBJECTIVE

Oral and oropharyngeal mucositis (OM) represents amultifactorialand complexinterplayof patient-, tumor-, and treatment-related factors. We aimed to build a predictive model for acute OM for locally advanced nasopharyngeal carcinoma (NPC) patients by combining clinical and dosimetric factors.

MATERIALS/METHODS

A series of consecutive NPC patients treated curatively with IMRT/VMAT + chemotherapy at 70 Gy (2-2.12 Gy/fr) was considered. For each patient, clinical- tumor- and treatment-related data were retrospectively collected. oral cavity (OC) and parotid glands (PG, considered as a single organ) were selected as organs-at-risk (OARs). Acute OM was assessed according to CTCAE v4.0 at baseline and weekly during RT. Two endpoints were considered: grade ≥3 and mean grade ≥1.5. DVHs were reduced to Equivalent Uniform Dose (EUD). Dosimetric and clinical/treatment features selected via LASSO were inserted into a multivariable logistic model. Goodness of fit was evaluated through Hosmer-Lemeshow test and calibration plot.

RESULTS

Data were collected for 132 patients. G ≥ 3 and mean G ≥ 1.5 OM were reported in 40 patients (30%). Analyses resulted in a 3-variables model for G ≥ 3 OM, including OC EUD with n = 0.05 (OR = 1.02), PG EUD with n = 1 (OR = 1.06), BMI ≥ 30 (OR = 3.8, for obese patients), and a single variable model for mean G ≥ 1.5 OM, i.e. OC EUD with n = 1 (mean dose) (OR = 1.07). Calibration was good in both cases.

CONCLUSION

OC mean dose was found to impact most on OM duration (mean G ≥ 1.5), while G ≥ 3 OM was associated to a synergic effect between PG mean dose and high dose received by small OC volumes, with BMI acting as a dose-modifying factor.

Integral dose based inverse optimization objective function promises lower toxicity in head-and-neck

PURPOSE

The voxels in a CT data sets contain density information. Besides its use in dose calculation density has no other application in modern radiotherapy treatment planning. This work introduces the use of density information by integral dose minimization in radiotherapy treatment planning for head-and-neck squamous cell carcinoma (HNSCC).

MATERIALS AND METHODS

Eighteen HNSCC cases were studied. For each case two intensity modulated radiotherapy (IMRT) plans were created: one based on dose-volume (DV) optimization, and one based on integral dose minimization (Energy hereafter) inverse optimization. The target objective functions in both optimization schemes were specified in terms of minimum, maximum, and uniform doses, while the organs at risk (OAR) objectives were specified in terms of DV- and Energy-objectives respectively. Commonly used dosimetric measures were applied to assess the performance of Energy-based optimization. In addition, generalized equivalent uniform doses (gEUDs) were evaluated. Statistical analyses were performed to estimate the performance of this novel inverse optimization paradigm.

RESULTS

Energy-based inverse optimization resulted in lower OAR doses for equivalent target doses and isodose coverage. The statistical tests showed dose reduction to the OARs with Energy-based optimization ranging from ∼2% to ∼15%.

CONCLUSIONS

Integral dose minimization based inverse optimization for HNSCC promises lower doses to nearby OARs. For comparable therapeutic effect the incorporation of density information into the optimization cost function allows reduction in the normal tissue doses and possibly in the risk and the severity of treatment related toxicities.

Image guidance and positioning accuracy in clinical practice: influence of positioning errors and imaging dose on the real dose distribution for head and neck cancer treatment

BACKGROUND

Modern radiotherapy offers the possibility of highly accurate tumor treatment. To benefit from this precision at its best, regular positioning verification is necessary. By the use of image-guided radiotherapy and the application of safety margins the influence of positioning inaccuracies can be counteracted. In this study the effect of additional imaging dose by set-up verification is compared with the effect of dose smearing by positioning inaccuracies for a collective of head-and-neck cancer patients.

METHODS

This study is based on treatment plans of 40 head-and-neck cancer patients. To evaluate the imaging dose several image guidance scenarios with different energies, techniques and frequencies were simulated and added to the original plan. The influence of the positioning inaccuracies was assessed by the use of real applied table shifts for positioning. The isocenters were shifted back appropriately to these values to simulate that no positioning correction had been performed. For the single fractions the shifted plans were summed considering three different scenarios: The summation of only shifted plans, the consideration of the original plan for the fractions with set-up verification, and the addition of the extra imaging dose to the latter. For both effects (additional imaging dose and dose smearing), plans were analyzed and compared considering target coverage, sparing of organs at risk (OAR) and normal tissue complication probability (NTCP).

RESULTS

Daily verification of the patient positioning using 3D imaging with MV energies result in non-negligible high doses. kV imaging has only marginal influence on plan quality, primarily related to sparing of organs at risk, even with daily 3D imaging. For this collective, sparing of organs at risk and NTCP are worse due to potential positioning errors.

CONCLUSION

Regular set-up verification is essential for precise radiation treatment. Relating to the additional dose, the use of kV modalities is uncritical for any frequency and technique. Dose smearing due to positioning errors for this collective mainly resulted in a decrease of OAR sparing. Target coverage also suffered from the positioning inaccuracies, especially for individual patients. Taking into account both examined effects the relevance of an extensive IGRT is clearly present, even at the expense of additional imaging dose and time expenditure.

Salvage reirradiation for local failure of prostate cancer after curative radiation therapy: Association of rectal toxicity with dose distribution and normal-tissue complication probability models

PURPOSE

This study aimed to assess the impact of radiation dose on rectal toxicity after salvage external beam radiation therapy (EBRT) with or without a brachytherapy boost for exclusive local failures after the primary EBRT for prostate cancer.

METHODS AND MATERIALS

Fourteen patients with no severe residual late toxicity after primary EBRT ± brachytherapy were reirradiated after a median time interval of 6.1 years. The median normalized total dose in 2 Gy fractions (NTD2Gy, α/β ratio = 1.5 Gy for prostate cancer cells) was 74 Gy at primary EBRT and 85.1 Gy at reirradiation. Rectal dose-volume histograms (converted to NTD2Gy_alpha/beta = 3 Gy) and the corresponding normal-tissue complication probability (NTCP) values for gastrointestinal (GI) toxicity were evaluated for 2 groups: High GI toxicity (grade ≥3) and low GI toxicity (grade ≤2).

RESULTS

The 5-year grade ≥3 GI toxicity-free survival rate was 57.1%. The median rectal V70Gy and maximum dose to 1 cm3 (D1ccrect) at primary EBRT were both predictive for grade ≥3 GI toxicity (9% vs 0%; P = .04 and 72.2 Gy vs 66.8 Gy; P < .01, respectively). When adding primary radiation therapy (RT) and reirradiation plans, the median D1ccrect was 139.8 Gy versus 126.7 Gy (P < .01) for high and low GI toxicity groups. NTCP >10% at primary RT was predictive for high GI toxicity at reirradiation (P < .05).

CONCLUSIONS

Even in the absence of residual toxicity after primary RT, rectal doses >70 Gy and NTCP >10% calculated for a first irradiation may be associated with a higher risk of developing high GI toxicity at reirradiation with a possible D1ccrect threshold of 130 Gy.

Impact of potentially variable RBE in liver proton therapy

Currently, the relative biological effectiveness (RBE) is assumed to be constant with a value of 1.1 in proton therapy. Although trends of RBE variations are well known, absolute values in patients are associated with considerable uncertainties. This study aims to evaluate the impact of a variable proton RBE in proton therapy liver trials using different fractionation schemes. Sixteen liver cancer cases were evaluated assuming two clinical schedules of 40 Gy/5 fractions and 58.05 Gy/15 fractions. The linear energy transfer (LET) and physical dose distribution in patients were simulated using Monte Carlo. The variable RBE distribution was calculated using a phenomenological model, considering the influence of the LET, fraction size and α/β value. Further, models to predict normal tissue complication probability (NTCP) and tumor control probability (TCP) were used to investigate potential RBE effects on outcome predictions. Applying the variable RBE model to the 5 and 15 fractions schedules results in an increase in mean fraction-size equivalent dose (FED) to the normal liver of 5.0% and 9.6% respectively. For patients with a mean FED to the normal liver larger than 29.8 Gy, this results in a non-negligible increase in the predicted NTCP of the normal liver averaging 11.6%, ranging from 2.7% to 25.6%. On the other hand, decrease in TCP was less than 5% for both fractionation regimens for all patients when assuming a variable RBE instead of constant. Consequently, the difference in TCP between the two fractionation schedules did not change significantly assuming a variable RBE while the impact on the NTCP difference was highly case specific. In addition, both the NTCP and TCP decrease with increasing α/β value for both fractionation schemes, with the decreases being more pronounced when using a variable RBE compared to using RBE  =  1.1. Assuming a constant RBE of 1.1 most likely overestimates the therapeutic ratio in proton therapy for liver cancer, predominantly due to underestimation of the RBE-weighted dose to the normal liver. The impact of applying a variable RBE (as compared to RBE  =  1.1) on the NTCP difference of the two fractionation regimens is case dependent. A variable RBE results in a slight increase in TCP difference. Variations in patient radiosensitivity increase when using a variable RBE.

Multivariate normal tissue complication probability modeling of vaginal late toxicity after brachytherapy for cervical cancer

PURPOSE

To explore the best variables combination for a predictive model of vaginal toxicity in cervical cancer patients undergoing brachytherapy (BT).

METHODS AND MATERIALS

Clinical and 3-dimensional dosimetric parameters were retrospectively extracted from an institutional database of consecutive patients undergoing intracavitary BT after external beam radiation therapy from 2006 to 2013 for a cervical cancer. A least absolute shrinkage and selection operator selection procedure in Cox's proportional hazards regression model was performed to select a set of relevant predictors for a multivariate normal tissue complication probability model of Grade ≥2 vaginal late toxicity. Outcomes reliability was internally assessed by bootstrap resampling method.

RESULTS

One hundred sixty-nine women were included in the present study with a median followup time of 3.8 years (interquartile range [IQR], 1.9-5.6 years). The 2 years and 5 years cumulative incidence rates of Grade ≥2 late vaginal toxicity were 19.9% and 27.5%, respectively. Among 31 metrics and six clinical factors extracted, the optimal model included two dosimetric variables: V70_Gy and D5_% (the percentage volume that received a dose greater or equal to 70 Gy and the minimum dose given to the hottest 5% volume, respectively). Area under the ROC curve at 2 and 5 years of followup were 0.85 and 0.91, respectively. Regarding internal validation, median area under the ROC curve of bootstrap predictions was 0.83 (IQR, 0.78-0.88) and 0.89 (IQR, 0.85-0.93) at 2 and 5 years of followup, respectively.

CONCLUSIONS

A multivariate normal tissue complication probability model for severe vaginal toxicity based on two dosimetric variables (V70_Gy and D5_%) provides reliable discrimination capability in a cohort of cervical cancer treated with external beam radiation therapy and BT.

Impact of dose calculation algorithms on the dosimetric and radiobiological indices for lung stereotactic body radiotherapy (SBRT) plans calculated using LQ–L model

Purpose

To investigate discrepancies in dose calculation algorithms used for lung stereotactic body radiotherapy (SBRT) plans.

Methods and materials

In total, 30 patients lung SBRT treatment plans, initially generated using BrainLab Pencil Beam (BL_PB) algorithm for 10 Gy×5 Fractions to the planning target volume (PTV) were included in the study. These plans were recalculated using BrainLab Monte Carlo (BL_MC), Eclipse AAA (EC_AAA), Eclipse Acuros XB (EC_AXB) and ADAC Pinnacle CCC (AP_CCC) algorithms. Dose volume histograms of PTV were used to calculate dosimetric and radiobiological quality indices, and equivalent dose to 2 Gy per fraction using linear-quadratic-linear model. The BL_MC algorithm is considered gold standard tool to compare PTV parameters and quality indices to investigate dose calculation discrepancies of abovementioned plans.

Results

BL_PB overestimates doses that may be due to inability of the algorithm to properly account for electron scattering and transport in inhomogeneous medium. Compared with BL_MCNO plans, the EC_AAA and EC_AXB yield lower homogeneity indices and overestimate the dose in the penumbra region, whereas AP_CCC plans were comparable for small PTV (≈8 cc) and had significant difference for large PTV.

Conclusion

BL_PB algorithm overestimates PTV doses than BL_MC calculated doses. The EC_AAA, EC_AXB and AP_CCC algorithms calculate doses within acceptable limits of radiotherapy dose delivery recommendations.

Dosimetric variations in calculation grid size in prostate VMAT: a dose-volume histogram analysis using the Gaussian error function

Background

Varying the calculation grid size can change the results of dose-volume and radiobiological parameters in a treatment plan, and therefore has an impact on the treatment planning quality assurance.

Purpose

This study investigated the dosimetric influence of the calculation grid size variation in the prostate volumetric modulated arc therapy (VMAT) plan.

Methods and materials

Dose distributions of 10 prostate VMAT plans were acquired using calculation grid sizes of 1–5 mm. Dose-volume histogram (DVH) analysis was carried out to determine the dose-volume variation corresponding to the grid size change using the Gaussian error function (GEF). At the same time, dose-volume points, dose-volume parameters and radiobiological parameters were calculated based on DVHs of targets and organs at risk (OARs) for each grid size.

Results

Comparing percentage variations of GEF parameters between the planning target volume (PTV) and clinical target volume (CTV), GEF parameters of the PTV were found varied more significantly than the CTV. This resulted in larger variations of dose-volume (%ΔCI=40·02 versus 13·55%, %ΔHI=12·45 versus 2·93% and %ΔGI=0·22 versus 0·06%) and radiobiological parameters (%ΔTCP=0·61 versus 0·25% and %ΔEUD=2·11 versus 0·26%) of the PTV compared with CTV. For OARs, the rectal wall showed a larger dose-volume variation than the rectum. However, similar dose-volume variation due to grid size change was not found in the bladder, bladder wall and femur.

Conclusions

Knowing the dosimetric variation in this study is important to the radiotherapy staff in the quality assurance for the prostate VMAT planning.

Radiobiological Models of Tissue Response to Radiation in Treatment Planning Systems

Aims

To present several biological concepts and models of tissue response to fractionated radiotherapy. To describe practical implementation of these models in three-dimensional treatment planning systems.

Methods

Models of cell survival, Equivalent Uniform Dose (EUD) and Tumor Control Probability (TCP) are discussed. These models are based on the target-cell hypothesis which assumes that response of organs and tissues to radiation therapy can be explained and mathematically described in terms of survival of the specific target-cells.

Results

Several formulae for deriving and calculating EUD and TCP for a given three-dimensional dose distribution are presented and discussed.

Conclusions

Biological models of tissue response to radiation, when used wisely, have a potential to be useful in radiation therapy treatment planning. The models can advance our understanding of the underlying biological mechanisms, and may help in designing new and better treatment strategies. They should be particularly useful in modern conformai radiotherapy where treatment strategy for each patient can be individualized and optimized according to patient characteristics and available technology of delivering sophisticated treatment plans.

Dosimetric, Radiobiological and Secondary Cancer Risk Evaluation in Head-and-Neck Three-dimensional Conformal Radiation Therapy, Intensity-Modulated Radiation Therapy, and Volumetric Modulated Arc Therapy: A Phantom Study

This analysis estimated secondary cancer risks after volumetric modulated arc therapy (VMAT) and compared those risks to the risks associated with other modalities of head-and-neck (H&N) radiotherapy. Images of H&N anthropomorphic phantom were acquired with a computed tomography scanner and exported via digital imaging and communications in medicine (DICOM) standards to a treatment planning system. Treatment plans were performed using a VMAT dual-arc technique, a nine-field intensity-modulated radiation therapy (IMRT) technique, and a four-field three-dimensional conformal therapy (3DCRT) technique. The prescription dose was 66.0 Gy for all three techniques, but to accommodate the range of dosimeter responses, we delivered a single dose of 6.60 Gy to the isocenter. The lifetime risk for secondary cancers was estimated according to National Council on Radiation Protection and Measurements (NCRP) Report 116. VMAT delivered the lowest maximum doses to esophagus (23 Gy), and normal brain (40 Gy). In comparison, maximum doses for 3DCRT were 74% and 40%, higher than those for VMAT for the esophagus, and normal brain, respectively. The normal tissue complication probability and equivalent uniform dose for the brain (2.1%, 0.9%, 0.8% and 3.8 Gy, 2.6 Gy, 2.3 Gy) and esophagus (4.2%, 0.7%, 0.4% and 3.7 Gy, 2.2 Gy, 1.8 Gy) were calculated for the 3DCRT, IMRT and VMAT respectively. Fractional esophagus OAR volumes receiving more than 20 Gy were 3.6% for VMAT, 23.6% for IMRT, and 100% for 3DCRT. The calculations for mean doses, NTCP, EUD and OAR volumes suggest that the risk of secondary cancer induction after VMAT is lower than after IMRT and 3DCRT.

Inter-institutional analysis demonstrates the importance of lower than previously anticipated dose regions to prevent late rectal bleeding following prostate radiotherapy

PURPOSE

To investigate whether inter-institutional cohort analysis uncovers more reliable dose-response relationships exemplified for late rectal bleeding (LRB) following prostate radiotherapy.

MATERIAL AND METHODS

Data from five institutions were used. Rectal dose-volume histograms (DVHs) for 989 patients treated with 3DCRT or IMRT to 70-86.4 Gy@1.8-2.0 Gy/fraction were obtained, and corrected for fractionation effects (α/β = 3 Gy). Cohorts with best-fit Lyman-Kutcher-Burman volume-effect parameter a were pooled after calibration adjustments of the available LRB definitions. In the pooled cohort, dose-response modeling (incorporating rectal dose and geometry, and patient characteristics) was conducted on a training cohort (70%) followed by final testing on the remaining 30%. Multivariate logistic regression was performed to build models with bootstrap stability.

RESULTS

Two cohorts with low bleeding rates (2%) were judged to be inconsistent with the remaining data, and were excluded. In the remaining pooled cohorts (n = 690; LRB rate = 12%), an optimal model was generated for 3DCRT using the minimum rectal dose and the absolute rectal volume receiving less than 55 Gy (AUC = 0.67; p = 0.0002; Hosmer-Lemeshow p-value, pHL = 0.59). The model performed nearly as well in the hold-out testing data (AUC = 0.71; p < 0.0001; pHL = 0.63), indicating a logistically shaped dose-response.

CONCLUSION

We have demonstrated the importance of integrating datasets from multiple institutions, thereby reducing the impact of intra-institutional dose-volume parameters explicitly correlated with prescription dose levels. This uncovered an unexpected emphasis on sparing of the low to intermediate rectal dose range in the etiology of late rectal bleeding following prostate radiotherapy.

Assessing the uncertainty in a normal tissue complication probability difference (∆NTCP): radiation-induced liver disease (RILD) in liver tumour patients treated with proton vs X-ray therapy

Modern radiotherapy technologies such as proton beam therapy (PBT) permit dose escalation to the tumour and minimize unnecessary doses to normal tissues. To achieve appropriate patient selection for PBT, a normal tissue complication probability (NTCP) model can be applied to estimate the risk of treatment-related toxicity relative to X-ray therapy (XRT). A methodology for estimating the difference in NTCP (∆NTCP), including its uncertainty as a function of dose to normal tissue, is described in this study using the Delta method, a statistical method for evaluating the variance of functions, considering the variance-covariance matrix. We used a virtual individual patient dataset of radiation-induced liver disease (RILD) in liver tumour patients who were treated with XRT as a study model. As an alternative option for individual patient data, dose-bin data, which consists of the number of patients who developed toxicity in each dose level/bin and the total number of patients in that dose level/bin, are useful for multi-institutional data sharing. It provides comparable accuracy with individual patient data when using the Delta method. With reliable NTCP models, the ∆NTCP with uncertainty might potentially guide the use of PBT; however, clinical validation and a cost-effectiveness study are needed to determine the appropriate ∆NTCP threshold.

The Technique of Intensity-Modulated Radiotherapy in the Treatment of Cholangiocarcinoma

Aims and background

Conventional radiotherapy in inoperable cholangiocarcinoma is limited by radiotolerance of the surrounding tissues. The aim of our dosimetric study was an evaluation of intensity-modulated radiotherapy in the treatment of inoperable bile duct carcinoma.

Methods

Four patients with inoperable cholangiocarcinoma treated by self-expandable stent placed to the biliary tree and radiotherapy were studied. The rotational technique, conformal 3D BOX technique and intensity-modulated radiotherapy plan were compared. Dose volume histograms and the normal tissue complication probability concept were used for comparison. The stent was used for target motion verification.

Results

The intensity-modulated radiotherapy plans showed favorable dose distribution in planning target volume and remarkable sparing of organs at risk.

Conclusions

The intensity-modulated radiotherapy technique in bile duct carcinomas deserves further research and clinical evaluation.

Dose-volume and radiobiological dependence on the calculation grid size in prostate VMAT planning

This study evaluated the effects of dose-volume and radiobiological dependence on the calculation grid size in prostate volumetric-modulated arc therapy (VMAT) planning. Ten patients with prostate cancer were selected for this retrospective treatment planning study. Prostate VMAT plans were created for the patients using the 6 MV photon beam produced by a Varian TrueBEAM linac with the calculation grid size equal to 1, 2, 2.5, 3, 4, and 5 mm. Dose-volume histograms (DVHs) of targets and organs at risk were generated for different grid sizes. We calculated the radiobiological parameters of the tumor control probability (TCP) of clinical target volume (CTV) and planning target volume (PTV), and the normal tissue complication probability (NTCP) of organs at risk (rectal wall, rectum, bladder wall, bladder, left femur, and right femur). The homogeneity, conformity, and gradient indexes of CTV and PTV were calculated for different grid sizes. The TCP of PTV was found decreasing with a rate of 0.06%/mm when the calculation grid size increased from 1 to 5 mm. On the other hand, both NTCPs of rectal wall and rectum were found decreasing with rates of 0.03%/mm and 0.05%/mm, respectively, with an increase of grid size. The homogeneity index of PTV increased with a rate of 0.57/mm of the calculation grid size, whereas the conformity index of PTV decreased with a rate of 0.0075/mm. The gradient index of PTV was found increasing with a rate equal to 0.05/mm. In prostate VMAT planning, variations of dose-volume and radiobiological parameters with calculation grid size on PTV, rectal wall, and rectum were more significant than those of CTV and other organs at risk such as bladder wall, bladder, and femurs. Results in this study are important in the treatment planning quality assurance when the calculation grid size is varied to compromise a shorter dose computing time.

Radiobiological analysis of stereotactic body radiation therapy for an evidence-based planning target volume of the lung using multiphase CT images obtained with a pneumatic abdominal compression apparatus: a case study

This study evaluated the efficiency of stereotactic body radiation therapy of lung (SBRT-Lung) in generating a treatment volume using conventional multiple-phase three-dimensional computed tomography (3D-CT) of a patient immobilized with pneumatic abdominal compression. The institutional protocol for SBRT-Lung using the RapidArc technique relied on a planning target volume (PTV) delineated using 3D-CT and accounted for linear and angular displacement of the tumor during respiratory movements. The efficiency of the institutional protocol was compared with that of a conventional method for PTV delineation based on radiobiological estimates, such as tumor control probability (TCP) and normal tissue complication probability (NTCP), evaluated using dose-volume parameters. Pneumatic abdominal compression improved the TCP by 15%. This novel protocol improved the TCP by 0.5% but reduced the NTCP for lung pneumonitis (0.2%) and rib fracture (1.0%). Beyond the observed variations in the patient's treatment setup, the institutional protocol yielded a significantly consistent TCP (p < 0.005). The successful clinical outcome of this case study corroborates predictions based on radiobiological evaluation and deserves validation through an increased number of patients.

New Clinical Features and Dosimetric Predictor Identification for Symptomatic Radiation Pneumonitis after Tangential Irradiation in Breast Cancer Patients

Background: Tangential irradiation is the most popular postoperative radiotherapy technique for breast cancer. However, irradiation has been related to symptomatic radiation pneumonitis (SRP), which decreases the quality of life of patients. This study investigated the clinical features and dosimetric parameters related to SRP of the ipsilateral lung to identify risk factors for SRP in breast cancer patients after three-dimensional conformal radiation therapy (3D-CRT) with tangential fields. Material and Methods: A total of 515 breast cancer patients were evaluated and divided into two groups: the local-regional irradiation group (259 patients) and the simple local irradiation group (256 patients). Clinical symptoms were registered and patient data collected. The relationship between the incidence of SRP and dosimetric parameters for the ipsilateral lung was assessed within 6 months after 3D-CRT. Dosimetric parameters were compared using t tests. The dosimetric predictors for SRP were estimated using a logistic regression model and receiver operating characteristic curve analysis. Results: In total, 19 patients (3.7%) developed grade 2 SRP. In the local-regional irradiation group, the probability of SRP in the lung body was greater than that in the lung apex (3.9% vs. 1.5%). V20 and V30 were independent predictors for SRP in the local-regional irradiation group (odds ratio = 1.152 and 1.439, both p = 0.030), whereas only V20 was an independent predictor of SRP in the simple local irradiation group (odds ratio = 1.351, p = 0.001). With 39.8% as the optimal threshold for V20 and 25.7% for V30 for local-regional irradiation, SRP could be predicted with an accuracy of 80.3% and 79.9%, a sensitivity of 61.5% and 69.2%, and a specificity of 81.3% and 80.5%, respectively. With 20.2% as the optimal V20 threshold for simple local irradiation, SRP could be predicted with an accuracy of 88.7%, a sensitivity of 83.3% and a specificity of 89.6%. Conclusions: SRP has become a rare complication with mild symptoms and occurs mainly in the lung body. V20 and V30 may be useful dosimetric predictors to evaluate SRP risk of the ipsilateral lung in breast cancer.

Modelling duodenum radiotherapy toxicity using cohort dose-volume-histogram data

Background and purpose

Gastro-intestinal toxicity is dose-limiting in abdominal radiotherapy and correlated with duodenum dose-volume parameters. We aimed to derive updated NTCP model parameters using published data and prospective radiotherapy quality-assured cohort data.

Material and methods

A systematic search identified publications providing duodenum dose-volume histogram (DVH) statistics for clinical studies of conventionally-fractionated radiotherapy. Values for the Lyman-Kutcher-Burman (LKB) NTCP model were derived through sum-squared-error minimisation and using leave-one-out cross-validation. Data were corrected for fraction size and weighted according to patient numbers, and the model refined using individual patient DVH data for two further cohorts from prospective clinical trials.

Results

Six studies with published DVH data were utilised, and with individual patient data included outcomes for 531 patients in total (median follow-up 16 months). Observed gastro-intestinal toxicity rates ranged from 0% to 14% (median 8%). LKB parameter values for unconstrained fit to published data were: n = 0.070, m = 0.46, TD₅₀₍₁₎ [Gy] = 183.8, while the values for the model incorporating the individual patient data were n = 0.193, m = 0.51, TD₅₀₍₁₎ [Gy] = 299.1.

Conclusions

LKB parameters derived using published data are shown to be consistent to those previously obtained using individual patient data, supporting a small volume-effect and dependence on exposure to high threshold dose.

Predictors of Liver Toxicity Following Stereotactic Body Radiation Therapy for Hepatocellular Carcinoma

PURPOSE

To identify risk factors associated with a decline in liver function after stereotactic body radiation therapy (SBRT) for hepatocellular carcinoma.

METHODS AND MATERIALS

Data were analyzed from patients with hepatocellular carcinoma treated on clinical trials of 6-fraction SBRT. Liver toxicity was defined as an increase in Child-Pugh (CP) score ≥2 three months after SBRT. Clinical factors, SBRT details, and liver dose-volume histogram (DVH) parameters were tested for association with toxicity using logistic regression. CP class B patients were analyzed separately.

RESULTS

Among CP class A patients, 101 were evaluable, with a baseline score of A5 (72%) or A6 (28%). Fifty-three percent had portal vein thrombus. The median liver volume was 1286 cc (range, 766-3967 cc), and the median prescribed dose was 36 Gy (range, 27-54 Gy). Toxicity was seen in 26 patients (26%). Thrombus, baseline CP of A6, and lower platelet count were associated with toxicity on univariate analysis, as were several liver DVH-based parameters. Absolute and spared liver volumes were not significant. On multivariate analysis for CP class A patients, significant associations were found for baseline CP score of A6 (odds ratio [OR], 4.85), lower platelet count (OR, 0.90; median, 108 × 10⁹/L vs 150 × 10⁹/L), higher mean liver dose (OR, 1.33; median, 16.9 Gy vs 14.7 Gy), and higher dose to 800 cc of liver (OR, 1.11; median, 14.3 Gy vs 6.0 Gy). With 13 CP-B7 patients included or when dose to 800 cc of liver was replaced with other DVH parameters (eg, dose to 700 or 900 cc of liver) in the multivariate analysis, effective volume and portal vein thrombus were associated with an increased risk.

CONCLUSIONS

Baseline CP scores and higher liver doses (eg, mean dose, effective volume, doses to 700-900 cc) were strongly associated with liver function decline 3 months after SBRT. A lower baseline platelet count and portal vein thrombus were also associated with an increased risk.

Rectal dose to prostate cancer patients treated with proton therapy with or without rectal spacer

The purpose of this study was to evaluate whether a spacer inserted in the prerectal space could reduce modeled rectal dose and toxicity rates for patients with prostate cancer treated in silico with pencil beam scanning (PBS) proton therapy. A total of 20 patients were included in this study who received photon therapy (12 with rectal spacer (DuraSeal™ gel) and 8 without). Two PBS treatment plans were retrospectively created for each patient using the following beam arrangements: (1) lateral-opposed (LAT) fields and (2) left and right anterior oblique (LAO/RAO) fields. Dose volume histograms (DVH) were generated for the prostate, rectum, bladder, and right and left femoral heads. The normal tissue complication probability (NTCP) for ≥grade 2 rectal toxicity was calculated using the Lyman-Kutcher-Burman model and compared between patients with and without the rectal spacer. A significantly lower mean rectal DVH was achieved in patients with rectal spacer compared to those without. For LAT plans, the mean rectal V70 with and without rectal spacer was 4.19 and 13.5%, respectively. For LAO/RAO plans, the mean rectal V70 with and without rectal spacer was 5.07 and 13.5%, respectively. No significant differences were found in any rectal dosimetric parameters between the LAT and the LAO/RAO plans generated with the rectal spacers. We found that ≥ 9 mm space resulted in a significant decrease in NTCP modeled for ≥grade 2 rectal toxicity. Rectal spacers can significantly decrease modeled rectal dose and predicted ≥grade 2 rectal toxicity in prostate cancer patients treated in silico with PBS. A minimum of 9 mm separation between the prostate and anterior rectal wall yields the largest benefit.

Dosimetric and radiobiological comparison of Cyberknife and Tomotherapy in stereotactic body radiotherapy for localized prostate cancer

BACKGROUND AND PURPOSE

As recent studies have suggested relatively low α/β for prostate cancer, the interest in hypofractionated stereotactic body radiotherapy (SBRT) for prostate cancer is rising. The aim of this study is to compare dosimetric results of Cyberknife (CK) with Tomotherapy (HT) in SBRT for localized prostate cancer. Furthermore, the radiobiologic consequences of heterogeneous dose distribution are also analyzed.

MATERIAL AND METHOD

A total of 12 cases of localized prostate cancer previously treated with SBRT were collected. Treatments had been planned and delivered using CK. Then HT plans were generated for comparison afterwards. The prescribed dose was 37.5Gy in 5 fractions. Dosimetric indices for target volumes and organs at risk (OAR) were compared. For radiobiological evaluation, generalized equivalent uniform dose (gEUD) and normal tissue complication probability (NTCP) were calculated and compared.

RESULT

Both CK and HT achieved target coverage while meeting OAR constraints adequately. HT plans resulted in better dose homogeneity (Homogeneity index: 1.04±0.01 vs. 1.21±0.01; p = 0.0022), target coverage (97.74±0.86% vs. 96.56±1.17%; p = 0.0076) and conformity (new vonformity index: 1.16±0.05 vs. 1.21±0.04; p = 0.0096). HT was shown to predict lower late rectal toxicity as compared to CK. Integral dose to body was also significantly lower in HT plans (46.59±6.44 Gy'L vs 57.05±11.68 Gy'L; p = 0.0029).

CONCLUSION

Based on physical dosimetry and radiobiologic considerations, HT may have advantages over CK, specifically in rectal sparing which could translate into clinical benefit of decreased late toxicities.

Impact of dose calculation models on radiotherapy outcomes and quality adjusted life years for lung cancer treatment: do we need to measure radiotherapy outcomes to tune the radiobiological parameters of a normal tissue complication probability model?

BACKGROUND

The equivalent uniform dose (EUD) radiobiological model can be applied for lung cancer treatment plans to estimate the tumor control probability (TCP) and the normal tissue complication probability (NTCP) using different dose calculation models. Then, based on the different calculated doses, the quality adjusted life years (QALY) score can be assessed versus the uncomplicated tumor control probability (UTCP) concept in order to predict the overall outcome of the different treatment plans.

METHODS

Nine lung cancer cases were included in this study. For the each patient, two treatments plans were generated. The doses were calculated respectively from pencil beam model, as pencil beam convolution (PBC) turning on 1D density correction with Modified Batho's (MB) method, and point kernel model as anisotropic analytical algorithm (AAA) using exactly the same prescribed dose, normalized to 100% at isocentre point inside the target and beam arrangements. The radiotherapy outcomes and QALY were compared. The bootstrap method was used to improve the 95% confidence intervals (95% CI) estimation. Wilcoxon paired test was used to calculate P value.

RESULTS

Compared to AAA considered as more realistic, the PBC_MB overestimated the TCP while underestimating NTCP, P<0.05. Thus the UTCP and the QALY score were also overestimated.

CONCLUSIONS

To correlate measured QALY's obtained from the follow-up of the patients with calculated QALY from DVH metrics, the more accurate dose calculation models should be first integrated in clinical use. Second, clinically measured outcomes are necessary to tune the parameters of the NTCP model used to link the treatment outcome with the QALY. Only after these two steps, the comparison and the ranking of different radiotherapy plans would be possible, avoiding over/under estimation of QALY and any other clinic-biological estimates.

Modeling Radiotherapy Induced Normal Tissue Complications: An Overview beyond Phenomenological Models

An overview of radiotherapy (RT) induced normal tissue complication probability (NTCP) models is presented. NTCP models based on empirical and mechanistic approaches that describe a specific radiation induced late effect proposed over time for conventional RT are reviewed with particular emphasis on their basic assumptions and related mathematical translation and their weak and strong points.

Investigation of whether in-room CT-based adaptive intracavitary brachytherapy for uterine cervical cancer is robust against interfractional location variations of organs and/or applicators

This study investigates whether in-room computed tomography (CT)-based adaptive treatment planning (ATP) is robust against interfractional location variations, namely, interfractional organ motions and/or applicator displacements, in 3D intracavitary brachytherapy (ICBT) for uterine cervical cancer. In ATP, the radiation treatment plans, which have been designed based on planning CT images (and/or MR images) acquired just before the treatments, are adaptively applied for each fraction, taking into account the interfractional location variations. 2D and 3D plans with ATP for 14 patients were simulated for 56 fractions at a prescribed dose of 600 cGy per fraction. The standard deviations (SDs) of location displacements (interfractional location variations) of the target and organs at risk (OARs) with 3D ATP were significantly smaller than those with 2D ATP (P < 0.05). The homogeneity index (HI), conformity index (CI) and tumor control probability (TCP) in 3D ATP were significantly higher for high-risk clinical target volumes than those in 2D ATP. The SDs of the HI, CI, TCP, bladder and rectum D_2cc, and the bladder and rectum normal tissue complication probability (NTCP) in 3D ATP were significantly smaller than those in 2D ATP. The results of this study suggest that the interfractional location variations give smaller impacts on the planning evaluation indices in 3D ATP than in 2D ATP. Therefore, the 3D plans with ATP are expected to be robust against interfractional location variations in each treatment fraction.

A treatment planning study of proton arc therapy for para-aortic lymph node tumors: dosimetric evaluation of conventional proton therapy, proton arc therapy, and intensity modulated radiotherapy

BACKGROUND

The purpose of this study is to evaluate the dosimetric benefits of a proton arc technique for treating tumors of the para-aortic lymph nodes (PALN).

METHOD

In nine patients, a proton arc therapy (PAT) technique was compared with intensity modulated radiation therapy (IMRT) and proton beam therapy (PBT) techniques with respect to the planning target volume (PTV) and organs at risk (OAR). PTV coverage, conformity index (CI), homogeneity index (HI) and OAR doses were compared. Organ-specific radiation induced cancer risks were estimated by applying organ equivalent dose (OED) and normal tissue complication probability (NTCP).

RESULTS

The PAT techniques showed better PTV coverage than IMRT and PBT plans. The CI obtained with PAT was 1.19 ± 0.02, which was significantly better than that for the IMRT techniques. The HI was lowest for the PAT plan and highest for IMRT. The dose to the OARs was always below the acceptable limits and comparable for all three techniques. OED results calculated based on a plateau dose-response model showed that the risk of secondary cancers in organs was much higher when IMRT or PBT were employed than when PAT was used. NTCPs of PAT to the stomach (0.29 %), small bowel (0.69 %) and liver (0.38 %) were substantially lower than those of IMRT and PBT.

CONCLUSION

This study demonstrates that there is a potential role for PAT as a commercialized instrument in the future to proton therapy.

New approach in lung cancer radiotherapy offers better normal tissue sparing

PURPOSE

Medical images are more than pictures. They contain additional quantitative information which can be interrogated, quantified, and utilized. Besides anatomical information computed tomography (CT) imaging data provide electron density information. Radiotherapy use of this density information is limited to its application only in dose calculations. The direct product of dose, density, and volume forms a quantity called integral dose. The integral dose delivered to a volume of interest is the total energy deposited in that volume. Here it is hypothesized that minimization of the integral dose is advantageous in radiotherapy planning. The purpose of this work is to study the incorporation of quantitative imaging information in radiotherapy inverse optimization through total energy minimization (Energy hereafter).

DESIGN

Twenty lung patient plans were studied. For each patient density was quantified on voxel-by-voxel basis through image gray value-to-density conversion curves. Energy-based objective function was used for inverse radiotherapy plan optimization. The obtained plans were evaluated in the light of current standard of care, based on dose-volume (DVH) optimization approach.

RESULTS

The statistical significance analyses of the results indicated that the doses to normal tissue were between 14% and 45% lower, when Energy-based optimization was used instead of DVH-based optimization.

CONCLUSION

Incorporation of quantitative imaging information, through CT derived density, in the optimization cost function allows reduction of dose to normal tissue for NSCLC cases. Energy-based radiotherapy plans result in lower normal tissue dose and potentially lower complication rates compared to standard of care.

Dosimetric and radiobiological characterizations of prostate intensity-modulated radiotherapy and volumetric-modulated arc therapy: A single-institution review of ninety cases

This study reviewed prostate volumetric-modulated arc therapy (VMAT) plans with intensity-modulated radiotherapy (IMRT) plans after prostate IMRT technique was replaced by VMAT in an institution. Characterizations of dosimetry and radiobiological variation in prostate were determined based on treatment plans of 40 prostate IMRT patients (planning target volume = 77.8–335 cm³) and 50 VMAT patients (planning target volume = 120–351 cm³) treated before and after 2013, respectively. Both IMRT and VMAT plans used the same dose-volume criteria in the inverse planning optimization. Dose-volume histogram, mean doses of target and normal tissues (rectum, bladder and femoral heads), dose-volume points (D_99% of planning target volume; D_30%, D_50%, V_{30 Gy} and V_{35 Gy} of rectum and bladder; D_5%, V_{14 Gy}, V_{22 Gy} of femoral heads), conformity index (CI), homogeneity index (HI), gradient index (GI), prostate tumor control probability (TCP), and rectal normal tissue complication probability (NTCP) based on the Lyman-Burman-Kutcher algorithm were calculated for each IMRT and VMAT plan. From our results, VMAT plan was found better due to its higher (1.05%) CI, lower (0.83%) HI and (0.75%) GI than IMRT. Comparing doses in normal tissues between IMRT and VMAT, it was found that IMRT mostly delivered higher doses of about 1.05% to the normal tissues than VMAT. Prostate TCP and rectal NTCP were found increased (1%) for VMAT than IMRT. It is seen that VMAT technique can decrease the dose-volume evaluation criteria for the normal tissues. Based on our dosimetric and radiobiological results in treatment plans, it is concluded that our VMAT implementation could produce comparable or slightly better target coverage and normal tissue sparing with a faster treatment time in prostate radiotherapy.

Assessing the shift of radiobiological metrics in lung radiotherapy plans using 2D gamma index

BACKGROUND

The purpose of this work is to investigate the 2D gamma (γ) maps to illustrate the change of radiobiological outcomes for lung radiotherapy plans and evaluate the correlation between tumor control probability (TCP), normal tissue complication probability (NTCP) with γ passing rates (γ-rates).

METHODS

Nine patients with lung cancer were used. The doses were calculated using Modified Batho method integrated with pencil beam convolution (MB-PBC) and anisotropic analytical algorithm (AAA) using the same beam arrangements and prescription dose. The TCP and NTCP were estimated, respectively, using equivalent uniform dose (EUD) model and Lyman-Kutcher-Burman (LKB) model. The correlation between ΔTCP or ΔNTCP with γ-rates, from 2%/2 and 3%/3 mm, were tested to explore the best correlation predicting the relevant γ criteria using Spearman's rank test (ρ). Wilcoxon paired test was used to calculate P value.

RESULTS

TCP value was significantly lower in the recalculated AAA plans as compared to MB plans. However, AAA predicted more NTCP on lung pneumonitis according to the LKB model and using relevant radiobiological parameters (n, m and TD50) for MB-PBC and AAA, with P=0.03. The data showed a weak correlation between radiobiological metrics with γ-rates or γ-mean, ρ<0.3.

CONCLUSIONS

AAA and MB yield different TCP values as well as NTCP for lung pneumonitis based on the LKB model parameters. Therefore, 2D γ-maps, generated with 2%/2 or 3%/3 mm, could illustrate visual information about the radiobiological changes. The information is useful to evaluate the clinical outcome of a radiotherapy treatment and to approve the treatment plan of the patient if the dose constraints are respected. On the other hand, the γ-maps tool can be used as quality assurance (QA) process to check the predicted TCP and NTCP from radiobiological models.

Patient Specific Characteristics Are an Important Factor That Determines the Risk of Acute Grade ≥ 2 Rectal Toxicity in Patients Treated for Prostate Cancer with IMRT and Daily Image Guidance Based on Implanted Gold Markers

AIM

To model acute rectal toxicity in Intensity Modulated Radiation Therapy (IMRT) for prostate cancer using dosimetry and patient specific characteristics.

METHODS

A database of 79 prostate cancer patients treated with image guided IMRT was used to fit parameters of Lyman-Kutcher-Burman (LKB) and logistic regression Normal Tissue Complications Probability (NTCP) models to acute grade ≥ 2 rectal toxicities. We used a univariate regression model to find the dosimetric index which was most correlated with toxicity and a multivariate logistic regression model with machine learning algorithm to integrate dosimetry with patient specific characteristics. We used Receiver Operating Characteristics (ROC) analysis and the area under the ROC curve (AUC) to quantify the predictive power of models.

RESULTS

Sixteen patients (20.3%) developed acute grade≥2 rectal toxicity. Our best estimate (95% confidence interval) of LKB model parameters for acute rectal toxicity are exponent n=0.13 (0.1-0.16), slope m=0.09 (0.08-0.11), and threshold dose TD50=56.8 (53.7-59.9) Gy. The best dosimetric indices in the univariate logistic regression NTCP model were D25% and V50Gy. The best AUC of dosimetry only modeling was 0.67 (0.54, 0.8). In the multivariate logistic regression two patient specific variables were particularly strongly correlated with acute rectal toxicity, the use of statin drugs and PSA level prior to IMRT, while two additional variables, age and diabetes were weakly correlated. The AUC of the logistic regression NTCP model improved to 0.88 (0.8, 0.96) when patient specific characteristics were included. In a group of 79 patients, 40 took Statins and 39 did not. Among patients who took statins, (4/40)=10% developed acute grade ≥2 rectal toxicity, compared to (12/39)=30.8% who did not take statins (p=0.03). The average and standard deviation of PSA distribution for patients with acute rectal toxicity was PSA_tox = 5.77 ± 2.27 and it was PSA_notox = 9.5 ± 7.8 for the remainder (p=0.01).

CONCLUSIONS

Patient specific characteristics strongly influence the likelihood of acute grade ≥ 2 rectal toxicity in radiation therapy for prostate cancer.

Anteriorly Oriented Beam Arrangements with Daily in Vivo Range Verification for Proton Therapy of Prostate Cancer: Rectal Toxicity Rates

Purpose

To model whether in vivo range verification could reduce high-grade rectal toxicity for patients with prostate cancer treated with pencil beam scanning proton therapy by allowing novel proton beam arrangements compared to standard lateral opposed beams.

Materials and Methods

Proton plans were generated for 8 patients with prostate cancer previously treated with photons by volumetric-modulated arc therapy (VMAT). The VMAT plans were generated by using a uniform 6-mm planning target volume (PTV) expansion. For the proton plans an additional distal margin (3.5% of beam range) was added to the uniform 6-mm PTV to account for range uncertainty, using 3 beam arrangements: (1) lateral opposed beams (LAT), (2) left and right anterior oblique beams (LAO/RAO), and (3) a single anterior-posterior beam (AP). Assuming use of in vivo range verification, plans were generated by using a reduced distal PTV and distal range uncertainty expansion (2 mm each) with AP (AP-2 mm) and LAO/RAO (LAO/RAO-2 mm) beam arrangements. Estimates of normal tissue complication probability (NTCP) for ≥grade 2 rectal bleeding were generated by using the Lyman-Kutcher-Burman model.

Results

Each proton and photon plan was able to achieve all prespecified rectal and bladder constraints. For the VMAT, LAT, AP, and LAO/RAO plans, estimated NTCP values for ≥grade 2 rectal bleeding were 0.19, 0.21, 0.24, and 0.2, respectively. For the AP-2 mm and LAO/RAO-2 mm plans, NTCP values were reduced to 0.11 and 0.1 with respect to ≥grade 2 rectal bleeding.

Conclusion

Presuming that in vivo range verification for pencil beam scanning proton therapy could localize the distal falloff of the Bragg peak to within 2 mm, novel beam arrangements (AP and LAO/RAO) may reduce the risk of serious rectal bleeding, compared to VMAT and LAT proton treatment techniques. These are achieved without an increase in modeled bladder complication rates.

Assessment of radiobiological metrics applied to patient-specific QA process of VMAT prostate treatments

VMAT is a powerful technique to deliver hypofractionated prostate treatments. The lack of correlations between usual 2D pretreatment QA results and the clinical impact of possible mistakes has allowed the development of 3D verification systems. Dose determination on patient anatomy has provided clinical predictive capability to patient-specific QA process. Dose-volume metrics, as evaluation criteria, should be replaced or complemented by radiobiological indices. These metrics can be incorporated into individualized QA extracting the information for response parameters (gEUD, TCP, NTCP) from DVHs. The aim of this study is to assess the role of two 3D verification systems dealing with radiobiological metrics applied to a prostate VMAT QA program. Radiobiological calculations were performed for AAPM TG-166 test cases. Maximum differences were 9.3% for gEUD, -1.3% for TCP, and 5.3% for NTCP calculations. Gamma tests and DVH-based comparisons were carried out for both systems in order to assess their performance in 3D dose determination for prostate treatments (high-, intermediate-, and low-risk, as well as prostate bed patients). Mean gamma passing rates for all structures were bet-ter than 92.0% and 99.1% for both 2%/2 mm and 3%/3 mm criteria. Maximum discrepancies were (2.4% ± 0.8%) and (6.2% ± 1.3%) for targets and normal tis-sues, respectively. Values for gEUD, TCP, and NTCP were extracted from TPS and compared to the results obtained with the two systems. Three models were used for TCP calculations (Poisson, sigmoidal, and Niemierko) and two models for NTCP determinations (LKB and Niemierko). The maximum mean difference for gEUD calculations was (4.7% ± 1.3%); for TCP, the maximum discrepancy was (-2.4% ± 1.1%); and NTCP comparisons led to a maximum deviation of (1.5% ± 0.5%). The potential usefulness of biological metrics in patient-specific QA has been explored. Both systems have been successfully assessed as potential tools for evaluating the clinical outcome of a radiotherapy treatment in the scope of pretreatment QA.

Urinary bladder dose-response relationships for patient-reported genitourinary morbidity domains following prostate cancer radiotherapy

BACKGROUND AND PURPOSE

Radiotherapy (RT) induced genitourinary (GU) morbidity is typically assessed by physicians as single symptoms or aggregated scores including symptoms from various domains. Here we apply a method to group patient-reported GU symptoms after RT for localized prostate cancer based on their interplay, and study how these relate to urinary bladder dose.

MATERIALS AND METHODS

Data were taken from two Scandinavian studies (N=207/276) including men treated with external-beam RT (EBRT) to 78/70Gy (2Gy/fraction; median time-to-follow-up: 3.6-6.4y). Within and across cohorts, bladder dose-volume parameters were tested as predictors for GU symptom domains identified from two study-specific questionnaires (35 questions on frequency, incontinence, obstruction, pain, urgency, and sensory symptoms) using univariate and multivariate logistic regression analysis (MVA) with 10-fold cross-validation. Performance was evaluated using Area Under the Receiver Operating Characteristic Curve (Az).

RESULTS

For the identified Incontinence (2-5 symptoms), Obstruction (3-5 symptoms), and Urgency (2-7 symptoms) domains, MVA demonstrated that bladder doses close to the prescription doses were the strongest predictors for Obstruction (Az: 0.53-0.57) and Urgency (Az: 0.60). For Obstruction, performance increased for the across cohort analysis (Az: 0.61-0.64).

CONCLUSIONS

Our identified patient-reported GU symptom domains suggest that high urinary bladder doses, and increased focus on both obstruction and urgency is likely to further add to the understanding of GU tract RT responses.

Radiobiological impact of dose calculation algorithms on biologically optimized IMRT lung stereotactic body radiation therapy plans

Background

The aim of this study is to evaluate the radiobiological impact of Acuros XB (AXB) vs. Anisotropic Analytic Algorithm (AAA) dose calculation algorithms in combined dose-volume and biological optimized IMRT plans of SBRT treatments for non-small-cell lung cancer (NSCLC) patients.

Methods

Twenty eight patients with NSCLC previously treated SBRT were re-planned using Varian Eclipse (V11) with combined dose-volume and biological optimization IMRT sliding window technique. The total dose prescribed to the PTV was 60 Gy with 12 Gy per fraction. The plans were initially optimized using AAA algorithm, and then were recomputed using AXB using the same MUs and MLC files to compare with the dose distribution of the original plans and assess the radiobiological as well as dosimetric impact of the two different dose algorithms. The Poisson Linear-Quadatric (PLQ) and Lyman-Kutcher-Burman (LKB) models were used for estimating the tumor control probability (TCP) and normal tissue complication probability (NTCP), respectively. The influence of the model parameter uncertainties on the TCP differences and the NTCP differences between AAA and AXB plans were studied by applying different sets of published model parameters. Patients were grouped into peripheral and centrally-located tumors to evaluate the impact of tumor location.

Results

PTV dose was lower in the re-calculated AXB plans, as compared to AAA plans. The median differences of PTV(D_95%) were 1.7 Gy (range: 0.3, 6.5 Gy) and 1.0 Gy (range: 0.6, 4.4 Gy) for peripheral tumors and centrally-located tumors, respectively. The median differences of PTV(mean) were 0.4 Gy (range: 0.0, 1.9 Gy) and 0.9 Gy (range: 0.0, 4.3 Gy) for peripheral tumors and centrally-located tumors, respectively. TCP was also found lower in AXB-recalculated plans compared with the AAA plans. The median (range) of the TCP differences for 30 month local control were 1.6 % (0.3 %, 5.8 %) for peripheral tumors and 1.3 % (0.5 %, 3.4 %) for centrally located tumors. The lower TCP is associated with the lower PTV coverage in AXB-recalculated plans. No obvious trend was observed between the calculation-resulted TCP differences and tumor size or location. AAA and AXB yield very similar NTCP on lung pneumonitis according to the LKB model estimation in the present study.

Conclusion

AAA apparently overestimates the PTV dose; the magnitude of resulting difference in calculated TCP was up to 5.8 % in our study. AAA and AXB yield very similar NTCP on lung pneumonitis based on the LKB model parameter sets we used in the present study.

Recent advancements in toxicity prediction following prostate cancer radiotherapy

In external beam radiotherapy for prostate cancer limiting toxicities for dose escalation are bladder and rectum toxicities. Normal tissue complication probability models aim at quantifying the risk of developping adverse events following radiotherapy. These models, originally proposed in the context of uniform irradiation, have evolved to implementations based on the state-of-the-art classification methods which are trained using empirical data. Recently, the use of image processing techniques combined with population analysis methods has led to a new generation of models to understand the risk of normal tissue complications following radiotherapy. This paper overviews those methods in the case of prostate cancer radiation therapy and propose some lines of future research.

Machine Learning Approaches for Predicting Radiation Therapy Outcomes: A Clinician's Perspective

Radiation oncology has always been deeply rooted in modeling, from the early days of isoeffect curves to the contemporary Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) initiative. In recent years, medical modeling for both prognostic and therapeutic purposes has exploded thanks to increasing availability of electronic data and genomics. One promising direction that medical modeling is moving toward is adopting the same machine learning methods used by companies such as Google and Facebook to combat disease. Broadly defined, machine learning is a branch of computer science that deals with making predictions from complex data through statistical models. These methods serve to uncover patterns in data and are actively used in areas such as speech recognition, handwriting recognition, face recognition, "spam" filtering (junk email), and targeted advertising. Although multiple radiation oncology research groups have shown the value of applied machine learning (ML), clinical adoption has been slow due to the high barrier to understanding these complex models by clinicians. Here, we present a review of the use of ML to predict radiation therapy outcomes from the clinician's point of view with the hope that it lowers the "barrier to entry" for those without formal training in ML. We begin by describing 7 principles that one should consider when evaluating (or creating) an ML model in radiation oncology. We next introduce 3 popular ML methods--logistic regression (LR), support vector machine (SVM), and artificial neural network (ANN)--and critique 3 seminal papers in the context of these principles. Although current studies are in exploratory stages, the overall methodology has progressively matured, and the field is ready for larger-scale further investigation.

Impact of spot size on plan quality of spot scanning proton radiosurgery for peripheral brain lesions

PURPOSE

To determine the plan quality of proton spot scanning (SS) radiosurgery as a function of spot size (in-air sigma) in comparison to x-ray radiosurgery for treating peripheral brain lesions.

METHODS

Single-field optimized (SFO) proton SS plans with sigma ranging from 1 to 8 mm, cone-based x-ray radiosurgery (Cone), and x-ray volumetric modulated arc therapy (VMAT) plans were generated for 11 patients. Plans were evaluated using secondary cancer risk and brain necrosis normal tissue complication probability (NTCP).

RESULTS

For all patients, secondary cancer is a negligible risk compared to brain necrosis NTCP. Secondary cancer risk was lower in proton SS plans than in photon plans regardless of spot size (p = 0.001). Brain necrosis NTCP increased monotonically from an average of 2.34/100 (range 0.42/100-4.49/100) to 6.05/100 (range 1.38/100-11.6/100) as sigma increased from 1 to 8 mm, compared to the average of 6.01/100 (range 0.82/100-11.5/100) for Cone and 5.22/100 (range 1.37/100-8.00/100) for VMAT. An in-air sigma less than 4.3 mm was required for proton SS plans to reduce NTCP over photon techniques for the cohort of patients studied with statistical significance (p = 0.0186). Proton SS plans with in-air sigma larger than 7.1 mm had significantly greater brain necrosis NTCP than photon techniques (p = 0.0322).

CONCLUSIONS

For treating peripheral brain lesions--where proton therapy would be expected to have the greatest depth-dose advantage over photon therapy--the lateral penumbra strongly impacts the SS plan quality relative to photon techniques: proton beamlet sigma at patient surface must be small (<7.1 mm for three-beam single-field optimized SS plans) in order to achieve comparable or smaller brain necrosis NTCP relative to photon radiosurgery techniques. Achieving such small in-air sigma values at low energy (<70 MeV) is a major technological challenge in commercially available proton therapy systems.

Radiobiological evaluation of intensity modulated radiation therapy treatments of patients with head and neck cancer: A dual-institutional study

In clinical practice, evaluation of clinical efficacy of treatment planning stems from the radiation oncologist's experience in accurately targeting tumors, while keeping minimal toxicity to various organs at risk (OAR) involved. A more objective, quantitative method may be raised by using radiobiological models. The purpose of this work is to evaluate the potential correlation of OAR-related toxicities to its radiobiologically estimated parameters in simultaneously integrated boost (SIB) intensity modulated radiation therapy (IMRT) plans of patients with head and neck tumors at two institutions. Lyman model for normal tissue complication probability (NTCP) and the Poisson model for tumor control probability (TCP) models were used in the Histogram Analysis in Radiation Therapy (HART) analysis. In this study, 33 patients with oropharyngeal primaries in the head and neck region were used to establish the correlation between NTCP values of (a) bilateral parotids with clinically observed rates of xerostomia, (b) esophagus with dysphagia, and (c) larynx with dysphagia. The results of the study indicated a strong correlation between the severity of xerostomia and dysphagia with Lyman NTCP of bilateral parotids and esophagus, respectively, but not with the larynx. In patients without complications, NTCP values of these organs were negligible. Using appropriate radiobiological models, the presence of a moderate to strong correlation between the severities of complications with NTCP of selected OARs suggested that the clinical outcome could be estimated prior to treatment.

Normal tissue complication probability modeling for cochlea constraints to avoid causing tinnitus after head-and-neck intensity-modulated radiation therapy

Background

Radiation-induced tinnitus is a side effect of radiotherapy in the inner ear for cancers of the head and neck. Effective dose constraints for protecting the cochlea are under-reported. The aim of this study is to determine the cochlea dose limitation to avoid causing tinnitus after head-and-neck cancer (HNC) intensity-modulated radiation therapy (IMRT).

Methods

In total 211 patients with HNC were included; the side effects of radiotherapy were investigated for 422 inner ears in the cohort. Forty-nine of the four hundred and twenty-two samples (11.6 %) developed grade 2+ tinnitus symptoms after IMRT, as diagnosed by a clinician. The Late Effects of Normal Tissues–Subjective, Objective, Management, Analytic (LENT-SOMA) criteria were used for tinnitus evaluation. The logistic and Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) models were used for the analyses.

Results

The NTCP-fitted parameters were TD₅₀ = 46.31 Gy (95 % CI, 41.46–52.50), γ₅₀ = 1.27 (95 % CI, 1.02–1.55), and TD₅₀ = 46.52 Gy (95 % CI, 41.91–53.43), m = 0.35 (95 % CI, 0.30–0.42) for the logistic and LKB models, respectively. The suggested guideline TD₂₀ for the tolerance dose to produce a 20 % complication rate within a specific period of time was TD₂₀ = 33.62 Gy (95 % CI, 30.15–38.27) (logistic) and TD₂₀ = 32.82 Gy (95 % CI, 29.58–37.69) (LKB).

Conclusions

To maintain the incidence of grade 2+ tinnitus toxicity <20 % in IMRT, we suggest that the mean dose to the cochlea should be <32 Gy. However, models should not be extrapolated to other patient populations without further verification and should first be confirmed before clinical implementation.

Relationships between dose to the gastro-intestinal tract and patient-reported symptom domains after radiotherapy for localized prostate cancer

BACKGROUND

Gastrointestinal (GI) morbidity after radiotherapy (RT) for prostate cancer is typically addressed by studying specific single symptoms. The aim of this study was to explore the interplay between domains of patient- reported outcomes (PROs) on GI morbidity, and to what extent these are explained by RT dose to the GI tract.

MATERIAL AND METHODS

The study included men from two Scandinavian studies (N = 211/277) who had undergone primary external beam radiotherapy (EBRT) for localized prostate cancer to 70-78 Gy (2 Gy/fraction). Factor analysis was applied to previously identified PRO-based symptom domains from two study-specific questionnaires. Number of questions: 43; median time to follow-up: 3.6-6.4 years) and dose-response outcome variables were defined from these domains. Dose/volume parameters of the anal sphincter (AS) or the rectum were tested as predictors for each outcome variable using logistic regression with 10-fold cross-validation. Performance was assessed using area under the receiver operating characteristic curve (Az) and model frequency.

RESULTS

Outcome variables from Defecation urgency (number of symptoms: 2-3), Fecal leakage (4-6), Mucous (4), and Pain (3-6) were defined. In both cohorts, intermediate rectal doses predicted Defecation urgency (mean Az: 0.53-0.54; Frequency: 70-75%), and near minimum and low AS doses predicted Fecal leakage (mean Az: 0.63-0.67; Frequency: 83-99%). In one cohort, high AS doses predicted Mucous (mean Az: 0.54; Frequency: 96%), whereas in the other, low AS doses and intermediate rectal doses predicted Pain (mean Az: 0.69; Frequency: 28-82%).

CONCLUSION

We have demonstrated that Defecation urgency, Fecal leakage, Mucous, and Pain following primary EBRT for localized prostate cancer primarily are predicted by intermediate rectal doses, low AS doses, high AS doses, or a combination of low AS and intermediate rectal doses, respectively. This suggests that there is a domain-specific dose-response for the GI tract. To reduce risk of GI morbidity, dose distributions of both the AS region and the rectum should, therefore, be considered when prescribing prostate cancer RT.

Multivariable normal-tissue complication modeling of acute esophageal toxicity in advanced stage non-small cell lung cancer patients treated with intensity-modulated (chemo-)radiotherapy

BACKGROUND AND PURPOSE

The majority of normal-tissue complication probability (NTCP) models for acute esophageal toxicity (AET) in advanced stage non-small cell lung cancer (AS-NSCLC) patients treated with (chemo-)radiotherapy are based on three-dimensional conformal radiotherapy (3D-CRT). Due to distinct dosimetric characteristics of intensity-modulated radiation therapy (IMRT), 3D-CRT based models need revision. We established a multivariable NTCP model for AET in 149 AS-NSCLC patients undergoing IMRT.

MATERIALS AND METHODS

An established model selection procedure was used to develop an NTCP model for Grade ⩾2 AET (53 patients) including clinical and esophageal dose-volume histogram parameters.

RESULTS

The NTCP model predicted an increased risk of Grade ⩾2 AET in case of: concurrent chemoradiotherapy (CCR) [adjusted odds ratio (OR) 14.08, 95% confidence interval (CI) 4.70-42.19; p<0.001], increasing mean esophageal dose [Dmean; OR 1.12 per Gy increase, 95% CI 1.06-1.19; p<0.001], female patients (OR 3.33, 95% CI 1.36-8.17; p=0.008), and ⩾cT3 (OR 2.7, 95% CI 1.12-6.50; p=0.026). The AUC was 0.82 and the model showed good calibration.

CONCLUSIONS

A multivariable NTCP model including CCR, Dmean, clinical tumor stage and gender predicts Grade ⩾2 AET after IMRT for AS-NSCLC. Prior to clinical introduction, the model needs validation in an independent patient cohort.

A dosimetric comparison of ultra-hypofractionated passively scattered proton radiotherapy and stereotactic body radiotherapy (SBRT) in the definitive treatment of localized prostate cancer

BACKGROUND

We compared target and normal tissue dosimetric indices between ultra-hypofractionated passively scattered proton radiotherapy and stereotactic body radiotherapy (SBRT) in the definitive treatment of localized prostate cancer.

MATERIAL AND METHODS

Ten patients were treated definitively for localized prostate cancer with SBRT to a dose of 36.25 Gy in 5 fractions prescribed to a volume encompassing the prostate only. Dose-volume constraints were applied to the rectum, bladder, penile bulb, femoral heads, and prostatic and membranous urethra. Three-field passively scattered proton plans were retrospectively generated using target volumes from the same patients. Dosimetric indices were compared between the SBRT and proton plans using the Wilcoxon signed rank test.

RESULTS

All dose constraints were achieved using both ultra-hypofractionated passively scattered proton and SBRT planning. Proton plans demonstrated significant improvement over SBRT in mean dose delivered to the penile bulb (5.2 CGE vs. 11.4 Gy; p=0.002), rectum (6.7 CGE vs. 10.6 Gy; p=0.002), and membranous urethra (32.2 CGE vs. 34.4 Gy; p=0.006) with improved target homogeneity resulting in a significant reduction in hot spots and volumes of tissue exposed to low doses of radiation. Compared to proton planning, SBRT planning resulted in significant improvement in target conformality with a mean index of 1.17 versus 1.72 (p=0.002), resulting in a dose reduction to the volume of bladder receiving more than 90% of the PD (V32.6, 7.5% vs. 15.9%; p=0.01) and mean dose to the left (7.1 Gy vs. 10.4 CGE; p=0.004) and right (4.0 Gy vs. 10.9 CGE; p=0.01) femoral heads.

CONCLUSION

Target and normal tissue dose constraints for ultra-hypofractionated definitive radiotherapy of localized prostate cancer are readily achieved using both CK SBRT and passively scattered proton-based therapy suggesting feasibility of either modality.

Dose-mass inverse optimization for minimally moving thoracic lesions

In the past decade, several different radiotherapy treatment plan evaluation and optimization schemes have been proposed as viable approaches, aiming for dose escalation or an increase of healthy tissue sparing. In particular, it has been argued that dose-mass plan evaluation and treatment plan optimization might be viable alternatives to the standard of care, which is realized through dose-volume evaluation and optimization. The purpose of this investigation is to apply dose-mass optimization to a cohort of lung cancer patients and compare the achievable healthy tissue sparing to that one achievable through dose-volume optimization. Fourteen non-small cell lung cancer (NSCLC) patient plans were studied retrospectively. The range of tumor motion was less than 0.5 cm and motion management in the treatment planning process was not considered. For each case, dose-volume (DV)-based and dose-mass (DM)-based optimization was performed. Nine-field step-and-shoot IMRT was used, with all of the optimization parameters kept the same between DV and DM optimizations. Commonly used dosimetric indices (DIs) such as dose to 1% the spinal cord volume, dose to 50% of the esophageal volume, and doses to 20 and 30% of healthy lung volumes were used for cross-comparison. Similarly, mass-based indices (MIs), such as doses to 20 and 30% of healthy lung masses, 1% of spinal cord mass, and 33% of heart mass, were also tallied. Statistical equivalence tests were performed to quantify the findings for the entire patient cohort. Both DV and DM plans for each case were normalized such that 95% of the planning target volume received the prescribed dose. DM optimization resulted in more organs at risk (OAR) sparing than DV optimization. The average sparing of cord, heart, and esophagus was 23, 4, and 6%, respectively. For the majority of the DIs, DM optimization resulted in lower lung doses. On average, the doses to 20 and 30% of healthy lung were lower by approximately 3 and 4%, whereas lung volumes receiving 2000 and 3000 cGy were lower by 3 and 2%, respectively. The behavior of MIs was very similar. The statistical analyses of the results again indicated better healthy anatomical structure sparing with DM optimization. The presented findings indicate that dose-mass-based optimization results in statistically significant OAR sparing as compared to dose-volume-based optimization for NSCLC. However, the sparing is case-dependent and it is not observed for all tallied dosimetric endpoints.

Contrasting analytical and data-driven frameworks for radiogenomic modeling of normal tissue toxicities in prostate cancer

BACKGROUND AND PURPOSE

We explore analytical and data-driven approaches to investigate the integration of genetic variations (single nucleotide polymorphisms [SNPs] and copy number variations [CNVs]) with dosimetric and clinical variables in modeling radiation-induced rectal bleeding (RB) and erectile dysfunction (ED) in prostate cancer patients.

MATERIALS AND METHODS

Sixty-two patients who underwent curative hypofractionated radiotherapy (66 Gy in 22 fractions) between 2002 and 2010 were retrospectively genotyped for CNV and SNP rs5489 in the xrcc1 DNA repair gene. Fifty-four patients had full dosimetric profiles. Two parallel modeling approaches were compared to assess the risk of severe RB (Grade⩾3) and ED (Grade⩾1); Maximum likelihood estimated generalized Lyman-Kutcher-Burman (LKB) and logistic regression. Statistical resampling based on cross-validation was used to evaluate model predictive power and generalizability to unseen data.

RESULTS

Integration of biological variables xrcc1 CNV and SNP improved the fit of the RB and ED analytical and data-driven models. Cross-validation of the generalized LKB models yielded increases in classification performance of 27.4% for RB and 14.6% for ED when xrcc1 CNV and SNP were included, respectively. Biological variables added to logistic regression modeling improved classification performance over standard dosimetric models by 33.5% for RB and 21.2% for ED models.

CONCLUSION

As a proof-of-concept, we demonstrated that the combination of genetic and dosimetric variables can provide significant improvement in NTCP prediction using analytical and data-driven approaches. The improvement in prediction performance was more pronounced in the data driven approaches. Moreover, we have shown that CNVs, in addition to SNPs, may be useful structural genetic variants in predicting radiation toxicities.

The radiobiology of hypofractionation

If the α/β ratio is high (e.g. 10 Gy) for tumour clonogen killing, but low (e.g. 3 Gy) for late normal tissue complications, then delivering external beam radiotherapy in a large number (20-30) of small (≈2 Gy) dose fractions should yield the highest 'therapeutic ratio'; this is demonstrated via the linear-quadratic model of cell killing. However, this 'conventional wisdom' is increasingly being challenged, partly by the success of stereotactic body radiotherapy (SBRT) or stereotactic ablative radiotherapy (SABR) extreme hypofractionation regimens of three to five large fractions for early stage non-small cell lung cancer and partly by indications that for certain tumours (prostate, breast) the α/β ratio may be of the same order or even lower than that characterising late complications. It is shown how highly conformal dose delivery combined with quasi-parallel normal tissue behaviour (n close to 1) enables 'safe' hypofractionation; this can be predicted by the (α/β)eff concept for normal tissues. Recent analyses of the clinical outcomes of non-small cell lung cancer radiotherapy covering 'conventional' hyper- to extreme hypofractionation (stereotactic ablative radiotherapy) regimens are consistent with linear-quadratic radiobiology, even at the largest fraction sizes, despite there being theoretical reasons to expect 'LQ violation' above a certain dose. Impairment of re-oxygenation between fractions and the very high (α/β) for hypoxic cells can complicate the picture regarding the analysis of clinical outcomes; it has also been suggested that vascular damage may play a role for very large dose fractions. Finally, the link between high values of (α/β)eff and normal-tissue sparing for quasi-parallel normal tissues, thereby favouring hypofractionation, may be particularly important for proton therapy, but more generally, improved conformality, achieved by whatever technique, can be translated into individualisation of both prescription dose and fraction number via the 'isotoxic' (iso-normal tissue complication probability) concept.