Evolution of motion uncertainty in rectal cancer: implications for adaptive radiotherapy

Evaluating primary tumor position variation for rectal cancer patients treated with long course radiotherapy

Objective.To quantify interfraction shape and positional variations of primary tumor volumes for rectal cancer patients receiving long course radiotherapy by comparing two quantification strategies: a center-of-mass (COM) method and a surface-based metric that captures local deformations.Approach.This study utilized repeat MRI scans before and during radiotherapy (RT) for rectal cancer to investigate the positional variation of the primary gross tumor volume (GTVp). Sixteen patients underwent six MRI exams, with the initial three before the RT course and the subsequent three at one, two, and four weeks into the RT course. GTVp's were delineated on 3D T2-weighted MRIs, and positional variation analyzed using both COM and point-based surface displacements against the initial scan. Surface displacements were quantified using a bidirectional local distance measure, analyzing 3D displacement vectors. Additionally, the study examined local right-left (RL) and anterior-posterior (AP) surface variations relative to tumor height in the rectum by mapping baseline GTVp volumes onto a reference rectum structure.Main results.Systematic error for COM measurements were 1.7, 1.3 and 2.0 mm for AP, RL, and cranial-caudal (CC) direction, respectively. Random errors were 2.1, 1.2 and 2.2 mm, while the GM errors were -0.3, 0.5 and -0.3 mm for AP, RL, and CC directions, respectively. An increase in systematic and random errors were observed when comparing 95th percentile surface displacements to the COM measurements, indicating local displacements which the COM did not detect. Additionally, a general tendency for higher-located tumors to experience larger left-right and AP surface variations were seen when evaluating the 95th percentile.Significance.COM-based analysis might underestimate local deformations. Consequently, surface-based methods might provide more robust estimations of systematic, random and group mean errors for planning target volume-margin calculation. The surface variations tend to increase for tumors located in the upper part of the rectum.

Investigation of changes in planning target volume and regression probability of rectal boost using in-silico cone-beam computed tomography-guided online-adaptive radiotherapy

Background and purpose

Radiotherapy boost to the primary tumour may enable organ preservation in locally advanced rectal cancer (LARC). This study evaluated cone-beam computed tomography (CBCT)-guided online-adaptive radiotherapy (ART) to reduce rectal boost planning target volume (PTVBoost) margins and allow dose escalation.

Materials and methods

Eleven LARC patients were included in this in silico study. Population-based PTVBoost margins were computed for non-adaptive and online-ART using van Herk's formalism. Dose/volume results were compared between: non-adaptive RT with a 25 x 2.16 Gy boost (Non-ART54Gy), ART with a 25 x 2.16 Gy boost (ART54Gy), and ART with an escalated boost of 25 x 2.4 Gy (ART60Gy). Tumour regression probability was compared between each plan using a dose-response model.

Results

PTVBoost margins for non-adaptive vs. online-ART were 14.2 vs. 3.3 mm in the antero-posterior, 5.0 vs. 3.2 mm in the left-right, and 12.3 vs. 8.7 mm in the supero-inferior axes. PTVBoost and pelvic lymph node PTV coverage (V95%) were significantly improved with ART54Gy and ART60Gy compared to Non-ART54Gy (p < 0.001). High-priority organ-at-risk constraints (priority 1&2) were violated in 26.8 % of cases for Non-ART54Gy, 21.2 % of cases for ART54Gy, and 20.8 % of cases for ART60Gy. Tumour regression probability was superior for ART60Gy (20.8 %) compared to ART54Gy (17.0 %, p < 0.001) and Non-ART54Gy (16.9 %, p < 0.001).

Conclusions

Online-ART significantly reduce rectal boost PTV margin. It allows better target volume coverage with a similar risk of radiation-induced toxicities, even when escalating the dose. Therefore, online-ART should be considered to perform dose-escalation in LARC patients with the objective of organ preservation.

Incorporating patient-specific prior clinical knowledge to improve clinical target volume auto-segmentation generalisability for online adaptive radiotherapy of rectal cancer: A multicenter validation

BACKGROUND & PURPOSE

Deep learning (DL) based auto-segmentation has shown to be beneficial for online adaptive radiotherapy (OART). However, auto-segmentation of clinical target volumes (CTV) is complex, as clinical interpretations are crucial in their definition. The resulting variation between clinicians and institutes hampers the generalizability of DL networks. In OART the CTV is delineated during treatment preparation which makes the clinician intent explicitly available during treatment. We studied whether multicenter generalisability improves when using this prior clinical knowledge, the pre-treatment delineation, as a patient-specific prior for DL models for online auto-segmentation of the mesorectal CTV.

MATERIAL & METHODS

We included intermediate risk or locally advanced rectal cancer patients from three centers. Patient-specific weight maps were created by combining the patient-specific CTV delineation on the pre-treatment scan with population-based variation of likely inter-fraction mesorectal CTV deformations. We trained two models to auto-segment the mesorectal CTV on in-house data, one with (MRI + prior) and one without (MRI-only) priors. Both models were applied to two external datasets. An external baseline model was trained without priors from scratch for one external center. Performance was evaluated on the DSC, surface Dice, 95HD and MSD.

RESULTS

For both external centers, the MRI + prior model outperformed the MRI-only model significantly on the segmentation metrics (p-values < 0.01). There was no significant difference between the external baseline model and the MRI + prior model.

CONCLUSION

Adding patient-specific weight maps makes the CTV segmentation model more robust to institutional preferences. Performance was comparable to a model trained locally from scratch. This makes this approach suitable for generalization to multiple centers.

Estimation of adaptive radiation therapy requirements for rectal cancer: a two-center study

BACKGROUND

Rectal cancer patients are potential beneficiaries of adaptive radiotherapy (ART) which demands considerable resources. Currently, there is no definite guidance on what kind of patients and when will benefit from ART. This study aimed to develop and validate a methodology for estimating ART requirements in rectal cancer before treatment course.

METHODS AND MATERIALS

This study involved 66 rectal cancer patients from center 1 and 27 patients from center 2. The ART requirements were evaluated by comparing 8 dose volume histogram (DVH) metrics of targets and organs at risk (OARs) between planning and treatment fractions. Tolerance ranges of deviation of DVH metrics were derived from 10 patients and applied to assess fractional variability. Eighteen features, encompassing diagnostic, dosimetric, and time-related information, were utilized to formulate a stepwise logistic regression model for fraction-level ART requirement estimation. The super parameters were determined through 5-fold cross-validation with 250 training fractions and the methodology was validated with 109 internal testing fractions and 134 external testing fractions.

RESULTS

The area under the curve (AUC) of training dataset was 0.74 (95% CI: 0.61 to 0.85), while in the internal and external testing, the AUC achieved 0.76 (95% CI: 0.60-0.90) and 0.68 (95% CI: 0.56-0.81). Using a best (or clinical applicable) cut-off value of 33.4% (11%), the predictive model achieved a sensitivity of 46.2% (69.2%) and specificity of 97.9% (68.7%). During the modeling, 5 features were retained: Homogeneity index (OR = 6.06, 95% CI: 2.93-14.8), planning target volume (OR = 1.77, 95% CI: 1.17-2.69), fraction dose (OR = 45.37, 95% CI: 5.74-469), accumulated dose (OR = 2.29, 95% CI: 1.35-4.14), and whether neoadjuvant chemoradiotherapy (OR > 1000).

CONCLUSION

ART requirements are associated with target volume, target dose homogeneity, fraction dose, dose accumulation and whether neoadjuvant radiotherapy. The predictive model exhibited the capability to predict fraction-level ART requirements.

Dealing with rectum motion during radiotherapy: How can we anticipate it?

Introduction

Intra- and inter-fraction rectum motion is important for pelvic radiotherapy (RT). This study assesses how RT session duration, the presence or the absence of an intra-rectal tumour, and the distance from the anorectal junction (ARJd) impact rectal motion.

Materials and methods

Analyses used cone-beam computed tomographies (CBCTs) from RT patients treated for rectal and prostate cancer. Three structures were evaluated: (1) the entire rectum in patients without a rectal tumour (RectumProstate); (2) the non-invaded portion (RectumRectum) and (3) the tumour-invaded portion (RectumTumour) in rectal cancer patients.Intrafraction motion was assessed using the Hausdorff distance 95% and the Mean distance-to-agreement between structures delineated on the first CBCT and the 2 subsequent CBCTs within a same RT session. Interfraction motion was quantified by comparing structures delineated on the planning-CT and the first CBCT of each session.Linear mixed model evaluated rectum motion in relation to time, tumour presence, and ARJd, respectively.

Results

We included 10 patients with and 10 without rectal cancer, collecting 385 CBCTs. A significant correlation (p < 0.05) between rectum motion and RT session duration was found. Intrafraction motion was significantly higher in prostate cancer patients (RectumProstate motion > RectumRectum and RectumTumour, p < 0.01). For interfraction motion, only the mean distance to agreement was significantly higher for RectumProstate (p < 0.05). Motion increased significantly with ARJd for all three structures (p < 0.001).

Conclusions

Session duration, absence of a tumour, and ARJd are associated with larger intra- and interfraction rectal motion. This highlights the need for tailored RT treatment, including online-adaptive RT, to manage intra- and interfraction variations. Rectal motion should be handled differently for patients with prostate cancer and those with rectal cancer.

On the trail of CBCT-guided adaptive rectal boost radiotherapy, does daily delineation require a radiation oncologist?

Introduction

Dose-escalation radiotherapy for rectal tumours is increasingly considered as a non-operative approach, with online-adaptive radiotherapy (oART) supporting this approach by correcting inter-fraction tumour position errors. However, using cone-beam computed tomography (CBCT)-guided oART requires daily target volume delineation by different operators, leading to inter-operator delineation variability and potential dosimetric issues. This study aims to compare and quantify the inter-operator and inter-professional delineation variability of the rectal boost volume on CBCT, including volumes by an automatically delineated oART treatment planning system.

Materials and methods

A rectal boost volume, defined as the primary tumour extended to the entire adjacent rectal wall, was delineated on 10 CBCTs from 5 patients by 15 operators: 4 expert radiation oncologists (ROs), 4 radiation therapists (RTTs) and 7 non-expert ROs. These contours were compared between the different professional groups. A comparison to the average volume of the group (ROs, RTTs, or non-expert ROs) with the lowest delineation variability was also performed for each individual volume including the volume automatically generated by an oART treatment planning system.

Results

Delineation variability was the highest in the superior (range: 2.3-6.0 mm), and inferior (2.3-12.4 mm) directions, compared to the left (0.2-4.4 mm), right (0.3-2.0 mm), anterior (0.1-2.9 mm), and posterior (0.5-4.0 mm) directions. Non-expert ROs, RTTs, and automatic oART volume showed similar ranges of delineation errors when compared to the expert ROs' volume, which was chosen as reference volume since this professional group showed the lowest variability.

Discussion

Expert ROs showed consistent results. Other professional groups exhibit similar variability, comparable to the automatic oART volume. Therefore, RTTs could safely perform the rectal boost delineation without non-expert ROs supervision in the absence of expert ROs during CBCT-based oART. Moreover, these findings provide quantitative data to compute accurate margins for the rectal boost planning target volume in a CBCT-guided oART workflow.

Individual lymph node position variation for rectal cancer patients treated with long course chemoradiotherapy

Background and purpose

Delivery of high precision radiotherapy lymph node boosts requires detailed information on the interfraction positional variation of individual lymph nodes. In this study we characterized interfraction positional shifts of suspected malignant lymph nodes for rectal cancer patients receiving long course radiotherapy. Furthermore, we investigated parameters which could affect the magnitude of the position variation.

Materials and Methods

Fourteen patients from a prospective clinical imaging study with a total of 61 suspected malignant lymph nodes in the mesorectum, presacral, and lateral regions, were included. The primary gross tumor volume (GTVp) and all suspected malignant lymph nodes were delineated on six magnetic resonance imaging scans per patient. Positional variation was calculated as systematic and random errors, based on shifts of center-of-mass, and estimated relative to either bony structures or the GTVp using a hierarchical linear mixed model.

Results

Depending on location and direction, systematic and random variations (relative to bony structures) were within 0.6-2.8 mm and 0.6-2.9 mm, respectively. Systematic and random variations increased when evaluating position relative to GTVp (median increase of 0.6 mm and 0.5 mm, respectively). Correlations with scan time-point and relative bladder volume were found in some directions.

Conclusions

Using linear mixed modeling, we estimated systematic and random positional variation for suspected malignant lymph nodes in rectal cancer patients treated with long course radiotherapy. Statistically significant correlations of the magnitude of the lymph node shifts were found related to scan time-point and relative bladder volume.

Evaluation of Deep Learning Clinical Target Volumes Auto-Contouring for Magnetic Resonance Imaging-Guided Online Adaptive Treatment of Rectal Cancer

Purpose

Segmentation of clinical target volumes (CTV) on medical images can be time-consuming and is prone to interobserver variation (IOV). This is a problem for online adaptive radiation therapy, where CTV segmentation must be performed every treatment fraction, leading to longer treatment times and logistic challenges. Deep learning (DL)-based auto-contouring has the potential to speed up CTV contouring, but its current clinical use is limited. One reason for this is that it can be time-consuming to verify the accuracy of CTV contours produced using auto-contouring, and there is a risk of bias being introduced. To be accepted by clinicians, auto-contouring must be trustworthy. Therefore, there is a need for a comprehensive commissioning framework when introducing DL-based auto-contouring in clinical practice. We present such a framework and apply it to an in-house developed DL model for auto-contouring of the CTV in rectal cancer patients treated with MRI-guided online adaptive radiation therapy.

Methods and Materials

The framework for evaluating DL-based auto-contouring consisted of 3 steps: (1) Quantitative evaluation of the model's performance and comparison with IOV; (2) Expert observations and corrections; and (3) Evaluation of the impact on expected volumetric target coverage. These steps were performed on independent data sets. The framework was applied to an in-house trained nnU-Net model, using the data of 44 rectal cancer patients treated at our institution.

Results

The framework established that the model's performance after expert corrections was comparable to IOV, and although the model introduced a bias, this had no relevant impact on clinical practice. Additionally, we found a substantial time gain without reducing quality as determined by volumetric target coverage.

Conclusions

Our framework provides a comprehensive evaluation of the performance and clinical usability of target auto-contouring models. Based on the results, we conclude that the model is eligible for clinical use.

Online Adaptive MRI-Guided Radiotherapy for Primary Tumor and Lymph Node Boosting in Rectal Cancer

The purpose of this study was to characterize the motion and define the required treatment margins of the pathological mesorectal lymph nodes (GTVln) for two online adaptive MRI-guided strategies for sequential boosting. Secondly, we determine the margins required for the primary gross tumor volume (GTVprim). Twenty-eight patients treated on a 1.5T MR-Linac were included in the study. On T2-weighted images for adaptation (MRIadapt) before and verification after irradiation (MRIpost) of five treatment fractions per patient, the GTVln and GTVprim were delineated. With online adaptive MRI-guided radiotherapy, daily plan adaptation can be performed through the use of two different strategies. In an adapt-to-shape (ATS) workflow the interfraction motion is effectively corrected by redelineation and the only relevant motion is intrafraction motion, while in an adapt-to-position (ATP) workflow the margin (for GTVln) is dominated by interfraction motion. The margin required for GTVprim will be identical to the ATS workflow, assuming each fraction would be perfectly matched on GTVprim. The intrafraction motion was calculated between MRIadapt and MRIpost for the GTVln and GTVprim separately. The interfraction motion of the GTVln was calculated with respect to the position of GTVprim, assuming each fraction would be perfectly matched on GTVprim. PTV margins were calculated for each strategy using the Van Herk recipe. For GTVln we randomly sampled the original dataset 20 times, with each subset containing a single randomly selected lymph node for each patient. The resulting margins for ATS ranged between 3 and 4 mm (LR), 3 and 5 mm (CC) and 5 and 6 mm (AP) based on the 20 randomly sampled datasets for GTVln. For ATP, the margins for GTVln were 10-12 mm in LR and AP and 16-19 mm in CC. The margins for ATS for GTVprim were 1.7 mm (LR), 4.7 mm (CC) and 3.2 mm anterior and 5.6 mm posterior. Daily delineation using ATS of both target volumes results in the smallest margins and is therefore recommended for safe dose escalation to the primary tumor and lymph nodes.

Improving the clinical workflow of a MR-Linac by dosimetric evaluation of synthetic CT

Adaptive radiotherapy performed on the daily magnetic resonance imaging (MRI) is an option to improve the treatment quality. In the adapt-to-shape workflow of 1.5-T MR-Linac, the contours of structures are adjusted on the basis of patient daily MRI, and the adapted plan is recalculated on the MRI-based synthetic computed tomography (syCT) generated by bulk density assignment. Because dosimetric accuracy of this strategy is a priority and requires evaluation, this study aims to explore the usefulness of adding an assessment of dosimetric errors associated with recalculation on syCT to the clinical workflow. Sixty-one patients, with various tumor sites, treated using a 1.5-T MR-Linac were included in this study. In Monaco V5.4, the target and organs at risk (OARs) were contoured, and a reference CT plan that contains information about the outlined contours, their average electron density (ED), and the priority of ED assignment was generated. To evaluate the dosimetric error of syCT caused by the inherent approximation within bulk density assignment, the reference CT plan was recalculated on the syCT obtained from the reference CT by forcing all contoured structures to their mean ED defined on the reference plan. The dose–volume histogram (DVH) and dose distribution of the CT and syCT plan were compared. The causes of dosimetric discrepancies were investigated, and the reference plan was reworked to minimize errors if needed. For 54 patients, gamma analysis of the dose distribution on syCT and CT show a median pass rate of 99.7% and 98.5% with the criteria of 3%/3 mm and 2%/2 mm, respectively. DVH difference of targets and OARs remained less than 1.5% or 1 Gy. For the remaining patients, factors (i.e., inappropriate ED assignments) influenced the dosimetric agreement of the syCT vs. CT reference DVH by up to 21%. The causes of the errors were promptly identified, and the DVH dosimetry was realigned except for two lung treatments for which a significant discrepancy remained. The recalculation on the syCT obtained from the planning CT is a powerful tool to assess and decrease the minimal error committed during the adaptive plan on the MRI-based syCT.

Time-course assessment of 3D-image distortion on the 1.5 T Marlin/Elekta Unity MR-LINAC

PURPOSE

The efficacy of MR-guided radiotherapy on a MR-LINAC (MR-L) is dependent on the geometric accuracy of its MR images over clinically relevant Fields-of-View (FOVs). Our objectives were to: evaluate gradient non-linearity (GNL) on the Elekta Unity MR-L across time via 76 weekly measurements of 3D-distortion over concentrically larger diameter spherical volumes (DSVs); quantify distortion measurement error; and assess the temporal stability of spatial distortion using statistical process control (SPC).

METHODS

MR-image distortion was assessed using a large-FOV 3D-phantom containing 1932 markers embedded in seven parallel plates, spaced 25 mm × 25 mm in- and 55 mm through-plane. Automatically analyzed T1 images yielded distortions in 200, 300, 400 and 500 mm concentric DSVs. Distortion measurement error was evaluated using median absolute difference analysis of imaging repeatability tests.

RESULTS

Over the measurement period absolute time-averaged distortion varied between: dr = 0.30 - 0.49 mm, 0.53 - 0.80 mm, 1.0 - 1.4 mm and 2.28 - 2.37 mm, for DSVs 200, 300, 400 and 500 mm at the 98^th percentile level. Repeatability tests showed that imaging/repositioning introduces negligible error: mean ≤ 0.02 mm (max ≤ 0.3 mm). SPC analysis showed image distortion was stable across all DSVs; however, noticeable changes in GNL were observed following servicing at the one-year mark.

CONCLUSIONS

Image distortion on the MR-L is in the sub-millimeter range for DSVs ≤ 300 mm and stable across time, with SPC analysis indicating all measurements remain within control for each DSV.

Effect of intrafraction adaptation on PTV margins for MRI guided online adaptive radiotherapy for rectal cancer

Purpose

To determine PTV margins for intrafraction motion in MRI-guided online adaptive radiotherapy for rectal cancer and the potential benefit of performing a 2nd adaptation prior to irradiation.

Methods

Thirty patients with rectal cancer received radiotherapy on a 1.5 T MR-Linac. On T2-weighted images for adaptation (MRI_adapt), verification prior to (MRI_ver) and after irradiation (MRI_post) of 5 treatment fractions per patient, the primary tumor GTV (GTV_prim) and mesorectum CTV (CTV_meso) were delineated. The structures on MRI_adapt were expanded to corresponding PTVs. We determined the required expansion margins such that on average over 5 fractions, 98% of CTV_meso and 95% of GTV_prim on MRI_post was covered in 90% of the patients. Furthermore, we studied the benefit of an additional adaptation, just prior to irradiation, by evaluating the coverage between the structures on MRI_ver and MRI_post. A threshold to assess the need for a secondary adaptation was determined by considering the overlap between MRI_adapt and MRI_ver.

Results

PTV margins for intrafraction motion without 2nd adaptation were 6.4 mm in the anterior direction and 4.0 mm in all other directions for CTV_meso and 5.0 mm isotropically for GTV_prim. A 2nd adaptation, applied for all fractions where the motion between MRI_adapt and MRI_ver exceeded 1 mm (36% of the fractions) would result in a reduction of the PTV_meso margin to 3.2 mm/2.0 mm. For PTV_prim a margin reduction to 3.5 mm is feasible when a 2nd adaptation is performed in fractions where the motion exceeded 4 mm (17% of the fractions).

Conclusion

We studied the potential benefit of intrafraction motion monitoring and a 2nd adaptation to reduce PTV margins in online adaptive MRIgRT in rectal cancer. Performing 2nd adaptations immediately after online replanning when motion exceeded 1 mm and 4 mm for CTV_meso and GTV_prim respectively, could result in a 30–50% margin reduction with limited reduction of dose to the bowel.

Reducing systematic errors due to deformation of organs at risk in radiotherapy

PURPOSE

In radiotherapy (RT), the planning CT (pCT) is commonly used to plan the full RT-course. Due to organ deformation and motion, the organ shapes seen at the pCT will not be identical to their shapes during RT. Any difference between the pCT organ shape and the organ's mean shape during RT will cause systematic errors. We propose to use statistical shrinkage estimation to reduce this error using only the pCT and the population mean shape computed from training data.

METHODS

The method was evaluated for the rectum in a cohort of 37 prostate cancer patients that had a pCT and 7-10 treatment CTs with rectum delineations. Deformable registration was performed both within-patient and between patients, resulting in point-to-point correspondence between all rectum shapes, which enabled us to compute a population mean rectum. Shrinkage estimates were found by combining the pCTs linearly with the population mean. The method was trained and evaluated using leave-one-out cross validation. The shrinkage estimates and the patient mean shapes were compared geometrically using the Dice similarity index (DSI), Hausdorff distance (HD), and bidirectional local distance. Clinical dose/volume histograms, equivalent uniform dose (EUD) and minimum dose to the hottest 5% volume (D5%) were compared for the shrinkage estimate and the pCT.

RESULTS

The method resulted in moderate but statistically significant increase in similarity to the patient mean shape over the pCT. On average, the HD was reduced from 15.6 to 13.4 mm, while the DSI was increased from 0.74 to 0.78. Significant reduction in the bias of volume estimates was found in the DVH-range of 52.5-65 Gy, where the bias was reduced from -1.3 to -0.2 percentage points, but no significant improvement was found in EUD or D5%, CONCLUSIONS: The results suggest that shrinkage estimation can reduce systematic errors due to organ deformations in RT. The method has potential to increase the accuracy in RT of deformable organs and can improve motion modeling.

Planning target volume margin assessment for online adaptive MR-guided dose-escalation in rectal cancer on a 1.5 T MR-Linac

PURPOSE

This study assessed the margins needed to cover tumor intrafraction motion during an MR-guided radiotherapy (MRgRT) dose-escalation strategy in intermediate risk rectal cancer.

METHODS

Fifteen patients with rectal cancer were treated with neoadjuvant short-course radiotherapy, 5x5 Gy, according to an online adaptive workflow on a 1.5 T MR-linac. Per patient, 26 3D T2 weighted MRIs were made; one reference scan preceding treatment and five scans per treatment fraction. The primary tumor was delineated on each scan as gross tumor volume (GTV). Target coverage margins were assessed by isotropically expanding the reference GTV until more than 95% of the voxels of the sequential GTVs were covered. A margin with a coverage probability threshold of 90% was defined as adequate. Intra- and interfraction margins to cope with the movement of the GTV in the period between scans were calculated to indicate the target volume margins. Furthermore, the margin needed to cover GTV movement was calculated for different time intervals.

RESULTS

The required margins to cover inter- and intrafraction GTV motion were 17 mm and 6 mm, respectively. Analysis based on time intervals between scans showed smaller margins were needed for adequate GTV coverage as time intervals became shorter, with a 4 mm margin required for a procedure of 15 min or less.

CONCLUSION

The shorter the treatment time, the smaller the margins needed to cover for the GTV movement during an online adaptive MRgRT dose-escalation strategy for intermediate risk rectal cancer. When time intervals between replanning and the end of dose delivery could be reduced to 15 min, a 4 mm margin would allow adequate target coverage.

Technical Note: End-to-end verification of an MR-Linac using a dynamic motion phantom

PURPOSE

MR-Linac integrates an MRI scanner and a linear accelerator to provide adaptive radiation treatment. Superior tissue contrast and real-time imaging can give the clinicians confidence to reduce the margins of the planning target volume (PTV). The purpose of this study was to verify the dosimetric accuracy of an MR-Linac system in treating a moving target and assess the error with different motion patterns and adaptation methods.

METHODS

We performed an end-to-end test for Elekta Unity (Elekta) using the 4D Dynamic Thorax Phantom (CIRS MRgRT 008Z), comparing the measured and planned dose. The moving phantom had four measurement locations in the tumor, liver, kidney, and spinal cord regions with a PTW30013 ion chamber. For seven different motion patterns, we first acquired simulation CT using a slow-scanning protocol, based on which we generated reference plans. The treatment technique was the standard intensity-modulated radiation therapy (IMRT). We tested both adaptation workflows: the Adapt-to-Position (ATP) and the Adapt-to-Shape (ATS). The three-dimensional (3D) distribution was measured using a diode array phantom (Sun Nuclear Inc.) to check the dose distribution accuracy as part of the routine QA process. We also performed end-to-end tests on a conventional Linac. Finally, we used SPSS Statistics 22.0 (Inc., Chicago, IL, USA) for data analysis.

RESULTS

All pretreatment reference plans and delivered plans had excellent QA results with a better than 95% passing rate of relative gamma analysis (2%/2 mm criteria). The adaptive planning for MR-Linac produced quality plans. The measured dose in the target agreed with the calculated dose.

CONCLUSIONS

The adaptive treatment on the MR-Linac system investigated met the expected performance with tumor motions. The outline of the target could be visualized and accurately contoured on the 3D MR for online planning. Under different motion patterns, the difference between the measured and calculated dose was acceptable clinically.

Inter- and intra-fractional stability of rectal gas in pelvic cancer patients during MRIgRT

PURPOSE

Due to the electron return effect (ERE) during magnetic resonance imaging guided radiotherapy (MRIgRT), rectal gas during pelvic treatments can result in hot spots of over-dosage in the rectal wall. Determining the clinical impact of this effect on rectal toxicity requires estimation of the amount and mobility (and stability) of rectal gas during treatment. We therefore investigated the amount of rectal gas and local inter- and intra-fractional changes of rectal gas in pelvic cancer patients.

METHODS

To estimate the volume of gas present at treatment planning, the rectal gas contents in the planning computed tomography (CT) scans of 124 bladder, 70 cervical and 2180 prostate cancer patients were calculated. To estimate inter- and intra-fractional variations in rectal gas, 174 and 131 T2-w MRIs for six cervical and eleven bladder cancer patients were used. These scans were acquired during four scan-sessions (~20-25 min each) at various time-points. Additionally, 258 T2-w MRIs of the first five prostate cancer patients treated using MRIgRT at our center, acquired during each fraction, were analyzed. Rectums were delineated on all scans. The area of gas within the rectum delineations was identified on each MRI slice using thresholding techniques. The area of gas on each slice of the rectum was used to calculate the inter- and intra-fractional group mean, systematic and random variations along the length of the rectum. The cumulative dose perturbation as a result of the gas was estimated. Two approaches were explored: accounting or not accounting for the gas at the start of the scan-session.

RESULTS

Intra-fractional variations in rectal gas are small compared to the absolute volume of rectal gas detected for all patient groups. That is, rectal gas is likely to remain stable for periods of 20-25 min. Larger volumes of gas and larger variations in gas volume were observed in bladder cancer patients compared with cervical and prostate cancer patients. For all patients, local cumulative dose perturbations per beam over an entire treatment in the order of 60 % were estimated when gas had not been accounted for in the daily adaption. The calculated dose perturbation over the whole treatment was dramatically reduced in all patients when accounting for the gas in the daily set-up image.

CONCLUSION

Rectal gas in pelvic cancer patients is likely to remain stable over the course of an MRIgRT fraction, and also likely to reappear in the same location in multiple fractions, and can therefore result in clinically relevant over-dosage in the rectal wall. The over-dosage is reduced when accounting for gas in the daily adaption.

Online adaptive radiotherapy compared to plan selection for rectal cancer: quantifying the benefit

Background

To compare online adaptive radiation therapy (ART) to a clinically implemented plan selection strategy (PS) with respect to dose to the organs at risk (OAR) for rectal cancer.

Methods

The first 20 patients treated with PS between May–September 2016 were included. This resulted in 10 short (SCRT) and 10 long (LCRT) course radiotherapy treatment schedules with a total of 300 Conebeam CT scans (CBCT). New dual arc VMAT plans were generated using auto-planning for both the online ART and PS strategy.

For each fraction bowel bag, bladder and mesorectum were delineated on daily Conebeam CTs. The dose distribution planned was used to calculate daily DVHs. Coverage of the CTV was calculated, as defined by the dose received by 99% of the CTV volume (D99%). The volume of normal tissue irradiated with 95% of the prescribed fraction dose was calculated by calculating the volume receiving 95% of the prescribed fraction or more dose minus the volume of the CTV. For each fraction the difference between the plan selection and online adaptive strategy of each DVH parameter was calculated, as well as the average difference per patient.

Results

Target coverage remained the same for online ART. The median volume of the normal tissue irradiated with 95% of the prescribed dose dropped from 642 cm3 (PS) to 237 cm3 (online-ART)(p < 0.001). Online ART reduced dose to the OARs for all tested dose levels for SCRT and LCRT (p < 0.001). For V15Gy of the bowel bag the median difference over all fractions of all patients was − 126 cm³ in LCRT, while the average difference per patient ranged from − 206 cm³ to − 40 cm³. For SCRT the median difference was − 62 cm³, while the range of the average difference per patient was − 105 cm3 to − 51 cm³.

For V15Gy of the bladder the median difference over all fractions of all patients was 26% in LCRT, while the average difference per patient ranged from − 34 to 12%. For SCRT the median difference of V95% was − 8%, while the range of the average difference per patient was − 29 to 0%.

Conclusions

Online ART for rectal cancer reduces dose the OARs significantly compared to a clinically implemented plan selection strategy, without compromising target coverage.

Trial registration

Medical Research Involving Human Subjects Act (WMO) does not apply to this study and was retrospectively approved by the Medical Ethics review Committee of the Academic Medical Center (W19_357 # 19.420; Amsterdam University Medical Centers, Location Academic Medical Center, Amsterdam, The Netherlands).

Daily adaptive radiotherapy for patients with prostate cancer using a high field MR-linac: Initial clinical experiences and assessment of delivered doses compared to a C-arm linac

INTRODUCTION

MR-guided adapted radiotherapy (MRgART) using a high field MR-linac has recently become available. We report the estimated delivered fractional dose of the first five prostate cancer patients treated at our centre using MRgART and compare this to C-Arm linac daily Image Guided Radiotherapy (IGRT).

METHODS

Patients were treated using adapted treatment plans shaped to their daily anatomy. The treatments were recalculated on an MR image acquired immediately prior to treatment delivery in order to estimate the delivered fractional dose. C-arm linac non-adapted VMAT treatment plans were recalculated on the same MR images to estimate the fractional dose that would have been delivered using conventional radiotherapy techniques using a daily IGRT protocol.

RESULTS

95% and 93% of mandatory target coverage objectives and organ at risk dose constraints were achieved by MRgART and C-arm linac delivered dose estimates, respectively. Both delivery techniques were estimated to have achieved 98% of mandatory Organ At Risk (OAR) dose constraints whereas for the target clinical goals, 86% and 80% were achieved by MRgART and C-arm linac delivered dose estimates.

CONCLUSIONS

Prostate MRgART can be delivered using the a high field MR-linac. Radiotherapy performed on a C-arm linac offers a good solution for prostate cancer patients who present with favourable anatomy at the time of reference imaging and demonstrate stable anatomy throughout the course of their treatment. For patients with critical OARs abutting target volumes on their reference image we have demonstrated the potential for a target dose coverage improvement for MRgART compared to C-arm linac treatment.

Cone-beam computed tomography for organ motion evaluation in locally advanced rectal cancer patients

PURPOSE

Due to a reported dose-response relationship in rectal cancer radiotherapy, a greater interest in dose intensification on small boost volume arises. Considering the need of an appropriate target movements evaluation, this retrospective study aimed to use cone-beam computed tomography (CBCT) for GTV and mesorectum organ motion (OM) evaluation, in locally advanced rectal cancer (LARC) patients treated with neoadjuvant chemo-radiotherapy, in prone and supine position.

METHODS

Thirty-two LARC patients were analyzed. GTV and mesorectum were delineated on MRI co-registrated with CT simulation. GTV and mesorectum OM was estimated on all CBCTs, performed during treatment, co-registrated with CT simulation. OM evaluation was obtained, as mean shift in left and right (L-R), postero-anterior (P-A) and cranio-caudal (Cr-C) directions. Volumes variability was calculated by DICE index.

RESULTS

A total of 296 CBCTs were analyzed. Mean shifts of the GTV and mesorectum in prone position were - 0.16 cm and 0.15 cm in L-R direction, 0.28 cm and - 0.40 cm in P-A direction, and 0.14 cm and - 0.21 cm, in Cr-C direction; for supine position the mean shifts of the GTV were - 0.10 cm and 0.17 cm in R-L direction, 0.26 cm and - 0.23 cm in A-P direction, 0.09 cm and - 0.11 cm in Cr-C direction. Mean DICE index for GTV and mesorectum was 0.74 and 0.86, in prone position, and 0.78 and 0.89 in supine position, respectively.

CONCLUSION

GTV and mesorectum OM was less than 4 mm in all directions in both positions, with a 1 mm less deviation in supine position. CBCTs resulted effective for OM assessment, and it could be an appropriate method for the implementation on an intensification treatment.

Feasibility of Gold Fiducial Markers as a Surrogate for Gross Tumor Volume Position in Image-Guided Radiation Therapy of Rectal Cancer

PURPOSE

To evaluate the feasibility of fiducial markers as a surrogate for gross tumor volume (GTV) position in image-guided radiation therapy of rectal cancer.

METHODS AND MATERIALS

We analyzed 35 fiducials in 19 patients with rectal cancer who received short-course radiation therapy or long-course chemoradiation therapy. Magnetic resonance imaging examinations were performed before and after the first week of radiation therapy, and daily pre- and postirradiation cone beam computed tomography scans were acquired in the first week of radiation therapy. Between the 2 magnetic resonance imaging examinations, the fiducial displacement relative to the center of gravity of the GTV (COG_GTV) and the COG_GTV displacement relative to bony anatomy were determined. Using the cone beam computed tomography scans, inter- and intrafraction fiducial displacement relative to bony anatomy were determined.

RESULTS

The systematic error of the fiducial displacement relative to the COG_GTV was 2.8, 2.4, and 4.2 mm in the left-right, anterior-posterior (AP), and craniocaudal (CC) directions, respectively. Large interfraction systematic errors of up to 8.0 mm and random errors up to 4.7 mm were found for COG_GTV and fiducial displacements relative to bony anatomy, mostly in the AP and CC directions. For tumors located in the mid and upper rectum, these errors were up to 9.4 mm (systematic) and 5.6 mm (random) compared with 4.9 mm and 2.9 mm for tumors in the lower rectum. Systematic and random errors of the intrafraction fiducial displacement relative to bony anatomy were ≤2.1 mm in all directions.

CONCLUSIONS

Large interfraction errors of the COG_GTV and the fiducials relative to bony anatomy were found. Therefore, despite the observed fiducial displacement relative to the COG_GTV, the use of fiducials as a surrogate for GTV position reduces the required margins in the AP and CC directions for a GTV boost using image-guided radiation therapy of rectal cancer. This reduction in margin may be larger in patients with tumors located in the mid and upper rectum compared with the lower rectum.

ReconSocket: a low-latency raw data streaming interface for real-time MRI-guided radiotherapy

With the recent advent of hybrid MRI-guided radiotherapy systems, continuous intra-fraction MR imaging for motion monitoring has become feasible. The ability to perform real-time custom image reconstructions is however often lacking. In this work we present a low-latency streaming solution, ReconSocket, which provides a real-time stream of k-space data from the magnetic resonance imaging (MRI) to custom reconstruction servers. We determined the performance of the data streaming by measuring the streaming latency (i.e. non-zero time delay due to data transfer and processing) and jitter (i.e. deviations from periodicity) using an ultra-fast 1D MRI acquisition of a moving phantom. Simultaneously, its position was recorded with near-zero time delay. The feasibility of low-latency custom reconstructions was tested by measuring the imaging latency (i.e. time delay between physical change and appearance of that change on the image) for several non-Cartesian 2D and 3D acquisitions using an in-house implemented reconstruction server. The measured streaming latency of the ReconSocket interface was [Formula: see text] ms. 98% of the incoming data packets arrived within a jitter range of 367 [Formula: see text]s. This shows that the ReconSocket interface can provide reliable real-time access to MRI data, acquired during the course of a MRI-guided radiotherapy fraction. The total imaging latency was measured to be 221 ms (2D) and 3889 ms (3D) for exemplary acquisitions, using the custom image reconstruction server. These imaging latencies are approximately equal to half of the temporal footprint (T _acq /2) of the respective 2D and 3D golden-angle radial sequences. For radial sequences, it was previously showed that T _acq /2 is the expected contribution of only the data acquisition to the total imaging latency. Indeed, the contribution of the non-Cartesian reconstruction to the total imaging latency was minor (<10%): 21 ms for 2D, 300 ms for 3D, indicating that the acquisition, i.e. the physical encoding of the image itself is the major contributor to the total imaging latency.

Adaptive radiotherapy: The Elekta Unity MR-linac concept

BACKGROUND AND PURPOSE

The promise of the MR-linac is that one can visualize all anatomical changes during the course of radiotherapy and hence adapt the treatment plan in order to always have the optimal treatment. Yet, there is a trade-off to be made between the time spent for adapting the treatment plan against the dosimetric gain. In this work, the various daily plan adaptation methods will be presented and applied on a variety of tumour sites. The aim is to provide an insight in the behavior of the state-of-the-art 1.5 T MRI guided on-line adaptive radiotherapy methods.

MATERIALS AND METHODS

To explore the different available plan adaptation workflows and methods, we have simulated online plan adaptation for five cases with varying levels of inter-fraction motion, regions of interest and target sizes: prostate, rectum, esophagus and lymph node oligometastases (single and multiple target). The plans were evaluated based on the clinical dose constraints and the optimization time was measured.

RESULTS

The time needed for plan adaptation ranged between 17 and 485 s. More advanced plan adaptation methods generally resulted in more plans that met the clinical dose criteria. Violations were often caused by insufficient PTV coverage or, for the multiple lymph node case, a too high dose to OAR in the vicinity of the PTV. With full online replanning it was possible to create plans that met all clinical dose constraints for all cases.

CONCLUSION

Daily full online replanning is the most robust adaptive planning method for Unity. It is feasible for specific sites in clinically acceptable times. Faster methods are available, but before applying these, the specific use cases should be explored dosimetrically.

MRI-based tumor inter-fraction motion statistics for rectal cancer boost radiotherapy

BACKGROUND

In patients diagnosed with rectal cancer, dose escalation is currently being investigated in a large number of studies. Since there is little known on gross tumor volume (GTV) inter-fraction motion for rectal cancer, a wide variety in margins is used. Purpose of this study is to quantify GTV inter-fraction motion statistics on different timescales and to give estimates of planning target volume (PTV) margins.

MATERIAL AND METHODS

Thirty-two patients, diagnosed with rectal cancer, were included. To investigate motion from week-to-week, 16 patients underwent a pretreatment and five weekly MRIs, prior to a radiotherapy (RT) fraction of the chemoradiotherapy treatment. To investigate motion from day-to-day, the remaining 16 patients underwent five daily MRIs before each fraction in one week of RT. GTV was delineated on all scans according to guidelines. Scans were aligned on bony anatomy with the first MRI. For both datasets separately, GTV inter-fraction motion was determined based on center-of-gravity displacement. Therefrom, systematic and random errors were determined in left/right (LR), anterior/posterior and cranial/caudal (CC) direction. PTV margin estimates were calculated and evaluated on GTV coverage.

RESULTS

Systematic and random errors were found in the range of 2.3-4.8 mm and 1.5-3.3 mm from week-to-week, and 1.8-4.5 mm and 1.8-4.0 mm from day-to-day, respectively. On both timescales, similar motion patterns were found; the most motion was observed in CC whilst the least motion was observed in LR. On the week-to-week data more systematic and less random motion was observed compared to the day-to-day data. Overall, only slight differences in margin estimates were found. Derived PTV margin estimates were found to give adequate GTV coverage.

CONCLUSION

GTV inter-fraction motion, on a week-to-week and day-to-day timescale, can be accounted for using motion statistics presented in this study.

Margin and PTV volume reduction using a population based library of plans strategy for rectal cancer radiotherapy

PURPOSE

Day-to-day shape variation in the rectum CTV results in considerable geometric uncertainties during rectal cancer radiotherapy. To ensure coverage a large CTV-to-PTV margin is required. The purpose of this study was to increase the accuracy of treatment delivery by building a population based library of planning CTVs for rectal cancer patients and to evaluate its potential for rectum PTV margin and PTV volume reduction.

METHODS

Analysis was done retrospectively on 33 early-stage rectal cancer patients with daily repeat CTs who received short-course pre-operative radiotherapy in 5 fractions of 5 Gy. We created signed distance maps from the planning rectum CTV to each of the repeat CTVs, from which we calculated the group mean, systematic and random error. The correlation between different regions of the rectum CTV was analyzed and used in combination with the distance maps to create the library of nine planning CTVs. For each of the repeat CTVs the best fitting CTV structure in the library was automatically selected defined by the plan that minimized the mean absolute distance between the repeat and library CTV. Residual distance maps were calculated from which a new PTV margin was constructed. Bootstrapping was performed on the margin difference to assess its significance.

RESULTS

Residual errors were found to decrease with the number of plans in the library, but adding more than five plans yields negligible further error reduction. Margin reduction of up to 50% was achieved at the upper-anterior site of the mesorectum. The average PTV volume decreased by 15.5% when a library is introduced.

CONCLUSIONS

A library of plans strategy for rectal cancer based on population statistics is feasible and results in a considerably reduced average rectum PTV volume compared to conventional radiotherapy.

[Image-guided radiotherapy contribution and patient setup for anorectal cancer treatment]

Intensity-modulated radiation therapy is recommended in anal squamous cell carcinoma treatment and is increasingly used in rectal cancer. It adapts the dose to target volumes, with a high doses gradient. Intensity-modulated radiation therapy allows to reduce toxicity to critical normal structures and to consider dose-escalation studies or systemic treatment intensification. Image-guided radiation therapy is a warrant of quality for intensity-modulated radiation therapy, especially for successful delivery of the dose as planned. There is no recommended international or national anorectal cancer image-guided radiation therapy protocol currently available. Dose-escalation trials or expert opinions about intensity-modulated/image-guided radiation therapy good practice guidelines recommend daily volumetric imaging throughout the treatment or during the five first fractions and weekly thereafter as a minimum. Image-guided radiation therapy allows to reduce margins related to patient setup errors. Internal margin, related to the internal organ motion, needs to be adapted according to short- or long-course radiotherapy, gender, rectal location; it can be higher than current recommended planning target volume margins, particularly in the upper and anterior part of mesorectum, which has the most significant movement. Image-guided radiation therapy based on volumetric imaging allows to take target volume shrinkage into account and to develop adaptive strategies, in particular for mesorectum shrinkage during rectal cancer treatment. Lastly, the emergence of new image-guided radiation therapy technologies including MRI (which plays a major role in pelvic tumours assessment and diagnosis) opens up interesting perspectives for adaptive radiotherapy, taking into account both organs' movements and tumour shrinkage.

Feasibility of MR-only proton dose calculations for prostate cancer radiotherapy using a commercial pseudo-CT generation method

A magnetic resonance (MR)-only radiotherapy workflow can reduce cost, radiation exposure and uncertainties introduced by CT-MRI registration. A crucial prerequisite is generating the so called pseudo-CT (pCT) images for accurate dose calculation and planning. Many pCT generation methods have been proposed in the scope of photon radiotherapy. This work aims at verifying for the first time whether a commercially available photon-oriented pCT generation method can be employed for accurate intensity-modulated proton therapy (IMPT) dose calculation. A retrospective study was conducted on ten prostate cancer patients. For pCT generation from MR images, a commercial solution for creating bulk-assigned pCTs, called MR for Attenuation Correction (MRCAT), was employed. The assigned pseudo-Hounsfield Unit (HU) values were adapted to yield an increased agreement to the reference CT in terms of proton range. Internal air cavities were copied from the CT to minimise inter-scan differences. CT- and MRCAT-based dose calculations for opposing beam IMPT plans were compared by gamma analysis and evaluation of clinically relevant target and organ at risk dose volume histogram (DVH) parameters. The proton range in beam's eye view (BEV) was compared using single field uniform dose (SFUD) plans. On average, a [Formula: see text] mm) gamma pass rate of 98.4% was obtained using a [Formula: see text] dose threshold after adaptation of the pseudo-HU values. Mean differences between CT- and MRCAT-based dose in the DVH parameters were below 1 Gy ([Formula: see text]). The median proton range difference was [Formula: see text] mm, with on average 96% of all BEV dose profiles showing a range agreement better than 3 mm. Results suggest that accurate MR-based proton dose calculation using an automatic commercial bulk-assignment pCT generation method, originally designed for photon radiotherapy, is feasible following adaptation of the assigned pseudo-HU values.

Motion management strategies and technical issues associated with stereotactic body radiotherapy of thoracic and upper abdominal tumors: A review from NRG oncology

The efficacy of stereotactic body radiotherapy (SBRT) has been well demonstrated. However, it presents unique challenges for accurate planning and delivery especially in the lungs and upper abdomen where respiratory motion can be significantly confounding accurate targeting and avoidance of normal tissues. In this paper, we review the current literature on SBRT for lung and upper abdominal tumors with particular emphasis on addressing respiratory motion and its affects. We provide recommendations on strategies to manage motion for different, patient-specific situations. Some of the recommendations will potentially be adopted to guide clinical trial protocols.