Diffusion MRI correlation with p16 status and prediction for tumor progression in locally advanced head and neck cancer

Abstract

Purpose:

To investigate p16 effects on diffusion image metrics and associations with tumor progression in patients with locally advanced head and neck cancers.

Methods:

Diffusion images pretreatment and after 20 Gy (2wk) of RT were analyzed in patients with cT4/N3 p16+ oropharynx cancer (OPSCC) (N=51) and locoregionally advanced head and neck squamous cell carcinoma (LAHNSCC) (N=28), enrolled onto a prospective adaptive RT trial. Mean ADC values, subvolumes with ADC <1.2 um²/ms (TV_LADC), and peak values of low (µ_L) and high (µ_H) components of ADC histograms in primary and total nodal gross tumor volumes were analyzed for prediction of freedom from local, distant, or any progression (FFLP, FFDP or FFLRDP) using multivariate Cox proportional-hazards model with clinical factors. P value with false discovery control <0.05 was considered as significant.

Results:

With a mean follow up of 36 months, 18 of LAHNSCC patients and 16 of p16+ OPSCC patients had progression. After adjusting for p16, small µ_L and ADC values, and large TV_LADC of primary tumors pre-RT were significantly associated with superior FFLRDP, FFLP and FFDP in the LAHNSCC (p<0.05), but no diffusion metrics were significant in p16+ oropharynx cancers. Post ad hoc analysis of the p16+ OPSCC only showed that large TV_LADC of the total nodal burden pre-RT was significantly associated with inferior FFDP (p=0.05).

Conclusion:

ADC metrics were associated with different progression patterns in the LAHNSCC and p16+ OPSCC, possibly explained by differences in cancer biology and morphology. A deep understanding of ADC metrics is warranted to establish imaging biomarkers for adaptive RT in HNSCC.

Introduction

Diffusion magnetic resonance imaging (MRI) measures water mobility in the tissue environment with high sensitivity to microstructures of cells and cell membrane permeability. Apparent diffusion coefficient (ADC), a commonly used diffusion imaging parameter, has been shown to be prognostic and predictive for outcomes in head and neck squamous cell carcinomas (HNSCC) (1–7). Diffusion MRI does not require a gadolinium (Gd) based contrast agent and can be obtained within a few minutes. The simplicity of ADC from acquisition to computation as well as its predictive value results in advantages of ADC as an emerging imaging biomarker for stratifying progression risk in patients with HNSCC during adaptive (de-escalation or intensification) radiation therapy (RT).

Recent studies have shown that pre-treatment ADC values in p16- HNSCC tumors are greater than in p16+ ones, although with varying levels of significance (4, 8–10). It is known that p16+ oropharynx cancers have improved outcomes compared to other locally advanced head and neck cancers (LAHNSCC) (11–13). To date, whether the prognostic or predictive value of ADC is affected by tumor biology and morphology differences between p16- and p16+ HNSCC has not been assessed. Furthermore, it has been shown that the ADC distribution in HNSCC is deviated from a Gaussian distribution, which motivates studies of skewness and kurtosis of ADC distributions (5, 9). Furthermore, the ADC distribution in HNSCC, its response to RT, and its association with progression have not been characterized to account for p16 status differences.

In this study, we aimed to quantitatively characterize the distribution of ADC and its changes during RT in poor prognosis locally advanced HNSCCs, including cT4 or N3 p16+ oropharyngeal squamous cell carcinoma (OPSCC) and LAHNSCC (composed of p16- disease and p16+ non-oropharyngeal sites). We compared parameters of ADC distributions and their association with progression patterns between p16+ OPSCC and LAHNSCC tumors. This investigation could provide insight in selection of ADC metrics for prediction of progression risk for locally advanced HNSCC.

Methods

Patients

Patients with locally advanced HNSCC who were enrolled in a randomized phase II clinical trial between March 2014 and January 2020 were included in this analysis (NCT02031250). This trial was approved by the Institutional Review Board of the University of Michigan. Written consent was obtained from all enrolled patients. The clinical trial results are reported elsewhere (14). In brief, eligibility included patients with 1) cT4/N3 (AJCC 8 stage III) p16+ OPSCC, 2) locally advanced (T3-4/N2-3) p16- oropharyngeal or p16+ non-oropharyngeal head and neck cancer planned to undergo definitive chemoradiation therapy (CRT). p16 status was evaluated by immunohistochemistry. The patients were randomized to a standard arm of RT (70 Gy in 35 fractions) or an experimental RT boost arm, both with concurrent weekly cisplatin (40mg/m2) or carboplatin (AUC=2) for cisplatin ineligibility. In the experimental arm, a union of 1) the persisting low blood volume (BV) [BV<7.64 ml/100g based upon a previous histogram analysis (15)] pre-RT to after 20 Gy and 2) persisting low ADC [ADC< 1.2 um²/ms based upon a histogram analysis of previous works (16, 17)] pre-RT to after 20 Gy received 2.5 Gy per fraction for the last 15 of 35 fractions for a total dose of 80 Gy in 35 fractions. If the union of persisting subvolumes pre-RT to after 20 Gy was less than 1 cm³, the patient was entered into an observation arm and treated by standard RT (70 Gy in 35 fractions).

Diffusion imaging acquisition

Patients underwent MRI scans pre-treatment (< 2 weeks prior to the initiation of definitive CRT) and at fraction 10 (20 Gy) per protocol. All diffusion weighted (DW) images as well as T2- weighted and post-Gd T1-weighted images were acquired on a 3T scanner (Skyra, Siemens Healthineers). All patients were scanned in the treatment position using an individual-patient immobilization 5-point mask and bite block or Aquaplast mold as required for treatment. DW images were acquired by either a 2D spin-echo single shot echo-planar pulse sequence or a readout segmentation of long variable echo-trains (RESOLVE) pulse sequence that reduced geometric distortion (18) with spatial resolution of ~1.2×1.2×4.8 mm and b-values of 50 and 800 s/mm². ADC maps were calculated from the two b-value DW images to mitigate the perfusion effect by using in-house software that was technically validated in a QIN collaborative project (19). Quality and geometric alignment of ADC maps were assessed and reported previously (16).

Diffusion image metrics

Quantitative diffusion image metrics were calculated in the gross tumor volumes (GTVs) contoured manually on post-Gd T1-weighted images by the treating attending HN radiation oncologist and reviewed by the trial principal investigator (MM). Each tumor including primary and treated nodal tumors was contoured individually. Considering the dramatic reduction of gross head and neck movement during scanning by individual-patient immobilization devices, ADC maps were reformatted to match voxel-by-voxel of post-Gd T1-weighted images and overlaid with the GTV. Gross necrosis regions in the GTV were excluded by thresholding ADC below 2.7 um²/ms that was 10% below the value of free water diffusion.

The mean ADC value and the subvolume of low ADC thresholded at 1.2 um²/ms in each GTV (TV_LADC) were calculated. Also, as a bimodal distribution of ADC values in the primary GTV was observed, suggesting two major populations, the histogram of ADC was approximated by two Gaussian functions and fitted after binned with a size of 0.1 um²/ms using a Simplex optimizer written in C++ (see Supplementary Figure 1A). The low and high ADC components in the GTV, referred as respective L and H, were described by their peak ADC values (µ), widths (σ), and amplitudes (A). If a single Gaussian-like distribution was observed, the single peak ADC component was considered as 50% of each to be the low and high components.

Statistical analysis

The p16 effects on diffusion imaging metrics and association with tumor progression were tested using Kruskal-Wallis test and Cox proportional-hazards model. We considered freedom from local progression (FFLP), freedom from locoregional progression (FFLRP), freedom from distant progression (FFDP), and freedom from locoregional and distant progression (FFLRDP). The time to progression was defined from the starting date of RT to the date of local, regional or distant progression, and censored at other progressions that were not targeted in the test, death, or last follow-up. As approximately 95% of patients had T4/N3 diseases, p16 status and RT boost as clinical factors were considered in analysis. Smoking status as a controversial clinical factor for predicting specific patterns of progression was not included to limit overfitting (20–22). Multivariate Cox model was used to assess the image metrics one at a time with clinical factors for prediction of tumor progression. Considering multiple comparisons in the analysis, p values were corrected with false discovery rate control (FDC). The adjusted p value with FDC <0.05 was considered as significant. All analyses are summarized in Supplementary Figure 2.

Results

Patients and outcomes

We examined the imaging characteristics of 79 patients (median age of 64 years, 51 p16+ OPSCC and 8 females) randomized in an adaptive RT boost trial, 40 on the standard arm and 39 on the experimental arm. The patient characteristics are provided in Table 1. The details of the trial and outcomes were provided elsewhere (14). In brief, the mean follow-up was 36 months, median 30 months (range 8-83 months) for patients without death, with minimum of 12 months follow-up in all patients without disease progression except one who was lost in follow-up. LAHNSCC disease included 50% p16- oropharynx, 11% larynx, 18% hypopharynx, 14% sinonasal, and 7% EBV negative nasopharynx, with known similar outcomes (23, 24). In the 28 LAHNSCC patients, 10 had no evidence of disease (NED), 11 patients had local failure (LF) in which 4 had regional failure (RF) and 7 had distant failure (DF) at the same time, 6 had DF only, and 1 had both RF and DF. Also, 9 of the 28 received 80 Gy RT adaptive boost. In 51 of the p16+ oropharynx patients, 34 had NED in which 13 received 70 Gy standard treatment (NED₇₀) and 21 had 80 Gy RT boost (NED₈₀), 7 had LF in which 1 had RF at the same time, 9 had DF in which 1 had RF at the same time, and 1 had RF only. There was no effect of RT boost on FFLRP in the LAHNSCC patients and an observable boost effect in the p16+ oropharynx patients, see Supplementary Figure 1B, suggesting that the effect of the boosting dose may depend upon p16 status.

Characteristics of imaging metrics of LAHNSCC and p16+ oropharynx cancers

Considering differences of tumor biology, outcome, and response to radiation boosting between the p16+ OPSCC and the LAHNSCC, diffusion metrics (namely TV_LADC, mean ADC and peak ADC values of low (µ_L) and high (µ_H) components) were characterized first.

In the LAHNSCC, the patients with NED had the smallest values of mean ADC values in both primary tumor and total nodal tumor volumes pre-RT and at 2wk, of µ_L and µ_H in primary tumors pre-RT and at 2wk, and of primary and total nodal GTVs pre-RT and at 2wk, compared to the patients with LF (with or without other progression) or DF only, but the differences between the three subgroups were not significant with FDC (p>0.1) (Supplementary Table 1). Averaged histograms of pre-RT ADC distributions of primary GTVs associated with NED, LF and DF only pre-RT and 2wk are plotted in Figure 1.

Figure 1

In the p16+ OPSCC, RT boost effect was observed in FFLRP (Supplementary Figure 1B), but also in diffusion metrics (Supplementary Table 2). Compared with LF and DF subgroups, the patients with NED₇₀ had the greatest mean ADC and the smallest TV_LADC in primary tumors pre-RT and at 2wk (p<0.05). The µ_L and µ_H values of primary tumors followed the same trend as the mean ADC. Compared with the NED₇₀ patients, the NED₈₀ patients had significantly low values of mean ADC of primary tumors pre-RT and at 2wk and µ_H pre-RT (0.004, 0.04, and 0.02, respectively), and had significantly large TV_LADC of primary tumor pre-RT and at 2wk (p<0.01, and 0.006, respectively), see Figure 2. This suggests that RT boost may overcome higher tumor cellularity or density. Figure 1 shows the averaged histograms of ADC distributions in primary GTVs of the subgroups pre-RT and 2wk.

Figure 2

Predictive values of ADC metrics and GTVs for progression

Considering the observed opposite trends in ADC metrics between the p16+ OPSCC and the LAHNSCC, an interaction effect between the diffusion metric and p16 status was suggested. RT boost effect was observed only in the p16+ OPSCC but not in the LAHNSCC, suggesting that RT boost effect interacted with p16 status. To avoid overfitting due to too many co-variables in the multivariate Cox model (e.g., boosting effect, the interactions between boost and the diffusion metrics as well as between boost and p16 status) and to best model the data applied to standard clinical practice, the 21 patients with NED₈₀ were excluded from the progression prediction models but 2 patients with LF and with boost were included. All patients with DF were included since there was no boost effect expected for FFDP. In the multivariate Cox models of progression prediction, p16 status, a diffusion metric (TV_LADC, mean ADC, µ_L, or µ_H) and the interaction of the diffusion metric and p16 status were considered.

After adjusting for p16 effect, Cox models for prediction of FFLP found significant effects of pre-RT TV_LADC, pre-RT mean ADC, and pre-RT µ_L of primary tumor for the LAHNSCC (p<0.04 with FDC), and effects of the interactions of pre-RT TV_LADC, pre-RT mean ADC, 2wk µ_L of primary tumor with p16 status (p<0.03 without FDC), but no significant effect of TV_LADC, mean ADC, µ_L or µ_H of primary tumor pre-RT or at 2wk for p16+ OPSCC (Table 2). After adjusting for p16 effect, Cox models for prediction of FFDP found no significant effects of the tested diffusion metrics of p16+ OPSCC. In LAHNSCC, the significant effects of pre-RT TV_LADC, pre-RT mean ADC, pre-RT µ_L and pre-RT µ_H of primary tumor for the LAHNSCC (p<0.04 with FDC), and an effect of the interaction of pre-RT TV_LADC of nodal tumor with p16 status (p<0.04 without FDC) were noted (Table 2). After adjusting for p16 effect, Cox models for prediction of FFLRDP found significant effects of pre-RT TV_LADC, pre-RT mean ADC, pre-RT µ_L and of pre-RT µ_H of primary tumor for the LAHNSCC (p<0.05 with FDC), and significant effects of the interactions of pre-RT TV_LADC, and pre-RT µ_L of primary tumor with p16 status (p<0.02 with FDC), but no significant effect of the tested diffusion metrics of primary or nodal tumor pre-RT or at 2wk for p16+ oropharynx cancers (Table 3). After adjusting for p16 status, primary GTV at 2wk predicted significantly for FFLRDF (p<0.03, HR=2.21, GTV > the median value (48.0 cm³)) but no primary or nodal GTV pre-RT or 2wk predicted for FFLP or FFDP.

Post ad hoc analysis confined to patients with p16+ OPSCC showed that only TV_LADC of the total nodal burden pre-RT was a significant predictor for FFDP (p=0.05, HR=8.02(0.99-65.2) for nodal TV_LADC >2.3 cm³), and primary GTVs pre-RT and at 2wk were significantly predictors for FFLRDF (p<0.02, HR=3.52(1.21-10.2) for GTV pre-RT > 57.1 cm³ and p<0.02, HR=3.90(1.30-11.7) GTV 2wk >48.0 cm³). Kaplan-Meier curves of FFLP, FFDP and FFLDP in the p16+ OPSCC and LAHNSCC are shown in Figure 3.

Figure 3

Discussion

In this study, we assessed quantitative diffusion metrics and ADC histograms for their associations with tumor progression in the patients with locally advanced p16+ OPSCC and LAHNSCC (95% of T4/N3) and enrolled on a randomized phase II trial of adaptive RT boost. With expected outcome differences between the advanced p16+ OPSCC and LAHNSCC, we found different ADC characteristics and associations with tumor progression between the two. Particularly, for the LAHNSCC, low ADC and large TV_LADC were associated with low risk of local and distant tumor progression, which could be interpreted as a lower ADC value associated with a less extent of stroma in tumor. However, for the p16+ OPSCC, high ADC, and small TV_LADC were associated with a trend of low risks of local and distant tumor progression, which may be attributed to tumor infiltrating lymphocytes. In analysis confined to p16+ OPSCC, we did see that nodal TV_LADC may be an imaging marker for distant progression, and primary GTV pre-RT and at 2 weeks during RT seem to be a stronger predictor for local, regional and distant progression albeit with the limitations of ad hoc analysis. As all, the ADC value appears to be affected substantially by biology and morphology of p16- and p16+ tumors as well as associated heterogeneity. Diffusion images, although easily acquired, could be affected by many biological, clinical and physical factors. Further investigations of these factor effects on diffusion images through pathologic correlation are needed to guide either radiation de-escalation or treatment intensification trials using ADC metrics.

ADC is sensitive to tumor microstructure. Histologically, p16- HNSCC is typically comprised of a keratinizing morphology with angulated nests of tumor cells, abundant cytoplasm, stromal desmoplasia, and central necrosis; p16+ oropharynx cancer is typically described as a nonkeratinizing morphology with the presence of a large amount of tumor infiltrating lymphocytes as well as small amounts of cytoplasm, central necrosis and stroma (25). These differing morphological features affect the ADC value and the ADC distribution in these tumors. For instance, it has been reported that the ADC value is positively correlated to the total percentage area of stroma and inversely correlated to the cell density in the HN tumors (26). When the tumor morphology variation manifests in the macroscopic level, i.e., in the millimeter range, an effect can be observed in the ADC value and distribution. It is plausible that in the p16- tumors, a low mean ADC value or low peak ADC values of the two components indicate a low total percentage of stroma, and thereby less protection of tumor cells from CRT by stroma and better tumor control (27–32). In the p16+ oropharynx tumors, tumor microstructure may be affected by tumor infiltrating lymphocytes (33, 34), where a large subvolume with low ADC in the tumor volume, low mean ADC or low µ_L value, could need to be treated with high radiation doses to have an improved local and regional tumor control. ADC is measured at the macroscopic level but affected by very different microscopic morphology and biology. All these challenge the analysis and interpretation of ADC and ADC changes in HNSCC.

It is not entirely clear how p16 status affects tumor ADC distribution changes in response to CRT. Although after 20 Gy of radiation both types of tumors show an increase in mean ADC values, the LAHNSCC tumors maintained the bi-distribution in a certain extent while the p16+ oropharynx tumors show a rapid normalization of the ADC distributions. Note that the histograms in Figure 1 are from the subgroup averages, which further remove individual variations. Nevertheless, a mean ADC in the tumor could over-simplify the heterogeneity of ADC distribution in HNSCC, particularly for the locally advanced tumors, and could obscure the different responses between p16- and p16+ tumors. A histogram analysis of the ADC bi-distribution or the low ADC component (e.g., subvolumes of low ADC defined by a threshold) could remove confounding effects to an extent and thereby increase the predictive power for progression. While ADC metrics in HNSCC are useful to differentiate tumor control vs progression, primary GTV in the p16+ oropharynx tumors seems to have sufficient power for prediction of local, regional and distant progression, but not in other LAHNSCC tumors.

This study has a few limitations. The number of patients in this study is still small. In line with standard clinical practice, our LAHNSCC patients included p16- oropharynx cancer as well as both p16+ and p16- non-oropharyngeal cancers. The potential biologic or clinical significance of p16 status outside the oropharynx is unclear with some suggestion that p16 positivity outside the oropharynx is less correlated with HPV positivity and may be driven by other molecular mechanisms (12, 13). However, the tumor ADC distribution could be different at different sites. The 21 p16+ oropharynx patients who received 80 Gy local RT boost and had no progression had large GTVs and ADC metrics for high risk of tumor progression. To avoid overfitting, the 21 patients were excluded from the progression analysis. Nevertheless, characterizing the ADC metrics in these two groups of patients according to current clinical practice further reveals differences between the two types of HNSCC, and has the potential to improve the power of the ADC metrics as a biomarker for assessment of tumor response and prediction of progression in HNSCC. Considering complexity of tumor biology and response to chemoradiation therapy, integrating multi-imaging biomarkers, including FDG PET and dynamic contrast enhanced MRI, as well as liquid biomarkers (e.g., circulating tumor human papilloma virus DNA) during the early course of CRT could improve prediction of tumor progression, which could provide sufficient time and guidance for individualized (intensified or de-intensified) adaptation of CRT and thereby improve outcomes of patients with locally advanced head and neck cancers.

References

• 1

  HatakenakaMNakamuraKYabuuchiHShioyamaYMatsuoYOhnishiKet al. Pretreatment apparent diffusion coefficient of the primary lesion correlates with local failure in head-and-neck cancer treated with chemoradiotherapy or radiotherapy. Int J Radiat Oncol Biol Phys (2011) 81(2):339–45. doi: 10.1016/j.ijrobp.2010.05.051

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  KimSLoevnerLQuonHShermanEWeinsteinGKilgerAet al. Diffusion-weighted magnetic resonance imaging for predicting and detecting early response to chemoradiation therapy of squamous cell carcinomas of the head and neck. Clin Cancer Res (2009) 15(3):986–94. doi: 10.1158/1078-0432.CCR-08-1287

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  LambrechtMVan CalsterBVandecaveyeVDe KeyzerFRoebbenIHermansRet al. Integrating pretreatment diffusion weighted MRI into a multivariable prognostic model for head and neck squamous cell carcinoma. Radiother Oncol (2014) 110(3):429–34. doi: 10.1016/j.radonc.2014.01.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  CaoYAryalMLiPLeeCSchipperMHawkinsPGet al. Predictive Values of MRI and PET Derived Quantitative Parameters for Patterns of Failure in Both p16+ and p16- High Risk Head and Neck Cancer. Front Oncol (2019) 9:1118. doi: 10.3389/fonc.2019.01118

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  KingADChowKKYuKHMoFKYeungDKYuanJet al. Head and neck squamous cell carcinoma: diagnostic performance of diffusion-weighted MR imaging for the prediction of treatment response. Radiology (2013) 266(2):531–8. doi: 10.1148/radiol.12120167

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  Xiao-pingYJingHFei-pingLYinHQiangLLanlanWet al. Intravoxel incoherent motion MRI for predicting early response to induction chemotherapy and chemoradiotherapy in patients with nasopharyngeal carcinoma. J Magn Reson Imaging (2016) 43(5):1179–90. doi: 10.1002/jmri.25075

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  HongJYaoYZhangYTangTZhangHBaoDet al. Value of magnetic resonance diffusion-weighted imaging for the prediction of radiosensitivity in nasopharyngeal carcinoma. Otolaryngol Head Neck Surg (2013) 149(5):707–13. doi: 10.1177/0194599813496537

  • CrossRef
  • Google Scholar

• 8

  DriessenJPvan BemmelAJvan KempenPMJanssenLMTerhaardCHPameijerFAet al. Correlation of human papillomavirus status with apparent diffusion coefficient of diffusion-weighted MRI in head and neck squamous cell carcinomas. Head Neck (2016) 38 Suppl 1:E613–8. doi: 10.1002/hed.24051

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  de PerrotTLenoirVDomingo AyllonMDulguerovNPusztaszeriMBeckerMet al. Apparent diffusion coefficient histograms of human papillomavirus-positive and human papillomavirus-negative head and neck squamous cell carcinoma: assessment of tumor heterogeneity and comparison with histopathology. AJNR Am J Neuroradiol (2017) 38(11):2153–60. doi: 10.3174/ajnr.A5370

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  SchoutenCSde GraafPBloemenaEWitteBIBraakhuisBJBrakenhoffRHet al. Quantitative diffusion-weighted MRI parameters and human papillomavirus status in oropharyngeal squamous cell carcinoma. AJNR Am J Neuroradiol (2015) 36(4):763–7. doi: 10.3174/ajnr.A4271

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  AngKKWheelerRWeberRRosenthalDINguyen-TanPFet al. Human papillomavirus and survival of patients with oropharyngeal cancer. N Engl J Med (2010) 363(1):24–35. doi: 10.1056/NEJMoa0912217

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FakhryCZhangQGillisonMLNguyen-TanPFRosenthalDIWeberRSet al. Validation of NRG oncology/RTOG-0129 risk groups for HPV-positive and HPV-negative oropharyngeal squamous cell cancer: Implications for risk-based therapeutic intensity trials. Cancer (2019) 125(12):2027–38. doi: 10.1002/cncr.32025

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  ChungCHZhangQKongCSHarrisJFertigEJHarariPMet al. p16 protein expression and human papillomavirus status as prognostic biomarkers of nonoropharyngeal head and neck squamous cell carcinoma. J Clin Oncol (2014) 32(35):3930–8. doi: 10.1200/JCO.2013.54.5228

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  MierzwaMLAryalMLeeCSchipperMVanTilMMoralesKet al. Randomized phase II study of physiologic mri-directed adaptive radiation boost in poor prognosis head and neck cancer. Clin Cancer Res (2022) CCR-22-1522. doi: 10.1158/1078-0432.CCR-22-1522

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  WangPPopovtzerAEisbruchACaoY. An approach to identify, from DCE MRI, significant subvolumes of tumors related to outcomes in advanced head-and-neck cancer. Med Phys (2012) 39(8):5277–85. doi: 10.1118/1.4737022

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  TengFAryalMLeeJLeeCShenXHawkinsPGet al. Adaptive boost target definition in high-risk head and neck cancer based on multi-imaging risk biomarkers. Int J Radiat Oncol Biol Phys (2018) 102(4):969–77. doi: 10.1016/j.ijrobp.2017.12.269

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HouwelingACWolfALVogelWVHamming-VriezeOvan Vliet-VroegindeweijCvan de KamerJBet al. FDG-PET and diffusion-weighted MRI in head-and-neck cancer patients: implications for dose painting. Radiother Oncol (2013) 106(2):250–4. doi: 10.1016/j.radonc.2013.01.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  ZhaoMLiuZShaYWangSYeXPanYet al. Readout-segmented echo-planar imaging in the evaluation of sinonasal lesions: A comprehensive comparison of image quality in single-shot echo-planar imaging. Magn Reson Imaging (2016) 34(2):166–72. doi: 10.1016/j.mri.2015.10.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  NewittDCMalyarenkoDChenevertTLQuarlesCCBellLFedorovAet al. Multisite concordance of apparent diffusion coefficient measurements across the NCI Quantitative Imaging Network. J Med Imaging (Bellingham) (2018) 5(1):011003. doi: 10.1117/1.JMI.5.1.011003

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  RosenBSWilkieJRSunYIbrahimMCasperKAMillerJEet al. CT and FDG-PET radiologic biomarkers in p16+ oropharyngeal squamous cell carcinoma patients treated with definitive chemoradiotherapy. Radiother Oncol (2021) 155:174–81. doi: 10.1016/j.radonc.2020.10.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HuangSHXuWWaldronJSiuLShenXTongLet al. Refining American Joint Committee on Cancer/Union for International Cancer Control TNM stage and prognostic groups for human papillomavirus-related oropharyngeal carcinomas. J Clin Oncol (2015) 33(8):836–45. doi: 10.1200/JCO.2014.58.6412

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  OlivaMHuangSHXuWSuJHansenARBratmanSVet al. Impact of cisplatin dose and smoking pack-years in human papillomavirus-positive oropharyngeal squamous cell carcinoma treated with chemoradiotherapy. Eur J Cancer (2019) 118:112–20. doi: 10.1016/j.ejca.2019.06.019

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  StenmarkMHMcHughJBSchipperMWallineHMKomarckCFengFYet al. Nonendemic HPV-positive nasopharyngeal carcinoma: association with poor prognosis. Int J Radiat Oncol Biol Phys (2014) 88(3):580–8. doi: 10.1016/j.ijrobp.2013.11.246

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  GallowayTJZhangQENguyen-TanPFRosenthalDISoulieresDFortinAet al. Prognostic value of p16 status on the development of a complete response in involved oropharynx cancer neck nodes after cisplatin-based chemoradiation: A secondary analysis of NRG oncology RTOG 0129. Int J Radiat Oncol Biol Phys (2016) 96(2):362–71. doi: 10.1016/j.ijrobp.2016.05.026

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  WestraWH. The morphologic profile of HPV-related head and neck squamous carcinoma: implications for diagnosis, prognosis, and clinical management. Head Neck Pathol (2012) 6 Suppl 1:S48–54. doi: 10.1007/s12105-012-0371-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  DriessenJPCaldas-MagalhaesJJanssenLMPameijerFAKooijNTerhaardCHet al. Diffusion-weighted MR imaging in laryngeal and hypopharyngeal carcinoma: association between apparent diffusion coefficient and histologic findings. Radiology (2014) 272(2):456–63. doi: 10.1148/radiol.14131173

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  NeesseAMichlPFreseKKFeigCCookNJacobetzMAet al. Stromal biology and therapy in pancreatic cancer. Gut (2011) 60(6):861–8. doi: 10.1136/gut.2010.226092

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  MuellerMMFusenigNE. Friends or foes - bipolar effects of the tumour stroma in cancer. Nat Rev Cancer (2004) 4(11):839–49. doi: 10.1038/nrc1477

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  WangKMaWWangJYuLZhangXWangZet al. Tumor-stroma ratio is an independent predictor for survival in esophageal squamous cell carcinoma. J Thorac Oncol (2012) 7(9):1457–61. doi: 10.1097/JTO.0b013e318260dfe8

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  de KruijfEMvan NesJGvan de VeldeCJPutterHSmitVTLiefersGJet al. Tumor-stroma ratio in the primary tumor is a prognostic factor in early breast cancer patients, especially in triple-negative carcinoma patients. Breast Cancer Res Treat (2011) 125(3):687–96. doi: 10.1007/s10549-010-0855-6

  • CrossRef
  • Google Scholar

• 31

  HuijbersATollenaarRAPeltGWZeestratenECDuttonSMcConkeyCCet al. The proportion of tumor-stroma as a strong prognosticator for stage II and III colon cancer patients: validation in the VICTOR trial. Ann Oncol (2013) 24(1):179–85. doi: 10.1093/annonc/mds246

  • CrossRef
  • Google Scholar

• 32

  WisemanBSWerbZ. Stromal effects on mammary gland development and breast cancer. Science (2002) 296(5570):1046–9. doi: 10.1126/science.1067431

  • CrossRef
  • Google Scholar

• 33

  WardMJThirdboroughSMMellowsTRileyCHarrisSSuchakKet al. Tumour-infiltrating lymphocytes predict for outcome in HPV-positive oropharyngeal cancer. Br J Cancer (2014) 110(2):489–500. doi: 10.1038/bjc.2013.639

  • CrossRef
  • Google Scholar

• 34

  ChakravarthyAHendersonSThirdboroughSMOttensmeierCHSuXLechnerMet al. Human papillomavirus drives tumor development throughout the head and neck: improved prognosis is associated with an immune response largely restricted to the oropharynx. J Clin Oncol (2016) 34(34):4132–41. doi: 10.1200/JCO.2016.68.2955

  • CrossRef
  • Google Scholar