Consensus guide on CT-based prediction of stopping-power ratio using a Hounsfield look-up table for proton therapy

BACKGROUND AND PURPOSE

Studies have shown large variations in stopping-power ratio (SPR) prediction from computed tomography (CT) across European proton centres. To standardise this process, a step-by-step guide on specifying a Hounsfield look-up table (HLUT) is presented here.

MATERIALS AND METHODS

The HLUT specification process is divided into six steps: Phantom setup, CT acquisition, CT number extraction, SPR determination, HLUT specification, and HLUT validation. Appropriate CT phantoms have a head- and body-sized part, with tissue-equivalent inserts in regard to X-ray and proton interactions. CT numbers are extracted from a region-of-interest covering the inner 70% of each insert in-plane and several axial CT slices in scan direction. For optimal HLUT specification, the SPR of phantom inserts is measured in a proton beam and the SPR of tabulated human tissues is computed stoichiometrically at 100 MeV. Including both phantom inserts and tabulated human tissues increases HLUT stability. Piecewise linear regressions are performed between CT numbers and SPRs for four tissue groups (lung, adipose, soft tissue, and bone) and then connected with straight lines. Finally, a thorough but simple validation is performed.

RESULTS

The best practices and individual challenges are explained comprehensively for each step. A well-defined strategy for specifying the connection points between the individual line segments of the HLUT is presented. The guide was tested exemplarily on three CT scanners from different vendors, proving its feasibility.

CONCLUSION

The presented step-by-step guide for CT-based HLUT specification with recommendations and examples can contribute to reduce inter-centre variations in SPR prediction.

Proton therapy in the treatment of hepatocellular carcinoma

Liver cancer represents one of the most common causes of death from cancer worldwide. Hepatocellular carcinoma (HCC) accounts for 90% of all primary liver cancers. Among local therapies, evidence regarding the use of radiation therapy is growing. Proton therapy currently represents the most advanced radiation therapy technique with unique physical properties which fit well with liver irradiation. Here, in this review, we aim to 1) illustrate the rationale for the use of proton therapy (PT) in the treatment of HCC, 2) discuss the technical challenges of advanced PT in this disease, 3) review the major clinical studies regarding the use of PT for HCC, and 4) analyze the potential developments and future directions of PT in this setting.

ESTRO-ACROP guideline on surface guided radiation therapy

Surface guidance systems enable patient positioning and motion monitoring without using ionising radiation. Surface Guided Radiation Therapy (SGRT) has therefore been widely adopted in radiation therapy in recent years, but guidelines on workflows and specific quality assurance (QA) are lacking. This ESTRO-ACROP guideline aims to give recommendations concerning SGRT roles and responsibilities and highlights common challenges and potential errors. Comprehensive guidelines for procurement, acceptance, commissioning, and QA of SGRT systems installed on computed tomography (CT) simulators, C-arm linacs, closed-bore linacs, and particle therapy treatment systems are presented that will help move to a consensus among SGRT users and facilitate a safe and efficient implementation and clinical application of SGRT.

Experimental assessment of inter-centre variation in stopping-power and range prediction in particle therapy

PURPOSE

Experimental assessment of inter-centre variation and absolute accuracy of stopping-power-ratio (SPR) prediction within 17 particle therapy centres of the European Particle Therapy Network.

MATERIAL AND METHODS

A head and body phantom with seventeen tissue-equivalent materials were scanned consecutively at the participating centres using their individual clinical CT scan protocol and translated into SPR with their in-house CT-number-to-SPR conversion. Inter-centre variation and absolute accuracy in SPR prediction were quantified for three tissue groups: lung, soft tissues and bones. The integral effect on range prediction for typical clinical beams traversing different tissues was determined for representative beam paths for the treatment of primary brain tumours as well as lung and prostate cancer.

RESULTS

An inter-centre variation in SPR prediction (2σ) of 8.7%, 6.3% and 1.5% relative to water was determined for bone, lung and soft-tissue surrogates in the head setup, respectively. Slightly smaller variations were observed in the body phantom (6.2%, 3.1%, 1.3%). This translated into inter-centre variation of integral range prediction (2σ) of 2.9%, 2.6% and 1.3% for typical beam paths of prostate-, lung- and primary brain-tumour treatments, respectively. The absolute error in range exceeded 2% in every fourth participating centre. The consideration of beam hardening and the execution of an independent HLUT validation had a positive effect, on average.

CONCLUSION

The large inter-centre variations in SPR and range prediction justify the currently clinically used margins accounting for range uncertainty, which are of the same magnitude as the inter-centre variation. This study underlines the necessity of higher standardisation in CT-number-to-SPR conversion.

Technical Note: CT calibration for proton treatment planning by cross-calibration with proton CT data

PURPOSE

This study explores the possibility of a new method for x-ray computed tomography (CT) calibration by means of cross-calibration with proton CT (pCT) data. The proposed method aims at a more accurate conversion of CT Hounsfield Units (HU) into proton stopping power ratio (SPR) relative to water to be used in proton-therapy treatment planning.

METHODS

X-ray CT scan was acquired on a synthetic anthropomorphic phantom, composed of different tissue equivalent materials (TEMs). A pCT apparatus was instead adopted to obtain a reference three-dimensional distribution of the phantom's SPR values. After rigid registration, the x-ray CT was artificially blurred to the same resolution of pCT. Then a scatter plot showing voxel-by-voxel SPR values as a function of HU was employed to link the two measurements and thus obtaining a cross-calibrated x-ray CT calibration curve. The cross-calibration was tested at treatment planning system and then compared with a conventional calibration based on exactly the same TEMs constituting the anthropomorphic phantom.

RESULTS

Cross-calibration provided an accurate SPR mapping, better than by conventional TEMs calibration. The dose distribution of single beams optimized on the reference SPR map was recomputed on cross-calibrated CT, showing, with respect to conventional calibration, minor deviation at the dose fall-off (lower than 1%).

CONCLUSIONS

The presented data demonstrated that, by means of reference pCT data, a heterogeneous phantom can be used for CT calibration, paving the way to the use of biological samples, with their accurate description of patients' tissues. This overcomes the limitations of conventional CT calibration requiring homogenous samples, only available by synthetic TEMs, which fail in accurately mimicking the properties of biological tissues. Once a heterogeneous biological sample is provided with its corresponding reference SPR maps, a cross-calibration procedure could be adopted by other PT centers, even when not equipped with a pCT system.

Clinical implementation of pencil beam scanning proton therapy for liver cancer with forced deep expiration breath hold

PURPOSE

To present our technique for liver cancer treatments with proton therapy in pencil beam scanning mode and to evaluate the impact of uncertainties on plan quality.

MATERIALS AND METHODS

Seventeen patients affected by liver cancer were included in this study. Patients were imaged and treated in forced breath-hold using the Active Breathing Coordinator system and monitored with an optical tracking system. Three simulation CTs were acquired to estimate the anatomical variability between breath-holds and generate an internal target volume (ITV). The treatment plans were optimized with a Single Field Optimization technique aimed at minimizing the use of range shifter. Plan robustness was tested simulating systematic range and setup uncertainties, as well as the interplay effect between breath-holds. The appropriateness of margin was further verified based on the actual positioning data acquired during treatment.

RESULTS

The dose distributions of the nominal plans achieved a satisfactory target coverage in 11 out of 17 patients, while in the remaining 6 D₉₅ to the PTV was affected by the constraint on mean liver dose. The constraints for all other organs at risk were always within tolerances. The interplay effect had a limited impact on the dose distributions: the worst case scenario showed a D₉₅ reduction in the ITV < 3.9 GyRBE and no OAR with D₁ > 105% of the prescription dose. The robustness analysis showed that for 13 out of 17 patients the ITV coverage in terms of D₉₅ was better than D₉₅ of the PTV in the nominal plan. For the remaining 4 patients, the maximum difference between ITV D₉₅ and PTV D₉₅ was ≤0.7% even for the largest simulated setup error and it was deemed clinically acceptable. Hot spots in the OARs were always lower than 105% of the prescription dose. Positioning images confirmed that the breath hold technique and the PTV margin were adequate to compensate for inter- and intra-breath-hold variations in liver position.

CONCLUSION

We designed and clinically applied a technique for the treatment of liver cancer with proton pencil beam scanning in forced deep expiration breath-hold. The initial data on plan robustness and patient positioning suggest that the choices in terms of planning technique and treatment margins are able to reach the desired balance between target coverage and organ at risk sparing.

Proton or photon radiosurgery for cardiac ablation of ventricular tachycardia? Breath and ECG gated robust optimization

PURPOSE

Ventricular tachycardia (VT) is a life-threatening heart disorder. The aim of this preliminary study is to assess the feasibility of stereotactic body radiation therapy (SBRT) photon and proton therapy (PT) plans for the treatment of VT, adopting robust optimization technique for both irradiation techniques.

METHODS

ECG gated CT images (in breath hold) were acquired for one patient. Conventional planning target volume (PTV) and robust optimized plans (25GyE in single fraction) were simulated for both photon (IMRT, 5 and 9 beams) and proton (SFO, 2 beams) plans. Robust optimized plans were obtained both for protons and photons considering in the optimization setup errors (5 mm in the three orthogonal directions), range (±3.5%) and the clinical target volume (CTV) motion due to heartbeat and breath-hold variability.

RESULTS

The photon robust optimization method, compared to PTV-based plans, showed a reduction in the average dose to the heart by about 25%; robust optimization allowed also reducing the mean dose to the left lung from 3.4. to 2.8 Gy for 9-beams configuration and from 4.1 to 2.9 Gy for 5-beams configuration. Robust optimization with protons, allowed further reducing the OAR doses: average dose to the heart and to the left lung decreased from 7.3 Gy to 5.2 GyE and from 2.9 Gy to 2.2 GyE, respectively.

CONCLUSIONS

Our study demonstrates the importance of the optimization technique adopted in the treatment planning system for VT treatment. It has been shown that robust optimization can significantly reduce the dose to healthy cardiac tissues and that PT further increases this gain.

Clinical results of active scanning proton therapy for primary liver tumors

BACKGROUND

Evidence for the efficacy of radiation therapy for primary liver cancer is growing. In this context, proton therapy (PT) can potentially improve the therapeutic ratio, as demonstrated by recent clinical studies. Here we report the first European clinical experience on the use of PT for primary liver cancer.

METHODS

All patients treated for primary liver cancer in our center entered the analysis. Patients were simulated during deep expiration breath-hold. A 15-fraction treatment schedule was adopted using active scanning PT. Clinical outcome and toxicity were retrospectively analyzed.

RESULTS

Between January 2018 and December 2019, 18 patients were treated. Fourteen patients had hepatocellular carcinoma (HCC), three patients had intrahepatic cholangiocarcinoma (ICC), and one patient had synchronous ICC-HCC. The Child-Pugh score was A5 in the majority of patients with HCC (71.4%). Median prescription dose was 58.05 Gy (range, 50.31-67.5). Median follow-up was 10 months (range, 1-19). The majority of deaths occurred from liver tumor progression. One-year overall survival (OS) was 63%. A significant correlation between worse OS and patient performance status, vascular invasion, and tumor stage was recorded. One-year local control was 90%. Toxicity was low, with a decrease in Child-Pugh score ⩾2 points detected in one patient. No cases of classic radiation-induced liver disease occurred.

CONCLUSIONS

Our initial results of active scanning PT for primary liver cancer demonstrated the feasibility, safety, and effectiveness of this advanced technique in this setting. The potential of the combination of PT with other locoregional therapies is under evaluation.

Clinical implementation in proton therapy of multi-field optimization by a hybrid method combining conventional PTV with robust optimization

To implement a robust multi-field optimization (MFO) technique compatible with the application of a Monte Carlo (MC) algorithm and to evaluate its robustness. Nine patients (three brain, five head-and-neck, one spine) underwent proton treatment generated by a novel robust MFO technique. A hybrid (hMFO) approach was implemented, planning dose coverage on isotropic PTV compensating for setup errors, whereas range calibration uncertainties are incorporated into PTV robust optimization process. hMFO was compared with single-field optimization (SFO) and full robust multi-field optimization (fMFO), both on the nominal plan and the worst-case scenarios assessed by robustness analysis. The SFO and the fMFO plans were normalized to hMFO on CTV to obtain iso-D95 coverage, and then the organs at risk (OARs) doses were compared. On the same OARs, in the normalized nominal plans the potential impact of variable relative biological effectiveness (RBE) was investigated. hMFO reduces the number of scenarios computed for robust optimization (from twenty-one in fMFO to three), making it practicable with the application of a MC algorithm. After normalizing on D95 CTV coverage, nominal hMFO plans were superior compared to SFO in terms of OARs sparing (p   <  0.01), without significant differences compared to fMFO. The improvement in OAR sparing with hMFO with respect to SFO was preserved in worst-case scenarios (p   <  0.01), confirming that hMFO is as robust as SFO to physical uncertainties, with no significant differences when compared to the worst case scenarios obtained by fMFO. The dose increase on OARs due to variable RBE was comparable to the increase due to physical uncertainties (i.e. 4-5 Gy(RBE)), but without significant differences between these techniques. hMFO allows improving plan quality with respect to SFO, with no significant differences with fMFO and without affecting robustness to setup, range and RBE uncertainties, making clinically feasible the application of MC-based robust optimization.

Accurate proton treatment planning for pencil beam crossing titanium fixation implants

PURPOSE

To present a planning strategy for proton pencil-beam scanning when titanium implants need to be crossed by the beam.

METHODS

We addressed three issues: the implementation of a CT calibration curve to assign to titanium the correct stopping power; the effect of artefacts on CT images and their reduction by a dedicated algorithm; the differences in dose computation depending on the dose engine, pencil-beam vs Monte-Carlo algorithms. We performed measurement tests on a simple cylinder phantom and on a real implant. These phantoms were irradiated with three geometries (single spots, uniform mono-energetic layer and uniform box), measuring the exit dose either by radio-chromic film or multi-layer ionization chamber. The procedure was then applied on two patients treated for chordoma.

RESULTS

We had to set in the calibration curve a mass density equal to 4.37 g/cm³ to saturated Hounsfield Units, in order to have the correct stopping power assigned to titanium in TPS. CT artefact reduction algorithm allowed a better reconstruction of the shape and size of the implant. Monte-Carlo resulted accurate in computing the dose distribution whereas the pencil-beam algorithm failed due to sharp density interfaces between titanium and the surrounding material. Finally, the treatment plans obtained on two patients showed the impact of the dose engine algorithm, with 10-20% differences between pencil-beam and Monte-Carlo in small regions distally to the titanium screws.

CONCLUSION

The described combination of CT calibration, artefacts reduction and Monte-Carlo computation provides a reliable methodology to compute dose in patients with titanium implants.

FLUKA particle therapy tool for Monte Carlo independent calculation of scanned proton and carbon ion beam therapy

While Monte Carlo (MC) codes are considered as the gold standard for dosimetric calculations, the availability of user friendly MC codes suited for particle therapy is limited. Based on the FLUKA MC code and its graphical user interface (GUI) Flair, we developed an easy-to-use tool which enables simple and reliable simulations for particle therapy. In this paper we provide an overview of functionalities of the tool and with the presented clinical, proton and carbon ion therapy examples we demonstrate its reliability and the usability in the clinical environment and show its flexibility for research purposes. The first, easy-to-use FLUKA MC platform for particle therapy with GUI functionalities allows a user with a minimal effort and reduced knowledge about MC details to apply MC at their facility and is expected to enhance the popularity of the MC for both research and clinical quality assurance and commissioning purposes.

Improvements in pencil beam scanning proton therapy dose calculation accuracy in brain tumor cases with a commercial Monte Carlo algorithm

A commercial Monte Carlo (MC) algorithm (RayStation version 6.0.024) for the treatment of brain tumors with pencil beam scanning (PBS) proton therapy is validated and compared via measurements and analytical calculations in clinically realistic scenarios. For the measurements a 2D ion chamber array detector (MatriXX PT) was placed underneath the following targets: (1) an anthropomorphic head phantom (with two different thicknesses) and (2) a biological sample (i.e. half a lamb's head). In addition, we compared the MC dose engine versus the RayStation pencil beam (PB) algorithm clinically implemented so far, in critical conditions such as superficial targets (i.e. in need of a range shifter (RS)), different air gaps, and gantry angles to simulate both orthogonal and tangential beam arrangements. For every plan the PB and MC dose calculations were compared to measurements using a gamma analysis metrics (3%, 3 mm). For the head phantom the gamma passing rate (GPR) was always  >96% and on average  >99% for the MC algorithm; the PB algorithm had a GPR of  ⩽90% for all the delivery configurations with a single slab (apart 95% GPR from the gantry of 0° and small air gap) and in the case of two slabs of the head phantom the GPR was  >95% only in the case of small air gaps for all three (0°, 45°, and 70°) simulated beam gantry angles. Overall the PB algorithm tends to overestimate the dose to the target (up to 25%) and underestimate the dose to the organ at risk (up to 30%). We found similar results (but a bit worse for the PB algorithm) for the two targets of the lamb's head where only two beam gantry angles were simulated. Our results suggest that in PBS proton therapy a range shifter (RS) needs to be used with caution when planning a treatment with an analytical algorithm due to potentially great discrepancies between the planned dose and the dose delivered to the patient, including in the case of brain tumors where this issue could be underestimated. Our results also suggest that a MC evaluation of the dose has to be performed every time the RS is used and, mostly, when it is used with large air gaps and beam directions tangential to the patient surface.

Evaluation of the water-equivalence of plastic materials in low- and high-energy clinical proton beams

The aim of this work was to evaluate the water-equivalence of new trial plastics designed specifically for light-ion beam dosimetry as well as commercially available plastics in clinical proton beams. The water-equivalence of materials was tested by computing a plastic-to-water conversion factor, [Formula: see text]. Trial materials were characterized experimentally in 60 MeV and 226 MeV un-modulated proton beams and the results were compared with Monte Carlo simulations using the FLUKA code. For the high-energy beam, a comparison between the trial plastics and various commercial plastics was also performed using FLUKA and Geant4 Monte Carlo codes. Experimental information was obtained from laterally integrated depth-dose ionization chamber measurements in water, with and without plastic slabs with variable thicknesses in front of the water phantom. Fluence correction factors, [Formula: see text], between water and various materials were also derived using the Monte Carlo method. For the 60 MeV proton beam, [Formula: see text] and [Formula: see text] factors were within 1% from unity for all trial plastics. For the 226 MeV proton beam, experimental [Formula: see text] values deviated from unity by a maximum of about 1% for the three trial plastics and experimental results showed no advantage regarding which of the plastics was the most equivalent to water. Different magnitudes of corrections were found between Geant4 and FLUKA for the various materials due mainly to the use of different nonelastic nuclear data. Nevertheless, for the 226 MeV proton beam, [Formula: see text] correction factors were within 2% from unity for all the materials. Considering the results from the two Monte Carlo codes, PMMA and trial plastic #3 had the smallest [Formula: see text] values, where maximum deviations from unity were 1%, however, PMMA range differed by 16% from that of water. Overall, [Formula: see text] factors were deviating more from unity than [Formula: see text] factors and could amount to a few percent for some materials.

Zero field PDD and TMR data for unflattened beams in conventional linacs: A tool for independent dose calculations

PURPOSE

To investigate the applicability of the formalism described in BJR supplement n.25 for Flattening Filter Free (FFF) beams in determining the zero-field tissue maximum ratio (TMR) for an independent calculation method of Percentage Depth Doses (PDDs) and relative dose factors (RDFs) at different experimental setups.

METHODS

Experimental PDDs for field size from 40×40cm² to 2×2cm² with Source Surface Distance (SSD) 100cm were acquired. The normalized peak scatter factor for each square field was obtained by fitting experimental RDFs in water and collimator factors (CFs) in air. Maximum log-likelihood methods were used to extract fit parameters in competing models and the Bayesian Information Criterion was used to select the best one. In different experimental setups additional RDFs and TPR₁₀²⁰s for field sizes other than reference field were measured and Monte Carlo simulations of PDDs at SSD 80cm were carried out to validate the results. PDD agreements were evaluated by gamma analysis.

RESULTS

The BJR formalism allowed to predict the PDDs obtained with MC within 2%/2mm at SSD 80cm from 100% down to 50% of the maximum dose. The agreement between experimental TPR₁₀²⁰s and RDFs values at SSD=90cm and BJR calculations were within 1% for field sizes greater than 5×5cm² while it was within 3% for fields down to 2×2cm².

CONCLUSIONS

BJR formalism can be used for FFF beams to predict PDD and RDF at different SSDs and can be used for independent MU calculations.

Investigation into the effects of high-Z nano materials in proton therapy

High-Z nano materials have been previously shown to increase the amount of dose deposition within the tumour due to an increase in secondary electrons. This study evaluates the effects of high-Z nano materials in combination with protons, and the impact of proton energy, nanoparticle material and concentration. These effects were studied in silico through Monte Carlo simulation and experimentally through a phantom study, with particular attention to macroscale changes to the Bragg peak in the presence of nanoparticles. Three nanoparticle materials were simulated (gold, silver and platinum) at three concentrations (0.01, 0.1 and 6.5 mg ml(-1)) at two clinical proton energies (60 and 226 MeV). Simulations were verified experimentally using Gafchromic film measurements of gold nanoparticles suspended in water at two available high concentrations (5.5 mg ml(-1) and 1.1 mg ml(-1)). A significant change to Bragg peak features was evident, where at 226 MeV and 6.5 mg ml(-1), simulations of gold showed a 4.7 mm longitudinal shift of the distal edge and experimentally at 5.5 mg ml(-1), a shift of 2.2 mm. Simulations showed this effect to be material dependent, where platinum having the highest physical density caused the greatest shift with increasing concentration. A dose enhancement of 6%  ±  0.05 and 5%  ±  0.15 (60 MeV and 226 MeV, respectively) was evident with gold at 6.5 mg ml(-1) to water alone, compared to the 21%  ±  0.53 observed experimentally as dose to film with 5.5 mg ml(-1) of gold nanoparticles suspended in water at 226 MeV. The introduction of nanoparticles has strong potential to enhance dose in proton therapy, however the changes to the Bragg peak distribution that occur with high concentrations need to be accounted for to ensure tumour coverage.

Characterization and validation of a Monte Carlo code for independent dose calculation in proton therapy treatments with pencil beam scanning

We propose a method of creating and validating a Monte Carlo (MC) model of a proton Pencil Beam Scanning (PBS) machine using only commissioning measurements and avoiding the nozzle modeling. Measurements with a scintillating screen coupled with a CCD camera, ionization chamber and a Faraday Cup were used to model the beam in TOPAS without using any machine parameter information but the virtual source distance from the isocenter. Then the model was validated on simple Spread Out Bragg Peaks (SOBP) delivered in water phantom and with six realistic clinical plans (many involving 3 or more fields) on an anthropomorphic phantom. In particular the behavior of the moveable Range Shifter (RS) feature was investigated and its modeling has been proposed. The gamma analysis (3%,3 mm) was used to compare MC, TPS (XiO-ELEKTA) and measured 2D dose distributions (using radiochromic film). The MC modeling proposed here shows good results in the validation phase, both for simple irradiation geometry (SOBP in water) and for modulated treatment fields (on anthropomorphic phantoms). In particular head lesions were investigated and both MC and TPS data were compared with measurements. Treatment plans with no RS always showed a very good agreement with both of them (γ-Passing Rate (PR)  >  95%). Treatment plans in which the RS was needed were also tested and validated. For these treatment plans MC results showed better agreement with measurements (γ-PR  >  93%) than the one coming from TPS (γ-PR  <  88%). This work shows how to simplify the MC modeling of a PBS machine for proton therapy treatments without accounting for any hardware components and proposes a more reliable RS modeling than the one implemented in our TPS. The validation process has shown how this code is a valid candidate for a completely independent treatment plan dose calculation algorithm. This makes the code an important future tool for the patient specific QA verification process.

[Allergic diseases caused by pesticides: 3 case reports]

Current pesticide exposure levels can induce a state of sensitization in exposed workers and only in exceptional situations cause severe intoxication. Three cases are reported--one suffering from urticaria/angio-edema, one from asthma and the other from oculo-rhinitis. Only one subject was atopical. The substances to which sensitization was found were Cynoxamil, Mancozeb, Thiophanate, Seccatutto for the first case, Paraquat for the second and Dodine for the third subject. The results suggest the need to extend allergological tests also to pesticides when workers, or subjects in any case exposed, show clinical signs with suspected immunological pathogenesis. Furthermore, in order to reduce the number of new cases of pesticide allergy, a more severe control also of subjects exposed to low doses is proposed.