Next Article in Journal
Thymoquinone Radiosensitizes Human Colorectal Cancer Cells in 2D and 3D Culture Models
Next Article in Special Issue
Radiation-Induced Esophagitis in Non-Small-Cell Lung Cancer Patients: Voxel-Based Analysis and NTCP Modeling
Previous Article in Journal
Cytokine-Induced Senescence in the Tumor Microenvironment and Its Effects on Anti-Tumor Immune Responses
Previous Article in Special Issue
A Novel and Automated Approach to Classify Radiation Induced Lung Tissue Damage on CT Scans
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Substantial Sparing of Organs at Risk with Modern Proton Therapy in Lung Cancer, but Altered Breathing Patterns Can Jeopardize Target Coverage

by
Camilla Grindeland Boer
1,
Kristine Fjellanger
1,2,*,
Inger Marie Sandvik
1,
Maren Ugland
1,
Grete May Engeseth
1,3 and
Liv Bolstad Hysing
1,2
1
Department of Oncology and Medical Physics, Haukeland University Hospital, 5021 Bergen, Norway
2
Institute of Physics and Technology, University of Bergen, 5007 Bergen, Norway
3
Department of Clinical Science, University of Bergen, 5021 Bergen, Norway
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(6), 1365; https://doi.org/10.3390/cancers14061365
Submission received: 31 January 2022 / Revised: 1 March 2022 / Accepted: 4 March 2022 / Published: 8 March 2022

Abstract

:

Simple Summary

Treatment of locally advanced non-small cell lung cancer (LA-NSCLC) is a fine balance between toxicity and cure. Modern proton therapy might offer a more gentle radiation treatment compared to state-of-the-art photon radiotherapy, but is also more susceptible to the influence of breathing motion and anatomical changes. In this study, the influence of such uncertainties on treatment delivery was thoroughly investigated. Modern proton therapy did indeed show potential to reduce the risk of toxicity for the heart and lungs. This potential was maintained under the influence of anatomical and delivery uncertainties. However, changes in breathing motion jeopardized the target dose distribution in a subset of patients. We therefore recommend imaging at onset or early in treatment to recognize these patients and adapt the treatment.

Abstract

Enhancing treatment of locally advanced non-small cell lung cancer (LA-NSCLC) by using pencil beam scanning proton therapy (PBS-PT) is attractive, but little knowledge exists on the effects of uncertainties occurring between the planning (Plan) and the start of treatment (Start). In this prospective simulation study, we investigated the clinical potential for PBS-PT under the influence of such uncertainties. Imaging with 4DCT at Plan and Start was carried out for 15 patients that received state-of-the-art intensity-modulated radiotherapy (IMRT). Three PBS-PT plans were created per patient: 3D robust single-field uniform dose (SFUD), 3D robust intensity-modulated proton therapy (IMPT), and 4D robust IMPT (4DIMPT). These were exposed to setup and range uncertainties and breathing motion at Plan, and changes in breathing motion and anatomy at Start. Target coverage and dose-volume parameters relevant for toxicity were compared. The organ at risk sparing at Plan was greatest with IMPT, followed by 4DIMPT, SFUD and IMRT, and persisted at Start. All plans met the preset criteria for target robustness at Plan. At Start, three patients had a lack of CTV coverage with PBS-PT. In conclusion, the clinical potential for heart and lung toxicity reduction with PBS-PT was substantial and persistent. Altered breathing patterns between Plan and Start jeopardized target coverage for all PBS-PT techniques.

1. Introduction

State-of-the-art treatment for inoperable locally advanced non-small cell lung cancer (LA-NSCLC) is concurrent chemotherapy and intensity-modulated photon radiotherapy (IMRT) to a dose of 60 Gy. Still, 5-year survival rates for stage III disease are only around 30%, and side effects from treatment are common and potentially fatal, limiting the possibility for dose escalation with IMRT [1,2].
Proton therapy (PT) has advantageous depth dose characteristics with the potential to reduce side effects and facilitate dose escalation in LA-NSCLC patients [3,4,5,6,7]. Although phase II clinical trials have been promising [8,9,10], PT showed no advantage over IMRT in a randomized trial by Liao and colleagues [11,12]. These early clinical trials have mainly applied passive scattering PT, but state-of-the-art PT uses pencil beam scanning (PBS), allowing more conformal dose distributions with lower doses to critical organs [3,13].
PBS-PT is, however, not straightforwardly delivered in the thoracic region due to the inherent sensitivity to uncertainties [14,15]. Much concern has been dedicated to how the breathing motion of the primary tumor can interplay with PBS-PT spot delivery [16,17]. PBS-PT has therefore mainly been offered to patients with limited breathing motion [14]. However, recent studies with PBS-PT confirm that interplay uncertainties are canceled out by fractionation—as for IMRT—and more attention should be focused on changes in breathing patterns and anatomy [15,18,19,20,21].
Various optimization techniques for PBS-PT exist, and it is believed that these are differently influenced by uncertainties [14,16]. With single-field uniform dose (SFUD), each field delivers a uniform dose to the entire target volume, while intensity-modulated proton therapy (IMPT) contains non-uniform dose distributions for the individual fields. 3D and 4D robust optimization can be applied to account for uncertainties due to patient setup and proton range as well as breathing motion, respectively [22,23,24].
In theory, IMPT has the potential to produce the most conformal treatment plans, while SFUD and 4D robust optimization are strategies to increase robustness [17]. There is, however, limited knowledge on how different robustly optimized PBS-PT techniques perform in practice for LA-NSCLC patients, as the planning CT is commonly used for both optimization and evaluation [16]. We see a need to investigate this in order to guide the use of PBS-PT, balancing organ at risk (OAR) sparing and robustness for both targets and OARs. Furthermore, little knowledge exists on the dosimetric advantages of PBS-PT compared to state-of-the-art IMRT at the start of treatment. Comparisons for, e.g., patient selection between protons and photons are also usually carried out on the planning scan, even though it is known that robustness towards changes matters [25]. A few studies have focused on the impact of anatomical changes occurring during the six weeks of treatment that should be handled by means of adaptive radiotherapy (ART) [26,27,28].
The purpose of this prospective simulation study was to compare 3D robust SFUD, 3D robust IMPT and 4D robust IMPT in terms of target coverage and OAR sparing under the influence of setup and range uncertainties, breathing motion and interplay at planning, as well as changes in the breathing motion pattern and anatomy from the planning to the start of treatment. Further, using the clinical IMRT plan as a reference, our objective was to evaluate if the potential for OAR sparing expected at planning was persistent at the start of treatment.

2. Materials and Methods

2.1. Patient Material and Clinical IMRT Planning

Fifteen consecutive patients with stage III NSCLC receiving radiochemotherapy with curative intent at Haukeland University Hospital in Bergen, Norway, in 2019–2020 were prospectively included in an in silico simulation study. All patients gave informed consent, and the study was approved by the regional committee for medical and health research ethics (protocol code 2019/749).
Imaging was performed on a Big Bore CT scanner (Philips Healthcare, Best, The Netherlands), using a Posirest-2 support device (Civco Radiotherapy, Coralville, IA, USA) for fixation in the supine position with arms resting above the head, and the Philips bellows device for registration of the breathing curve.
4DCTs with 10 respiratory phases and deep inspiration breath hold (DIBH) CTs were acquired at planning (Plan) and at the start of treatment (fraction 2 or 3; Start). The average intensity projection (AIP) of the 4DCT was used for delineation and treatment planning. Gross tumor volumes (GTVs) for the primary tumor and lymph nodes were defined on the AIP, based on a diagnostic CT with intravenous contrast, an FDG-PET-CT, and biopsy of mediastinal lymph nodes. To define the internal GTVs (IGTVs), each 4DCT phase was blended with the AIP, and the structure was expanded to include the GTV positions on all phases. Exceptions from this were three patients treated in DIBH due to lung dose exceeding the constraints or large tumor motion blurring the 4DCT. In these cases, IGTV delineation included the GTV on three consecutive planning DIBH scans. A clinical target volume (CTV) was created using a 5 mm margin from the IGTV, without extending into uninvolved organs such as bone, heart, esophagus and major vessels. The GTVs and CTV were deformably mapped to all phases of the 4DCT and later used in 4DIMPT optimization. For clinical IMRT planning, a planning target volume (PTV) with 5 mm margin from the CTV was used. Target delineation was performed by the same oncologist (I.M.S.) on all Plan and Start scans. The lungs, heart, esophagus, spinal canal and, if relevant, the brachial plexus were delineated according to RTOG guidelines [29].
Clinical treatment planning was performed in Eclipse v. 15.6 (Varian Medical Systems, Palo Alto, CA, USA). All patients received IMRT with a prescribed dose of 60 or 66 Gy in 2 Gy fractions, depending on lung function, lung dose and proximity of the brachial plexus to the PTV. The beam configuration was adjusted to fit the anatomy of each patient, mainly using six beams and avoiding entry through the contralateral lung. For the PTV, D98% > 95% of prescribed dose was required, and the maximum dose in the plan should be <107%. Dose constraints for OARs are listed in Table A1. The Acuros External Beam algorithm was used for dose calculation, and the plans were normalized to the median dose in the PTV.
The motion amplitude of the primary tumor at Plan and Start was evaluated in Eclipse, using deformable mapping of the primary tumor GTV from the AIP to each breathing phase of the 4DCT and measuring the motion of the GTV center of mass in all directions.
The AIP of the 4DCT acquired at Start was rigidly matched to the AIP at Plan using six degrees of freedom and a volume of interest covering the PTV, as well as skeletal structures and the body contour in proximity to the PTV.

2.2. Proton Therapy Planning

Proton planning was performed in RayStation v. 8B (RaySearch Laboratories, Stockholm, Sweden). To ensure high plan quality, all plans were made by an experienced planning expert within photon therapy and comparative proton planning (C.G.B.) and reviewed by an experienced medical physicist (M.U.). The 4DCT phases were deformably registered to their respective AIP, and the deformed target volumes (GTVs and CTVs) and OARs were mapped onto each phase. For the AIP scan, a density override representative for tumor tissue (1.06 g/cm3, ~40 HU) was used for all plans for the IGTV. For the 4DCT phases, the original density values were applied (i.e., no density override).
For each patient, three PBS-PT treatment plans were created on the Plan AIP using different optimization techniques: SFUD, 4DIMPT and IMPT. 3D robust optimization according to the minimax approach with setup uncertainty of 5 mm in each direction and 3.5% range uncertainty (21 scenarios) was used for SFUD and IMPT [22]. 4D robust optimization, applying the same settings for setup and range uncertainty on all 4DCT breathing phases (231 scenarios), was used for 4DIMPT [23]. In 3D and 4D robust optimization, the reference plan is evaluated in each uncertainty scenario, and in each iteration, the scenario with the currently worst objective value is improved. 3D and 4D robust optimization were applied for the CTV, and for the spinal canal if close to the CTV. Rescanning methods were not used.
Each plan had two (10 patients) or three (5 patients) coplanar fields with gantry angles carefully selected with regard to the patient anatomy and the distance between the beam entry and the CTV (Figure A1). For each patient, beam angles were individually selected, and the same field setup was used in the three PT plans. Range shifters of 4 cm or 7.5 cm were used for all fields, and the air gaps were 5–12 cm from the body contour depending on beam angles and risk of collision. Sigma of spot sizes in air at isocenter (without range shifter) were 3.7 to 7.2 mm depending on energy.
The same prescription as in the clinical plan was used, applying a relative biological effectiveness of 1.1 for protons. A generic IBA beam model was used for planning, and a Monte Carlo algorithm was used for dose calculation (using 0.5% statistical uncertainty).

2.3. Robustness Evaluation

An overview of the acquired image data and the robustness evaluation is shown in Figure 1. Robustness towards setup and range variations (Plan S/R) was evaluated on the Plan AIP using combined isocenter shifts of 2.9 mm in 3 directions simultaneously (corresponding to 5 mm isotropic shifts) and 3.5% range uncertainty (16 scenarios).
Robustness towards breathing motion at Plan (Plan CT0/50) was evaluated by recalculating all PT plans on the extreme breathing phases of the 4DCT: CT0 (maximum inspiration) and CT50 (maximum expiration). In addition, interplay evaluation was performed at Plan (Plan Interplay) using a script provided by RaySearch. The 10 breathing phases of the 4DCT were in turn used as the starting phase for treatment delivery, and the spots were distributed on the CTs of the different phases based on delivery time and breathing cycle length. Constant breathing periods of five seconds were used for all patients. The dose on each phase was calculated and mapped to the reference image (AIP), where the total dose was calculated. This resulted in 10 different interplay dose distributions depending on which phase delivery started in. For all robustness simulations, reported values represent the worst-case scenario for each parameter.
All PT plans were also recalculated on the AIP (Start), CT0 and CT50 (Start CT0/50) of the Start 4DCT to evaluate robustness towards changes in breathing motion and anatomy that can occur between planning and onset of treatment.

2.4. Dosimetric Evaluation

Dose distributions at Plan were compared using D98% and D2% for the CTV, as well as the homogeneity index HI = (D2% − D98%)/D50%, and the conformity index CI = (TVRI/TV) × (TVRI/VRI), where TV is the target volume, TVRI is the target volume covered by the reference (95%) isodose and VRI is the volume of the reference isodose [30]. For healthy tissue and OARs, the following parameters relevant for toxicity were evaluated: D2cc for the patient body, Dmean, V5Gy and V20Gy for the lungs, Dmean and V30Gy for the heart, Dmean for the esophagus and Dmax for the spinal canal.
For OARs, the planning criteria (Table A1) were also required in robustness evaluation. In addition, the D2cc to the patient body should be <107% of the prescribed dose. CTV D98 > 95% and CTV D2% < 107% were required in setup and range and extreme phase evaluation as well as in the Start recalculations. The interplay effect is expected to cause under- and overdosage in the tumor and OARs that average out during fractionated treatment. Ensuring at least 1.8 Gy per fraction, i.e., CTV D98 > 90%, and CTV D2% and body D2cc < 110% were considered acceptable in interplay evaluation.
A structured overview of the various evaluations and criteria is shown in Table A2. Initially, we present the target coverage and OAR sparing for the various techniques at Plan. Thereafter, we investigate the robustness of the target dose and OAR doses, respectively. For evaluation of the actual clinical potential of proton therapy compared to photon therapy, we lastly compare target coverage and OAR sparing at Start.

2.5. Statistical Analysis

Statistical analyses were performed in IBM SPSS Statistics (IBM Corp., Armonk, NY, USA). Friedman’s test (non-parametric two-way analysis of variance by ranks) was used for comparison of the different techniques. Bonferroni correction was applied to adjust the p-value for multiple testing in post hoc analysis. A significance level of 0.05 was used.

3. Results

3.1. Patient Characteristics and Breathing Motion

The median CTV volume was 137 cc (range 66–435 cc). Nine patients had disease stage IIIA, five IIIB and one IIIC. Primary tumor positions were left upper (3), left lower (4), right upper (4) and right lower (3) lobe. One patient only had mediastinal lymph nodes and one only had a primary tumor; the rest had both primary tumor and lymph nodes included in the target volume. The prescribed dose was 60 Gy for 6 patients and 66 Gy for 9 patients.
The breathing motion of the primary tumor was largest in the cranio–caudal direction, with a median amplitude of 4 mm and a maximum of 15 mm in the planning 4DCTs. Large variability in breathing motion was observed between patients. Six patients had a motion amplitude >5 mm, all in the cranio–caudal direction. Median breathing motion amplitudes were similar at Plan and Start (Table A3), and for most patients, the change in amplitude from Plan to Start was ≤2 mm in all directions. Three patients had a larger change in amplitude in the cranio–caudal direction (−6 mm, +4 mm and −3 mm), and these were also the three patients with the largest breathing motion amplitudes at Plan.

3.2. Target Coverage and OAR Sparing at Plan

All treatment plans achieved the required CTV D98% > 95% and D2% < 107% of the prescribed dose at Plan (Table 1). The median CTV D98% was, however, significantly higher for IMRT than for all PT techniques. The median PTV D98% in the IMRT plans was 95.7% (range 94.6–97.0%). Healthy tissue and OAR doses were lower for all proton techniques than for IMRT (Table 1, Figure 2). The only exception was the D2cc of the body, where IMRT gave the lowest dose. Among the PT techniques, significant differences were found between SFUD and IMPT in Dmean for the lungs and esophagus and V20Gy for the lungs, all in favor of IMPT. The mean rank was the worst with IMRT and the best with IMPT for all of the evaluated OAR parameters.

3.3. Target Dose Robustness at Plan and Start

All IMPT and 4DIMPT plans achieved the criteria for D98% and D2% for the CTV on setup and range evaluation and extreme phase evaluation at Plan (Figure 3). One SFUD plan narrowly failed with a D98% of 94.8% on CT0. In interplay evaluations, all plans fulfilled the goal of D98% > 90%. CTV D2% slightly exceeded 107% in interplay evaluations of three plans, two of which were SFUD and one was 4DIMPT. Thus, all PBS-PT techniques had satisfying target robustness at Plan. The results from the extensive robustness evaluation at Plan and Start for CTV D98% are shown in Figure 3, and a summary of D98% and D2% values for all evaluations of all proton techniques are listed in the Appendix A (Table A4).
For the Start AIP recalculation, the CTV D98% was above 95% of the prescribed dose for 13/15 patients with SFUD and 4DIMPT and 12/15 with IMPT (Figure 3). In general, the CTV D98% decreased in extreme-phase evaluations, but for 11/15 patients, it was still above 95%, independent of the optimization technique. The differences in the median CTV D98% between the PT techniques were small for both AIP and CT0/50 recalculations at Start (Figure 3, Table A4). It was, however, statistically significant between 4DIMPT and IMPT on CT0/50, in disfavor of IMPT. D2% was similar and <107% for all plans on all scans.
One of the patients that stood out with insufficient CTV coverage at Start (AIP) was patient 3, with D98% of 87.2%, 91.7% and 82.7% for SFUD, 4DIMPT and IMPT, respectively. This patient had a change in breathing pattern between Plan and Start, causing the CTV in the mediastinum to expand 15 mm caudally and 3 mm cranially (Figure A2). Similar changes were seen for patients 11 and 15. The CTV coverage at Start for these patients was not sufficient with any PT optimization technique (Figure 3). For one of the patients with insufficient and one of the patients with sufficient CTV coverage at Start, dose distributions for all PT techniques are shown in Figure A3. For patient 8, the low CTV D98% at the Start CT0/50 was likely caused by delineation uncertainty. The IMRT plans for patients 3 and 11 were planned and recalculated on DIBH CTs. The values for IMRT and PT techniques can therefore not be directly compared.

3.4. OAR Dose Robustness at Plan and Start

In setup and range and extreme-phase evaluations at Plan, both IMPT and 4DIMPT achieved the constraints for OARs (Table A1) in all plans. One of the SFUD plans failed in the setup and range evaluation, exceeding the Dmax criterion for the spinal canal with 52.0 Gy in the worst-case scenario. In interplay evaluation, 26 out of 45 plans had a D2cc to the body >107%; however, only three plans exceeded 110% of the prescribed dose. Two of these were IMPT plans, and one was 4DIMPT. The OAR constraints were met for all patients and all techniques in interplay evaluation.
Relevant dose-volume parameters for OARs at Plan and Start for IMRT and all PT techniques are shown in Figure 2. The pattern of OAR sparing with PT compared to IMRT persisted at Start. Median changes in dose-volume parameters from Plan to Start were 6% or lower for all parameters and all techniques (Table A5). Nevertheless, large variations between patients in the relative change of dose-volume parameters (ranging from −58% to 103%) from Plan to Start were seen for individual patients with all techniques. For most of the patients, constraints were still achieved for all OARs. For one patient, the esophagus shifted towards the CTV, causing a ~30% increase in mean dose to above 35 Gy for all techniques. Hotspots (D2cc > 107%) to the healthy tissue occurred at Start with one IMRT plan and two SFUD plans.

3.5. Target Coverage and OAR Sparing at Start

Out of the PBS-PT techniques, IMPT showed the greatest potential for toxicity reduction. A comparison of all 105 OAR dose-volume parameters calculated at Start resulted in the best mean rank for IMPT (1.51), followed by 4DIMPT (2.06), SFUD (2.56) and IMRT (3.87), with all pairwise comparisons being significant. Figure 4 shows the per-patient advantage of IMPT in the sparing of OAR mean doses, as well as the price to pay in target coverage. The latter was, however, only significantly different between IMPT and IMRT, probably influenced by the use of DIBH for IMRT.
Substantial dose reductions were achieved with IMPT compared to state-of-the-art IMRT for the lungs, heart and spinal canal (Figure 4 and Figure 5). For the lungs, the median Dmean was reduced from 13.7 to 9.6 Gy, V5Gy from 55.1 to 28.4% and V20Gy from 23.4 to 18.6% with IMPT compared to IMRT. The median heart Dmean was reduced from 8.2 to 3.0 Gy, with Dmean < 10 Gy for all patients with IMPT and 10/15 with IMRT. The median heart V30Gy was reduced from 8.3 to 3.6%, the median esophagus Dmean was reduced from 20.1 to 18.1 Gy and the median spinal canal Dmax was reduced from 45.5 to 32.7 Gy. All differences were statistically significant, and among the 105 individual parameters compared, 102 were in favor of IMPT, 2 were in favor of IMRT, and 1 was tied.

4. Discussion

This study shows that the potential for OAR sparing with PBS-PT compared to state-of-the-art IMRT was substantial and persistent from the planning to the start of treatment. Among the various optimization techniques, IMPT spared OARs the most. There were surprisingly small differences between the PBS-PT techniques in the response to various uncertainties, but IMPT was slightly less robust towards breathing motion than 4DIMPT. All techniques were acceptable with respect to robustness evaluations at Plan, including interplay, and also at Start for the majority of patients. However, all robust optimization techniques failed to account for changes in breathing motion patterns occurring in three patients, causing unacceptable coverage of the mediastinal lymph nodes.
Given strategies to recognize patients with altered breathing motion and account for the lack of target robustness in these patients, we believe robustly optimized IMPT and 4DIMPT can reduce the risk of both radiation pneumonitis and heart toxicity compared to IMRT. Lung Dmean, V5Gy and V20Gy were all significantly reduced with IMPT and 4DIMPT, and these parameters have previously been correlated to the probability of radiation pneumonitis [31]. Interestingly, this could potentially be a key to better outcome as well, since patients with radiation pneumonitis have been excluded from adjuvant treatment with immune checkpoint inhibitors [32]. Note that SFUD did not reduce lung doses compared to IMRT in our study, and therefore we would not expect any reduced risk of pneumonitis with SFUD. The reductions in heart dose seen with all PBS-PT techniques are also likely clinically relevant. Atkins et al. showed that a heart Dmean > 10 Gy significantly increased the risk of mortality in LA-NSCLC [33]. In our study, a mean dose to the heart below 10 Gy was achieved for 10/15 patients with IMRT and for all patients across all PBS-PT techniques.
Sparing of the spinal canal beyond the max dose constraint is not expected to give a clinical benefit in itself. However, with the large reduction seen with all PBS-PT techniques compared to IMRT, less effort must be spent on this highly prioritized constraint in the optimization, possibly giving room for the considerable dose reduction seen for other OARs.
Mean doses to the esophagus were slightly reduced with IMPT compared to IMRT in our study. It is unknown whether this would lead to a reduction in esophagitis [34], especially since there are additional uncertainties in elevated LET that were not considered in the current study. The esophagus is highly mobile and often located in close proximity to the target volume, and can move into the high-dose region. This was the case for one of the patients in our study. In a recent clinical dose-escalation study (including 47 patients with stage III NSCLC) by Iwata et al., ART was used to monitor the position of the esophagus and adjust treatment accordingly if needed [8]. Dose-escalated PT was well tolerated in this phase II study, with no grade ≥3 radiation pneumonitis and one case of acute grade 3 esophagitis. Additionally, the 5-year overall survival of 59% (probably influenced by combination with immunotherapy) shows promise. This study mainly used passive scattering PT, although some patients with small tumor motion had single-field optimized spot-scanning plans.
Recently, Ribeiro et al. published a comprehensive robustness analysis, including weekly imaging during treatment, for 10 stage III NSCLC patients with small to moderate tumor motion, showing the feasibility of PBS-PT in the majority of patients [27]. Our study strengthens these findings by confirming the results in an independent patient group with larger motion variability. Inoue et al. also investigated the robustness of 3D robustly optimized IMPT in stage III NSCLC [35]. They reported a limited impact of setup and range uncertainties, breathing motion and interplay effects on the dose distribution when using properly selected robust optimization parameters. This is in line with our analysis for the planning scan.
A strength of our study was the prospective study design with repeated imaging at the start of treatment. At this time point, we expected a small probability of anatomical changes in need of ART, based on experience from photon therapy [36]. A CT at fraction 2 or 3 was therefore chosen for robustness evaluation, as it would reveal if any of the optimization techniques were particularly sensitive towards interfractional variations such as changes in breathing pattern or positioning of the patient.
Indeed our results show that none of the optimization techniques for PBS-PT were able to handle substantial changes in the breathing pattern. With current robust optimization methods, it is therefore important to verify dose delivery at the onset of treatment. Adaptive protocols in PT are commonly based on weekly 4DCTs, starting at the end of the first treatment week, but imaging at the onset of treatment could recognize these patients earlier. Importantly, the observed target under-dosage was mainly located in the mediastinal lymph nodes (and not the primary tumor), which are hard to locate on, e.g., CBCT. A possibility is to use the carina as a surrogate structure in addition to the diaphragm, as done by Møller et al. in their ART protocol [36]. The carina position has been shown to correlate better with lung volume than, e.g., diaphragm position [37]. Alternative strategies to avoid dose degradation due to breathing motion changes could be respiratory gating or breath-hold strategies. Although images in DIBH were acquired in the current study and used clinically in IMRT for three patients, analysis of PBS-PT in DIBH was beyond the scope of the current study.
Both Ribeiro et al. and Hoffmann et al. reported that altered shoulder position caused a loss in robustness when evaluating dose during treatment [27,28]. This was not observed in our study, but in principle, this could also occur from the planning to the start of treatment and should be kept in mind when evaluating robustness at the onset of treatment. Our study design was limited to observing changes between the planning and start of treatment, and hence anatomical changes such as atelectasis or pleural effusion were not observed. Such changes can occur during treatment and largely impact the delivered dose, but they are well known and can be corrected for by existing adaptive protocols [28]. The novelty of our study lies in focusing on uncertainties that so far have received less attention. We have shown that these are neither handled by current robust optimization techniques nor adaptive protocols.
Regarding the comparison of different PBS-PT optimization techniques, Ribeiro et al. compared 3D and 4D robustly optimized IMPT plans with layered rescanning in their study [27]. Similar to us, they found only small differences in robustness between the techniques. However, IMPT was (somewhat surprisingly) slightly more robust than 4DIMPT in their study, while we found the opposite. This might be explained by the difference in the use of density override for the target. In the study by Ribeiro et al., density override was only used for the IMPT plans, while we used it on the AIP for both techniques. This is an example of one out of several technical details that might influence robustness; rescanning is another [14]. Indeed, with the use of rescanning, the uncertainties due to the interplay effect could be limited even further than reported here. Liao and colleagues have pointed out the importance of treatment planning experience in PT for NSCLC [11]. In addition to comprehensive treatment planning guidelines, solutions for automated treatment planning could be useful to ensure the high plan quality needed in PBS-PT for LA-NSCLC [38].
The number of treatment fields could also influence the robustness of the PT plans. In this study, two fields were used for ten patients and three fields for five patients. On the one hand, adding a third field could increase the robustness, as the dose contribution is divided between more treatment angles, and changes in anatomy affecting one of the fields have a lesser impact on the dose distribution. However, some issues came with increasing the number of fields. For some patients, finding a third, robust angle could be difficult due to, e.g., arm position or large breasts or fat folds, where it was preferred to avoid beam entry due to positioning uncertainty. In the 4D optimization, splitting the fields with field-specific targets was not possible, so the fields had to be able to contribute to both the primary tumor and the lymph node volumes, giving some limitations for robust angles because of the surrounding anatomy. Hence, the requirement for the field setup to work for all optimization techniques was a limitation in this study.
Another limitation of the current study is the low number of included patients. Despite this, there was a large variation in tumor size and position, and breathing motion ranged from negligible to substantial. These parameters also varied among the patients that failed the robustness criteria at Start. Finally, the 4D optimization and extreme phase and interplay evaluations performed in this study required deformable image registration and mapping of contours to each phase of the 4DCT. As delineation was performed on the AIP as a part of the clinical routine, the contours were mapped from the AIP to each phase. Due to blurring of the edges, the GTV on the AIP may be slightly larger than in reality, and the plans may therefore be slightly more robust than if delineation had been performed on one of the phase images.

5. Conclusions

The potential of IMPT and 4DIMPT for reducing heart and lung toxicity in the treatment of LA-NSCLC was substantial and persistent at Start. SFUD only showed potential for reduced heart toxicity. All proton optimization techniques responded similarly to uncertainties and were sufficiently robust towards setup and range uncertainties as well as interplay at Plan, and for the majority of patients in recalculations at Start. Altered breathing patterns between Plan and Start jeopardized target coverage for all PBS-PT techniques. Adaptive protocols for free-breathing PBS-PT should include imaging at onset of or early in treatment, and possibly a surrogate for visualization of the mediastinal target. Given such strategies to recognize patients with altered breathing patterns, we believe there is great potential for PBS-PT to improve the treatment of LA-NSCLC.

Author Contributions

Conceptualization, C.G.B., I.M.S., M.U. and L.B.H.; methodology, C.G.B., M.U., G.M.E. and L.B.H.; software, C.G.B. and M.U.; validation, K.F., I.M.S. and M.U.; formal analysis, C.G.B., K.F. and L.B.H.; investigation, C.G.B., K.F., I.M.S., M.U., G.M.E. and L.B.H.; resources, L.B.H.; data curation, K.F., I.M.S. and L.B.H.; writing—original draft preparation, C.G.B., K.F. and L.B.H.; writing—review and editing, C.G.B., K.F., I.M.S., M.U., G.M.E. and L.B.H.; visualization, C.G.B., K.F. and L.B.H.; supervision, G.M.E. and L.B.H.; project administration, L.B.H.; funding acquisition, C.G.B., K.F. and L.B.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from the Trond Mohn Foundation (grant number BFS2017TMT07) and Helse Vest RHF (grant number F-12505).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the regional committee for medical and health research ethics in Western Norway (protocol code 2019/749, 21 June 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons as they are part of an ongoing study.

Acknowledgments

The authors are thankful to Sara Pilskog for scientific discussions and clinical personnel at Haukeland University Hospital for data collection.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Planning dose constraints for organs at risk.
Table A1. Planning dose constraints for organs at risk.
OrganDose Constraint
LungsV5Gy < 65%
V20Gy < 35%
Dmean < 20 Gy
EsophagusDmean < 34 Gy
HeartV30Gy < 40%
Spinal canalDmax < 50 Gy
Brachial plexusDmax < 66 Gy
Table A2. Overview of the performed robustness evaluations, and the criteria used in each evaluation. The criteria in parentheses were used for interplay evaluation. * OAR constraints are shown in Table A1.
Table A2. Overview of the performed robustness evaluations, and the criteria used in each evaluation. The criteria in parentheses were used for interplay evaluation. * OAR constraints are shown in Table A1.
Evaluations ofEvaluations at Plan Evaluations at StartEvaluation Criteria
Target coverage and OAR sparing at PlanPlan AIP CTV D98% > 95%
CTV D2% < 107%
Body D2cc < 107%
OARs within constraints *
Target dose robustnessPlan S/R
Plan CT0/50
Plan Interplay
Start AIP
Start CT0/50
CTV D98% > 95% (90%)
CTV D2% < 107% (110%)
OAR dose robustnessPlan S/R
Plan CT0/50
Plan Interplay
Start AIP
Start CT0/50
Body D2cc < 107% (110%)
OARs within constraints *
Target coverage and OAR sparing at Start Start AIPCTV D98% > 95%
CTV D2% < 107%
Body D2cc < 107%
OARs within constraints *
Table A3. Breathing motion amplitudes of the primary tumor in the Plan and Start 4DCTs. Median value and range are given for each direction.
Table A3. Breathing motion amplitudes of the primary tumor in the Plan and Start 4DCTs. Median value and range are given for each direction.
DirectionMotion—PlanMotion—Start
x (left–right)1 mm (0–5)1 mm (0–6)
y (anterior–posterior)2 mm (1–4)2 mm (1–4)
z (cranio–caudal)4 mm (1–15)4 mm (1–13)
Table A4. Robust evaluation of CTV dose, showing median values and range.
Table A4. Robust evaluation of CTV dose, showing median values and range.
ParameterEvaluationSFUD4DIMPTIMPT
CTV D98% (%)Plan CT0/5097.5 (94.8–98.8)98.0 (97.0–99.0)97.8 (96.5–98.5)
Plan setup/range 97.0 (95.6–97.7)97.2 (95.6–97.7)97.3 (95.1–98.1)
Plan interplay 93.6 (90.7–97.5)95.0 (93.5–97.0)94.6 (91.8–97.0)
Start97.6 (85.3–98.8)97.6 (86.2–98.7)97.5 (81.8–98.2)
Start CT0/5097.4 (43.2–98.8)97.4 (61.3–98.9)97.1 (35.6–98.4)
CTV D2% (%) Plan CT0/50102.2 (101.2–103.2)102.5 (102.0–103.5)102.5 (102.0–103.9)
Plan setup/range 103.3 (102.2–104.5)103.2 (102.3–104.2)103.3 (102.5–105.5)
Plan interplay 104.9 (102.7–107.6)105.3 (103.2–107.5)105.5 (104.1–106.7)
Start102.6 (101.4–104.4)102.8 (102.0–103.7)102.6 (101.8–103.7)
Start CT0/50102.3 (101.1–104.9)102.5 (102.0–104.0)102.5 (101.8–103.8)
Table A5. Difference in OAR parameters from Plan to Start, relative to the Plan value. Median and range are shown for each technique. Positive values indicate a higher dose and negative values indicate a lower dose at Start.
Table A5. Difference in OAR parameters from Plan to Start, relative to the Plan value. Median and range are shown for each technique. Positive values indicate a higher dose and negative values indicate a lower dose at Start.
StructureDifference, Plan vs. Start
IMRTSFUD4DIMPTIMPT
Lungs Dmean (Gy)1 (−6/20)0 (−11/20)0 (−11/33)0 (−12/32)
Lungs V5Gy (%)−1 (−5/7)−1 (−11/26)−2 (−11/26)−1 (−12/28)
Lungs V20Gy (%)1 (−7/29)3 (−13/30)2 (−10/30)3 (−11/31)
Heart Dmean (Gy)−4 (−13/68)−6 (−55/70)−4 (−47/78)−6 (−51/78)
Heart V30Gy (%)0 (−29/162)0 (−58/84)0 (−42/97)0 (−50/103)
Esophagus Dmean (Gy)−3 (−17/28)−5 (−27/30)−4 (−31/31)−4 (−31/35)
Spinal canal Dmax (Gy)−1 (−9/6)1 (−14/6)0 (−16/9)1 (−21/8)
Body D2cc (Gy)0 (−2/6)0 (−1/4)0 (0/1)0 (−3/1)
Figure A1. Common beam configurations for different tumor locations: (a) frontal oblique tumor with oblique and frontal fields, (b) nodes in the frontal part of the mediastinum with two slightly frontal oblique fields, (c) dorsal tumor with dorsal oblique and dorsal fields, and (d) large tumor near the brachial plexus, with slightly frontal oblique, dorsal and dorsal oblique fields.
Figure A1. Common beam configurations for different tumor locations: (a) frontal oblique tumor with oblique and frontal fields, (b) nodes in the frontal part of the mediastinum with two slightly frontal oblique fields, (c) dorsal tumor with dorsal oblique and dorsal fields, and (d) large tumor near the brachial plexus, with slightly frontal oblique, dorsal and dorsal oblique fields.
Cancers 14 01365 g0a1
Figure A2. SFUD plan for patient 3. Deeper breathing and increased breathing motion causes the CTV (pink) in the mediastinum to extend 15 mm more caudally and 3 mm more cranially on the Start AIP (right) compared to the Plan AIP (left). The yellow cross shows a reference position in both images. The 95% isodose (green) does not cover the caudal area of the CTV at Start. Notably, the diaphragm position had also changed in this patient, and this change would be detectable on CBCT imaging.
Figure A2. SFUD plan for patient 3. Deeper breathing and increased breathing motion causes the CTV (pink) in the mediastinum to extend 15 mm more caudally and 3 mm more cranially on the Start AIP (right) compared to the Plan AIP (left). The yellow cross shows a reference position in both images. The 95% isodose (green) does not cover the caudal area of the CTV at Start. Notably, the diaphragm position had also changed in this patient, and this change would be detectable on CBCT imaging.
Cancers 14 01365 g0a2
Figure A3. SFUD (left), 4DIMPT (middle), and IMPT (right) plans recalculated on the Start AIP for (a) one of the patients that failed the robustness criteria (patient 11) and (b) one of the patients with satisfactory target coverage (patient 10). The CTV is delineated in pink. Doses higher than 95% of the prescribed dose are shown, according to the color scale in the top right corner.
Figure A3. SFUD (left), 4DIMPT (middle), and IMPT (right) plans recalculated on the Start AIP for (a) one of the patients that failed the robustness criteria (patient 11) and (b) one of the patients with satisfactory target coverage (patient 10). The CTV is delineated in pink. Doses higher than 95% of the prescribed dose are shown, according to the color scale in the top right corner.
Cancers 14 01365 g0a3

References

  1. Bradley, J.D.; Hu, C.; Komaki, R.R.; Masters, G.A.; Blumenschein, G.R.; Schild, S.E.; Bogart, J.A.; Forster, K.M.; Magliocco, A.M.; Kavadi, V.S.; et al. Long-Term Results of NRG Oncology RTOG 0617: Standard-versus High-Dose Chemoradiotherapy with or without Cetuximab for Unresectable Stage III Non–Small-Cell Lung Cancer. J. Clin. Oncol. 2020, 38, 706–714. [Google Scholar] [CrossRef] [PubMed]
  2. Ma, L.; Men, Y.; Feng, L.; Kang, J.; Sun, X.; Yuan, M.; Jiang, W.; Hui, Z. A Current Review of Dose-Escalated Radiotherapy in Locally Advanced Non-Small Cell Lung Cancer. Radiol. Oncol. 2019, 53, 6–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gjyshi, O.; Liao, Z. Proton Therapy for Locally Advanced Non-Small Cell Lung Cancer. Br. J. Radiol. 2020, 93, 20190378. [Google Scholar] [CrossRef] [PubMed]
  4. Ramella, S.; D’Angelillo, R.M. Proton Beam or Photon Beam Radiotherapy in the Treatment of Non-Small-Cell Lung Cancer. Lancet Oncol. 2020, 21, 873–875. [Google Scholar] [CrossRef]
  5. Chang, J.Y.; Zhang, X.; Wang, X.; Kang, Y.; Riley, B.; Bilton, S.; Mohan, R.; Komaki, R.; Cox, J.D. Significant Reduction of Normal Tissue Dose by Proton Radiotherapy Compared with Three-Dimensional Conformal or Intensity-Modulated Radiation Therapy in Stage I or Stage III Non–Small-Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2006, 65, 1087–1096. [Google Scholar] [CrossRef]
  6. Zhang, X.; Li, Y.; Pan, X.; Xiaoqiang, L.; Mohan, R.; Komaki, R.; Cox, J.D.; Chang, J.Y. Intensity-Modulated Proton Therapy Reduces the Dose to Normal Tissue Compared with Intensity-Modulated Radiation Therapy or Passive Scattering Proton Therapy and Enables Individualized Radical Radiotherapy for Extensive Stage IIIB Non-Small-Cell Lung Cancer: A Virtual Clinical Study. Int. J. Radiat. Oncol. Biol. Phys. 2010, 77, 357–366. [Google Scholar] [CrossRef] [Green Version]
  7. Kim, N.; Myoung Noh, J.; Lee, W.; Park, B.; Park, H.; Young Park, J.; Pyo, H. Proton Beam Therapy Reduces the Risk of Severe Radiation-Induced Lymphopenia during Chemoradiotherapy for Locally Advanced Non-Small Cell Lung Cancer: A Comparative Analysis of Proton versus Photon Therapy. Radiother. Oncol. 2021, 156, 166–173. [Google Scholar] [CrossRef] [PubMed]
  8. Iwata, H.; Akita, K.; Yamaba, Y.; Kunii, E.; Takakuwa, O.; Yoshihara, M.; Hattori, Y.; Nakajima, K.; Hayashi, K.; Toshito, T.; et al. Concurrent Chemo-Proton Therapy Using Adaptive Planning for Unresectable Stage 3 Non-Small Cell Lung Cancer: A Phase 2 Study. Int. J. Radiat. Oncol. Biol. Phys. 2021, 109, 1359–1367. [Google Scholar] [CrossRef]
  9. Chang, J.Y.; Verma, V.; Li, M.; Zhang, W.; Komaki, R.; Lu, C.; Allen, P.K.; Liao, Z.; Welsh, J.; Lin, S.H.; et al. Proton Beam Radiotherapy and Concurrent Chemotherapy for Unresectable Stage III Non–Small Cell Lung Cancer: Final Results of a Phase 2 Study. JAMA Oncol. 2017, 3, e172032. [Google Scholar] [CrossRef] [PubMed]
  10. Nguyen, Q.-N.; Ly, N.B.; Komaki, R.; Levy, L.B.; Gomez, D.R.; Chang, J.Y.; Allen, P.K.; Mehran, R.J.; Lu, C.; Gillin, M.; et al. Long-Term Outcomes after Proton Therapy, with Concurrent Chemotherapy, for Stage II–III Inoperable Non-Small Cell Lung Cancer. Radiother. Oncol. 2015, 115, 367–372. [Google Scholar] [CrossRef] [Green Version]
  11. Liao, Z.; Lee, J.J.; Komaki, R.; Gomez, D.R.; O’Reilly, M.S.; Fossella, F.V.; Blumenschein, G.R.; Heymach, J.V.; Vaporciyan, A.A.; Swisher, S.G.; et al. Bayesian Adaptive Randomization Trial of Passive Scattering Proton Therapy and Intensity-Modulated Photon Radiotherapy for Locally Advanced Non–Small-Cell Lung Cancer. J. Clin. Oncol. 2018, 36, 1813–1822. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, P.; Xu, T.; Gomez, D.R.; Deng, W.; Wei, X.; Elhalawani, H.; Jin, H.; Guan, F.; Mirkovic, D.; Xu, Y.; et al. Patterns of Local-Regional Failure after Intensity Modulated Radiation Therapy or Passive Scattering Proton Therapy with Concurrent Chemotherapy for Non-Small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2019, 103, 123–131. [Google Scholar] [CrossRef]
  13. Gjyshi, O.; Xu, T.; Elhammali, A.; Boyce-Fappiano, D.; Chun, S.G.; Gandhi, S.; Lee, P.; Chen, A.B.; Lin, S.H.; Chang, J.Y.; et al. Toxicity and Survival after Intensity-Modulated Proton Therapy versus Passive Scattering Proton Therapy for NSCLC. J. Thorac. Oncol. 2021, 16, 269–277. [Google Scholar] [CrossRef] [PubMed]
  14. Chang, J.Y.; Zhang, X.; Knopf, A.; Li, H.; Mori, S.; Dong, L.; Lu, H.-M.; Liu, W.; Badiyan, S.N.; Both, S.; et al. Consensus Guidelines for Implementing Pencil-Beam Scanning Proton Therapy for Thoracic Malignancies on Behalf of the PTCOG Thoracic and Lymphoma Subcommittee. Int. J. Radiat. Oncol. Biol. Phys. 2017, 99, 41–50. [Google Scholar] [CrossRef] [PubMed]
  15. Meijers, A.; Knopf, A.-C.; Crijns, A.P.G.; Ubbels, J.F.; Niezink, A.G.H.; Langendijk, J.A.; Wijsman, R.; Both, S. Evaluation of Interplay and Organ Motion Effects by Means of 4D Dose Reconstruction and Accumulation. Radiother. Oncol. 2020, 150, 268–274. [Google Scholar] [CrossRef] [PubMed]
  16. Jie, A.W.; Marignol, L. Pro-Con of Proton: Dosimetric Advantages of Intensity-Modulation over Passive Scatter for Thoracic Malignancies. Tech. Innov. Patient Support Radiat. Oncol. 2020, 15, 37–46. [Google Scholar] [CrossRef] [PubMed]
  17. Paganetti, H.; Grassberger, C.; Sharp, G.C. Physics of Particle Beam and Hypofractionated Beam Delivery in NSCLC. Semin. Radiat. Oncol. 2021, 31, 162–169. [Google Scholar] [CrossRef]
  18. Selvaraj, J.; Uzan, J.; Baker, C.; Nahum, A. 4D Radiobiological Modelling of the Interplay Effect in Conventionally and Hypofractionated Lung Tumour IMRT. Br. J. Radiol. 2015, 88, 20140372. [Google Scholar] [CrossRef] [Green Version]
  19. Den Otter, L.A.; Anakotta, R.M.; Weessies, M.; Roos, C.T.G.; Sijtsema, N.M.; Muijs, C.T.; Dieters, M.; Wijsman, R.; Troost, E.G.C.; Richter, C.; et al. Investigation of Inter-fraction Target Motion Variations in the Context of Pencil Beam Scanned Proton Therapy in Non-small Cell Lung Cancer Patients. Med. Phys. 2020, 47, 3835–3844. [Google Scholar] [CrossRef] [PubMed]
  20. Møller, D.S.; Khalil, A.A.; Knap, M.M.; Hoffmann, L. Adaptive Radiotherapy of Lung Cancer Patients with Pleural Effusion or Atelectasis. Radiother. Oncol. 2014, 110, 517–522. [Google Scholar] [CrossRef]
  21. Kwint, M.; Conijn, S.; Schaake, E.; Knegjens, J.; Rossi, M.; Remeijer, P.; Sonke, J.-J.; Belderbos, J. Intra Thoracic Anatomical Changes in Lung Cancer Patients during the Course of Radiotherapy. Radiother. Oncol. 2014, 113, 392–397. [Google Scholar] [CrossRef] [PubMed]
  22. Fredriksson, A.; Forsgren, A.; Hårdemark, B. Minimax Optimization for Handling Range and Setup Uncertainties in Proton Therapy: Minimax Optimization for Handling Uncertainties in Proton Therapy. Med. Phys. 2011, 38, 1672–1684. [Google Scholar] [CrossRef] [PubMed]
  23. Engwall, E.; Fredriksson, A.; Glimelius, L. 4D Robust Optimization Including Uncertainties in Time Structures Can Reduce the Interplay Effect in Proton Pencil Beam Scanning Radiation Therapy. Med. Phys. 2018, 45, 4020–4029. [Google Scholar] [CrossRef] [Green Version]
  24. Ge, S.; Wang, X.; Liao, Z.; Zhang, L.; Sahoo, N.; Yang, J.; Guan, F.; Mohan, R. Potential for Improvements in Robustness and Optimality of Intensity-Modulated Proton Therapy for Lung Cancer with 4-Dimensional Robust Optimization. Cancers 2019, 11, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Teoh, S.; Fiorini, F.; George, B.; Vallis, K.A.; Van den Heuvel, F. Proton vs Photon: A Model-Based Approach to Patient Selection for Reduction of Cardiac Toxicity in Locally Advanced Lung Cancer. Radiother. Oncol. 2020, 152, 151–162. [Google Scholar] [CrossRef] [Green Version]
  26. Chen, M.; Yang, J.; Liao, Z.; Chen, J.; Xu, C.; He, X.; Zhang, X.; Zhu, R.X.; Li, H. Anatomic Change over the Course of Treatment for Non–Small Cell Lung Cancer Patients and Its Impact on Intensity-Modulated Radiation Therapy and Passive-Scattering Proton Therapy Deliveries. Radiat. Oncol 2020, 15, 55. [Google Scholar] [CrossRef]
  27. Ribeiro, C.O.; Visser, S.; Korevaar, E.W.; Sijtsema, N.M.; Anakotta, R.M.; Dieters, M.; Both, S.; Langendijk, J.A.; Wijsman, R.; Muijs, C.T.; et al. Towards the Clinical Implementation of Intensity-Modulated Proton Therapy for Thoracic Indications with Moderate Motion: Robust Optimised Plan Evaluation by Means of Patient and Machine Specific Information. Radiother. Oncol. 2021, 157, 210–218. [Google Scholar] [CrossRef] [PubMed]
  28. Hoffmann, L.; Alber, M.; Jensen, M.F.; Holt, M.I.; Møller, D.S. Adaptation Is Mandatory for Intensity Modulated Proton Therapy of Advanced Lung Cancer to Ensure Target Coverage. Radiother. Oncol. 2017, 122, 400–405. [Google Scholar] [CrossRef]
  29. Kong, F.-M.; Ritter, T.; Quint, D.J.; Senan, S.; Gaspar, L.E.; Komaki, R.U.; Hurkmans, C.W.; Timmerman, R.; Bezjak, A.; Bradley, J.D.; et al. Consideration of Dose Limits for Organs at Risk of Thoracic Radiotherapy: Atlas for Lung, Proximal Bronchial Tree, Esophagus, Spinal Cord, Ribs, and Brachial Plexus. Int. J. Radiat. Oncol. Biol. Phys. 2011, 81, 1442–1457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Feuvret, L.; Noël, G.; Mazeron, J.-J.; Bey, P. Conformity Index: A Review. Int. J. Radiat. Oncol. Biol. Phys. 2006, 64, 333–342. [Google Scholar] [CrossRef] [PubMed]
  31. Marks, L.B.; Bentzen, S.M.; Deasy, J.O.; Kong, F.-M.; Bradley, J.D.; Vogelius, I.S.; El Naqa, I.; Hubbs, J.L.; Lebesque, J.V.; Timmerman, R.D.; et al. Radiation Dose–Volume Effects in the Lung. Int. J. Radiat. Oncol. Biol. Phys. 2010, 76, S70–S76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Antonia, S.J.; Villegas, A.; Daniel, D.; Vicente, D.; Murakami, S.; Hui, R.; Kurata, T.; Chiappori, A.; Lee, K.H.; de Wit, M.; et al. Overall Survival with Durvalumab after Chemoradiotherapy in Stage III NSCLC. N. Engl. J. Med. 2018, 379, 2342–2350. [Google Scholar] [CrossRef]
  33. Atkins, K.M.; Rawal, B.; Chaunzwa, T.L.; Lamba, N.; Bitterman, D.S.; Williams, C.L.; Kozono, D.E.; Baldini, E.H.; Chen, A.B.; Nguyen, P.L.; et al. Cardiac Radiation Dose, Cardiac Disease, and Mortality in Patients with Lung Cancer. J. Am. Coll. Cardiol. 2019, 73, 2976–2987. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, Z.; Chen, M.; Sun, J.; Jiang, S.; Wang, L.; Wang, X.; Sahoo, N.; Gunn, G.B.; Frank, S.J.; Nguyen, Q.-N.; et al. Lyman–Kutcher–Burman Normal Tissue Complication Probability Modeling for Radiation-Induced Esophagitis in Non-Small Cell Lung Cancer Patients Receiving Proton Radiotherapy. Radiother. Oncol. 2020, 146, 200–204. [Google Scholar] [CrossRef] [PubMed]
  35. Inoue, T.; Widder, J.; van Dijk, L.V.; Takegawa, H.; Koizumi, M.; Takashina, M.; Usui, K.; Kurokawa, C.; Sugimoto, S.; Saito, A.I.; et al. Limited Impact of Setup and Range Uncertainties, Breathing Motion, and Interplay Effects in Robustly Optimized Intensity Modulated Proton Therapy for Stage III Non-Small Cell Lung Cancer. Int. J. Radiat. Oncol. Biol. Phys. 2016, 96, 661–669. [Google Scholar] [CrossRef] [PubMed]
  36. Møller, D.S.; Holt, M.I.; Alber, M.; Tvilum, M.; Khalil, A.A.; Knap, M.M.; Hoffmann, L. Adaptive Radiotherapy for Advanced Lung Cancer Ensures Target Coverage and Decreases Lung Dose. Radiother. Oncol. 2016, 121, 32–38. [Google Scholar] [CrossRef] [PubMed]
  37. Van der Weide, L.; van Sörnsen de Koste, J.R.; Lagerwaard, F.J.; Vincent, A.; van Triest, B.; Slotman, B.J.; Senan, S. Analysis of Carina Position as Surrogate Marker for Delivering Phase-Gated Radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 2008, 71, 1111–1117. [Google Scholar] [CrossRef]
  38. Kouwenberg, J.; Penninkhof, J.; Habraken, S.; Zindler, J.; Hoogeman, M.; Heijmen, B. Model Based Patient Pre-Selection for Intensity-Modulated Proton Therapy (IMPT) Using Automated Treatment Planning and Machine Learning. Radiother. Oncol. 2021, 158, 224–229. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overview of the acquired image data and the robustness evaluation.
Figure 1. Overview of the acquired image data and the robustness evaluation.
Cancers 14 01365 g001
Figure 2. OAR and normal tissue doses at Plan (blue) and Start (red) for the different optimization techniques; panel (ah) shows each of the evaluated parameters. Purple horizontal lines indicate planning constraints in cases where the value is included in the plot. Boxplots show the median (line), mean (cross) and spread, with outliers as dots outside the box. IMRT = intensity-modulated radiotherapy, SFUD = single-field uniform dose, 4DIMPT = 4D robustly optimized intensity-modulated proton therapy, IMPT = 3D robustly optimized intensity-modulated proton therapy.
Figure 2. OAR and normal tissue doses at Plan (blue) and Start (red) for the different optimization techniques; panel (ah) shows each of the evaluated parameters. Purple horizontal lines indicate planning constraints in cases where the value is included in the plot. Boxplots show the median (line), mean (cross) and spread, with outliers as dots outside the box. IMRT = intensity-modulated radiotherapy, SFUD = single-field uniform dose, 4DIMPT = 4D robustly optimized intensity-modulated proton therapy, IMPT = 3D robustly optimized intensity-modulated proton therapy.
Cancers 14 01365 g002
Figure 3. Robust evaluation on the Plan (blue) and Start (red) CTs for each optimization technique. The acceptance criteria were D98% > 90% for interplay evaluation and >95% for other robust evaluations. In cases where the criteria were not met, the patient number is given next to the observation in the figure. DIBH CT was used in IMRT for patients 3 and 11, hence planning and recalculation were not performed on the same scans as for the PT plans. CT0/50 includes two observations (both extreme phases) per patient. Boxplots show the median (line), mean (cross) and spread, with outliers as dots outside the box. Plan = planning CT, Start = start of treatment CT, CT0/50 = extreme phase evaluation, S/R = setup and range evaluation.
Figure 3. Robust evaluation on the Plan (blue) and Start (red) CTs for each optimization technique. The acceptance criteria were D98% > 90% for interplay evaluation and >95% for other robust evaluations. In cases where the criteria were not met, the patient number is given next to the observation in the figure. DIBH CT was used in IMRT for patients 3 and 11, hence planning and recalculation were not performed on the same scans as for the PT plans. CT0/50 includes two observations (both extreme phases) per patient. Boxplots show the median (line), mean (cross) and spread, with outliers as dots outside the box. Plan = planning CT, Start = start of treatment CT, CT0/50 = extreme phase evaluation, S/R = setup and range evaluation.
Cancers 14 01365 g003
Figure 4. Absolute difference in mean OAR dose and CTV D98% per patient between IMPT and SFUD (a), IMPT and 4DIMPT (b) and IMPT and IMRT (c) as calculated on the AIP at Start. Negative values are always in favor of IMPT (for the OARs, the other technique is subtracted from IMPT, while for the CTV, IMPT is subtracted from the other technique). The patients are sorted according to increasing breathing motion. Patient 1 had no primary tumor, patients 2–7 had tumor motion amplitude >0.5 cm and patients 8–15 had tumor motion amplitude <0.5 cm. Notably, the patients with a large advantage in CTV D98% for SFUD and 4DIMPT (patients 3, 11 and 15) are the same patients where changes in breathing pattern and anatomy deteriorated CTV coverage for all proton techniques. Note also that three of the patients (2, 3 and 11) received IMRT in DIBH and cannot be directly compared in (c).
Figure 4. Absolute difference in mean OAR dose and CTV D98% per patient between IMPT and SFUD (a), IMPT and 4DIMPT (b) and IMPT and IMRT (c) as calculated on the AIP at Start. Negative values are always in favor of IMPT (for the OARs, the other technique is subtracted from IMPT, while for the CTV, IMPT is subtracted from the other technique). The patients are sorted according to increasing breathing motion. Patient 1 had no primary tumor, patients 2–7 had tumor motion amplitude >0.5 cm and patients 8–15 had tumor motion amplitude <0.5 cm. Notably, the patients with a large advantage in CTV D98% for SFUD and 4DIMPT (patients 3, 11 and 15) are the same patients where changes in breathing pattern and anatomy deteriorated CTV coverage for all proton techniques. Note also that three of the patients (2, 3 and 11) received IMRT in DIBH and cannot be directly compared in (c).
Cancers 14 01365 g004
Figure 5. (a) Population average DVH for IMRT and IMPT plans recalculated at Start. Six patients had CTV_60 and nine had CTV_66. (b) Dose distribution for patient 7 with IMRT (left) and IMPT (right), recalculated at Start. The spread of the low and medium doses results in higher mean doses to the lungs (blue), esophagus (green) and heart (yellow) with IMRT. The CTV is delineated in pink, and the spinal canal in cyan.
Figure 5. (a) Population average DVH for IMRT and IMPT plans recalculated at Start. Six patients had CTV_60 and nine had CTV_66. (b) Dose distribution for patient 7 with IMRT (left) and IMPT (right), recalculated at Start. The spread of the low and medium doses results in higher mean doses to the lungs (blue), esophagus (green) and heart (yellow) with IMRT. The CTV is delineated in pink, and the spinal canal in cyan.
Cancers 14 01365 g005
Table 1. Dose-volume parameters for the target, patient body, and OARs for the different optimization techniques, evaluated on the planning scan. The mean rank (obtained by Friedman’s test) for each technique regarding each evaluated parameter is also shown. A mean rank of 1 would mean that this was the best plan for all patients, while a mean rank of 4 means that this was the worst plan for all patients. Statistically significant differences (p < 0.05) compared to the IMRT plan are shown in bold.
Table 1. Dose-volume parameters for the target, patient body, and OARs for the different optimization techniques, evaluated on the planning scan. The mean rank (obtained by Friedman’s test) for each technique regarding each evaluated parameter is also shown. A mean rank of 1 would mean that this was the best plan for all patients, while a mean rank of 4 means that this was the worst plan for all patients. Statistically significant differences (p < 0.05) compared to the IMRT plan are shown in bold.
IMRTSFUD4DIMPTIMPT
Median (Range)Mean RankMedian (Range)Mean RankMedian (Range)Mean RankMedian (Range)Mean Rank
CTV D98% (%)98.6 (98.2–98.8)1.2098 (97.4–98.8)2.9398 (97.1–98.9)2.9398 (97.1–98.7)2.93
CTV D2% (%)102.2 (101.7–102.7)1.93102.5 (101.5–103.0)2.13102.6 (102–103.4)2.93102.4 (102.2–103.9)3.00
CTV CI0.41 (0.33–0.62)1.700.34 (0.26–0.53)3.330.34 (0.25–0.64)2.500.34 (0.26–0.61)2.47
CTV HI0.036 (0.030–0.044)1.270.046 (0.027–0.051)2.530.047 (0.038–0.059)3.130.044 (0.036–0.068)3.07
Body D2cc (Gy)67.5 (61.3–68.4)1.2068.3 (61.1–69.0)2.4768.1 (62.1–69.7)3.0768.5 (61.9–69.4)3.27
Lungs Dmean (Gy)13.6 (6.6–16.8)3.9310.2 (4.5–14.5) *2.879.6 (4.1–13.2)1.879.2 (4.5–13.1)1.33
Lungs V5Gy (%)54.9 (31.8–63.2)4.0028.3 (14.0–41.9)2.4727.5 (13.1–41.6)2.0027.6 (13.7–40.9)1.53
Lungs V20Gy (%)22.6 (9.6–30.6)3.8019.0 (9.0–29.0) *2.9318.0 (8.2–25.1)1.9317.6 (8.8–25.5)1.33
Heart Dmean (Gy)8.1 (0.9- 20.7)4.002.8 (0.5–10.1)2.532.8 (0.5–10.0)1.932.8 (0.5–10.0)1.53
Heart V30Gy (%)8.6 (0.0–28.2)3.903.2 (0.0–13.7)2.303.3 (0.9–13.9)2.373.3 (0.0–13.7)1.43
Esophagus Dmean (Gy)20.4 (14.1–32.2)3.3320.0 (10.0–30.8) *3.0019.5 (8.9–31.8)2.2720.0 (9.0–29.9)1.40
Spinal Canal Dmax (Gy)46.0 (31.8–53.9)4.0032.3 (18.8–41.1)2.4031.4 (20.6–40.3)2.0032.5 (12.8–39.6)1.60
* Statistically significant difference between SFUD and IMPT.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boer, C.G.; Fjellanger, K.; Sandvik, I.M.; Ugland, M.; Engeseth, G.M.; Hysing, L.B. Substantial Sparing of Organs at Risk with Modern Proton Therapy in Lung Cancer, but Altered Breathing Patterns Can Jeopardize Target Coverage. Cancers 2022, 14, 1365. https://doi.org/10.3390/cancers14061365

AMA Style

Boer CG, Fjellanger K, Sandvik IM, Ugland M, Engeseth GM, Hysing LB. Substantial Sparing of Organs at Risk with Modern Proton Therapy in Lung Cancer, but Altered Breathing Patterns Can Jeopardize Target Coverage. Cancers. 2022; 14(6):1365. https://doi.org/10.3390/cancers14061365

Chicago/Turabian Style

Boer, Camilla Grindeland, Kristine Fjellanger, Inger Marie Sandvik, Maren Ugland, Grete May Engeseth, and Liv Bolstad Hysing. 2022. "Substantial Sparing of Organs at Risk with Modern Proton Therapy in Lung Cancer, but Altered Breathing Patterns Can Jeopardize Target Coverage" Cancers 14, no. 6: 1365. https://doi.org/10.3390/cancers14061365

APA Style

Boer, C. G., Fjellanger, K., Sandvik, I. M., Ugland, M., Engeseth, G. M., & Hysing, L. B. (2022). Substantial Sparing of Organs at Risk with Modern Proton Therapy in Lung Cancer, but Altered Breathing Patterns Can Jeopardize Target Coverage. Cancers, 14(6), 1365. https://doi.org/10.3390/cancers14061365

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop