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.

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.

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.

Deep Seated Tumour Treatments With Electrons of High Energy Delivered at FLASH Rates: The Example of Prostate Cancer.

Different therapies are adopted for the treatment of deep seated tumours in combination or as an alternative to surgical removal or chemotherapy: radiotherapy with photons (RT), particle therapy (PT) with protons or even heavier ions like 12C, are now available in clinical centres. In addition to these irradiation modalities, the use of Very High Energy Electron (VHEE) beams (100-200 MeV) has been suggested in the past, but the diffusion of that technique was delayed due to the needed space and budget, with respect to standard photon devices. These disadvantages were not paired by an increased therapeutic efficacy, at least when comparing to proton or carbon ion beams. In this contribution we investigate how recent developments in electron beam therapy could reshape the treatments of deep seated tumours. In this respect we carefully explored the application of VHEE beams to the prostate cancer, a well-known and studied example of deep seated tumour currently treated with high efficacy both using RT and PT. The VHEE Treatment Planning System was obtained by means of an accurate Monte Carlo (MC) simulation of the electrons interactions with the patient body. A simple model of the FLASH effect (healthy tissues sparing at ultra-high dose rates), has been introduced and the results have been compared with conventional RT. The study demonstrates that VHEE beams, even in absence of a significant FLASH effect and with a reduced energy range (70-130 MeV) with respect to implementations already explored in literature, could be a good alternative to standard RT, even in the framework of technological developments that are nowadays affordable.

Addendum: Measurement of charged particle yields from PMMA irradiated by a 220 MeV/u 12C beam.

In this paper we report the re-analysis of the data published in (Piersanti et al. 2014) documenting the charged secondary particles production induced by the interaction of a 220 MeV/u 12C ion beam impinging on a polymethyl methacrylate (PMMA) target, measured in 2012 at the GSI facility in Darmstadt (Germany). This re-analysis takes into account the inhomogeneous light response of the LYSO crystal in the experimental setup measured in a subsequent experiment (2014) performed in the Heidelberg Ion- Beam Therapy Center. A better description of the detector and re-calculation of the geometrical efficiencies have been implemented as well, based on an improved approach that accounts also for the energy dependence of the emission spectrum. The new analysis has small effect on the total secondary charged flux, but has an impact on the production yield and emission velocity distributions of the different particle species (protons, deuterons and tritons) at different angles with respect to the beam direction (60° and 90°). All these observables indeed depend on the particle identification algorithms and hence on the LYSO detector energy response. The results of the data re-analysis presented here are intended to supersede and replace the results published in (Piersanti et al. 2014).

Benchmarking Geant4 hadronic models for prompt-γ monitoring in carbon ion therapy.

PURPOSE: The real-time monitoring of the spread-out Bragg peak would allow the planned dose delivered during treatment to be directly verified, but this poses a major challenge in modern ion beam therapy. A possible method to achieve this goal is to exploit the production of secondary particles by the nuclear reactions of the beam with the patient and correlate their emission profile to the planned target volume position. In this study, we present both the production rate and energy spectra of the prompt-γ produced by the interactions of the 12 C ion beam with a polymethyl methacrylate (PMMA) target. We also assess three different Monte Carlo models for prompt-γ simulation based on our experimental data. METHODS: The experiment was carried out at the GSI Helmholtz Centre for Heavy Ion Research, Darmstadt, Germany with a 220 MeV/u 12 C ions beam impinging on a 5× 5× 20 cm3 polymethyl methacrylate beam stopping target, with the prompt-γ being detected by a hexagonally-shaped barium fluoride scintillator with a circumscribed radius of 5.4 cm and a length of 14 cm, placed at 60° and 90° with respect to the beam direction. Monte Carlo simulations were carried out with three different hadronic models from the Geant4 code: binary ion cascade (BIC), quantum molecular dynamics (QMD), and Liege intranuclear cascade (INCL++ ). RESULTS: An experimental prompt-γ yield of 1.06 × 10-2  sr-1 was measured at 90°. A good agreement was observed between the shapes of the experimental and simulated energy spectra, especially with the INCL++ physics list. The prompt-γ yield obtained with this physics list was compatible with the measurement within 2σ, with a relative difference of 26% on average. BIC and QMD physics lists proved to be less accurate than INCL++ , with the difference between the measured and simulated yields exceeding 100%. The differences between the three physics lists were ascribed to important discrepancies between the models of the physical processes producing prompt-γ emissions. CONCLUSION: In conclusion, this study provides prompt-γ yield values in agreement with previously published results for different carbon ions energies. This work demonstrates that the INCL++ physics list from Geant4 is more accurate than BIC and QMD to reproduce prompt-γ emission properties.

Design of a new tracking device for on-line beam range monitor in carbon therapy.

Charged particle therapy is a technique for cancer treatment that exploits hadron beams, mostly protons and carbon ions. A critical issue is the monitoring of the beam range so to check the correct dose deposition to the tumor and surrounding tissues. The design of a new tracking device for beam range real-time monitoring in pencil beam carbon ion therapy is presented. The proposed device tracks secondary charged particles produced by beam interactions in the patient tissue and exploits the correlation of the charged particle emission profile with the spatial dose deposition and the Bragg peak position. The detector, currently under construction, uses the information provided by 12 layers of scintillating fibers followed by a plastic scintillator and a pixelated Lutetium Fine Silicate (LFS) crystal calorimeter. An algorithm to account and correct for emission profile distortion due to charged secondaries absorption inside the patient tissue is also proposed. Finally detector reconstruction efficiency for charged particle emission profile is evaluated using a Monte Carlo simulation considering a quasi-realistic case of a non-homogenous phantom.

Monitoring of Hadrontherapy Treatments by Means of Charged Particle Detection.

The interaction of the incoming beam radiation with the patient body in hadrontherapy treatments produces secondary charged and neutral particles, whose detection can be used for monitoring purposes and to perform an on-line check of beam particle range. In the context of ion-therapy with active scanning, charged particles are potentially attractive since they can be easily tracked with a high efficiency, in presence of a relatively low background contamination. In order to verify the possibility of exploiting this approach for in-beam monitoring in ion-therapy, and to guide the design of specific detectors, both simulations and experimental tests are being performed with ion beams impinging on simple homogeneous tissue-like targets (PMMA). From these studies, a resolution of the order of few millimeters on the single track has been proven to be sufficient to exploit charged particle tracking for monitoring purposes, preserving the precision achievable on longitudinal shape. The results obtained so far show that the measurement of charged particles can be successfully implemented in a technology capable of monitoring both the dose profile and the position of the Bragg peak inside the target and finally lead to the design of a novel profile detector. Crucial aspects to be considered are the detector positioning, to be optimized in order to maximize the available statistics, and the capability of accounting for the multiple scattering interactions undergone by the charged fragments along their exit path from the patient body. The experimental results collected up to now are also valuable for the validation of Monte Carlo simulation software tools and their implementation in Treatment Planning Software packages.

Time Evolution of DOTATOC Uptake in Neuroendocrine Tumors in View of a Possible Application of Radioguided Surgery with ß- Decay.

UNLABELLED: A novel radioguided surgery (RGS) technique exploiting ß- radiation has been proposed. To develop such a technique, a suitable radiotracer able to deliver a ß- emitter to the tumor has to be identified. A first candidate is represented by 90Y-labeled DOTATOC, a compound commonly used today for peptide radioreceptor therapy. The application of this ß- RGS to neuroendocrine tumors (NET) requires study of the uptake of DOTATOC and its time evolution both in tumors and in healthy organs and evaluation of the corresponding performance of the technique. METHODS: Uptake by lesions and healthy organs (kidneys, spleen, liver and healthy muscle) was estimated on 177Lu-DOTATOC SPECT/CT scans of 15 patients affected by NET with different localizations, treated at IRCCS-Arcispedale Santa Maria Nuova, Reggio Emilia, Italy. For each patient, SPECT/CT images, acquired at 0.5, 4, 20, 40, and 70 h after injection, were studied. For each lesion, the tumor-to-nontumor ratio (TNR) with respect to all healthy organs and its time evolution were studied. A subset of patients showing hepatic lesions was selected, and the TNR with respect to the nearby healthy tissue was calculated. By means of a Monte Carlo simulation of the probe for ß- RGS, the activity that is to be administered for a successful detection was estimated lesion-by-lesion. RESULTS: Uptake of DOTATOC on NETs maximized at about 24 h after injection. The cases of hepatic lesions showed a TNR with respect to the tumor margins compatible with the application of ß- RGS. In particular, 0.1-mL residuals are expected to be detectable within 1 s with 5% false-negative and 1% false-positive by administering the patient as little as 1 MBq/kg. CONCLUSION: The balance between tumor uptake and metabolic washout in healthy tissue causes the TNR to increase with time, reaching its maximum after 24 h, and this characteristic can be exploited when a radiotracer with a long half-life, such as 90Y, is used. In particular, if 90Y-DOTATOC is used with liver NET metastases, the proposed RGS technique is believed to be feasible by injecting an activity that is one third of that commonly used for PET imaging.

Toward radioguided surgery with ß- decays: uptake of a somatostatin analogue, DOTATOC, in meningioma and high-grade glioma.

UNLABELLED: A novel radioguided surgery (RGS) technique for cerebral tumors using ß(-) radiation is being developed. Checking for a radiotracer that can deliver a ß(-) emitter to the tumor is a fundamental step in the deployment of such a technique. This paper reports a study of the uptake of (90)Y-DOTATOC in meningiomas and high-grade gliomas (HGGs) and a feasibility study of the RGS technique in these types of tumor. Estimates were performed assuming the use of a ß(-) probe under development with a sensitive area 2.55 mm in radius to detect 0.1-mL residuals. METHODS: Uptake and background from healthy tissues were estimated on (68)Ga-DOTATOC PET scans of 11 meningioma patients and 12 HGG patients. A dedicated statistical analysis of the DICOM images was developed and validated. The feasibility study was performed using full simulation of emission and detection of the radiation, accounting for the measured uptake and background rate. RESULTS: All meningioma patients but one with an atypical extracranial tumor showed high uptake of DOTATOC. In terms of feasibility of the RGS technique, we estimated that by administering a 3 MBq/kg activity of radiotracer, the time needed to detect a 0.1-mL remnant with 5% false-negative and 1% false-positive rates is less than 1 s. Actually, to achieve a detection time of 1 s the required activities to administer were as low as 0.2-0.5 MBq/kg in many patients. In HGGs, the uptake was lower than in meningiomas, but the tumor-to-nontumor ratio was higher than 4, which implies that the tracer can still be effective for RGS. It was estimated that by administering 3 mBq/kg of radiotracer, the time needed to detect a 0.1-mL remnant is less than 6 s, with the exception of the only oligodendroma in the sample. CONCLUSION: Uptake of (90)Y-DOTATOC in meningiomas was high in all studied patients. Uptake in HGGs was significantly worse than in meningiomas but was still acceptable for RGS, particularly if further research and development are done to improve the performance of the ß(-) probe.