Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 15 February 2022
Sec. Head and Neck Cancer

Radiomic Analysis of Tumour Heterogeneity Using MRI in Head and Neck Cancer Following Chemoradiotherapy: A Feasibility Study

Amrita Guha,,Amrita Guha1,2,3Mustafa AnjariMustafa Anjari4Gary Cook,Gary Cook3,5Vicky Goh,Vicky Goh3,4Steve Connor,,*Steve Connor3,4,6*
  • 1Department of Radio-Diagnosis, Tata Memorial Hospital, Mumbai, India
  • 2Training School Complex, Homi Bhabha National Institute, Mumbai, India
  • 3School of Biomedical Engineering & Imaging Sciences, Faculty of Life Sciences & Medicine, King’s College London, London, United Kingdom
  • 4Department of Radiology, Guy’s and St Thomas’ Hospital, London, United Kingdom
  • 5King’s College London & Guy’s and St Thomas’ Positron Emission Tomography (PET) Centre, London, United Kingdom
  • 6Department of Neuroradiology, King’s College Hospital, London, United Kingdom

Objectives: To evaluate interval changes in heterogeneity on diffusion-weighted apparent diffusion coefficient (ADC) maps and T1-weighted post-gadolinium (T1w post gad) MRI in head and neck carcinoma (HNSCC), with and without chemo-radiotherapy (CRT) response.

Methods: This prospective observational cohort study included 24 participants (20 men, age 62.9 ± 8.8 years) with stage III and IV HNSCC. The primary tumour (n = 23) and largest lymph node (n = 22) dimensions, histogram parameters and grey-level co-occurrence matrix (GLCM) parameters were measured on ADC maps and T1w post gad sequences, performed pretreatment and 6 and 12 weeks post CRT. The 2-year treatment response at primary and nodal sites was recorded. The Wilcoxon signed-rank test was used to compare interval changes in parameters after stratifying for treatment response and failure (p < 0.001 statistical significance).

Results: 23/23 primary tumours and 18/22 nodes responded to CRT at 2 years. Responding HNSCC demonstrated a significant interval change in ADC histogram parameters (kurtosis, coefficient of variation, entropy, energy for primary tumour; kurtosis for nodes) and T1w post gad GLCM (entropy and contrast in the primary tumour and nodes) by 6 weeks post CRT (p < 0.001). Lymph nodes with treatment failure did not demonstrate an interval alteration in heterogeneity parameters.

Conclusions: ADC maps and T1w post gad MRI demonstrate the evolution of heterogeneity parameters in successfully treated HNSCC by 6 weeks post CRT; however, this is not observed in lymph nodes failing treatment.

Advances in Knowledge: Early reduction in heterogeneity is demonstrated on MRI when HNSCC responds to CRT.

Introduction

Head and neck squamous cell carcinoma (HNSCC) is the seventh most common cancer worldwide (1). Chemo-radiotherapy (CRT) provides the best opportunity for cure in advanced-stage HNSCC (2); however, tumour resistance or insufficient therapy may result in treatment failure in more than 30% of patients (3). Whilst earlier detection of residual viable tumour allows for salvage surgery and improved survival (4), it is currently challenging to evaluate this with clinical examination and cross-sectional imaging, due to the presence of posttreatment tissue distortion (5). Metabolic imaging with 18F-fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) (6) may overcome some of the difficulties in interpretation with conventional CT and MRI in this clinical context, but it is generally delayed for at least 12 weeks, due to earlier false positives from post-CRT inflammation (7).

Assessing changes in MRI signal heterogeneity within tumours (8, 9) may better reflect residual disease. There are limited studies addressing the application of such analysis of imaging features in the prediction of treatment response (10, 11), and there are no data on their early posttreatment evolution in HNSCC. It is well accepted that HNSCCs display a marked “heterogeneity” on histology with variation in proliferation and cellular differentiation within different regions of the tumours (12, 13). We hypothesised that alterations in tumour heterogeneity following CRT would be reflected by diffusion-weighted and post-gadolinium T1w MRI and that measurements of signal heterogeneity may augment standard response assessment based on size criteria alone.

Therefore, our primary objective was to evaluate the interval changes in signal heterogeneity on diffusion-weighted (DWI) apparent diffusion coefficient (ADC) maps and T1-weighted post-gadolinium (T1w post gad) MRI within the primary tumour and largest metastatic lymph node at 6 and 12 weeks following CRT. Our secondary objectives were to evaluate how interval changes vary according to treatment response and to compare the interval changes in parameters to those of conventional size criteria.

Materials and Methods

Study Design

Participants were recruited for a prospective single-centre cohort observational study (ISRCTN58327080; Research Ethics Committee approval 13/LO/1876) and provided informed consent for participation.

Participants

Participants were eligible if there was histologically confirmed stage III or IV primary squamous cell carcinoma of the head and neck (HNSCC) without distant metastatic disease, a 1-cm2 area of measurable primary tumour and/or metastatic locoregional node on the basis of standard clinico-radiological staging, and where curative CRT was planned. Exclusion criteria were prior chemo- or radiotherapy; ECOG performance status >2; lack of capacity to provide informed consent; and known contrast agent allergy or renal impairment.

Treatment

Intensity-modulated radiotherapy consisting of 70 Gy in 35 fractions was delivered, 2 Gy per fraction once daily, 5 days a week. Concomitant intravenous cisplatin at a dose of 35 mg/m2 every 7 days, starting on day 1 of radiotherapy, was administered to all participants.

MRI Imaging

Participants underwent MRI before treatment and at 6 and 12 weeks after completion of chemoradiotherapy. Imaging was performed on a 1.5-Tesla MRI system (MAGNETOM Aera, Siemens Healthcare, Erlangen, Germany) using a 20-channel phased-array surface neck coil and included T2-weighted, diffusion-weighted imaging (DWI) and post gad T1w sequences. The MRI acquisition protocol is summarised in Supplemental Table 1.

MRI Processing and Analysis

Two sets of freehand regions of interest (ROI) were delineated by a radiologist (AG, 7 years of experience) on both DWI (b = 800) and post gad T1w images with reference to the other sequences. ROIs were placed at the site of the measurable primary tumour and/or largest pathological lymph node on the pretreatment, 6-week and 12-week post-treatment MRIs. ROIs were placed on multiple sections to encompass the whole lesion depicted as high signal on the DWI b-800 or high signal on post gad T1-w sequences (Figure 1). Internal areas of necrosis (non-enhancement and high diffusion signal, respectively) were included, but areas of peri-tumoural inflammation (adjacent avid gadolinium enhancement and high signal on T2w images) were excluded. When a focus of increased signal relative to background was not evident on the post-CRT imaging, this was defined as non-measurable, and a representative ROI (6-mm diameter/28 mm2) was placed at the site of the pretreatment ROI. For DWI, the ROI was transferred onto the ADC maps calculated from the b = 100 and b = 800 values to derive ADC parameters. In patients where it was challenging to draw ROIs on the posttreatment scans, the T2W and T1 post-gadolinium images were also correlated with ensuring the accuracy of delineation.

FIGURE 1
www.frontiersin.org

Figure 1 Workflow of the study: ROIs were placed on the post gad T1w sequences: primary tumour (A), single largest pathological lymph node (B) and on the ADC map: primary tumour (C). The radiomics first- and second-order parameters were derived and are pictorially depicted here by a sample of the shaded surface intensity histogram (D). These analyses were then repeated at the 6- and 12-week MRI studies.

Primary tumour and nodal histogram and grey-level co-occurrence matrix (GLCM) parameters were extracted from the ROIs using an in-house radiomics software platform, validated through the International Biomarker Standardization Initiative (14). Analysis was restricted to these parameters to reflect the limitations of the 2D acquisition and in-plane voxel resolution. For example, GLCM analysis was only applied to post gad T1w images since the larger voxel size (0.9 × 0.9 × 5 mm) of the diffusion acquisition resulted in insufficient voxels for analysis.

The mean and median gadolinium-enhanced T1w signal and ADC value were also recorded as well as the lesion long axis length and volume on the post gad T1w images.

Clinical Data and Treatment Outcome

Patient demographics, TN status, human papillomavirus (HPV) status, tumour site and subsite tumour stage were recorded. Follow-up included a standard-of-care 12-week 18F-FDG PET/CT study which was used to guide management and clinical assessment at 1 year and 2 years following completion of chemo-radiotherapy. The response of primary tumour and largest lymph node at 2 years was determined by cytological or histological confirmation (biopsy or resection) or by serial progression on imaging follow-up.

Statistical Analyses

SPSS statistical software Version 21.0 was used for analysis (IBM SPSS Statistics for Windows, Version 21.0. Armonk, NY: IBM Corp. Released 2012).

Pearson’s correlation was performed as a method of parameter reduction amongst the initial 21 GLCM parameters, excluding those with a high correlation coefficient (>0.75). This left 4 GLCM parameters for further analysis (Supplemental Table 2).

Descriptive data were collated for pretreatment, 6-week and 12-week post-CRT histogram and the selected GLCM parameters. This was performed separately for primary tumour and lymph node values and according to whether there was treatment response at 2 years.

As variables were not normally distributed, the Wilcoxon signed-rank test was performed to assess the interval changes in heterogeneity parameters. A p-value of p < 0.001 was considered statistically significant to account for multiple testing.

Results

Participants

The participant consort flow diagram is demonstrated in Figure 2. Of the 24 participants (20 men, 4 women, mean age 62.9 +_8.79 years), there were 5 (20%) patients with stage III disease and 19 (79%) patients with stage IV disease. Participant characteristics including primary site, HPV status and TN staging are summarised in Table 1. There were 23 primary tumours and 22 lymph nodes analysed.

FIGURE 2
www.frontiersin.org

Figure 2 Patient flowchart.

TABLE 1
www.frontiersin.org

Table 1 Participant characteristics of the HNSCC cohort.

At 6 weeks, there were 19/23 non-measurable primary tumours and 10/22 non-measurable nodes whilst at 12 weeks there were 23/23 non-measurable primary tumours and 18/22 non-measurable nodes. All primary tumours and 18/22 lymph nodes responded to CRT at 2 years. The 4 cases of lymph node recurrence corresponded to the same largest lymph node site as that undergoing analysis. Patients with non-responding lymph nodes underwent salvage neck dissection following initial cytological assessment.

Interval Changes Following CRT in the Primary Tumour

All primary tumours responded to CRT at 2 years following treatment.

The primary tumour showed a significant mean decrease of 71% (6.3+_0.8 cm) in diameter and 98% (1.45+_0.5 cm3) in volume at 6 weeks and a further reduction between 6 and 12 weeks (p <.0001).

There was a significant increase in mean values of ADCmean and ADCmedian (729 to 1,019 × 10-6 mm2/s; 40% and 692 to 1,062 × 10-6 mm2/s; 53%) by 6 weeks post CRT, but there was no further significant increase between 6 and 12 weeks. Kurtosis (-42%), coefficient of variation (CoV) (-46%) and entropy (-44%) significantly decreased, and energy significantly (+238%) increased by 6 weeks of treatment (p <.0001), again without further significant alterations between 6 and 12 weeks post CRT (Table 2).

TABLE 2
www.frontiersin.org

Table 2 Pretreatment, 6-week post-CRT and 12-week post-CRT primary tumour ADC and post-gadolinium T1w histogram parameters and their interval changes with treatment (all with treatment success n = 23).

None of the histogram parameters showed any significant changes on post gad T1w images at 6 or 12 weeks. GLCM contrast significantly increased (+106%) and GLCM entropy (-36%) significantly decreased by 6 weeks, with no further significant change between 6 and 12 weeks post CRT (Table 2).

Interval Changes Following CRT in the Lymph Nodes

There were 18/22 lymph nodes that responded to CRT at 2 years following treatment.

The responding nodes showed a significant decrease in the mean diameter of 68% (1.4+_0.3) and volume of 95% (8+_1.4) at 6 weeks and a further reduction between 6 and 12 weeks (p <.0001).

There was no significant increase in mean values of ADCmean and ADCmedian by 6 or 12 weeks post CRT.

ADC kurtosis showed a significant mean decrease (-35%) at 6 weeks. By 12 weeks, energy (+170%) and CoV (+280%) showed a significant mean increase compared to the pretreatment values whilst there was a statistically significant mean reduction in entropy (-39%).

None of the post gad T1w histogram parameters showed significant change with treatment. However, mean GLCM contrast significantly increased and mean GLCM entropy significantly decreased (+61% and -20%) by 6 weeks post CRT with no further significant alteration between 6 and 12 weeks (Table 3).

TABLE 3
www.frontiersin.org

Table 3 Pretreatment, 6-week post-CRT and 12-week post-CRT ADC and post-gadolinium T1w texture parameters and their interval changes in metastatic lymph nodes of participants with treatment response (n = 18).

The non-responding metastatic lymph nodes did not demonstrate a significant change in any ADC or post gadolinium T1w histogram or GLCM parameters between pretreatment and post-CRT studies. They also did not show a significant decrease in diameter or volume (Table 4).

TABLE 4
www.frontiersin.org

Table 4 Pretreatment, 6-week post-CRT and 12-week post-CRT ADC and post-gadolinium T1w histogram parameters and their interval changes in metastatic lymph nodes of participants with treatment failure (n = 4).

Discussion

Summary of Findings

Interval changes in signal heterogeneity have been evaluated to provide insight into the effect of CRT on the primary tumour and lymph nodes in HNSCC, and to identify whether heterogeneity changes might provide an earlier indication of treatment response. Changes in selected histogram parameters were observed by 6 weeks on ADC maps in responding primary tumours (kurtosis, CoV, entropy, energy) and responding lymph nodes (kurtosis). Whilst there were no post CRT changes in histogram parameters on post gad T1w images, there were significant decreases in GLCM entropy for both successfully treated primary tumours and lymph nodes. The non-responding lymph nodes did not demonstrate any of these interval changes in heterogeneity. Interval reduction in size was also demonstrated by 6 weeks in responding lymph nodes but not in non-responding lymph nodes.

ADC

Previous studies have evaluated the role of diffusion-weighted MRI in predicting HNSCC treatment response post CRT, and most have found that a greater rise in ADCmean from the pretreatment to the intra-treatment (1519) or posttreatment MRI (20, 21) predicts treatment outcomes. Similarly, our study demonstrated that primary tumour ADCmean increased by 6 weeks post CRT. However, interval changes in ADCmean did not differ between metastatic lymph nodes that did and did not respond to treatment. This may be due to our inclusion of necrotic areas in the ROIs, which could have led to an overestimation of baseline ADCmean.

King et al. found that ADC histogram analysis could predict intra-treatment response at the primary site as early as 2 weeks after the start of treatment, with primary tumours showing higher histogram skewness (20). de Parrot et al. demonstrated that the ADC histogram kurtosis ratio (22) was useful in differentiating histologic grades of HNSCC at b = 2,000 s/mm2. A general shift in skewness to the left means that the tail of the ADC curves got heavier toward the negative and signifies a greater asymmetry of the data posttreatment. Similar findings were reported by Forouton et al. (23), who demonstrated negative skewness and greater asymmetry following therapy in the treatment group when compared to the control group. Whilst our data showed a non-significant increase in skewness post CRT, the histogram parameters of kurtosis, coefficient of variation and energy were found to demonstrate a significant increase. This was demonstrated by 6-week post CRT for primary tumour and 12 weeks post CRT for lymph nodes. The underlying histopathological basis for this is yet to be studied; however, it appears there is a trend toward greater homogeneity of the tumour with treatment.

Post-Gadolinium T1w MRI

There are limited data on the interval changes in macroscopic T1w gadolinium enhancement following CRT for HNSCC. A qualitative grading (24) showed that changes in the pattern of necrosis on T1w gadolinium enhancement imaging could predict the outcome. The GLCM is a tabulation of how often different combinations of grey levels occur between neighbouring pixels in an image (2123). Reinert et al. (25) found that GLCM entropy significantly decreased, whereas GLCM uniformity significantly increased (p   < .001) after therapy in a study assessing response of Hodgkin lymphoma to chemotherapy. Whilst none of the histogram parameters changed significantly post CRT in our cohort, GLCM contrast increased whilst GLCM entropy decreased with treatment in both the primary tumour and lymph node by 6 weeks post CRT. It is of interest that changes were observed in similar parameters in both metastatic lymph nodes (largely non-enhancing necrosis in 16/18 pretreatment) and primary tumours (homogenously enhancing in 22/23 pretreatment) considering the marked difference in pretreatment enhancement patterns.

Differences Between the Primary Tumour and Nodes

The interval changes in ADC-based histogram parameters differed between the primary tumours and metastatic lymph nodes, with earlier changes being observed with primary tumours. Although there are no previous comparable data available for ADC-based histogram parameters, a potential disparity between the timing and nature of posttreatment ADC-related biomarkers between the primary tumour and nodes has been established in previous studies. In particular, these have demonstrated that the ADCmean acquired from nodal sites is able to predict treatment outcomes, whereas ADCs acquired from primary sites do not (26), whilst predictive threshold ADC measurements differ between primary tumour and nodal sites (27).

Comparison With Primary Tumour and Nodal Size and Volume Change

The assessment of morphological features such as changes in tumour dimensions has provided prognostic stratification (20, 24, 28, 29) in the posttreatment setting, but they have variable diagnostic accuracy in determining treatment response (30, 31). A significant reduction in lymph node dimensions at 6 weeks post CRT was only shown in those with successful treatment. However, it is of interest that the difference in magnitude between responding and non-responding lymph node dimensions was less than that demonstrated with some heterogeneity parameters. The decrease in diameter and volume of the responding lymph nodes at 6 weeks post CRT was 69% and 95% as compared to 50% and 90% with the non-responding lymph nodes. Thus, the addition of histogram parameters could perhaps aid in increasing the discriminatory power of size criteria in assessing posttreatment changes, even as early as 6 weeks posttreatment.

Limitations

It is acknowledged that there are some limitations to this study.

Firstly, it should be noted that the cohort was largely represented by HPV-positive oropharyngeal cancer. HPV-OPC is a clinically and histologically distinct form of HNSCC which has a better outcome irrespective of treatment choice. Since it exhibits particular histopathological features such as indistinct cell borders and comedo-necrosis, there should be caution in extrapolation of the results to a wider HNSCC population.

Secondly, there was only a small sample of non-responders to CRT with only four metastatic lymph nodes and no primary tumours which did not respond to treatment. The small sample of non-responding lymph nodes precluded a direct statistical comparison of parameters between responding and non-responding lymph nodes. It also raises the possibility that interval changes in histogram and GLCM parameters in the non-responding lymph nodes were not detected due to a type 2 error. Future studies should be conducted with a larger HPV negative cohort and a greater proportion of non-responders such that these texture parameters can be evaluated for their ability to predict treatment response.

Finally, it could be argued that the outcomes were influenced by our method of ROI analysis which was focused on whole volumes of the primary tumour or lymph node. Both volume and area-based ROI data are variably presented in the literature. However, there are limited data on how this impacts on histogram or GLCM analysis. Generous sampling of the whole volume may better reflect tumour heterogeneity. However, smaller ROIs benefit from the ability to exclude macroscopic necrosis. It is possible that sampling certain parameters such as ADCmean may have revealed significant changes if confined to enhancing “viable tumour”, particularly with respect to the largely necrotic lymph nodes.

Conclusions

This exploratory study indicates that significant post-CRT interval changes in histogram and GLCM parameters are demonstrated on ADC maps and T1w post gad images as early as 6 weeks following successful treatment for stage III and IV HNSCC. The selection of appropriate ADC map histogram parameters may differ, depending on whether biomarkers for primary tumour or lymph node response are being evaluated. A small sample of non-responding metastatic lymph nodes did not demonstrate such interval evolution of these heterogeneity parameters; however, standard size criteria were also able to predict treatment failure in these patients.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Ethics Statement

Institutional approval from the Research Ethics Committee (REC reference 13/LO/1876) and informed consent was obtained from all participants. The study conforms to recognized standards of the Declaration of Helsinki.

Author Contributions

Conceptualization, SC, VG; Methodology, SC, VG, GC; Investigation, SC, AG, MA; Data Curation, AG, SC; Writing - Original Draft, AG, SC; Writing - Review & Editing, VG, GC; Supervision, SC, VG; Project Administration, SC: Funding Acquisition, SC.

Funding

This research was funded in whole, or in part, by the Wellcome Trust [203148/Z/16/Z]. For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission. Authors acknowledge funding support from Wellcome/Engineering and Physical Sciences Research Council Centre for Medical Engineering at King’s College London (WT 203148/Z/16/Z); National Institute for Health Research Biomedical Research Centre at Guy’s & St Thomas’ Hospitals and King’s College London; Cancer Research UK National Cancer Imaging Translational Accelerator (A27066); the UK Research & Innovation London Medical Imaging and Artificial Intelligence Centre. Authors also acknowledge funding suport from Guy’s and St Thomas’ Hospital Charity (ref EFT130501) and the Royal College of Radiologists: Kodak Radiology Fund Research Bursary.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.784693/full#supplementary-material

References

1. Ferlay J, Soerjomataram I, Dikshit R, Esers S, Mathers C, Rebelo M, et al. Cancer Incidence and Mortality Worldwide: Sources, Methods and Major Patterns in GLOBOCAN 2012. Int J Cancer (2015) 136(5):E359–86. doi: 10.1002/ijc.29210

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Bernier J, Pfister DG, Cooper JS. Adjuvant Chemo- and Radiotherapy for Poor Prognosis Head and Neck Squamous Cell Carcinomas. Crit Rev Oncol Hematol (2005) 56(3):353–64. doi: 10.1016/j.critrevonc.2005.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hauswald H, Simon C, Hecht S, Debus J, Lindel K. Long-Term Outcome and Patterns of Failure in Patients With Advanced Head and Neck Cancer. Radiat Oncol (2011) 6:70. doi: 10.1186/1748-717X-6-70

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mabanta SR, Mendenhall WM, Stringer SP, Cassisi NJ. Salvage Treatment for Neck Recurrence After Irradiation Alone for Head and Neck Squamous Cell Carcinoma With Clinically Positive Neck Nodes. Head Neck (1999) 21(7):591–4. doi: 10.1002/(SICI)1097-0347(199910)21:7<591::AID-HED1>3.0.CO;2-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Varoquaux A, Rager O, Dulguerov P, Burkhardt K, Ailianou A, Becker M. Diffusion-Weighted and PET/MR Imaging After Radiation Therapy for Malignant Head and Neck Tumors. RadioGraphics (2015) 35(5):1502–27. doi: 10.1148/rg.2015140029

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Subramaniam RM, Truong M, Peller P, Sakai O, Mercier G. Fluorodeoxyglucose–Positron-Emission Tomography Imaging of Head and Neck Squamous Cell Cancer. Am J Neuroradiol (2010) 31(4):598–604. doi: 10.3174/ajnr.A1760

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gupta T, Master Z, Kannan S, Agarwal JP, Ghsoh-Laskar S, Rangarajan V, et al. Diagnostic Performance of Post-Treatment FDG PET or FDG PET/CT Imaging in Head and Neck Cancer: A Systematic Review and Meta-Analysis. Eur J Nucl Med Mol Imaging (2011) 38(11):2083. doi: 10.1007/s00259-011-1893-y

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Just N. Improving Tumour Heterogeneity MRI Assessment With Histograms. Br J Cancer (2014) 111(12):2205–13. doi: 10.1038/bjc.2014.512

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kim YJ, Kim SH, Lee AW, Jin M-S, Kang BJ, Song BJ. Histogram Analysis of Apparent Diffusion Coefficients After Neoadjuvant Chemotherapy in Breast Cancer. Jpn J Radiol (2016) 34(10):657–66. doi: 10.1007/s11604-016-0570-2

PubMed Abstract | CrossRef Full Text | Google Scholar

10. King AD, Chow K-K, Yu K-H, Mo FKF, Yeung DKW, Yuan J, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ren J-L, Yuan Y, Li X-X, Shi Y-Q, Tao X-F. Histogram Analysis of Apparent Diffusion Coefficient Maps in the Prognosis of Patients With Locally Advanced Head and Neck Squamous Cell Carcinoma: Comparison of Different Region of Interest Selection Methods. Eur J Radiol (2018) 106:7–13. doi: 10.1016/j.ejrad.2018.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Jacob R, Welkoborsky HJ, Mann WJ, Höfken F, Dienes HP, Freije JE. Heterogeneity of Squamous Cell Carcinomas of the Head and Neck-Analysis of Tumor Biologic Factors and Proliferation Rates. Laryngoscope (1996) 106(9):1170–5. doi: 10.1097/00005537-199609000-00023

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Veneroni S, Silvestrini R, Costa A, Salvatori P, Faranda A, Molinari R. Biological Indicators of Survival in Patients Treated by Surgery for Squamous Cell Carcinoma of the Oral Cavity and Oropharynx. Oral Oncol (1997) 33(6):408–13. doi: 10.1016/S0964-1955(97)00036-5

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zwanenburg A, Vallières M, Abdalah MA, Aerts HJWL, Andrearczyk V, Apte A, et al. The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Image-Based Phenotyping. Radiology (2020) 295(2):328–38. doi: 10.1148/radiol.2020191145

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kim S, Loevner L, Quon H, Sherman E, Weinstein G, Kilger A, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Marzi S, Piludu F, Sanguineti G, Marucci L, Farneti A, Terrenato I, et al. The Prediction of the Treatment Response of Cervical Nodes Using Intravoxel Incoherent Motion Diffusion-Weighted Imaging. Eur J Radiol (2017) 92:93–102. doi: 10.1016/j.ejrad.2017.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Berrak S, Chawla S, Kim S, Quon H, Sherman E, Loevner LA, et al. Diffusion Weighted Imaging in Predicting Progression Free Survival in Patients With Squamous Cell Carcinomas of the Head and Neck Treated With Induction Chemotherapy. Acad Radiol (2011) 18(10):1225–32. doi: 10.1016/j.acra.2011.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Matoba M, Tuji H, Shimode Y, Toyoda I, Kuginuki Y, Miwa K, et al. Fractional Change in Apparent Diffusion Coefficient as an Imaging Biomarker for Predicting Treatment Response in Head and Neck Cancer Treated With Chemoradiotherapy. Am J Neuroradiol (2014) 35(2):379–85. doi: 10.3174/ajnr.A3706

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Galbán CJ, Mukherji SK, Chenevert TL, Meyer CR, Hamstra DA, Bland PH, et al. A Feasibility Study of Parametric Response Map Analysis of Diffusion-Weighted Magnetic Resonance Imaging Scans of Head and Neck Cancer Patients for Providing Early Detection of Therapeutic Efficacy. Transl Oncol (2009) 2(3):184–90. doi: 10.1593/tlo.09175

PubMed Abstract | CrossRef Full Text | Google Scholar

20. King AD, Mo FKF, Yu K-H, Yeung DKW, Zhou H, Bhatia KS, et al. Squamous Cell Carcinoma of the Head and Neck: Diffusion-Weighted MR Imaging for Prediction and Monitoring of Treatment Response. Eur Radiol (2010) 20(9):2213–20. doi: 10.1007/s00330-010-1769-8

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Vandecaveye V, Dirix P, De Keyzer F, Op de Beeck K, Vander Poorten V, Hauben E, et al. Diffusion-Weighted Magnetic Resonance Imaging Early After Chemoradiotherapy to Monitor Treatment Response in Head-and-Neck Squamous Cell Carcinoma. Int J Radiat Oncol (2012) 82(3):1098–107. doi: 10.1016/j.ijrobp.2011.02.044

CrossRef Full Text | Google Scholar

22. de Perrot T, Lenoir V, Domingo Ayllón M, Dulguerov N, Pusztaszeri M, Becker M. 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. Am J Neuroradiol (2017) 38(11):2153–60. doi: 10.3174/ajnr.A5370

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Foroutan P, Kreahling JM, Morse DL, Grove O, Lloyd MC, Reed D, et al. Diffusion MRI and Novel Texture Analysis in Osteosarcoma Xenotransplants Predicts Response to Anti-Checkpoint Therapy. PloS One (2013) 8(12):e82875. doi: 10.1371/journal.pone.0082875

PubMed Abstract | CrossRef Full Text | Google Scholar

24. King AD, Yu K-H, Mo FKF, Law BKH, Yeun TWC, Bhatia KS, et al. Cervical Nodal Metastases From Head and Neck Squamous Cell Carcinoma: MRI Criteria for Treatment Assessment. Head Neck (2016) 38(S1):E1598–604. doi: 10.1002/hed.24285

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Reinert CP, Wanek L, Bösmüller H, Federmann B, Fritz J, Sokler M, et al. Computed Tomography Texture Analysis for Assessment of Chemotherapy Response of Hodgkin Lymphoma. Med (Baltimore) (2020) 99(7):e19146–6. doi: 10.1097/MD.0000000000019146

CrossRef Full Text | Google Scholar

26. Noij DP, Pouwels PJW, Ljumanovic R, Knol DL, Doornaert P, de Bree R, et al. Predictive Value of Diffusion-Weighted Imaging Without and With Including Contrast-Enhanced Magnetic Resonance Imaging in Image Analysis of Head and Neck Squamous Cell Carcinoma. Eur J Radiol (2015) 84(1):108–16. doi: 10.1016/j.ejrad.2014.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wong KH, Panek R, Bhide SA, Nutting CM, Harrington KJ, Newbold KL. The Emerging Potential of Magnetic Resonance Imaging in Personalizing Radiotherapy for Head and Neck Cancer: An Oncologist’s Perspective. Br J Radiol (2017) 90(1071):20160768. doi: 10.1259/bjr.20160768

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Bhatia KSS, King AD, Yu K-H, Vlantis AC, Tse GM, Mo FKF, et al. Does Primary Tumour Volumetry Performed Early in the Course of Definitive Concomitant Chemoradiotherapy for Head and Neck Squamous Cell Carcinoma Improve Prediction of Primary Site Outcome? Br J Radiol (2010) 83(995):964–70. doi: 10.1259/bjr/27631720

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Labadie RF, Yarbrough WG, Weissler MC, Pillsbury HC, Mukherji SK. Nodal Volume Reduction After Concurrent Chemo- and Radiotherapy: Correlation Between Initial CT and Histopathologic Findings. Am J Neuroradiol (2000) 21(2):310–4.

PubMed Abstract | Google Scholar

30. Clavel S, Charron M-P, Bélair M, Delouya G. The Role of Computed Tomography in the Management of the Neck After Chemoradiotherapy in Patients With Head-and-Neck Cancer. Int J Radiat Oncol Biol Phys (2012) 82(2):567–73. doi: 10.1016/j.ijrobp.2010.11.066

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Nishimura G, Matsuda H, Taguchi T, Takahashi M, Komatsu M, Sano D, et al. Treatment Evaluation of Metastatic Lymph Nodes After Concurrent Chemoradiotherapy in Patients With Head and Neck Squamous Cell Carcinoma. Anticancer Res (2012) 32(2):595–600.

PubMed Abstract | Google Scholar

Keywords: radiomics, heterogeneity, head and neck cancer, chemoradiotherapy, texture analysis, magnetic resonance imaging

Citation: Guha A, Anjari M, Cook G, Goh V and Connor S (2022) Radiomic Analysis of Tumour Heterogeneity Using MRI in Head and Neck Cancer Following Chemoradiotherapy: A Feasibility Study. Front. Oncol. 12:784693. doi: 10.3389/fonc.2022.784693

Received: 28 September 2021; Accepted: 13 January 2022;
Published: 15 February 2022.

Edited by:

Shreerang Bhide, Institute of Cancer Research (ICR), United Kingdom

Reviewed by:

Derfel Ap Dafydd, Royal Marsden Hospital, United Kingdom
Avraham Eisbruch, University of Michigan, United States

Copyright © 2022 Guha, Anjari, Cook, Goh and Connor. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Steve Connor, steve.connor@kcl.ac.uk

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.