In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO.

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam's Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

First application of the BIANCA biophysical model to carbon-ion patient cases.

Objective.The main objective of this work consists of applying, for the first time, the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model to the RBE calculation for C-ion cancer patients, and comparing the outcomes with those obtained by the LEM I model, which is applied in clinics. Indeed, the continuous development of heavy-ion cancer therapy requires modelling of biological effects of ion beams on tumours and normal tissues. The relative biological effectiveness (RBE) of heavy ions is higher than that of protons, with a significant variation along the beam path. Therefore, it requires a precise modelling, especially for the pencil-beam scanning technique. Currently, two radiobiological models, LEM I and MKM, are in use for heavy ions in scanned pencil-beam facilities.Approach.Utilizing an interface with the FLUKA Particle Therapy Tool, BIANCA was applied to re-calculate the RBE-weighted dose distribution for carbon-ion treatment of three patients (chordoma, head-and-neck and prostate) previously irradiated at CNAO, where radiobiological optimization was based on LEM I. The predictions obtained by BIANCA were based either on chordoma cell survival (RBEsurv), or on dicentric aberrations in peripheral blood lymphocytes (RBEab), which are indicators of late normal tissue damage, including secondary tumours. The simulation outcomes were then compared with those provided by LEM I.Main results.While in the target and in the entrance channel BIANCA predictions were lower than those obtained by LEM I, the two models provided very similar results in the considered OAR. The observed differences betweenRBEsurvandRBEab(which were also dependent on fractional dose and LET) suggest that in normal tissues the information on cell survival should be integrated by information more closely related to the induction of late damage, such as chromosome aberrations.Significance.This work showed that BIANCA is suitable for treatment plan optimization in ion-beam therapy, especially considering that it can predict both cell survival and chromosome aberrations and has previously shown good agreement with carbon-ion experimental data.

Radiobiological damage by space radiation: extension of the BIANCA model to heavy ions up to iron, and pilot application to cosmic ray exposure.

Space research seems to be object of a renewed interest, also considering that human missions to the Moon, and possibly Mars, are being planned. Among the risks affecting such missions, astronauts' exposure to space radiation is a major concern. In this work, the question of the evaluation of biological damage by Galactic Cosmic Rays (GCR) was addressed by a biophysical model called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), which simulates the induction of cell death and chromosome aberrations by different ions. While previously BIANCA has been validated for calculating cell death along hadrontherapy beams up to oxygen, herein the approach was extended up to Fe ions. Specifically, experimental survival curves available in literature for V79 cells irradiated by Si-, Ne-, Ar- and Fe-ions were reproduced, and a reference radiobiological database describing V79 cell survival as a function of ion type (1 ⩽Z⩽ 26), energy and dose was constructed. Analogous databases were generated for Chinese hamster ovary hamster cells and human skin fibroblasts, finding good agreement between simulations and data. Concerning chromosome aberrations, which are regarded as radiation risk biomarkers, dicentric data in human lymphocytes irradiated by heavy ions up to iron were reproduced, and a radiobiological database allowing calculation of lymphocyte dicentric yields as a function of dose, ion type (1 ⩽Z⩽ 26) and energy was constructed. Following interface between BIANCA and the FLUKA Monte Carlo transport code, a feasibility study was performed to calculate the relative biological effectiveness (RBE) of different GCR spectrum components, for both dicentrics and cell death. Fe-ions, although representing only 10% of the total absorbed dose, were found to be responsible for about 35%-40% of the RBE-weighted dose. Interestingly, the RBE for dicentrics was higher than that for cell survival. More generally, this work shows that BIANCA can calculate RBE values for cell death and lymphocyte dicentrics not only for ion therapy, but also for space radiation.

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations.

PURPOSE: In-beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy-to-interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in-beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment. METHODS: We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel-wise two-tailed statistical tests of the simulated PET images, resembling the voxel-based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change. RESULTS: The VBM resembling method has been successfully applied to the simulated in-beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two-tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment. CONCLUSIONS: The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in-beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Healthy Tissue Damage Following Cancer Ion Therapy: A Radiobiological Database Predicting Lymphocyte Chromosome Aberrations Based on the BIANCA Biophysical Model.

Chromosome aberrations are widely considered among the best biomarkers of radiation health risk due to their relationship with late cancer incidence. In particular, aberrations in peripheral blood lymphocytes (PBL) can be regarded as indicators of hematologic toxicity, which is a major limiting factor of radiotherapy total dose. In this framework, a radiobiological database describing the induction of PBL dicentrics as a function of ion type and energy was developed by means of the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model, which has been previously applied to predict the effectiveness of therapeutic-like ion beams at killing tumour cells. This database was then read by the FLUKA Monte Carlo transport code, thus allowing us to calculate the Relative Biological Effectiveness (RBE) for dicentric induction along therapeutic C-ion beams. A comparison with previous results showed that, while in the higher-dose regions (e.g., the Spread-Out Bragg Peak, SOBP), the RBE for dicentrics was lower than that for cell survival. In the lower-dose regions (e.g., the fragmentation tail), the opposite trend was observed. This work suggests that, at least for some irradiation scenarios, calculating the biological effectiveness of a hadrontherapy beam solely based on the RBE for cell survival may lead to an underestimation of the risk of (late) damage to healthy tissues. More generally, following this work, BIANCA has gained the capability of providing RBE predictions not only for cell killing, but also for healthy tissue damage.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach.

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments' contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/ß values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments' contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

Monitoring Carbon Ion Beams Transverse Position Detecting Charged Secondary Fragments: Results From Patient Treatment Performed at CNAO.

Particle therapy in which deep seated tumours are treated using 12C ions (Carbon Ions RadioTherapy or CIRT) exploits the high conformity in the dose release, the high relative biological effectiveness and low oxygen enhancement ratio of such projectiles. The advantages of CIRT are driving a rapid increase in the number of centres that are trying to implement such technique. To fully profit from the ballistic precision achievable in delivering the dose to the target volume an online range verification system would be needed, but currently missing. The 12C ions beams range could only be monitored by looking at the secondary radiation emitted by the primary beam interaction with the patient tissues and no technical solution capable of the needed precision has been adopted in the clinical centres yet. The detection of charged secondary fragments, mainly protons, emitted by the patient is a promising approach, and is currently being explored in clinical trials at CNAO. Charged particles are easy to detect and can be back-tracked to the emission point with high efficiency in an almost background-free environment. These fragments are the product of projectiles fragmentation, and are hence mainly produced along the beam path inside the patient. This experimental signature can be used to monitor the beam position in the plane orthogonal to its flight direction, providing an online feedback to the beam transverse position monitor chambers used in the clinical centres. This information could be used to cross-check, validate and calibrate, whenever needed, the information provided by the ion chambers already implemented in most clinical centres as beam control detectors. In this paper we study the feasibility of such strategy in the clinical routine, analysing the data collected during the clinical trial performed at the CNAO facility on patients treated using 12C ions and monitored using the Dose Profiler (DP) detector developed within the INSIDE project. On the basis of the data collected monitoring three patients, the technique potential and limitations will be discussed.