Locoregional control of non-small cell lung cancer in relation to automated early assessment of tumor regression on cone beam computed tomography

Early radiologic and metabolic tumour response assessment during combined chemo-radiotherapy for locally advanced NSCLC

Background

The role of early treatment response for patients with locally advanced non-small cell lung cancer (LA-NSCLC) treated with concurrent chemo-radiotherapy (cCRT) is unclear. The study aims to investigate the predictive value of response to induction chemotherapy (iCX) and the correlation with pattern of failure (PoF).

Materials and methods

Patients with LA-NSCLC treated with cCRT were included for analyses (n = 276). Target delineations were registered from radiotherapy planning PET/CT to diagnostic PET/CT, in between which patients received iCX. Volume, sphericity, and SUVpeak were extracted from each scan. First site of failure was categorised as loco-regional (LR), distant (DM), or simultaneous LR+M (LR+M). Fine and Gray models for PoF were performed: a baseline model (including performance status (PS), stage, and histology), an image model for squamous cell carcinoma (SCC), and an image model for non-SCC. Parameters included PS, volume (VOL) of tumour, VOL of lymph nodes, ΔVOL, sphericity, SUVpeak, ΔSUVpeak, and oligometastatic disease.

Results

Median follow-up was 7.6 years. SCC had higher sub-distribution hazard ratio (sHR) for LRF (sHR = 2.771 [1.577:4.87], p < 0.01) and decreased sHR for DM (sHR = 0.247 [0.125:0.485], p  <  0.01). For both image models, high diagnostic SUVpeak increased risk of LRF (sHR = 1.059 [1.05:1.106], p < 0.01 for SCC, sHR = 1.12 [1.03:1.21], p < 0.01 for non-SCC). Patients with SCC and less decrease in VOL had higher sHR for DM (sHR = 1.025[1.001:1.048] pr. % increase, p = 0.038).

Conclusion

Poor response in disease volume was correlated with higher sHR of DM for SCC, no other clear correlation of response and PoF was observed. Histology significantly correlated with PoF with SCC prone to LRF and non-SCC prone to DM as first site of failure. High SUVpeak at diagnosis increased the risk of LRF for both histologies.

Prediction of SBRT response in liver cancer by combining original and delta cone-beam CT radiomics: a pilot study

This study aims to explore the feasibility of utilizing a combination of original and delta cone-beam CT (CBCT) radiomics for predicting treatment response in liver tumors undergoing stereotactic body radiation therapy (SBRT). A total of 49 patients are included in this study, with 36 receiving 5-fraction SBRT, 3 receiving 4-fraction SBRT, and 10 receiving 3-fraction SBRT. The CBCT and planning CT images from liver cancer patients who underwent SBRT are collected to extract overall 547 radiomics features. The CBCT features which are reproducible and interchangeable with pCT are selected for modeling analysis. The delta features between fractions are calculated to depict tumor change. The patients with 4-fraction SBRT are only used for screening robust features. In patients receiving 5-fraction SBRT, the predictive ability of both original and delta CBCT features for two-level treatment response (local efficacy vs. local non-efficacy; complete response (CR) vs. partial response (PR)) is assessed by utilizing multivariable logistic regression with leave-one-out cross-validation. Additionally, univariate analysis is conducted to validate the capability of CBCT features in identifying local efficacy in patients receiving 3-fraction SBRT. In patients receiving 5-fraction SBRT, the combined models incorporating original and delta CBCT radiomics features demonstrate higher area under the curve (AUC) values compared to models using either original or delta features alone for both classification tasks. The AUC values for predicting local efficacy vs. local non-efficacy are 0.58 for original features, 0.82 for delta features, and 0.90 for combined features. For distinguishing PR from CR, the respective AUC values for original, delta and combined features are 0.79, 0.80, and 0.89. In patients receiving 3-fraction SBRT, eight valuable CBCT radiomics features are identified for predicting local efficacy. The combination of original and delta radiomics derived from fractionated CBCT images in liver cancer patients undergoing SBRT shows promise in providing comprehensive information for predicting treatment response.

A Novel Method for Evaluating Early Tumor Response Based on Daily CBCT Images for Lung SBRT

BACKGROUND

We aimed to develop a new tumor response assessment method for lung SBRT.

METHODS

In total, 132 lung cancer patients with 134 tumors who received SBRT treatment with daily CBCT were included in this study. The information about tumor size (area), contrast (contrast-to-noise ratio (CNR)), and density/attenuation (μ) was derived from the CBCT images for the first and the last fractions. The ratios of tumor area, CNR, and μ (RA, RCNR, Rμ) between the last and first fractions were calculated for comparison. The product of the three rations was defined as a new parameter (R) for assessment. Tumor response was independently assessed by a radiologist based on a comprehensive analysis of the CBCT images.

RESULTS

R ranged from 0.27 to 1.67 with a mean value of 0.95. Based on the radiologic assessment results, a receiver operation characteristic (ROC) curve with the area under the curve (AUC) of 95% was obtained and the optimal cutoff value (RC) was determined as 1.1. The results based on RC achieved a 94% accuracy, 94% specificity, and 90% sensitivity.

CONCLUSION

The results show that R was correlated with early tumor response to lung SBRT and that using R for evaluating tumor response to SBRT would be viable and efficient.

Identification of reproducible radiomic features from on-board volumetric images: A multi-institutional phantom study

BACKGROUND

Radiomics analysis using on-board volumetric images has attracted research attention as a method for predicting prognosis during treatment; however, the lack of standardization is still one of the main concerns.

PURPOSE

This study investigated the factors that influence the reproducibility of radiomic features extracted from on-board volumetric images using an anthropomorphic radiomics phantom. Furthermore, a phantom experiment was conducted with different treatment machines from multiple institutions as external validation to identify reproducible radiomic features.

METHODS

The phantom was designed to be 35 × 20 × 20 cm with eight types of heterogeneous spheres (⌀ = 1, 2, and 3 cm). On-board volumetric images were acquired using 15 treatment machines from eight institutions. Of these, kilovoltage cone-beam computed tomography (kV-CBCT) image data acquired from four treatment machines at one institution were used as an internal evaluation dataset to explore the reproducibility of radiomic features. The remaining image data, including kV-CBCT, megavoltage-CBCT (MV-CBCT), and megavoltage computed tomography (MV-CT) provided by seven different institutions (11 treatment machines), were used as an external validation dataset. A total of 1,302 radiomic features, including 18 first-order, 75 texture, 465 (i.e., 93 × 5) Laplacian of Gaussian (LoG) filter-based, and 744 (i.e., 93 × 8) wavelet filter-based features, were extracted within the spheres. The intraclass correlation coefficient (ICC) was calculated to explore feature repeatability and reproducibility using an internal evaluation dataset. Subsequently, the coefficient of variation (COV) was calculated to validate the feature variability of external institutions. An absolute ICC exceeding 0.85 or COV under 5% was considered indicative of a highly reproducible feature.

RESULTS

For internal evaluation, ICC analysis showed that the median percentage of radiomic features with high repeatability was 95.2%. The ICC analysis indicated that the median percentages of highly reproducible features for inter-tube current, reconstruction algorithm, and treatment machine were decreased by 20.8%, 29.2%, and 33.3%, respectively. For external validation, the COV analysis showed that the median percentage of reproducible features was 31.5%. A total of 16 features, including nine LoG filter-based and seven wavelet filter-based features, were indicated as highly reproducible features. The gray-level run-length matrix (GLRLM) was classified as containing the most frequent features (N = 8), followed by the gray-level dependence matrix (N = 7) and gray-level co-occurrence matrix (N = 1) features.

CONCLUSIONS

We developed the standard phantom for radiomics analysis of kV-CBCT, MV-CBCT, and MV-CT images. With this phantom, we revealed that the differences in the treatment machine and image reconstruction algorithm reduce the reproducibility of radiomic features from on-board volumetric images. Specifically, the most reproducible features for external validation were LoG or wavelet filter-based GLRLM features. However, the acceptability of the identified features should be examined in advance at each institution before applying the findings to prognosis prediction.

Preliminary exploration of response the course of radiotherapy for stage III non-small cell lung cancer based on longitudinal CT radiomics features

Purpose

Explore the longitudinal CT-based radiomics to demonstrate the changing trend of radiotherapy response and to determine at which point after the onset of treatment radiomics exhibit the greatest change for stage III NSCLC patients.

Methods and materials

Ten stage III NSCLC patients in line with inclusion criteria were enrolled retrospectively, each of whom received radiotherapy or concurrent chemo-radiotherapy and performed eight series of follow-up CT imaging. Longitudinal radiomics were extracted on region of interest from the eight registered images, then two steps were conducted to select significant features as indicators of tumor change: 1) stable features were selected by Kendall rank correlation; 2) texture feature types with a steadily changing trend were retained and intensity features with stable change trends were selected to represent the large number of them. Next, the trend and rate of tumor change were analyzed using the Delta method and Curve-fitting method. Finally, the statistics in the distribution of stable features in patients were calculated.

Results

675 stable features were selected from a total number of 1371 radiomics features, then 12 texture features types were retained and three intensity features were chosen to represent their own category. Among the final selected feature types, it was found that the two time points were weeks 1 and 3 with the higher rate of change. One patient had very few stable tumor features out of a total of 101 features, and the rate of change of features of another patient was conspicuously higher than the average level with number of 301 features.

Conclusion

The longitudinal CT radiomics could demonstrate the change trend of tumor and at which point exhibit the greatest change during radiotherapy, and potentially be used for treatment decisions concerning adaptive radiotherapy.

Predicting response to radiotherapy of head and neck squamous cell carcinoma using radiomics from cone-beam CT images

INTRODUCTION

Radiotherapy (RT) for head and neck cancer is now guided by cone-beam computed tomography (CBCT). We aim to identify a CBCT radiomic signature predictive of progression to RT.

MATERIAL AND METHODS

A cohort of 93 patients was split into training (n = 60) and testing (n = 33) sets. A total of 88 features were extracted from the gross tumor volume (GTV) on each CBCT. Receiver operating characteristic (ROC) curves were used to determine the power of each feature at each week of treatment to predict progression to radio(chemo)therapy. Only features with AUC > 0.65 at each week were pre-selected. Absolute differences were calculated between features from each weekly CBCT and baseline CBCT1 images. The smallest detectable change (C = 1.96 × SD, SD being the standard deviation of differences between feature values calculated on CBCT1 and CBCTn) with its confidence interval (95% confidence interval [CI]) was determined for each feature. The features for which the change was larger than C for at least 5% of patients were then selected. A radiomics-based model was built at the time-point that showed the highest AUC and compared with models relying on clinical variables.

RESULTS

Seven features had an AUC > 0.65 at each week, and six exhibited a change larger than the predefined CI 95%. After exclusion of inter-correlated features, only one parameter remains, Coarseness. Among clinical variable, only hemoglobin value was significant. AUC for predicting the treatment response were 0.78 (p = .006), 0.85 (p < .001), and 0.99 (p < .001) for clinical, CBCT4-radiomics (Coarseness) and clinical + radiomics based models respectively. The mean AUC of this last model on a 5-fold cross-validation was 0.80 (±0.09). On the testing cohort, the best prediction was given by the combined model (balanced accuracy [BAcc] 0.67 , p < .001).

CONCLUSIONS

We described a feature selection methodology for delta-radiomics that is able to select reproducible features which are informative due to their change during treatment. A selected delta radiomics feature may improve clinical-based prediction models.

Identification of modes of tumor regression in non-small cell lung cancer patients during radiotherapy

PURPOSE

Observed gross tumor volume (GTV) shrinkage during radiotherapy (RT) raises the question of whether to adapt treatment to changes observed on the acquired images. In the literature, two modes of tumor regression have been described: elastic and non-elastic. These modes of tumor regression will affect the safety of treatment adaptation. This study applies a novel approach, using routine cone-beam computed tomography (CBCT) and deformable image registration to automatically distinguish between elastic and non-elastic tumor regression.

METHODS

In this retrospective study, 150 locally advanced non-small cell lung cancer patients treated with 55 Gray of radiotherapy were included. First, the two modes of tumor regression were simulated. For each mode of tumor regression, one timepoint was simulated. Based on the results of simulated data, the approach used for analysis in real patients was developed. CBCTs were non-rigidly registered to the baseline CBCT using a cubic B-spline algorithm, NiftyReg. Next, the Jacobian determinants were computed from the deformation vector fields. To capture local volume changes, 10 Jacobian values were sampled perpendicular to the surface of the GTV, across the lung-tumor boundary. From the simulated data, we can distinguish elastic from non-elastic tumor regression by comparing the Jacobian values samples between 5 and 12.5 mm inside and 5 and 12.5 mm outside the planning GTV. Finally, morphometric results were compared between tumors of different histologies.

RESULTS

Most patients (92.3%) in our cohort showed stable disease in the first week of treatment and non-elastic shrinkage in the later weeks of treatment. At week 2, 125 patients (88%) showed stable disease, three patients (2.1%) disease progression, and 11 patients (8%) regression. By treatment completion, 91 patients (64%) had stable disease, one patient (0.7%) progression and 46 patients (32%) regression. A slight difference in the mode of tumor change was observed between tumors of different histologies.

CONCLUSION

Our novel approach shows that it may be possible to automatically quantify and identify global changes in lung cancer patients during RT, using routine CBCT images. Our results show that different regions of the tumor change in different ways. Therefore, careful consideration should be taken when adapting RT.

Liquid Biopsies for Molecular Biology-Based Radiotherapy

Molecular alterations drive cancer initiation and evolution during development and in response to therapy. Radiotherapy is one of the most commonly employed cancer treatment modalities, but radiobiologic approaches for personalizing therapy based on tumor biology and individual risks remain to be defined. In recent years, analysis of circulating nucleic acids has emerged as a non-invasive approach to leverage tumor molecular abnormalities as biomarkers of prognosis and treatment response. Here, we evaluate the roles of circulating tumor DNA and related analyses as powerful tools for precision radiotherapy. We highlight emerging work advancing liquid biopsies beyond biomarker studies into translational research investigating tumor clonal evolution and acquired resistance.

Analysing the possibility of utilizing CBCT radiomics as an independent modality: a phantom study

Aim:

To verify if computed tomography (CT) radiomics were reproducible by cone beam CT (CBCT) radiomics by using Catphan^® 504.

Materials and Methods:

Catphan^® 504 was imaged using the default IGRT OBI CBCT imaging protocols and CT scanner. Seven known density image regions of the phantom were segmented and image feature was extracted by Imaging Biomarker Explorer (IBEX) software. The 49 selected features from four feature categories were analyzed by considering each region of interest (ROI) segment as individual image set. Correlation was studies using interclass correlation coefficient (ICC) and Pearson’s correlation coefficient.

Results:

The ICC of the three feature categories, namely intensity, GLCM, and GLRLM was significant (p-value<0.05) in comparison with CT, while the ICC of the fourth feature category, NID, was no significant. The average absolute Pearson’s correlation coefficient from the features of the images was as follows: CT: r=0.679±0.257, CBC_Thead: r=0.707±0.231, CBCT_thorax: r=0.643±0.260, and CBCT_pelvis: r=0.594±0.276.

Conclusion:

It seems that the various densities of Catphan^® 504 ROI image segments of the CT radiomics are reproducible with CBCT radiomics and CBCT radiomics can be used as an independent modality.

The dynamics and prognostic value of FDG PET-metrics in weekly monitoring of (chemo)radiotherapy for NSCLC

PURPOSE

To test if the relative change in FDG-PET SUV_max over the course of treatment was associated with disease progression and overall survival. Additionally, the prognostic values of other first-order PET-metric changes were investigated.

METHODS

The study included 38 patients with stage II-III NSCLC, who underwent concurrent chemoradiotherapy. Patients received two pre-treatment FDG-PET scans and four during-treatment scans at weekly intervals. SUV_max was normalized to the start of treatment and analyzed using linear regression. Linear regression coefficients of other first order PET-metrics were grouped according to dissimilarity. Associations to patient outcome were analyzed using Cox hazard ratio.

RESULTS

Twenty-eight patients satisfied the criteria for analysis. All PET-metrics demonstrated a strong linear correlation with time during treatment [median R-range: -0.87: -0.97]. No strong associations (p > 0.10) were found for the relative slope of SUV_max to patient outcomes. Other first-order metrics did correlate with outcome but the single imaging time-point maximizing the association of PET response with outcome varied per PET metric and outcome parameter.

CONCLUSION

All investigated FDG PET metrics linearly decreased during treatment. Relative change in SUV_max was not associated to patient outcome while several other first order PET-metrics were related to patient outcome. A single optimal imaging time-point could not be identified.

A method for predictive modeling of tumor regression for lung adaptive radiotherapy

PURPOSE

The purpose of this work is to create a decision support methodology to predict when patients undergoing radiotherapy treatment for locally advanced lung cancer would potentially benefit from adaptive radiotherapy. The proposed methodology seeks to eliminate the manual subjective review by developing an automated statistical learning model to predict when tumor regression would trigger implementation of adaptive radiotherapy based on quantified anatomic changes observed in individual patients on-treatment cone beam computed tomographies (CTs). This proposed process seeks to improve the efficacy and efficiency of both the existing manual and automated adaptive review processes for locally advanced stage III lung cancer.

METHODS

A predictive algorithm was developed as a decision support tool to determine the potential utility of mid-treatment adaptive radiotherapy based on anatomic changes observed on 1158 daily CBCT images across 43 patients. The anatomic changes on each axial slice within specified regions-of-interest were quantified into a single value utilizing imaging similarity criteria comparing the daily CBCT to the initial simulation CT. The range of the quantified metrics for each fraction across all axial slices are reduced to specified quantiles, which are used as the predictive input to train a logistic regression algorithm. A "ground-truth" of the need for adaptive radiotherapy based on tumor regression was evaluated systematically on each of the daily CBCTs and used as the classifier in the logistic regression algorithm. Accuracy of the predictive model was assessed utilizing both a tenfold cross validation and an independent validation dataset, with the sensitivity, specificity, and fractional accuracy compared to the ground-truth.

RESULTS

The sensitivity and specificity for the individual daily fractions ranged from 87.9%-94.3% and 91.9%-98.6% for a probability threshold of 0.2-0.5, respectively. The corresponding average treatment fraction difference between the model predictions and assessed ART "ground-truth" ranged from -2.25 to -0.07 fractions, with the model predictions consistently predicting the potential need for ART earlier in the treatment course. By initially utilizing a lower probability threshold, the higher sensitivity minimizes the chance of false negative by alerting the clinician to review a higher number of questionable cases.

CONCLUSIONS

The proposed methodology accurately predicted the first fraction at which individual patients may benefit from ART based on quantified anatomic changes observed in the on-treatment volumetric imaging. The generalizability of the proposed method has potential to expand to additional modes of adaptive radiotherapy for lung cancer patients with observed underlying anatomic changes.

A Pilot Study Examining the Prognostic Utility of Tumor Shrinkage on Cone-Beam Computed Tomography (CBCT) for Stage III Locally Advanced Non-Small Cell Lung Cancer Patients Treated with Definitive Chemoradiation

There has been growing interest in utilizing information from cone-beam computed tomography (CBCT) to help guide both treatment delivery and prognosis. In this assessment of locally advanced unresectable stage III non-small cell lung cancer (NSCLC) treated with definitive chemoradiation, we aimed to determine the survival advantage associated with using CBCT to measure tumor regression. Patient, tumor, and treatment characteristics were collected. The serial tumor shrinkage for each patient was determined from tumor volume contours on weekly CBCTs. Survival analysis was performed using the Kaplan–Meier technique and a Cox proportional hazards model. At least two-thirds of patients had a tumor volume reduction of at least 5% after each week of chemoradiation. A weekly reduction in tumor volume of 5% or greater seen on the CBCT images during radiation therapy was significantly associated with improved overall survival, which remained significant when adjusted for age, histology, grade, and T- and N-stages (p = 0.0036). Additionally, the presence of N3 disease was associated with a five-fold increased risk of recurrence (p = 0.0006) and a nearly three-fold increased risk of death (p = 0.053) compared with N0–N2 disease. Tumor volume shrinkage observed in the CBCT images during definitive chemoradiation holds promise as a prognostic indicator of stage III NSCLC, especially given its affordability, availability, and applicability. Further evaluation in a prospective fashion is warranted to validate the tumor volume shrinkage and its clinical utility.

The Predictive Value of Tumor Volume and Its Change on Short-Term Outcome for Esophageal Squamous Cell Carcinoma Treated With Radiotherapy or Chemoradiotherapy

Objectives

To investigate the tumor volume and its change on short-term outcome in esophageal squamous cell carcinoma (ESCC) patients who underwent definitive radiotherapy or chemoradiotherapy.

Methods and Materials

All data were retrospectively collected from 418 ESCC patients who received radiotherapy or chemoradiotherapy at our institution between 2015 and 2019. Short-term outcome using the treatment response evaluation was assessed according to the RECIST 1.1. The tumor volume change rate (TVCR) was defined as follows: TVCR = {1 - [gross tumor volume (GTV) at shrinking irradiation field planning)]/(GTV at the initial treatment planning)} ×100%. Chi square test was used to compare the clinic characteristics in different TVCR groups, and the difference between initial GTV (GTVi) and shrinking GTV (GTVs) was compared using Wilcoxon's sign rank test. Logistic regression analysis and Spearman correlation was performed.

Results

There was a significant decrease in GTVi compared to GTVs (P < 0.001). In univariate analysis, age, cT-stage, TNM stage, treatment modality, GTVi, and TVCR were associated with short-term outcome (all P < 0.05). In multivariate analysis, gender and TVCR were statistically significant (P = 0.010, <0.001) with short-term outcome, and the combined predictive value of gender and TVCR exceeded that of TVCR (AUC, 0.876 vs 0.855).

Conclusions

TVCR could serve to forecast short-term outcome of radiotherapy or chemoradiotherapy in ESCC. It was of great significance to guide the individualized treatment of ESCC.

Identification of patterns of tumour change measured on CBCT images in NSCLC patients during radiotherapy

In this study, we propose a novel approach to investigate changes in the visible tumour and surrounding tissues with the aim of identifying patterns of tumour change during radiotherapy (RT) without segmentation on the follow-up images. On-treatment cone-beam computed tomography (CBCT) images of 240 non-small cell lung cancer (NSCLC) patients who received 55 Gy of RT were included. CBCTs were automatically aligned onto planning computed tomography (planning CT) scan using a two-step rigid registration process. To explore density changes across the lung-tumour boundary, eight shells confined to the shape of the gross tumour volume (GTV) were created. The shells extended 6 mm inside and outside of the GTV border, and each shell is 1.5 mm thick. After applying intensity correction on CBCTs, the mean intensity was extracted from each shell across all CBCTs. Thereafter, linear fits were created, indicating density change over time in each shell during treatment. The slopes of all eight shells were clustered to explore patterns in the slopes that show how tumours change. Seven clusters were obtained, 97% of the patients were clustered into three groups. After visual inspection, we found that these clusters represented patients with little or no density change, progression and regression. For the three groups, the survival curves were not significantly different between the groups, p-value = 0.51. However, the results show that definite patterns of tumour change exist, suggesting that it may be possible to identify patterns of tumour changes from on-treatment CBCT images.

Geometric and Dosimetric Evaluation of a Commercially Available Auto-segmentation Tool for Gross Tumour Volume Delineation in Locally Advanced Non-small Cell Lung Cancer: a Feasibility Study

AIMS

To quantify the reliability of a commercially available auto-segmentation tool in locally advanced non-small cell lung cancer using serial four-dimensional computed tomography (4DCT) scans during conventionally fractionated radiotherapy.

MATERIALS AND METHODS

Eight patients with serial 4DCT scans (n = 44) acquired over the course of radiotherapy were assessed. Each 4DCT had a physician-defined primary tumour manual contour (MC). An auto-contour (AC) and a user-adjusted auto-contour (UA-AC) were created for each scan. Geometric agreement of the AC and the UA-AC to the MC was assessed using the dice similarity coefficient (DSC), the centre of mass (COM) shift from the MC and the structure volume difference from the MC. Bland Altman analysis was carried out to assess agreement between contouring methods. Dosimetric reliability was assessed by comparison of planning target volume dose coverage on the MC and UA-AC. The time trend analysis of the geometric accuracy measures from the initial planning scan through to the final scan for each patient was evaluated using a Wilcoxon signed ranks test to assess the reliability of the UA-AC over the duration of radiotherapy.

RESULTS

User adjustment significantly improved all geometric comparison metrics over the AC alone. Improved agreement was observed in smaller tumours not abutting normal soft tissue and median values for geometric comparisons to the MC for DSC, tumour volume difference and COM offset were 0.80 (range 0.49-0.89), 0.8 cm³ (range 0.0-5.9 cm³) and 0.16 cm (range 0.09-0.69 cm), respectively. There were no significant differences in dose metrics measured from the MC and the UA-AC after Bonferroni correction. Variation in geometric agreement between the MC and the UA-AC were observed over the course of radiotherapy with both DSC (P = 0.035) and COM shift from the MC (ns) worsening. The median tumour volume difference from the MC improved at the later time point.

CONCLUSIONS

These findings suggest that the UA-AC can produce geometrically and dosimetrically acceptable contours for appropriately selected patients with non-small cell lung cancer. Larger studies are required to confirm the findings.

Investigation of inter-fraction target motion variations in the context of pencil beam scanned proton therapy in non-small cell lung cancer patients

Purpose

For locally advanced‐stage non‐small cell lung cancer (NSCLC), inter‐fraction target motion variations during the whole time span of a fractionated treatment course are assessed in a large and representative patient cohort. The primary objective is to develop a suitable motion monitoring strategy for pencil beam scanning proton therapy (PBS‐PT) treatments of NSCLC patients during free breathing.

Methods

Weekly 4D computed tomography (4DCT; 41 patients) and daily 4D cone beam computed tomography (4DCBCT; 10 of 41 patients) scans were analyzed for a fully fractionated treatment course. Gross tumor volumes (GTVs) were contoured and the 3D displacement vectors of the centroid positions were compared for all scans. Furthermore, motion amplitude variations in different lung segments were statistically analyzed. The dosimetric impact of target motion variations and target motion assessment was investigated in exemplary patient cases.

Results

The median observed centroid motion was 3.4 mm (range: 0.2–12.4 mm) with an average variation of 2.2 mm (range: 0.1–8.8 mm). Ten of 32 patients (31.3%) with an initial motion <5 mm increased beyond a 5‐mm motion amplitude during the treatment course. Motion observed in the 4DCBCT scans deviated on average 1.5 mm (range: 0.0–6.0 mm) from the motion observed in the 4DCTs. Larger motion variations for one example patient compromised treatment plan robustness while no dosimetric influence was seen due to motion assessment biases in another example case.

Conclusions

Target motion variations were investigated during the course of radiotherapy for NSCLC patients. Patients with initial GTV motion amplitudes of < 2 mm can be assumed to be stable in motion during the treatment course. For treatments of NSCLC patients who exhibit motion amplitudes of > 2 mm, 4DCBCT should be considered for motion monitoring due to substantial motion variations observed.

Cone-beam CT radiomics features might improve the prediction of lung toxicity after SBRT in stage I NSCLC patients

Background

Stereotactic body radiotherapy (SBRT) is the standard care for inoperable early stage non‐small cell lung cancer (NSCLC). The purpose of our study was to investigate whether a prediction model based on cone‐beam CT (CBCT) plus pretreatment CT radiomics features could improve the prediction of tumor control and lung toxicity after SBRT in comparison to a model based on pretreatment CT radiomics features alone.

Methods

A total of 34 cases of stage I NSCLC patients who received SBRT were included in the study. The pretreatment planning CT and serial CBCT radiomics features were analyzed using the imaging biomarker explorer (IBEX) software platform. Multivariate logistic regression was conducted for the association between progression‐free survival (PFS), lung toxicity and features. The predictive capabilities of the models based on CBCT and CT features were compared using receiver operating characteristic (ROC) curves.

Results

Five CBCT features and two planning CT features were correlated with disease progression. Six CBCT features and two planning CT features were related to lung injury. The ROC curves indicated that the model based on the CBCT plus planning CT features might be better than the model based on the planning CT features in predicting lung injury. The other ROC curves indicated that the model based on the planning CT features was similar to the model based on the CBCT plus planning CT features in predicting disease progression.

Conclusions

Both pretreatment CT and CBCT radiomics features could predict disease progression and lung injury. A model with CBCT plus pretreatment CT radiomics features might improve the prediction of lung toxicity in comparison with a model with pretreatment CT features alone.

Key points

Significant findings of the study: A model with cone‐beam CT radiomics features plus pre‐treatment CT radiomics features might improve the prediction of lung toxicity after SBRT in stage I NSCLC patients.

What this study adds: In the prediction of PFS and lung toxicity in early‐stage NSCLC patients treated with SBRT, CBCT radiomics could be another effective method.

Diffusion MRI outlined viable tumour volume beats GTV in intra-treatment stratification of outcome

BACKGROUND AND PURPOSE

In radiotherapy, treatment response is generally evaluated many weeks after end of the treatment course. If the treatment outcome could be predicted during radiotherapy better tumour control could be achieved through timely adaptation of the treatment strategy. In this study intra-treatment change based on the diffusion MRI outlined viable tumour volume (VTV) was assessed and compared to the standard GTV to study their outcome prediction capacity.

MATERIALS AND METHODS

Thirty-eight brain metastases from twenty-one cancer patients were analysed in this prospective trial. Diffusion and structural MRI was acquired on a 1 T machine before, during, and at follow-up 2-3 months after radiotherapy. The VTV was defined as a region with high cellularity using high b-value diffusion MRI scans. Further, the diffusivity of the VTV was derived as the apparent diffusion coefficient (ADC). Treatment outcome was determined using RECIST defined bounds in the T1W MRI follow-up scan. Longitudinal statistical analysis was performed using a linear mixed effect model.

RESULTS

The GTV declined in both responding and non-responding (significantly) tumours with inseparable rates during radiotherapy. The VTV volume fraction reduced significantly in the responding tumours only. The ADC of the VTV increased significantly in responding metastases whereas it decreased in non-responding metastases. Furthermore, no association between baseline tumour size or primary disease and outcome was observed.

CONCLUSION

GTV size change during radiotherapy is not a reliable predictor of outcome in brain metastases. On the other hand, change in the volume fraction of VTV and diffusivity of VTV shows ability to stratify treatment outcome.

Tumor regression during radiotherapy for non-small cell lung cancer patients using cone-beam computed tomography images

PURPOSE

Previous literature has reported contradicting results regarding the relationship between tumor volume changes during radiotherapy treatment for non-small cell lung cancer (NSCLC) patients and locoregional recurrence-free rate or overall survival. The aim of this study is to validate the results from a previous study by using a different volume extraction procedure and evaluating an external validation dataset.

METHODS

For two datasets of 94 and 141 NSCLC patients, gross tumor volumes were determined manually to investigate the relationship between tumor volume regression and locoregional control using Kaplan-Meier curves. For both datasets, different subgroups of patients based on histology and chemotherapy regimens were also investigated. For the first dataset (n = 94), automatically determined tumor volumes were available from a previously published study to further compare their correlation with updated clinical data.

RESULTS

A total of 70 out of 94 patients were classified into the same group as in the previous publication, splitting the dataset based on median tumor regression calculated by the two volume extraction methods. Non-adenocarcinoma patients receiving concurrent chemotherapy with large tumor regression show reduced locoregional recurrence-free rates in both datasets (p < 0.05 in dataset 2). For dataset 2, the opposite behavior is observed for patients not receiving chemotherapy, which was significant for overall survival (p = 0.01) but non-significant for locoregional recurrence-free rate (p = 0.13).

CONCLUSION

The tumor regression pattern observed during radiotherapy is not only influenced by irradiation but depends largely on the delivered chemotherapy schedule, so it follows that the relationship between patient outcome and the degree of tumor regression is also largely determined by the chemotherapy schedule. This analysis shows that the relationship between tumor regression and outcome is complex, and indicates factors that could explain previously reported contradicting findings. This, in turn, will help guide future studies to fully understand the relationship between tumor regression and outcome.

Anatomical Adaptation-Early Clinical Evidence of Benefit and Future Needs in Lung Cancer

Definitive treatment of locally advanced non-small-cell lung cancer with radiation is challenging. During the course of treatment, anatomical changes such as tumor regression, tumor displacement/deformation, pleural effusion, and/or atelectasis can result in a deviation of the administered radiation dose from the intended prescribed treatment and thereby worsen local control and toxicity. Adaptive radiotherapy can help correct for these changes and can be generally categorized into 3 philosophical paradigms: (1) maintenance of prescribed dose to the initially defined target volume; (2) dose reduction to healthy organs while maintaining initial prescribed dose to a regressing tumor volume; or (3) dose escalation to a regressing tumor volume with isotoxicity to healthy organs. Numerous single institution studies have investigated these methods, and results from large prospective clinical trials will hopefully provide consensus on the method, utility, and efficacy of implementing adaptive radiation therapy (ART) in a clinical setting. Additional development into standardization and automation of the ART workflow, specifically in identifying when ART is warranted and in reducing the manual clinical effort needed to produce an adaptive plan, will be paramount to making ART feasible for the broader radiation therapy community.

Longitudinal radiomics of cone-beam CT images from non-small cell lung cancer patients: Evaluation of the added prognostic value for overall survival and locoregional recurrence

•

Longitudinal CBCT radiomics does not show added prognostic information for NSCLC.

•

A CT-radiomics model could be validated in an external validation dataset.

•

Dataset heterogeneity and small cohort sizes could cause poor validation performance.

Background and purpose

The prognostic value of radiomics for non-small cell lung cancer (NSCLC) patients has been investigated for images acquired prior to treatment, but no prognostic model has been developed that includes the change of radiomic features during treatment. Therefore, the aim of this study was to investigate the potential added prognostic value of a longitudinal radiomics approach using cone-beam computed tomography (CBCT) for NSCLC patients.

Materials and methods

This retrospective study includes a training dataset of 141 stage I–IV NSCLC patients and three external validation datasets of 94, 61 and 41 patients, all treated with curative intended (chemo)radiotherapy. The change of radiomic features extracted from CBCT images was summarized as the slope of a linear regression. The CBCT slope-features and CT-extracted features were used as input for a Cox proportional hazards model. Moreover, prognostic performance of clinical parameters was investigated for overall survival and locoregional recurrence. Model performances were assessed using the Kaplan–Meier curves and c-index.

Results

The radiomics model contained only CT-derived features and reached a c-index of 0.63 for overall survival and could be validated on the first validation dataset. No model for locoregional recurrence could be developed that validated on the validation datasets. The clinical parameters model could not be validated for either overall survival or locoregional recurrence.

Conclusion

In this study we could not confirm our hypothesis that longitudinal CBCT-extracted radiomic features contribute to improved prognostic information. Moreover, performance of baseline radiomic features or clinical parameters was poor, probably affected by heterogeneity within and between datasets.

Incidence and evolution of imaging changes on cone-beam CT during and after radical radiotherapy for non-small cell lung cancer

BACKGROUND AND PURPOSE

Cone beam CT (CBCT) is used to improve accuracy of radical radiotherapy by adjusting treatment to the observed imaging changes. To ensure appropriate adjustment, image interpretation should precede any changes to treatment delivery. This study provides the methodology for image interpretation and the frequency and evolution of the changes in patients undergoing radical radiotherapy for localised and locally advanced non-small cell lung cancer (NSCLC).

PATIENTS AND METHODS

From December 2012 to December 2014, 250 patients with localised and locally advanced NSCLC had 2462 chest CBCT scans during the course of fractionated radical radiotherapy (RT) (3-5 daily CBCTs in the first week followed by at least weekly imaging, mean 9.5 per patient, range 1-21). All CBCT images were reviewed describing changes and their evolution using diagnostic imaging definitions and validated by an independent chest radiologist.

RESULTS

During radical RT for NSCLC 328 imaging changes were identified on CBCT in 180 (72%) patients; 104 (32%) had reduction and 41 (13%) increase in tumour size; 48 (15%) had changes in consolidations contiguous to the primary lesion, 26 (8%) non-contiguous consolidations, 43 (13%) changes in tumour cavitation, 36 (11%) pleural effusion and 30 (9%) changes in atelectasis. In 105 patients imaging changes were noted in continuity with the treated tumour of which only 41 (39%) represented tumour enlargement; others included new or enlarging adjacent consolidation (34%), and new or enlarging atelectasis (19%). The changes evolved during treatment.

CONCLUSION

Imaging changes on CBCT include real and apparent changes in tumour size and parenchymal changes which evolve during treatment. Correct image interpretation, particularly when occurring adjacent to the tumour, is essential prior to adjustment to treatment delivery.

Greater reduction in mid-treatment FDG-PET volume may be associated with worse survival in non-small cell lung cancer

BACKGROUND AND PURPOSE

This study tested the hypotheses that 1) changes in mid-treatment fluorodeoxyglucose (FDG)-positron emission tomography (PET) parameters are predictive of overall survival (OS) and 2) mid-treatment FDG-PET-adapted treatment has the potential to improve survival in patients with non-small cell lung cancer (NSCLC).

MATERIAL AND METHODS

Patients with stage I-III NSCLC requiring daily fractionated radiation were eligible. FDG-PET-CT scans were obtained prior to and mid-treatment with radiotherapy at 40-50 Gy. The normalized maximum standardized uptake value (NSUVmax), normalized mean SUV (NSUVmean), PET-metabolic tumor volume (MTV), total lesion glycolysis (TLG), and computed tomography-based gross tumor volume (CT-GTV) were consistently measured for all patients. The primary study endpoint was OS.

RESULTS

The study is comprised of 102 patients who received 3-dimensional conformal radiotherapy, among whom 30 patients who received mid-treatment PET-adapted dose escalation radiotherapy. All PET-CT parameters decreased significantly (P < 0.001) mid-treatment, with greater reductions in FDG-volumetric parameters compared to FDG-activity factors. Mid-treatment changes in MTV (P = 0.053) and TLG (P = 0.021) were associated with OS, while changes in NSUVmax, NSUVmean, and CT-GTV were not (all Ps>0.1). Patients receiving conventional radiation (60-70 Gy) with MTV reductions greater than the mean had a median survival of 14 months, compared to those with MTV reductions less than the mean who had a median survival of 22 months. By contrast, patients receiving mid-treatment PET-adapted radiation with MTV reductions greater than the mean had a median survival of 33 months, compared to those with MTV reductions less than the mean who had a median survival of 19 months. Overall, PET-adapted treatment resulted in a 19% better 5-year survival than conventional radiation.

CONCLUSION

Changes in mid-treatment PET-volumetric parameters were significantly associated with survival in NSCLC. A greater reduction in the mid-treatment MTV was associated with worse survival in patients treated with standard radiation, but with better survival in patients who received mid-treatment PET-adapted treatment.

[Image-guided radiotherapy in lung cancer]

Image-guided radiotherapy takes place at every step of the treatment in lung cancer, from treatment planning, with fusion imaging, to daily in-room repositioning. Managing tumoral and surrounding thoracic structures motion has been allowed since the routine use of 4D computed tomography (4DCT). The integration of respiratory motion has been made with "passive" techniques based on reconstruction images from 4DCT planning, or "active" techniques adapted to the patient's breathing. Daily repositioning is based on regular images, weekly or daily, low (kV) or high (MV) energy. MRI and functional imaging also play an important part in lung cancer radiation and open the way for adaptative radiotherapy.

CALIPER: A deformable image registration algorithm for large geometric changes during radiotherapy for locally advanced non-small cell lung cancer

PURPOSE

Locally advanced non-small cell lung cancer (NSCLC) patients may experience dramatic changes in anatomy during radiotherapy and could benefit from adaptive radiotherapy (ART). Deformable image registration (DIR) is necessary to accurately accumulate dose during plan adaptation, but current algorithms perform poorly in the presence of large geometric changes, namely atelectasis resolution. The goal of this work was to develop a DIR framework, named Consistent Anatomy in Lung Parametric imagE Registration (CALIPER), to handle large geometric changes in the thorax.

METHODS

Registrations were performed on pairs of baseline and mid-treatment CT datasets of NSCLC patients presenting with atelectasis at the start of treatment. Pairs were classified based on atelectasis volume change as either full, partial, or no resolution. The evaluated registration algorithms consisted of several combinations of a hybrid intensity- and feature-based similarity cost function to investigate the ability to simultaneously match healthy lung parenchyma and adjacent atelectasis. These components of the cost function included a mass-preserving intensity cost in the lung parenchyma, use of filters to enhance vascular structures in the lung parenchyma, manually delineated lung lobes as labels, and several intensity cost functions to model atelectasis change. Registration error was quantified with landmark-based target registration error and post-registration alignment of atelectatic lobes.

RESULTS

The registrations using both lobe labels and vasculature enhancement in addition to intensity of the CT images were found to have the highest accuracy. Of these registrations, the mean (SD) of mean landmark error across patients was 2.50 (1.16) mm, 2.80 (0.70) mm, and 2.04 (0.13) mm for no change, partial resolution, and full atelectasis resolution, respectively. The mean (SD) atelectatic lobe Dice similarity coefficient was 0.91 (0.08), 0.90 (0.08), and 0.89 (0.04), respectively, for the same groups. Registration accuracy was comparable to healthy lung registrations of current state-of-the-art algorithms reported in literature.

CONCLUSIONS

The CALIPER algorithm developed in this work achieves accurate image registration for challenging cases involving large geometric and topological changes in NSCLC patients, a requirement for enabling ART in this patient group.

Is heterogeneity in stage 3 non-small cell lung cancer obscuring the potential benefits of dose-escalated concurrent chemo-radiotherapy in clinical trials?

The current standard of care for the management of inoperable stage 3 non-small cell lung cancer (NSCLC) is concurrent chemoradiotherapy (cCRT) using radiotherapy dose-fractionation and chemotherapy regimens that were established 3 decades ago. In an attempt to improve the chances of long-term control from cCRT, dose-escalation of the radiotherapy dose was assessed in the RTOG 0617 randomised control study comparing the standard 60 Gy in 30 fractions with a high-dose arm receiving 74 Gy in 37 fractions. Following the publication of this trial the thoracic oncology community were surprised to learn that there was worse survival in the dose-escalated arm and that for now the standard of care must remain with the lower dose. In this article we review the RTOG 0617 paper with subsequent analyses and studies to explore why the use of dose-escalated cCRT in stage 3 NSCLC has not shown the benefits that were expected. The overarching theme of this opinion piece is how heterogeneity between stage 3 NSCLC cases in terms of patient, tumour, and clinical factors may obscure the potential benefits of dose-escalation by causing imbalances in the arms of studies such as RTOG 0617. We also examine recent advances in the staging, management, and technological delivery of radiotherapy in NSCLC and how these may be employed to optimise cCRT trials in the future and ensure that any potential benefits of dose-escalation can be detected.

The Value of CBCT-based Tumor Density and Volume Variations in Prediction of Early Response to Chemoradiation Therapy in Advanced NSCLC

The correlations between early responses and the variations in physical density and primary tumor volume (TV) according to cone-beam computed tomography (CBCT) during chemoradiotherapy for non-small cell lung cancer (NSCLC) patients were investigated. 54 patients with inoperable and locally advanced NSCLC were included in this study. The CT numbers (CTN) and TV were measured on each of the seven observation points. The changes in the mean CTN values and the variation ratios of TV during the treatment course were analysed and correlated with the clinical outcomes, as evaluated by the RECIST criteria. For patients who responded to treatment, the CTN and TV change ratio decreased by 28.44 ± 13.12 HU and 32.01% (range, 8.46-61.67%); these values were significantly higher than those in the non-responding patients, with 19.63 ± 8.67 HU and 23.20% (range, -15.57-38%) (p = 0.016, p = 0.048), respectively. The area under curve for the combination of CTN and TV was larger than either alone (AUC = 0.751, p = 0.002). The differences between response and non-response were most significant between Fraction 10 and Fraction 15 for CTN changes and between Fraction 5 and Fraction 10 for the TV regression ratio. The changes in CTN and TV obtained from CBCT images have the potential capability to predict an early response of NSCLC.

Feature selection methodology for longitudinal cone-beam CT radiomics

BACKGROUND

Cone-beam CT (CBCT) scans are typically acquired daily for positioning verification of non-small cell lung cancer (NSCLC) patients. Quantitative information, derived using radiomics, can potentially contribute to (early) treatment adaptation. The aims of this study were to (1) describe and investigate a methodology for feature selection of a longitudinal radiomics approach (2) investigate which time-point during treatment is potentially useful for early treatment response assessment.

MATERIAL AND METHODS

For 90 NSCLC patients CBCT scans of the first two fractions of treatment (considered as 'test-retest' scans) were analyzed, as well as weekly CBCT images. One hundred and sixteen radiomic features were extracted from the GTV of all scans and subsequently absolute and relative differences were calculated between weekly CBCT images and the CBCT of the first fraction. Test-retest scans were used to determine the smallest detectable change (C = 1.96 * SD) allowing for feature selection by choosing a minimum number of patients for which a feature should change more than 'C' to be considered as relevant. Analysis of which features change at which moment during treatment was used to investigate which time-point is potentially relevant to extract longitudinal radiomics information for early treatment response assessment.

RESULTS

A total of six absolute delta features changed for at least ten patients at week 2 of treatment and increased to 61 at week 3, 79 at week 4 and 85 at week 5. There was 93% overlap between features selected at week 3 and the other weeks.

CONCLUSIONS

This study describes a feature selection methodology for longitudinal radiomics that is able to select reproducible delta radiomics features that are informative due to their change during treatment, which can potentially be used for treatment decisions concerning adaptive radiotherapy. Nonetheless, the prognostic value of the selected delta radiomic features should be investigated in future studies.

Impact of beam angle choice on pencil beam scanning breath-hold proton therapy for lung lesions

INTRODUCTION

The breath-hold technique inter alia has been suggested to mitigate the detrimental effect of motion on pencil beam scanned (PBS) proton therapy dose distributions. The aim of this study was to evaluate the robustness of incident proton beam angles to day-to-day anatomical variations in breath-hold.

MATERIALS AND METHODS

Single field PBS plans at five degrees increments in the transversal plane were made and water-equivalent path lengths (WEPLs) were derived on the planning breath-hold CT (BHCT) for 30 patients diagnosed with locally-advanced non-small cell lung cancer (NSCLC), early stage NSCLC or lung metastasis. Our treatment planning system was subsequently used to recalculate the plans and derive WEPL on a BHCT scan acquired at the end of the treatment. Changes to the V_95%, D₉₅ and mean target dose were evaluated.

RESULTS

The difference in WEPL as a function of the beam angle was highly patient specific, with a median of 3.3 mm (range: 0.0-41.1 mm). Slightly larger WEPL differences were located around the lateral or lateral anterior/posterior beam angles. Linear models revealed that changes in dose were associated to the changes in WEPL and the tumor baseline shift (p < 0.05).

CONCLUSIONS

WEPL changes and tumor baseline shift can serve as reasonable surrogates for dosimetric uncertainty of the target coverage and are well-suited for routine evaluation of plan robustness. The two lateral beam angles are not recommended to use for PBS proton therapy of lung cancer patients treated in breath-hold, due to the poor robustness for several of the patients evaluated.

Statistical tools for the temporal analysis and classification of lung lesions

BACKGROUND AND OBJECTIVE

Lung cancer remains one of the most common cancers globally. Temporal evaluation is an important tool for analyzing the malignant behavior of lesions during treatment, or of indeterminate lesions that may be benign. This work proposes a methodology for the analysis, quantification, and visualization of small (local) and large (global) changes in lung lesions. In addition, we extract textural features for the classification of lesions as benign or malignant.

METHODS

We employ the statistical concept of uncertainty to associate each voxel of a lesion to a probability that changes occur in the lesion over time. We employ the Jensen divergence and hypothesis test locally to verify voxel-to-voxel changes, and globally to capture changes in lesion volumes.

RESULTS

For the local hypothesis test, we determine that the change in density varies by between 3.84 and 40.01% of the lesion volume in a public database of malignant lesions under treatment, and by between 5.76 and 35.43% in a private database of benign lung nodules. From the texture analysis of regions in which the density changes occur, we are able to discriminate lung lesions with an accuracy of 98.41%, which shows that these changes can indicate the true nature of the lesion.

CONCLUSION

In addition to the visual aspects of the density changes occurring in the lesions over time, we quantify these changes and analyze the entire set using volumetry. In the case of malignant lesions, large b-divergence values are associated with major changes in lesion volume. In addition, this occurs when the change in volume is small but is associated with significant changes in density, as indicated by the histogram divergence. For benign lesions, the methodology shows that even in cases where the change in volume is small, a change of density occurs. This proves that even in lesions that appear stable, a change in density occurs.

A longitudinal four-dimensional computed tomography and cone beam computed tomography dataset for image-guided radiation therapy research in lung cancer

PURPOSE

To describe in detail a dataset consisting of serial four-dimensional computed tomography (4DCT) and 4D cone beam CT (4DCBCT) images acquired during chemoradiotherapy of 20 locally advanced, nonsmall cell lung cancer patients we have collected at our institution and shared publicly with the research community.

ACQUISITION AND VALIDATION METHODS

As part of an NCI-sponsored research study 82 4DCT and 507 4DCBCT images were acquired in a population of 20 locally advanced nonsmall cell lung cancer patients undergoing radiation therapy. All subjects underwent concurrent radiochemotherapy to a total dose of 59.4-70.2 Gy using daily 1.8 or 2 Gy fractions. Audio-visual biofeedback was used to minimize breathing irregularity during all fractions, including acquisition of all 4DCT and 4DCBCT acquisitions in all subjects. Target, organs at risk, and implanted fiducial markers were delineated by a physician in the 4DCT images. Image coordinate system origins between 4DCT and 4DCBCT were manipulated in such a way that the images can be used to simulate initial patient setup in the treatment position. 4DCT images were acquired on a 16-slice helical CT simulator with 10 breathing phases and 3 mm slice thickness during simulation. In 13 of the 20 subjects, 4DCTs were also acquired on the same scanner weekly during therapy. Every day, 4DCBCT images were acquired on a commercial onboard CBCT scanner. An optically tracked external surrogate was synchronized with CBCT acquisition so that each CBCT projection was time stamped with the surrogate respiratory signal through in-house software and hardware tools. Approximately 2500 projections were acquired over a period of 8-10 minutes in half-fan mode with the half bow-tie filter. Using the external surrogate, the CBCT projections were sorted into 10 breathing phases and reconstructed with an in-house FDK reconstruction algorithm. Errors in respiration sorting, reconstruction, and acquisition were carefully identified and corrected.

DATA FORMAT AND USAGE NOTES

4DCT and 4DCBCT images are available in DICOM format and structures through DICOM-RT RTSTRUCT format. All data are stored in the Cancer Imaging Archive (TCIA, http://www.cancerimagingarchive.net/) as collection 4D-Lung and are publicly available.

DISCUSSION

Due to high temporal frequency sampling, redundant (4DCT and 4DCBCT) data at similar timepoints, oversampled 4DCBCT, and fiducial markers, this dataset can support studies in image-guided and image-guided adaptive radiotherapy, assessment of 4D voxel trajectory variability, and development and validation of new tools for image registration and motion management.

Modelling and Bayesian adaptive prediction of individual patients' tumour volume change during radiotherapy

The aim of this study is to develop a mathematical modelling method that can predict individual patients’ response to radiotherapy, in terms of tumour volume change during the treatment. The main concept is to start from a population-average model, which is subsequently updated from an individual’s tumour volume measurement. The model becomes increasingly personalized and so too does the prediction it produces. This idea of adaptive prediction was realised by using a Bayesian approach for updating the model parameters. The feasibility of the developed method was demonstrated on the data from 25 non-small cell lung cancer patients treated with helical tomotherapy, during which tumour volume was measured from daily imaging as part of the image-guided radiotherapy. The method could provide useful information for adaptive treatment planning and dose scheduling based on the patient’s personalised response.

Time-Resolved Versus Integrated Transit Planar Dosimetry for Volumetric Modulated Arc Therapy

Purpose:

It is desirable that dosimetric deviations during radiation treatments are detected. Integrated transit planar dosimetry is commonly used to evaluate external beam treatments such as volumetric-modulated arc therapy. This work focuses on patient geometry changes which result in differences between the planned and the delivered radiation dose. Integrated transit planar dosimetry will average out some deviations. Novel time-resolved transit planar dosimetry compares the delivered dose of volumetric-modulated arc therapy to the planned dose at various time points. Four patient cases are shown where time-resolved transit planar dosimetry detects patient geometry changes during treatment.

Methods:

A control point to control point comparison between the planned dose and the treatment dose of volumetric-modulated arc therapy beams is calculated using the planning computed tomography and the kV cone-beam computed tomography of the day and evaluated with a time-resolved γ function. Results were computed for 4 patients treated with volumetric-modulated arc therapy, each showing an anatomical change: pleural effusion, rectal gas pockets, and tumor regression.

Results:

In all cases, the geometrical change was detected by time-resolved transit planar dosimetry, whereas integrated transit planar dosimetry showed minor or no indication of the dose discrepancy. Both tumor regression cases were detected earlier in the treatment with time-resolved planar dosimetry in comparison to integrated transit planar dosimetry. The pleural effusion and the gas pocket were detected exclusively with time-resolved transit planar dosimetry.

Conclusions:

Clinical cases were presented in this proof-of-principle study in which integrated transit planar dosimetry did not detect dosimetrically relevant deviations to the same extent time-resolved transit planar dosimetry was able to. Time-resolved transit planar dosimetry also provides results that can be presented as a function of arc delivery angle allowing easier interpretation compared to integrated transit planar dosimetry.

Open source deformable image registration system for treatment planning and recurrence CT scans : Validation in the head and neck region

BACKGROUND

Clinical application of deformable registration (DIR) of medical images remains limited due to sparse validation of DIR methods in specific situations, e. g. in case of cancer recurrences. In this study the accuracy of DIR for registration of planning CT (pCT) and recurrence CT (rCT) images of head and neck squamous cell carcinoma (HNSCC) patients was evaluated.

PATIENTS AND MATERIALS

Twenty patients treated with definitive IMRT for HNSCC in 2010-2012 were included. For each patient, a pCT and an rCT scan were used. Median interval between the scans was 8.5 months. One observer manually contoured eight anatomical regions-of-interest (ROI) twice on pCT and once on rCT.

METHODS

pCT and rCT images were deformably registered using the open source software elastix. Mean surface distance (MSD) and Dice similarity coefficient (DSC) between contours were used for validation of DIR. A measure for delineation uncertainty was estimated by assessing MSD from the re-delineations of the same ROI on pCT. DIR and manual contouring uncertainties were correlated with tissue volume and rigidity.

RESULTS

MSD varied 1-3 mm for different ROIs for DIR and 1-1.5 mm for re-delineated ROIs performed on pCT. DSC for DIR varied between 0.58 and 0.79 for soft tissues and was 0.79 or higher for bony structures, and correlated with the volumes of ROIs (r = 0.5, p < 0.001) and tissue rigidity (r = 0.54, p < 0.001).

CONCLUSION

DIR using elastix in HNSCC on planning and recurrence CT scans is feasible; an uncertainty of the method is close to the voxel size length of the planning CT images.

Cone-beam computed tomography in lung stereotactic ablative radiation therapy: predictive parameters of early response

OBJECTIVE:

To analyze lung lesion volume variations by contouring on cone-beam CT (CBCT) images to evaluate the early predictive parameters of stereotactic ablative radiation therapy (SABR) treatment response.

METHODS:

The prescribed dose of SABR was varied according to the tumour site (central or peripheral) and maximum diameter of the lesions by using a strategy of risk-adapted dose prescription with a dose range between 48 and 70 Gy in 3-10 consecutive fractions. For the purpose of the analysis, the gross tumour volume (GTV) was recontoured for each patient at first and last CBCT using two lung levels/windows: (a) -600/1000 HU and (b) -1000/250 HU. Univariate analysis was performed to evaluate a correlation between lung lesion variations on CBCT using the two levels/windows and treatment response 6 months after SABR. Independent variables were the number of fractions, time between initial and final fraction, biologically effective dose and pre-SABR GTV. Cut points of lesion volume reduction were evaluated to determine the correlation with complete response 6 months after SABR.

RESULTS:

41 lung lesions were evaluated. 82 lung lesions were recontoured for each CBCT level/window. A lung lesion shrinkage of at least 20% was revealed to be statistically related to complete response 6 months after SABR for both the CBCT levels/windows used. The probability of complete response ranged between six and eight times higher in respect to CBCT levels/windows -600/1000 HU and -1000/250 HU, respectively, compared with patients without a lesion shrinkage of 20% at the last session of SABR.

CONCLUSION:

According to current findings, a lung lesion shrinkage of at least 20% at the last session of SABR could be predictable of complete response 6 months thereafter. Further investigations about this topic are needed.

ADVANCES IN KNOWLEDGE:

Prediction of the early tumour response could be useful to personalize imaging restaging after the completion of SABR or to incorporate additional therapies in case of poor responders to improve clinical outcomes.

Voxel-based comparative analysis of lung lesions in CT for therapeutic purposes

Lung cancer remains as one of the most incident types of cancer throughout the world. Temporal evaluation has become a very useful tool when one wishes to analyze some malignancy-indicating behavior. The objective of the present work is to detect changes in the local densities of lung lesions over time (follow-up analysis). From the detected changes, local information as well as extent region of changes can complement the studies regarding the malignant or benign nature of the lesion. Based on this idea, we attempt to use techniques that allow the observation of changes in the lesion over time, based on remote sensing techniques which highlight changes occurring in the environment. The techniques used were the image differencing, image rationing, median filtering, image regression and the fuzzy XOR operator. Based on the global measurement of change percentage in the density, we found density variations which were considered significant in a range from 2.22 to 36.57 % of the volume of the lesion. The results achieved are promising since, besides the visual aspects of the changes in density of the lung lesion over time, we managed to quantify these changes and compare them by volumetric analysis, a more commonly used technique for analysis of changes in lung lesions.

Impact of early tumor reduction on outcome differs by histological subtype in stage III non-small-cell lung cancer treated with definitive radiotherapy

BACKGROUND

We retrospectively investigated the impact on survival of early tumor reduction during definitive radiotherapy for inoperable stage III non-small cell lung cancer (NSCLC) patients, according to their histological subtypes.

METHODS

Between November 2006 and December 2012, 152 consecutive patients with inoperable stage III NSCLC who underwent definitive radiotherapy were reviewed retrospectively. Forty-one patients were excluded for not satisfying the inclusion criteria. Forty-five (40.5 %) and 48 (43.2 %) patients were diagnosed with squamous cell carcinoma (SQC) and adenocarcinoma (ADC), respectively. The tumor reduction rate (TRR) was defined as follows: TRR = 1-[gross tumor volume (GTV) on computed tomography at shrinking irradiation field planning]/(GTV on computed tomography at the initial treatment planning). The Cox proportional hazard model was used to identify significant prognostic factors for overall survival (OS) and progression-free survival (PFS).

RESULTS

We evaluated 111 patients, with a median follow-up time of 52.2 months in surviving patients. The median TRR was 45.9 %. In all patients, there were significant associations between TRR and PFS (P = 0.036) on multivariate analysis, although TRR had no correlation with OS (P = 0.141). With respect to histological subtype, multivariate analyses revealed that a higher TRR showed significant associations with better OS and PFS in the SQC group (P = 0.013 and 0.040, respectively). In contrast, a higher TRR was associated with poorer OS in the ADC group (P = 0.030); there was no association between TRR and PFS.

CONCLUSION

We found that a higher TRR is a promising prognostic factor for better survival and disease control in SQC patients.

Automatic planning of head and neck treatment plans

Treatment planning is time‐consuming and the outcome depends on the person performing the optimization. A system that automates treatment planning could potentially reduce the manual time required for optimization and could also provide a method to reduce the variation between persons performing radiation dose planning (dosimetrist) and potentially improve the overall plan quality. This study evaluates the performance of the Auto‐Planning module that has recently become clinically available in the Pinnacle³ radiation therapy treatment planning system. Twenty‐six clinically delivered head and neck treatment plans were reoptimized with the Auto‐Planning module. Comparison of the two types of treatment plans were performed using DVH metrics and a blinded clinical evaluation by two senior radiation oncologists using a scale from one to six. Both evaluations investigated dose coverage of target and dose to healthy tissues. Auto‐Planning was able to produce clinically acceptable treatment plans in all 26 cases. Target coverages in the two types of plans were similar, but automatically generated plans had less irradiation of healthy tissue. In 94% of the evaluations, the autoplans scored at least as high as the previously delivered clinical plans. For all patients, the Auto‐Planning tool produced clinically acceptable head and neck treatment plans without any manual intervention, except for the initial target and OAR delineations. The main benefit of the method is the likely improvement in the overall treatment quality since consistent, high‐quality plans are generated which even can be further optimized, if necessary. This makes it possible for the dosimetrist to focus more time on difficult dose planning goals and to spend less time on the more tedious parts of the planning process.

PACS number: 87.55.de