Abstract

Background. Our purpose was to elucidate possible correlations between histogram parameters derived from dynamic contrast-enhanced MRI (DCE-MRI) with several histopathological features in head and neck squamous cell carcinomas (HNSCC). Methods. Thirty patients with primary HNSCC were prospectively acquired. Histogram analysis was derived from the DCE-MRI parameters: Ktrans, Kep, and Ve. Additionally, in all cases, expression of human papilloma virus (p16) hypoxia-inducible factor-1-alpha (Hif1-alpha), vascular endothelial growth factor (VEGF), epidermal growth factor receptor (EGFR), and tumor suppressor protein p53 were estimated. Results. Kep kurtosis was significantly higher in p16 tumors, and Ve min was significantly lower in p16 tumors compared to the p16 negative tumors. In the overall sample, Kep entropy correlated well with EGFR expression (, ). In p16 positive carcinomas, Ktrans max correlated with VEGF expression (, ), Ktrans kurtosis correlated with Hif1-alpha expression (, ), and Ktrans entropy correlated with EGFR expression (, ). Regarding Kep parameters, mode correlated with VEGF expression (, ), and entropy correlated with Hif1-alpha expression (, ). In p16 negative carcinomas, Kep mode correlated with Her2 expression (, ), Ve max correlated with p53 expression (, ), and Ve p10 correlated with EGFR expression (, ). Conclusion. DCE-MRI can reflect several histopathological features in HNSCC. Associations between DCE-MRI and histopathology in HNSCC depend on p16 status. Kep kurtosis and Ve min can differentiate p16 positive and p16 negative carcinomas.

1. Introduction

Head and neck squamous cell carcinoma (HNSCC) is a frequently occurring malignancy [1]. Previously, the role of imaging modalities was to locate the primary tumor and detect infiltration of bordering body structures and distant metastasis [2]. However, modern imaging modalities can also provide valuable information regarding tumor microstructure and might be able to predict several histopathological features in tumors [3, 4].

Dynamic contrast-enhanced MRI (DCE-MRI) is a functional imaging technique, which is able to assess tumor vascularization by measurement of sequential changes of signal intensity over time after contrast media application [5, 6]. In DCE-MRI, quantitative parameters like Ktrans (volume transfer constant in min−1), Ve (volume fraction of the extravascular extracellular space which is dimensionless), and Kep (rate constant in min−1) can be obtained [6].

Previous reports suggested that DCE-MRI can reflect tumor vessel density [6]. However, besides perfusion, DCE-MRI is also linked to cellularity, as well as to proliferation index [7, 8]. Furthermore, it has been shown that DCE-MRI can predict survival and treatment response to radiochemotherapy in HNSCC [5, 911]. Additionally, it can predict tumor recurrence [12] and metastatic spread [13]. Besides the prognostic information, DCE-MRI can also aid in discrimination between benign and malignant head and neck tumors [14].

Histogram analysis is used to analyze radiological images. By using this technique, every voxel of a region of interest (ROI) is issued into a histogram. Thereby, a broad spectrum of new parameters can be estimated: minimum, mean, maximum, median, mode, percentiles, kurtosis, skewness, and entropy. According to the literature, heterogeneity of the histogram might also display heterogeneity of the tumor [15].

Several histopathological parameters play an important role in HNSCC. For example, p16 expression, associated with human papilloma virus, is one of the most important prognostic factors in HNSCC [16]. Other parameters, such as vascular endothelial growth factor (VEGF), hypoxia-inducible factor-1-alpha (Hif1-alpha), epidermal growth factor receptor (EGFR), and tumor suppressor protein p53 expression, are also of prognostic relevance and might aid in treatment response prediction in HNSCC [17, 18]. Presumably, imaging might also be able to reflect these expression profiles, especially by using the more advanced histogram-based analysis. Recently, a first promising study identified statistical differences between p16 positive and p16 negative carcinomas using histogram-based parameters derived from diffusion-weighted imaging [19]. Previously, only two studies analyzed relationships between DCE-MRI and histopathological parameters like the proliferation index Ki 67 and/or tumor cellularity in HNSCC using conventional ROI-based analysis [7, 20]. Presumably, histogram-based DCE parameters may show more associations with histopathology.

Therefore, the aim of this study was to estimate whole lesion histogram parameters derived from DCE-MRI and to elucidate possible correlations with several clinically relevant histopathological features in HNSCC.

2. Materials and Methods

This prospective study was approved by the institutional review board (Ethics committee of the University of Leipzig, study codes 180-2007, 201-10-12072010, and 341-15-05102015). All methods were performed in accordance with the relevant guidelines and regulations. All patients gave their written informed consent.

2.1. Patients

For this study, 30 patients (22 men and 8 women; mean age 57.0 ± 10.6 years; range 33–77 years) with histopathological proven primary HNSCC were included into the present study. Different tumor localizations were identified: the oropharynx in 46.7% of cases, tongue in 23.3%, hypopharynx in 10%, larynx in 16.7%, and nasopharynx in 3.3% of cases. There were T3 staged cancers in 33.3% and T4 in 40% cases and only 26.7% with T1 and T2 cancers. 90% of cases were nodal positive and 10% of patients without any nodal metastases. Well and moderately differentiated tumors were identified in 36.7% of patients and poorly differentiated in 63.3%. All patients did not receive any form of cancer treatment before the investigation.

2.2. DCE-MRI

In all patients, dynamic contrast-enhanced (DCE) imaging was performed using T1w DCE sequences according to a imaging protocol, as reported previously (TR/TE 2.47/0.97 ms, flip angle 8°, voxel size 1.2 × 1.0 × 5.0 mm, and slice thickness 5 mm) [7, 21]. The sequence included forty scans at 6 seconds. The contrast application of 0.1 mmol gadobutrol per kg of bodyweight (Gadovist®, Bayer Healthcare, Leverkusen, Germany) started after the fifth scan with a rate of 3 ml per second (Spectris Solaris, Medrad, Bayer Healthcare, Leverkusen, Germany). The acquired images were further analyzed with Tissue 4D (Siemens Medical Systems, Erlangen, Germany), which uses a population-based technique for the arterial input function (AIF). The AIF was modelled to the gadolinium dose and according to the biexponential model of Tofts and Kermode. Finally, Ktrans, Ve, and Kep were calculated (for exemplary parameter images, see Figures 1 and 2).

2.3. Histogram Analysis

The acquired DCE-MRI data were processed with a Matlab-based application (Mathworks, Natick, MA, USA). On the Ktrans, Kep, and Ve maps, a volume of interest was drawn inside the tumor boundary using all slices with visible tumor areas and thus providing a whole lesion measurement. All measures were performed by one experienced author (AS, 15 years of general radiological experience). The following parameters were estimated for Ktrans, Kep, and Ve: mean, maximum, minimum, median, mode, 10th, 25th, 75th, and 90th percentiles, as well as kurtosis, skewness, and entropy.

2.4. Histopathological Findings

In every patient, the diagnosis was confirmed by tumor biopsy. The histological specimens were deparaffinized, rehydrated, and cut into 5 μm slices. Moreover, the histological slices were stained by the epidermal growth factor receptor (EGFR, EMERGO Europe, clone 111.6, dilution 1 : 30), vascular endothelial growth factor (VEGF, EMERGO Europe, clone VG1, dilution 1 : 20), tumor suppressor protein p53 (DakoCytomation, Glostrup, Denmark; clone DO-7, dilution 1 : 100), hypoxia-inducible factor-1 (Hif1-alpha) (Biocare Medical, 60 Berry Dr Pacheco, CA 94553; clone EP1215Y, dilution 1 : 100), and p16 (p16 expression, CINtec Histology, Roche, Germany), as performed in our previous study [22].

Pannoramic microscope scanner (Pannoramic SCAN, 3DHISTECH Ltd., Budapest, Hungary) with Carl Zeiss objectives up to 41x bright field magnification by default was used to digitalize all specimens. In the used bottom-up technique, the whole sample was acquired at a high resolution. All slides were analyzed with Pannoramic Viewer 1.15.4 (open source software, 3D HISTECH Ltd., Budapest, Hungary), and three representative images with a magnification of ×200 were extracted from each patient.

The histopathological images were further investigated by using the ImageJ software 1.48v (National Institutes of Health Image program). The tumors were divided according to the p16 status.

Finally, expression of EGFR, VEGF, HIF1-alpha, and p53 (Figures 1 and 2) was semiautomatically estimated as a sum of stained areas (in µm2) by using a brightness threshold. Figure 1 displays a p16 negative, and Figure 2 shows a p16 positive carcinoma.

2.5. Statistical Analysis

Statistical analysis was performed using GraphPad Prism (GraphPad Software, La Jolla, CA, USA). Collected data were evaluated by means of descriptive statistics.

Spearman's correlation coefficient (ρ) was used to analyze associations between investigated imaging and histopathology parameters. Mann–Whitney U test was used for discrimination between p16 groups. values below 0.05 were considered statistically significant.

3. Results

There were 10 (33.3%) p16 negative and 20 (66.7%) p16 positive tumors. Kep kurtosis was significantly higher in p16 tumors, and Ve min was significantly lower in p16 positive tumors compared to the p16 negative tumors, and , respectively (Figure 3).

In the overall sample, the correlation analysis revealed only one statistically significant correlation between Kep entropy and EGFR expression (ρ = 0.38, ) (Figure 4).

In the p16 positive carcinomas, Ktrans max correlated with VEGF expression (ρ = 0.46, ), Ktrans kurtosis correlated with Hif1-alpha expression (ρ = 0.46, ) and Ktrans entropy correlated with EGFR expression (ρ = 0.50, ). Regarding Kep parameters, mode correlated with VEGF expression (ρ = 0.51, ), and entropy correlated with Hif1-alpha expression (ρ = 0.47, ). None of the Ve values were associated with the analyzed histochemical parameters.

In the p16 negative group, the following associations could be identified: Kep mode correlated with Her2 expression (ρ = −0.72, ), Ve max correlated with p53 expression (ρ = −0.80, ), and Ve p10 correlated with EGFR expression (ρ = 0.68, ).

4. Discussion

This present study identified statistically significant associations between histogram parameters derived from DCE-MRI and different histopathological features in HNSCC. Furthermore, it showed that these relationships depended on the p16 status.

There is increasing evidence that MRI, especially using functional imaging modalities, is able to reflect tumor microstructure and to predict tumor behavior [3, 7, 8, 20]. It is widely acknowledged that DCE-MRI is associated with vascularity in tissues, especially with microvessel density as the most investigated parameter. For example, significant associations between DCE-MRI and microvessel density have been reported in experimental [23] as well as in clinical investigations [7, 24, 25].

Notably, it has been shown that different DCE parameters might also reflect different aspects of tumor microstructure [7]. So, Ve might also be strongly associated with cellularity because it reflects the amount of extracellular space, as it was exemplarily shown in a glioma model [8]. This might be one reason for the different correlations identified in the present study.

Several studies elucidated possible correlations between imaging and histopathology in HNSCC. For example, it has been shown that diffusion-weighted imaging (DWI) correlated with Ki 67 expression as well with nucleic areas [3, 26]. In another study, Ktrans correlated inversely with Ki 67 expression (r = −0.62), whereas Ve tended to correlate with the cell count [7]. Furthermore, Jansen et al. showed that Kep correlated statistically significant with VEGF expression (r = 0.808) [20].

In the present study, Kep mode correlated with VEGF expression in p16 positive patients. Interestingly, also Ktrans max correlated in a similar fashion with VEGF expression. Furthermore, Ktrans max also showed a significant association with Hif1-alpha. Presumably, the maximum value of Ktrans may reflect tumor areas with the highest vessel density. Therefore, the observed correlation between Ktrans max and expression of VEGF is logical. Our results are in agreement with some previous reports. For example, in gliomas, also a positive correlation between VEGF and Ktrans was observed [2729].

However, some studies did not find significant associations between DCE-MRI and histopathology. For example, in breast cancer, no correlations between histogram parameters derived from DCE-MRI and VEGF expression could be identified [24].

Rasmussen et al. found associations between standardized uptake values (SUV) derived from positron emission tomography (PET) with 2-deoxy-2-[18F]fluoro-D-glucose (FDG) and histopathology in HNSCC [30]. There were negative correlations for Bcl-2 and p16 and positive with β-tubulin-1 index. Moreover, in another study, SUV was only associated with VEGF expression, whereas no association was found for GLUT-1, Ki 67, P53, CD68, Hif1-alpha, and CD31 [31]. Our results indicate that DCE-MRI might be more sensitive than FDG PET for prediction of histopathological features.

It is believed that the histogram-based analysis of radiological images can better reflect tumor than conventional ROI-based analysis [15]. For example, it was shown that histogram analysis of DCE and DWI can identify more correlations between parameters of these imaging modalities [32].

The present study showed that kurtosis values derived from Kep and Ve min were significantly different in p16 positive compared to p16 negative tumors. This novel finding might be caused by several underlying tissue characteristics. In a recent study by de Perrot et al., histogram analysis derived from the ADC map was used to differentiate between p16 positive and p16 negative HNSCSS [19]. Ve is a parameter, which might be related to ADC values and cellularity [8, 30]. Interestingly, Ve min that represents voxels with the lowest extracellular space, and, presumably, areas with the highest cell density, was lower in p16 positive lesions. This finding may suggest that p16 positive tumors may show a higher cell density than p16 negative tumors. In the study by de Perrot et al., also kurtosis derived from ADC maps could distinguish p16 positive and p16 negative carcinomas [19].

These findings might be related to several causes. As reported previously, p16 positive cancers were more often nonkeratinizing and had a high Ki 67 expression [19]. Moreover, expression profiles of p16 positive and p16 negative cancers might differ significantly emphasizing their different tumor behavior. So, it was shown that expression of Eps8 is different in these subtypes of HNSCC [33]. This EGFR substrate contributes to the carcinogenesis and might be involved in invasiveness in HNSCC [31]. Interestingly, the expression of Eps8 correlated with the tumor stage and p16 status but not with anatomical localization of tumors [33]. Moreover, the expression of other histopathological parameters such as EGFR, VEGF, and NOTCH1 differ between p16 positive and negative tumors, which suggest differences in tumor angiogenesis in these entities [34]. This might be also a reason for the identified influence of p16 expression on association between imaging and histopathology.

Furthermore, it is known that p16 expression is one of the most important prognostic factors in HNSCC with a more favorable outcome for p16 positive cancers [16]. The other investigated histopathological features are also of clinical importance. So, EGFR is involved in the regulation of many cellular pathways, including cell proliferation, apoptosis, and cellular differentiation [35]. It was identified that EGFR expression is a good prognostic parameter in HNSCC [35, 36]. Furthermore, p53 regulates the activity of pathways, which lead to cell cycle arrest, senescence, or apoptosis [37]. Another parameter, namely, VEGF predicts outcome in HNSCC. VEGF overexpression has been reported as a poor indicator for patients with head and neck cancer [38]. Finally, Hif1-alpha characterizes cellular responses to hypoxic stress and is related to the neoangiogenesis [39]. Overexpression of Hif1-alpha was also significantly associated with poor survival in HNSCC [39]. Therefore, the possibility to characterize HNSCC based on imaging is very important. The identified associations between DCE-MRI parameters and several histopathological markers can be used in clinical practice.

There are several limitations of this study to address. Firstly, our patient sample size is small yet good comparable to similar studies. Secondly, we performed a whole tumor measurement for the DCE-MRI images, whereas the histopathology was investigated only on a small part of the tumor, which might limit our correlation results. Further prospective studies are needed to confirm our preliminary results.

In conclusion, the present study identified statistically significant correlations between histogram parameters derived from DCE-MRI and expression of VEGF, EGFR, p53, and Hif1-alpha in HNSCC. Associations between DCE-MRI and histopathology in HNSCC depend on the p16 status. Furthermore, Kep kurtosis and Ve minimum can differentiate p16 positive and p16 negative carcinomas.

Data Availability

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

Ethical Approval

The study was approved by the institutional review board of the University of Leipzig. All procedures performed in the study were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Conceptualization was done by Alexey Surov. Data curation was performed by Leonard Leifels, Alexey Surov, and Hans Jonas Meyer. Formal analysis was conducted by Leonard Leifels, Hans Jonas Meyer, and Anne Kathrin Höhn. Investigation was performed by Leonard Leifels, Alexey Surov, and Anne Kathrin Höhn. Methodology was contributed by Hans Jonas Meyer and Alexey Surov. Project administration was done by Alexey Surov. Resources were contributed by Alexey Surov and Hans Jonas Meyer. Software was contributed by Hans Jonas Meyer and Alexey Surov. Supervision was done by Alexey Surov. Validation was performed by Alexey Surov and Hans Jonas Meyer. Visualization was performed by Alexey Surov. Writing of the original draft was performed by Hans Jonas Meyer. Writing in terms of review and editing was performed by Leonard Leifels, Alexey Surov, Hans Jonas Meyer, and Anne Kathrin Höhn.