The sensitivity of radiobiological models in carbon ion radiotherapy (CIRT) and its consequences on the clinical treatment plan: Differences between LEM and MKM models.

PURPOSE: Carbon ion radiotherapy (CIRT) relies on relative biological effectiveness (RBE)-weighted dose calculations. Japanese clinics predominantly use the microdosimetric kinetic model (MKM), while European centers utilize the local effect model (LEM). Despite both models estimating RBE-distributions in tissue, their physical and mathematical assumptions differ, leading to significant disparities in RBE-weighted doses. Several European clinics adopted Japanese treatment schedules, necessitating adjustments in dose prescriptions and organ at risk (OAR) constraints. In the context of these two clinically used standards for RBE-weighted dose estimation, the objective of this study was to highlight specific scenarios for which the translations between models diverge, as shortcomings between them can influence clinical decisions. METHODS: Our aim was to discuss planning strategies minimizing those discrepancies, ultimately striving for more accurate and robust treatments. Evaluations were conducted in a virtual water phantom and patient CT-geometry, optimizing LEM RBE-weighted dose first and recomputing MKM thereafter. Dose-averaged linear energy transfer (LETd) distributions were also assessed. RESULTS: Results demonstrate how various parameters influence LEM/MKM translation. Similar LEM-dose distributions lead to markedly different MKM-dose distributions and variations in LETd. Generally, a homogeneous LEM RBE-weighted dose aligns with lower MKM values in most of the target volume. Nevertheless, paradoxical MKM hotspots may emerge (at the end of the range), potentially influencing clinical outcomes. Therefore, translation between models requires great caution. CONCLUSIONS: Understanding the relationship between these two clinical standards enables combining European and Japanese based experiences. The implementation of optimal planning strategies ensures the safety and acceptability of the clinical plan for both models and therefore enhances plan robustness from the RBE-weighted dose and LETd distribution point of view. This study emphasizes the importance of optimal planning strategies and the need for comprehensive CIRT plan quality assessment tools. In situations where simultaneous LEM and MKM computation capabilities are lacking, it can provide guidance in plan design, ultimately contributing to enhanced CIRT outcomes.

Commissioning and clinical implementation of an independent dose calculation system for scanned proton beams.

PURPOSE: Experimental patient-specific QA (PSQA) is a time and resource-intensive process, with a poor sensitivity in detecting errors. Radiation therapy facilities aim to substitute it by means of independent dose calculation (IDC) in combination with a comprehensive beam delivery QA program. This paper reports on the commissioning of the IDC software tool myQA iON (IBA Dosimetry) for proton therapy and its clinical implementation at the MedAustron Ion Therapy Center. METHODS: The IDC commissioning work included the validation of the beam model, the implementation and validation of clinical CT protocols, and the evaluation of patient treatment data. Dose difference maps, gamma index distributions, and pass rates (GPR) have been reviewed. The performance of the IDC tool has been assessed and clinical workflows, simulation settings, and GPR tolerances have been defined. RESULTS: Beam model validation showed agreement of ranges within ± 0.2 mm, Bragg-Peak widths within ± 0.1 mm, and spot sizes at various air gaps within ± 5% compared to physical measurements. Simulated dose in 2D reference fields deviated by -0.3% ± 0.5%, while 3D dose distributions differed by 1.8% on average to measurements. Validation of the CT calibration resulted in systematic differences of 2.0% between IDC and experimental data for tissue like samples. GPRs of 99.4 ± 0.6% were found for head, head and neck, and pediatric CT protocols on a 2%/2 mm gamma criterion. GPRs for the adult abdomen protocol were at 98.9% on average with 3%/3 mm. Root causes of GPR outliers, for example, implants were identified and evaluated. CONCLUSION: IDC has been successfully commissioned and integrated into the MedAustron clinical workflow for protons in 2021. IDC has been stepwise and safely substituting experimental PSQA since February 2021. The initial reduction of proton experimental PSQA was about 25% and reached up to 90% after 1 year.

The trends and significance of SSTR PET/CT added to MRI in follow-up imaging of low-grade meningioma treated with fractionated proton therapy.

PURPOSE: Overexpression of the somatostatin receptor (SSTR) has led to adoption of SSTR PET/CT for diagnosis and radiotherapy planning in meningioma, but data on SSTR expression during follow-up remain scarce. We investigated PET/CT quantifiers of SSTR tracers in WHO grade I meningioma following fractionated proton beam therapy (PBT) compared to standard response assessment with MRI. METHODS: Twenty-two patients diagnosed with low-grade meningioma treated by PBT were included. Follow-up included clinical visits, MRI, and [68Ga]Ga-DOTATOC PET/CT scans. Radiologic tumor response, MRI and PET volume (VMRI and VPET), maximum and mean standardied uptake value (SUVmax/SUVmean), total lesion activity (TLA), and heterogeneity index (HI) were evaluated. RESULTS: Median follow-up was 35.3 months (range: 6.4-47.9). Nineteen patients (86.4%, p = 0.0009) showed a decrease of SUVmax between baseline and first follow-up PET/CT (median: -24%, range: -53% to +89%) and in 81.8% of all cases, the SUVmax, SUVmean, and TLA at last follow-up were eventually lower than at baseline (p = 0.0043). Ambiguous trends without significance between the timepoints analyzed were observed for VPET. HI increased between baseline and last follow-up in 75% of cases (p = 0.024). All patients remained radiologically and clinically stable. Median VMRI decreased by -9.3% (range 0-32.5%, p < 0.0001) between baseline and last follow-up. CONCLUSION: PET/CT in follow-up of irradiated meningioma showed an early trend towards decreased binding of SSTR-specific tracers following radiation and MRI demonstrated consistently stable or decreasing tumor volume. Translational research is needed to clarify the underlying biology of the subsequent increase in SSTR PET quantifiers.

Implementation of Sphinx/Lynx as daily QA equipment for scanned proton and carbon ion beams.

PURPOSE: Reporting on the first implementation of a proton dedicated commercial device (IBA Sphinx/Lynx) for daily Quality Assurance (QA) of scanned proton and carbon ion beams. METHODS: Daily QA trendlines over more than 3 years for protons and more than 2 years for carbon ions have been acquired. Key daily QA parameters were reviewed, namely the spot size and position, beam range, Bragg peak width, coincidence (between beam and imaging system isocenters), homogeneity and dose. RESULTS: The performance of the QA equipment for protons and carbon ions was evaluated. Daily QA trendlines allowed us to detect machine performance drifts and changes. The definition of tolerances and action levels is provided and compared with levels used in the literature. CONCLUSION: The device has been successfully implemented for routine daily QA activities in a dual particle therapy facility for more than 2 years. It improved the efficiency of daily QA and provides a comprehensive QA process.

Harmonization of proton treatment planning for head and neck cancer using pencil beam scanning: first report of the IPACS collaboration group.

Background and purpose: A collaborative network between proton therapy (PT) centres in Trento in Italy, Poland, Austria, Czech Republic and Sweden (IPACS) was founded to implement trials and harmonize PT. This is the first report of IPACS with the aim to show the level of harmonization that can be achieved for proton therapy planning of head and neck (sino-nasal) cancer.Methods: CT-data sets of five patients were included. During several face-to-face and online meetings, a common treatment planning protocol was developed. Each centre used its own treatment planning system (TPS) and planning approach with some restrictions specified in the treatment planning protocol. In addition, volumetric modulated arc therapy (VMAT) photon plans were created.Results: For CTV1, the average Dmedian was 59.3 ± 2.4 Gy(RBE) for protons and 58.8 ± 2.0 Gy(RBE) for VMAT (aim was 56 Gy(RBE)). For CTV2, the average Dmedian was 71.2 ± 1.0 Gy(RBE) for protons and 70.6 ± 0.4 Gy(RBE) for VMAT (aim was 70 Gy(RBE)). The average D2% for the spinal cord was 25.1 ± 8.5 Gy(RBE) for protons and 47.6 ± 1.4 Gy(RBE) for VMAT. The average D2% for chiasm was 46.5 ± 4.4 Gy(RBE) for protons and 50.8 ± 1.4 Gy(RBE) for VMAT, respectively. Robust evaluation was performed and showed the least robust plans for plans with a low number of beams.Discussion: In conclusion, several influences on harmonization were identified: adherence/interpretation to/of the protocol, available technology, experience in treatment planning and use of different beam arrangements. In future, all OARs that should be included in the optimization need to be specified in order to further harmonize treatment planning.

Individually optimized stereotactic radiotherapy for pancreatic head tumors: A planning feasibility study.

AIM: Aim of this study was to perform a planning feasibility analysis of a 3-level dose prescription using an IMRT-SIB technique. BACKGROUND: Radiation therapy of locally advanced pancreatic cancer should administer a minimum dose to the duodenum and a very high dose to the vascular infiltration areas to improve the possibility of a radical resection. MATERIALS AND METHODS: Fifteen patients with pancreatic head adenocarcinoma and vascular involvement were included. The duodenal PTV (PTVd) was defined as the GTV overlapping the duodenal PRV. Vascular CTV (CTVv) was defined as the surface of contact or infiltration between the tumor and vessel plus a 5 mm margin. Vascular PTV (PTVv) was considered as the CTVv plus an anisotropic margin. The tumor PTV (PTVt) was defined as the GTV plus a margin including the PTVv and excluding the PTVd. The following doses were prescribed: 30 Gy (6 Gy/fraction) to PTVd, 37.5 Gy (7.5 Gy/fraction) to PTVt, and 45 Gy (9 Gy/fraction) to PTVv, respectively. Treatment was planned with an IMRT technique. RESULTS: The primary end-point (PTVv Dmean > 90%) was achieved in all patients. PTVv D98% > 90% was achieved in 6 patients (40%). OARs constraints were achieved in all patients. CONCLUSIONS: Although the PTVv D95% > 95% objective was achieved only in 40% of patients, the study showed that in 100% of patients it was possible to administer a strongly differentiated mean/median dose. Prospective trials based on clinical application of this strategy seem to be justified in selected patients without overlap between PTVd and PTVv.

ART for head and neck patients: On the difference between VMAT and IMPT.

UNLABELLED: Anatomical changes in the head-and-neck (H&N) region during the course of treatment can cause deteriorated dose distributions. Different replanning strategies were investigated for volumetric modulated arc therapy (VMAT) and intensity-modulated proton therapy (IMPT). MATERIAL AND METHODS: For six H&N patients two repeated computed tomography (CT) and magnetic resonance (MR) (CT1/MR1 at week 2 and CT2/MR2 at week 4) scans were acquired additionally to the initial planning CT/MR. Organs-at-risk (OARs) and three targets (CTV70Gy, CTV63Gy, CTV56Gy) were delineated on MRs and transferred to respective CT data set. Simultaneously integrated boost plans were created using VMAT (two arcs) and IMPT (four beams). To assess the need of replanning the initial VMAT and IMPT plans were recalculated on repeated CTs. Furthermore, VMAT and IMPT plans were replanned on the repeated CTs. A Demon algorithm was used for deformable registration of the repeated CTs with the initial CT and utilized for dose accumulation. Total dose estimations were performed to compare ART versus standard treatment strategies. RESULTS: Dosimetric evaluation of recalculated plans on CT1 and CT2 showed increasing OAR doses for both, VMAT and IMPT. The target coverage of recalculated VMAT plans was considered acceptable in three cases, while for all IMPT plans it dropped. Adaptation of the treatment reduced D2% for brainstem by 6.7 Gy for VMAT and by 8 Gy for IMPT, for particular patients. These D2% reductions were reaching 9 Gy and 14 Gy for the spinal cord. ART improved target dose homogeneity, especially for protons, i.e. D2% decreased by up to 8 Gy while D98% increased by 1.2 Gy. CONCLUSION: ART showed benefits for both modalities. However, as IMPT is more conformal, the magnitude of dosimetric changes was more pronounced compared to VMAT. Large anatomic variations had a severe impact on treatment plan quality for both VMAT and IMPT. ART is justified in those cases irrespective of treatment modalities.

Sacral-Nerve-Sparing Planning Strategy in Pelvic Sarcomas/Chordomas Treated with Carbon-Ion Radiotherapy.

To minimize radiation-induced lumbosacral neuropathy (RILSN), we employed sacral-nerve-sparing optimized carbon-ion therapy strategy (SNSo-CIRT) in treating 35 patients with pelvic sarcomas/chordomas. Plans were optimized using Local Effect Model-I (LEM-I), prescribed DRBE|LEM-I|D50% (median dose to HD-PTV) = 73.6 (70.4-76.8) Gy (RBE)/16 fractions. Sacral nerves were contoured between L5-S3 levels. DRBE|LEM-I to 5% of sacral nerves-to-spare (outside HD-CTV) (DRBE|LEM-I|D5%) were restricted to <69 Gy (RBE). The median follow-up was 25 months (range of 2-53). Three patients (9%) developed late RILSN (≥G3) after an average period of 8 months post-CIRT. The RILSN-free survival at 2 years was 91% (CI, 81-100). With SNSo-CIRT, DRBE|LEM-I|D5% for sacral nerves-to-spare = 66.9 ± 1.9 Gy (RBE), maintaining DRBE|LEM-I to 98% of HD-CTV (DRBE|LEM-I|D98%) = 70 ± 3.6 Gy (RBE). Two-year OS and LC were 100% and 93% (CI, 84-100), respectively. LETd and DRBE with modified-microdosimetric kinetic model (mMKM) were recomputed retrospectively. DRBE|LEM-I and DRBE|mMKM were similar, but DRBE-filtered-LETd was higher in sacral nerves-to-spare in patients with RILSN than those without. At DRBE|LEM-I cutoff = 64 Gy (RBE), 2-year RILSN-free survival was 100% in patients with <12% of sacral nerves-to-spare voxels receiving LETd > 55 keV/µm than 75% (CI, 54-100) in those with ≥12% of voxels (p < 0.05). DRBE-filtered-LETd holds promise for the SNSo-CIRT strategy but requires longer follow-up for validation.

Investigation on the physical dose filtered by linear energy transfer for treatment plan evaluation in carbon ion therapy.

BACKGROUND: Large tumor size has been reported as a predicting factor for inferior clinical outcome in carbon ion radiotherapy (CIRT). Besides the clinical factors accompanied with such tumors, larger tumors receive typically more low linear energy transfer (LET) contributions than small ones which may be the underlying physical cause. Although dose averaged LET is often used as a single parameter descriptor to quantify the beam quality, there is no evidence that this parameter is the optimal clinical predictor for the complex mixed radiation fields in CIRT. PURPOSE: Purpose of this study was to investigate on a novel dosimetric quantity, namely high-LET-dose ( D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the physical dose filtered based on an LET threshold) as a single parameter estimator to differentiate between carbon ion treatment plans (cTP) with a small and large tumor volume. METHODS: Ten cTPs with a planning target volume, PTV ≥ 500 cm 3 $\mathrm{PTV}\ge {500}\,{{\rm cm}^{3}}$ (large) and nine with a PTV < 500 cm 3 $\mathrm{PTV}<{500}\,{{\rm cm}^{3}}$ (small) were selected for this study. To find a reasonable LET threshold ( L thr $\textrm {L}_{\textrm {thr}}$ ) that results in a significant difference in terms of D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the voxel based normalized high-LET-dose ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) distribution in the clinical target volume (CTV) was studied on a subset (12 out of 19 cTPs) for 18 LET thresholds, using standard distribution descriptors (mean, variance and skewness). The classical dose volume histogram concept was used to evaluate the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ distributions within the target of all 19 cTPs at the before determined L thr $\textrm {L}_{\textrm {thr}}$ . Statistical significance of the difference between the two groups in terms of mean D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ volume histogram parameters was evaluated by means of (two-sided) t-test or Mann-Whitney-U-test. In addition, the minimum target coverage at the above determined L thr $\textrm {L}_{\textrm {thr}}$ was compared and validated against three other thresholds to verify its potential in differentiation between small and large volume tumors. RESULTS: An L thr $\textrm {L}_{\textrm {thr}}$ of approximately 30 keV / µ m ${30}\,{\rm keV/}\umu {\rm m}$ was found to be a reasonable threshold to classify the two groups. At this threshold, the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ were significantly larger ( p < 0.05 $p<0.05$ ) in small CTVs. For the small tumor group, the near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ (and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) in the CTV were in average 9.3 ± 1.5 Gy $9.3\pm {1.5}\,{\rm Gy}$ (0.31 ± 0.08) and 13.6 ± 1.6 Gy $13.6\pm {1.6}\,{\rm Gy}$ (0.46 ± 0.06), respectively. For the large tumors, these parameters were 6.6 ± 0.2 Gy $6.6\pm {0.2}\,{\rm Gy}$ (0.20 ± 0.01) and 8.6 ± 0.4 Gy $8.6\pm {0.4}\,{\rm Gy}$ (0.28 ± 0.02). The difference between the two groups in terms of mean near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) was 2.7 Gy (11%) and 5.0 Gy (18%), respectively. CONCLUSIONS: The feasibility of high-LET-dose based evaluation was shown in this study where a lower D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ was found in cTPs with a large tumor size. Further investigation is needed to draw clinical conclusions. The proposed methodology in this work can be utilized for future high-LET-dose based studies.

Validation of a deep-learning segmentation model for adult and pediatric head and neck radiotherapy in different patient positions.

Background and purpose: Autocontouring for radiotherapy has the potential to significantly save time and reduce interobserver variability. We aimed to assess the performance of a commercial autocontouring model for head and neck (H&N) patients in eight orientations relevant to particle therapy with fixed beam lines, focusing on validation and implementation for routine clinical use. Materials and methods: Autocontouring was performed on sixteen organs at risk (OARs) for 98 adult and pediatric patients with 137 H&N CT scans in eight orientations. A geometric comparison of the autocontours and manual segmentations was performed using the Hausdorff Distance 95th percentile, Dice Similarity Coefficient (DSC) and surface DSC and compared to interobserver variability where available. Additional qualitative scoring and dose-volume-histogram (DVH) parameters analyses were performed for twenty patients in two positions, consisting of scoring on a 0-3 scale based on clinical usability and comparing the mean (Dmean) and near-maximum (D2%) dose, respectively. Results: For the geometric analysis, the model performance in head-first-supine straight and hyperextended orientations was in the same range as the interobserver variability. HD95, DSC and surface DSC was heterogeneous in other orientations. No significant geometric differences were found between pediatric and adult autocontours. The qualitative scoring yielded a median score of ≥ 2 for 13/16 OARs while 7/32 DVH parameters were significantly different. Conclusions: For head-first-supine straight and hyperextended scans, we found that 13/16 OAR autocontours were suited for use in daily clinical practice and subsequently implemented. Further development is needed for other patient orientations before implementation.

Dose averaged linear energy transfer optimization for large sacral chordomas in carbon ion therapy.

BACKGROUND: Carbon ion beams are well accepted as densely ionizing radiation with a high linear energy transfer (LET). However, the current clinical practice does not fully exploit the highest possible dose-averaged LET (LETd) and, consequently, the biological potential in the target. This aspect becomes worse in larger tumors for which inferior clinical outcomes and corresponding lower LETd was reported. PURPOSE: The vicinity to critical organs in general and the inferior overall survival reported for larger sacral chordomas treated with carbon ion radiotherapy (CIRT), makes the treatment of such tumors challenging. In this work it was aimed to increase the LETd in large volume tumors while maintaining the relative biological effectiveness (RBE)-weighted dose, utilizing the LETd optimization functions of a commercial treatment planning system (TPS). METHODS: Ten reference sequential boost carbon ion treatment plans, designed to mimic clinical plans for large sacral chordoma tumors, were generated. High dose clinical target volumes (CTV-HD) larger than 250 cm 3 $250 \,{\rm cm}^{3}$ were considered as large targets. The total RBE-weighted median dose prescription with the local effect model (LEM) was D RBE , 50 % = 73.6 Gy $\textrm {D}_{\rm RBE, 50\%}=73.6 \,{\rm Gy}$ in 16 fractions (nine to low dose and seven to high dose planning target volume). No LETd optimization was performed in the reference plans, while LETd optimized plans used the minimum LETd (Lmin) optimization function in RayStation 2023B. Three different Lmin values were investigated and specified for the seven boost fractions: L min = 60 keV / µ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ , L min = 80 keV / µ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ and L min = 100 keV / µ m $\textrm {L}_{\rm min}=100 \,{\rm keV}/{\umu }{\rm m}$ . To compare the LETd optimized against reference plans, LETd and RBE-weighted dose based goals similar to and less strict than clinical ones were specified for the target. The goals for the organs at risk (OAR) remained unchanged. Robustness evaluation was studied for eight scenarios ( ± 3.5 % $\pm 3.5\%$ range uncertainty and ± 3 mm $\pm 3 \,{\rm mm}$ setup uncertainty along the main three axes). RESULTS: The optimization method with L min = 60 keV / µ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ resulted in an optimal LETd distribution with an average increase of LET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (and LET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) in the CTV-HD by 8.9 ± 1.5 keV / µ m $8.9\pm 1.5 \,{\rm keV}/{\umu }{\rm m}$ ( 27 % $27\%$ ) (and 6.9 ± 1.3 keV / µ m $6.9\pm 1.3 \,{\rm keV}/{\umu }{\rm m}$ ( 17 % $17\%$ )), without significant difference in the RBE-weighted dose. By allowing ± 5 % $\pm 5\%$ over- and under-dosage in the target, the LET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (and LET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) can be increased by 11.3 ± 1.2 keV / µ m $11.3\pm 1.2 \,{\rm keV}/{\umu }{\rm m}$ ( 34 % $34\%$ ) (and 11.7 ± 3.4 keV / µ m $11.7\pm 3.4 \,{\rm keV}/{\umu }{\rm m}$ ( 29 % $29\%$ )), using the optimization parameters L min = 80 keV / µ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ . The pass rate for the OAR goals in the LETd optimized plans was in the same level as the reference plans. LETd optimization lead to less robust plans compared to reference plans. CONCLUSIONS: Compared to conventionally optimized treatment plans, the LETd in the target was increased while maintaining the RBE-weighted dose using TPS LETd optimization functionalities. Regularly assessing RBE-weighted dose robustness and acquiring more in-room images remain crucial and inevitable aspects during treatment.

Effects of nuclear interaction corrections and trichrome fragment spectra modelling on dose and linear energy transfer distributions in carbon ion radiotherapy.

Background and Purpose: Nuclear interaction correction (NIC) and trichrome fragment spectra modelling improve relative biological effectiveness-weighted dose (DRBE) and dose-averaged linear energy transfer (LETd) calculation for carbon ions. The effect of those novel approaches on the clinical dose and LET distributions was investigated. Materials and Methods: The effect of the NIC and trichrome algorithm was assessed, creating single beam plans for a virtual water phantom with standard settings and NIC + trichrome corrections. Reference DRBE and LETd distributions were simulated using FLUKA version 2021.2.9. Thirty clinically applied scanned carbon ion treatment plans were recalculated applying NIC, trichrome and NIC + trichrome corrections, using the LEM low dose approximation and compared to clinical plans (base RS). Four treatment sites were analysed: six prostate adenocarcinoma, ten head and neck, nine locally advanced pancreatic adenocarcinoma and five sacral chordoma. The FLUKA and clinical plans were compared in terms of DRBE deviations for D98%, D50%, D2% for the clinical target volume (CTV) and D50% in ring-like dose regions retrieved from isodose curves in base RS plans. Additionally, region-based median LETd deviations and global gamma parameters were evaluated. Results: Dose deviations comparing base RS and evaluation plans were within ± 1% supported by γ-pass rates over 97% for all cases. No significant LETd deviations were reported in the CTV, but significant median LETd deviations were up to 80% for very low dose regions. Conclusion: Our results showed improved accuracy of the predicted DRBE and LETd. Considering clinically relevant constraints, no significant modifications of clinical protocols are expected with the introduction of NIC + trichrome.

Prospective Analysis of Radiation-Induced Contrast Enhancement and Health-Related Quality of Life After Proton Therapy for Central Nervous System and Skull Base Tumors.

PURPOSE: Intracerebral radiation-induced contrast enhancement (RICE) can occur after photon as well as proton beam therapy (PBT). This study evaluated the incidence, characteristics, and risk factors of RICE after PBT delivered to, or in direct proximity to, the brain and its effect on health-related quality of life (HRQoL). METHODS AND MATERIALS: Four hundred twenty-one patients treated with pencil beam scanning PBT between 2017 and 2021 were included. Follow-up included clinical evaluation and contrast-enhanced magnetic resonance imaging at 3, 6, and 12 months after treatment completion and annually thereafter. RICE was graded according to Common Terminology Criteria for Adverse Events version 4, and HRQoL parameters were assessed via European Organisation for Research and Treatment of Cancer Quality of Life Questionnaire (EORTC QLQ)-C30 questionnaires. RESULTS: The median follow-up was 24 months (range, 6-54), and median dose to 1% relative volume of noninvolved central nervous system (D1%CNS) was 54.3 Gy relative biologic effectiveness (RBE; range, 30-76 Gy RBE). The cumulative RICE incidence was 15% (n = 63), of which 10.5% (n = 44) were grade 1, 3.1% (n = 13) were grade 2, and 1.4% (n = 6) were grade 3. No grade 4 or 5 events were observed. Twenty-six of 63 RICE (41.3%) had resolved at the latest follow-up. The median onset after PBT and duration of RICE in patients in whom the lesions resolved were 11.8 and 9.0 months, respectively. On multivariable analysis, D1%CNS > 57.6 Gy RBE, previous in-field radiation, and diabetes mellitus were identified as significant risk factors for RICE development. Previous radiation was the only factor influencing the risk of symptomatic RICE. After PBT, general HRQoL parameters were not compromised. In a matched cohort analysis of 54/50 patients with and without RICE, no differences in global health score or functional and symptom scales were seen. CONCLUSIONS: The overall incidence of clinically relevant RICE after PBT is very low and has no significant negative effect on long-term patient QoL.

Patterns of practice of image guided particle therapy for cranio-spinal irradiation: A site specific multi-institutional survey of European Particle Therapy Network.

PURPOSE: To investigate the current practice patterns in image-guided particle therapy (IGPT) for cranio-spinal irradiation (CSI). METHODS: A multi-institutional survey was distributed to European particle therapy centres to analyse all aspects of IGPT. Based on the survey results, a Delphi consensus analysis was developed to define minimum requirements and optimal workflow for clinical practice. The centres participating in the institutional survey were invited to join the Delphi process. RESULTS: Eleven centres participated in the survey. Imaging for treatment planning was rather similar among the centres with Computed Tomography (CT) being the main modality. For positioning verification, 2D IGPT was more commonly used than 3D IGPT. Two centres performed routinely imaging for plan adaptation, by the rest ad hoc. Eight centres participated in the Delphi consensus analysis. The full consensus was reached on the use of CT imaging without contrast for treatment planning and the role of magnetic resonance imaging (MRI) in target and organs-at-risk delineation. There was an agreement on the necessity to perform patient position verification and correction before each isocentre. The most important outcome was the clear need for standardization and harmonization of the workflow. CONCLUSION: There were differences in CSI IGPT clinical practice among the European particle therapy centres. Moreover, the optimal workflow as identified by experts was not yet reached. There is a strong need for consensus guidelines. The state-of-the-art imaging technology and protocols need to be implemented into clinical practice to improve the quality of IGPT for CSI.

Robustness of IMPT treatment plans with respect to inter-fractional set-up uncertainties: impact of various beam arrangements for cranial targets.

UNLABELLED: In the current study IMPT plan robustness was evaluated with respect to inter-fractional patient positioning for various beam arrangements and two tumor indications in the cranial region. MATERIAL AND METHODS: For 14 patients suffering from tumors in the cranial region [skull base (SB; n = 7) and paranasal sinus (PS; n = 7)] the CTV and OARs were delineated. A safety margin of 3 mm was applied to the CTV. A prescribed dose of 2 GyE was planned via three beam arrangements (α, ß, γ). Beam arrangement α consisted of lateral opposed fields for both tumor groups while beam arrangement ß was optimized according to respective tumor and OAR locations, using two beams only. Beam arrangement γ applied four beams in the SB group and three beams in the PS group. Dose distributions were recalculated subjected to virtual patient translations along the major anatomical axes. The following dosimetric indices were evaluated and compared to original plans: target coverage (TC), target dose homogeneity (HI), CTV median and average dose (D(median), D(mean)). For OARs near maximum dose and average dose (D2%, D(mean)) were evaluated. RESULTS: Dose distributions were distorted after introducing shifts. In the SB group, TC and HI were significantly different for caudal, cranial and anterior shifts for all beam arrangements. For PS patients, all but right shifts differed significantly from the original plans for all beam arrangements, although clinical relevance was not reached for arrangement γ (ΔTC < 1.5%). In general, beam arrangement γ exhibited the least spread of data regarding target indices and was consequently considered the most robust. Dosimetric parameters regarding the brainstem were mostly influenced by shifts along the anterio-posterior axis. CONCLUSION: For cranial IMPT, set-up uncertainties may lead to pronounced deterioration of dose distributions. According to our investigations, multi-beam arrangements were dosimetrically more robust and hence preferable over two beam arrangements.

Real-time 2D/3D registration using kV-MV image pairs for tumor motion tracking in image guided radiotherapy.

Intra-fractional respiratory motion during radiotherapy leads to a larger planning target volume (PTV). Real-time tumor motion tracking by two-dimensional (2D)/3D registration using on-board kilo-voltage (kV) imaging can allow for a reduction of the PTV though motion along the imaging beam axis cannot be resolved using only one projection image. We present a retrospective patient study investigating the impact of paired portal mega-voltage (MV) and kV images on registration accuracy. Material and methods. We used data from 10 patients suffering from non-small cell lung cancer (NSCLC) undergoing stereotactic body radiation therapy (SBRT) lung treatment. For each patient we acquired a planning computed tomography (CT) and sequences of kV and MV images during treatment. We compared the accuracy of motion tracking in six degrees-of-freedom (DOF) using the anterior-posterior (AP) kV sequence or the sequence of kV-MV image pairs. Results. Motion along cranial-caudal direction could accurately be extracted when using only the kV sequence but in AP direction we obtained large errors. When using kV-MV pairs, the average error was reduced from 2.9 mm to 1.5 mm and the motion along AP was successfully extracted. Mean registration time was 188 ms. Conclusion. Our evaluation shows that using kV-MV image pairs leads to improved motion extraction in six DOF and is suitable for real-time tumor motion tracking with a conventional LINAC.

State-of-the-art and potential of experimental microdosimetry in ion-beam therapy.

In radiotherapy, radiation-quality should be an expression of the biological and physical characteristics of ionizing radiation such as spatial distribution of ionization or energy deposition. Linear energy transfer (LET) and lineal energy (y) are two descriptors used to quantify the radiation quality. These two quantities are connected and exhibit similar features. In ion-beam therapy (IBT), lineal energy can be measured with microdosimeters, which are specifically designed to cope with the high fluence of particles in clinical beams, while the quantification of LET is generally based on calculations. In pre-clinical studies, microdosimetric spectra are used for the indirect determination of relative biological effectiveness (RBE), e.g., using the microdosimetric kinetic model (MKM) or biophysical response functions. In this context it is important to consider saturation effects, which occur when the highest values of y become less biologically relevant compared to the relative contribution they make to the physical dose. Recent clinical data suggests that local tumor control and normal tissue effects can be linked to macroscopic and microscopic dosimetry parameters. In particular, positive clinical outcomes have been correlated to the highest LET values in the density distribution, and there is no evident link to the saturation discussed above. A systematic collection of microdosimetric information in combination with clinical data in retrospective studies may clarify the role of radiation quality at the highest LET. In the clinical setting, microdosimetry is not widely used yet, despite its potential to be linked with LET by experimentally-determined y values. Through this connection, both play an important role in complex therapy techniques such as intensity modulated particle therapy (IMPT), LET-painting and multi-ion optimization. This review summarizes the current state of microdosimetry for IBT and its potential, as well as research and development needed to make experimental microdosimetry a mature procedure in a clinical context.

Patient Breathing Motion and Delivery Specifics Influencing the Robustness of a Proton Pancreas Irradiation.

Motion compensation strategies in particle therapy depend on the anatomy, motion amplitude and underlying beam delivery technology. This retrospective study on pancreas patients with small moving tumours analysed existing treatment concepts and serves as a basis for future treatment strategies for patients with larger motion amplitudes as well as the transition towards carbon ion treatments. The dose distributions of 17 hypofractionated proton treatment plans were analysed using 4D dose tracking (4DDT). The recalculation of clinical treatment plans employing robust optimisation for mitigating different organ fillings was performed on phased-based 4D computed tomography (4DCT) data considering the accelerator (pulsed scanned pencil beams delivered by a synchrotron) and the breathing-time structure. The analysis confirmed the robustness of the included treatment plans concerning the interplay of beam and organ motion. The median deterioration of D50% (ΔD50%) for the clinical target volume (CTV) and the planning target volume (PTV) was below 2%, while the only outlier was observed for ΔD98% with -35.1%. The average gamma pass rate over all treatment plans (2%/ 2 mm) was 88.8% ± 8.3, while treatment plans for motion amplitudes larger than 1 mm performed worse. For organs at risk (OARs), the median ΔD2% was below 3%, but for single patients, essential changes, e.g., up to 160% for the stomach were observed. The hypofractionated proton treatment for pancreas patients based on robust treatment plan optimisation and 2 to 4 horizontal and vertical beams showed to be robust against intra-fractional movements up to 3.7 mm. It could be demonstrated that the patient's orientation did not influence the motion sensitivity. The identified outliers showed the need for continuous 4DDT calculations in clinical practice to identify patient cases with more significant deviations.

Possibilities and challenges when using synthetic computed tomography in an adaptive carbon-ion treatment workflow.

BACKGROUND AND PURPOSE: Anatomical surveillance during ion-beam therapy is the basis for an effective tumor treatment and optimal organ at risk (OAR) sparing. Synthetic computed tomography (sCT) based on magnetic resonance imaging (MRI) can replace the X-ray based planning CT (X-rayCT) in photon radiotherapy and improve the workflow efficiency without additional imaging dose. The extension to carbon-ion radiotherapy is highly challenging; complex patient positioning, unique anatomical situations, distinct horizontal and vertical beam incidence directions, and limited training data are only few problems. This study gives insight into the possibilities and challenges of using sCTs in carbon-ion therapy. MATERIALS AND METHODS: For head and neck patients immobilised with thermoplastic masks 30 clinically applied actively scanned carbon-ion treatment plans on 15 CTs comprising 60 beams were analyzed. Those treatment plans were re-calculated on MRI based sCTs which were created employing a 3D U-Net. Dose differences and carbon-ion spot displacements between sCT and X-rayCT were evaluated on a patient specific basis. RESULTS: Spot displacement analysis showed a peak displacement by 0.2 cm caused by the immobilisation mask not measurable with the MRI. 95.7% of all spot displacements were located within 1 cm. For the clinical target volume (CTV) the median D50% agreed within -0.2% (-1.3 to 1.4%), while the median D0.01cc differed up to 4.2% (-1.3 to 25.3%) comparing the dose distribution on the X-rayCT and the sCT. OAR deviations depended strongly on the position and the dose gradient. For three patients no deterioration of the OAR parameters was observed. Other patients showed large deteriorations, e.g. for one patient D2% of the chiasm differed by 28.1%. CONCLUSION: The usage of sCTs opens several new questions, concluding that we are not ready yet for an MR-only workflow in carbon-ion therapy, as envisaged in photon therapy. Although omitting the X-rayCT seems unfavourable in the case of carbon-ion therapy, an sCT could be advantageous for monitoring, re-planning, and adaptation.

Planning Strategy to Optimize the Dose-Averaged LET Distribution in Large Pelvic Sarcomas/Chordomas Treated with Carbon-Ion Radiotherapy.

To improve outcomes in large sarcomas/chordomas treated with CIRT, there has been recent interest in LET optimization. We evaluated 22 pelvic sarcoma/chordoma patients treated with CIRT [large: HD-CTV ≥ 250 cm3 (n = 9), small: HD-CTV < 250 cm3 (n = 13)], DRBE|LEM-I = 73.6 (70.4-73.6) Gy (RBE)/16 fractions, using the local effect model-I (LEM-I) optimization and modified-microdosimetric kinetic model (mMKM) recomputation. We observed that to improve high-LETd distribution in large tumors, at least 27 cm3 (low-LETd region) of HD-CTV should receive LETd of ≥33 keV/µm (p < 0.05). Hence, LETd optimization using 'distal patching' was explored in a treatment planning setting (not implemented clinically yet). Distal-patching structures were created to stop beams 1-2 cm beyond the HD-PTV-midplane. These plans were reoptimized and DRBE|LEM-I, DRBE|mMKM, and LETd were recomputed. Distal patching increased (a) LETd50% in HD-CTV (from 38 ± 3.4 keV/µm to 47 ± 8.1 keV/µm), (b) LETdmin in low-LETd regions of the HD-CTV (from 32 ± 2.3 keV/µm to 36.2 ± 3.6 keV/µm), (c) the GTV fraction receiving LETd of ≥50 keV/µm, (from <10% to >50%) and (d) the high-LETd component in the central region of the GTV, without significant compromise in DRBE distribution. However, distal patching is sensitive to setup/range uncertainties, and efforts to ascertain robustness are underway, before routine clinical implementation.