Assessment of gene expressions from squamous cell carcinoma of the head and neck to predict radiochemotherapy-related xerostomia and dysphagia.

PURPOSE: We tested the hypothesis that gene expressions from biopsies of locally advanced head and neck squamous cell carcinoma (HNSCC) patients can supplement dose-volume parameters to predict dysphagia and xerostomia following primary radiochemotherapy (RCTx). MATERIAL AND METHODS: A panel of 178 genes previously related to radiochemosensitivity of HNSCC was considered for nanoString analysis based on tumour biopsies of 90 patients with locally advanced HNSCC treated by primary RCTx. Dose-volume parameters were extracted from the parotid, submandibular glands, oral cavity, larynx, buccal mucosa, and lips. Normal tissue complication probability (NTCP) models were developed for acute, late, and for the improvement of xerostomia grade ≥2 and dysphagia grade ≥3 using a cross-validation-based least absolute shrinkage and selection operator (LASSO) approach combined with stepwise logistic regression for feature selection. The final signatures were included in a logistic regression model with optimism correction. Performance was assessed by the area under the receiver operating characteristic curve (AUC). RESULTS: NTCP models for acute and late xerostomia and the improvement of dysphagia resulted in optimism-corrected AUC values of 0.84, 0.76, and 0.70, respectively. The minimum dose to the contralateral parotid was selected for both acute and late xerostomia and the minimum dose to the larynx was selected for dysphagia improvement. For the xerostomia endpoints, the following gene expressions were selected: RPA2 (cellular response to DNA damage), TCF3 (salivary gland cells development), GBE1 (glycogen storage and regulation), and MAPK3 (regulation of cellular processes). No gene expression features were selected for the prediction of dysphagia. CONCLUSION: This hypothesis-generating study showed the potential of improving NTCP models using gene expression data for HNSCC patients. The presented models require independent validation before potential application in clinical practice.

Experimental validation of 4D log file-based proton dose reconstruction for interplay assessment considering amplitude-sorted 4DCTs.

PURPOSE: The unpredictable interplay between dynamic proton therapy delivery and target motion in the thorax can lead to severe dose distortions. A fraction-wise four-dimensional (4D) dose reconstruction workflow allows for the assessment of the applied dose after patient treatment while considering the actual beam delivery sequence extracted from machine log files, the recorded breathing pattern and the geometric information from a 4D computed tomography scan (4DCT). Such an algorithm capable of accounting for amplitude-sorted 4DCTs was implemented and its accuracy as well as its sensitivity to input parameter variations was experimentally evaluated. METHODS: An anthropomorphic thorax phantom with a movable insert containing a target surrogate and a radiochromic film was irradiated with a monoenergetic field for various 1D target motion forms (sin, sin4 ) and peak-to-peak amplitudes (5/10/15/20/30 mm). The measured characteristic film dose distributions were compared to the respective sections in the 4D reconstructed doses using a 2D γ-analysis (3 mm, 3%); γ-pass rates were derived for different dose grid resolutions (1 mm/3 mm) and deformable image registrations (DIR, automatic/manual) applied during the 4D dose reconstruction process. In an additional analysis, the sensitivity of reconstructed dose distributions against potential asynchronous timing of the motion and machine log files was investigated for both a monoenergetic field and more realistic 4D robustly optimized fields by artificially introduced offsets of ±1/5/25/50/250 ms. The resulting dose distributions with asynchronized log files were compared to those with synchronized log files by means of a 3D γ-analysis (1 mm, 1%) and the evaluation of absolute dose differences. RESULTS: The induced characteristic interplay patterns on the films were well reproduced by the 4D dose reconstruction with 2D γ-pass rates ≥95% for almost all cases with motion magnitudes ≤15 mm. In general, the 2D γ-pass rates showed a significant decrease for larger motion amplitudes and increase when using a finer dose grid resolution but were not affected by the choice of motion form (sin, sin4 ). There was also a trend, though not statistically significant, toward the manually defined DIR for better quality of the reconstructed dose distributions in the area imaged by the film. The 4D dose reconstruction results for the monoenergetic as well as the 4D robustly optimized fields were robust against small asynchronies between motion and machine log files of up to 5 ms, which is in the order of potential network latencies. CONCLUSIONS: We have implemented a 4D log file-based proton dose reconstruction that accounts for amplitude-sorted 4DCTs. Its accuracy was proven to be clinically acceptable for target motion magnitudes of up to 15 mm. Particular attention should be paid to the synchronization of the log file generating systems as the reconstructed dose distribution may vary with log file asynchronies larger than those caused by realistic network delays.

Sample-size calculation for preclinical dose-response experiments using heterogeneous tumour models.

BACKGROUND: In preclinical radio-oncological research, local tumour control is considered the most relevant endpoint as it reflects the inactivation of cancer stem cells. Preclinical tumour-control assays may compare dose-response curves between different radiotherapy strategies, e.g., assessing additional targeted drugs and immunotherapeutic interventions, or between different radiation modalities. To mimic the biological heterogeneity of human tumour populations and to accommodate for approaches of personalized oncology, preclinical studies are increasingly performed combining larger panels of tumour models. For designing the study protocols and to obtain reliable results, prospective sample-size planning has to be developed that accounts for such heterogeneous cohorts. METHODS: A Monte-Carlo-based method was developed to estimate the sample size of a comparative 1:1 two-arm prospective tumour-control assay. Based on repeated logistic regression analysis, pre-defined dose levels, assumptions on the dose-response curves of the included tumour models and on the dose-modifying factors (DMF), the power is calculated for a given number of animals per dose group. RESULTS: Two applications are presented: (i) For a simple tumour-control assay with the head and neck squamous cell carcinoma (HNSCC) model FaDu, 10 animals would be required for each of 7 dose levels per arm to reveal a DMF of 1.25 with a power of 0.82 without drop out (total: 140 animals). (ii) In a more complex experiment combining six different lung tumour models to a heterogeneous population, 21 animals would be required for each of 11 dose levels per arm to reveal a DMF of 1.25 with a power of 0.81 without drop out (total: 462 animals). Analyzing the heterogeneous cohort reduces the required animal number by more than 50% compared to six individual tumour-control assays. CONCLUSION: An approach for estimating the required animal number for comparative tumour-control assays in a heterogeneous population is presented, allowing also the inclusion of different treatments as a personalized approach in the experimental arm. The software is publicly available and can be applied to plan comparisons of sigmoidal dose-response curves based on logistic regression.