Relating pre-treatment non-Gaussian intravoxel incoherent motion diffusion-weighted imaging to human papillomavirus status and response in oropharyngeal carcinoma

Background and purpose

Diffusion-weighted imaging (DWI) is a promising technique for response assessment in head-and-neck cancer. Recently, we optimized Non-Gaussian Intravoxel Incoherent Motion Imaging (NG-IVIM), an extension of the conventional apparent diffusion coefficient (ADC) model, for the head and neck. In the current study, we describe the first application in a group of patients with human papillomavirus (HPV)-positive and HPV-negative oropharyngeal squamous cell carcinoma. The aim of this study was to relate ADC and NG-IVIM DWI parameters to HPV status and clinical treatment response.

Materials and methods

Thirty-six patients (18 HPV-positive, 18 HPV-negative) were prospectively included. Presence of progressive disease was scored within one year. The mean pre-treatment ADC and NG-IVIM parameters in the gross tumor volume were compared between HPV-positive and HPV-negative patients. In HPV-negative patients, ADC and NG-IVIM parameters were compared between patients with and without progressive disease.

Results

ADC, the NG-IVIM diffusion coefficient D, and perfusion fraction f were significantly higher, while pseudo-diffusion coefficient D* and kurtosis K were significantly lower in the HPV-negative compared to HPV-positive patients. In the HPV-negative group, a significantly lower D was found for patients with progressive disease compared to complete responders. No relation with ADC was observed.

Conclusion

The results of our single-center study suggest that ADC is related to HPV status, but not an independent response predictor. The NG-IVIM parameter D, however, was independently associated to response in the HPV-negative group. Noteworthy in the opposite direction as previously thought based on ADC.

Accounting for fractionation and heterogeneous dose distributions in the modelling of osteoradionecrosis in oropharyngeal carcinoma treatment

BACKGROUND AND PURPOSE

Osteoradionecrosis (ORN) of the mandible is a severe complication following radiotherapy (RT). With a renewed interest in hypofractionation for head and neck radiotherapy, more information concerning ORN development after high fraction doses is important. The aim of this explorative study was to develop a model for ORN risk prediction applicable across different fractionation schemes using Equivalent Uniform Doses (EUD).

MATERIAL AND METHODS

We performed a retrospective cohort study in 334 oropharyngeal squamous cell carcinoma (OPSCC) patients treated with either a hypofractionated Stereotactic Body Radiation Therapy (HF-SBRT) boost or conventional Intensity Modulated Radiation Therapy (IMRT). ORN was scored with the CTCAE v5.0. HF-SBRT and IMRT dose distributions were converted into equivalent dose in 2 Gy fractions (α/β = 0.85 Gy) and analyzed using EUD. The parameter a that led to an EUD that best discriminated patients with and without grade ≥ 2 ORN was selected. Patient and treatment-related risk factors of ORN were analyzed with uni- and multivariable regression analysis.

RESULTS

A total of 32 patients (9.6%) developed ORN grade ≥ 2. An EUD(a = 8) best discriminated between ORN and non-ORN (AUC = 0.71). In multivariable regression, pre-RT extractions (SHR = 2.34; p = 0.012), mandibular volume (SHR = 1.04; p = 0.003), and the EUD(a = 8) (SHR = 1.14; p < 0.001) were significantly associated with ORN.

CONCLUSION

Risk models for ORN based on conventional DVH parameters cannot be directly applied to HF-SBRT fractionation schemes and dose distributions. However, after correcting for fractionation and non-uniform dose distributions using EUD, a single model can distinguish between ORN and non-ORN after conventionally fractionated radiotherapy and hypofractionated boost treatments.

Development of a local dose-response relationship for osteoradionecrosis within the mandible

PURPOSE

Osteoradionecrosis (ORN) of the mandible is a severe complication following radiotherapy of the head and neck, but not all regions of the mandible may be equally at risk. Therefore our goal was to explore a local dose response relationship for subregions of the mandible.

MATERIALS AND METHODS

All oropharyngeal cancer patients treated at our hospital between 2009 and 2016 were reviewed. Follow-up was cut-off at 3 years. For patients that developed ORN, the ORN volume was delineated on the planning CT. Each mandible was divided into 16 volumes of interest (VOIs) based on the location of the dental elements and the presence of ORN in each was scored. Generalized estimating equations were used to build a model for the probability of developing ORN in an element VOI.

RESULTS

Of the 219 included patients, 22 developed ORN in 89 element VOIs. Mean dose to the element VOI (odds ratio (OR) = 1.05 per Gy, 95% confidence interval (CI): (1.04,1.07)), pre-radiotherapy extractions of an element ipsilateral to element of interest (OR = 2.81, 95% CI: (1.12,7.05)), and smoking at start of radiotherapy (OR = 3.37, 95% CI: (1.29,8.78)) were significantly associated with an increased probability of ORN in the VOI.

CONCLUSION

The developed dose-response model indicates that the probability of ORN varies within the mandible and strongly depends on the local dose, the location of extractions, and smoking.

Multidelay pseudocontinuous arterial spin labeling to measure blood flow in the head and neck

Perfusion MRI is promising for the assessment, prediction, and monitoring of radiation toxicity in organs at risk in head and neck cancer. Arterial spin labeling (ASL) may be an attractive alternative for conventional perfusion MRI, that does not require the administration of contrast agents. However, currently, little is known about the characteristics and performance of ASL in healthy tissues in the head and neck region. Therefore, the purpose of this study was to optimize and evaluate multidelay pseudocontinuous ASL (pCASL) for the head and neck region and to explore nominal values and measurement repeatability for the blood flow (BF), and the transit time and T1 values needed for BF quantification in healthy tissues. Twenty healthy volunteers underwent a scan session consisting of four repeats of multidelay pCASL (postlabel delays: 1000, 1632, 2479 ms). Regions of interest were defined in the parotid glands, submandibular glands, tonsils, and the cerebellum (as a reference). Nominal values of BF were calculated as the average over four repeats per volunteer. The repeatability coefficient and within-subject coefficient of repeatability (wCV) of BF were calculated. The effect of T1 (map vs. cohort average) and transit time correction on BF was investigated. The mean BF (± SE) was 55.7 ± 3.1 ml/100 g/min for the parotid glands, 41.2 ± 2.8 ml/100 g/min for the submandibular glands, and 32.3 ± 2.2 ml/100 g/min for the tonsils. The best repeatability was found in the parotid glands (wCV = 13.3%-16.1%), followed by the submandibular glands and tonsils (wCV = 20.0%-24.6%). On average, the effect of T1 and transit time correction on BF was limited, although substantial bias occurred in individual acquisitions. In conclusion, we demonstrated the feasibility of BF measurements in the head and neck region using multidelay pCASL and reported on nominal BF values, BF repeatability, the effect of T1, and transit time in various tissues in the head and neck region.

Post radiation mucosal ulcer risk after a hypofractionated stereotactic boost and conventional fractionated radiotherapy for oropharyngeal carcinoma

BACKGROUND/PURPOSE

Post radiation mucosal ulcers (PRMU) after treatment for oropharyngeal squamous cell carcinoma (OPSCC) can have a huge negative impact on patients' quality of life, but little is known concerning risk factors and the impact of fraction size. Therefore, the goal of this study was to determine the pattern of PRMU development and to identify risk factors after a hypofractionated stereotactic body radiotherapy boost (SBRT) compared to conventionally fractionated radiotherapy for OPSCC.

MATERIAL AND METHODS

We performed a retrospective cohort study (N = 332) of OPSCC patients with ≥ 1-year disease-free survival, treated with 46 Gy Intensity Modulated Radiotherapy (IMRT) (2 Gy fractions) followed by either an SBRT boost of 16.5 Gy (5.5 Gy fractions) (N = 180), or 24 Gy IMRT (2 Gy fractions) (N = 152). PRMU (grade ≥ 2) was scored when observed > three months after the last radiotherapy (RT) fraction (CTCAE v5.0). Potential risk factors were analyzed with Cox regression models using death as competing risk. Dose at the PRMU site was calculated by projecting delineated PRMU on the planning CT.

RESULTS

All cases of PRMU (N = 64) occurred within 24 months; all were grade 2. The cumulative incidence at 2 years in the SBRT boost group was 26% (N = 46) vs. 12% (N = 18) for conventional fractionation (p = 0.003). Most PRMU developed within nine months (N = 48). PRMU occurring > nine months (N = 16) were mainly observed in the SBRT boost group (N = 15). Sex (p = 0.048), acute tube feeding (p = < 0.001), tumor subsite tonsil (p = 0.001), and N stage (p = 0.017) were associated with PRMU risk at multivariable regression in the hypofractionated SBRT boost group. All 25 delineated PRMU were located within the high dose regions.

CONCLUSION

The risk of PRMU should be included in the cost benefit analysis when considering future research using a hypofractionated SBRT boost for OPSCC patients.

The COMPLETE trial: HolistiC early respOnse assessMent for oroPharyngeaL cancEr paTiEnts; Protocol for an observational study

INTRODUCTION

The locoregional failure (LRF) rate in human papilloma virus (HPV)-negative oropharyngeal squamous cell carcinoma (OPSCC) remains disappointingly high and toxicity is substantial. Response prediction prior to or early during treatment would provide opportunities for personalised treatment. Currently, there are no accurate predictive models available for correct OPSCC patient selection. Apparently, the pivotal driving forces that determine how a OPSCC responds to treatment, have yet to be elucidated. Therefore, the holistiC early respOnse assessMent for oroPharyngeaL cancer paTiEnts study focuses on a holistic approach to gain insight in novel potential prognostic biomarkers, acquired before and early during treatment, to predict response to treatment in HPV-negative patients with OPSCC.

METHODS AND ANALYSIS

This single-centre prospective observational study investigates 60 HPV-negative patients with OPSCC scheduled for primary radiotherapy (RT) with cisplatin or cetuximab, according to current clinical practice. A holistic approach will be used that aims to map the macroscopic (with Intra Voxel Incoherent Motion Diffusion Kurtosis Imaging (IVIM-DKI); before, during, and 3 months after RT), microscopic (with biopsies of the primary tumour acquired before treatment and irradiated ex vivo to assess radiosensitivity), and molecular landscape (with circulating tumour DNA (ctDNA) analysed before, during and 3 months after treatment). The main end point is locoregional control (LRC) 2 years after treatment. The primary objective is to determine whether a relative change in the mean of the diffusion coefficient D (an IVIM-DKI parameter) in the primary tumour early during treatment, improves the performance of a predictive model consisting of tumour volume only, for 2 years LRC after treatment. The secondary objectives investigate the potential of other IVIM-DKI parameters, ex vivo sensitivity characteristics, ctDNA, and combinations thereof as potential novel prognostic markers.

ETHICS AND DISSEMINATION

The study was approved by the Medical Ethical Committee of Erasmus Medical Center. The main results of the trial will be presented in international meetings and medical journals.

TRIAL REGISTRATION NUMBER

NL8458.

Comparison of 81mKrypton and 99mTc-Technegas for ventilation single-photon emission computed tomography in severe chronic obstructive pulmonary disease

INTRODUCTION

Ventilation and perfusion single-photon emission computed tomography combined with computed tomography (SPECT/CT) is a powerful tool to assess the state of the lungs in chronic obstructive pulmonary disease (COPD). 81mKrypton is a gaseous ventilation tracer and distributes similarly to air, but is not widely available and relatively expensive. 99mTc-Technegas is cheaper and has wider availability, but is an aerosol, which may deposit in hot spots as the severity of COPD increases. In this study, 81mKrypton and 99mTc-Technegas were compared quantitatively in patients with severe COPD.

METHODS

The penetration ratio, the heterogeneity index (with and without band filtering for relevant clinical sizes) and hot spot appearance were assessed in eleven patients with severe COPD that underwent simultaneous dual-isotope ventilation SPECT/CT with both 99mTc-Technegas and 81mKrypton.

RESULTS

Significant differences were found in the penetration ratio for the medium energy general purpose (MEGP) collimators, but not for the low energy general purpose (LEGP) collimators. The difference in the overall and the band filtered heterogeneity index was significant in most cases. All patients suffered from 99mTc-Technegas hot spots in at least one lung. Comparison of MEGP 81mKrypton and LEGP Technegas scans revealed similar results as the comparison for the MEGP collimators.

CONCLUSION

Caution should be taken when replacing 81mKrypton with 99mTc-Technegas as a ventilation tracer in patients with severe COPD as there are significant differences in the distribution of the tracers over the lungs. Furthermore, this patient group is prone to 99mTc-Technegas hot spots and might need additional scanning if hot spots severely hamper image interpretation.

An optimal acquisition and post-processing pipeline for hybrid IVIM-DKI in head and neck

Purpose

To optimize the diffusion‐weighting b values and postprocessing pipeline for hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region.

Methods

Optimized diffusion‐weighting b value sets ranging between 5 and 30 b values were constructed by optimizing the Cramér‐Rao lower bound of the hybrid intravoxel incoherent motion diffusion kurtosis imaging model. With this model, the perfusion fraction, pseudodiffusion coefficient, diffusion coefficient, and kurtosis were estimated. Sixteen volunteers were scanned with a reference b value set and 3 repeats of the optimized sets, of which 1 with volunteers swallowing on purpose.

The effects of (1) b value optimization and number of b values, (2) registration type (none vs. intervolume vs. intra‐ and intervolume registration), and (3) manual swallowing artifact rejection on the parameter precision were assessed.

Results

The SD was higher in the reference set for perfusion fraction, diffusion coefficient, and kurtosis by a factor of 1.7, 1.5, and 2.3 compared to the optimized set, respectively. A smaller SD (factor 0.7) was seen in pseudodiffusion coefficient. The sets containing 15, 20, and 30 b values had comparable repeatability in all parameters, except pseudodiffusion coefficient, for which set size 30 was worse. Equal repeatability for the registration approaches was seen in all parameters of interest. Swallowing artifact rejection removed the bias when present.

Conclusion

To achieve optimal hybrid intravoxel incoherent motion diffusion kurtosis imaging in the head and neck region, b value optimization and swallowing artifact image rejection are beneficial. The optimized set of 15 b values yielded the optimal protocol efficiency, with a precision comparable to larger b value sets and a 50% reduction in scan time.