Biologically effective dose distribution based on the linear quadratic model and its clinical relevance

AAPM Task Group Report 267: A joint AAPM GEC-ESTRO report on biophysical models and tools for the planning and evaluation of brachytherapy

Brachytherapy utilizes a multitude of radioactive sources and treatment techniques that often exhibit widely different spatial and temporal dose delivery patterns. Biophysical models, capable of modeling the key interacting effects of dose delivery patterns with the underlying cellular processes of the irradiated tissues, can be a potentially useful tool for elucidating the radiobiological effects of complex brachytherapy dose delivery patterns and for comparing their relative clinical effectiveness. While the biophysical models have been used largely in research settings by experts, it has also been used increasingly by clinical medical physicists over the last two decades. A good understanding of the potentials and limitations of the biophysical models and their intended use is critically important in the widespread use of these models. To facilitate meaningful and consistent use of biophysical models in brachytherapy, Task Group 267 (TG-267) was formed jointly with the American Association of Physics in Medicine (AAPM) and The Groupe Européen de Curiethérapie and the European Society for Radiotherapy & Oncology (GEC-ESTRO) to review the existing biophysical models, model parameters, and their use in selected brachytherapy modalities and to develop practice guidelines for clinical medical physicists regarding the selection, use, and interpretation of biophysical models. The report provides an overview of the clinical background and the rationale for the development of biophysical models in radiation oncology and, particularly, in brachytherapy; a summary of the results of literature review of the existing biophysical models that have been used in brachytherapy; a focused discussion of the applications of relevant biophysical models for five selected brachytherapy modalities; and the task group recommendations on the use, reporting, and implementation of biophysical models for brachytherapy treatment planning and evaluation. The report concludes with discussions on the challenges and opportunities in using biophysical models for brachytherapy and with an outlook for future developments.

The impact of brachytherapy boost for anal canal cancers in the era of de-escalation treatments

PURPOSE

To analyze clinical outcomes of high-dose-rate (HDR) interstitial brachytherapy boost (ISBT) after external beam radiation therapy (EBRT) or chemoradiotherapy (CRT) for the treatment of anal canal cancers (ACC).

METHODS AND MATERIALS

A total of 78 patients with ACC were treated at our institution by ISBT. Local Control (LC), disease-free survival (DFS), overall survival (OS), colostomy-free survival (CFS) and toxicity rates were analyzed.

RESULTS

With a median followup (FU) of 59.8 months (95% CI [55.8-64.2]), six (7.7%) local recurrences with 2 patients (2.6%) having persistent disease at 3 months were observed. The 5-year rate of LC for the entire population was 92% [83-96%]. The 5-year DFS rate was 86% [76-93%]. The 5-year OS was 96% [88-99%]. In the univariate analysis, chemotherapy was significantly associated with morbidity grade ≥2. Late digestive toxicity grade ≥3 was reported in 8.9% patients, 1 patient underwent colostomy due to toxicity. The 5-year CFS rate was 88% [79-94%].

CONCLUSIONS

HDR interstitial brachytherapy boost provide excellent rates of tumor control and colostomy-free survival with a favorable profile of GI toxicity. Continence in anal cancer survivors is a challenge and the boost technique must be discussed in a multidisciplinary approach as part of de-escalation treatments.

High-Dose Once-Daily Thoracic Radiotherapy in Limited-Stage Small-Cell Lung Cancer: CALGB 30610 (Alliance)/RTOG 0538

PURPOSE

Although level 1 evidence supports 45-Gy twice-daily radiotherapy as standard for limited-stage small-cell lung cancer, most patients receive higher-dose once-daily regimens in clinical practice. Whether increasing radiotherapy dose improves outcomes remains to be prospectively demonstrated.

METHODS

This phase III trial, CALGB 30610/RTOG 0538 (ClinicalTrials.gov identifier: NCT00632853), was conducted in two stages. In the first stage, patients with limited-stage disease were randomly assigned to receive 45-Gy twice-daily, 70-Gy once-daily, or 61.2-Gy concomitant-boost radiotherapy, starting with either the first or second (of four total) chemotherapy cycles. In the second stage, allocation to the 61.2-Gy arm was discontinued following planned interim toxicity analysis, and the study continued with two remaining arms. The primary end point was overall survival (OS) in the intention-to-treat population.

RESULTS

Trial accrual opened on March 15, 2008, and closed on December 1, 2019. All patients randomly assigned to 45-Gy twice-daily (n = 313) or 70-Gy once-daily radiotherapy (n = 325) are included in this analysis. After a median follow-up of 4.7 years, OS was not improved on the once-daily arm (hazard ratio for death, 0.94; 95% CI, 0.76 to 1.17; P = .594). Median survival is 28.5 months for twice-daily treatment, and 30.1 months for once-daily treatment, with 5-year OS of 29% and 32%, respectively. Treatment was tolerable, and the frequency of severe adverse events, including esophageal and pulmonary toxicity, was similar on both arms.

CONCLUSION

Although 45-Gy twice-daily radiotherapy remains the standard of care, this study provides the most robust information available to help guide the choice of thoracic radiotherapy regimen for patients with limited-stage small-cell lung cancer.

The Effect of Atm Loss on Radiosensitivity of a Primary Mouse Model of Pten-Deleted Brainstem Glioma

Diffuse midline gliomas arise in the brainstem and other midline brain structures and cause a large proportion of childhood brain tumor deaths. Radiation therapy is the most effective treatment option, but these tumors ultimately progress. Inhibition of the phosphoinositide-3-kinase (PI3K)-like kinase, ataxia-telangiectasia mutated (ATM), which orchestrates the cellular response to radiation-induced DNA damage, may enhance the efficacy of radiation therapy. Diffuse midline gliomas in the brainstem contain loss-of-function mutations in the tumor suppressor PTEN, or functionally similar alterations in the phosphoinositide-3-kinase (PI3K) pathway, at moderate frequency. Here, we sought to determine if ATM inactivation could radiosensitize a primary mouse model of brainstem glioma driven by Pten loss. Using Cre/loxP recombinase technology and the RCAS/TVA retroviral gene delivery system, we established a mouse model of brainstem glioma driven by Pten deletion. We find that Pten-null brainstem gliomas are relatively radiosensitive at baseline. In addition, we show that deletion of Atm in the tumor cells does not extend survival of mice bearing Pten-null brainstem gliomas after focal brain irradiation. These results characterize a novel primary mouse model of PTEN-mutated brainstem glioma and provide insights into the mechanism of radiosensitization by ATM deletion, which may guide the design of future clinical trials.

Mechanisms and Review of Clinical Evidence of Variations in Relative Biological Effectiveness in Proton Therapy

Proton therapy is increasingly being used as a radiation therapy modality. There is uncertainty about the biological effectiveness of protons relative to photon therapies as it depends on several physical and biological parameters. Radiation oncology currently applies a constant and generic value for the relative biological effectiveness (RBE) of 1.1, which was chosen conservatively to ensure tumor coverage. The use of a constant value has been challenged particularly when considering normal tissue constraints. Potential variations in RBE have been assessed in several published reviews but have mostly focused on data from clonogenic cell survival experiments with unclear relevance for clinical proton therapy. The goal of this review is to put in vitro findings in relation to clinical observations. Relevant in vivo pathways determining RBE for tumors and normal tissues are outlined, including not only damage to tumor cells and parenchyma but also vascular damage and immune response. Furthermore, the current clinical evidence of varying RBE is reviewed. The assessment can serve as guidance for treatment planning, personalized dose prescriptions, and outcome analysis.

Accelerated hyper-versus normofractionated radiochemotherapy with temozolomide in patients with glioblastoma: a multicenter retrospective analysis

Background and purpose

The standard treatment of glioblastoma patients consists of surgery followed by normofractionated radiotherapy (NFRT) with concomitant and adjuvant temozolomide chemotherapy. Whether accelerated hyperfractionated radiotherapy (HFRT) yields comparable results to NFRT in combination with temozolomide has only sparsely been investigated. The objective of this study was to compare NFRT with HFRT in a multicenter analysis.

Materials and methods

A total of 484 glioblastoma patients from four centers were retrospectively pooled and analyzed. Three-hundred-ten and 174 patients had been treated with NFRT (30 × 1.8 Gy or 30 × 2 Gy) and HFRT (37 × 1.6 Gy or 30 × 1.8 Gy twice/day), respectively. The primary outcome of interest was overall survival (OS) which was correlated with patient-, tumor- and treatment-related variables via univariable and multivariable Cox frailty models. For multivariable modeling, missing covariates were imputed using multiple imputation by chained equations, and a sensitivity analysis was performed on the complete-cases-only dataset.

Results

After a median follow-up of 15.7 months (range 0.8–88.6 months), median OS was 16.9 months (15.0–18.7 months) in the NFRT group and 14.9 months (13.2–17.3 months) in the HFRT group (p = 0.26). In multivariable frailty regression, better performance status, gross-total versus not gross-total resection, MGMT hypermethylation, IDH mutation, smaller planning target volume and salvage therapy were significantly associated with longer OS (all p < 0.01). Treatment differences (HFRT versus NFRT) had no significant effect on OS in either univariable or multivariable analysis.

Conclusions

Since HFRT with temozolomide was not associated with worse OS, we assume HFRT to be a potential option for patients wishing to shorten their treatment time.

Tumor-immune ecosystem dynamics define an individual Radiation Immune Score to predict pan-cancer radiocurability

Radiotherapy efficacy is the result of radiation-mediated cytotoxicity coupled with stimulation of antitumor immune responses. We develop an in silico 3-dimensional agent-based model of diverse tumor-immune ecosystems (TIES) represented as anti- or pro-tumor immune phenotypes. We validate the model in 10,469 patients across 31 tumor types by demonstrating that clinically detected tumors have pro-tumor TIES. We then quantify the likelihood radiation induces antitumor TIES shifts toward immune-mediated tumor elimination by developing the individual Radiation Immune Score (iRIS). We show iRIS distribution across 31 tumor types is consistent with the clinical effectiveness of radiotherapy, and in combination with a molecular radiosensitivity index (RSI) combines to predict pan-cancer radiocurability. We show that iRIS correlates with local control and survival in a separate cohort of 59 lung cancer patients treated with radiation. In combination, iRIS and RSI predict radiation-induced TIES shifts in individual patients and identify candidates for radiation de-escalation and treatment escalation. This is the first clinically and biologically validated computational model to simulate and predict pan-cancer response and outcomes via the perturbation of the TIES by radiotherapy.

The effects of PSA kinetics on the outcome of hypofractionated salvage radiotherapy for biochemical recurrence of prostate cancer after prostatectomy

The feasibility and efficacy of hypofractionated salvage radiotherapy (HS-RT) for prostate cancer (PC) with biochemical recurrence (BR) after prostatectomy, and the usefulness of prostate-specific antigen (PSA) kinetics as a predictor of BR, were evaluated in 38 patients who received HS-RT without androgen deprivation therapy between May 2009 and January 2017. Their median age, PSA level and PSA doubling time (PSA-DT) at the start of HS-RT were 68 (53-74) years, 0.28 (0.20-0.79) ng/ml and 7.7 (2.3-38.5) months, respectively. A total dose of 60 Gy in 20 fractions (three times a week) was three-dimensionally delivered to the prostate bed. After a median follow-up of 62 (30-100) months, 19 (50%) patients developed a second BR after HS-RT, but only 1 patient died before the last follow-up. The 5-year overall survival and BR-free survival rates were 97.1 and 47.4%, respectively. Late grade 2 gastrointestinal and genitourinary morbidities were observed in 0 and 5 (13%) patients, respectively. The PSA level as well as pathological T-stage and surgical margin status were regarded as significant predictive factors for a second BR by multivariate analysis. BR developed within 6 months after HS-RT in 11 (85%) of 13 patients with a PSA-DT < 10 months compared with 1 (17%) of 6 with a PSA-DT ≥ 10 months (median time to BR: 3 vs 14 months, P < 0.05). Despite the small number of patients, our HS-RT protocol seems feasible, and PSA kinetics may be useful for predicting the risk of BR and determining the appropriate follow-up schedule.

Radiotherapy for local recurrence of non-small-cell lung cancer after lobectomy and lymph node dissection-can local recurrence be radically cured by radiation?

BACKGROUND

There is no standard therapeutic approach for local recurrence of non-small cell lung cancer (NSCLC) after complete resection. We investigated the outcomes of radiotherapy (RT) for patients with local recurrence.

METHODS

We reviewed 46 patients who underwent curative-intent RT for local recurrence after lobectomy or pneumonectomy accompanied with mediastinal lymph node dissection between 2002 and 2014. We analyzed overall survival (OS), progression-free survival (PFS), local control, tumour response and the re-recurrence pattern.

RESULTS

Among the 46 patients, 16 received concurrent chemotherapy. The median follow-up period was 48 months. The response rate was 91%. The 5-year OS and local control rates were 47.9 and 65.3%, respectively, and the 5-year PFS rate was 22.8%. Female sex and complete response to radiation were favourable prognostic factors. Of the 33 patients with recurrence after radiation, 32 (97%) had distant metastasis.

CONCLUSIONS

Although RT for local recurrence has high efficacy, distant relapse after radiation remains a major issue. Therefore, combination systemic therapy for local recurrence at any site should be further investigated. Since it is difficult to achieve a radical cure for local recurrence using RT, further study, for the administration of post-operative adjuvant therapy, is recommended.

Radiation Dose Escalation is Crucial in Anti-CTLA-4 Antibody Therapy to Enhance Local and Distant Antitumor Effect in Murine Osteosarcoma

We previously reported that a combination of 10 Gy of X-ray irradiation and dual immune checkpoint blockade with anti-CTLA-4 (C4) and anti-PD-L1 antibodies produced a significant shrinkage of irradiated and unirradiated tumors (abscopal effect) and prolonged overall survival. However, the optimal radiation delivery regimen combined with single immune checkpoint blockade of C4 for inducing a maximum systemic antitumor response still remains unclear, particularly for patients with osteosarcoma. We used syngeneic C3H mice that were subcutaneously injected with LM8 osteosarcoma cells into both legs. C4 was administered three times, and one side of the tumor was irradiated by X-ray beams. The optimal radiation dose required to induce the abscopal effect was explored with a focus on the induction of the type-I interferon pathway. Radiation delivered in a single fraction of 10 Gy, 4.5 Gy × 3 fractions (fx), and 2 Gy × 8 fx with C4 failed to produce significant inhibition of unirradiated tumor growth compared with monotherapy with C4. Dose escalation to 16 Gy in a single fraction, or the equivalent hypofractionated dose of 8 Gy × 3 fx, which significantly increased secretion of IFN-β in vitro, produced a dramatic regression of both irradiated and unirradiated tumors and prolonged overall survival in combination with C4. Furthermore, irradiation at 16 Gy in both a single fraction and 8 Gy × 3 fx diminished regulatory T cells in the unirradiated tumor microenvironment. These results suggest that total dose escalation of radiation is crucial in C4 therapy to enhance the antitumor response in both local and distant tumors and prolonged overall survival regardless of fractionation for osteosarcoma.

Immunologic Consequences of Sequencing Cancer Radiotherapy and Surgery

PURPOSE

Early-stage cancers are routinely treated with surgery followed by radiotherapy (SR). Radiotherapy before surgery (RS) has been widely ignored for some cancers. We evaluate overall survival (OS) and disease-free survival (DFS) with SR and RS for different cancer types and simulate the plausibility of RS- and SR-induced antitumor immunity contributing to outcomes.

MATERIALS AND METHODS

We analyzed a SEER data set of early-stage cancers treated with SR or RS. OS and DFS were calculated for cancers with sufficient numbers for statistical power (cancers of lung and bronchus, esophagus, rectum, cervix uteri, corpus uteri, and breast). We simulated the immunologic consequences of SR, RS, and radiotherapy alone in a mathematical model of tumor-immune interactions.

RESULTS

RS improved OS for cancers with low 20-year survival rates (lung: hazard ratio [HR], 0.88; P = .046) and improved DFS for cancers with higher survival (breast: HR = 0.64; P < .001). For rectal cancer, with intermediate 20-year survival, RS improved both OS (HR = 0.89; P = .006) and DFS (HR = 0.86; P = .04). Model simulations suggested that RS could increase OS by eliminating cancer for a broader range of model parameters and radiotherapy-induced antitumor immunity compared with SR for selected parameter combinations. This could create an immune memory that may explain increased DFS after RS for certain cancers.

CONCLUSION

Study results suggest plausibility that radiation to the bulk of the tumor could induce a more robust immune response and better harness the synergy of radiotherapy and antitumor immunity than postsurgical radiation to the tumor bed. This exploratory study provides motivation for prospective evaluation of immune activation of RS versus SR in controlled clinical studies.

Mathematical modeling of radiotherapy cancer treatment using Caputo fractional derivative

BACKGROUND

This paper presents a mathematical model that simulates a radiotherapy cancer treatment process. The model takes into consideration two important radiobiological factors, which are repair and repopulation of cells. The model was used to simulate the fractionated treatment process of six patients. The results gave the population changes in the cells and the final volumes of the normal and cancer cells.

METHOD

The model was formulated by integrating the Caputo fractional derivative with the previous cancer treatment model. Thereafter, the linear-quadratic with the repopulation model was coupled into the model to account for the cells' population decay due to radiation. The treatment process was then simulated with numerical variables, numerical parameters, and radiation parameters. The numerical parameters which included the proliferation coefficients of the cells, competition coefficients of the cells, and the perturbation constant of the normal cells were obtained from previous literature. The radiation and numerical parameters were obtained from reported clinical data of six patients treated with radiotherapy. The patients had tumor volumes of 24.1cm³, 17.4cm³, 28.4cm³, 18.8cm³, 30.6cm³, and 12.6cm³ with fractionated doses of 2 Gy for the first two patients and 1.8 Gy for the other four. The initial tumor volumes were used to obtain initial populations of cells after which the treatment process was simulated in MATLAB. Subsequently, a global sensitivity analysis was done to corroborate the model with clinical data. Finally, 96 radiation protocols were simulated by using the biologically effective dose formula. These protocols were used to obtain a regression equation connecting the value of the Caputo fractional derivative with the fractionated dose.

RESULTS

The final tumor volumes, from the results of the simulations, were 3.58cm³, 8.61cm³, 5.68cm³, 4.36cm³, 5.75cm³, and 6.12cm³, while those of the normal cells were 23.87cm³, 17.29cm³, 28.17cm³, 18.68cm³, 30.33cm³, and 12.55cm³. The sensitivity analysis showed that the most sensitive model factors were the value of the Caputo fractional derivative and the proliferation coefficient of the cancer cells. Lastly, the obtained regression equation accounted for 99.14% of the prediction.

CONCLUSION

The model can simulate a cancer treatment process and predict the results of other radiation protocols.

Numerical simulation of normal and cancer cells' populations with fractional derivative under radiotherapy

Background This paper presents a numerical simulation of normal and cancer cells' population dynamics during radiotherapy. The model used for the simulation was the improved cancer treatment model with radiotherapy. The model simulated the population changes during a fractionated cancer treatment process. The results gave the final populations of the cells, which provided the final volumes of the tumor and normal cells. Method The improved model was obtained by integrating the previous cancer treatment model with the Caputo fractional derivative. In addition, the cells' population decay due to radiation was accounted for by coupling the linear-quadratic model into the improved model. The simulation of the treatment process was done with numerical variables, numerical parameters, and radiation parameters. The numerical variables include the populations of the cells and the time of treatment. The numerical parameters were the model factors which included the proliferation rates of cells, competition coefficients of cells, and perturbation constant for normal cells. The radiation parameters were clinical data based on the treatment procedure. The numerical parameters were obtained from the previous literature while the numerical variables and radiation parameters, which were clinical data, were obtained from reported data of four cancer patients treated with radiotherapy. The four cancer patients had tumor volumes of 28.4 cm³, 18.8 cm³, 30.6 cm³, and 12.6 cm³ and were treated with different treatment plans and a fractionated dose of 1.8 Gy each. The initial populations of cells were obtained by using the tumor volumes. The computer simulations were done with MATLAB. Results The final volumes of the tumors, from the results of the simulations, were 5.67 cm³, 4.36 cm³, 5.74 cm³, and 6.15 cm³ while the normal cells' volumes were 28.17 cm³, 18.68 cm³, 30.34 cm³, and 12.54 cm³. The powers of the derivatives were 0.16774, 0.16557, 0.16835, and 0.16. A variance-based sensitivity analysis was done to corroborate the model with the clinical data. The result showed that the most sensitive factors were the power of the derivative and the cancer cells' proliferation rate. Conclusion The model provided information concerning the status of treatments and can also predict outcomes of other treatment plans.

High-dose-rate vs. low-dose-rate interstitial brachytherapy boost for anal canal cancers

PURPOSE

The purpose of this study was to analyze and compare clinical outcomes of low-dose-rate (LDR) and high-dose-rate (HDR) interstitial brachytherapy boost (ISBT) after EBRT or radio chemotherapy for the treatment of anal canal cancers.

METHODS AND MATERIALS

One hundred patients with anal canal cancers were treated at our institution by ISBT [LDR (n = 50); HDR (n = 50)]. Chronic toxicity rates, local control, disease-free survival, overall survival, and colostomy-free survival of the two different dose-rate brachytherapy modalities were analyzed and compared.

RESULTS

With a median followup of 42.2 months (95% CI, [34.5-48.8]), 9 (9% [4.8-16.2%]) local recurrences were observed, 4 (8% [3.2-18.8%]) in LDR vs. 5 (10% [4.4-21.4%]) in HDR group (odds ratio [OR] = 1.28 [0.32-5.07], p = 0.73). The 5-year rate of local control for the entire population was 90% [81-95%], 93% [79-98%] vs. 86% [69-94%] for LDR and HDR, respectively (p = 0.38). The 5-year disease-free survival rate for all patients was 82% [71-90%], 88% [73-95%] vs. 72% [44-88%] for LDR and HDR, respectively (p = 0.21). The 5-year overall survival rate for global population was 94% [84-98%], with no significant differences between LDR (97% [79-100%]) and HDR (93% [80-98%]) (p = 0.27). The 5-year colostomy-free survival rate was 92% [83-96%], respectively, 95% [83-99%] vs. 86% [69-94%] for LDR and HDR (p = 0.21). Significant differences were found in terms of chronic toxicity rates, with 28 (56% [42.3-68.8%]) patients concerned in low-dose-rate brachytherapy vs. 17 (34% [22.4-47.9%]) in high-dose-rate brachytherapy (OR = 0.40 [0.18-0.91], p = 0.03).

CONCLUSIONS

Local recurrence rates were comparable between both groups; HDR brachytherapy seem to have a better toxicity profile. Our data confirmed the finding that HDR can be used to safely administer ISBT without increasing chronic toxicity.

Comparison of hyper- and hypofractionated radiation schemes with IMRT technique in small cell lung cancer: Clinical outcomes and the introduction of extended LQ and TCP models

PURPOSE

To evaluate the outcomes of 45 Gy/15 fractions/once-daily and 45 Gy/30 fractions/twice-daily radiation schemes utilizing intensity-modulated radiation therapy (IMRT) in extensive stage small cell lung cancer (SCLC), and to build up a new radiobiological model for tumor control probability (TCP) considering multiple biological effects.

METHODS

Fifty-eight consecutive patients diagnosed with extensive stage SCLC, treated with chemotherapy and chest irradiation, were retrospectively reviewed. Thirty-seven received hyperfractionated IMRT (Hyper-IMRT, 45 Gy/30 fractions/twice-daily) and 21 received hypofractionated IMRT (Hypo-IMRT, 45 Gy/15 fractions/once-daily). Local progression-free survival (LPFS) and overall survival (OS) were calculated and compared. An extended linear-quadratic (LQ) model, LQRG, incorporating cell repair, redistribution, reoxygenation, regrowth and Gompertzian tumor growth was created based on the clinical data. The TCP model was reformulated to predict LPFS. The classical LQ and TCP models were compared with the new models. Akaike information criterion (AIC) was used to assess the quality of the models.

RESULTS

The 2-year LPFS (34.1% vs 27.9%, p = 0.44) and OS (76.9% vs 76.9%, p = 0.26) were similar between Hyper- and Hypo-IMRT patients. According to the LQRG model, the α/β calculated was 9.2 (95% confidence interval: 8.7-9.9) Gy after optimization. The average absolute and relative fitting errors for LPFS were 9.1% and 18.7% for Hyper-IMRT, and 8.8% and 16.2% for Hypo-IMRT of the new TCP model, compared with 29.1% and 62.3% for Hyper-IMRT, and 30.7% and 65.3% for Hypo-IMRT of the classical model.

CONCLUSIONS

Hypo- and Hyper-IMRT resulted in comparable local control in the chest irradiation of extensive stage SCLC. The LQRG model has better performance in predicting the TCP (or LPFS) of the two schemes.

Appropriate radiation dose for symptomatic relief and local control in patients with adult T cell leukemia/lymphoma

Adult T-cell leukemia/lymphoma (ATL) is an aggressive peripheral T-cell neoplasm that occurs only in patients with human T-cell leukemia virus type 1. No large study or randomized trial investigating radiotherapy (RT) for ATL has been performed. We retrospectively reviewed 55 courses of RT for 41 consecutive patients with ATL who underwent RT between 2000 and 2016 at our institutions. The results showed that RT for local ATL lesions can achieve symptomatic improvement in 92% of cases. Local remission, either complete remission (CR) or partial response (PR), was achieved in 100% of the patients (CR: 89%, PR: 11%) with ≥40 Gy irradiation. CR or PR was achieved in 71% (CR: 29%, PR: 43%) with 30-39 Gy and in 73% (CR: 6.7%, PR: 67%) with ≤29 Gy irradiation. The mean total radiation dose in the CR and PR groups differed significantly (38 vs 25 Gy, P = 0.0002). The maximum acute toxicity was Grade 0-2 in all patients, except for one patient experienced Grade 3 radiation dermatitis. In-field relapses occurred in 36% of patients, and the frequency of in-field relapses was 11%, 30% and 71% among those who achieved CR, PR and SD, respectively. All 9 patients who received total skin irradiation experienced cutaneous relapses, with a median of 63 days (range, 7-210 days). Almost all (39 of 41) patients with ATL experienced out-of-field progression after RT. In conclusion, RT was confirmed to be effective and safe for palliative treatment of local ATL lesions.

Re-irradiation volumetric modulated arc therapy optimization based on cumulative biologically effective dose objectives

The objective of this note is to introduce a clinical tool that generates ideal base plan dose distributions to enable re‐irradiation volumetric modulated arc therapy (VMAT) optimization based on cumulative biological effective dose objectives for specific organs at risk (OARs). The tool is demonstrated with a lung cancer case that required re‐irradiation at our clinic. First, previous treatment dose is deformed onto the retreatment computed tomography (CT) using commercial software. Then, the in‐house Matlab tool alters the deformed previous dose using radiobiological concepts on a voxel‐by‐voxel manner to generate an ideal base plan dose distribution. Ideal base plans that were generated using the in‐house Matlab tool were compatible with the Varian Eclipse™ treatment planning system. The tool enabled optimization of VMAT re‐irradiation plans using cumulative dose limits for OARs and all OAR cumulative dose objectives were met on the first optimization for the recurrent lung cancer case tested.

Accelerated Intensity-Modulated External Radiotherapy as a New Technical Approach to Treat the Index Quadrant after Conserving Surgery in Early Breast Cancer: A Preliminary Study

Purpose

We present a novel technical approach to treat the index quadrant after conserving surgery in patients with early breast cancer and study its clinical feasibility.

Methods and material

Patients selected for the study, after breast conserving surgery with histologically verified breast carcinoma, signed a full informed consent to intensity-modulated external radiotherapy of the partial breast. Treatment was performed with the 6 MV beam from one of the Elekta Precise LINAC units installed in our Department. The prescribed dose was 30 Gy in 5 fractions in 10 days.

Results

Acute toxicity was minimal. No skin changes were noted during treatment or during the first 6 months after radiotherapy treatment.

Conclusions

Accelerated partial breast irradiation using intensity-modulated external radiotherapy is technically feasible. We think the approach will give good results in terms of local control, toxicity and quality of life, at the same time sparing resources for the patient and health care system.

A feasibility study on adaptive 18F-FDG-PET-guided radiotherapy for recurrent and second primary head and neck cancer in the previously irradiated territory

PURPOSE

To evaluate feasibility, disease control, survival, and toxicity after adaptive ¹⁸F-fluorodeoxyglucose (FDG) positron emisson tomography (PET) guided radiotherapy in patients with recurrent and second primary head and neck squamous cell carcinoma.

METHODS

A prospective trial investigated the feasibility of adaptive intensity modulated radiotherapy (IMRT) ± concomitant cetuximab in 10 patients. The primary endpoint was achieving a 2-year survival free of grade >3 toxicity in ≥30% of patients. Three treatment plans based on 3 PET/CT scans were consecutively delivered in 6 weeks. The range of dose painting was 66.0-85.0 Gy in the dose-painted tumoral volumes in 30 fractions.

RESULTS

Two-year locoregional and distant control rates were 38 and 76%, respectively. Overall and disease-free survival at 2 years was 20%. No grade 4 or 5 acute toxicity was observed in any of the patients, except for arterial mucosal hemorrhage in 1 patient. Three months after radiotherapy, grade 4 dysphagia and mucosal wound healing problems were observed in 1/7 and 1/6 of patients, respectively. Grade 5 toxicity (fatal bleeding) was seen in 2 patients, at 3.8 and 4.1 months of follow-up. Data on 2‑year toxicity could only be assessed in 1 of the 2 surviving patients, in whom grade 4 mucosal wound healing problems were observed; no other grade >3 toxicity was observed. In this respect, a 30% 2‑year survival free of grade >3 toxicity will not be achieved.

CONCLUSIONS

Adaptive PET-guided reirradiation is feasible. However, due to slow accrual and treatment results that seemed inconsistent with achieving the primary endpoint, the trial was stopped early.

Late mucosal ulcers in dose-escalated adaptive dose-painting treatments for head-and-neck cancer

BACKGROUND

To identify predictive factors for the development of late grade 4 mucosal ulcers in adaptive dose-escalated treatments for head-and-neck cancer.

MATERIAL AND METHODS

Patient data of four dose-escalated three-phase adaptive dose-painting by numbers (DPBN) clinical trials were analyzed in this study. Correlations between the development of late grade 4 ulcers and factors related with the treatment, disease characteristics and the patient were investigated. Dosimetrical thresholds were searched among the highest doses received by 1.75 cm³ (D_1.75cc) of the primary gross tumor volume (GTV_T) and the corresponding normalized isoeffective dose (NID2_1.75cc, with a reference dose of 2Gy/fraction and α/β of 3 Gy).

RESULTS

From 39 studied patients, nine developed late grade 4 mucosal ulcers. The continuation to either smoke or drink alcohol after therapy was the factor that showed a strong (eight out of nine patients) association with the occurrence of grade 4 ulcers. Six of the patients who continued to smoke or/and drink had D_1.75cc and NID2_1.75cc above 84 Gy and 95.5 Gy, respectively. Seven of the patients with grade 4 had the dose levels above these thresholds, but even if the D_1.75cc threshold was significant in the prediction of late grade 4 ulcers, it could not be considered as the only contributing factor.

CONCLUSIONS

The search for patterns provided strong reasons to apply a dosimetrical threshold for the peak-dose volume of 1.75 cm³ as a preventive measure for late grade 4 mucosal ulcers. Also, patients that continue to smoke or drink alcohol after therapy have increased risk to develop late mucosal ulcers.

Development of 3D biological effective dose distribution software program

Purpose

To develop a software program to convert physical dose distribution into biological effective dose (BED).

Methods

The MATLAB-based BED distribution software program was designed to import the radiotherapy treatment plan from the computer treatment planning system and to convert the physical dose distribution into the BED distribution. The BED calculation was based on the linear-quadratic-linear model (LQ-L model). Besides radiobiological parameters, other specific data could be fed in through the panel. The accuracy of the program was verified by comparing the BED distribution with manual calculation.

Results

This software program was able to import the radiotherapy treatment plans and pull out pixel-wised physical dose for BED calculation, and display the isoBED lines on the computed tomographic (CT) image. The verification of BED dose distribution was performed in both phantom and clinical cases. It revealed that there were no differences between the program and manual BED calculations.

Conclusion

It is feasible and practical to use this in-house BED distribution software program in clinical practices and research work. However, it should be used with caution as the validity of the program depends on the accuracy of the published biological parameters.

Clinical analysis of the approximate, 3-dimensional, biological effective dose equation in multiphase treatment plans

A multiphase, approximate biological effective dose (BED_A) equation was introduced because most treatment planning systems (TPS) are incapable of calculating the true BED (BED_T). This work investigates the accuracy and precision of the multiphase BED_A relative to the BED_T in clinical cases. Ten patients with head and neck cancer and 10 patients with prostate cancer were studied using their treatment plans from Pinnacle³ 9.2 (Philips Medical, Fitchburg, WI). The organs at risk (OARs) that were studied are the normal brain, left and right optic nerves, optic chiasm, spinal cord, brainstem, bladder, and rectum. BED_A and BED_T distributions were calculated using MATLAB 2010b (MathWorks, Natick, MA) and analyzed on a voxel basis for percent error, percent error volume histograms (PEVHs), Pearson correlation coefficient, and Bland-Altman analysis. The maximum BED values that were calculated using the BED_A and BED_T methods were also analyzed. BED_A was found to always underestimate BED_T. The accuracy and precision of BED_A distributions varied between the organs: for optic chiasm and brainstem, 50% of the patients had an overall BED_A percent error of <1%; for left and right optic nerves, rectum, and bladder, 60% to 70% of the patients had an overall BED_A percent error of <1%; and for normal brain and spinal cord, 80% of the patients had an overall BED_A percent error of <1%. BED_A distributions had maximum errors ranging from 2% to 11%, with the 11% error occurring for bladder. BED_A produced much more accurate maximum BED values with adjacent organs such as normal brain, bladder, and rectum. This study has shown that BED_A can calculate BED distributions with acceptable accuracy under certain circumstances. However, its consistency and accuracy strongly depend on the dose distributions of the different treatment phases. One should be cautious when using BED_A.

Importance of deformable image registration and biological dose summation in planning of radiotherapy retreatments

The purpose of this study is to evaluate the effect of nonrigid and fractionation-corrected dose summation on total doses in radiotherapy and to demonstrate the benefits of such dose summation on clinical decision-making for planning of retreatments. Dose summation of organs at risk (OARs) was investigated for 3 clinical cases with need of retreatment to the same site: head and neck, brain, and mediastinum. Three different summation methods over old and new radiotherapy treatment plans are presented and compared: (1) rigid raw sum with rigid registration of the planning images and direct dose summing; (2) deformable raw sum with deformable image registration and direct dose summing; and (3) deformable biological sum with deformable registration and takes into account the dose per fraction in biological manner in certain critical organs. In 2 cases, a user-defined dose downscaling is applied to take into account the time between the treatments and the healing from the radiation-induced effects. Of the 3 summation methods presented, the deformable biological sum was considered to offer the most biologically plausible account of the treatment. There were remarkable differences between near-maximum doses (D0.1cc) and dose-volume histogram (DVH) curves for OARs between different summation methods. The differences between deformable raw sum and rigid raw sum D0.1cc doses are in the range from -8 Gy to 2 Gy. Similarly, the deviation was from -14 Gy to 5 Gy for the deformable biological sum compared with the rigid raw sum. These differences come from incorrect summation of doses in the rigid raw sum case, and from the dose per fraction effect in biological summation. We conclude that computing the 3-dimensional deformable biological summation could be a valuable tool for treating patients with complex retreatments. It has the potential to assist the oncologist in refining plans for maximally curative doses while respecting appropriate tissue tolerances.

Quality assessment of delineation and dose planning of early breast cancer patients included in the randomized Skagen Trial 1

BACKGROUND AND PURPOSE

To report on a Quality assessment (QA) of Skagen Trial 1, exploring hypofractionation for breast cancer patients with indication for regional nodal radiotherapy.

MATERIAL AND METHODS

Deviations from protocol regarding target volume delineations and dose parameters (Dmin, Dmax, D98%, D95% and D2%) from randomly selected dose plans were assessed. Target volume delineation according to ESTRO guidelines was obtained through atlas based automated segmentation and centrally approved as gold standard (GS). Dice similarity scores (DSC) with original delineations were measured. Dose parameters measured in the two delineations were reported to assess their dosimetric outcome.

RESULTS

Assessment included 88 plans from 12 centres in 4 countries. DSC showed high agreement in contouring, 99% and 96% of the patients had a complete delineation of target volumes and organs at risk. No deviations in the dosimetric outcome were found in 76% of the patients, 82% and 95% of the patients had successful coverage of breast/chestwall and CTVn_L2-4-interpectoral. Dosimetric outcomes of original delineation and GS were comparable.

CONCLUSIONS

QA showed high protocol compliance and adequate dose coverage in most patients. Inter-observer variability in contouring was low. Dose parameters were in harmony with protocol regardless original or GS segmentation.

A graphical user interface (GUI) toolkit for the calculation of three-dimensional (3D) multi-phase biological effective dose (BED) distributions including statistical analyses

A toolkit has been developed for calculating the 3-dimensional biological effective dose (BED) distributions in multi-phase, external beam radiotherapy treatments such as those applied in liver stereotactic body radiation therapy (SBRT) and in multi-prescription treatments. This toolkit also provides a wide range of statistical results related to dose and BED distributions. MATLAB 2010a, version 7.10 was used to create this GUI toolkit. The input data consist of the dose distribution matrices, organ contour coordinates, and treatment planning parameters from the treatment planning system (TPS). The toolkit has the capability of calculating the multi-phase BED distributions using different formulas (denoted as true and approximate). Following the calculations of the BED distributions, the dose and BED distributions can be viewed in different projections (e.g. coronal, sagittal and transverse). The different elements of this toolkit are presented and the important steps for the execution of its calculations are illustrated. The toolkit is applied on brain, head & neck and prostate cancer patients, who received primary and boost phases in order to demonstrate its capability in calculating BED distributions, as well as measuring the inaccuracy and imprecision of the approximate BED distributions. Finally, the clinical situations in which the use of the present toolkit would have a significant clinical impact are indicated.

Effect of Irradiation on Tissue Penetration Depth of Doxorubicin after Pressurized Intra-Peritoneal Aerosol Chemotherapy (PIPAC) in a Novel Ex-Vivo Model

Background: This study was performed to assess the impact of irradiation on the tissue penetration depth of doxorubicin delivered during Pressurized Intra-Peritoneal Aerosol Chemotherapy (PIPAC).

Methods: Fresh post mortem swine peritoneum was cut into 10 proportional sections. Except for 2 control samples, all received irradiation with 1, 2, 7 and 14 Gy, respectively. Four samples received PIPAC 15 minutes after irradiation and 4 other after 24 hours. Doxorubicin was aerosolized in an ex-vivo PIPAC model at 12 mmHg/36°C. In-tissue doxorubicin penetration was measured using fluorescence microscopy on frozen thin sections.

Results: Doxorubicin penetration after PIPAC (15 minutes after irradiation) was 476 ± 74 µm for the control sample, 450 ± 45µm after 1 Gy (p > 0.05), 438 ± 29 µm after 2 Gy (p > 0.05), 396 ± 32 µm after 7 Gy (p = 0.005) and 284 ± 57 after 14 Gy irradiation (p < 0.001). The doxorubicin penetration after PIPAC (24 hours after irradiation) was 428 ± 77 µm for the control sample, 393 ± 41 µm after 1 Gy (p > 0.05), 379 ± 56 µm after 2 Gy (p > 0.05), 352 ± 53 µm after 7 Gy (p = 0.008) and 345 ± 53 after 14 Gy irradiation (p = 0.001).

Conclusions: Higher (fractional) radiation dose might reduce the tissue penetration depth of doxorubicin in our ex-vivo model. However, irradiation with lower (fractional) radiation dose does not affect the tissue penetration negatively. Further studies are warranted to investigate if irradiation can be used safely as chemopotenting agent for patients with peritoneal metastases treated with PIPAC.

Comparison of breast sequential and simultaneous integrated boost using the biologically effective dose volume histogram (BEDVH)

Purpose

A method is presented to radiobiologically compare sequential (SEQ) and simultaneously integrated boost (SIB) breast radiotherapy.

Methods

The method is based on identically prescribed biologically effective dose (iso-BED) which was achieved by different prescribed doses due to different fractionation schemes. It is performed by converting the calculated three-dimensional dose distribution to the corresponding BED distribution taking into consideration the different number of fractions for generic α/β ratios. A cumulative BED volume histogram (BEDVH) is then derived from the BED distribution and is compared for the two delivery schemes. Ten breast cancer patients (4 right-sided and 6 left-sided) were investigated. Two tangential intensity modulated whole breast beams with two other oblique (with different gantry angles) beams for the boost volume were used. The boost and the breast target volumes with either α/β = 10 or 3 Gy, and ipsi-lateral and contra-lateral lungs, heart, and contra-lateral breast as organs at risk (OARs) with α/β = 3 Gy were compared.

Results

Based on the BEDVH comparisons, the use of SIB reduced the biological breast mean dose by about 3 %, the ipsi-lateral lung and heart by about 10 %, and contra-lateral breast and lung by about 7 %.

Conclusion

BED based comparisons should always be used in comparing plans that have different fraction sizes. SIB schemes are dosimetrically more advantageous than SEQ in breast target volume and OARs for equal prescribed BEDs for breast and boost.

The Role of Hypofractionated Radiation Therapy with Photons, Protons, and Heavy Ions for Treating Extracranial Lesions

Traditionally, the ability to deliver large doses of ionizing radiation to a tumor has been limited by radiation-induced toxicity to normal surrounding tissues. This was the initial impetus for the development of conventionally fractionated radiation therapy, where large volumes of healthy tissue received radiation and were allowed the time to repair the radiation damage. However, advances in radiation delivery techniques and image guidance have allowed for more ablative doses of radiation to be delivered in a very accurate, conformal, and safe manner with shortened fractionation schemes. Hypofractionated regimens with photons have already transformed how certain tumor types are treated with radiation therapy. Additionally, hypofractionation is able to deliver a complete course of ablative radiation therapy over a shorter period of time compared to conventional fractionation regimens making treatment more convenient to the patient and potentially more cost-effective. Recently, there has been an increased interest in proton therapy because of the potential further improvement in dose distributions achievable due to their unique physical characteristics. Furthermore, with heavier ions the dose conformality is increased and, in addition, there is potentially a higher biological effectiveness compared to protons and photons. Due to the properties mentioned above, charged particle therapy has already become an attractive modality to further investigate the role of hypofractionation in the treatment of various tumors. This review will discuss the rationale and evolution of hypofractionated radiation therapy, the reported clinical success with initially photon and then charged particle modalities, and further potential implementation into treatment regimens going forward.

A patient-specific computational model of hypoxia-modulated radiation resistance in glioblastoma using 18F-FMISO-PET

Glioblastoma multiforme (GBM) is a highly invasive primary brain tumour that has poor prognosis despite aggressive treatment. A hallmark of these tumours is diffuse invasion into the surrounding brain, necessitating a multi-modal treatment approach, including surgery, radiation and chemotherapy. We have previously demonstrated the ability of our model to predict radiographic response immediately following radiation therapy in individual GBM patients using a simplified geometry of the brain and theoretical radiation dose. Using only two pre-treatment magnetic resonance imaging scans, we calculate net rates of proliferation and invasion as well as radiation sensitivity for a patient's disease. Here, we present the application of our clinically targeted modelling approach to a single glioblastoma patient as a demonstration of our method. We apply our model in the full three-dimensional architecture of the brain to quantify the effects of regional resistance to radiation owing to hypoxia in vivo determined by [(18)F]-fluoromisonidazole positron emission tomography (FMISO-PET) and the patient-specific three-dimensional radiation treatment plan. Incorporation of hypoxia into our model with FMISO-PET increases the model-data agreement by an order of magnitude. This improvement was robust to our definition of hypoxia or the degree of radiation resistance quantified with the FMISO-PET image and our computational model, respectively. This work demonstrates a useful application of patient-specific modelling in personalized medicine and how mathematical modelling has the potential to unify multi-modality imaging and radiation treatment planning.

Practical aspects and uncertainty analysis of biological effective dose (BED) regarding its three-dimensional calculation in multiphase radiotherapy treatment plans

PURPOSE

There is a growing interest in the radiation oncology community to use the biological effective dose (BED) rather than the physical dose (PD) in treatment plan evaluation and optimization due to its stronger correlation with radiobiological effects. Radiotherapy patients may receive treatments involving a single only phase or multiple phases (e.g., primary and boost). Since most treatment planning systems cannot calculate the analytical BED distribution in multiphase treatments, an approximate multiphase BED expression, which is based on the total physical dose distribution, has been used. The purpose of this paper is to reveal the mathematical properties of the approximate BED formulation, relative to the true BED.

METHODS

The mathematical properties of the approximate multiphase BED equation are analyzed and evaluated. In order to better understand the accuracy of the approximate multiphase BED equation, the true multiphase BED equation was derived and the mathematical differences between the true and approximate multiphase BED equations were determined. The magnitude of its inaccuracies under common clinical circumstances was also studied. All calculations were performed on a voxel-by-voxel basis using the three-dimensional dose matrices.

RESULTS

Results showed that the approximate multiphase BED equation is accurate only when the dose-per-fractions (DPFs) in both the first and second phases are equal, which occur when the dose distribution does not significantly change between the phases. In the case of heterogeneous dose distributions, which significantly vary between the phases, there are fewer occurrences of equal DPFs and hence the inaccuracy of the approximate multiphase BED is greater. These characteristics are usually seen in the dose distributions being delivered to organs at risk rather than to targets.

CONCLUSIONS

The finding of this study indicates that the true multiphase BED equation should be implemented in the treatment planning systems due to the inconsistent accuracy of the approximate multiphase BED equation in most of the clinical situations.

Bringing the heavy: carbon ion therapy in the radiobiological and clinical context

Radiotherapy for the treatment of cancer is undergoing an evolution, shifting to the use of heavier ion species. For a plethora of malignancies, current radiotherapy using photons or protons yields marginal benefits in local control and survival. One hypothesis is that these malignancies have acquired, or are inherently radioresistant to low LET radiation. In the last decade, carbon ion radiotherapy facilities have slowly been constructed in Europe and Asia, demonstrating favorable results for many of the malignancies that do poorly with conventional radiotherapy. However, from a radiobiological perspective, much of how this modality works in overcoming radioresistance, and extending local control and survival are not yet fully understood. In this review, we will explain from a radiobiological perspective how carbon ion radiotherapy can overcome the classical and recently postulated contributors of radioresistance (α/β ratio, hypoxia, cell proliferation, the tumor microenvironment and metabolism, and cancer stem cells). Furthermore, we will make recommendations on the important factors to consider, such as anatomical location, in the future design and implementation of clinical trials. With the existing data available we believe that the expansion of carbon ion facilities into the United States is warranted.

Estimating dose painting effects in radiotherapy: a mathematical model

Tumor heterogeneity is widely considered to be a determinant factor in tumor progression and in particular in its recurrence after therapy. Unfortunately, current medical techniques are unable to deduce clinically relevant information about tumor heterogeneity by means of non-invasive methods. As a consequence, when radiotherapy is used as a treatment of choice, radiation dosimetries are prescribed under the assumption that the malignancy targeted is of a homogeneous nature. In this work we discuss the effects of different radiation dose distributions on heterogeneous tumors by means of an individual cell-based model. To that end, a case is considered where two tumor cell phenotypes are present, which we assume to strongly differ in their respective cell cycle duration and radiosensitivity properties. We show herein that, as a result of such differences, the spatial distribution of the corresponding phenotypes, whence the resulting tumor heterogeneity can be predicted as growth proceeds. In particular, we show that if we start from a situation where a majority of ordinary cancer cells (CCs) and a minority of cancer stem cells (CSCs) are randomly distributed, and we assume that the length of CSC cycle is significantly longer than that of CCs, then CSCs become concentrated at an inner region as tumor grows. As a consequence we obtain that if CSCs are assumed to be more resistant to radiation than CCs, heterogeneous dosimetries can be selected to enhance tumor control by boosting radiation in the region occupied by the more radioresistant tumor cell phenotype. It is also shown that, when compared with homogeneous dose distributions as those being currently delivered in clinical practice, such heterogeneous radiation dosimetries fare always better than their homogeneous counterparts. Finally, limitations to our assumptions and their resulting clinical implications will be discussed.

Critical dose and toxicity index of organs at risk in radiotherapy: analyzing the calculated effects of modified dose fractionation in non-small cell lung cancer

To increase the efficacy of radiotherapy for non-small cell lung cancer (NSCLC), many schemes of dose fractionation were assessed by a new "toxicity index" (I), which allows one to choose the fractionation schedules that produce less toxic treatments. Thirty-two patients affected by non resectable NSCLC were treated by standard 3-dimensional conformal radiotherapy (3DCRT) with a strategy of limited treated volume. Computed tomography datasets were employed to re plan by simultaneous integrated boost intensity-modulated radiotherapy (IMRT). The dose distributions from plans were used to test various schemes of dose fractionation, in 3DCRT as well as in IMRT, by transforming the dose-volume histogram (DVH) into a biological equivalent DVH (BDVH) and by varying the overall treatment time. The BDVHs were obtained through the toxicity index, which was defined for each of the organs at risk (OAR) by a linear quadratic model keeping an equivalent radiobiological effect on the target volume. The less toxic fractionation consisted in a severe/moderate hyper fractionation for the volume including the primary tumor and lymph nodes, followed by a hypofractionation for the reduced volume of the primary tumor. The 3DCRT and IMRT resulted, respectively, in 4.7% and 4.3% of dose sparing for the spinal cord, without significant changes for the combined-lungs toxicity (p < 0.001). Schedules with reduced overall treatment time (accelerated fractionations) led to a 12.5% dose sparing for the spinal cord (7.5% in IMRT), 8.3% dose sparing for V20 in the combined lungs (5.5% in IMRT), and also significant dose sparing for all the other OARs (p < 0.001). The toxicity index allows to choose fractionation schedules with reduced toxicity for all the OARs and equivalent radiobiological effect for the tumor in 3DCRT, as well as in IMRT, treatments of NSCLC.

Analysis of fractionation correction methodologies for multiple phase treatment plans in radiation therapy

PURPOSE

Radiation therapy is often delivered by multiple sequential treatment plans. For an accurate radiobiological evaluation of the overall treatment, fractionation corrections to each dose distribution must be applied before summing the three-dimensional dose matrix of each plan since the simpler approach of performing the fractionation correction to the total dose-volume histograms, obtained by the arithmetical sum of the different plans, becomes inaccurate for more heterogeneous dose patterns. In this study, the differences between these two fractionation correction methods, named here as exact (corrected before) and approximate (after summation), respectively, are assessed for different cancer types.

METHODS

Prostate, breast, and head and neck (HN) tumor patients were selected to quantify the differences between two fractionation correction methods (the exact vs the approximate). For each cancer type, two different treatment plans were developed using uniform (CRT) and intensity modulated beams (IMRT), respectively. The responses of the target and normal tissue were calculated using the Poisson linear-quadratic-time model and the relative seriality model, respectively. All treatments were radiobiologically evaluated and compared using the complication-free tumor control probability (P+), the biologically effective uniform dose (D) together with common dosimetric criteria.

RESULTS

For the prostate cancer patient, an underestimation of around 14%-15% in P+ was obtained when the fractionation correction was applied after summation compared to the exact approach due to significant biological and dosimetric variations obtained between the two fractionation correction methods in the involved lymph nodes. For the breast cancer patient, an underestimation of around 3%-4% in the maximum dose in the heart was obtained. Despite the dosimetric differences in this organ, no significant variations were obtained in treatment outcome. For the HN tumor patient, an underestimation of about 5% in treatment outcome was obtained for the CRT plan as a result of an underestimation of the planning target volume control probability by about 10%. An underestimation of about 6% in the complication probability of the right parotid was also obtained. For all the other organs at risk, dosimetric differences of up to 4% were obtained but with no significant impact in the expected clinical outcome. However, for the IMRT plan, an overestimation in P+ of 4.3% was obtained mainly due to an underestimation of the complication probability of the left and right parotids (2.9% and 5.8%, respectively).

CONCLUSIONS

The use of the exact fractionation correction method, which is applying fractionation correction on the separate dose distributions of a multiple phase treatment before their summation was found to have a significant expected clinical impact. For regions of interest that are irradiated with very heterogeneous dose distributions and significantly different doses per fraction in the different treatment phases, the exact fractionation correction method needs to be applied since a significant underestimation of the true patient outcome can be introduced otherwise.

Radiation dose responses for chemoradiation therapy of pancreatic cancer: an analysis of compiled clinical data using biophysical models

PURPOSE

We analyzed recent clinical data obtained from chemoradiation of unresectable, locally advanced pancreatic cancer (LAPC) in order to examine possible benefits from radiation therapy dose escalation.

METHODS AND MATERIALS

A modified linear quadratic model was used to fit clinical tumor response and survival data of chemoradiation treatments for LAPC reported from 20 institutions. Biophysical radiosensitivity parameters were extracted from the fits.

RESULTS

Examination of the clinical data demonstrated an enhancement in tumor response with higher irradiation dose, an important clinical result for palliation and quality of life. Little indication of improvement in 1-year survival with increased radiation dose was observed. Possible dose escalation schemes are proposed based on calculations of the biologically effective dose required for a 50% tumor response rate.

CONCLUSIONS

Based on the evaluation of tumor response data, the escalation of radiation dose presents potential clinical benefits which when combined with normal tissue complication analyses may result in improved treatment outcome for locally advanced pancreatic cancer patients.

Mesenchymal stem cells show radioresistance in vivo

Irradiation impacts on the viability and differentiation capacity of tissue-borne mesenchymal stem cells (MSC), which play a pivotal role in bone regeneration. As a consequence of radiotherapy, bones may develop osteoradionecrosis. When irradiating human bone-derived MSC in vitro with increasing doses, the cells’ self-renewal capabilities were greatly reduced. Mitotically stalled cells were still capable of differentiating into osteoblasts and pre-adipocytes. As a large animal model comparable to the clinical situation, pig mandibles were subjected to fractionized radiation of 2 χ 9 Gy within 1 week. This treatment mimics that of a standardized clinical treatment regimen of head and neck cancer patients irradiated 30 χ 2 Gy. In the pig model, fractures which had been irradiated, showed delayed osseous healing. When isolating MSC at different time points post-irradiation, no significant changes regarding proliferation capacity and osteogenic differentiation potential became apparent. Therefore, pig mandibles were irradiated with a single dose of either 9 or 18 Gy in vivo, and MSC were isolated immediately afterwards. No significant differences between the untreated and 9 Gy irradiated bone with respect to proliferation and osteogenic differentiation were unveiled. Yet, cells isolated from 18 Gy irradiated specimens exhibited a reduced osteogenic differentiation capacity, and during the first 2 weeks proliferation rates were greatly diminished. Thereafter, cells recovered and showed normal proliferation behaviour. These findings imply that MSC can effectively cope with irradiation up to high doses in vivo. This finding should thus be implemented in future therapeutic concepts to protect regenerating tissue from radiation consequences.

Providing a fast conversion of total dose to biological effective dose (BED) for hybrid seed brachytherapy

Optimization of permanent seed implant brachytherapy plans for treatment of prostate cancer should be based on biological effective dose (BED) distributions, since dose does not accurately represent biological effects between different types of sources. Currently, biological optimization for these plans is not feasible due to the amount of time necessary to calculate the BED distribution. This study provides a fast calculation method, based on the total dose, to calculate the BED distribution. Distributions of various numbers of hybrid seeds were used to calculate total dose distributions, as well as BED distributions. Hybrid seeds are a mixture of different isotopes (in this study ¹²⁵I and ¹⁰³Pd). Three ratios of hybrid seeds were investigated: 25/75, 50/50, and 75/25. The total dose and BED value from each voxel were coupled together to produce graphs of total dose vs. BED. Equations were then derived from these graphs. The study investigated four types of tissue: bladder, rectum, prostate, and other normal tissue. Equations were derived from the total dose – BED correspondence. Accuracy of conversion from total dose to BED was within 2 Gy; however, accuracy of conversion was found to be better for high total dose regions as compared to lower dose regions. The method introduced in this paper allows one to perform fast conversion of total dose to BED for brachytherapy using hybrid seeds, which makes the BED‐based plan optimization practical. The method defined here can be extended to other ratios, as well as other tissues that are affected by permanent seed implant brachytherapy (i.e., breast).

PACS number: 87.55.de

Parotid gland sparing with helical tomotherapy in head-and-neck cancer

PURPOSE

This study evaluated the ability of helical tomotherapy to spare the function of the parotid glands in patients with head-and-neck cancer by analyzing dose-volume histograms, salivary gland scintigraphy, and quality of life assessment.

METHODS AND MATERIALS

Data from 76 consecutive patients treated with helical tomotherapy (Hi-Art Tomotherapy) at the University Hospital Brussel were analyzed. During planning, priority was given to planning target volume (PTV) coverage: ≥ 95% of the dose must be delivered to ≥ 95% of the PTV. Elective nodal regions received 54 Gy (1.8 Gy/fraction). A dose of 70.5 Gy (2.35 Gy/fraction) was prescribed to the primary tumor and pathologic lymph nodes (simultaneous integrated boost scheme). Objective scoring of salivary excretion was performed by salivary gland scintigraphy. Subjective scoring of salivary gland function was evaluated by the European Organization for Research and Treatment of Cancer quality of life questionnaires Quality of Life Questionnaire-C30 (QLQ-C30) and Quality of Life Questionnaire-Head and Neck 35 (H&N35).

RESULTS

Analysis of dose-volume histograms (DVHs) showed excellent coverage of the PTV. The volume of PTV receiving 95% of the prescribed dose (V95%) was 99.4 (range, 96.3-99.9). DVH analysis of parotid gland showed a median value of the mean parotid dose of 32.1 Gy (range, 17.5-70.3 Gy). The median parotid volume receiving a dose <26 Gy was 51.2%. Quality of life evaluation demonstrated an initial deterioration of almost all scales and items in QLQ-C30 and QLQ-H&N35. Most items improved in time, and some reached baseline values 18 months after treatment.

CONCLUSION

DVH analysis, scintigraphic evaluation of parotid function, and quality of life assessment of our patient group showed that helical tomotherapy makes it possible to preserve parotid gland function without compromising disease control. We recommend mean parotid doses of <34 Gy and doses <26 Gy to a maximum 47% of the parotid volume as planning goals. Intensity-modulated radiotherapy should be considered as standard treatment in patients with head-and-neck cancer.

Range uncertainty in proton therapy due to variable biological effectiveness

Traditionally, dose in proton radiotherapy is prescribed as Gy(RBE) by scaling up the physical dose by 10%. The relative biological effectiveness (RBE) of protons is considered to vary with dose-averaged linear energy transfer (LET(d)), dose (d) and (α/β)(x). The increase of RBE with depth causes a shift of the falloff of the beam, i.e. a change of the beam range. The magnitude of this shift will depend on dose and (α/β)(x). The aim of this project was to quantify the dependence of the range shift on these parameters. Three double-scattered beams of different ranges incident on a computational phantom consisting of different regions of interest (ROIs) were used. Each ROI was assigned with (α/β)(x) values between 0.5 and 20 Gy. The distribution of LET(d) within each ROI was obtained from a Monte Carlo simulation. The LET(d) distribution depends on the beam energy and thus its nominal range. The RBE values within the ROIs were calculated for doses between 1 and 15 Gy using an in-house developed biophysical model. Dose-volume histograms of the RBE-weighted doses were extracted for each ROI for a 'fixed RBE' (RBE = 1.1) and a 'variable RBE' (RBE = f (d, α/β, LET(d))), and the percentage difference in range was obtained from the difference of the percentage volumes at the distal 80% of the dose. Range differences in normal tissue ((α/β)(x) = 3 Gy) of the order of 3-2 mm were obtained, respectively, for a shallow (physical range 4.8 cm) and a deep (physical range 12.8 cm) beam, when a dose of 1 Gy normalized to the mid-SOBP was delivered. As the dose increased to 15 Gy, the variable RBE decreases below 1.1 which induces ranges of about 1 mm shorter than those obtained with an RBE of 1.1. The shift in the range of an SOBP when comparing biological dose distributions obtained with a fixed or a variable RBE was quantified as a function of dose, (α/β)(x) and physical range (as a surrogate of the initial beam energy). The shift increases with the physical range but decreases with increasing dose or (α/β)(x). The results of our study allow a quantitative consideration of RBE-caused range uncertainties as a function of treatment site and dose in treatment planning.

ICRP publication 118: ICRP statement on tissue reactions and early and late effects of radiation in normal tissues and organs--threshold doses for tissue reactions in a radiation protection context

This report provides a review of early and late effects of radiation in normal tissues and organs with respect to radiation protection. It was instigated following a recommendation in Publication 103 (ICRP, 2007), and it provides updated estimates of 'practical' threshold doses for tissue injury defined at the level of 1% incidence. Estimates are given for morbidity and mortality endpoints in all organ systems following acute, fractionated, or chronic exposure. The organ systems comprise the haematopoietic, immune, reproductive, circulatory, respiratory, musculoskeletal, endocrine, and nervous systems; the digestive and urinary tracts; the skin; and the eye. Particular attention is paid to circulatory disease and cataracts because of recent evidence of higher incidences of injury than expected after lower doses; hence, threshold doses appear to be lower than previously considered. This is largely because of the increasing incidences with increasing times after exposure. In the context of protection, it is the threshold doses for very long follow-up times that are the most relevant for workers and the public; for example, the atomic bomb survivors with 40-50years of follow-up. Radiotherapy data generally apply for shorter follow-up times because of competing causes of death in cancer patients, and hence the risks of radiation-induced circulatory disease at those earlier times are lower. A variety of biological response modifiers have been used to help reduce late reactions in many tissues. These include antioxidants, radical scavengers, inhibitors of apoptosis, anti-inflammatory drugs, angiotensin-converting enzyme inhibitors, growth factors, and cytokines. In many cases, these give dose modification factors of 1.1-1.2, and in a few cases 1.5-2, indicating the potential for increasing threshold doses in known exposure cases. In contrast, there are agents that enhance radiation responses, notably other cytotoxic agents such as antimetabolites, alkylating agents, anti-angiogenic drugs, and antibiotics, as well as genetic and comorbidity factors. Most tissues show a sparing effect of dose fractionation, so that total doses for a given endpoint are higher if the dose is fractionated rather than when given as a single dose. However, for reactions manifesting very late after low total doses, particularly for cataracts and circulatory disease, it appears that the rate of dose delivery does not modify the low incidence. This implies that the injury in these cases and at these low dose levels is caused by single-hit irreparable-type events. For these two tissues, a threshold dose of 0.5Gy is proposed herein for practical purposes, irrespective of the rate of dose delivery, and future studies may elucidate this judgement further.

First results of the randomised UK FAST Trial of radiotherapy hypofractionation for treatment of early breast cancer (CRUKE/04/015)

BACKGROUND AND PURPOSE

Randomised trials testing 15- or 16-fraction regimens of adjuvant radiotherapy in women with early breast cancer have reported favourable outcomes compared with standard fractionation. To evaluate hypofractionation further, two 5-fraction schedules delivering 1 fraction per week have been tested against a 25-fraction regimen.

MATERIALS AND METHODS

Women aged ⩾50years with node negative early breast cancer were randomly assigned after microscopic complete tumour resection to 50Gy in 25 fractions versus 28.5 or 30Gy in 5 once-weekly fractions of 5.7 or 6.0Gy, respectively, to the whole breast. The primary endpoint was 2-year change in photographic breast appearance.

RESULTS

Nine hundred and fifteen women were recruited from 2004 to 2007. Seven hundred and twenty-nine patients had 2-year photographic assessments. Risk ratios for mild/marked change were 1.70 (95% CI 1.26-2.29, p<0.001) for 30Gy and 1.15 (0.82-1.60, p=0.489) for 28.5Gy versus 50Gy. Three-year rates of physician-assessed moderate/marked adverse effects in the breast were 17.3% (13.3-22.3%, p<0.001) for 30Gy and 11.1% (7.9-15.6%, p=0.18) for 28.5Gy compared with 9.5% (6.5-13.7%) after 50Gy. With a median follow-up in survivors of 37.3months, 2 local tumour relapses and 23 deaths have occurred.

CONCLUSIONS

At 3years median follow-up, 28.5Gy in 5 fractions is comparable to 50Gy in 25 fractions, and significantly milder than 30Gy in 5 fractions, in terms of adverse effects in the breast.

Response-probability volume histograms and iso-probability of response charts in treatment plan evaluation

PURPOSE

This study aims at demonstrating a new method for treatment plan evaluation and comparison based on the radiobiological response of individual voxels. This is performed by applying them on three different cancer types and treatment plans of different conformalities. Furthermore, their usefulness is examined in conjunction with traditionally applied radiobiological and dosimetric treatment plan evaluation criteria.

METHODS

Three different cancer types (head and neck, breast and prostate) were selected to quantify the benefits of the proposed treatment plan evaluation method. In each case, conventional conformal radiotherapy (CRT) and intensity modulated radiotherapy (IMRT) treatment configurations were planned. Iso-probability of response charts was produced by calculating the response probability in every voxel using the linear-quadratic-Poisson model and the dose-response parameters of the corresponding structure to which this voxel belongs. The overall probabilities of target and normal tissue responses were calculated using the Poisson and the relative seriality models, respectively. The 3D dose distribution converted to a 2 Gy fractionation, D2(GY) and iso-BED distributions are also shown and compared with the proposed methodology. Response-probability volume histograms (RVH) were derived and compared with common dose volume histograms (DVH). The different dose distributions were also compared using the complication-free tumor control probability, P+, the biologically effective uniform dose, D, and common dosimetric criteria.

RESULTS

3D Iso-probability of response distributions is very useful for plan evaluation since their visual information focuses on the doses that are likely to have a larger clinical effect in that particular organ. The graphical display becomes independent of the prescription dose highlighting the local radiation therapy effect in each voxel without the loss of important spatial information. For example, due to the exponential nature of the Poisson distribution, cold spots in the target volumes or hot spots in the normal tissues are much easier to be identified. Response-volume histograms, as DVH, can also be derived and used for plan comparison. RVH are advantageous since by incorporating the radiobiological properties of each voxel they summarize the 3D distribution into 2D without the loss of relevant information. Thus, more clinically relevant radiobiological objectives and constraints could be defined and used in treatment planning optimization. These measures become increasingly important when dose distributions need to be designed according to the microscopic biological properties of tumor and normal tissues.

CONCLUSIONS

The proposed methods do not aim to replace quantifiers like the probabilities of total tissue response, which ultimately are the quantities of interest to evaluate treatment success. However, iso-probability of response charts and response-probability volume histograms illustrates more clearly the difference in effectiveness between different treatment plans than the information provided by alternative dosimetric data. The use of 3D iso-probability of response distributions could serve as a good descriptor of the effectiveness of a dose distribution indicating primarily the regions in a tissue that dominate its response.

Analysis of acute radiation-induced esophagitis in non-small-cell lung cancer patients using the Lyman NTCP model

PURPOSE

To analyze acute esophagitis (AE) in a Chinese population receiving 3D conformal radiotherapy (3DCRT) for non-small cell lung cancer (NSCLC), combined or not with chemotherapy (CT), using the Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) model.

MATERIALS AND METHODS

157 Chinese patients (pts) presented with NSCLC received 3DCRT: alone (34 pts) or combined with sequential CT (59 pts) (group 1) or with concomitant CT (64 pts) (group 2). Parameters (TD(50), n, and m) of the LKB NTCP model predicting for>grade 2 AE (RTOG grading) were identified using maximum likelihood analysis. Univariate and multivariate analyses using a binary regression logistic model were performed to identify patient, tumor and dosimetric predictors of AE.

RESULTS

Grade 2 or 3 AE occurred in 24% and 52% of pts in group 1 and 2, respectively (p<0.001). For the 93 group 1 pts, the fitted LKB model parameters were: m=0.15, n=0.29 and TD(50)=46 Gy. For the 64 group 2 pts, the parameters were: m=0.42, n=0.09 and TD(50)=36 Gy. In multivariate analysis, the only significant predictors of AE were: NTCP (p<0.001) and V(50), as continuous variable (RR=1.03, p=0.03) or being more than a threshold value of 11% (RR=3.6, p=0.009).

CONCLUSIONS

A LKB NTCP model has been established to predict AE in a Chinese population, receiving thoracic RT, alone or combined with CT. The parameters of the models appear slightly different than the previous one described in Western countries, with a lower volume effect for Chinese patients.