A Framework for Assessing the Effect of Cardiac and Respiratory Motion for Stereotactic Arrhythmia Radioablation Using a Digital Phantom With a 17-Segment Model: A STOPSTORM.eu Consortium Study

PURPOSE

The optimal motion management strategy for patients receiving stereotactic arrhythmia radioablation (STAR) for the treatment of ventricular tachycardia (VT) is not fully known. We developed a framework using a digital phantom to simulate cardiorespiratory motion in combination with different motion management strategies to gain insight into the effect of cardiorespiratory motion on STAR.

METHODS AND MATERIALS

The 4-dimensional (4D) extended cardiac-torso (XCAT) phantom was expanded with the 17-segment left ventricular (LV) model, which allowed placement of STAR targets in standardized ventricular regions. Cardiac- and respiratory-binned 4D computed tomography (CT) scans were simulated for free-breathing, reduced free-breathing, respiratory-gating, and breath-hold scenarios. Respiratory motion of the heart was set to population-averaged values of patients with VT: 6, 2, and 1 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction was adjusted by reducing LV ejection fraction to 35%. Target displacement was evaluated for all segments using envelopes encompassing the cardiorespiratory motion. Envelopes incorporating only the diastole plus respiratory motion were created to simulate the scenario where cardiac motion is not fully captured on 4D respiratory CT scans used for radiation therapy planning.

RESULTS

The average volume of the 17 segments was 6 cm3 (1-9 cm3). Cardiac contraction-relaxation resulted in maximum segment (centroid) motion of 4, 6, and 3.5 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction-relaxation resulted in a motion envelope increase of 49% (24%-79%) compared with individual segment volumes, whereas envelopes increased by 126% (79%-167%) if respiratory motion also was considered. Envelopes incorporating only the diastole and respiration motion covered on average 68% to 75% of the motion envelope.

CONCLUSIONS

The developed LV-segmental XCAT framework showed that free-wall regions display the most cardiorespiratory displacement. Our framework supports the optimization of STAR by evaluating the effect of (cardio)respiratory motion and motion management strategies for patients with VT.

Evaluation of a cone-beam computed tomography system calibrated for accurate radiotherapy dose calculation

Background and purpose

Dose calculation on cone-beam computed tomography (CBCT) images has been less accurate than on computed tomography (CT) images due to lower image quality and discrepancies in CT numbers for CBCT. As increasing interest arises in offline and online re-planning, dose calculation accuracy was evaluated for a novel CBCT imager integrated into a ring gantry treatment machine.

Materials and methods

The new CBCT system allowed fast image acquisition (5.9 s) by using new hardware, including a large-size flat panel detector, and incorporated image-processing algorithms with iterative reconstruction techniques, leading to accurate CT numbers allowing dose calculation. In this study, CBCT- and CT-based dose calculations were compared based on three anthropomorphic phantoms, after CBCT-to-mass-density calibration was performed. Six plans were created on the CT scans covering various target locations and complexities, followed by CBCT to CT registrations, copying of contours, and re-calculation of the plans on the CBCT scans. Dose-volume histogram metrics for target volumes and organs-at-risk (OARs) were evaluated, and global gamma analyses were performed.

Results

Target coverage differences were consistently below 1.2 %, demonstrating the agreement between CT and re-calculated CBCT dose distributions. Differences in Dmean for OARs were below 0.5 Gy for all plans, except for three OARs, which were below 0.8 Gy (<1.1 %). All plans had a 3 %/1mm gamma pass rate > 97 %.

Conclusions

This study demonstrated comparable results between dose calculations performed on CBCT and CT acquisitions. The new CBCT system with enhanced image quality and CT number accuracy opens possibilities for off-line and on-line re-planning.

Clinical benefit of range uncertainty reduction in proton treatment planning based on dual-energy CT for neuro-oncological patients

OBJECTIVE

Several studies have shown that dual-energy CT (DECT) can lead to improved accuracy for proton range estimation. This study investigated the clinical benefit of reduced range uncertainty, enabled by DECT, in robust optimisation for neuro-oncological patients.

METHODS

DECT scans for 27 neuro-oncological patients were included. Commercial software was applied to create stopping-power ratio (SPR) maps based on the DECT scan. Two plans were robustly optimised on the SPR map, keeping the beam and plan settings identical to the clinical plan. One plan was robustly optimised and evaluated with a range uncertainty of 3% (as used clinically; denoted 3%-plan); the second plan applied a range uncertainty of 2% (2%-plan). Both plans were clinical acceptable and optimal. The dose-volume histogram parameters were compared between the two plans. Two experienced neuro-radiation oncologists determined the relevant dose difference for each organ-at-risk (OAR). Moreover, the OAR toxicity levels were assessed.

RESULTS

For 24 patients, a dose reduction >0.5/1 Gy (relevant dose difference depending on the OAR) was seen in one or more OARs for the 2%-plan; e.g. for brainstem D0.03cc in 10 patients, and hippocampus D40% in 6 patients. Furthermore, 12 patients had a reduction in toxicity level for one or two OARs, showing a clear benefit for the patient.

CONCLUSION

Robust optimisation with reduced range uncertainty allows for reduction of OAR toxicity, providing a rationale for clinical implementation. Based on these results, we have clinically introduced DECT-based proton treatment planning for neuro-oncological patients, accompanied with a reduced range uncertainty of 2%.

ADVANCES IN KNOWLEDGE

This study shows the clinical benefit of range uncertainty reduction from 3% to 2% in robustly optimised proton plans. A dose reduction to one or more OARs was seen for 89% of the patients, and 44% of the patients had an expected toxicity level decrease.

PET/CT radiomics for prediction of hyperprogression in metastatic melanoma patients treated with immune checkpoint inhibitors

Purpose

This study evaluated pretreatment 2[18F]fluoro-2-deoxy-D-glucose (FDG)-PET/CT-based radiomic signatures for prediction of hyperprogression in metastatic melanoma patients treated with immune checkpoint inhibition (ICI).

Material and method

Fifty-six consecutive metastatic melanoma patients treated with ICI and available imaging were included in the study and 330 metastatic lesions were individually, fully segmented on pre-treatment CT and FDG-PET imaging. Lesion hyperprogression (HPL) was defined as lesion progression according to RECIST 1.1 and doubling of tumor growth rate. Patient hyperprogression (PD-HPD) was defined as progressive disease (PD) according to RECIST 1.1 and presence of at least one HPL. Patient survival was evaluated with Kaplan-Meier curves. Mortality risk of PD-HPD status was assessed by estimation of hazard ratio (HR). Furthermore, we assessed with Fisher test and Mann-Whitney U test if demographic or treatment parameters were different between PD-HPD and the remaining patients. Pre-treatment PET/CT-based radiomic signatures were used to build models predicting HPL at three months after start of treatment. The models were internally validated with nested cross-validation. The performance metric was the area under receiver operating characteristic curve (AUC).

Results

PD-HPD patients constituted 57.1% of all PD patients. PD-HPD was negatively related to patient overall survival with HR=8.52 (95%CI 3.47-20.94). Sixty-nine lesions (20.9%) were identified as progressing at 3 months. Twenty-nine of these lesions were classified as hyperprogressive, thereby showing a HPL rate of 8.8%. CT-based, PET-based, and PET/CT-based models predicting HPL at three months after the start of treatment achieved testing AUC of 0.703 +/- 0.054, 0.516 +/- 0.061, and 0.704 +/- 0.070, respectively. The best performing models relied mostly on CT-based histogram features.

Conclusions

FDG-PET/CT-based radiomic signatures yield potential for pretreatment prediction of lesion hyperprogression, which may contribute to reducing the risk of delayed treatment adaptation in metastatic melanoma patients treated with ICI.

Transferability of radiomic signatures from experimental to human interstitial lung disease

Background

Interstitial lung disease (ILD) defines a group of parenchymal lung disorders, characterized by fibrosis as their common final pathophysiological stage. To improve diagnosis and treatment of ILD, there is a need for repetitive non-invasive characterization of lung tissue by quantitative parameters. In this study, we investigated whether CT image patterns found in mice with bleomycin induced lung fibrosis can be translated as prognostic factors to human patients diagnosed with ILD.

Methods

Bleomycin was used to induce lung fibrosis in mice (n_control = 36, n_experimental = 55). The patient cohort consisted of 98 systemic sclerosis (SSc) patients (n_ILD = 65). Radiomic features (n_histogram = 17, n_texture = 137) were extracted from microCT (mice) and HRCT (patients) images. Predictive performance of the models was evaluated with the area under the receiver-operating characteristic curve (AUC). First, predictive performance of individual features was examined and compared between murine and patient data sets. Second, multivariate models predicting ILD were trained on murine data and tested on patient data. Additionally, the models were reoptimized on patient data to reduce the influence of the domain shift on the performance scores.

Results

Predictive power of individual features in terms of AUC was highly correlated between mice and patients (r = 0.86). A model based only on mean image intensity in the lung scored AUC = 0.921 ± 0.048 in mice and AUC = 0.774 (CI95% 0.677-0.859) in patients. The best radiomic model based on three radiomic features scored AUC = 0.994 ± 0.013 in mice and validated with AUC = 0.832 (CI95% 0.745-0.907) in patients. However, reoptimization of the model weights in the patient cohort allowed to increase the model’s performance to AUC = 0.912 ± 0.058.

Conclusion

Radiomic signatures of experimental ILD derived from microCT scans translated to HRCT of humans with SSc-ILD. We showed that the experimental model of BLM-induced ILD is a promising system to test radiomic models for later application and validation in human cohorts.

Robustness of radiomic features in magnetic resonance imaging for patients with glioblastoma: multi-center study

Background and purpose

Radiomics offers great potential in improving diagnosis and treatment for patients with glioblastoma multiforme. However, in order to implement radiomics in clinical routine, the features used for prognostic modelling need to be stable. This comprises significant challenge in multi-center studies. The aim of this study was to evaluate the impact of different image normalization methods on MRI features robustness in multi-center study.

Methods

Radiomics stability was checked on magnetic resonance images of eleven patients. The images were acquired in two different hospitals using contrast-enhanced T1 sequences. The images were normalized using one of five investigated approaches including grey-level discretization, histogram matching and z-score. Then, radiomic features were extracted and features stability was evaluated using intra-class correlation coefficients. In the second part of the study, improvement in the prognostic performance of features was tested on 60 patients derived from publicly available dataset.

Results

Depending on the normalization scheme, the percentage of stable features varied from 3.4% to 8%. The histogram matching based on the tumor region showed the highest amount of the stable features (113/1404); while normalization using fixed bin size resulted in 48 stable features. The histogram matching also led to better prognostic value (median c-index increase of 0.065) comparing to non-normalized images.

Conclusions

MRI normalization plays an important role in radiomics. Appropriate normalization helps to select robust features, which can be used for prognostic modelling in multicenter studies. In our study, histogram matching based on tumor region improved both stability of radiomic features and their prognostic value.

Improved Survival Prediction by Combining Radiological Imaging and S-100B Levels Into a Multivariate Model in Metastatic Melanoma Patients Treated With Immune Checkpoint Inhibition

Purpose

We explored imaging and blood bio-markers for survival prediction in a cohort of patients with metastatic melanoma treated with immune checkpoint inhibition.

Materials and Methods

94 consecutive metastatic melanoma patients treated with immune checkpoint inhibition were included into this study. PET/CT imaging was available at baseline (Tp0), 3 months (Tp1) and 6 months (Tp2) after start of immunotherapy. Radiological response at Tp2 was evaluated using iRECIST. Total tumor burden (TB) at each time-point was measured and relative change of TB compared to baseline was calculated. LDH, CRP and S-100B were also analyzed. Cox proportional hazards model and logistic regression were used for survival analysis.

Results

iRECIST at Tp2 was significantly associated with overall survival (OS) with C-index=0.68. TB at baseline was not associated with OS, whereas TB at Tp1 and Tp2 provided similar predictive power with C-index of 0.67 and 0.71, respectively. Appearance of new metastatic lesions during follow-up was an independent prognostic factor (C-index=0.73). Elevated LDH and S-100B ratios at Tp2 were significantly associated with worse OS: C-index=0.73 for LDH and 0.73 for S-100B. Correlation of LDH with TB was weak (r=0.34). A multivariate model including TB change, S-100B, and appearance of new lesions showed the best predictive performance with C-index=0.83.

Conclusion

Our analysis shows only a weak correlation between LDH and TB. Additionally, baseline TB was not a prognostic factor in our cohort. A multivariate model combining early blood and imaging biomarkers achieved the best predictive power with regard to survival, outperforming iRECIST.

A 2.5D convolutional neural network for HPV prediction in advanced oropharyngeal cancer

BACKGROUND

Infection with human papilloma virus (HPV) is one of the most relevant prognostic factors in advanced oropharyngeal cancer (OPC) treatment. In this study we aimed to assess the diagnostic accuracy of a deep learning-based method for HPV status prediction in computed tomography (CT) images of advanced OPC.

METHOD

An internal dataset and three public collections were employed (internal: n = 151, HNC1: n = 451; HNC2: n = 80; HNC3: n = 110). Internal and HNC1 datasets were used for training, whereas HNC2 and HNC3 collections were used as external test cohorts. All CT scans were resampled to a 2 mm³ resolution and a sub-volume of 72x72x72 pixels was cropped on each scan, centered around the tumor. Then, a 2.5D input of size 72x72x3 pixels was assembled by selecting the 2D slice containing the largest tumor area along the axial, sagittal and coronal planes, respectively. The convolutional neural network employed consisted of the first 5 modules of the Xception model and a small classification network. Ten-fold cross-validation was applied to evaluate training performance. At test time, soft majority voting was used to predict HPV status.

RESULTS

A final training mean [range] area under the curve (AUC) of 0.84 [0.76-0.89], accuracy of 0.76 [0.64-0.83] and F1-score of 0.74 [0.62-0.83] were achieved. AUC/accuracy/F1-score values of 0.83/0.75/0.69 and 0.88/0.79/0.68 were achieved on the HNC2 and HNC3 test sets, respectively.

CONCLUSION

Deep learning was successfully applied and validated in two external cohorts to predict HPV status in CT images of advanced OPC, proving its potential as a support tool in cancer precision medicine.

Quantification of the spatial distribution of primary tumors in the lung to develop new prognostic biomarkers for locally advanced NSCLC

The anatomical location and extent of primary lung tumors have shown prognostic value for overall survival (OS). However, its manual assessment is prone to interobserver variability. This study aims to use data driven identification of image characteristics for OS in locally advanced non-small cell lung cancer (NSCLC) patients. Five stage IIIA/IIIB NSCLC patient cohorts were retrospectively collected. Patients were treated either with radiochemotherapy (RCT): RCT1* (n = 107), RCT2 (n = 95), RCT3 (n = 37) or with surgery combined with radiotherapy or chemotherapy: S1* (n = 135), S2 (n = 55). Based on a deformable image registration (MIM Vista, 6.9.2.), an in-house developed software transferred each primary tumor to the CT scan of a reference patient while maintaining the original tumor shape. A frequency-weighted cumulative status map was created for both exploratory cohorts (indicated with an asterisk), where the spatial extent of the tumor was uni-labeled with 2 years OS. For the exploratory cohorts, a permutation test with random assignment of patient status was performed to identify regions with statistically significant worse OS, referred to as decreased survival areas (DSA). The minimal Euclidean distance between primary tumor to DSA was extracted from the independent cohorts (negative distance in case of overlap). To account for the tumor volume, the distance was scaled with the radius of the volume-equivalent sphere. For the S1 cohort, DSA were located at the right main bronchus whereas for the RCT1 cohort they further extended in cranio-caudal direction. In the independent cohorts, the model based on distance to DSA achieved performance: AUC_RCT2 [95% CI] = 0.67 [0.55–0.78] and AUC_RCT3 = 0.59 [0.39–0.79] for RCT patients, but showed bad performance for surgery cohort (AUC_S2 = 0.52 [0.30–0.74]). Shorter distance to DSA was associated with worse outcome (p = 0.0074). In conclusion, this explanatory analysis quantifies the value of primary tumor location for OS prediction based on cumulative status maps. Shorter distance of primary tumor to a high-risk region was associated with worse prognosis in the RCT cohort.

Preselection of robust radiomic features does not improve outcome modelling in non-small cell lung cancer based on clinical routine FDG-PET imaging

Background

Radiomics is a promising tool for identifying imaging-based biomarkers. Radiomics-based models are often trained on single-institution datasets; however, multi-centre imaging datasets are preferred for external generalizability owing to the influence of inter-institutional scanning differences and acquisition settings. The study aim was to determine the value of preselection of robust radiomic features in routine clinical positron emission tomography (PET) images to predict clinical outcomes in locally advanced non-small cell lung cancer (NSCLC).

Methods

A total of 1404 primary tumour radiomic features were extracted from pre-treatment [¹⁸F]fluorodeoxyglucose (FDG)-PET scans of stage IIIA/N2 or IIIB NSCLC patients using a training cohort (n = 79; prospective Swiss multi-centre randomized phase III trial SAKK 16/00; 16 centres) and an internal validation cohort (n = 31; single centre). Robustness studies investigating delineation variation, attenuation correction and motion were performed (intraclass correlation coefficient threshold > 0.9). Two 12-/24-month event-free survival (EFS) and overall survival (OS) logistic regression models were trained using standardized imaging: (1) with robust features alone and (2) with all available features. Models were then validated using fivefold cross-validation, and validation on a separate single-centre dataset. Model performance was assessed using area under the receiver operating characteristic curve (AUC).

Results

Robustness studies identified 179 stable features (13%), with 25% stable features for 3D versus 4D acquisition, 31% for attenuation correction and 78% for delineation. Univariable analysis found no significant robust features predicting 12-/24-month EFS and 12-month OS (p value > 0.076). Prognostic models without robust preselection performed well for 12-month EFS in training (AUC = 0.73) and validation (AUC = 0.74). Patient stratification into two risk groups based on 12-month EFS was significant for training (p value = 0.02) and validation cohorts (p value = 0.03).

Conclusions

A PET-based radiomics model using a standardized, multi-centre dataset to predict EFS in locally advanced NSCLC was successfully established and validated with good performance. Prediction models with robust feature preselection were unsuccessful, indicating the need for a standardized imaging protocol.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13550-021-00809-3.

Targeting Treatment Resistance in Head and Neck Squamous Cell Carcinoma - Proof of Concept for CT Radiomics-Based Identification of Resistant Sub-Volumes

Purpose

Radiomics has already been proposed as a prognostic biomarker in head and neck cancer (HNSCC). However, its predictive power in radiotherapy has not yet been studied. Here, we investigated a local radiomics approach to distinguish between tumor sub-volumes with different levels of radiosensitivity as a possible target for radiation dose intensification.

Materials and Methods

Of 40 patients (n=28 training and n=12 validation) with biopsy confirmed locally recurrent HNSCC, pretreatment contrast-enhanced CT images were registered with follow-up PET/CT imaging allowing identification of controlled (GTVcontrol) vs non-controlled (GTVrec) tumor sub-volumes on pretreatment imaging. A bi-regional model was built using radiomic features extracted from pretreatment CT in the GTVrec and GTVcontrol to differentiate between those regions. Additionally, concept of local radiomics was implemented to perform detection task. The original tumor volume was divided into sub-volumes with no prior information on the location of recurrence. Radiomic features from those sub-volumes were then used to detect recurrent sub-volumes using multivariable logistic regression.

Results

Radiomic features extracted from non-controlled regions differed significantly from those in controlled regions (training AUC = 0.79 CI 95% 0.66 - 0.91 and validation AUC = 0.88 CI 95% 0.72 – 1.00). Local radiomics analysis allowed efficient detection of non-controlled sub-volumes both in the training AUC = 0.66 (CI 95% 0.56 – 0.75) and validation cohort 0.70 (CI 95% 0.53 – 0.86), however performance of this model was inferior to bi-regional model. Both models indicated that sub-volumes characterized by higher heterogeneity were linked to tumor recurrence.

Conclusion

Local radiomics is able to detect sub-volumes with decreased radiosensitivity, associated with location of tumor recurrence in HNSCC in the pre-treatment CT imaging. This proof of concept study, indicates that local CT radiomics can be used as predictive biomarker in radiotherapy and potential target for dose intensification.

Radiomic Analysis to Predict Outcome in Recurrent Glioblastoma Based on Multi-Center MR Imaging From the Prospective DIRECTOR Trial

Background

Based on promising results from radiomic approaches to predict O⁶-methylguanine DNA methyltransferase promoter methylation status (MGMT status) and clinical outcome in patients with newly diagnosed glioblastoma, the current study aimed to evaluate radiomics in recurrent glioblastoma patients.

Methods

Pre-treatment MR-imaging data of 69 patients enrolled into the DIRECTOR trial in recurrent glioblastoma served as a training cohort, and 49 independent patients formed an external validation cohort. Contrast-enhancing tumor and peritumoral volumes were segmented on MR images. 180 radiomic features were extracted after application of two MR intensity normalization techniques: fixed number of bins and linear rescaling. Radiomic feature selection was performed via principal component analysis, and multivariable models were trained to predict MGMT status, progression-free survival from first salvage therapy, referred to herein as PFS₂, and overall survival (OS). The prognostic power of models was quantified with concordance index (CI) for survival data and area under receiver operating characteristic curve (AUC) for the MGMT status.

Results

We established and validated a radiomic model to predict MGMT status using linear intensity interpolation and considering features extracted from gadolinium-enhanced T1-weighted MRI (training AUC = 0.670, validation AUC = 0.673). Additionally, models predicting PFS₂ and OS were found for the training cohort but were not confirmed in our validation cohort.

Conclusions

A radiomic model for prediction of MGMT promoter methylation status from tumor texture features in patients with recurrent glioblastoma was successfully established, providing a non-invasive approach to anticipate patient's response to chemotherapy if biopsy cannot be performed. The radiomic approach to predict PFS₂ and OS failed.

Impact of CT convolution kernel on robustness of radiomic features for different lung diseases and tissue types

OBJECTIVES

In this study, we aimed to assess the impact of different CT reconstruction kernels on the stability of radiomic features and the transferability between different diseases and tissue types. Three lung diseases were evaluated, i.e. non-small cell lung cancer (NSCLC), malignant pleural mesothelioma (MPM) and interstitial lung disease related to systemic sclerosis (SSc-ILD) as well as four different tissue types, i.e. primary tumor, largest involved lymph node ipsilateral and contralateral lung.

METHODS

Pre-treatment non-contrast enhanced CT scans from 23 NSCLC, 10 MPM and 12 SSc-ILD patients were collected retrospectively. For each patient, CT scans were reconstructed using smooth and sharp kernel in filtered back projection. The regions of interest (ROIs) were contoured on the smooth kernel-based CT and transferred to the sharp kernel-based CT. The voxels were resized to the largest voxel dimension of each cohort. In total, 1386 features were analyzed. Feature stability was assessed using the intraclass correlation coefficient. Features above the stability threshold >0.9 were considered stable.

RESULTS

We observed a strong impact of the reconstruction method on stability of the features (at maximum 26% of the 1386 features were stable). Intensity features were the most stable followed by texture and wavelet features. The wavelet features showed a positive correlation between percentage of stable features and size of the ROI (R2 = 0.79, p = 0.005). Lymph node radiomics showed poorest stability (<10%) and lung radiomics the largest stability (26%). Robustness analysis done on the contralateral lung could to a large extent be transferred to the ipsilateral lung, and the overlap of stable lung features between different lung diseases was more than 50%. However, results of robustness studies cannot be transferred between tissue types, which was investigated in NSCLC and MPM patients; the overlap of stable features for lymph node and lung, as well as for primary tumor and lymph node was very small in both disease types.

CONCLUSION

The robustness of radiomic features is strongly affected by different reconstruction kernels. The effect is largely influenced by the tissue type and less by the disease type.

ADVANCES IN KNOWLEDGE

The study presents to our knowledge the most complete analysis on the impact of convolution kernel on the robustness of CT-based radiomics for four relevant tissue types in three different lung diseases. .

Comparison of beam segment versus full plan re-optimization in daily magnetic resonance imaging-guided online-adaptive radiotherapy

The optimal approach for magnetic resonance imaging-guided online adaptive radiotherapy is currently unknown and needs to consider patient on-couch time constraints. The aim of this study was to compare two different plan optimization approaches. The comparison was performed in 238 clinically applied online-adapted treatment plans from 55 patients, in which the approach of re-optimization was selected based on the physician’s choice. For 33 patients where both optimization approaches were used at least once, the median treatment planning dose metrics of both target and organ at risk differed less than 1%. Therefore, we concluded that beam segment weight optimization was chosen adequately for most patients without compromising plan quality.

Radiomics Feature Activation Maps as a New Tool for Signature Interpretability

Introduction

In the field of personalized medicine, radiomics has shown its potential to support treatment decisions. However, the limited feature interpretability hampers its introduction into the clinics. Here, we propose a new methodology to create radiomics feature activation maps, which allows to identify the spatial-anatomical locations responsible for signature activation based on local radiomics. The feasibility of this technique will be studied for histological subtype differentiation (adenocarcinoma versus squamous cell carcinoma) in non-small cell lung cancer (NSCLC) using computed tomography (CT) radiomics.

Materials and Methods

Pre-treatment CT scans were collected from a multi-centric Swiss trial (training, n=73, IIIA/N2 NSCLC, SAKK 16/00) and an independent cohort (validation, n=32, IIIA/N2/IIIB NSCLC). Based on the gross tumor volume (GTV), four peritumoral region of interests (ROI) were defined: lung_exterior (expansion into the lung), iso_exterior (expansion into lung and soft tissue), gradient (GTV border region), GTV+Rim (GTV and iso_exterior). For each ROI, 154 radiomic features were extracted using an in-house developed software implementation (Z-Rad, Python v2.7.14). Features robust against delineation variability served as an input for a multivariate logistic regression analysis. Model performance was quantified using the area under the receiver operating characteristic curve (AUC) and verified using five-fold cross validation and internal validation. Local radiomic features were extracted from the GTV+Rim ROI using non-overlapping 3x3x3 voxel patches previously marked as GTV or rim. A binary activation map was created for each patient using the median global feature value from the training. The ratios of activated/non-activated patches of GTV and rim regions were compared between histological subtypes (Wilcoxon test).

Results

Iso_exterior, gradient, GTV+Rim showed good performances for histological subtype prediction (AUC_training=0.68-0.72 and AUC_validation=0.73-0.74) whereas GTV and lung_exterior models failed validation. GTV+Rim model feature activation maps showed that local texture feature distribution differed significantly between histological subtypes in the rim (p=0.0481) but not in the GTV (p=0.461).

Conclusion

In this exploratory study, radiomics-based prediction of NSCLC histological subtypes was predominantly based on the peritumoral region indicating that radiomics activation maps can be useful for tracing back the spatial location of regions responsible for signature activation.

Dosimetric and geometric end-to-end accuracy of a magnetic resonance guided linear accelerator

The introduction of real-time imaging by magnetic resonance guided linear accelerators (MR-Linacs) enabled adaptive treatments and gating on the tumor position. Different end-to-end tests monitored the accuracy of our MR-Linac during the first year of clinical operation. We report on the stability of these tests covering a static, adaptive and gating workflow. Film measurements showed gamma passing rates of 96.4% ± 3.4% for the static tests (five measurements) and for the two adaptive tests 98.9% and 99.99%, respectively (criterion 2%/2mm). The gated point dose measurements in the breathing phantom were 2.7% lower than in the static phantom.

Treatment plan quality during online adaptive re-planning

BACKGROUND

Online adaptive radiotherapy is intended to prevent plan degradation caused by inter-fractional tumor volume and shape changes, but time limitations make online re-planning challenging. The aim of this study was to compare the quality of online-adapted plans to their respective reference treatment plans.

METHODS

Fifty-two patients treated on a ViewRay MRIdian Linac were included in this retrospective study. In total 238 online-adapted plans were analyzed, which were optimized with either changing of the segment weights (n = 85) or full re-optimization (n = 153). Five different treatment sites were evaluated: prostate, abdomen, liver, lung and pelvis. Dosimetric parameters of gross tumor volume (GTV), planning target volume (PTV), 2 cm ring around the PTV and organs at risk (OARs) were considered. The Wilcoxon signed-rank test was used to assess differences between online-adapted and reference treatment plans, p < 0.05 was considered significant.

RESULTS

The average duration of the online adaptation, consisting of contour editing, plan optimization and quality assurance (QA), was 24 ± 6 min. The GTV was slightly larger (average ± SD: 1.9% ± 9.0%) in the adapted plans than in the reference plans (p < 0.001). GTV-D95% exhibited no significant changes when considering all plans, but GTV-D2% increased by 0.40% ± 1.5% on average (p < 0.001). There was a very small yet significant decrease in GTV-coverage for the abdomen plans. The ring Dmean increased on average by 1.0% ± 3.6% considering all plans (p < 0.001). There was a significant reduction of the dose to the rectum of 4.7% ± 16% on average (p < 0.001) for prostate plans.

CONCLUSIONS

Dosimetric quality of online-adapted plans was comparable to reference treatment plans and OAR dose was either comparable or decreased, depending on treatment site. However, dose spillage was slightly increased.

FDG PET versus CT radiomics to predict outcome in malignant pleural mesothelioma patients

BACKGROUND

Careful selection of malignant pleural mesothelioma (MPM) patients for curative treatment is of highest importance, as the multimodal treatment regimen is challenging for patients and harbors a high risk of substantial toxicity. Radiomics-a quantitative method for image analysis-has shown its prognostic ability in different tumor entities and could therefore play an important role in optimizing patient selection for radical cancer treatment. So far, radiomics as a prognostic tool in MPM was not investigated.

MATERIALS AND METHODS

This study is based on 72 MPM patients treated with surgery in a curative intent at our institution between 2009 and 2017. Pre-treatment Fluorine-18 fluorodeoxyglucose (FDG) PET and CT scans were used for radiomics outcome modeling. After extraction of 1404 CT and 1410 FDG PET features from each image, a preselection by principal component analysis was performed to include only robust, non-redundant features for the cox regression to predict the progression-free survival (PFS) and the overall survival (OS). Results were validated on a separate cohort. Additionally, SUVmax and SUVmean, and volume were tested for their prognostic ability for PFS and OS.

RESULTS

For the PFS a concordance index (c-index) of 0.67 (95% CI 0.52-0.82) and 0.66 (95% CI 0.57-0.78) for the training cohort (n = 36) and internal validation cohort (n = 36), respectively, were obtained for the PET radiomics model. The PFS advantage of the low-risk group translated also into an OS advantage. On CT images, no radiomics model could be trained. SUV max and SUV mean were also not prognostic in terms of PFS and OS.

CONCLUSION

We were able to build a successful FDG PET radiomics model for the prediction of PFS in MPM. Radiomics could serve as a tool to aid clinical decision support systems for treatment of MPM in future.

Consolidation cetuximab after concurrent triplet radiochemotherapy+cetuximab in patients with advanced head and neck cancer: A randomized phase II study

BACKGROUND AND PURPOSE

Preclinical data suggest that cetuximab should be continued after end of concurrent radiotherapy+cetuximab due to its efficacy against residual tumor cells in the irradiated tumor bed. Based on this concept the phase II add-on cetuximab (AOC) study was designed.

MATERIALS AND METHODS

Altogether 63 patients with advanced head and neck cancer were treated with radiochemotherapy (70 Gy, cisplatin 40 mg/m² weekly) in combination with concurrent cetuximab (loading dose 400 mg/m², then 250 mg/m² weekly). Thereafter patients were randomized to cetuximab consolidation (500 mg/m² biweekly × 6) or no further treatment. The primary endpoint was the 2-year locoregional control (LRC) rate. As translational research endpoints serum markers were analyzed before and during treatment and CT-based quantitative image analysis (radiomics) was performed.

RESULTS

Median follow-up was 24 months. The 2-year LRC rates were 67.9% and 67.7% in the treatment arms with and without consolidation cetuximab, respectively. Higher than median levels of three serum markers were negatively associated with the 2-year LRC rate in the overall patient cohort: Osteopontin, IL8 and FasL2 (p ≤ 0.05). A radiomics model consisting of two radiomics features could be built showing that higher entropy and higher complexity of tumor Hounsfield unit distribution indicates worse LRC (concordance index 0.66). No correlation was found between biological and imaging markers.

CONCLUSIONS

There was no evidence that consolidation cetuximab would improve the 2-year LRC rate. Prognostic biological and imaging markers could be identified for the overall patient cohort. Studies with larger patient numbers are needed to correlate biological and imaging markers.

The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Image-based Phenotyping

Background Radiomic features may quantify characteristics present in medical imaging. However, the lack of standardized definitions and validated reference values have hampered clinical use. Purpose To standardize a set of 174 radiomic features. Materials and Methods Radiomic features were assessed in three phases. In phase I, 487 features were derived from the basic set of 174 features. Twenty-five research teams with unique radiomics software implementations computed feature values directly from a digital phantom, without any additional image processing. In phase II, 15 teams computed values for 1347 derived features using a CT image of a patient with lung cancer and predefined image processing configurations. In both phases, consensus among the teams on the validity of tentative reference values was measured through the frequency of the modal value and classified as follows: less than three matches, weak; three to five matches, moderate; six to nine matches, strong; 10 or more matches, very strong. In the final phase (phase III), a public data set of multimodality images (CT, fluorine 18 fluorodeoxyglucose PET, and T1-weighted MRI) from 51 patients with soft-tissue sarcoma was used to prospectively assess reproducibility of standardized features. Results Consensus on reference values was initially weak for 232 of 302 features (76.8%) at phase I and 703 of 1075 features (65.4%) at phase II. At the final iteration, weak consensus remained for only two of 487 features (0.4%) at phase I and 19 of 1347 features (1.4%) at phase II. Strong or better consensus was achieved for 463 of 487 features (95.1%) at phase I and 1220 of 1347 features (90.6%) at phase II. Overall, 169 of 174 features were standardized in the first two phases. In the final validation phase (phase III), most of the 169 standardized features could be excellently reproduced (166 with CT; 164 with PET; and 164 with MRI). Conclusion A set of 169 radiomics features was standardized, which enabled verification and calibration of different radiomics software. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Kuhl and Truhn in this issue.

Radiomics, Tumor Volume, and Blood Biomarkers for Early Prediction of Pseudoprogression in Patients with Metastatic Melanoma Treated with Immune Checkpoint Inhibition

PURPOSE

We assessed the predictive potential of positron emission tomography (PET)/CT-based radiomics, lesion volume, and routine blood markers for early differentiation of pseudoprogression from true progression at 3 months.

EXPERIMENTAL DESIGN

112 patients with metastatic melanoma treated with immune checkpoint inhibition were included in our study. Median follow-up duration was 22 months. 716 metastases were segmented individually on CT and 2[18F]fluoro-2-deoxy-D-glucose (FDG)-PET imaging at three timepoints: baseline (TP0), 3 months (TP1), and 6 months (TP2). Response was defined on a lesion-individual level (RECIST 1.1) and retrospectively correlated with FDG-PET/CT radiomic features and the blood markers LDH/S100. Seven multivariate prediction model classes were generated.

RESULTS

Two-year (median) overall survival, progression-free survival, and immune progression-free survival were 69% (not reached), 24% (6 months), and 42% (16 months), respectively. At 3 months, 106 (16%) lesions had progressed, of which 30 (5%) were identified as pseudoprogression at 6 months. Patients with pseudoprogressive lesions and without true progressive lesions had a similar outcome to responding patients and a significantly better 2-year overall survival of 100% (30 months), compared with 15% (10 months) in patients with true progressions/without pseudoprogression (P = 0.002). Patients with mixed progressive/pseudoprogressive lesions were in between at 53% (25 months). The blood prediction model (LDH+S100) achieved an AUC = 0.71. Higher LDH/S100 values indicated a low chance of pseudoprogression. Volume-based models: AUC = 0.72 (TP1) and AUC = 0.80 (delta-volume between TP0/TP1). Radiomics models (including/excluding volume-related features): AUC = 0.79/0.78. Combined blood/volume model: AUC = 0.79. Combined blood/radiomics model (including volume-related features): AUC = 0.78. The combined blood/radiomics model (excluding volume-related features) performed best: AUC = 0.82.

CONCLUSIONS

Noninvasive PET/CT-based radiomics, especially in combination with blood parameters, are promising biomarkers for early differentiation of pseudoprogression, potentially avoiding added toxicity or delayed treatment switch.

Privacy-preserving distributed learning of radiomics to predict overall survival and HPV status in head and neck cancer

A major challenge in radiomics is assembling data from multiple centers. Sharing data between hospitals is restricted by legal and ethical regulations. Distributed learning is a technique, enabling training models on multicenter data without data leaving the hospitals ("privacy-preserving" distributed learning). This study tested feasibility of distributed learning of radiomics data for prediction of two year overall survival and HPV status in head and neck cancer (HNC) patients. Pretreatment CT images were collected from 1174 HNC patients in 6 different cohorts. 981 radiomic features were extracted using Z-Rad software implementation. Hierarchical clustering was performed to preselect features. Classification was done using logistic regression. In the validation dataset, the receiver operating characteristics (ROC) were compared between the models trained in the centralized and distributed manner. No difference in ROC was observed with respect to feature selection. The logistic regression coefficients were identical between the methods (absolute difference <10-7). In comparison of the full workflow (feature selection and classification), no significant difference in ROC was found between centralized and distributed models for both studied endpoints (DeLong p > 0.05). In conclusion, both feature selection and classification are feasible in a distributed manner using radiomics data, which opens new possibility for training more reliable radiomics models.

Systematic Review of Privacy-Preserving Distributed Machine Learning From Federated Databases in Health Care

Big data for health care is one of the potential solutions to deal with the numerous challenges of health care, such as rising cost, aging population, precision medicine, universal health coverage, and the increase of noncommunicable diseases. However, data centralization for big data raises privacy and regulatory concerns.Covered topics include (1) an introduction to privacy of patient data and distributed learning as a potential solution to preserving these data, a description of the legal context for patient data research, and a definition of machine/deep learning concepts; (2) a presentation of the adopted review protocol; (3) a presentation of the search results; and (4) a discussion of the findings, limitations of the review, and future perspectives.Distributed learning from federated databases makes data centralization unnecessary. Distributed algorithms iteratively analyze separate databases, essentially sharing research questions and answers between databases instead of sharing the data. In other words, one can learn from separate and isolated datasets without patient data ever leaving the individual clinical institutes.Distributed learning promises great potential to facilitate big data for medical application, in particular for international consortiums. Our purpose is to review the major implementations of distributed learning in health care.

CT radiomics and PET radiomics: ready for clinical implementation?

INTRODUCTION

Today, rapid technical and clinical developments result in an increasing number of treatment options for oncological diseases. Thus, decision support systems are needed to offer the right treatment to the right patient. Imaging biomarkers hold great promise in patient-individual treatment guidance. Routinely performed for diagnosis and staging, imaging datasets are expected to hold more information than used in the clinical practice. Radiomics describes the extraction of a large number of meaningful quantitative features from medical images, such as computed tomography (CT) and positron emission tomography (PET). Due to the non-invasive nature and ability to capture 3D image-based heterogeneity, radiomic features are potential surrogate markers of the cancer phenotype. Several radiomic studies are published per day, owing to encouraging results of many radiomics-based patient outcome models. Despite this comparably large number of studies, radiomics is mainly studied in proof of principle concept. Hence, a translation of radiomics from a hot topic research field into an essential clinical decision-making tool is lacking, but of high clinical interest.

EVIDENCE ACQUISITION

Herein, we present a literature review addressing the clinical evidence of CT and PET radiomics. An extensive literature review was conducted in PubMed, including papers on robustness and clinical applications.

EVIDENCE SYNTHESIS

We summarize image-modality related influences on the robustness of radiomic features and provide an overview of clinical evidence reported in the literature. Today, more evidence has been provided for CT imaging, however, PET imaging offers the promise of direct imaging of biological processes and functions. We provide a summary of future research directions, which needs to be addressed in order to successfully introduce radiomics into clinical medicine. In comparison to CT, more focus should be directed towards harmonization of PET acquisition and reconstruction protocols, which is important for transferable modelling.

CONCLUSIONS

Both CT and PET radiomics are promising pre-treatment and intra-treatment biomarkers for outcome prediction. Most studies are performed in retrospective setting, however their validation in prospective data collections is ongoing.

Delta-radiomics for prediction of pseudoprogression in malignant melanoma treated with immune checkpoint inhibition

9575

Background: Distinguishing progressive disease (PD) from pseudoprogression (PP) in patients treated with immune-checkpoint inhibition (ICI) is challenging and usually requires confirmation follow-up imaging or invasive diagnostic techniques. This project aimed to identify predictive radiomic signatures for PP from CT imaging. Methods: The response to ICI of 105 metastatic melanoma patients with 645 metastases was retrospectively correlated with radiomic signatures (172 total features). All metastatic lesions were delineated at 3 time points: prior to ICI (t0), at 3 (t1) and 6 months (t2). Response was defined individually for each metastasis using RECIST 1.1, comparing baseline t0 to t2. Three prediction models for PP were built: CT radiomics at t0 and t1, as well as the relative difference between both t0 and t1 (delta-radiomics). Results: Median follow-up was 18 months and 2-year OS and PFS were 72% and 25%, respectively. Median OS: not reached, median PFS: 6 months. Response per lesion at t1: 13% complete remission (CR), 19% partial remission (PR), 52% stable disease (SD) and 16% PD. At t2: 16% CR, 31% PR, 38% SD and 15% PD. 106 progressive lesions were identified at t1, of which, 26 changed to SD, 1 to CR and 3 to PR at t2, resulting in 30 PPs (4.7%). Metastasis location significantly influenced response rates but was not associated with PP (p = 0.4). Lung metastases had significantly higher response rates than soft tissue (p < 0.001), liver (p < 0.001) and bone metastases (p = 0.008). Univariate analysis followed by removal of correlated features revealed no significant radiomic features associated with PP at t0. One independent feature was identified at t1 (AUC 0.74), while delta-radiomics was the best performing approach, identifying four independent features (AUC 0.72 to 0.81). A final multivariate delta radiomics logistic regression model was generated and internally validated, achieving an AUC of 0.81 (± 0.11, 10-fold cross-validation). Conclusions: Metastasis location significantly influenced response rates and CT-based delta-radiomics is a promising biomarker for early differentiation between pseudoprogression and true progression in metastatic melanoma patients treated with ICI.

Interchangeability of radiomic features between [18F]-FDG PET/CT and [18F]-FDG PET/MR

PURPOSE

Radiomics is a promising tool for identification of new prognostic biomarkers. However, image reconstruction settings and test-retest variability may influence the absolute values of radiomic features. Unstable radiomic features cannot be used as reliable biomarkers. PET/MR is becoming increasingly available and often replaces PET/CT for different indications. The aim of this study was to quantify to what extend [18F]-FDG PET/CT radiomics models can be transferred to [18F]-FDG PET/MR and thereby to investigate the feasibility of combined PET/CT-PET/MR models. For this purpose, we compared PET radiomic features calculated on PET/MR and PET/CT and on a 4D-gated PET/MR dataset to select radiomic features that are robust to attenuation correction differences and test-retest variability, respectively.

METHODS

Two cohorts of patients with lung lesions were studied. In the first cohort (n = 10), inhale and exhale phases of a 4D [18F]-FDG PET/MR (4DPETMR) scan were used as a surrogate for a test-retest dataset. In the second cohort (n = 9), patients underwent first an [18F]-FDG PET/MR scan (SIGNA PET/MR, GE Healthcare, Waukesha) followed by an [18F]-FDG PET/CT scan (Discovery 690, GE Healthcare) with a delay of 33 ± 5 min (PETCT-PETMR). Lesions were segmented on inhale and exhale 4D-PET phases and on the individual PET scans from PET/CT and PET/MR with two semi-automated methods (gradient-based and threshold-based). The scan resolution was 2.73 × 2.73 × 3.27 mm and 2.34 × 2.34 × 2.78 mm for the PET/CT and PET/MR, respectively. In total, 1355 radiomic features were calculated, i.e., shape (n = 18), intensity (n = 17), texture (n = 136), and wavelet (n = 1184). The intraclass correlation coefficient (ICC) was calculated to compare the radiomic features of the 4DPETMR (ICC(1,1)) and PETCT-PETMR (ICC(3,1)) datasets. An ICC > 0.9 was considered stable among both types of PET scans.

RESULTS AND CONCLUSION

The 4DPETMR showed highest stability for shape, intensity, and texture (>80%) and lower stability for wavelet features (40%). Gradient-based method showed higher stability compared to threshold-based method except from shape features. In PETCT-PETMR, more than 61% of shape and intensity features were stable for both segmentation methods. However, a reduced stability was observed for texture (50%) and wavelet (<30%) features. More wavelet features were robust in the smoothed images (low-pass filtering) compared to images with emphasized heterogeneity (high-pass filtering). Comparing stable features of both investigations, highest agreement was found for intensity and lower agreement for shape, texture, and wavelet features. Only 53.6% of stable texture features in 4DPETMR were also stable in PETCT-PETMR, and even less in case of wavelet features (40.4%). Approximately 16.9% (texture) and 43.2% (wavelet) of stable PETCT-PETMR features are unstable in 4DPETMR. To conclude, shape and intensity features were robust when comparing two types of [18F]-FDG PET scans (PET/CT and PET/MR). Reduced stability was observed for texture and wavelet features. We identified multiple origins of instability of radiomic features, such as attenuation correction differences, different uptake times, and spatial resolution. This needs to be considered when models based on PET/CT are transferred PET/MR models or when combined models are used.

Body motion during dynamic couch tracking with healthy volunteers

In precision radiotherapy, the intrafractional motion can cause a considerable uncertainty of the location of the tumor to be treated. An established approach is the expansion of the target volume to account for the motion. An alternative approach is couch-tracking, in which the patient is continually moved to compensate the intrafractional motion. However, couch-tracking itself might induce uncertainty of the patient's body position, because the body is non-rigid. One hundred healthy volunteers were positioned supine on a robotic couch. Optical markers were placed on the torso of the volunteers as well as on the couch, and their positions were tracked with an optical surface measurement system. Using these markers, the uncertainty of the body position relative to the couch position was estimated while the couch was static or moving. Over the included 83 healthy volunteers, the median of the uncertainty increased by 0.8 mm (SI), 0.4 mm (LR) and 0.4 mm (AP) when the couch moved. Couch motion was found to increase the uncertainty of the body position relative to the couch. However, this uncertainty is one order of magnitude smaller than the intrafractional tumor motion amplitudes to be compensated. Therefore, even with body motion present, the couch-tracking approach is a viable option. The study was registered at ClinicalTrials.gov (NCT02820532) and the Swiss national clinical trials portal (SNCTP000001878).

Influence of inter-observer delineation variability on radiomics stability in different tumor sites

BACKGROUND

Radiomics is a promising methodology for quantitative analysis and description of radiological images using advanced mathematics and statistics. Tumor delineation, which is still often done manually, is an essential step in radiomics, however, inter-observer variability is a well-known uncertainty in radiation oncology. This study investigated the impact of inter-observer variability (IOV) in manual tumor delineation on the reliability of radiomic features (RF).

METHODS

Three different tumor types (head and neck squamous cell carcinoma (HNSCC), malignant pleural mesothelioma (MPM) and non-small cell lung cancer (NSCLC)) were included. For each site, eleven individual tumors were contoured on CT scans by three experienced radiation oncologists. Dice coefficients (DC) were calculated for quantification of delineation variability. RF were calculated with an in-house developed software implementation, which comprises 1404 features: shape (n = 18), histogram (n = 17), texture (n = 137) and wavelet (n = 1232). The IOV of RF was studied using the intraclass correlation coefficient (ICC). An ICC >0.8 indicates a good reproducibility. For the stable RF, an average linkage hierarchical clustering was performed to identify classes of uncorrelated features.

RESULTS

Median DC was high for NSCLC (0.86, range 0.57-0.90) and HNSCC (0.72, 0.21-0.89), whereas it was low for MPM (0.26, 0-0.9) indicating substantial IOV. Stability rate of RF correlated with DC and depended on tumor site, showing a high stability in NSCLC (90% of total parameters), acceptable stability in HNSCC (59% of total parameters) and low stability in MPM (36% of total parameters). Shape features showed the weakest stability across all tumor types. Hierarchical clustering revealed 14 groups of correlated and stable features for NSCLC and 6 groups for both HNSCC and MPM.

CONCLUSION

Inter-observer delineation variability has a relevant influence on radiomics analysis and is strongly influenced by tumor type. This leads to a reduced number of suitable imaging features.

Development and validation of a radiomic signature to predict HPV (p16) status from standard CT imaging: a multicenter study

OBJECTIVES

Human papillomavirus (HPV) positive oropharyngeal cancer (oropharyngeal squamous cell carcinoma, OPSCC) is biologically and clinically different from HPV negative OPSCC. Here, we evaluate the use of a radiomic approach to identify the HPV status of OPSCC.

METHODS

Four independent cohorts, totaling 778 OPSCC patients with HPV determined by p16 were collected. We randomly assigned 80% of all data for model training (N = 628) and 20% for validation (N = 150). On the pre-treatment CT images, 902 radiomic features were calculated from the gross tumor volume. Multivariable modeling was performed using least absolute shrinkage and selection operator. To assess the impact of CT artifacts in predicting HPV (p16), a model was developed on all training data (Mall) and on the artifact-free subset of training data (Mno art). Models were validated on all validation data (Vall), and the subgroups with (Vart) and without (Vno art) artifacts. Kaplan-Meier survival analysis was performed to compare HPV status based on p16 and radiomic model predictions.

RESULTS

The area under the receiver operator curve for Mall and Mno art ranged between 0.70 and 0.80 and was not significantly different for all validation data sets. There was a consistent and significant split between survival curves with HPV status determined by p16 [p = 0.007; hazard ratio (HR): 0.46], Mall (p = 0.036; HR: 0.55) and Mno art (p = 0.027; HR: 0.49).

CONCLUSION

This study provides proof of concept that molecular information can be derived from standard medical images and shows potential for radiomics as imaging biomarker of HPV status. Advances in knowledge: Radiomics has the potential to identify clinically relevant molecular phenotypes.

Exploratory Radiomics in Computed Tomography Perfusion of Prostate Cancer

BACKGROUND/AIM

An evaluation if radiomic features of CT perfusion (CTP) can predict tumor grade and aggressiveness in prostate cancer was performed.

MATERIALS AND METHODS

Forty-seven patients had biopsy-confirmed prostate cancer, and received a CTP. Blood volume (BV), blood flow (BF) and mean transit time (MTT) maps were derived and 1,701 radiomic features were determined per patient. Regression models were built to estimate post-surgical Gleason score (GS), microvessel density (MVD) and distinguish between the different risk groups.

RESULTS

Six out of the 47 patients had to be excluded from further analysis. A weak relationship between postsurgical GS and one radiomic parameter was found (R²=0.21, p=0.01). The same parameter combined with MTT inter-quartile range was prognostic for the risk group categorisation (AUC=0.81). Two different radiomic parameters were able to distinguish between low-intermediate risk and high-intermediate risk (AUC=0.77). Four parameters correlated with MVD (R2=0.53, p<0.02).

CONCLUSION

This exploratory study shows the potential of radiomics to classify prostate cancer.

Unconscious physiological response of healthy volunteers to dynamic respiration-synchronized couch motion

Background

Intrafractional motion can be a substantial uncertainty in precision radiotherapy. Conventionally, the target volume is expanded to account for the motion. Couch-tracking is an alternative, where the patient is moved to compensate for the tumor motion. However, the couch motion may influence the patient’s stress and respiration behavior decreasing the couch-tracking effectiveness.

Methods

In total, 100 volunteers were positioned supine on a robotic couch, which moved dynamically and respiration synchronized. During the measurement, the skin conductivity, the heartrate, and the gaze location were measured indicating the volunteer’s stress. Volunteers rated the subjective motion sickness using a questionnaire. The measurement alternated between static and tracking segments (three cycles), each 1 min long.

Results

The respiration amplitude showed no significant difference between tracking and static segments, but decreased significantly from the first to the last tracking segment (p < 0.0001). The respiration frequency differed significantly between tracking and static segments (p < 0.0001), but not between the first and the last tracking segment. The physiological parameters and the questionnaire showed mild signals of stress and motion sickness.

Conclusion

Generally, people tolerated the couch motions. The interaction between couch motion and the patient’s breathing pattern should be considered for a clinical implementation.

Trial registration

The study was registered at ClinicalTrials.gov (NCT02820532) and the Swiss national clinical trials portal (SNCTP000001878) on June 20, 2016.

Electronic supplementary material

The online version of this article (10.1186/s13014-017-0925-6) contains supplementary material, which is available to authorized users.