R. Ricotti

3D-printed applicators for high dose rate brachytherapy: Dosimetric assessment at different infill percentage

PURPOSE Dosimetric assessment of high dose rate (HDR) brachytherapy applicators, printed in 3D with acrylonitrile butadiene styrene (ABS) at different infill percentage. MATERIALS AND METHODS A low-cost, desktop, 3D printer (Hamlet 3DX100, Hamlet, Dublin, IE) was used for manufacturing simple HDR applicators, reproducing typical geometries in brachytherapy: cylindrical (common in vaginal treatment) and flat configurations (generally used to treat superficial lesions). Printer accuracy was investigated through physical measurements. The dosimetric consequences of varying the applicators density by tuning the printing infill percentage were analysed experimentally by measuring depth dose profiles and superficial dose distribution with Gafchromic EBT3 films (International Specialty Products, Wayne, NJ). Dose distributions were compared to those obtained with a commercial superficial applicator. RESULTS Measured printing accuracy was within 0.5mm. Dose attenuation was not sensitive to the density of the material. Surface dose distribution comparison of the 3D printed flat applicators with respect to the commercial superficial applicator showed an overall passing rate greater than 94% for gamma analysis with 3% dose difference criteria, 3mm distance-to-agreement criteria and 10% dose threshold. CONCLUSION Low-cost 3D printers are a promising solution for the customization of the HDR brachytherapy applicators. However, further assessment of 3D printing techniques and regulatory materials approval are required for clinical application.

Geometric and dosimetric accuracy and imaging dose of the real-time tumour tracking system of a gimbal mounted linac

PURPOSE To suggest a comprehensive testing scheme to evaluate the geometric and dosimetric accuracy and the imaging dose of the VERO dynamic tumour tracking (DTT) for its clinical implementation. METHODS Geometric accuracy was evaluated for gantry 0° and 90° in terms of prediction (EP), mechanical (EM) and tracking (ET) errors for sinusoidal patterns with 10 and 20 mm amplitudes, 2-6 s periods and phase shift up to 1 s and for 3 patient patterns. The automatic 4D model update was investigated simulating changes in the breathing pattern during treatment. Dosimetric accuracy was evaluated with gafchromic films irradiated in static and moving phantom with and without DTT. The entrance skin dose (ESD) was assessed using a solid state detector and gafchromic films. RESULTS The RMS of EP, EM, and ET were up to 0.8, 0.5 and 0.9 mm for all non phased-shifted motion patterns while for the phased-shifted ones, EP and ET increased to 2.2 and 2.6 mm. Up to 4 updates are necessary to restore a good correlation model, according to type of change. For 100 kVp and 1 mA s X-ray beam, the ESD per portal due to 20 s fluoroscopy was 16.6 mGy, while treatment verification at a frequency of 1 Hz contributed with 4.2 mGy/min. CONCLUSIONS The proposed testing scheme highlighted that the VERO DTT system tracks a moving target with high accuracy. The automatic update of the 4D model is a powerful tool to guarantee the accuracy of tracking without increasing the imaging dose.

Dosimetric characterization of 3D printed bolus at different infill percentage for external photon beam radiotherapy

BACKGROUND AND PURPOSE 3D printing is rapidly evolving and further assessment of materials and technique is required for clinical applications. We evaluated 3D printed boluses with acrylonitrile butadiene styrene (ABS) and polylactide (PLA) at different infill percentage. MATERIAL AND METHODS A low-cost 3D printer was used. The influence of the air inclusion within the 3D printed boluses was assessed thoroughly both with treatment planning system (TPS) and with physical measurements. For each bolus, two treatment plans were calculated with Monte Carlo algorithm, considering the computed tomography (CT) scan of the 3D printed bolus or modelling the 3D printed bolus as a virtual bolus structure with a homogeneous density. Depth dose measurements were performed with Gafchromic films. RESULTS High infill percentage corresponds to high density and high homogeneity within bolus material. The approximation of the bolus in the TPS as a homogeneous material is satisfying for infill percentages greater than 20%. Measurements performed with PLA boluses are more comparable to the TPS calculated profiles. For boluses printed at 40% and 60% infill, the discrepancies between calculated and measured dose distribution are within 5%. CONCLUSIONS 3D printing technology allows modulating the shift of the build-up region by tuning the infill percentage of the 3D printed bolus in order to improve superficial target coverage.

From technological advances to biological understanding: The main steps toward high-precision RT in breast cancer

Radiotherapy improves local control in breast cancer (BC) patients which increases overall survival in the long term. Improvements in treatment planning and delivery and a greater understanding of BC behaviour have laid the groundwork for high-precision radiotherapy, which is bound to further improve the therapeutic index. Precise identification of target volumes, better coverage and dose homogeneity have had a positive impact on toxicity and local control. The conformity of treatment dose due to three-dimensional radiotherapy and new techniques such as intensity modulated radiotherapy makes it possible to spare surrounding normal tissue. The widespread use of dose-volume constraints and histograms have increased awareness of toxicity. Real time image guidance has improved geometric precision and accuracy, together with the implementation of quality assurance programs. Advances in the precision of radiotherapy is also based on the choice of the appropriate fractionation and approach. Adaptive radiotherapy is not only a technical concept, but is also a biological concept based on the knowledge that different types of BC have distinctive patterns of locoregional spread. A greater understanding of cancer biology helps in choosing the treatment best suited to a particular situation. Biomarkers predictive of response play a crucial role. The combination of radiotherapy with molecular targeted therapies may enhance radiosensitivity, thus increasing the cytotoxic effects and improving treatment response. The appropriateness of an alternative fractionation, partial breast irradiation, dose escalating/de-escalating approaches, the extent of nodal irradiation have been examined for all the BC subtypes. The broadened concept of adaptive radiotherapy is vital to high-precision treatments.

Multimodal image registration for the identification of dominant intraprostatic lesion in high-precision radiotherapy treatments

PURPOSE The integration of CT and multiparametric MRI (mpMRI) is a challenging task in high-precision radiotherapy for prostate cancer. A simple methodology for multimodal deformable image registration (DIR) of prostate cancer patients is presented. METHODS CT and mpMRI of 10 patients were considered. Organs at risk and prostate were contoured on both scans. The dominant intraprostatic lesion was additionally delineated on MRI. After a preliminary rigid image registration, the voxel intensity of all the segmented structures in both scans except the prostate was increased by a specific amount (a constant additional value, A), in order to enhance the contrast of the main organs influencing its position and shape. 70 couples of scans were obtained by varying A from 0 to 800 and they were subsequently non-rigidly registered. Quantities derived from image analysis and contour statistics were considered for the tuning of the best performing A. RESULTS A = 200 resulted the minimum enhancement value required to obtain statistically significant superior registration results. Mean centre of mass distance between corresponding structures decreases from 7.4 mm in rigid registration to 5.3 mm in DIR without enhancement (DIR-0) and to 2.7 mm in DIR with A = 200 (DIR-200). Mean contour distance was 2.5, 1.9 and 0.67 mm in rigid registration, DIR-0 and DIR-200, respectively. In DIR-200 mean contours overlap increases of +13 and +24% with respect to DIR-0 and rigid registration, respectively. CONCLUSION Contour propagation according to the vector field resulting from DIR-200 allows the delineation of dominant intraprostatic lesion on CT scan and its use for high-precision radiotherapy treatment planning. Advances in knowledge: We investigated the application of a B-spline, mutual information-based multimodal DIR coupled with a simple, patient-unspecific but efficient contrast enhancement procedure in the pelvic body area, thus obtaining a robust and accurate methodology to transfer the functional information deriving from mpMRI onto a planning CT reference volume.

Intra-fraction respiratory motion and baseline drift during breast Helical Tomotherapy

BACKGROUND AND PURPOSE To investigate the intra-fraction breast motion during long-lasting treatments of breast cancer with Helical Tomotherapy by means of an optical tracking system. MATERIALS AND METHODS A set of seven radio-transparent passive markers was placed on the thoraco-abdominal surface of twenty breast cancer patients and tracked by an infrared tracking system. A continuous non-invasive monitoring of intra-fraction motion from patient setup verification and correction to the end of radiation delivery was thus obtained. The measured displacements were analysed in terms of cyclic respiratory motion and slow baseline drift. RESULTS The average monitoring time per patient was 15.57min. The breathing amplitude of the chest was less than 2mm, on average, along all anatomical directions. The baseline drift of the body led to more significant setup uncertainties than the respiratory motion. The main intra-fraction baseline drifts were in posterior and inferior directions and occurred within the first eight minutes of monitoring. Considering the intra-fraction motion only, the resultant clinical-to-planning target volume safety margins are highly patient-specific and largely anisotropic. CONCLUSION The non-respiratory motion occurring during prolonged treatments induces notable uncertainties. Non-invasive continuous monitoring of patient setup variations including baseline drifts is recommended in order to minimize dosimetric deviations, which might jeopardize the therapeutic ratio between target coverage and the sparing of organs at risk.

4D ML reconstruction as a tool for volumetric PET-based treatment verification in ion beam radiotherapy

PURPOSE An innovative strategy to improve the sensitivity of positron emission tomography (PET)-based treatment verification in ion beam radiotherapy is proposed. METHODS Low counting statistics PET images acquired during or shortly after the treatment (Measured PET) and a Monte Carlo estimate of the same PET images derived from the treatment plan (Expected PET) are considered as two frames of a 4D dataset. A 4D maximum likelihood reconstruction strategy was adapted to iteratively estimate the annihilation events distribution in a reference frame and the deformation motion fields that map it in the Expected PET and Measured PET frames. The outputs generated by the proposed strategy are as follows: (1) an estimate of the Measured PET with an image quality comparable to the Expected PET and (2) an estimate of the motion field mapping Expected PET to Measured PET. The details of the algorithm are presented and the strategy is preliminarily tested on analytically simulated datasets. RESULTS The algorithm demonstrates (1) robustness against noise, even in the worst conditions where 1.5 × 10(4) true coincidences and a random fraction of 73% are simulated; (2) a proper sensitivity to different kind and grade of mismatches ranging between 1 and 10 mm; (3) robustness against bias due to incorrect washout modeling in the Monte Carlo simulation up to 1/3 of the original signal amplitude; and (4) an ability to describe the mismatch even in presence of complex annihilation distributions such as those induced by two perpendicular superimposed ion fields. CONCLUSIONS The promising results obtained in this work suggest the applicability of the method as a quantification tool for PET-based treatment verification in ion beam radiotherapy. An extensive assessment of the proposed strategy on real treatment verification data is planned.

Variability in axillary lymph node delineation for breast cancer radiotherapy in presence of guidelines on a multi-institutional platform

Abstract Aim: To quantify the variability between radiation oncologists (ROs) when outlining axillary nodes in breast cancer. Material and methods: For each participating center, three ROs with different levels of expertise, i.e., junior (J), senior (S) and expert (E), contoured axillary nodal levels (L1, L2, L3 and L4) on the CT images of three different patients (P) of an increasing degree of anatomical complexity (from P1 to P2 to P3), according to contouring guidelines. Consensus contours were generated using the simultaneous truth and performance level estimation (STAPLE) method. Results: Fifteen centers and 42 ROs participated. Overall, the median Dice similarity coefficient was 0.66. Statistically significant differences were observed according to the level of expertise (better agreement for J and E, worse for S); the axillary level (better agreement for L1 and L4, worse for L3); the patient (better agreement for P1, worse for P3). Statistically significant differences in contouring were found in 18% of the inter-center comparison. Less than a half of the centers could claim to have a good agreement between the internal ROs. Conclusions: The overall intra-institute and inter-institute agreement was moderate. Central lymph-node levels were the most critical and variability increased as the complexity of the patient’s anatomy increased. These findings might have an effect on the interpretation of results from multicenter and even mono-institute studies.

Multi atlas based segmentation: should we prefer the best atlas group over the group of best atlases?

Multi atlas based segmentation (MABS) uses a database of atlas images, and an atlas selection process is used to choose an atlas subset for registration and voting. In the current state of the art, atlases are chosen according to a similarity criterion between the target subject and each atlas in the database. In this paper, we propose a new concept for atlas selection that relies on selecting the best performing group of atlases rather than the group of highest scoring individual atlases. Experiments were performed using CT images of 50 patients, with contours of brainstem and parotid glands. The dataset was randomly split into two groups: 20 volumes were used as an atlas database and 30 served as target subjects for testing. Classic oracle selection, where atlases are chosen by the highest dice similarity coefficient (DSC) with the target, was performed. This was compared to oracle group selection, where all the combinations of atlas subgroups were considered and scored by computing DSC with the target subject. Subsequently, convolutional neural networks were designed to predict the best group of atlases. The results were also compared with the selection strategy based on normalized mutual information (NMI). Oracle group was proven to be significantly better than classic oracle selection (p  <  10-5). Atlas group selection led to a median  ±  interquartile DSC of 0.740  ±  0.084, 0.718  ±  0.086 and 0.670  ±  0.097 for brainstem and left/right parotid glands respectively, outperforming NMI selection 0.676  ±  0.113, 0.632  ±  0.104 and 0.606  ±  0.118 (p  <  0.001) as well as classic oracle selection. The implemented methodology is a proof of principle that selecting the atlases by considering the performance of the entire group of atlases instead of each single atlas leads to higher segmentation accuracy, being even better then current oracle strategy. This finding opens a new discussion about the most appropriate atlas selection criterion for MABS.

First clinical investigation of a 4D maximum likelihood reconstruction for 4D PET-based treatment verification in ion beam therapy

BACKGROUND AND PURPOSE In clinical applications of Positron Emission Tomography (PET)-based treatment verification in ion beam therapy (PT-PET), detection and interpretation of inconsistencies between Measured PET and Expected PET are mostly limited by Measured PET noise, due to low count statistics, and by Expected PET bias, especially due to inaccurate washout modelling in off-line implementations. In this work, a recently proposed 4D Maximum Likelihood (ML) reconstruction algorithm which considers Measured PET and Expected PET as two different motion phases of a 4D dataset is assessed on clinical 4D PET-CT datasets acquired after carbon ion therapy. MATERIAL AND METHODS The 4D ML reconstruction algorithm estimates: (1) Measured PET of enhanced image quality with respect to the conventional Measured PET, thanks to the exploitation of Expected PET; (2) the deformation field mapping the Expected PET onto the Measured PET as a measure of the occurred displacements. RESULTS Results demonstrate the desired sensitivity to inconsistencies due to breathing motion and/or setup modification, robustness to noise in different count statistics scenarios, but a limited sensitivity to Expected PET washout inaccuracy. CONCLUSIONS The 4D ML reconstruction algorithm supports clinical 4D PT-PET in ion beam therapy. The limited sensitivity to washout inaccuracy can be detected and potentially overcome.