Influence of different dose calculation algorithms on the estimate of NTCP for lung complications

Evaluation of Tumor Control and Normal Tissue Complication Probabilities in Patients Receiving Comprehensive Nodal Irradiation for Left-Sided Breast Cancer

Women with left-sided breast cancer receiving adjuvant radiotherapy have increased incidence of cardiac mortality due to ischemic heart disease; to date, no threshold dose for late cardiac/pulmonary morbidity or mortality has been established. We investigated the likelihood of cardiac death and radiation pneumonitis in women with left-sided breast cancer who received comprehensive lymph node irradiation. The differences in dosimetric parameters between free-breathing (FB) and deep inspiration breath hold (DIBH) techniques were also addressed. Based on NTCP calculations, the probability of cardiac death was significantly reduced with the DIBH compared to the FB technique (p < 0.001). The risk of radiation pneumonitis was not clinically significant. There was no difference in coverage between FB and DIBH plans. Doses to healthy structures were significantly lower in DIBH plan than in FB plan for V20, V30, and ipsilateral total lung volume. Inspiratory gating reduces the dose absorbed by the heart without compromising the target range, thus reducing the likelihood of cardiac death.

A systematic review of the normal tissue complication probability models and parameters: Head and neck cancers treated with conformal radiotherapy

This systematic review study aims to provide comprehensive data on different radiobiological models, parameters, and endpoints used for calculating the normal tissue complication probability (NTCP) based on clinical data from head and neck cancer patients treated with conformal radiotherapy. A systematic literature search was carried out according to the PRISMA guideline for the identification of relevant publications in six electronic databases of Embase, PubMed, Scopus, and Google Scholar to July 2022 using specific keywords in the paper's title and abstract. The initial search resulted in 1368 articles for all organs for the review article about the NTCP parameters. One hundred and seventy-eight articles were accepted for all organs with complete parameters for the mentioned models and finally, 20 head and neck cancer articles were accepted for review. Analysis of the studies shows that the Lyman-Kutcher-Burman (LKB) model properly links the NTCP curve parameters to the postradiotherapy endpoints. In the LKB model for esophagus, the minimum, and maximum corresponding parameters were reported as TD50  = 2.61 Gy with grade ≥3 radiation-induced esophagitis endpoints as the minimum TD50 and TD50  = 68 Gy as the maximum ones. nmin  = 0.06, nmax  = 1.04, mmin  = 0.1, and mmax  = 0.65, respectively. Unfortunately, there was not a wide range of published articles on other organs at risk like ear or cauda equina except Burman et al. (Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys Ther. 1991;21:123-135). Findings suggest that the validation of different radiobiological models and their corresponding parameters need to be investigated in vivo and in vitro for developing a more accurate NTCP model to be used for radiotherapy treatment planning optimization.

Impact of delineation errors on the estimated organ at risk dose and of dose errors on the normal tissue complication probability model

BACKGROUND

Normal tissue complication probability (NTCP) models are often based on doses retrieved from delineated volumes. For retrospective dose-response studies focusing on organs that have not been delineated historically, automatic segmentation might be considered. However, automatic segmentation risks generating considerable delineation errors and knowledge regarding how these errors impact the estimated organ dose is important. Furthermore, organ-at-risk (OAR) dose uncertainties cannot be eliminated and might affect the resulting NTCP model. Therefore, it is also of interest to study how OAR dose errors impact the NTCP modeling results.

PURPOSE

To investigate how random delineation errors of the proximal bronchial tree, heart, and esophagus impact the estimated OAR dose, and to investigate how random errors in the doses used for dose-response modeling affect the estimated NTCPs.

METHODS

We investigated the impact of random delineation errors on the estimated OAR dose using the treatment plans of 39 patients treated with conventionally fractionated radiation therapy of non-small-cell lung cancer. Study-specific reference structures were defined by manually contouring the proximal bronchial tree, heart and esophagus. For each patient and organ, 120 reshaped structures were created by introducing random shifts and margins to the entire reference structure. The mean and near-maximum dose to the reference and reshaped structures were compared. In a separate investigation, the impact of random dose errors on the NTCP model was studied performing dose-response modeling with study sets containing treatment outcomes and OAR doses with and without introduced errors. Universal patient populations with defined population risks, dose-response relationships and distributions of OAR doses were used as ground truth. From such a universal population, we randomly sampled data sets consisting of OAR dose and treatment outcome into reference populations. Study sets of different sizes were created by repeatedly introducing errors to the OAR doses of each reference population. The NTCP models generated with dose errors were compared to the reference NTCP model of the corresponding reference population.

RESULTS

A total of 14 040 reshaped structures with random delineation errors were created. The delineation errors resulted in systematic mean dose errors of less than 1% of the prescribed dose (PD). Mean dose differences above 15% of PD and near-maximum doses differences above 25% of PD were observed for 211 and 457 reshaped structures, respectively. Introducing random errors to OAR doses used for dose-response modeling resulted in systematic underestimations of the median NTCP. For all investigated scenarios, the median differences in NTCP were within 0.1 percentage points (p.p.) when comparing different study sizes.

CONCLUSIONS

Introducing random delineation errors to the proximal bronchial tree, heart and esophagus resulted in mean dose and near-maximum dose differences above 15% and 25% of PD, respectively. We did not observe an association between the dose level and the magnitude of the dose errors. For the scenarios investigated in this study, introducing random errors to OAR doses used for dose-response modeling resulted in systematic underestimations of the median NTCP for reference risks higher than the universal population risk. The median NTCP underestimation was similar for different study sizes, all within 0.1 p.p.

Impact of respiratory motion in dosimetric and clinical advantages for adjuvant left-sided breast radiotherapy

We investigated the organ-sparing effect of the deep inspiration breath hold (DIBH) technique among different levels of lung expansion for left-side breast radiotherapy. This retrospective study enrolled 30 patients who received adjuvant left breast radiotherapy after breast-conserving surgery (BCS). Simulation scans of both DIBH and deep breathing four-dimensional computed tomography (4DCT) were acquired, and three treatment plans were generated for each patient. One plan was based on the DIBH images, and the other two plans were based on the mid-lung expansion (ME) and initial lung expansion (IE) phases retrieved from 4DCT data sets. Dosimetric comparisons and normal tissue complication probability (NTCP) models were conducted. We used image registration for displacement analysis and sought potential factors related to the dose benefit of DIBH. The DIBH plans resulted significantly lower doses to the heart, left ventricle (LV) and left anterior descending coronary artery (LAD), including the high- to low-dose areas, followed by the ME plans and IE plans (p < 0.05). DIBH reduced the risk of long-term cardiac mortality by 40% and radiation pneumonitis of the left lung by 37.96% compared with the IE plans (p < 0.001). The reduction in the mean dose to the heart and LV significantly correlated with anterior displacement of the left lung. The DIBH technique is a feasible tool to provide dosimetric and clinical advantages for adjuvant left-sided breast radiotherapy. Breathing pattern and the level of lung expansion seem to play an important role.

Analysis of early respiratory-related mortality after radiation therapy of non-small-cell lung cancer: feasibility of automatic data extraction for dose-response studies

Purpose: To examine the feasibility of automatic data extraction from clinical radiation therapy (RT) databases at four hospitals to investigate the impact of mean lung dose (MLD) and age on the risk of early respiratory-related death and early overall death for patients treated with RT for non-small-cell lung cancer (NSCLC).Material and methods: We included adult patients with NSCLC receiving curatively intended RT between 2002 and 2017 at four hospitals. A script was developed to automatically extract RT-related data. The cause of death for patients deceased within 180 days of the start of RT was retrospectively assessed. Using logistic regression, the risks of respiratory-related death and of overall death within 90 and 180 days were investigated using MLD and age as variables.Results: Altogether, 1785 patients were included in the analysis of early overall mortality and 1655 of early respiratory-related mortality. The respiratory-related mortalities within 90 and 180 days were 0.9% (15/1655) and 3.6% (60/1655). The overall mortalities within 90 and 180 days were 2.5% (45/1785) and 10.6% (190/1785). Higher MLD and older age were associated with an increased risk of respiratory-related death within 180 days and overall death within 90 and 180 days (all p<.05). For example, the risk of respiratory-related death within 180 days and their 95% confidence interval for patients aged 65 and 75 years with MLDs of 20 Gy was according to our logistic model 3.8% (2.6-5.0%) and 7.7% (5.5-10%), respectively.Conclusions: Automatic data extraction was successfully used to pool data from four hospitals. MLD and age were associated with the risk of respiratory-related death within 180 days of the start of RT and with overall death within 90 and 180 days. A model quantifying the risk of respiratory-related death within 180 days was formulated.

Knowledge-based planning for oesophageal cancers using a model trained with plans from a different treatment planning system

Background: This study aimed to evaluate knowledge-based volume modulated arc therapy (VMAT) plans for oesophageal cancers using a model trained with plans optimised with a different treatment planning system (TPS) and to compare lung dose sparing in two TPSs, Eclipse and RayStation.Materials and methods: A total of 64 patients with stage I-III oesophageal cancers were treated using hybrid VMAT (H-VMAT) plans optimised using RayStation. Among them, 40 plans were used for training the model for knowledge-based planning (KBP) in RapidPlan. The remaining 24 plans were recalculated using RapidPlan to validate the KBP model. H-VMAT plans calculated using RapidPlan were compared with H-VMAT plans optimised using RayStation with respect to planning target volume doses, lung doses, and modulation complexity.Results: In the lung, there were significant differences between the volume ratios receiving doses in excess of 5, 10, and 20 Gy (V₅, V₁₀, and V₂₀). The V₅ for the lung with H-VMAT plans optimised using RapidPlan was significantly higher than that of H-VMAT plans optimised using RayStation (p < .01), with a mean difference of 10%. Compared to H-VMAT plans optimised using RayStation, the V₁₀ and V₂₀ for the lung were significantly lower with H-VMAT plans optimised using RapidPlan (p = .04 and p = .02), with differences exceeding 1.0%. In terms of modulation complexity, the change in beam output at each control point was more constant with H-VMAT plans optimised using RapidPlan than with H-VMAT plans optimised using RayStation. The range of the change with H-VMAT plans optimised using RapidPlan was one third that of H-VMAT plans optimised using RayStation.Conclusion: Two optimisers in Eclipse and RayStation had different dosimetric performance in lung sparing and modulation complexity. RapidPlan could not improve low lung doses, however, it provided an appreciate intermediated doses compared to plans optimised with RayStation.

A treatment planning study comparing IMRT techniques and cyber knife for stereotactic body radiotherapy of low-risk prostate carcinoma

PURPOSE

Comparing radiation treatment plans by using the same safety margins and dose objectives for all techniques, to ascertain the optimal radiation technique for the stereotactic body radiotherapy (SBRT) of low-risk prostate cancer.

MATERIAL AND METHODS

Treatment plans for 27 randomly selected patients were compared using intensity-modulated (IMRT) techniques as Sliding Window (SW), volumetric modulated arc therapy (VMAT), and helical tomotherapy (HT), as well as Cyber Knife (CK) system. The target dose was calculated to 36.25 Gy delivered in five fractions over 1 week. Dosimetric indices for target volume and organs at risk (OAR) as well as normal tissue complication probability (NTCP) of late rectal and urinary bladder toxicities were analyzed.

RESULTS

The CK provided lower homogeneity in the target volume, but higher values for most of the conformity indices compared to the IMRT approaches. The SW demonstrated superior rectum sparing at medium-to-high dose range (V18 Gy - V32.6 Gy) compared to other techniques (p < 0.05). The whole urinary bladder experienced the best shielding by SW and VMAT at the medium dose (V18 Gy, p < 0.05 versus CK), however we obtained no relevant differences between techniques at the high dose. Generally, the CK demonstrated significantly superior rectum and bladder exposure at V18 Gy as compared to HT, SW, and VMAT. For the rectum, mean NTCP values were significantly superior for HT (NTCP = 2.3%, p < 0.05), and for urinary bladder, the NTCP showed no significant advantages for any technique.

CONCLUSION

No absolute dosimetric advantage was revealed to choose between CK or IMRT techniques for the SBRT of low-grade prostate cancer. Using the same safety margins and dose objectives, IMRT techniques demonstrated superior sparing of the rectum and bladder at a medium dose compared to CK. Comparing different IMRT approaches SW displayed superior rectum sparing at a medium-to-high dose range, whereas both SW and RA revealed superior bladder sparing compared to HT. HT demonstrated a significantly lower NTCP outcome compared to CK or IMRT techniques regarding the rectum. Radiation plans can be optimized further by an individual modification of dose objectives independent of the treatment plan strategy.

Comparison of gamma index based on dosimetric error and clinically relevant dose-volume index based on three-dimensional dose prediction in breast intensity-modulated radiation therapy

Background

Measurement-guided dose reconstruction has lately attracted significant attention because it can predict the delivered patient dose distribution. Although the treatment planning system (TPS) uses sophisticated algorithm to calculate the dose distribution, the calculation accuracy depends on the particular TPS used. This study aimed to investigate the relationship between the gamma passing rate (GPR) and the clinically relevant dose–volume index based on the predicted 3D patient dose distribution derived from two TPSs (XiO, RayStation).

Methods

Twenty-one breast intensity-modulated radiation therapy plans were inversely optimized using XiO. With the same plans, both TPSs calculated the planned dose distribution. We conducted per-beam measurements on the coronal plane using a 2D array detector and analyzed the difference in 2D GPRs between the measured and planned doses by commercial software. Using in-house software, we calculated the predicted 3D patient dose distribution and derived the predicted 3D GPR, the predicted per-organ 3D GPR, and the predicted clinically relevant dose–volume indices [dose–volume histogram metrics and the value of the tumor-control probability/normal tissue complication probability of the planning target volume and organs at risk]. The results derived from XiO were compared with those from RayStation.

Results

While the mean 2D GPRs derived from both TPSs were 98.1% (XiO) and 100% (RayStation), the mean predicted 3D GPRs of ipsilateral lung (73.3% [XiO] and 85.9% [RayStation]; p < 0.001) had no correlation with 2D GPRs under the 3% global/3 mm criterion. Besides, this significant difference in terms of referenced TPS between XiO and RayStation could be explained by the fact that the error of predicted V_5Gy of ipsilateral lung derived from XiO (29.6%) was significantly larger than that derived from RayStation (− 0.2%; p < 0.001).

Conclusions

GPR is useful as a patient quality assurance to detect dosimetric errors; however, it does not necessarily contain detailed information on errors. Using the predicted clinically relevant dose–volume indices, the clinical interpretation of dosimetric errors can be obtained. We conclude that a clinically relevant dose–volume index based on the predicted 3D patient dose distribution could add to the clinical and biological considerations in the GPR, if we can guarantee the dose calculation accuracy of referenced TPS.

Estimation of the risk for radiation-induced liver disease following photon- or proton-beam radiosurgery of liver metastases

BACKGROUND

Radiotherapy of liver metastases is commonly being performed with photon-beam based stereotactic body radiation therapy (SBRT). The high risk for radiation-induced liver disease (RILD) is a limiting factor in these treatments. The use of proton-beam based SBRT could potentially improve the sparing of the healthy part of the liver. The aim of this study was to use estimations of normal tissue complication probability (NTCP) to identify liver-metastases patients that could benefit from being treated with intensity-modulated proton therapy (IMPT), based on the reduction of the risk for RILD.

METHODS

Ten liver metastases patients, previously treated with photon-beam based SBRT, were retrospectively planned with IMPT. A CTV-based robust optimisation (accounting for setup and range uncertainties), combined with a PTV-based conventional optimisation, was performed. A robustness criterion was defined for the CTV (V95% > 98% for at least 10 of the 12 simulated scenarios). The NTCP was estimated for different endpoints using the Lyman-Kutcher-Burman model. The ΔNTCP (NTCPIMPT - NTCPSBRT) for RILD was registered for each patient. The patients for which the NTCP (RILD) < 5% were also identified. A generic relative biological effectiveness of 1.1 was assumed for the proton beams.

RESULTS

For all patients, the objectives set for the PTV and the robustness criterion set for the CTV were fulfilled with the IMPT plans. An improved sparing of the healthy part of the liver, right kidney, lungs, spinal cord and the skin was achieved with the IMPT plans, compared to the SBRT plans. Mean liver doses larger than the threshold value of 32 Gy led to NTCP values for RILD exceeding 5% (7 patients with SBRT and 3 patients with the IMPT plans). ΔNTCP values (RILD) ranging between - 98% and - 17% (7 patients) and between 0 and 2% (3 patients), were calculated.

CONCLUSIONS

In this study, liver metastases patients that could benefit from being treated with IMPT, based on the NTCP reductions, were identified. The clinical implementation of such a model-based approach to select liver metastases patients to proton therapy needs to be made with caution while considering the uncertainties involved in the NTCP estimations.

Dosimetrical and radiobiological approach to manage the dosimetric shift in the transition of dose calculation algorithm in radiation oncology: how to improve high quality treatment and avoid unexpected outcomes?

BACKGROUND

For a given prescribed dose of radiotherapy, with the successive generations of dose calculation algorithms, more monitor units (MUs) are generally needed. This is due to the implementation of successive improvements in dose calculation: better heterogeneity correction and more accurate estimation of secondary electron transport contribution. More recently, there is the possibility to report the dose-to-medium, physically more accurate compared to the dose-to-water as the reference one. This last point is a recent concern and the main focus of this study.

METHODS

In this paper, we propose steps for a general analysis procedure to estimate the dosimetric alterations, and the potential clinical changes, between a reference algorithm and a new one. This includes dosimetric parameters, gamma index, radiobiology indices based on equivalent uniform dose concept and statistics with bootstrap simulation. Finally, we provide a general recommendation on the clinical use of new algorithms regarding the dose prescription or dose limits to the organs at risks.

RESULTS

The dosimetrical and radiobiological data showed a significant effect, which might exceed 5-10%, of the calculation method on the dose the distribution and clinical outcomes for lung cancer patients. Wilcoxon signed rank paired comparisons indicated that the delivered dose in MUs was significantly increased (> 2%) using more advanced dose calculation methods as compared to the reference one.

CONCLUSION

This paper illustrates and explains the use of dosimetrical, radiobiologcal and statistical tests for dosimetric comparisons in radiotherapy. The change of dose calculation algorithm may induce a dosimetric shift, which has to be evaluated by the physicists and the oncologists. This includes the impact on tumor control and on the risk of toxicity based on normal tissue dose constraints. In fact, the alteration in dose distribution makes it hard to keep exactly the same tumor control probability along with the same normal tissue complication probability.

Image guidance in proton therapy for lung cancer

Proton therapy is a promising but challenging treatment modality for the management of lung cancer. The technical challenges are due to respiratory motion, low dose tolerance of adjacent normal tissue and tissue density heterogeneity. Different imaging modalities are applied at various steps of lung proton therapy to provide information on target definition, target motion, proton range, patient setup and treatment outcome assessment. Imaging data is used to guide treatment design, treatment delivery, and treatment adaptation to ensure the treatment goal is achieved. This review article will summarize and compare various imaging techniques that can be used in every step of lung proton therapy to address these challenges.

Using gEUD based plan analysis method to evaluate proton vs. photon plans for lung cancer radiation therapy

The goal of this study was to exam the efficacy of current DVH based clinical guidelines draw from photon experience for lung cancer radiation therapy on proton therapy. Comparison proton plans and IMRT plans were generated for 10 lung patients treated in our proton facility. A gEUD based plan evaluation method was developed for plan evaluation. This evaluation method used normal lung gEUD(a) curve in which the model parameter "a" was sampled from the literature reported value. For all patients, the proton plans delivered lower normal lung V5 Gy with similar V20 Gy and similar target coverage. Based on current clinical guidelines, proton plans were ranked superior to IMRT plans for all 10 patients. However, the proton and IMRT normal lung gEUD(a) curves crossed for 8 patients within the tested range of "a", which means there was a possibility that proton plan would be worse than IMRT plan for lung sparing. A concept of deficiency index (DI) was introduced to quantify the probability of proton plans doing worse than IMRT plans. By applying threshold on DI, four patients' proton plan was ranked inferior to the IMRT plan. Meanwhile if a threshold to the location of curve crossing was applied, 6 patients' proton plan was ranked inferior to the IMRT plan. The contradictory ranking results between the current clinical guidelines and the gEUD(a) curve analysis demonstrated there is potential pitfalls by applying photon experience directly to the proton world. A comprehensive plan evaluation based on radio-biological models should be carried out to decide if a lung patient would really be benefit from proton therapy.

Intensity-modulated radiation therapy versus volumetric-modulated arc therapy in non-small cell lung cancer: assessing the risk of radiation pneumonitis

Purpose

This study aimed to compare intensity-modulated radiation therapy (IMRT) and volumetric-modulated arc therapy (VMAT) regarding plan quality and healthy lung sparing, in stage III non-small cell lung cancer (NSCLC) patients.

Materials and methods

The plans of 60 patients were allocated either to the IMRT (n=30) or the VMAT (n=30) group. The dose prescribed to the planning target volume (PTV) was evaluated at the 95% level and the mean lung dose (MLD) and the healthy lung receiving 5, 10 and 20 Gy (V5, V10and V20, respectively) were analysed. The normal tissue complication probability (NTCP) for radiation pneumonitis was calculated with the Lyman–Kutcher–Burman model.

Results

Both techniques achieved comparable results for target coverage (V95%=97·87 versus 97·18%,p>0·05) and homogeneity. The MLD (15·57 versus 16·98 Gy,p>0·05), V5(60·35 versus 67·25%,p>0·05) and V10(45·22 versus 53·14%,p=0·011) were lower for IMRT, whereas VMAT reduced V20(26·44 versus 25·90%,p>0·05). The NTCP for radiation pneumonitis was higher for VMAT, but no statistical significance was observed (11·07 versus 12·75,p>0·05).

Conclusion

Both techniques seemed suitable for NSCLC treatment, but IMRT presented better results regarding lung sparing thus being beneficial in reducing the risk of radiation-induced pneumonitis.

Normal lung sparing Tomotherapy technique in stage III lung cancer

PURPOSE

Radiation pneumonitis (RP) has been a challenging obstacle in treating stage III lung cancer patients. Beam angle optimization (BAO) technique for Tomotherapy was developed to reduce the normal lung dose for stage III non-small cell lung cancer (NSCLC). Comparative analyses on plan quality by 3 different Intensity-modulated radiation therapy (IMRT) methods with BAO were done.

MATERIALS AND METHODS

Ten consecutive stage IIIB NSCLC patients receiving linac-based static IMRT (L-IMRT) with total 66 Gy in 33 fractions to the PTV were selected. Two additional Tomotherapy-based IMRT plans (helical beam (TH-IMRT) and static beam (TD-IMRT)) were generated on each patient. To reduce the normal lung dose, Beam angles were optimized by using complete and directional block functions in Tomotherapy based on knowledge based statistical analysis. Plan quality was compared with target coverage, normal organ sparing capability, and normal tissue complication probability (NTCP). Actual beam delivery times and risk of RP related with planning target volume (PTV) were also evaluated.

RESULTS

The best PTV coverage measured by conformity index and homogeneity index was achievable by TH-IMRT (0.82 and 1.06), followed by TD-IMRT (0.81 and 1.07) and L-IMRT (0.75 and 1.08). Mean lung dose was the lowest in TH-IMRT plan followed by TD-IMRT and L-IMRT, all of which were ≤20 Gy. TH-IMRT plan could significantly lower the lung volumes receiving low to medium dose levels: V5~30 when compared to L-IMRT plan; and V5~20 when compared to TD-IMRT plan, respectively. TD-IMRT plan was significantly better than L-IMRT with respects to V20 and V30 and there was no significant difference with respect to V40 among three plans. The NTCP of the lung was the lowest in TH-IMRT plan, followed by TD-IMRT and L-IMRT (6.42% vs. 6.53% vs. 8.11%). Beam delivery time was the shortest in TD-IMRT plan followed by L-IMRT. As PTV length increased, NTCP and Mean lung dose proportionally increased significantly in all three plans.

CONCLUSION

Advantageous profiles by TH-IMRT could be achieved by BAO by complete and directional block functions. Current observation could help radiation oncologists to make wise selection of IMRT method for stage IIIB NSCLC.

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.

Evaluation of the radiobiological gamma index with motion interplay in tangential IMRT breast treatment

The purpose of this study was to evaluate the impact of the motion interplay effect in early-stage left-sided breast cancer intensity-modulated radiation therapy (IMRT), incorporating the radiobiological gamma index (RGI). The IMRT dosimetry for various breathing amplitudes and cycles was investigated in 10 patients. The predicted dose was calculated using the convolution of segmented measured doses. The physical gamma index (PGI) of the planning target volume (PTV) and the organs at risk (OAR) was calculated by comparing the original with the predicted dose distributions. The RGI was calculated from the PGI using the tumor control probability (TCP) and the normal tissue complication probability (NTCP). The predicted mean dose and the generalized equivalent uniform dose (gEUD) to the target with various breathing amplitudes were lower than the original dose (P < 0.01). The predicted mean dose and gEUD to the OARs with motion were higher than for the original dose to the OARs (P < 0.01). However, the predicted data did not differ significantly between the various breathing cycles for either the PTV or the OARs. The mean RGI gamma passing rate for the PTV was higher than that for the PGI (P < 0.01), and for OARs, the RGI values were higher than those for the PGI (P < 0.01). The gamma passing rates of the RGI for the target and the OARs other than the contralateral lung differed significantly from those of the PGI under organ motion. Provided an NTCP value <0.05 is considered acceptable, it may be possible, by taking breathing motion into consideration, to escalate the dose to achieve the PTV coverage without compromising the TCP.

Correlation between pneumonitis risk in radiation oncology and lung density measured with X-ray computed tomography

BACKGROUND

The risk of toxicity with radiation oncology for lung cancer limits the maximal radiation dose that can be delivered to thoracic tumors. This study aims at investigating the correlation between normal tissue complication probability (NTCP) and physical lung density by analyzing the computed tomography (CT) scan imaging used for radiotherapy dose planning.

METHODS

Data from CT of lung cancer patients (n=10), treated with three dimensional radiotherapy, were selected for this study. The dose was calculated using analytical anisotropic algorithm (AAA). Dose volume histograms (DVH) for healthy lung (lung excluding targets) were calculated. The NTCP for lung radiation induced pneumonitis was computed using initial radiobiological parameters from Lyman-Kutcher and Burman (LKB) model and readjusted parameters for AAA, with α/β=3. The correlation coefficient "rho" was calculated using Spearman's rank test. The bootstrap method was used to estimate the 95% confidence interval (95% CI). Wilcoxon paired test was used to calculate P values.

RESULTS

Bootstrapping simulation revealed significant difference between NTCP computed with the initial radiobiological parameters and that computed with the parameters readjusted for AAA (P=0.03). The results of simulations based on 1,000 replications showed no correlation for NTCP with density, with "rho" <0.3.

CONCLUSIONS

For a given set of patients, we assessed the correlation between NTCP and lung density using bootstrap analysis. The lack of correlation could result either from a very accurate dose calculation, by AAA, whatever the lung density yielding a NTCP result only dependant of the dose and not any more of the density; or to the very limited range of natural variation of relative electronic density (0.15 to 0.20) observed in this small series of patients. Another important parameter is the bootstrap simulation with 1,000 random samplings may have underestimated the correlation, since the initial data (n=10) showed a weak correlation.

Quantitative comparison of dose distribution in radiotherapy plans using 2D gamma maps and X-ray computed tomography

BACKGROUND

The advanced dose calculation algorithms implemented in treatment planning system (TPS) have remarkably improved the accuracy of dose calculation especially the modeling of electrons transport in the low density medium. The purpose of this study is to evaluate the use of 2D gamma (γ) index to quantify and evaluate the impact of the calculation of electrons transport on dose distribution for lung radiotherapy.

METHODS

X-ray computed tomography images were used to calculate the dose for twelve radiotherapy treatment plans. The doses were originally calculated with Modified Batho (MB) 1D density correction method, and recalculated with anisotropic analytical algorithm (AAA), using the same prescribed dose. Dose parameters derived from dose volume histograms (DVH) and target coverage indices were compared. To compare dose distribution, 2D γ-index was applied, ranging from 1%/1 mm to 6%/6 mm. The results were displayed using γ-maps in 2D. Correlation between DVH metrics and γ passing rates was tested using Spearman's rank test and Wilcoxon paired test to calculate P values.

RESULTS

the plans generated with AAA predicted more heterogeneous dose distribution inside the target, with P<0.05. However, MB overestimated the dose predicting more coverage of the target by the prescribed dose. The γ analysis showed that the difference between MB and AAA could reach up to ±10%. The 2D γ-maps illustrated that AAA predicted more dose to organs at risks, as well as lower dose to the target compared to MB.

CONCLUSIONS

Taking into account of the electrons transport on radiotherapy plans showed a significant impact on delivered dose and dose distribution. When considering the AAA represent the true cumulative dose, a readjusting of the prescribed dose and an optimization to protect the organs at risks should be taken in consideration in order to obtain the better clinical outcome.

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.

Normal tissue considerations and dose-volume constraints in the moderately hypofractionated treatment of non-small cell lung cancer

Hypofractionated radiation therapy (RT) regimes in non-small cell lung cancer (NSCLC) have become increasingly popular with a number of international trials currently underway. The majority of the dose-volume-constraints (DVCs) published in the literature refer to conventional 2Gy per fraction deliveries. Here relevant organs-at-risk (OARs) are identified and available dose-volume constraint data discussed and summarised for moderately hypofractionated NSCLC regimes. The OARs examined include lung, brachial plexus, heart, oesophagus, airway and spinal cord. Where available the toxicity rates are also reported with all data summarised tabulated to aid its use in the clinic.

Comparison of doses and NTCP to risk organs with enhanced inspiration gating and free breathing for left-sided breast cancer radiotherapy using the AAA algorithm

Background

The purpose of this study was to investigate the potential dose reduction to the heart, left anterior descending (LAD) coronary artery and the ipsilateral lung for patients treated with tangential and locoregional radiotherapy for left-sided breast cancer with enhanced inspiration gating (EIG) compared to free breathing (FB) using the AAA algorithm. The radiobiological implication of such dose sparing was also investigated.

Methods

Thirty-two patients, who received tangential or locoregional adjuvant radiotherapy with EIG for left-sided breast cancer, were retrospectively enrolled in this study. Each patient was CT-scanned during FB and EIG. Similar treatment plans, with comparable target coverage, were created in the two CT-sets using the AAA algorithm. Further, the probability of radiation induced cardiac mortality and pneumonitis were calculated using NTCP models.

Results

For tangential treatment, the median V_25Gy for the heart and LAD was decreased for EIG from 2.2% to 0.2% and 40.2% to 0.1% (p < 0.001), respectively, whereas there was no significant difference in V_20Gy for the ipsilateral lung (p = 0.109). For locoregional treatment, the median V_25Gy for the heart and LAD was decreased for EIG from 3.3% to 0.2% and 51.4% to 5.1% (p < 0.001), respectively, and the median ipsilateral lung V_20Gy decreased from 27.0% for FB to 21.5% (p = 0.020) for EIG. The median excess cardiac mortality probability decreased from 0.49% for FB to 0.02% for EIG (p < 0.001) for tangential treatment and from 0.75% to 0.02% (p < 0.001) for locoregional treatment. There was no significant difference in risk of radiation pneumonitis for tangential treatment (p = 0.179) whereas it decreased for locoregional treatment from 6.82% for FB to 3.17% for EIG (p = 0.004).

Conclusions

In this study the AAA algorithm was used for dose calculation to the heart, LAD and left lung when comparing the EIG and FB techniques for tangential and locoregional radiotherapy of breast cancer patients. The results support the dose and NTCP reductions reported in previous studies where dose calculations were performed using the pencil beam algorithm.