Advances in Imaging for HPV-Related Oropharyngeal Cancer: Applications to Radiation Oncology

Determinants of the level of circulating-tumor HPV16 DNA in patients with HPV-associated oropharyngeal cancer at the time of diagnosis

Circulating tumor HPV DNA (ctHPV16) assessed in liquid biopsy may be used as a marker of cancer in patients with HPV-associated oropharyngeal cancer (HPV + OPC). Factors influencing the initial ctHPV16 quantity are not well recognized. In this study we aimed to establish what factors are related to the level of ctHPV16 at the time of diagnosis. 51 patients (37 men and 14 women, median age of 57 years old) with HPV + OPC prior to definitive treatment were included. ctHPV16 was measured by qPCR. Tumor and nodal staging were assessed according to AJCC8. Blood derived factors included squamous cell carcinoma antigen (SCC-Ag), serum soluble fragment of cytokeratin 19 (CYFRA 21-1), C-reactive protein (CRP), albumin level (Alb), neutrophils (Neut), thrombocytes (Plt) and lymphocyte (Lym) count, Neut/Lym ratio were assessed. The volumes of the primary tumor (TV) and involved lymph nodes (NV) were calculated using MRI, CT or PET-CT scans. Data were analysed using parametric and nonparametric methods. Variables for multivariable linear regression analysis were chosen based on the results from univariable analysis (correlation, univariable regression and difference). There were 9 (18%), 10 (19%) and 32 (63%) patients who had TV and NV assessed in MRI, CT or PET respectively. Primary tumor neither as T-stage nor TV was related to ctHPV16 level. Significant differences in the ctHPV16 between patients with high vs low pain (P = 0.038), NV (P = 0.023), TV + NV (P = 0.018), CYFRA 21-1 (P = 0.002), CRP (P = 0.019), and N1 vs N3 (P = 0.044) were observed. ctHPV16 was significantly associated with CYFRA 21-1 (P = 0.017), N stage (P = 0.005), NV (P = 0.009), TV + NV (P = 0.002), CRP (P = 0.019), and pain (P = 0.038). In univariable linear regression analysis the same variables predicted ctHPV16 level. In multivariable analyses, CYFRA 21-1 and CRP (both as categorical variables) were predictors of ctHPV16 level even above NV. ctHPV16 at presentation is driven by tumor volume measured mostly by N. CYFRA 21-1 and CRP are additional factors related to ctHPV16 prior to the treatment.

Imaging of human papilloma virus (HPV) related oropharynx tumour: what we know to date

The tumours of head and neck district are around 3% of all malignancies and squamous cell carcinoma is the most frequent histotype, with rapid increase during the last two decades because of the increment of the infection due to human papilloma virus (HPV). Even if the gold standard for the diagnosis is histological examination, including the detection of viral DNA and transcription products, imaging plays a fundamental role in the detection and staging of HPV + tumours, in order to assess the primary tumour, to establish the extent of disease and for follow-up. The main diagnostic tools are Computed Tomography (CT), Positron Emission Tomography-Computed Tomography (PET-CT) and Magnetic Resonance Imaging (MRI), but also Ultrasound (US) and the use of innovative techniques such as Radiomics have an important role. Aim of our review is to illustrate the main imaging features of HPV + tumours of the oropharynx, in US, CT and MRI imaging. In particular, we will outline the main limitations and strengths of the various imaging techniques, the main uses in the diagnosis, staging and follow-up of disease and the fundamental differential diagnoses of this type of tumour. Finally, we will focus on the innovative technique of texture analysis, which is increasingly gaining importance as a diagnostic tool in aid of the radiologist.

Longitudinal diffusion and volumetric kinetics of head and neck cancer magnetic resonance on a 1.5 T MR-linear accelerator hybrid system: A prospective R-IDEAL stage 2a imaging biomarker characterization/pre-qualification study

Objectives

We aim to characterize the serial quantitative apparent diffusion coefficient (ADC) changes of the target disease volume using diffusion-weighted imaging (DWI) acquired weekly during radiation therapy (RT) on a 1.5 T MR-Linac and correlate these changes with tumor response and oncologic outcomes for head and neck squamous cell carcinoma (HNSCC) patients as part of a programmatic R-IDEAL biomarker characterization effort.

Methods

Thirty patients with HNSCC who received curative-intent RT at MD Anderson Cancer Center, were included. Baseline and weekly MRI were obtained, and various ADC parameters were extracted from the regions of interest (ROIs). Baseline and weekly ADC parameters were correlated with response during and after RT, and the recurrence using the Mann-Whitney U test. The Wilcoxon signed-rank test was used to compare the weekly ADC versus baseline values. Weekly volumetric changes (Δvolume) for each ROI were correlated with ΔADC using Spearman's Rho test. Recursive partitioning analysis (RPA) identified the optimal ΔADC threshold associated with different oncologic outcomes.

Results

There was a significant rise in all ADC parameters at different time points of RT compared to baseline for both gross primary disease (GTV-P) and gross nodal disease volumes (GTV-N). The increased ADC values for GTV-P were statistically significant only for primary tumors achieving complete remission (CR) during RT. RPA identified GTV-P ΔADC 5th percentile > 13% at the mid-RT as the most significant parameter associated with primary tumors' CR during RT (p < 0.001). There was a significant decrease in residual volume of both GTV-P & GTV-N throughout the course of RT. A significant negative correlation between mean ΔADC and Δvolume for GTV-P at the 3rd and 4th week of RT was detected (r = -0.39, p = 0.044 & r = -0.45, p = 0.019, respectively).

Conclusion

Assessment of ADC kinetics at regular intervals throughout RT seems to be correlated with RT response. Further studies with larger cohorts and multi-institutional data are needed for validation of ΔADC as a model for prediction of response to RT.

Prospective evaluation of in vivo and phantom repeatability and reproducibility of diffusion-weighted MRI sequences on 1.5 T MRI-linear accelerator (MR-Linac) and MR simulator devices for head and neck cancers

INTRODUCTION

Diffusion-weighted imaging (DWI) on MRI-linear accelerator (MR-linac) systems can potentially be used for monitoring treatment response and adaptive radiotherapy in head and neck cancers (HNC) but requires extensive validation. We performed technical validation to compare six total DWI sequences on an MR-linac and MR simulator (MR sim) in patients, volunteers, and phantoms.

METHODS

Ten human papillomavirus-positive oropharyngeal cancer patients and ten healthy volunteers underwent DWI on a 1.5 T MR-linac with three DWI sequences: echo planar imaging (EPI), split acquisition of fast spin echo signals (SPLICE), and turbo spin echo (TSE). Volunteers were also imaged on a 1.5 T MR sim with three sequences: EPI, BLADE (vendor tradename), and readout segmentation of long variable echo trains (RESOLVE). Participants underwent two scan sessions per device and two repeats of each sequence per session. Repeatability and reproducibility within-subject coefficient of variation (wCV) of mean ADC were calculated for tumors and lymph nodes (patients) and parotid glands (volunteers). ADC bias, repeatability/reproducibility metrics, SNR, and geometric distortion were quantified using a phantom.

RESULTS

In vivo repeatability/reproducibility wCV for parotids were 5.41%/6.72%, 3.83%/8.80%, 5.66%/10.03%, 3.44%/5.70%, 5.04%/5.66%, 4.23%/7.36% for EPIMR-linac, SPLICE, TSE, EPIMR sim, BLADE, RESOLVE. Repeatability/reproducibility wCV for EPIMR-linac, SPLICE, TSE were 9.64%/10.28%, 7.84%/8.96%, 7.60%/11.68% for tumors and 7.80%/9.95%, 7.23%/8.48%, 10.82%/10.44% for nodes. All sequences except TSE had phantom ADC biases within ± 0.1x10-3 mm2/s for most vials (EPIMR-linac, SPLICE, and BLADE had 2, 3, and 1 vials out of 13 with larger biases, respectively). SNR of b = 0 images was 87.3, 180.5, 161.3, 171.0, 171.9, 130.2 for EPIMR-linac, SPLICE, TSE, EPIMR sim, BLADE, RESOLVE.

CONCLUSION

MR-linac DWI sequences demonstrated near-comparable performance to MR sim sequences and warrant further clinical validation for treatment response assessment in HNC.

Magnetic resonance imaging in naso-oropharyngeal carcinoma: role of texture analysis in the assessment of response to radiochemotherapy, a preliminary study

OBJECTIVE

Identifying MRI texture parameters able to distinguish inflammation, fibrosis, and residual cancer in patients with naso-oropharynx carcinoma after radiochemotherapy (RT-CHT).

MATERIAL AND METHODS

In this single-centre, observational, retrospective study, texture analysis was performed on ADC maps and post-gadolinium T1 images of patients with histological diagnosis of naso-oropharyngeal carcinoma treated with RT-CHT. An initial cohort of 99 patients was selected; 57 of them were later excluded. The final cohort of 42 patients was divided into 3 groups (inflammation, fibrosis, and residual cancer) according to MRI, 18F-FDG-PET/CT performed 3-4 months after RT-CHT, and biopsy. Pre-RT-CHT lesions and the corresponding anatomic area post-RT-CHT were segmented with 3D slicer software from which 107 textural features were derived. T-Student and Wilcoxon signed-rank tests were performed, and features with p-value < 0.01 were considered statistically significant. Cut-off values-obtained by ROC curves-to discriminate post-RT-CHT non-tumoural changes from residual cancer were calculated for the parameters statistically associated to the diseased status at follow-up.

RESULTS

Two features-Energy and Grey Level Non-Uniformity-were statistically significant on T1 images in the comparison between 'positive' (residual cancer) and 'negative' patients (inflammation and fibrosis). Energy was also found to be statistically significant in both patients with fibrosis and residual cancer. Grey Level Non-Uniformity was significant in the differentiation between residual cancer and inflammation. Five features were statistically significant on ADC maps in the differentiation between 'positive' and 'negative' patients. The reduction in values of such features between pre- and post-RT-CHT was correlated with a good response to therapy.

CONCLUSIONS

Texture analysis on post-gadolinium T1 images and ADC maps can differentiate residual cancer from fibrosis and inflammation in early follow-up of naso-oropharyngeal carcinoma treated with RT-CHT.

Longitudinal diffusion and volumetric kinetics of head and neck cancer magnetic resonance on a 1.5T MR-Linear accelerator hybrid system: A prospective R-IDEAL Stage 2a imaging biomarker characterization/ pre-qualification study

Objectives

We aim to characterize the serial quantitative apparent diffusion coefficient (ADC) changes of the target disease volume using diffusion-weighted imaging (DWI) acquired weekly during radiation therapy (RT) on a 1.5T MR-Linac and correlate these changes with tumor response and oncologic outcomes for head and neck squamous cell carcinoma (HNSCC) patients as part of a programmatic R-IDEAL biomarker characterization effort.

Methods

Thirty patients with pathologically confirmed HNSCC who received curative-intent RT at the University of Texas MD Anderson Cancer Center, were included in this prospective study. Baseline and weekly Magnetic resonance imaging (MRI) (weeks 1-6) were obtained, and various ADC parameters (mean, 5 th , 10 th , 20 th , 30 th , 40 th , 50 th , 60 th , 70 th , 80 th , 90 th and 95 th percentile) were extracted from the target regions of interest (ROIs). Baseline and weekly ADC parameters were correlated with response during RT, loco-regional control, and the development of recurrence using the Mann-Whitney U test. The Wilcoxon signed-rank test was used to compare the weekly ADC versus baseline values. Weekly volumetric changes (Δvolume) for each ROI were correlated with ΔADC using Spearman's Rho test. Recursive partitioning analysis (RPA) was performed to identify the optimal ΔADC threshold associated with different oncologic outcomes.

Results

There was an overall significant rise in all ADC parameters during different time points of RT compared to baseline values for both gross primary disease volume (GTV-P) and gross nodal disease volumes (GTV-N). The increased ADC values for GTV-P were statistically significant only for primary tumors achieving complete remission (CR) during RT. RPA identified GTV-P ΔADC 5 th percentile >13% at the 3 rd week of RT as the most significant parameter associated with CR for primary tumor during RT (p <0.001). Baseline ADC parameters for GTV-P and GTV-N didn't significantly correlate with response to RT or other oncologic outcomes. There was a significant decrease in residual volume of both GTV-P & GTV-N throughout the course of RT. Additionally, a significant negative correlation between mean ΔADC and Δvolume for GTV-P at the 3 rd and 4 th week of RT was detected (r = -0.39, p = 0.044 & r = -0.45, p = 0.019, respectively).

Conclusion

Assessment of ADC kinetics at regular intervals throughout RT seems to be correlated with RT response. Further studies with larger cohorts and multi-institutional data are needed for validation of ΔADC as a model for prediction of response to RT.

Application of simultaneous uncertainty quantification for image segmentation with probabilistic deep learning: Performance benchmarking of oropharyngeal cancer target delineation as a use-case

Background

Oropharyngeal cancer (OPC) is a widespread disease, with radiotherapy being a core treatment modality. Manual segmentation of the primary gross tumor volume (GTVp) is currently employed for OPC radiotherapy planning, but is subject to significant interobserver variability. Deep learning (DL) approaches have shown promise in automating GTVp segmentation, but comparative (auto)confidence metrics of these models predictions has not been well-explored. Quantifying instance-specific DL model uncertainty is crucial to improving clinician trust and facilitating broad clinical implementation. Therefore, in this study, probabilistic DL models for GTVp auto-segmentation were developed using large-scale PET/CT datasets, and various uncertainty auto-estimation methods were systematically investigated and benchmarked.

Methods

We utilized the publicly available 2021 HECKTOR Challenge training dataset with 224 co-registered PET/CT scans of OPC patients with corresponding GTVp segmentations as a development set. A separate set of 67 co-registered PET/CT scans of OPC patients with corresponding GTVp segmentations was used for external validation. Two approximate Bayesian deep learning methods, the MC Dropout Ensemble and Deep Ensemble, both with five submodels, were evaluated for GTVp segmentation and uncertainty performance. The segmentation performance was evaluated using the volumetric Dice similarity coefficient (DSC), mean surface distance (MSD), and Hausdorff distance at 95% (95HD). The uncertainty was evaluated using four measures from literature: coefficient of variation (CV), structure expected entropy, structure predictive entropy, and structure mutual information, and additionally with our novel Dice-risk measure. The utility of uncertainty information was evaluated with the accuracy of uncertainty-based segmentation performance prediction using the Accuracy vs Uncertainty (AvU) metric, and by examining the linear correlation between uncertainty estimates and DSC. In addition, batch-based and instance-based referral processes were examined, where the patients with high uncertainty were rejected from the set. In the batch referral process, the area under the referral curve with DSC (R-DSC AUC) was used for evaluation, whereas in the instance referral process, the DSC at various uncertainty thresholds were examined.

Results

Both models behaved similarly in terms of the segmentation performance and uncertainty estimation. Specifically, the MC Dropout Ensemble had 0.776 DSC, 1.703 mm MSD, and 5.385 mm 95HD. The Deep Ensemble had 0.767 DSC, 1.717 mm MSD, and 5.477 mm 95HD. The uncertainty measure with the highest DSC correlation was structure predictive entropy with correlation coefficients of 0.699 and 0.692 for the MC Dropout Ensemble and the Deep Ensemble, respectively. The highest AvU value was 0.866 for both models. The best performing uncertainty measure for both models was the CV which had R-DSC AUC of 0.783 and 0.782 for the MC Dropout Ensemble and Deep Ensemble, respectively. With referring patients based on uncertainty thresholds from 0.85 validation DSC for all uncertainty measures, on average the DSC improved from the full dataset by 4.7% and 5.0% while referring 21.8% and 22% patients for MC Dropout Ensemble and Deep Ensemble, respectively.

Conclusion

We found that many of the investigated methods provide overall similar but distinct utility in terms of predicting segmentation quality and referral performance. These findings are a critical first-step towards more widespread implementation of uncertainty quantification in OPC GTVp segmentation.

Advanced Techniques in Head and Neck Cancer Imaging: Guide to Precision Cancer Management

Precision treatment requires precision imaging. With the advent of various advanced techniques in head and neck cancer treatment, imaging has become an integral part of the multidisciplinary approach to head and neck cancer care from diagnosis to staging and also plays a vital role in response evaluation in various tumors. Conventional anatomic imaging (CT scan, MRI, ultrasound) remains basic and focuses on defining the anatomical extent of the disease and its spread. Accurate assessment of the biological behavior of tumors, including tumor cellularity, growth, and response evaluation, is evolving with recent advances in molecular, functional, and hybrid/multiplex imaging. Integration of these various advanced diagnostic imaging and nonimaging methods aids understanding of cancer pathophysiology and provides a more comprehensive evaluation in this era of precision treatment. Here we discuss the current status of various advanced imaging techniques and their applications in head and neck cancer imaging.

Segmentation stability of human head and neck cancer medical images for radiotherapy applications under de-identification conditions: Benchmarking data sharing and artificial intelligence use-cases

Background

Demand for head and neck cancer (HNC) radiotherapy data in algorithmic development has prompted increased image dataset sharing. Medical images must comply with data protection requirements so that re-use is enabled without disclosing patient identifiers. Defacing, i.e., the removal of facial features from images, is often considered a reasonable compromise between data protection and re-usability for neuroimaging data. While defacing tools have been developed by the neuroimaging community, their acceptability for radiotherapy applications have not been explored. Therefore, this study systematically investigated the impact of available defacing algorithms on HNC organs at risk (OARs).

Methods

A publicly available dataset of magnetic resonance imaging scans for 55 HNC patients with eight segmented OARs (bilateral submandibular glands, parotid glands, level II neck lymph nodes, level III neck lymph nodes) was utilized. Eight publicly available defacing algorithms were investigated: afni_refacer, DeepDefacer, defacer, fsl_deface, mask_face, mri_deface, pydeface, and quickshear. Using a subset of scans where defacing succeeded (N=29), a 5-fold cross-validation 3D U-net based OAR auto-segmentation model was utilized to perform two main experiments: 1.) comparing original and defaced data for training when evaluated on original data; 2.) using original data for training and comparing the model evaluation on original and defaced data. Models were primarily assessed using the Dice similarity coefficient (DSC).

Results

Most defacing methods were unable to produce any usable images for evaluation, while mask_face, fsl_deface, and pydeface were unable to remove the face for 29%, 18%, and 24% of subjects, respectively. When using the original data for evaluation, the composite OAR DSC was statistically higher (p ≤ 0.05) for the model trained with the original data with a DSC of 0.760 compared to the mask_face, fsl_deface, and pydeface models with DSCs of 0.742, 0.736, and 0.449, respectively. Moreover, the model trained with original data had decreased performance (p ≤ 0.05) when evaluated on the defaced data with DSCs of 0.673, 0.693, and 0.406 for mask_face, fsl_deface, and pydeface, respectively.

Conclusion

Defacing algorithms may have a significant impact on HNC OAR auto-segmentation model training and testing. This work highlights the need for further development of HNC-specific image anonymization methods.

Multifactorial Model Based on DWI-Radiomics to Determine HPV Status in Oropharyngeal Squamous Cell Carcinoma

Background: Oropharyngeal squamous cell carcinoma (OPSCC) associated with human papillomavirus (HPV) has higher rates of locoregional control and a better prognosis than HPV-negative OPSCC. These differences are due to some unique biological characteristics that are also visible through advanced imaging modalities. We investigated the ability of a multifactorial model based on both clinical factors and diffusion-weighted imaging (DWI) to determine the HPV status in OPSCC. Methods: The apparent diffusion coefficient (ADC) and the perfusion-free tissue diffusion coefficient D were derived from DWI, both in the primary tumor (PT) and lymph node (LN). First- and second-order radiomic features were extracted from ADC and D maps. Different families of machine learning (ML) algorithms were trained on our dataset using five-fold cross-validation. Results: A cohort of 144 patients was evaluated retrospectively, which was divided into a training set (n = 95) and a validation set (n = 49). The 50th percentile of DPT, the inverse difference moment of ADCLN, smoke habits, and tumor subsite (tonsil versus base of the tongue) were the most relevant predictors. Conclusions: DWI-based radiomics, together with patient-related parameters, allowed us to obtain good diagnostic accuracies in differentiating HPV-positive from HPV-negative patients. A substantial decrease in predictive power was observed in the validation cohort, underscoring the need for further analyses on a larger sample size.

Evaluation of deep learning-based multiparametric MRI oropharyngeal primary tumor auto-segmentation and investigation of input channel effects: Results from a prospective imaging registry

BACKGROUND/PURPOSE

Oropharyngeal cancer (OPC) primary gross tumor volume (GTVp) segmentation is crucial for radiotherapy. Multiparametric MRI (mpMRI) is increasingly used for OPC adaptive radiotherapy but relies on manual segmentation. Therefore, we constructed mpMRI deep learning (DL) OPC GTVp auto-segmentation models and determined the impact of input channels on segmentation performance.

MATERIALS/METHODS

GTVp ground truth segmentations were manually generated for 30 OPC patients from a clinical trial. We evaluated five mpMRI input channels (T2, T1, ADC, Ktrans, Ve). 3D Residual U-net models were developed and assessed using leave-one-out cross-validation. A baseline T2 model was compared to mpMRI models (T2 + T1, T2 + ADC, T2 + Ktrans, T2 + Ve, all five channels [ALL]) primarily using the Dice similarity coefficient (DSC). False-negative DSC (FND), false-positive DSC, sensitivity, positive predictive value, surface DSC, Hausdorff distance (HD), 95% HD, and mean surface distance were also assessed. For the best model, ground truth and DL-generated segmentations were compared through a blinded Turing test using three physician observers.

RESULTS

Models yielded mean DSCs from 0.71 ± 0.12 (ALL) to 0.73 ± 0.12 (T2 + T1). Compared to the T2 model, performance was significantly improved for FND, sensitivity, surface DSC, HD, and 95% HD for the T2 + T1 model (p < 0.05) and for FND for the T2 + Ve and ALL models (p < 0.05). No model demonstrated significant correlations between tumor size and DSC (p > 0.05). Most models demonstrated significant correlations between tumor size and HD or Surface DSC (p < 0.05), except those that included ADC or Ve as input channels (p > 0.05). On average, there were no significant differences between ground truth and DL-generated segmentations for all observers (p > 0.05).

CONCLUSION

DL using mpMRI provides reasonably accurate segmentations of OPC GTVp that may be comparable to ground truth segmentations generated by clinical experts. Incorporating additional mpMRI channels may increase the performance of FND, sensitivity, surface DSC, HD, and 95% HD, and improve model robustness to tumor size.