Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study.

GOAL: In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes. METHODS/MATERIALS: We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans. RESULTS: Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields. CONCLUSION: The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study.

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

Is motherhood still possible after pelvic carbon ion radiotherapy? A promising combined fertility-preservation approach.

INTRODUCTION: Preserving the endocrine and reproductive function in young female cancer patients undergoing pelvic radiation is a significant challenge. While the photon beam radiation's adverse effects on the uterus and ovaries are well established, the impact of pelvic carbon ion radiotherapy on women's reproductive function is largely unexplored. Strategies such as oocyte cryopreservation and ovarian transposition are commonly recommended for safeguarding future fertility. METHODS: This study presents a pioneering case of successful pregnancy after carbon ion radiotherapy for locally advanced sacral chondrosarcoma. RESULTS: A multidisciplinary approach facilitated the displacement of ovaries and uterus before carbon ion radiotherapy, resulting in the preservation of endocrine and reproductive function. CONCLUSION: The patient achieved optimal oncological response and delivered a healthy infant following the completion of cancer treatment.

The Role of Carbon Ion Therapy in the Changing Oncology Landscape-A Narrative Review of the Literature and the Decade of Carbon Ion Experience at the Italian National Center for Oncological Hadrontherapy.

BACKGROUND: Currently, 13 Asian and European facilities deliver carbon ion radiotherapy (CIRT) for preclinical and clinical activity, and, to date, 55 clinical studies including CIRT for adult and paediatric solid neoplasms have been registered. The National Center for Oncological Hadrontherapy (CNAO) is the only Italian facility able to accelerate both protons and carbon ions for oncological treatment and research. METHODS: To summarise and critically evaluate state-of-the-art knowledge on the application of carbon ion radiotherapy in oncological settings, the authors conducted a literature search till December 2022 in the following electronic databases: PubMed, Web of Science, MEDLINE, Google Scholar, and Cochrane. The results of 68 studies are reported using a narrative approach, highlighting CNAO's clinical activity over the last 10 years of CIRT. RESULTS: The ballistic and radiobiological hallmarks of CIRT make it an effective option in several rare, radioresistant, and difficult-to-treat tumours. CNAO has made a significant contribution to the advancement of knowledge on CIRT delivery in selected tumour types. CONCLUSIONS: After an initial ramp-up period, CNAO has progressively honed its clinical, technical, and dosimetric skills. Growing engagement with national and international networks and research groups for complex cancers has led to increasingly targeted patient selection for CIRT and lowered barriers to facility access.

Refining critical structure contouring in STereotactic Arrhythmia Radioablation (STAR): Benchmark results and consensus guidelines from the STOPSTORM.eu consortium.

BACKGROUND AND PURPOSE: In patients with recurrent ventricular tachycardia (VT), STereotactic Arrhythmia Radioablation (STAR) shows promising results. The STOPSTORM.eu consortium was established to investigate and harmonise STAR treatment in Europe. The primary goals of this benchmark study were to standardise contouring of organs at risk (OAR) for STAR, including detailed substructures of the heart, and accredit each participating centre. MATERIALS AND METHODS: Centres within the STOPSTORM.eu consortium were asked to delineate 31 OAR in three STAR cases. Delineation was reviewed by the consortium expert panel and after a dedicated workshop feedback and accreditation was provided to all participants. Further quantitative analysis was performed by calculating DICE similarity coefficients (DSC), median distance to agreement (MDA), and 95th percentile distance to agreement (HD95). RESULTS: Twenty centres participated in this study. Based on DSC, MDA and HD95, the delineations of well-known OAR in radiotherapy were similar, such as lungs (median DSC = 0.96, median MDA = 0.1 mm and median HD95 = 1.1 mm) and aorta (median DSC = 0.90, median MDA = 0.1 mm and median HD95 = 1.5 mm). Some centres did not include the gastro-oesophageal junction, leading to differences in stomach and oesophagus delineations. For cardiac substructures, such as chambers (median DSC = 0.83, median MDA = 0.2 mm and median HD95 = 0.5 mm), valves (median DSC = 0.16, median MDA = 4.6 mm and median HD95 = 16.0 mm), coronary arteries (median DSC = 0.4, median MDA = 0.7 mm and median HD95 = 8.3 mm) and the sinoatrial and atrioventricular nodes (median DSC = 0.29, median MDA = 4.4 mm and median HD95 = 11.4 mm), deviations between centres occurred more frequently. After the dedicated workshop all centres were accredited and contouring consensus guidelines for STAR were established. CONCLUSION: This STOPSTORM multi-centre critical structure contouring benchmark study showed high agreement for standard radiotherapy OAR. However, for cardiac substructures larger disagreement in contouring occurred, which may have significant impact on STAR treatment planning and dosimetry evaluation. To standardize OAR contouring, consensus guidelines for critical structure contouring in STAR were established.

STereotactic Arrhythmia Radioablation (STAR): the Standardized Treatment and Outcome Platform for Stereotactic Therapy Of Re-entrant tachycardia by a Multidisciplinary consortium (STOPSTORM.eu) and review of current patterns of STAR practice in Europe.

The EU Horizon 2020 Framework-funded Standardized Treatment and Outcome Platform for Stereotactic Therapy Of Re-entrant tachycardia by a Multidisciplinary (STOPSTORM) consortium has been established as a large research network for investigating STereotactic Arrhythmia Radioablation (STAR) for ventricular tachycardia (VT). The aim is to provide a pooled treatment database to evaluate patterns of practice and outcomes of STAR and finally to harmonize STAR within Europe. The consortium comprises 31 clinical and research institutions. The project is divided into nine work packages (WPs): (i) observational cohort; (ii) standardization and harmonization of target delineation; (iii) harmonized prospective cohort; (iv) quality assurance (QA); (v) analysis and evaluation; (vi, ix) ethics and regulations; and (vii, viii) project coordination and dissemination. To provide a review of current clinical STAR practice in Europe, a comprehensive questionnaire was performed at project start. The STOPSTORM Institutions' experience in VT catheter ablation (83% ≥ 20 ann.) and stereotactic body radiotherapy (59% > 200 ann.) was adequate, and 84 STAR treatments were performed until project launch, while 8/22 centres already recruited VT patients in national clinical trials. The majority currently base their target definition on mapping during VT (96%) and/or pace mapping (75%), reduced voltage areas (63%), or late ventricular potentials (75%) during sinus rhythm. The majority currently apply a single-fraction dose of 25 Gy while planning techniques and dose prescription methods vary greatly. The current clinical STAR practice in the STOPSTORM consortium highlights potential areas of optimization and harmonization for substrate mapping, target delineation, motion management, dosimetry, and QA, which will be addressed in the various WPs.

Immune checkpoint inhibitors and Carbon iON radiotherapy In solid Cancers with stable disease (ICONIC).

ICONIC is a multicenter, open-label, nonrandomized phase II clinical trial aiming to assess the feasibility and clinical activity of the addition of carbon ion radiotherapy to immune checkpoint inhibitors in cancer patients who have obtained disease stability with pembrolizumab administered as per standard-of-care. The primary end point is objective response rate, and the secondary end points are safety, survival and disease control rate. Translational research is an exploratory aim. The planned sample size is 27 patients. The study combination will be considered worth investigating if at least four objective responses are observed. If the null hypothesis is rejected, ICONIC will be the first proof of concept of the feasibility and clinical activity of the addition of carbon ion radiotherapy to immune checkpoint inhibitors in oncology.

ICONIC is a multicenter, open-label, nonrandomized, phase II clinical trial aiming to evaluate the feasibility and clinical activity of the addition of carbon ion radiotherapy to immune checkpoint inhibitors in cancer patients who have obtained disease stability with pembrolizumab administered as per standard-of-care. Considering that no clinical trials have been conducted thus far to assess the safety of the association between immune checkpoint inhibitors and carbon ion radiotherapy, the current clinical study will provide controlled data about the safety of this unprecedented therapeutic combination. Clinical Trial Registration: NCT05229614 (ClinicalTrials.gov).

Synthetic CT in Carbon Ion Radiotherapy of the Abdominal Site.

The generation of synthetic CT for carbon ion radiotherapy (CIRT) applications is challenging, since high accuracy is required in treatment planning and delivery, especially in an anatomical site as complex as the abdomen. Thirty-nine abdominal MRI-CT volume pairs were collected and a three-channel cGAN (accounting for air, bones, soft tissues) was used to generate sCTs. The network was tested on five held-out MRI volumes for two scenarios: (i) a CT-based segmentation of the MRI channels, to assess the quality of sCTs and (ii) an MRI manual segmentation, to simulate an MRI-only treatment scenario. The sCTs were evaluated by means of similarity metrics (e.g., mean absolute error, MAE) and geometrical criteria (e.g., dice coefficient). Recalculated CIRT plans were evaluated through dose volume histogram, gamma analysis and range shift analysis. The CT-based test set presented optimal MAE on bones (86.03 ± 10.76 HU), soft tissues (55.39 ± 3.41 HU) and air (54.42 ± 11.48 HU). Higher values were obtained from the MRI-only test set (MAEBONE = 154.87 ± 22.90 HU). The global gamma pass rate reached 94.88 ± 4.9% with 3%/3 mm, while the range shift reached a median (IQR) of 0.98 (3.64) mm. The three-channel cGAN can generate acceptable abdominal sCTs and allow for CIRT dose recalculations comparable to the clinical plans.

The role of multiple anatomical scenarios in plan optimization for carbon ion radiotherapy of pancreatic cancer.

PURPOSE /OBJECTIVE: To quantify benefits of robust optimization on multiple 4DCT acquisitions combined with off-line treatment adaptation for neoadjuvant carbon ion therapy (CIRT) of pancreatic cancer. MATERIAL/METHODS: For 10 previously treated patients, 4DCTs were acquired around -15 (CTPlan), -5 (RE1), -1 (RE2) and +6 (RE3) days from RT start. Treatment plans were newly optimized to a dose prescription of 38.4 Gy(RBE) (8 fractions) with a constraint of 38 Gy(RBE) to 1% of the gastrointestinal organs at risk volume (D1%). Three strategies were tested: (A) robust optimization on CTPlan maximum exhale (0Ex) with 3 mm set-up, 3% range uncertainty, including 30%-inhale; (B) addition of the RE1-0Ex scenario; (C) plan recalculation at each REi and adaptation (RPi) according to deviation thresholds from clinical goals. The cumulative variation of target coverage and GI-OARs doses was evaluated. Duodenum contours of all 4DCTs of each patient were registered on CTPlan-0Ex. The capacity of pre-RT acquisitions to predict duodenum position was investigated by computing the intersection of contours at CTplan, RE1, or their union, with respect to subsequent 4DCTs and the CTV, coupled with increasing margin. RESULTS: (A) No recalculation exceeded the D1% constraint. (B) The inclusion of RE1-0Ex in the optimization problem improved inter-fraction robustness on a patient-specific basis, but was non-significant on average. (C) Half of the plans would be re-optimized to recover target coverage and/or minimize duodenum dose, at least once. A significant difference was observed between pre-RT duodenum contours when intersecting subsequent contours, either with a margin expansion. CONCLUSION: Anatomical variations highlighted at multiple REi proved that a fast and efficient online adaptation is essential to optimize treatment quality of CIRT for pancreatic cancer.

The impact of particle radiotherapy on the functioning of cardiac implantable electronic devices: a systematic review of in vitro and in vivo studies according to PICO criteria.

The number of oncological patients who may benefit from proton beam radiotherapy (PBT) or carbon ion radiotherapy (CIRT), overall referred to as particle radiotherapy (RT), is expected to strongly increase in the next future, as well as the number of cardiological patients requiring cardiac implantable electronic devices (CIEDs). The management of patients with a CIED requiring particle RT deserves peculiar attention compared to those undergoing conventional photon beam RT, mostly due to the potential generation of secondary neutrons by particle beams interactions. Current consensus documents recommend managing these patients as being at intermediate/high risk of RT-induced device malfunctioning regardless of the dose on the CIED and the beam delivery method used, despite the last one significantly affects secondary neutrons generation (very limited neutrons production with active scanning as opposed to the passive scattering technique). The key issues for the current review were expressed in four questions according to the Population, Intervention, Control, Outcome criteria. Three in vitro and five in vivo studies were included. Based on the available data, PBT and CIRT with active scanning have a limited potential to interfere with CIED that has only emerged from in vitro study so far, while a significant potential for neutron-related, not severe, CIED malfunctions (resets) was consistently reported in both clinical and in vitro studies with passive scattering.

Role of diffusion-weighted MRI in recurrent rectal cancer treated with carbon ion radiotherapy.

Aim: To evaluate the association between pretreatment diffusion-weighted MRI (DW-MRI) and 12-month radiological response in locally recurrent rectal cancer treated with carbon ion radiotherapy. Methods: Histogram analysis was performed on pretreatment DW-MRI for patients re-irradiated with carbon ion radiotherapy for local recurrence of rectal cancer. Results: A total of 17 patients were enrolled in the study. Pretreatment DW-MRI b-value of 1000 s/mm2 (b1000) and apparent diffusion coefficient (ADC) lesion median values for 1-year nonresponders (six patients) and responders (11 patients) demonstrated a median (interquartile of median values) of 62.5 (23.9) and 34.0 (13.0) and 953.0 (277.0) and 942.5 (339.0) µm2/s, respectively. All b1000 histogram features (h-features) and ADC h-kurtosis showed statistically significant differences, whereas only b1000 h-median, b1000 h-interquartile range and ADC h-kurtosis demonstrated remarkable diagnostic accuracy. Conclusion: DW-MRI showed promising results in predicting carbon ion radiotherapy outcome in local recurrence of rectal cancer, particularly with regard to b1000 h-median, b1000 h-interquartile range and ADC h-kurtosis.

Carbon ion radiotherapy is a form of advanced radiotherapy that is especially suitable for radioresistant and/or difficult-to-irradiate tumors. In case of recurrence of rectal cancer after pelvic photon beam radiotherapy, carbon ion radiotherapy may be an option. In this study, the authors looked at the potential role of specific MRI sequences performed before treatment to predict response to carbon ion radiotherapy. If confirmed in a larger prospective cohort, the findings of this study may drive clinical decisions toward a more tumor- and patient-tailored therapeutic approach.

Cardiotoxicity model-based patient selection for Hodgkin lymphoma proton therapy.

INTRODUCTION: Hodgkin lymphoma (HL) is a highly curable hematological malignancy. Consolidation radiation therapy techniques have made significant progresses to improve organ-at-risk sparing in order to reduce late radiation-induced toxicity. Recent technical breakthroughs notably include intensity modulated proton therapy (IMPT), which has demonstrated a major dosimetric benefit at the cardiac level for mediastinal HL patients. However, its implementation in clinical practice is still challenging, notably due to the limited access to proton therapy facilities. In this context, the purpose of this study was to estimate the benefit of IMPT for HL proton therapy for diverse cardiac adverse events and to propose a general frame for mediastinal HL patient selection strategy for IMPT based on cardiotoxicity reduction, patient clinical factors, and IMPT treatment availability. MATERIAL AND METHODS: This retrospective dosimetric study included 30 mediastinal HL patients treated with VMAT. IMPT plans were generated on the initial simulation scans. Dose to the heart, to the left ventricle and to the valves were retrieved to calculate the relative risk (RR) of ischemic heart disease (IHD), congestive heart failure (CHF) and valvular disease (VD). Composite relative risk reduction (cRRR) of late cardiotoxicity, between VMAT and IMPT, were calculated as the weighted mean of relative risk reduction for IHD, CHF and VD, calculated across a wide range of cardiovascular risk factor combinations. The proportion of mediastinal HL patients who could benefit from IMPT was estimated in European countries, based on the country population and on the number of active gantries, to propose country-specific cRRR thresholds for patient selection. RESULTS: Compared with VMAT, IMPT significantly reduced average mean doses to the heart (2.36 Gy vs 0.99 Gy, p < 0.01), to the left ventricle (0.67 Gy vs 0.03, p < 0.01) and to the valves (1.29 Gy vs. 0.06, p < 0.01). For a HL patient without cardiovascular risk factor other than anthracycline-based chemotherapy, the relative risks of late cardiovascular complications were significantly lower after IMPT compared with VMAT for ischemic heart disease (1.07 vs 1.17, p < 0.01), for congestive heart failure (2.84 vs. 3.00, p < 0.01), and for valvular disease (1.01 vs. 1.06, p < 0.01). The median cRRR of cardiovascular adverse events with IMPT was 4.8%, ranging between 0.1% and 30.5%, depending on the extent of radiation fields and on the considered cardiovascular risk factors. The estimated proportion of HL patients currently treatable with IMPT in European countries with proton therapy facilities ranged between 8.0% and 100% depending on the country, corresponding to cRRR thresholds ranging from 24.0% to 0.0%. CONCLUSION: While a statistically significant clinical benefit is theoretically expected for ischemic heart disease, cardiac heart failure and valvular disease for mediastinal HL patients with IMPT, the overall cardiotoxicity risk reduction is notable only for a minority of patients. In the context of limited IMPT availability, this study proposed a general model-based selection approach for mediastinal HL patient based on calculated cardiotoxicity reduction, taking into consideration patient clinical characteristics and IMPT facility availability.

Case Report: Treatment Planning Study to Demonstrate Feasibility of Transthoracic Ultrasound Guidance to Facilitate Ventricular Tachycardia Ablation With Protons.

Background: Cardiac arrhythmias, such as ventricular tachycardia, are disruptions in the normal cardiac function that originate from problems in the electrical conduction of signals inside the heart. Recently, a non-invasive treatment option based on external photon or proton beam irradiation has been used to ablate the arrhythmogenic structures. Especially in proton therapy, based on its steep dose gradient, it is crucial to monitor the motion of the heart in order to ensure that the radiation dose is delivered to the correct location. Transthoracic ultrasound imaging has the potential to provide guidance during this treatment delivery. However, it has to be noted that the presence of an ultrasound probe on the chest of the patient introduces constraints on usable beam angles for both protons and photon treatments. This case report investigates the possibility to generate a clinically acceptable proton treatment plan while the ultrasound probe is present on the chest of the patient. Case: A treatment plan study was performed based on a 4D cardiac-gated computed tomography scan of a 55 year-old male patient suffering from refractory ventricular tachycardia who underwent cardiac radioablation. A proton therapy treatment plan was generated for the actual treatment target in presence of an ultrasound probe on the chest of this patient. The clinical acceptability of the generated plan was confirmed by evaluating standard target dose-volume metrics, dose to organs-at-risk and target dose conformity and homogeneity. Conclusion: The generation of a clinically acceptable proton therapy treatment plan for cardiac radioablation of ventricular tachycardia could be performed in the presence of an ultrasound probe on the chest of the patient. These results establish a basis and justification for continued research and product development for ultrasound-guided cardiac radioablation.

A randomized controlled phase III study comparing hadrontherapy with carbon ions versus conventional radiotherapy - including photon and proton therapy - for the treatment of radioresistant tumors: the ETOILE trial.

BACKGROUND: Some cancers such as sarcomas (bone and soft tissue sarcomas) and adenoid cystic carcinomas are considered as radioresistant to low linear energy transfer radiation (including photons and protons) and may therefore beneficiate from a carbon ion therapy. Despite encouraging results obtained in phase I/II trials compared to historical data with photons, the spread of carbon ions has been limited mainly because of the absence of randomized medical data. The French health authorities stressed the importance of having randomized data for carbon ion therapy. METHODS: The ETOILE study is a multicenter prospective randomized phase III trial comparing carbon ion therapy to either advanced photon or proton radiotherapy for inoperable or macroscopically incompletely resected (R2) radioresistant cancers including sarcomas and adenoid cystic carcinomas. In the experimental arm, carbon ion therapy will be performed at the National Center for Oncological Hadrontherapy (CNAO) in Pavia, Italy. In the control arm, photon or proton radiotherapy will be carried out in referent centers in France. The primary endpoint is progression-free survival (PFS). Secondary endpoints are overall survival and local control, toxicity profile, and quality of life. In addition, a prospective health-economic study and a radiobiological analysis will be conducted. To demonstrate an absolute improvement in the 5-year PFS rate of 20% in favor of carbon ion therapy, 250 patients have to be included in the study. DISCUSSION: So far, no clinical study of phase III has demonstrated the superiority of carbon ion therapy compared to conventional radiotherapy, including proton therapy, for the treatment of radioresistant tumors. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT02838602 . Date of registration: July 20, 2016. The posted information will be updated as needed to reflect protocol amendments and study progress.

Ablative Radiotherapy (ART) for Locally Advanced Pancreatic Cancer (LAPC): Toward a New Paradigm?

Locally advanced pancreatic cancer (LAPC) represents a major urgency in oncology. Due to the massive involvement of the peripancreatic vessels, a curative-intent surgery is generally precluded. Historically, LAPC has been an indication for palliative systemic therapy. In recent years, with the introduction of intensive multi-agent chemotherapy regimens and aggressive surgical approaches, the survival of LAPC patients has significantly improved. In this complex and rapidly evolving scenario, the role of radiotherapy is still debated. The use of standard-dose conventional fractionated radiotherapy in LAPC has led to unsatisfactory oncological outcomes. However, technological advances in radiation therapy over recent years have definitively changed this paradigm. The use of ablative doses of radiotherapy, in association with image-guidance, respiratory organ-motion management, and adaptive protocols, has led to unprecedented results in terms of local control and survival. In this overview, principles, clinical applications, and current pitfalls of ablative radiotherapy (ART) as an emerging treatment option for LAPC are discussed.

Time-resolved MRI for off-line treatment robustness evaluation in carbon-ion radiotherapy of pancreatic cancer.

PURPOSE: In this study, we investigate the use of magnetic resonance imaging (MRI) for the clinical evaluation of gating treatment robustness in carbon-ion radiotherapy (CIRT) of pancreatic cancer. Indeed, MRI allows radiation-free repeated scans and fast dynamic sequences for time-resolved (TR) imaging (cine-MRI), providing information on inter- and intra-fraction cycle-to-cycle variations of respiratory motion. MRI can therefore support treatment planning and verification, overcoming the limitations of the current clinical standard, that is, four-dimensional computed tomography (4DCT), which describes an "average" breathing cycle neglecting breathing motion variability. METHODS: We integrated a technique to generate a virtual CT (vCT) from 3D MRI with a method for 3D reconstruction from 2D cine-MRI, to produce TR vCTs for dose recalculations. For eight patients, the method allowed evaluating inter-fraction variations at end-exhale and intra-fraction cycle-to-cycle variability within the gating window in terms of tumor displacement and dose to the target and organs at risk. RESULTS: The median inter-fraction tumor motion was in the range 3.33-12.16 mm, but the target coverage was robust (-0.4% median D95% variation). Concerning cycle-to-cycle variations, the gating technique was effective in limiting tumor displacement (1.35 mm median gating motion) and corresponding dose variations (-3.9% median D95% variation). The larger exposure of organs at risk (duodenum and stomach) was caused by inter-fraction motion, whereas intra-fraction cycle-to-cycle dose variations were limited. CONCLUSIONS: This study proposed a method for the generation of TR vCTs from MRI, which enabled an off-line evaluation of gating treatment robustness and suggested its feasibility to support treatment planning of pancreatic tumors in CIRT.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors.

PURPOSE: To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario. MATERIALS AND METHODS: The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed. RESULTS: For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D95% and gastro-intestinal OARs D2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs. CONCLUSIONS: The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

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.