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.

Analysis of tumour microstructure estimation from conventional diffusion MRI and application to skull-base chordoma. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:3761-3764. [PMID: 34892054 DOI: 10.1109/embc46164.2021.9630129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Skull-base chordoma (SBC) is a rare tumour whose molecular and radiological characteristics are still being investigated. In neuro-oncology microstructural imaging techniques, like diffusion-weighted MRI (DW-MRI), have been widely investigated, with the apparent diffusion coefficient (ADC) being one of the most used DW-MRI parameters due to its ease of acquisition and computation. ADC is a potential biomarker without a clear link to microstructure. The aim of this work was to derive microstructural information from conventional ADC, showing its potential for the characterisation of skull-base chordomas. Sixteen patients affected by SBC, who underwent conventional DW-MRI were retrospectively selected. From mono-exponential fits of DW-MRI, ADC maps were estimated using different sets of b-values. DW-MRI signals were simulated from synthetic substrates , which mimic the cellular packing of a tumour tissue with well-defined microstructural features. Starting from a published method, an error-driven procedure was evaluated to improve the estimates of microstructural parameters obtained through the simulated signals. A quantitative description of the tumour microstructure was then obtained from the DW-MRI images. This allowed successfully differentiating patients according to histologically-verified cell proliferation information.Clinical Relevance - The impact on cancer management derives from the expected improvement of radiation treatment quality tailored to a patient-specific non-invasive description of tumour microstructure.

Dosimetric Impact of Inter-Fraction Anatomical Changes in Carbon Ion Boost Treatment for High-Risk Prostate Cancer (AIRC IG 14300)

Rectum and bladder volumes play an important role in the dose distribution reproducibility in prostate cancer adenocarcinoma (PCa) radiotherapy, especially for particle therapy, where density variation can strongly affect the dose distribution. We investigated the reliability and reproducibility of our image-guided radiotherapy (IGRT) and treatment planning protocol for carbon ion radiotherapy (CIRT) within the phase II mixed beam study (AIRC IG 14300) for the treatment of high-risk PCa. In order to calculate the daily dose distribution, a set of synthetic computed tomography (sCT) images was generated from the cone beam computed tomography (CBCT) images acquired in each treatment session. Planning target volume (PTV) together with rectum and bladder volume variation was evaluated with sCT dose-volume histogram (DVH) metric deviations from the planning values. The correlations between the bladder and rectum volumes, and the corresponding DVH metrics, were also assessed. No significant difference in the bladder, rectum, and PTV median volumes between the planning computed tomography (pCT) and the sCT was found. In addition, no significant difference was assessed when comparing the average DVHs and median DVH metrics between pCT and sCT. Dose deviations determined by bladder and rectum filling variations demonstrated that dose distributions were reproducible in terms of both target coverage and organs at risk (OARs) sparing.

Image-based shading correction for narrow-FOV truncated pelvic CBCT with deep convolutional neural networks and transfer learning

Purpose: Cone beam computed tomography (CBCT) is a standard solution for in‐room image guidance for radiation therapy. It is used to evaluate and compensate for anatomopathological changes between the dose delivery plan and the fraction delivery day. CBCT is a fast and versatile solution, but it suffers from drawbacks like low contrast and requires proper calibration to derive density values. Although these limitations are even more prominent with in‐room customized CBCT systems, strategies based on deep learning have shown potential in improving image quality. As such, this article presents a method based on a convolutional neural network and a novel two‐step supervised training based on the transfer learning paradigm for shading correction in CBCT volumes with narrow field of view (FOV) acquired with an ad hoc in‐room system.

Methods: We designed a U‐Net convolutional neural network, trained on axial slices of corresponding CT/CBCT couples. To improve the generalization capability of the network, we exploited two‐stage learning using two distinct data sets. At first, the network weights were trained using synthetic CBCT scans generated from a public data set, and then only the deepest layers of the network were trained again with real‐world clinical data to fine‐tune the weights. Synthetic data were generated according to real data acquisition parameters. The network takes a single grayscale volume as input and outputs the same volume with corrected shading and improved HU values.

Results: Evaluation was carried out with a leave‐one‐out cross‐validation, computed on 18 unique CT/CBCT pairs from six different patients from a real‐world dataset. Comparing original CBCT to CT and improved CBCT to CT, we obtained an average improvement of 6 dB on peak signal‐to‐noise ratio (PSNR), +2% on structural similarity index measure (SSIM). The median interquartile range (IQR) Hounsfield unit (HU) difference between CBCT and CT improved from 161.37 (162.54) HU to 49.41 (66.70) HU. Region of interest (ROI)‐based HU difference was narrowed by 75% in the spongy bone (femoral head), 89% in the bladder, 85% for fat, and 83% for muscle. The improvement in contrast‐to‐noise ratio for these ROIs was about 67%.

Conclusions: We demonstrated that shading correction obtaining CT‐compatible data from narrow‐FOV CBCTs acquired with a customized in‐room system is possible. Moreover, the transfer learning approach proved particularly beneficial for such a shading correction approach.

Dosimetric effect of variable rectum and sigmoid colon filling during carbon ion radiotherapy to sacral chordoma

PURPOSE

Carbon ion radiotherapy (CIRT) is sensitive to anatomical density variations. We examined the dosimetric effect of variable intestinal filling condition during CIRT to ten sacral chordoma patients.

METHODS

For each patient, eight virtual computed tomography scans (vCTs) were generated by varying the density distribution within the rectum and the sigmoid in the planning computed tomography (pCT) with a density override approach mimicking a heterogeneous combination of gas and feces. Totally full and empty intestinal preparations were modelled. In addition, five different intestinal filling conditions were modelled by a mixed density pattern derived from two combined and weighted Gaussian distributions simulating gas and feces respectively. Finally, a patient-specific mixing proportion was estimated by evaluating the daily amount of gas detected in the cone beam computed tomography (CBCT). Dose distribution was recalculated on each vCT and dose volume histograms (DVHs) were examined.

RESULTS

No target coverage degradation was observed at different vCTs. Rectum and sigma dose degradation ranged respectively between: [-6.7; 21.6]GyE and [-0.7; 15.4]GyE for D_50%; [-377.4; 1197.9] and [-95.2; 1027.5] for AUC; [-1.2; 10.7]GyE and [-2.6; 21.5]GyE for D_1%.

CONCLUSIONS

Variation of intestinal density can greatly influence the penetration depth of charged particle and might compromise dose distribution. In particular cases, with large clinical target volume in very close proximity to rectum and sigmoid colon, it is appropriate to evaluate the amount of gas present in the daily CBCT images even if it is totally included in the reference planning structures.

Potential role of functional imaging in predicting outcome for patients treated with carbon ion therapy: a review

OBJECTIVE

Carbon ion radiation therapy (CIRT) is an emerging radiation technique with advantageous physical and radiobiologic properties compared to conventional radiotherapy (RT) providing better response in case of radioresistant and hypoxic tumors. Our aim is to critically review if functional imaging techniques could play a role in predicting outcome of CIRT-treated tumors, as already proven for conventional RT.

METHODS

14 studies, concerning Magnetic resonance imaging (MRI) and Positron Emission Tomography (PET), were selected after a comprehensive search on multiple electronic databases from January 2000 to March 2020.

RESULTS

MRI studies (n = 5) focused on diffusion-weighted MRI and, even though quantitative parameters were the same in all studies (apparent diffusion coefficient, ADC), results were not univocal, probably due to different imaging acquisition protocols and tumoral histology. For PET studies (n = 9), different tracers were used such as [¹⁸F]FDG and other uncommon tracers ([¹¹C]MET, [¹⁸F]FLT), with a relevant heterogeneity regarding parameters used for outcome assessment.

CONCLUSION

No conclusion can be drawn on the predictive value of functional imaging in CIRT-treated tumors. A standardization of image acquisition, multi-institutional large trials and external validations are needed in order to establish the prognostic value of functional imaging in CIRT and to guide clinical practice.

ADVANCES IN KNOWLEDGE

Emerging studies focused on functional imaging's role in predicting CIRT outcome. Due to the heterogeneity of images acquisition and studies, results are conflicting and prospective large studies with imaging standardized protocol are needed.

Non-invasive recognition of eye torsion through optical imaging of the iris pattern in ocular proton therapy

The introduction of non-invasive imaging techniques such as MRI imaging for treatment planning and optical eye tracking for in-room eye localization would obviate the requirement of clips implantation for many patients undergoing ocular proton therapy. This study specifically addresses the issue of torsional eye movement detection during patient positioning. Non-invasive detection of eye torsion is performed by measuring the iris pattern rotations using a beams eye view optical camera. When handling images of patients to be treated using proton therapy, a number of additional challenges are encountered, such as changing eye position, pupil dilatation and illumination. A method is proposed to address these extra challenges while also compensating for the effect of cornea distortion in eye torsion computation. The accuracy of the proposed algorithm was evaluated against corresponding measurement of eye torsion using the clips configuration measured on x-ray images. This study involves twenty patients who received ocular proton therapy at Paul Scherrer Institute and it is covered by ethical approval (EKNZ 2019-01987).

Is age rating enough to investigate changes in breathing motion pattern associated with aging of physically active women?

The most common way to analyze the effect of aging on breathing is to divide subjects into age groups. However, in addition to the fact that there is no consensus in the literature regarding age group division, such design critically influences the interpretation of the effects attributed to aging. Thus, this study aimed to investigate the feasibility to distinguish different age groups from the 3D kinematic variables of breathing motion (i.e., markers' coordinate as a function of time allowing the calculation of compartmental volume variations) and to analyze whether the aging could influence these variables. Seventy-three physically active women aged 19-80 years performed quiet breathing and vital capacity maneuvers. To record the thoracoabdominal breathing motion, the 3D coordinates of 32 retroreflective markers positioned on the trunk were used to estimate the volume variation of the superior thorax, inferior thorax, and abdomen. The percentage of contribution and the correlation coefficient were calculated to analyze the breathing motion pattern from the estimated volumes. The k-means cluster analysis was performed to analyze the age group classification. Linear regression was performed to investigate whether age can predict changes in the breathing motion pattern. The results showed that physically active women could not be classified into age groups from breathing motion. Despite significant p values of the linear regression, the high variability of the data suggested that age itself is not enough to predict the changes in breathing motion pattern when non-sedentary women are considered.

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 ¹²C 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 ¹²C 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 ¹²C 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.

Investigating DWI changes in white matter of meningioma patients treated with proton therapy

PURPOSE

To evaluate changes in diffusion and perfusion-related properties of white matter (WM) induced by proton therapy, which is capable of a greater dose sparing to organs at risk with respect to conventional X-ray radiotherapy, and to eventually expose early manifestations of delayed neuro-toxicities.

METHODS

Apparent diffusion coefficient (ADC) and IVIM parameters (D, D* and f) were estimated from diffusion-weighted MRI (DWI) in 46 patients affected by meningioma and treated with proton therapy. The impact on changes in diffusion and perfusion-related WM properties of dose and time, as well as the influence of demographic and pre-treatment clinical information, were investigated through linear mixed-effects models.

RESULTS

Decreasing trends in ADC and D were found for WM regions hit by medium-high (30-40 Gy(RBE)) and high (>40 Gy(RBE)) doses, which are compatible with diffusion restriction due to radiation-induced cellular injury. Significant influence of dose and time on median ADC changes were observed. Also, D* showed a significant dependency on dose, whereas f consistently showed no dependency on dose and time. Age, gender and surgery extent were also found to affect changes in ADC.

CONCLUSIONS

These results overall agree with those from studies conducted on cohorts of mixed proton and X-ray radiotherapy patients. Future work should focus on relating our findings with clinical information of co-morbidities and thus exploiting such or more advanced imaging data to build normal tissue complication probability models to better integrate clinical and dose information.

An MRI framework for respiratory motion modelling validation

Introduction

Respiratory motion models establish a correspondence between respiratory‐correlated (RC) 4‐dimensional (4D) imaging and respiratory surrogates, to estimate time‐resolved (TR) 3D breathing motion. To evaluate the performance of motion models on real patient data, a validation framework based on magnetic resonance imaging (MRI) is proposed, entailing the use of RC 4DMRI to build the model, and on both (i) TR 2D cine‐MRI and (ii) additional 4DMRI data for testing intra‐/inter‐fraction breathing motion variability.

Methods

Repeated MRI data were acquired in 7 patients with abdominal lesions. The considered model relied on deformable image registration (DIR) for building the model and compensating for inter‐fraction baseline variations. Both 2D and 3D validation were performed, by comparing model estimations with the ground truth 2D cine‐MRI and 4DMRI respiratory phases, respectively.

Results

The median DIR error was comparable to the voxel size (1.33 × 1.33 × 5 mm³), with higher values in the presence of large inter‐fraction motion (median value: 2.97 mm). In the 2D validation, the median estimation error on anatomical landmarks’ position resulted below 4 mm in every scenario, whereas in the 3D validation it was 1.33 mm and 4.21 mm when testing intra‐ and inter‐fraction motion, respectively. The range of motion described in the cine‐MRI was comparable to the motion of the building 4DMRI, being always above the estimation error. Overall, the model performance was dependent on DIR error, presenting reduced accuracy when inter‐fraction baseline variations occurred.

Conclusions

Results suggest the potential of the proposed framework in evaluating global motion models for organ motion management in MRI‐guided radiotherapy.

Development and validation of a new set-up simulator dedicated to ocular proton therapy at CNAO

An Eye Tracking System (ETS) is used at CNAO for providing a stable and reproducible ocular proton therapy (OPT) set-up, featuring a fixation light (FL) and monitoring stereo-cameras embedded in a rigid case. The aim of this work is to propose an ETS set-up simulation algorithm, that automatically provides the FL positioning in space, according to patient-specific gaze direction and avoiding interferences with patient, beam and collimator. Two configurations are provided: one in the CT room for acquiring images required for treatment planning with the patient lying on a couch, and one related to the treatment room with the patient sitting in front of the beam. Algorithm validation was performed reproducing ETS simulation (CT) and treatment (room) set-up for 30 patients previously treated at CNAO. The positioning accuracy of the device was quantified through a set of 14 control points applied to the ETS case and localizable both in the CT volume and in room X-ray images. Differences between the position of ETS reference points estimated by the algorithm and those measured by imaging systems are reported. The corresponding gaze direction deviation is on average 0.2° polar and 0.3° azimuth for positioning in CT room and 0.1° polar and 0.4° azimuth in the treatment room. The simulation algorithm was embedded in a clinically usable software application, which we assessed as capable of ensuring ETS positioning with an average accuracy of 2 mm in CT room and 1.5 mm in treatment room, corresponding to gaze direction deviations consistently lower than 1°.

Porcine lung phantom-based validation of estimated 4D-MRI using orthogonal cine imaging for low-field MR-Linacs

Real-time motion monitoring of lung tumors with low-field magnetic resonance imaging-guided linear accelerators (MR-Linacs) is currently limited to sagittal 2D cine magnetic resonance imaging (MRI). To provide input data for improved intrafractional and interfractional adaptive radiotherapy, the 4D anatomy has to be inferred from data with lower dimensionality. The purpose of this study was to experimentally validate a previously proposed propagation method that provides continuous time-resolved estimated 4D-MRI based on orthogonal cine MRI for a low-field MR-Linac. Ex vivo porcine lungs were injected with artificial nodules and mounted in a dedicated phantom that allows for the simulation of periodic and reproducible breathing motion. The phantom was scanned with a research version of a commercial 0.35 T MR-Linac. Respiratory-correlated 4D-MRI were reconstructed and served as ground truth images. Series of interleaved orthogonal slices in sagittal and coronal orientation, intersecting the injected targets, were acquired at 7.3 Hz. Estimated 4D-MRI at 3.65 Hz were created in post-processing using the propagation method and compared to the ground truth 4D-MRI. Eight datasets at different breathing frequencies and motion amplitudes were acquired for three porcine lungs. The overall median (95[Formula: see text] percentile) deviation between ground truth and estimated deformation vector fields was 2.3 mm (5.7 mm), corresponding to 0.7 (1.6) times the in-plane imaging resolution (3.5 × 3.5 mm²). Median (95[Formula: see text] percentile) estimated nodule position errors were 1.5 mm (3.8 mm) for nodules intersected by orthogonal slices and 2.1 mm (7.1 mm) for nodules located more than 2 cm away from either of the orthogonal slices. The estimation error depended on the breathing phase, the motion amplitude and the location of the estimated position with respect to the orthogonal slices. By using the propagation method, the 4D motion within the porcine lung phantom could be accurately and robustly estimated. The method could provide valuable information for treatment planning, real-time motion monitoring, treatment adaptation, and post-treatment evaluation of MR-guided radiotherapy treatments.

Improving the characterization of meningioma microstructure in proton therapy from conventional apparent diffusion coefficient measurements using Monte Carlo simulations of diffusion MRI

PURPOSE

Proton therapy could benefit from noninvasively gaining tumor microstructure information, at both planning and monitoring stages. The anatomical location of brain tumors, such as meningiomas, often hinders the recovery of such information from histopathology, and conventional noninvasive imaging biomarkers, like the apparent diffusion coefficient (ADC) from diffusion-weighted MRI (DW-MRI), are nonspecific. The aim of this study was to retrieve discriminative microstructural markers from conventional ADC for meningiomas treated with proton therapy. These markers were employed for tumor grading and tumor response assessment.

METHODS

DW-MRIs from patients affected by meningioma and enrolled in proton therapy were collected before (n = 35) and 3 months after (n = 25) treatment. For the latter group, the risk of an adverse outcome was inferred by their clinical history. Using Monte Carlo methods, DW-MRI signals were simulated from packings of synthetic cells built with well-defined geometrical and diffusion properties. Patients' ADC was modeled as a weighted sum of selected simulated signals. The weights that best described a patient's ADC were determined through an optimization procedure and used to estimate a set of markers of tumor microstructure: diffusion coefficient (D), volume fraction (vf), and radius (R). Apparent cellularity (ρ_app ) was estimated from vf and R for an easier clinical interpretability. Differences between meningothelial and atypical subtypes, and low- and high-grade meningiomas were assessed with nonparametric statistical tests, whereas sensitivity and specificity with ROC analyses. Similar analyses were performed for patients showing low or high risk of an adverse outcome to preliminary evaluate response to treatment.

RESULTS

Significant (P < 0.05) differences in median ADC, D, vf, R, and ρ_app values were found when comparing meningiomas' subtypes and grades. ROC analyses showed that estimated microstructural parameters reached higher specificity than ADC for subtyping (0.93 for D and vf vs 0.80 for ADC) and grading (0.75 for R vs 0.67 for ADC). High- and low-risk patients showed significant differences in ADC and microstructural parameters. The skewness of ρ_app was the parameter with highest AUC (0.90) and sensitivity (0.75).

CONCLUSIONS

Matching measured with simulated ADC yielded a set of potential imaging markers for meningiomas grading and response monitoring in proton therapy, showing higher specificity than conventional ADC. These markers can provide discriminative information about spatial patterns of tumor microstructure implying important advantages for patient-specific proton therapy workflows.

Radiomics and Dosiomics for Predicting Local Control after Carbon-Ion Radiotherapy in Skull-Base Chordoma

Skull-base chordoma (SBC) can be treated with carbon ion radiotherapy (CIRT) to improve local control (LC). The study aimed to explore the role of multi-parametric radiomic, dosiomic and clinical features as prognostic factors for LC in SBC patients undergoing CIRT. Before CIRT, 57 patients underwent MR and CT imaging, from which tumour contours and dose maps were obtained. MRI and CT-based radiomic, and dosiomic features were selected and fed to two survival models, singularly or by combining them with clinical factors. Adverse LC was given by in-field recurrence or tumour progression. The dataset was split in development and test sets and the models' performance evaluated using the concordance index (C-index). Patients were then assigned a low- or high-risk score. Survival curves were estimated, and risk groups compared through log-rank tests (after Bonferroni correction α = 0.0083). The best performing models were built on features describing tumour shape and dosiomic heterogeneity (median/interquartile range validation C-index: 0.80/024 and 0.79/0.26), followed by combined (0.73/0.30 and 0.75/0.27) and CT-based models (0.77/0.24 and 0.64/0.28). Dosiomic and combined models could consistently stratify patients in two significantly different groups. Dosiomic and multi-parametric radiomic features showed to be promising prognostic factors for LC in SBC treated with CIRT.

Integration of spatial distortion effects in a 4D computational phantom for simulation studies in extra-cranial MRI-guided radiation therapy: Initial results

PURPOSE

Spatial distortions in magnetic resonance imaging (MRI) are mainly caused by inhomogeneities of the static magnetic field, nonlinearities in the applied gradients, and tissue-specific magnetic susceptibility variations. These factors may significantly alter the geometrical accuracy of the reconstructed MR image, thus questioning the reliability of MRI for guidance in image-guided radiation therapy. In this work, we quantified MRI spatial distortions and created a quantitative model where different sources of distortions can be separated. The generated model was then integrated into a four-dimensional (4D) computational phantom for simulation studies in MRI-guided radiation therapy at extra-cranial sites.

METHODS

A geometrical spatial distortion phantom was designed in four modules embedding laser-cut PMMA grids, providing 3520 landmarks in a field of view of (345 × 260 × 480) mm³ . The construction accuracy of the phantom was verified experimentally. Two fast MRI sequences for extra-cranial imaging at 1.5 T were investigated, considering axial slices acquired with online distortion correction, in order to mimic practical use in MRI-guided radiotherapy. Distortions were separated into their sources by acquisition of images with gradient polarity reversal and dedicated susceptibility calculations. Such a separation yielded a quantitative spatial distortion model to be used for MR imaging simulations. Finally, the obtained spatial distortion model was embedded into an anthropomorphic 4D computational phantom, providing registered virtual CT/MR images where spatial distortions in MRI acquisition can be simulated.

RESULTS

The manufacturing accuracy of the geometrical distortion phantom was quantified to be within 0.2 mm in the grid planes and 0.5 mm in depth, including thickness variations and bending effects of individual grids. Residual spatial distortions after MRI distortion correction were strongly influenced by the applied correction mode, with larger effects in the trans-axial direction. In the axial plane, gradient nonlinearities caused the main distortions, with values up to 3 mm in a 1.5 T magnet, whereas static field and susceptibility effects were below 1 mm. The integration in the 4D anthropomorphic computational phantom highlighted that deformations can be severe in the region of the thoracic diaphragm, especially when using axial imaging with 2D distortion correction. Adaptation of the phantom based on patient-specific measurements was also verified, aiming at increased realism in the simulation.

CONCLUSIONS

The implemented framework provides an integrated approach for MRI spatial distortion modeling, where different sources of distortion can be quantified in time-dependent geometries. The computational phantom represents a valuable platform to study motion management strategies in extra-cranial MRI-guided radiotherapy, where the effects of spatial distortions can be modeled on synthetic images in a virtual environment.

Inter-fractional monitoring of [Formula: see text]C ions treatments: results from a clinical trial at the CNAO facility

The high dose conformity and healthy tissue sparing achievable in Particle Therapy when using C ions calls for safety factors in treatment planning, to prevent the tumor under-dosage related to the possible occurrence of inter-fractional morphological changes during a treatment. This limitation could be overcome by a range monitor, still missing in clinical routine, capable of providing on-line feedback. The Dose Profiler (DP) is a detector developed within the INnovative Solution for In-beam Dosimetry in hadronthErapy (INSIDE) collaboration for the monitoring of carbon ion treatments at the CNAO facility (Centro Nazionale di Adroterapia Oncologica) exploiting the detection of charged secondary fragments that escape from the patient. The DP capability to detect inter-fractional changes is demonstrated by comparing the obtained fragment emission maps in different fractions of the treatments enrolled in the first ever clinical trial of such a monitoring system, performed at CNAO. The case of a CNAO patient that underwent a significant morphological change is presented in detail, focusing on the implications that can be drawn for the achievable inter-fractional monitoring DP sensitivity in real clinical conditions. The results have been cross-checked against a simulation study.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

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 D_95% 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 D_95% 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 D_95% and gastro-intestinal OARs D_2% 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.

Multi-parametric qualitative and quantitative MRI assessment as predictor of histological grading in previously treated meningiomas

PURPOSE

Meningiomas are mainly benign tumors, though a considerable proportion shows aggressive behaviors histologically consistent with atypia/anaplasia. Histopathological grading is usually assessed through invasive procedures, which is not always feasible due to the inaccessibility of the lesion or to treatment contraindications. Therefore, we propose a multi-parametric MRI assessment as a predictor of meningioma histopathological grading.

METHODS

Seventy-three patients with 74 histologically proven and previously treated meningiomas were retrospectively enrolled (42 WHO I, 24 WHO II, 8 WHO III) and studied with MRI including T2 TSE, FLAIR, Gradient Echo, DWI, and pre- and post-contrast T1 sequences. Lesion masks were segmented on post-contrast T1 sequences and rigidly registered to ADC maps to extract quantitative parameters from conventional DWI and intravoxel incoherent motion model assessing tumor perfusion. Two expert neuroradiologists assessed morphological features of meningiomas with semi-quantitative scores.

RESULTS

Univariate analysis showed different distributions (p < 0.05) of quantitative diffusion parameters (Wilcoxon rank-sum test) and morphological features (Pearson's chi-square; Fisher's exact test) among meningiomas grouped in low-grade (WHO I) and higher grade forms (WHO II/III); the only exception consisted of the tumor-brain interface. A multivariate logistic regression, combining all parameters showing statistical significance in the univariate analysis, allowed discrimination between the groups of meningiomas with high sensitivity (0.968) and specificity (0.925). Heterogeneous contrast enhancement and low ADC were the best independent predictors of atypia and anaplasia.

CONCLUSION

Our multi-parametric MRI assessment showed high sensitivity and specificity in predicting histological grading of meningiomas. Such an assessment may be clinically useful in characterizing lesions without histological diagnosis. Key points • When surgery and biopsy are not feasible, parameters obtained from both conventional and diffusion-weighted MRI can predict atypia and anaplasia in meningiomas with high sensitivity and specificity. • Low ADC values and heterogeneous contrast enhancement are the best predictors of higher grade meningioma.

External validation of radiomics-based predictive models in low-dose CT screening for early lung cancer diagnosis

PURPOSE

Low-dose CT screening allows early lung cancer detection, but is affected by frequent false positive results, inter/intra observer variation and uncertain diagnoses of lung nodules. Radiomics-based models have recently been introduced to overcome these issues, but limitations in demonstrating their generalizability on independent datasets are slowing their introduction to clinic. The aim of this study is to evaluate two radiomics-based models to classify malignant pulmonary nodules in low-dose CT screening, and to externally validate them on an independent cohort. The effect of a radiomics features harmonization technique is also investigated to evaluate its impact on the classification of lung nodules from a multicenter data.

METHODS

Pulmonary nodules from two independent cohorts were considered in this study; the first cohort (110 subjects, 113 nodules) was used to train prediction models, and the second cohort (72 nodules) to externally validate them. Literature-based radiomics features were extracted and, after feature selection, used as predictive variables in models for malignancy identification. An in-house prediction model based on artificial neural network (ANN) was implemented and evaluated, along with an alternative model from the literature, based on a support vector machine (SVM) classifier coupled with a least absolute shrinkage and selection operator (LASSO). External validation was performed on the second cohort to evaluate models' generalization ability. Additionally, the impact of the Combat harmonization method was investigated to compensate for multicenter datasets variabilities. A new training of the models based on harmonized features was performed on the first cohort, then tested separately on the harmonized and non-harmonized features of the second cohort.

RESULTS

Preliminary results showed a good accuracy of the investigated models in distinguishing benign from malignant pulmonary nodules with both sets of radiomics features (i.e., non-harmonized and harmonized). The performance of the models, quantified in terms of Area Under the Curve (AUC), was > 0.89 in the training set and > 0.82 in the external validation set for all the investigated scenarios, outperforming the clinical standard (AUC of 0.76). Slightly higher performance was observed for the SVM-LASSO model than the ANN in the external dataset, although they did not result significantly different. For both harmonized and non-harmonized features, no statistical difference was found between Receiver operating characteristic (ROC) curves related to training and test set for both models.

CONCLUSIONS

Although no significant improvements were observed when applying the Combat harmonization method, both in-house and literature-based models were able to classify lung nodules with good generalization to an independent dataset, thus showing their potential as tools for clinical decision-making in lung cancer screening.