In-beam PET treatment monitoring of carbon therapy patients: Results of a clinical trial at CNAO.

PURPOSE: Carbon ion therapy treatments can be monitored non-invasively with in-beam Positron Emission Tomography (PET). At CNAO the INSIDE in-beam PET scanner has been used in a clinical trial (NCT03662373) to monitor cancer treatments with proton and carbon therapy. In this work we present the analysis results of carbon therapy data, acquired during the first phase of the clinical trial, analyzing data of nine patients treated at CNAO for various malignant tumors in the head-and-neck region. MATERIALS AND METHODS: The patient group contained two patients requiring replanning, and seven patients without replanning, based on established protocols. For each patient the PET images acquired along the course of treatment were compared with a reference, applying two analysis methods: the beam-eye-view (BEV) method and the γ-index analysis. Time trends in several parameters were investigated, as well as the agreement with control CTs, if available. RESULTS: Regarding the BEV-method, the average sigma value σ was 3.7 mm of range difference distributions for patients without changes (sensitivity of the INSIDE detector). The 3D-information obtained from the BEV analysis was partly in agreement with what was observed in the control CT. The data quality and quantity was insufficient for a definite interpretation of the time trends. CONCLUSION: We analyzed carbon therapy data acquired with the INSIDE in-beam PET detector using two analysis methods. The data allowed to evaluate sensitivity of the INSIDE detector for carbon therapy and to make several recommendations for the future.

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.

A Data-Driven Fragmentation Model for Carbon Therapy GPU-Accelerated Monte-Carlo Dose Recalculation.

The advent of Graphics Processing Units (GPU) has prompted the development of Monte Carlo (MC) algorithms that can significantly reduce the simulation time with respect to standard MC algorithms based on Central Processing Unit (CPU) hardware. The possibility to evaluate a complete treatment plan within minutes, instead of hours, paves the way for many clinical applications where the time-factor is important. FRED (Fast paRticle thErapy Dose evaluator) is a software that exploits the GPU power to recalculate and optimise ion beam treatment plans. The main goal when developing the FRED physics model was to balance accuracy, calculation time and GPU execution guidelines. Nowadays, FRED is already used as a quality assurance tool in Maastricht and Krakow proton clinical centers and as a research tool in several clinical and research centers across Europe. Lately the core software has been updated including a model of carbon ions interactions with matter. The implementation is phenomenological and based on carbon fragmentation data currently available. The model has been tested against the MC FLUKA software, commonly used in particle therapy, and a good agreement was found. In this paper, the new FRED data-driven model for carbon ion fragmentation will be presented together with the validation tests against the FLUKA MC software. The results will be discussed in the context of FRED clinical applications to 12C ions treatment planning.

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.

Radioguided surgery with ß- radiation in pancreatic Neuroendocrine Tumors: a feasibility study.

The possibility to use ß- decaying isotopes for radioguided surgery (RGS) has been recently proposed, and first promising tests on ex-vivo samples of Meningioma and intestinal Neuroendocrine Tumor (NET) have been published. This paper reports a study of the uptake of 68Ga-DOTATOC in pancreatic NETs (pNETs) in order to assess the feasibility of a new RGS approach using 90Y-DOTATOC. Tumor and healthy pancreas uptakes were estimated from 68Ga-DOTATOC PET/CT scans of 30 patients with pNETs. From the obtained SUVs (Standardised Uptake Value) and TNRs (Tumor Non tumor Ratio), an analysis algorithm relying on a Monte Carlo simulation of the detector has been applied to evaluate the performances of the proposed technique. Almost all considered patients resulted to be compatible with the application of ß--RGS assuming to administer 1.5 MBq/kg of activity of 90Y-DOTATOC 24 h before surgery, and a sampling time of few seconds. In just 2 cases the technique would have required a mildly increased amount of activity or of sampling time. Despite a high physiological uptake of 68Ga-DOTATOC in the healthy pancreas, the proposed RGS technique promises to be effective. This approach allows RGS to find application also in pancreatic diseases, where traditional techniques are not viable.

Radioguided surgery with ß radiation: a novel application with Ga68.

Radio Guided Surgery is a technique helping the surgeon in the resection of tumors: a radiolabeled tracer is administered to the patient before surgery and then the surgeon evaluates the completeness of the resection with a handheld detector sensitive to emitted radiation. Established methods rely on γ emitting tracers coupled with γ detecting probes. The efficacy of this technique is however hindered by the high penetration of γ radiation, limiting its applicability to low background conditions. To overtake such limitations, a novel approach to RGS has been proposed, relying on ß- emitting isotopes together with a dedicated ß probe. This technique has been proved to be effective in first ex-vivo trials. We discuss in this paper the possibility to extend its application cases to 68Ga, a ß+ emitting isotope widely used today in nuclear medicine. To this aim, a retrospective study on 45 prostatic cancer patients was performed, analysing their 68Ga-PSMA PET images to asses if the molecule uptake is enough to apply this technique. Despite the expected variability both in terms of SUV (median 4.1, IQR 3.0-6.1) and TNR (median 9.4, IQR 5.2-14.6), the majority of cases have been found to be compatible with ß-RGS with reasonable injected activity and probing time (5 s).

Feasibility of beta-particle radioguided surgery for a variety of "nuclear medicine" radionuclides.

PURPOSE: Beta-particle radioguided tumor resection may potentially overcome the limitations of conventional gamma-ray guided surgery by eliminating, or at least minimizing, the confounding effect of counts contributed by activity in adjacent normal tissues. The current study evaluates the clinical feasibility of this approach for a variety of radionuclides. Nowadays, the only ß- radioisotope suited to radioguided surgery is 90Y. Here, we study the ß- probe prototype capability to different radionuclides chosen among those used in nuclear medicine. METHODS: The counting efficiency of our probe prototype was evaluated for sources of electrons and photons of different energies. Such measurements were used to benchmark the Monte Carlo (MC) simulation of the probe behavior, especially the parameters related to the simulation of the optical photon propagation in the scintillation crystal. Then, the MC simulation was used to derive the signal and the background we would measure from a small tumor embedded in the patient body if one of the selected radionuclides is used. RESULTS: Based on the criterion of detectability of a 0.1 ml tumor for a counting interval of 1 s and an administered activity of 3 MBq/kg, the current probe yields a detectable signal over a wide range of Standard Uptake Values (SUVs) and tumor-to-non-tumor activity-concentration ratios (TNRs) for 31Si, 32P, 97Zr, and 188Re. Although efficient counting of 83Br, 133I, and 153Sm proved somewhat more problematic, the foregoing criterion can be satisfied for these isotopes as well for sufficiently high SUVs and TNRs.

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.