Automated tracking of morphologic changes in weekly magnetic resonance imaging during head and neck radiotherapy.

BACKGROUND AND PURPOSE: Anatomic changes during head and neck radiotherapy can impact dose delivery, necessitate adaptive replanning, and indicate patient-specific response to treatment. We have developed an automated system to track these changes through longitudinal MRI scans to aid identification and clinical intervention. The purpose of this article is to describe this tracking system and present results from an initial cohort of patients. MATERIALS AND METHODS: The Automated Watchdog in Adaptive Radiotherapy Environment (AWARE) was developed to process longitudinal MRI data for radiotherapy patients. AWARE automatically identifies and collects weekly scans, propagates radiotherapy planning structures, computes structure changes over time, and reports important trends to the clinical team. AWARE also incorporates manual structure review and revision from clinical experts and dynamically updates tracking statistics when necessary. AWARE was applied to patients receiving weekly T2-weighted MRI scans during head and neck radiotherapy. Changes in nodal gross tumor volume (GTV) and parotid gland delineations were tracked over time to assess changes during treatment and identify early indicators of treatment response. RESULTS: N = 91 patients were tracked and analyzed in this study. Nodal GTVs and parotids both shrunk considerably throughout treatment (-9.7 ± 7.7% and -3.7 ± 3.3% per week, respectively). Ipsilateral parotids shrunk significantly faster than contralateral (-4.3 ± 3.1% vs. -2.9 ± 3.3% per week, p = 0.005) and increased in distance from GTVs over time (+2.7 ± 7.2% per week, p < 1 × 10-5 ). Automatic structure propagations agreed well with manual revisions (Dice = 0.88 ± 0.09 for parotids and 0.80 ± 0.15 for GTVs), but for GTVs the agreement degraded 4-5 weeks after the start of treatment. Changes in GTV volume observed by AWARE as early as one week into treatment were predictive of large changes later in the course (AUC = 0.79). CONCLUSION: AWARE automatically identified longitudinal changes in GTV and parotid volumes during radiotherapy. Results suggest that this system may be useful for identifying rapidly responding patients as early as one week into treatment.

Accuracy of surface-guided patient setup for conventional radiotherapy of brain and nasopharynx cancer.

PURPOSE: To evaluate the accuracy of surface-guided radiotherapy (SGRT) in cranial patient setup by direct comparison between optical surface imaging (OSI) and cone-beam computed tomography (CBCT), before applying SGRT-only setup for conventional radiotherapy of brain and nasopharynx cancer. METHODS AND MATERIALS: Using CBCT as reference, SGRT setup accuracy was examined based on 269 patients (415 treatments) treated with frameless cranial stereotactic radiosurgery (SRS) during 2018-2019. Patients were immobilized in customized head molds and open-face masks and monitored using OSI during treatment. The facial skin area in planning CT was used as OSI region of interest (ROI) for automatic surface alignment and the skull was used as the landmark for automatic CBCT/CT registration. A 6 degrees of freedom (6DOF) couch was used. Immediately after CBCT setup, an OSI verification image was captured, recording the SGRT setup differences. These differences were analyzed in 6DOFs and as a function of isocenter positions away from the anterior surface to assess OSI-ROI bias. The SGRT in-room setup time was estimated and compared with CBCT and orthogonal 2D kilovoltage (2DkV) setups. RESULTS: The SGRT setup difference (magnitude) is found to be 1.0 ± 2.5 mm and 0.1˚±1.4˚ on average among 415 treatments and within 5 mm/3˚ with greater than 95% confidence level (P < 0.001). Outliers were observed for very-posterior isocenters: 15 differences (3.6%) are >5.0mm and 9 (2.2%) are >3.0˚. The setup differences show minor correlations (|r| < 0.45) between translational and rotational DOFs and a minor increasing trend (<1.0 mm) in the anterior-to-posterior direction. The SGRT setup time is 0.8 ± 0.3 min, much shorter than CBCT (5 ± 2 min) and 2DkV (2 ± 1 min) setups. CONCLUSION: This study demonstrates that SGRT has sufficient accuracy for fast in-room patient setup and allows real-time motion monitoring for beam holding during treatment, potentially useful to guide radiotherapy of brain and nasopharynx cancer with standard fractionation.

Illustrated instructions for mechanical quality assurance of a medical linear accelerator.

PURPOSE: The purpose of this study was to develop and test a set of illustrated instructions for effective training for mechanical quality assurance (QA) of medical linear accelerators (linac). METHODS: Illustrated instructions were created for mechanical QA and underwent several steps of review, testing, and refinement. Eleven testers with no recent QA experience were then recruited from our radiotherapy department (one student, two computational scientists, and eight dosimetrists). This group was selected because they have experience of radiation therapy but no preconceived ideas about how to do QA. The following parameters were progressively decalibrated on a Varian C-series linac: Group A = gantry angle, ceiling laser position, X1 jaw position, couch longitudinal position, physical graticule position (five testers); Group B = Group A + wall laser position, couch lateral and vertical position, collimator angle (three testers); Group C = Group B + couch angle, wall laser angle, and optical distance indicator (three testers). Testers were taught how to use the linac and then used the instructions to try to identify these errors. An experienced physicist observed each session, giving support on machine operation as necessary. RESULTS: Testers were able to follow the instructions. They determined gantry, collimator, and couch angle errors within 0.4°, 0.3°, and 0.9° of the actual changed values, respectively. Laser positions were determined within 1 mm and jaw positions within 2 mm. Couch position errors were determined within 2 mm and 3 mm for lateral/longitudinal and vertical errors, respectively. Accessory-positioning errors were determined within 1 mm. Optical distance indicator errors were determined within 2 mm when comparing with distance sticks and 6 mm when using blocks, indicating that distance sticks should be the preferred approach for inexperienced staff. CONCLUSIONS: Inexperienced users were able to follow these instructions and catch errors within the criteria suggested by AAPM TG-142 for linacs used for intensity-modulated radiation therapy. These instructions are, therefore, suitable for QA training.

Initial clinical experience with ArcCHECK for IMRT/VMAT QA.

Many devices designed for the purpose of performing patient-specific IMRT/VMAT QA are commercially available. In this work we report our experience and initial clinical results with the ArcCHECK. The ArcCHECK consists of a cylindrical array of diode detectors measuring entry and exit doses. The measured result is a cumulative dose displayed as a 2D matrix. The detector array requires both an absolute dose calibration, and a calibration of the detector response, relative to each other. In addition to the calibrations suggested by the manufacturer, various tests were performed in order to assess its stability and performance prior to clinical introduction. Tests of uniformity, linearity, and repetition rate dependence of the detector response were conducted and described in this work. Following initial test-ing, the ArcCHECK device was introduced in the clinic for routine patient-specific IMRT QA. The clinical results from one year of use were collected and analyzed. The gamma pass rates at the 3%/3 mm criterion were reported for 3,116 cases that included both IMRT and VMAT treatment plans delivered on 18 linear accelera-tors. The gamma pass rates were categorized based on the treatment site, treatment technique, type of MLCs, operator, ArcCHECK device, and LINAC model. We recorded the percent of failures at the clinically acceptable threshold of 90%. In addition, we calculated the threshold that encompasses two standard deviations (2 SD) (95%) of QAs (T95) for each category investigated. The commissioning measurements demonstrated that the device performed as expected. The uniformity of the detector response to a constant field arc delivery showed a 1% standard deviation from the mean. The variation in dose with changing repetition rate was within 1 cGy of the mean, while the measured dose showed a linear relation with delivered MUs. Our initial patient QA results showed that, at the clinically selected passing criterion, 4.5% of cases failed. On average T95 was 91%, rang-ing from 73% for gynecological sites to 96.5% for central nervous system sites. There are statistically significant differences in passing rates between IMRT and VMAT, high-definition (HD) and non-HD MLCs, and different LINAC models (p-values << 0.001). An additional investigation into the failing QAs and a com-parison with ion-chamber measurements reveals that the differences observed in the passing rates between the different studied factors can be largely explained by the field size dependence of the device. Based on our initial experience with the ArcCHECK, our passing rates are, on average, consistent with values reported in the AAPM TG-119. However, the significant variations between QAs that were observed based on the size of the treatment fields may need to be corrected to improve the specificity and sensitivity of the device.

Artificial-intelligence-driven measurements of brain metastases' response to SRS compare favorably with current manual standards of assessment.

Background: Evaluation of treatment response for brain metastases (BMs) following stereotactic radiosurgery (SRS) becomes complex as the number of treated BMs increases. This study uses artificial intelligence (AI) to track BMs after SRS and validates its output compared with manual measurements. Methods: Patients with BMs who received at least one course of SRS and followed up with MRI scans were retrospectively identified. A tool for automated detection, segmentation, and tracking of intracranial metastases on longitudinal imaging, MEtastasis Tracking with Repeated Observations (METRO), was applied to the dataset. The longest three-dimensional (3D) diameter identified with METRO was compared with manual measurements of maximum axial BM diameter, and their correlation was analyzed. Change in size of the measured BM identified with METRO after SRS treatment was used to classify BMs as responding, or not responding, to treatment, and its accuracy was determined relative to manual measurements. Results: From 71 patients, 176 BMs were identified and measured with METRO and manual methods. Based on a one-to-one correlation analysis, the correlation coefficient was R2 = 0.76 (P = .0001). Using modified BM response classifications of BM change in size, the longest 3D diameter data identified with METRO had a sensitivity of 0.72 and a specificity of 0.95 in identifying lesions that responded to SRS, when using manual axial diameter measurements as the ground truth. Conclusions: Using AI to automatically measure and track BM volumes following SRS treatment, this study showed a strong correlation between AI-driven measurements and the current clinically used method: manual axial diameter measurements.

Multi-modality imaging parameters that predict rapid tumor regression in head and neck radiotherapy.

Background and purpose: Volume regression during radiotherapy can indicate patient-specific treatment response. We aimed to identify pre-treatment multimodality imaging (MMI) metrics from positron emission tomography (PET), magnetic resonance imaging (MRI), and computed tomography (CT) that predict rapid tumor regression during radiotherapy in human papilloma virus (HPV) associated oropharyngeal carcinoma. Materials and methods: Pre-treatment FDG PET-CT, diffusion-weighted MRI (DW-MRI), and intra-treatment (at 1, 2, and 3 weeks) MRI were acquired in 72 patients undergoing chemoradiation therapy for HPV+ oropharyngeal carcinoma. Nodal gross tumor volumes were delineated on longitudinal images to measure intra-treatment volume changes. Pre-treatment PET standardized uptake value (SUV), CT Hounsfield Unit (HU), and non-gaussian intravoxel incoherent motion DW-MRI metrics were computed and correlated with volume changes. Intercorrelations between MMI metrics were also assessed using network analysis. Validation was carried out on a separate cohort (N = 64) for FDG PET-CT. Results: Significant correlations with volume loss were observed for baseline FDG SUVmean (Spearman ρ = 0.46, p < 0.001), CT HUmean (ρ = 0.38, p = 0.001), and DW-MRI diffusion coefficient, Dmean (ρ = -0.39, p < 0.001). Network analysis revealed 41 intercorrelations between MMI and volume loss metrics, but SUVmean remained a statistically significant predictor of volume loss in multivariate linear regression (p = 0.01). Significant correlations were also observed for SUVmean in the validation cohort in both primary (ρ = 0.30, p = 0.02) and nodal (ρ = 0.31, p = 0.02) tumors. Conclusions: Multiple pre-treatment imaging metrics were correlated with rapid nodal gross tumor volume loss during radiotherapy. FDG-PET SUV in particular exhibited significant correlations with volume regression across the two cohorts and in multivariate analysis.

Clinical Experience With an Offline Adaptive Radiation Therapy Head and Neck Program: Dosimetric Benefits and Opportunities for Patient Selection.

PURPOSE: The objective of this study was to develop a linear accelerator (LINAC)-based adaptive radiation therapy (ART) workflow for the head and neck that is informed by automated image tracking to identify major anatomic changes warranting adaptation. In this study, we report our initial clinical experience with the program and an investigation into potential trigger signals for ART. METHODS AND MATERIALS: Offline ART was systematically performed on patients receiving radiation therapy for head and neck cancer on C-arm LINACs. Adaptations were performed at a single time point during treatment with resimulation approximately 3 weeks into treatment. Throughout treatment, all patients were tracked using an automated image tracking system called the Automated Watchdog for Adaptive Radiotherapy Environment (AWARE). AWARE measures volumetric changes in gross tumor volumes (GTVs) and selected normal tissues via cone beam computed tomography scans and deformable registration. The benefit of ART was determined by comparing adaptive plan dosimetry and normal tissue complication probabilities against the initial plans recalculated on resimulation computed tomography scans. Dosimetric differences were then correlated with AWARE-measured volume changes to identify patient-specific triggers for ART. Candidate trigger variables were evaluated using receiver operator characteristic analysis. RESULTS: In total, 46 patients received ART in this study. Among these patients, we observed a significant decrease in dose to the submandibular glands (mean ± standard deviation: -219.2 ± 291.2 cGy, P < 10-5), parotids (-68.2 ± 197.7 cGy, P = .001), and oral cavity (-238.7 ± 206.7 cGy, P < 10-5) with the adaptive plan. Normal tissue complication probabilities for xerostomia computed from mean parotid doses also decreased significantly with the adaptive plans (P = .008). We also observed systematic intratreatment volume reductions (ΔV) for GTVs and normal tissues. Candidate triggers were identified that predicted significant improvement with ART, including parotid ΔV = 7%, neck ΔV = 2%, and nodal GTV ΔV = 29%. CONCLUSIONS: Systematic offline head and neck ART was successfully deployed on conventional LINACs and reduced doses to critical salivary structures and the oral cavity. Automated cone beam computed tomography tracking provided information regarding anatomic changes that may aid patient-specific triggering for ART.

Radiation Therapy for HPV-Positive Oropharyngeal Squamous Cell Carcinoma: An ASTRO Clinical Practice Guideline.

PURPOSE: Human Papilloma Virus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC) is a distinct disease from other head and neck tumors. This guideline provides evidence-based recommendations on the critical decisions in its curative treatment, including both definitive and postoperative radiation therapy (RT) management. METHODS: ASTRO convened a task force to address 5 key questions on the use of RT for management of HPV-associated OPSCC. These questions included indications for definitive and postoperative RT and chemoradiation; dose-fractionation regimens and treatment volumes; preferred RT techniques and normal tissue considerations; and posttreatment management decisions. The task force did not address indications for primary surgery versus RT. Recommendations were based on a systematic literature review and created using a predefined consensus-building methodology and system for grading evidence quality and recommendation strength. RESULTS: Concurrent cisplatin is recommended for patients receiving definitive RT with T3-4 disease and/or 1 node >3 cm, or multiple nodes. For similar patients who are ineligible for cisplatin, concurrent cetuximab, carboplatin/5-fluorouracil, or taxane-based systemic therapy are conditionally recommended. In the postoperative setting, RT with concurrent cisplatin (either schedule) is recommended for positive surgical margins or extranodal extension. Postoperative RT alone is recommended for pT3-4 disease, >2 nodes, or a single node >3 cm. Observation is conditionally recommended for pT1-2 disease and a single node ≤3 cm without other risk factors. For patients treated with definitive RT with concurrent systemic therapy, 7000 cGy in 33 to 35 fractions is recommended, and for patients receiving postoperative RT without positive surgical margins and extranodal extension, 5600 to 6000 cGy is recommended. For all patients receiving RT, intensity modulated RT over 3-dimensional techniques with reduction in dose to critical organs at risk (including salivary and swallowing structures) is recommended. Reassessment with positron emission tomography-computed tomography is recommended approximately 3 months after definitive RT/chemoradiation, and neck dissection is recommended for convincing evidence of residual disease; for equivocal positron emission tomography-computed tomography findings, either neck dissection or repeat imaging is recommended. CONCLUSIONS: The role and practice of RT continues to evolve for HPV-associated OPSCC, and these guidelines inform best clinical practice based on the available evidence.

Automatically tracking brain metastases after stereotactic radiosurgery.

Background and purpose: Patients with brain metastases (BMs) are surviving longer and returning for multiple courses of stereotactic radiosurgery. BMs are monitored after radiation with follow-up magnetic resonance (MR) imaging every 2-3 months. This study investigated whether it is possible to automatically track BMs on longitudinal imaging and quantify the tumor response after radiotherapy. Methods: The METRO process (MEtastasis Tracking with Repeated Observations was developed to automatically process patient data and track BMs. A longitudinal intrapatient registration method for T1 MR post-Gd was conceived and validated on 20 patients. Detections and volumetric measurements of BMs were obtained from a deep learning model. BM tracking was validated on 32 separate patients by comparing results with manual measurements of BM response and radiologists' assessments of new BMs. Linear regression and residual analysis were used to assess accuracy in determining tumor response and size change. Results: A total of 123 irradiated BMs and 38 new BMs were successfully tracked. 66 irradiated BMs were visible on follow-up imaging 3-9 months after radiotherapy. Comparing their longest diameter changes measured manually vs. METRO, the Pearson correlation coefficient was 0.88 (p < 0.001); the mean residual error was -8 ± 17%. The mean registration error was 1.5 ± 0.2 mm. Conclusions: Automatic, longitudinal tracking of BMs using deep learning methods is feasible. In particular, the software system METRO fulfills a need to automatically track and quantify volumetric changes of BMs prior to, and in response to, radiation therapy.

A Potential Pitfall and Clinical Solutions in Surface-Guided Deep Inspiration Breath Hold Radiation Therapy for Left-Sided Breast Cancer.

Purpose: Deep inspiration breath hold (DIBH) is an effective technique to spare the heart in treating left-sided breast cancer. Surface-guided radiation therapy (SGRT) is increasingly applied in DIBH setup and motion monitoring. Patient-specific breathing behavior, either thoracically driven or abdominally driven (A-DIBH), should be unaltered, online identified, and monitored accordingly to ensure reproducible heart-sparing treatment. Methods and Materials: Sixty patients with left-sided breast cancer treated with SGRT were analyzed: 20 A-DIBH patients with vertical chest elevation (VCE ≤ 5 mm) were prospectively identified, and 40 control patients were retrospectively and randomly selected for comparison. At simulation, both free-breathing (FB) and DIBH computed tomography (CT) were acquired, guided by a motion surrogate placed around the xiphoid process. For SGRT treatment setups, the region of interest (ROI) was defined on the CT chest surface, and the surrogate-based setup was a backup. For all 60 patients, the VCE was measured as the average of the FB-to-DIBH elevations at the breast and xiphoid process, together with abdominal elevation. In the 40-patient control group, A-DIBH patients (VCE ≤ 5 mm) were identified. Of the 20 A-DIBH patients, 10 were treated with volumetric modulated arc therapy plans, and 10 patients were treated with tangent plans. Clinical DIBH plans were recalculated on FB CT to compare maximum dose (DMax), 5% of the maximum dose (D5%), mean dose (DMean), and V30Gy, V20Gy, and V5Gy of the heart and lungs and their significance. Results: In the 20 A-DIBH patients, VCE = 3 ± 2 mm, surrogate motion (9 ± 6 mm), and abdomen motion of 14 ± 5 mm are found. Heart dose reduction from FB to DIBH is significant (P < .01): ∆DMax = -8.4 ± 9.8 Gy, ∆D5% = -2.4 ± 4.4 Gy, and ∆DMean = -0.6 ± 0.9 Gy. Six out of 40 control patients (15%) are found to have VCE ≤ 5 mm. Conclusions: A-DIBH (VCE ≤ 5 mm) patient population is significant (15%), and they should be identified in the SGRT workflow and monitored accordingly. A new abdominal ROI or an abdominal surrogate should be used instead of the conventional chest-only ROI. Patient-specific DIBH should be preserved for higher reproducibility to ensure heart sparing.

3-D fiducial motion tracking using limited MV projections in arc therapy.

PURPOSE: In-treatment fiducial tracking has recently received attention as a method for improving treatment accuracy, dose conformity, and sparing of healthy tissue. 3-D fiducial localization in arc-radiotherapy remains challenging due to the motion of targets and the complexity of arc deliveries. We propose a novel statistical method for estimating 3-D fiducial motion using limited 2-D megavoltage (MV) projections. METHODS: 3-D fiducial motion was estimated by a maximum a posteriori (MAP) approach to integrating information of fiducial projections with prior knowledge of target motion. To obtain the imaging geometries, short sequences of MV projections were selected in which fiducials were continuously visible. The MAP algorithm estimated the 3-D motion by maximizing the probability of displacement of fiducials in the sequences. Prior knowledge of target motion from a large statistical sample was built into the model to enhance the accuracy of estimation. In the case that a motion prior was unavailable, the algorithm can be simplified to the maximum likelihood (ML) approach. To compare tracking performance, a multiprojection geometric method was also presented by extending the typical two-project ion geometric estimation approach. The algorithms were evaluated using clinical prostate motion traces, and the performance was measured in quality indices and statistical evaluation. RESULTS: The results showed that the MAP method significantly outperforms the geometric method in all measures. In our simulations, the MAP method achieved an accuracy of less than 1 mm RMS error using only five continuous projections, whereas the geometric method required 15 projections to achieve a similar result. CONCLUSIONS: The approach presented can accurately estimate tumor motion using a limited number of continuous projections. The MAP motion estimation is superior to both the ML and geometric estimation methods.

Automatic segmentation of brain metastases using T1 magnetic resonance and computed tomography images.

An increasing number of patients with multiple brain metastases are being treated with stereotactic radiosurgery (SRS). Manually identifying and contouring all metastatic lesions is difficult and time-consuming, and a potential source of variability. Hence, we developed a 3D deep learning approach for segmenting brain metastases on MR and CT images. Five-hundred eleven patients treated with SRS were retrospectively identified for this study. Prior to radiotherapy, the patients were imaged with 3D T1 spoiled-gradient MR post-Gd (T1 + C) and contrast-enhanced CT (CECT), which were co-registered by a treatment planner. The gross tumor volume contours, authored by the attending radiation oncologist, were taken as the ground truth. There were 3 ± 4 metastases per patient, with volume up to 57 ml. We produced a multi-stage model that automatically performs brain extraction, followed by detection and segmentation of brain metastases using co-registered T1 + C and CECT. Augmented data from 80% of these patients were used to train modified 3D V-Net convolutional neural networks for this task. We combined a normalized boundary loss function with soft Dice loss to improve the model optimization, and employed gradient accumulation to stabilize the training. The average Dice similarity coefficient (DSC) for brain extraction was 0.975 ± 0.002 (95% CI). The detection sensitivity per metastasis was 90% (329/367), with moderate dependence on metastasis size. Averaged across 102 test patients, our approach had metastasis detection sensitivity 95 ± 3%, 2.4 ± 0.5 false positives, DSC of 0.76 ± 0.03, and 95th-percentile Hausdorff distance of 2.5 ± 0.3 mm (95% CIs). The volumes of automatic and manual segmentations were strongly correlated for metastases of volume up to 20 ml (r=0.97,p<0.001). This work expounds a fully 3D deep learning approach capable of automatically detecting and segmenting brain metastases using co-registered T1 + C and CECT.

Deep learning auto-segmentation and automated treatment planning for trismus risk reduction in head and neck cancer radiotherapy.

BACKGROUND AND PURPOSE: Reducing trismus in radiotherapy for head and neck cancer (HNC) is important. Automated deep learning (DL) segmentation and automated planning was used to introduce new and rarely segmented masticatory structures to study if trismus risk could be decreased. MATERIALS AND METHODS: Auto-segmentation was based on purpose-built DL, and automated planning used our in-house system, ECHO. Treatment plans for ten HNC patients, treated with 2 Gy × 35 fractions, were optimized (ECHO0). Six manually segmented OARs were replaced with DL auto-segmentations and the plans re-optimized (ECHO1). In a third set of plans, mean doses for auto-segmented ipsilateral masseter and medial pterygoid (MIMean, MPIMean), derived from a trismus risk model, were implemented as dose-volume objectives (ECHO2). Clinical dose-volume criteria were compared between the two scenarios (ECHO0 vs. ECHO1; ECHO1 vs. ECHO2; Wilcoxon signed-rank test; significance: p < 0.01). RESULTS: Small systematic differences were observed between the doses to the six auto-segmented OARs and their manual counterparts (median: ECHO1 = 6.2 (range: 0.4, 21) Gy vs. ECHO0 = 6.6 (range: 0.3, 22) Gy; p = 0.007), and the ECHO1 plans provided improved normal tissue sparing across a larger dose-volume range. Only in the ECHO2 plans, all patients fulfilled both MIMean and MPIMean criteria. The population median MIMean and MPIMean were considerably lower than those suggested by the trismus model (ECHO0: MIMean = 13 Gy vs. ≤42 Gy; MPIMean = 29 Gy vs. ≤68 Gy). CONCLUSIONS: Automated treatment planning can efficiently incorporate new structures from DL auto-segmentation, which results in trismus risk sparing without deteriorating treatment plan quality. Auto-planning and deep learning auto-segmentation together provide a powerful platform to further improve treatment planning.

Use of a realistic breathing lung phantom to evaluate dose delivery errors.

PURPOSE: To compare the effect of respiration-induced motion on delivered dose (the interplay effect) for different treatment techniques under realistic clinical conditions. METHODS: A flexible resin tumor model was created using rapid prototyping techniques based on a computed tomography (CT) image of an actual tumor. Twenty micro-MOSFETs were inserted into the tumor model and the tumor model was inserted into an anthropomorphic breathing phantom. Phantom motion was programed using the motion trajectory of an actual patient. A four-dimensional CT image was obtained and several treatment plans were created using different treatment techniques and planning systems: Conformal (Eclipse), step-and-shoot intensity-modulated radiation therapy (IMRT) (Pinnacle), step-and-shoot IMRT (XiO), dynamic IMRT (Eclipse), complex dynamic IMRT (Eclipse), hybrid IMRT [60% conformal, 40% dynamic IMRT (Eclipse)], volume-modulated are therapy (VMAT) [single-arc (Eclipse)], VMAT [double-arc (Eclipse)], and complex VMAT (Eclipse). The complex plans were created by artificially pushing the optimizer to give complex multileaf collimator sequences. Each IMRT field was irradiated five times and each VMAT field was irradiated ten times, with each irradiation starting at a random point in the respiratory cycle. The effect of fractionation was calculated by randomly summing the measured doses. The maximum deviation for each measurement point per fraction and the probability that 95% of the model tumor had dose deviations less than 2% and 5% were calculated as a function of the number of fractions. Tumor control probabilities for each treatment plan were calculated and compared. RESULTS: After five fractions, measured dose deviations were less than 2% for more than 95% of measurement points within the tumor model for all plans, except the complex dynamic IMRT, step-and-shoot IMRT (XiO), complex VMAT, and single-arc VMAT plans. Reducing the dose rate of the complex IMRT plans from 600 to 200 MU/min reduced the dose deviations to less than 2%. Dose deviations were less than 5% after five fractions for all plans, except the complex single-arc VMAT plan. CONCLUSIONS: Rapid prototyping techniques can be used to create realistic tumor models. For most treatment techniques, the dose deviations averaged out after several fractions. Treatments with unusually complicated multileaf collimator sequences had larger dose deviations. For IMRT treatments, dose deviations can be reduced by reducing the dose rate. For VMAT treatments, using two arcs instead of one is effective for reducing dose deviations.

Predictive dose accumulation for HN adaptive radiotherapy.

During radiation therapy (RT) of head and neck (HN) cancer, the shape and volume of the parotid glands (PG) may change significantly, resulting in clinically relevant deviations of delivered dose from the planning dose. Early and accurate longitudinal prediction of PG anatomical changes during the RT can be valuable to inform decisions on plan adaptation. We developed a deep neural network for longitudinal predictions using the displacement fields (DFs) between the planning computed tomography (pCT) and weekly cone beam computed tomography (CBCT). Sixty-three HN patients treated with volumetric modulated arc were retrospectively studied. We calculated DFs between pCT and week 1-3 CBCT by B-spline and Demon deformable image registration (DIR). The resultant DFs were subsequently used as input to our novel network to predict the week 4 to 6 DFs for generating predicted weekly PG contours and weekly dose distributions. For evaluation, we measured dice similarity (DICE), and the uncertainty of accumulated dose. Moreover, we compared the detection accuracies of candidates for adaptive radiotherapy (ART) when the trigger criteria were mean dose difference more than 10%, 7.5%, and 5%, respectively. The DICE of ipsilateral/contralateral PG at week 4 to 6 using the prediction model trained with B-spline were 0.81 [Formula: see text] 0.07/0.81 [Formula: see text] 0.04 (week 4), 0.79 [Formula: see text] 0.06/0.81 [Formula: see text] 0.05 (week 5) and 0.78 [Formula: see text] 0.06/0.82 [Formula: see text] (week 6). The DICE with the Demons model were 0.78 [Formula: see text] 0.08/0.82 [Formula: see text] 0.03 (week 4), 0.77 [Formula: see text] 0.07/0.82 [Formula: see text] 0.04 (week 5) and 0.75 [Formula: see text] 0.07/0.82 [Formula: see text] 0.02 (week 6). The dose volume histogram (DVH) analysis with the predicted accumulated dose showed the feasibility of predicting dose uncertainty due to the PG anatomical changes. The AUC of ART candidate detection with our predictive model was over 0.90. In conclusion, the proposed network was able to predict future anatomical changes and dose uncertainty of PGs with clinically acceptable accuracy, and hence can be readily integrated into the ART workflow.

Differences between planned and delivered dose for head and neck cancer, and their consequences for normal tissue complication probability and treatment adaptation.

BACKGROUND AND PURPOSE: Anatomical changes induce differences between planned and delivered dose. Adaptive radiotherapy (ART) may reduce these differences but the optimal implementation is insufficiently clear. The aims of this study were to quantify the difference between planned and delivered dose in HNC patients, assess the consequential difference in normal tissue complication probability (ΔNTCP) and to explore the value of ΔNTCP as an objective selection strategy for ART. MATERIALS AND METHODS: For 52 patients, daily doses were accumulated to estimate the delivered dose. The difference from planned dose was analyzed for CTVs and 9 organs-at-risk (OAR). ΔNTCP was calculated for xerostomia, dysphagia, parotid gland dysfunction and tube feeding dependency at 6 months. ART was deemed necessary if ΔNTCP was >5%. The positive predicted value (PPV) was calculated for identification of ART-patients by clinical judgement, and ΔNTCP at fraction 10 and 15. RESULTS: ΔNTCP >5% was seen five times for dysphagia and twice for the other toxicities. Only 5/9 patients with any ΔNTCP >5% clinically received ART, although ART had been done for 13/52 patients (PPV: 0.38). PPV was 0.86 and 0.75 for accumulated dose at fraction 10 and 15, respectively, using a ΔNTCP cut-off for the allocation of ART of 5%. Using other ΔNTCP cut-offs did not substantially improve PPV. With this cut-off the negative predictive value was 0.93 for ΔNTCP method of fraction 10 and fraction 15, and 0.90 for clinical judgement. CONCLUSION: To identify patients accurately for ART, NTCP calculations based on the dose differences between planned and delivered dose at fraction 10 are superior to clinical judgement.

A PET Radiomics Model to Predict Refractory Mediastinal Hodgkin Lymphoma.

First-order radiomic features, such as metabolic tumor volume (MTV) and total lesion glycolysis (TLG), are associated with disease progression in early-stage classical Hodgkin lymphoma (HL). We hypothesized that a model incorporating first- and second-order radiomic features would more accurately predict outcome than MTV or TLG alone. We assessed whether radiomic features extracted from baseline PET scans predicted relapsed or refractory disease status in a cohort of 251 patients with stage I-II HL who were managed at a tertiary cancer center. Models were developed and tested using a machine-learning algorithm. Features extracted from mediastinal sites were highly predictive of primary refractory disease. A model incorporating 5 of the most predictive features had an area under the curve (AUC) of 95.2% and total error rate of 1.8%. By comparison, the AUC was 78% for both MTV and TLG and was 65% for maximum standardize uptake value (SUVmax). Furthermore, among the patients with refractory mediastinal disease, our model distinguished those who were successfully salvaged from those who ultimately died of HL. We conclude that our PET radiomic model may improve upfront stratification of early-stage HL patients with mediastinal disease and thus contribute to risk-adapted, individualized management.

The development and testing of a digital PET phantom for the evaluation of tumor volume segmentation techniques.

Methods for accurate tumor volume segmentation of positron emission tomography (PET) images have been under investigation in recent years partly as a result of the increased use of PET in radiation treatment planning (RTP). None of the developed automated or semiautomated segmentation methods, however, has been shown reliable enough to be regarded as the standard. One reason for this is that there is no source of well characterized and reliable test data for evaluating such techniques. The authors have constructed a digital tumor phantom to address this need. The phantom was created using the Zubal phantom as input to the SimSET software used for PET simulations. Synthetic tumors were placed in the lung of the Zubal phantom to provide the targets for segmentation. The authors concentrated on the lung, since much of the interest to include PET in RTP is for nonsmall cell lung cancer. Several tests were performed on the phantom to ensure its close resemblance to clinical PET scans. The authors measured statistical quantities to compare image intensity distributions from regions-of-interest (ROIs) placed in the liver, the lungs, and tumors in phantom and clinical reconstructions. Using ROIs they also made measurements of autocorrelation functions to ensure the image texture is similar in clinical and phantom data. The authors also compared the intensity profile and appearance of real and simulated uniform activity spheres within uniform background. These measurements, along with visual comparisons of the phantom with clinical scans, indicate that the simulated phantom mimics reality quite well. Finally, they investigate and quantify the relationship between the threshold required to segment a tumor and the inhomogeneity of the tumor's image intensity distribution. The tests and various measurements performed in this study demonstrate how the phantom can offer a reliable way of testing and investigating tumor volume segmentation in PET.

Deep-Inspiration Breath-Hold Intensity Modulated Radiation Therapy to the Mediastinum for Lymphoma Patients: Setup Uncertainties and Margins.

PURPOSE: Patient setup for treating large target volumes can be challenging. In the present study, we measured the local uncertainties in the treatment of mediastinal lymphoma and investigated the need for region-specific planning target volume (PTV) margins. METHODS AND MATERIALS: The data from 30 patients who had undergone radiation therapy for mediastinal lymphoma were retrospectively analyzed. A computed tomography (CT)-on-rails (CTOR) system in the treatment room was used for daily image guidance. The total PTV was split into 6 regions: neck, supraclavicular fossa, axilla, mediastinum, upper heart, and lower heart. The total PTV and the 6 local regions were separately aligned to the planning CT scan using automatic rigid registration. The residual local errors using 3 setup strategies were investigated: no image guidance, CTOR setup to total PTV, and simulated cone beam CT setup to the mediastinum. Errors were recorded in the anteroposterior, superoinferior, and right-left directions separately. Using the residual error calculations, the margins required to cover 95% of the clinical target volume for 90% of the patients was estimated. RESULTS: For each patient, 12 to 21 days of daily CTOR data were available for analysis. The residual errors for the total PTV and mediastinum setups were both smaller than those with no image guidance. The lower heart region had more uncertainty with all 3 setup strategies. Margin analysis revealed that the magnitude of the margin is dependent on the imaging strategy, direction, and local region inside the PTV. Margins >7 mm are necessary to account for uncertainty in the neck, lower heart, and axilla regions even under daily CT guidance. CONCLUSIONS: Setup uncertainties in the mediastinum are not uniform and are dependent on target location and imaging strategy. However, with the appropriate margin, we can target regions that might not be visualized with the available on-board imager system.

Evaluation of diffusion weighted imaging for tumor delineation in head-and-neck radiotherapy by comparison with automatically segmented 18F-fluorodeoxyglucose positron emission tomography.

BACKGROUND AND PURPOSE: Diffusion weighted (DW) MRI may facilitate target volume delineation for head-and-neck (HN) radiation treatment planning. In this study we assessed the use of a dedicated, geometrically accurate, DW-MRI sequence for target volume delineation. The delineations were compared with semi-automatic segmentations on 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET) images and evaluated for interobserver variation. METHODS AND MATERIALS: Fifteen HN cancer patients underwent both DW-MRI and FDG-PET for RT treatment planning. Target delineation on DW-MRI was performed by three observers, while for PET a semi-automatic segmentation was performed using a Gaussian mixture model. For interobserver variation and intermodality variation, volumes, overlap metrics and Hausdorff distances were calculated from the delineations. RESULTS: The median volumes delineated by the three observers on DW-MRI were 10.8, 10.5 and 9.0 cm3 respectively, and was larger than the median PET volume (8.0 cm3). The median conformity index of DW-MRI for interobserver variation was 0.73 (range 0.38-0.80). Compared to PET, the delineations on DW-MRI by the three observers showed a median dice similarity coefficient of 0.71, 0.69 and 0.72 respectively. The mean Hausdorff distance was small with median (range) distances between PET and DW-MRI of 2.3 (1.5-6.8), 2.5 (1.6-6.9) and 2.0 (1.35-7.6) mm respectively. Over all patients, the median 95th percentile distances were 6.0 (3.0-13.4), 6.6 (4.0-24.0) and 5.3 (3.4-26.0) mm. CONCLUSION: Using a dedicated DW-MRI sequence, target volumes could be defined with good interobserver agreement and a good overlap with PET. Target volume delineation using DW-MRI is promising in head-and-neck radiotherapy, combined with other modalities, it can lead to more precise target volume delineation.