Next Article in Journal
Radiogenomic Analysis of Oncological Data: A Technical Survey
Next Article in Special Issue
Diagnostic Limitation of Fine-Needle Aspiration (FNA) on Indeterminate Thyroid Nodules Can Be Partially Overcome by Preoperative Molecular Analysis: Assessment of RET/PTC1 Rearrangement in BRAF and RAS Wild-Type Routine Air-Dried FNA Specimens
Previous Article in Journal
Potential Coagulation Factor-Driven Pro-Inflammatory Responses in Ovarian Cancer Tissues Associated with Insufficient O2 and Plasma Supply
Previous Article in Special Issue
Molecular Signature of Indeterminate Thyroid Lesions: Current Methods to Improve Fine Needle Aspiration Cytology (FNAC) Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Histogram Analysis of Diffusion Weighted Imaging at 3T is Useful for Prediction of Lymphatic Metastatic Spread, Proliferative Activity, and Cellularity in Thyroid Cancer

1
Department for Neuroradiology, University Hospital Leipzig, Leipzig 04103, Germany
2
Department for Diagnostic and Interventional Radiology, University Hospital Leipzig, Leipzig 04103, Germany
3
Department for Ophthalmology, University Hospital Leipzig, Leipzig 04103, Germany
4
Department for Diagnostic and Interventional Neuroradiology, Katharinenhospital Stuttgart, Stuttgart 70174, Germany
5
Institute for Pathology, University Hospital Halle-Wittenberg, Martin-Luther-University Halle-Wittenberg, Halle 06112, Germany
6
Institute for Pathology, University Hospital Leipzig, Leipzig 04103, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(4), 821; https://doi.org/10.3390/ijms18040821
Submission received: 9 March 2017 / Revised: 9 April 2017 / Accepted: 10 April 2017 / Published: 12 April 2017
(This article belongs to the Special Issue Current Knowledge in Thyroid Cancer—From Bench to Bedside)

Abstract

:
Pre-surgical diffusion weighted imaging (DWI) is increasingly important in the context of thyroid cancer for identification of the optimal treatment strategy. It has exemplarily been shown that DWI at 3T can distinguish undifferentiated from well-differentiated thyroid carcinoma, which has decisive implications for the magnitude of surgery. This study used DWI histogram analysis of whole tumor apparent diffusion coefficient (ADC) maps. The primary aim was to discriminate thyroid carcinomas which had already gained the capacity to metastasize lymphatically from those not yet being able to spread via the lymphatic system. The secondary aim was to reflect prognostically important tumor-biological features like cellularity and proliferative activity with ADC histogram analysis. Fifteen patients with follicular-cell derived thyroid cancer were enrolled. Lymph node status, extent of infiltration of surrounding tissue, and Ki-67 and p53 expression were assessed in these patients. DWI was obtained in a 3T system using b values of 0, 400, and 800 s/mm2. Whole tumor ADC volumes were analyzed using a histogram-based approach. Several ADC parameters showed significant correlations with immunohistopathological parameters. Most importantly, ADC histogram skewness and ADC histogram kurtosis were able to differentiate between nodal negative and nodal positive thyroid carcinoma. Conclusions: histogram analysis of whole ADC tumor volumes has the potential to provide valuable information on tumor biology in thyroid carcinoma. However, further studies are warranted.

1. Introduction

The incidence of thyroid cancer, being the most abundant endocrine malignancy, is rapidly increasing [1]. The vast majority of thyroid neoplasms is follicular cell-derived and subsumed under the umbrella categories of papillary thyroid cancer, follicular thyroid cancer, poorly differentiated thyroid cancer, and anaplastic thyroid cancer [2]. Although the overall five-year survival rates of thyroid cancer are 94% in women and 85% in men [3], certain entities of the disease are perpetually associated with poor outcomes (for example the tall cell variant of papillary thyroid cancer and undifferentiated thyroid cancer [1]). Some of the differentiated entities—most of all papillary thyroid cancer variants—frequently metastasize locally via the lymphatic system [4], and resultant local recurrence is not an uncommon scenario [5], leading to significant morbidity.
A variety of therapeutic options is available for thyroid cancer [6], but surgery still remains the predominant treatment [7]. Radical surgery is the most important form of therapy for undifferentiated thyroid cancer [8], and surgical treatment of significant nodal disease in well differentiated thyroid cancer is widely accepted to be associated with improved outcomes in terms of survival and recurrence rates [9]. Nonetheless, extensive surgery in this specific context carries a high risk of therapy-related morbidity like phrenic nerve palsy, brachial plexus palsy, cranial nerve injury, chyle leak, and pneumothorax [10].
Considering the broad spectrum of aggressiveness in thyroid cancer and the resulting necessity for customized treatment, employing presurgical imaging is of great importance, as it allows the thyroid surgeon to identify disease subtypes being associated with increased mortality and morbidity such as metastasizing and undifferentiated thyroid cancer.
Diffusion-weighted magnetic resonance imaging (DWI) has the potential to reveal tumor architectural details like cellular density and proliferative activity in different malignant entities [11,12]. Using a standard echo-planar imaging (EPI) technique, DWI has the capability to differentiate between malignant and benign thyroid nodules [13]. Furthermore, DWI can distinguish manifestations of papillary thyroid cancer with extra-glandular growth from those confined to the thyroid [14]. Using a RESOLVE sequence (which is less prone to susceptibility and motion-induced phase artifacts, has less T2* blurring and provides higher resolution than standard EPI DWI, [15]) in a 3T scanner, DWI even has the capability to distinguish between differentiated and undifferentiated subtypes of thyroid carcinoma [16].
However, in the clinical setting, obtained DWI data is commonly analyzed using a two-dimensional region of interest in the slice of the apparent diffusion coefficient (ADC) map representing the maximum diameter of the tumor. This approach does not account for the regularly encountered heterogeneity of whole tumors and certainly does not reflect the complex micro-architectural properties of malignantly transformed tissue.
An enhanced approach using every voxel of the tumor to compute a histogram of intensity levels could help to further increase prediction of histological features of tumors by magnetic resonance imaging (MRI) [17]. This way, the magnitude of tumor heterogeneity probably is revealed in a fashion superior to the commonly used two-dimensional method [17].
To the best of the authors’ knowledge, only one study used ADC histogram analysis in thyroid cancer to differentiate benign from malignant nodules and furthermore reveal extra-thyroidal growth of papillary thyroid cancer [18]. So far, no studies demonstrated predictability of lymph node involvement by ADC histogram analysis of the primary tumor. Therefore, the primary aim of this study was to investigate the potential of ADC histogram analysis (including percentiles, entropy, skewness, and kurtosis) on data obtained with RESOLVE DWI to distinguish between nodal-negative and nodal-positive thyroid cancer. The discriminability of metastatic from non-metastatic thyroid cancer is of great clinical importance. Hence, this study investigated a promising translational approach that might have the potential to significantly increase the value of clinical-oncological imaging. The secondary aim was to correlate ADC histogram parameters with expression of important prognostic markers like p53 and Ki-67. Last, it aimed to compare our findings with the results of previous studies, which investigated the potential of DWI to predict histopathological features in thyroid cancer.

2. Results

2.1. Diffusion Weighted Imaging and Immunohistopathology of Thyroid Carcinoma

For reasons of clarity and comprehensibility, results of MRI and histopathology were organized in tables. Figure 1 shows MRI findings of a patient with follicular thyroid carcinoma, presenting as heterogeneous enlargement of the right thyroid lobe. The corresponding immunohistological images are shown in Figure 2. The calculated DWI parameters of all investigated thyroid carcinomas are summarized in Table 1 and the corresponding histopathological data is given in Table 2.

2.2. Correlation Analysis

Table 3 displays results of the correlation analysis between immunohistopathological parameters and ADC fractions as well as histogram related parameters. Correlation analysis identified the following, significant correlations: ADCmean with p53 (r = 0.548, p = 0.034), ADCmax with Ki67 (r = −0.646, p = 0.009) and p53 (r = 0.645, p = 0.009), ADCp75 with p53 (r = 0.537, p = 0.025), ADCp90 with Ki67 (r = −0.568, p = 0.027) and p53 (r = 0.588, p = 0.021), ADCmedian with p53 (r = 0.556, p = 0.032), ADCmodus with p53 (r = 0.534, p = 0.040), and kurtosis with cell count (r = −0.571, p = 0.026). Figure 3 summarizes the significant correlations graphically and displays them as dot plots.

2.3. Group Comparisons

Histogram analysis derived ADC values are compared between the nodal negative and the nodal positive group in Figure 4. Levene’s Test revealed homoscedasticity for the nodal-negative and the nodal-positive group only regarding ADCskewness (p = 0.015). For all remaining ADC derived histogram parameters, Levene’s Test showed heterogeneity of variance when comparing the nodal-negative and the nodal-positive group. Hence, group comparisons were performed using unpaired t-test for ADCskewness and Mann-Whitney-U Test for all remaining parameters. The corresponding p-values are given in Table 4. Statistically significant differences were only identified for skewness (p = 0.031) and kurtosis (p = 0.028). No other significant differences or trends were delineable when comparing thyroid carcinoma patients with restricted vs. advanced infiltration pattern (results not presented).

3. Discussion

This study aimed to investigate the potential of 3T RESOLVE DWI using an ADC histogram analysis approach to distinguish between limited and advanced thyroid cancer with reference to the status of lymphatic metastatic dissemination. To the author’s best knowledge, this work is the first to show differences in ADC histogram parameters between nodal-positive and nodal-negative thyroid cancer.
In detail, skewness and kurtosis of the ADC histograms were significantly increased in nodal-positive compared to nodal-negative thyroid cancer. This finding corresponds to previous studies in other malignant tumors, exemplarily clear cell renal cell carcinoma, and rectal cancer, which revealed that increased skewness of ADC histograms is associated with a more advanced disease stage [19,20]. Furthermore, an increase in ADC histogram skewness was observed in patients suffering from recurrent high grade glioma who showed disease progress under anti-proliferative chemotherapy, indicating ongoing proliferation of glioma cells within the tumor [21]. The association between changes in ADC values and altered cellularity in tumors is a well-known phenomenon [22]. Considering this, the findings of the aforementioned studies and our results we hypothesize that the process of lymphatic metastatic spread of thyroid cancer is linked to profound changes in the tissue microarchitecture, related to proliferation of distinct tumor cell clusters and subsequent migration via the lymphatic system, which finds its reflection in corresponding changes of the ADC histogram.
Additionally, this study found significant correlations between ADC histogram analysis derived values of thyroid cancer and corresponding immune-reactivity for p53. p53 has great importance as tumor suppressor and controls cell fate via induction of apoptosis, cell cycle arrest and senescence [23]. Under normal conditions, p53 remains undetectable for its rapid proteasomic degradation [23]. In thyroid cancer, p53 has been used as prognostic marker being associated with favorable outcome [24,25]. ADC mean, ADC max, ADC median, ADC modus, ADC p75 and ADC p90 correlated significantly with p53 expression. In general, increased ADC values of tumors have been shown to be associated with good therapeutic responses [26]. It was thereupon concluded that increased ADC values of thyroid cancer—in consent with previously published work—indicate a favorable prognosis. Furthermore, a clear inverse correlation of ADC max and ADC p90 with Ki-67 expression was identified. Ki-67 is a nuclear protein strictly associated with cell division and widely used in the clinical routine to assess proliferative activity [27].
Increased proliferation of cells, as indicated by increased expression of Ki-67, consecutively decreases the corresponding extracellular space in a given volume of tissue and thereupon reduces water diffusibility, which is reflected by decreased ADC values [22]. Thus—in accordance to other malignancies [11,28,29]—decreased ADC values are associated with an increased proliferation rate within thyroid cancer tissue.
This study furthermore identified a significant inverse correlation between cell count and kurtosis. Only few studies investigated the potential of ADC kurtosis to reflect histological properties, for example Chandarana and colleagues were able to differentiate clear cell from papillary subtype of renal cell cancer by means of ADC kurtosis [30]. It is therefore concluded that ADC histogram kurtosis provides additional insight in tumor-architectural details, but further studies are necessary to validate this finding in order to further elaborate the significance of this parameter. Conventionally, ADCmean and ADCmin were used to investigate histopathological features like cellularity of tumors in vivo [22]. However, classical ADC parameters like ADCmean and ADCmin are strongly scanner-dependent and cannot be used to compare patients investigated in different MRI devices without normalization. In contrary, histogram parameters estimate characteristics of the ADC distribution, which is not scanner-dependent like the absolute ADC values. Therefore, ADC derived histogram parameters (skewness, entropy, kurtosis) might be superior when investigating histopathological features in vivo using more than one MRI scanner in a singular study.
This study suffers from few limitations. The major limitation is the small number of patients included in this study. Furthermore, this study did not include all clinically relevant subtypes of thyroid cancer, exemplarily medullary thyroid carcinomas were not investigated. Therefore, future works including greater cohorts with different histopathological subtypes have to confirm these findings and further elucidate the relationship between histopathological findings and ADC alterations. Also, ADC histogram analysis was performed by a single, experienced reader. The suitability of histogram analysis for the clinical routine necessitates assessment of inter-reader and intra-reader variability including readers with different levels of experience. A future work needs to investigate these phenomena in a larger cohort.
ADC histogram analysis can provide more detailed information on diffusion characteristics of tumors than commonly obtained ADC parameters. For example, a previously published study demonstrated that common ADC parameters (mean, max, and min) did not reflect histopathological features like cellularity and proliferative activity in thyroid carcinoma [16]. In contrast, this study demonstrated that certain ADC histogram parameters reflect distinct histopathological features very well. Although it has proven to be a very sensitive tool for detection of microstructural changes, the specificity of ADC histogram parameters for the underlying histological changes is unclear. Characteristic changes of ADC histogram parameters in different tumor entities might be related to very different histological changes. Therefore, the significance of ADC histogram analysis should be investigated in a tumor-specific manner.

4. Materials and Methods

This retrospective study was approved (No. 2014-99) by the local research ethics committee of the Martin-Luther-University Halle-Wittenberg.

4.1. Patients

The radiological database for thyroid carcinoma was reviewed. In total, 20 patients were identified, but only 15 patients with histopathologically confirmed thyroid carcinoma had received proper DWI (using the RESOLVE sequence) and were therefore enrolled in our study.
The patient group was comprised of one male and 14 female patients. The mean age was 67 years (with a standard deviation of 12.9 years). The distribution of histopathological subtypes was as follows; follicular thyroid carcinoma: n = 4, papillary thyroid carcinoma: n = 5, anaplastic thyroid carcinoma: n = 6. Five patients were diagnosed with nodal negative thyroid cancer, and 10 patients had pathologically confirmed lymph node metastases. One patient was diagnosed with distant metastatic disease (pulmonary and pleural manifestation). Infiltration pattern ranged from restriction to the thyroid gland to advanced infiltration including infiltration of the trachea, esophagus, and internal jugular vein. An overview of demographic, clinical and pathological information is given in Table 5.

4.2. MRI

MRI of the neck was performed for all patients using a 3T device (Magnetom Skyra, Siemens, Erlangen, Germany). The imaging protocol included the following sequences:
  • axial T2 weighted (T2w) turbo spin echo (TSE) sequence (TR/TE: 4000/69, flip angle: 150°, slice thickness: 4 mm, acquisition matrix: 200 × 222, field of view: 100 mm);
  • axial T1 weighted (T1w) turbo spin echo (TSE ) sequences (TR/TE: 765/9.5, flip angle: 150°, slice thickness: 5 mm, acquisition matrix: 200 × 222, field of view: 100 mm) before and after intravenous application of contrast medium (gadopentate dimeglumine, Magnevist®, Bayer Schering Pharma, Leverkusen, Germany);
  • axial DWI (readout-segmented, multi-shot EPI sequence; TR/TE: 5400/69, flip angle 180°, slice thickness: 4 mm, acquisition matrix: 200 × 222, field of view: 100 mm) with b values of 0, 400 and 800 s/mm2. ADC maps were generated automatically by the implemented software package and analyzed as described previously [28].
All images were available in digital form and were analyzed by an experienced radiologist without knowledge of the histopathological diagnosis on a PACS workstation (Centricity PACS, GE Medical Systems, Milwaukee, WI, USA). Figure 1 shows a representative axial T2 weighted image of follicular thyroid carcinoma and corresponding axial ADC images of the whole tumor, which were used for histogram analysis (also displayed in Figure 1).

4.3. Histogram Analysis of ADC Values

DWI data was transferred in DICOM format and processed offline with a custom-made Matlab-based application (The Mathworks, Natick, MA, USA) on a standard windows operated system. The ADC maps were displayed within a graphical user interface (GUI) that enables the reader to scroll through the slices and draw a volume of interest (VOI) at the tumor’s boundary. The VOI was created by manually drawing regions of interest (ROIs) along the margin of the tumor using all slices displaying the tumor (whole lesion measure). All measures were performed by one author (AS). The ROIs were modified in the GUI and saved (in Matlab-specific format) for later processing. After setting the ROIs, the following parameters were calculated and given in a spreadsheet format: ROI volume (cm3), mean (ADCmean), maximum (ADCmax), minimum (ADCmin), median (ADCmedian), modus (ADCmodus), and the following percentils: 10th (ADCp10), 25th (ADCp25), 75th (ADCp75), and 90th (ADCp90). Additionally, histogram-based characteristics of the VOI—kurtosis, skewness, and entropy—were computed. All calculations were performed using in-build Matlab functions.

4.4. Histopathology and Immunohistochemistry

All thyroid carcinomas were surgically resected and histopathologically analysed. In every case, the proliferation index was estimated on Ki-67 antigen stained specimens using MIB-1 monoclonal antibody (DakoCytomation, Glostrup, Denmark) as reported previously [31]. Furthermore, p53 index was estimated using monoclonal antibody p57, clone DO-7 (DakoCytomation). Two high power fields (0.16 mm2 per field, ×400) were analysed. The area with the highest number of positive nuclei was selected. Figure 2 exemplarily shows Ki-67 and p53 immunostaining of a follicular thyroid carcinoma. Additionally, cellular density was calculated for each tumor as average cell count per five high power fields (×400). Furthermore, average nuclear area and total nuclear area were estimated using ImageJ package 1.48v (National Institute of Health, Bethesda, MD, USA) as described previously [11]. All histopathological sections were analysed using a research microscope Jenalumar equipped with a Diagnostic instruments camera 4.2 (Zeiss, Jena, Germany).

4.5. Statistical Analysis

Statistical analysis was performed using IBM SPSS 23™ (SPSS Inc., Chicago, IL, USA). Collected data was first evaluated by means of descriptive statistics. Correlative analysis was then performed using Spearman’s correlation coefficient in order to analyze associations between histogram analysis derived values of ADC and (immuno-) histopathological parameters. Subsequently, Levene’s Test for homogeneity of variance was performed to assess the equality of variances of ADC derived histogram parameters between different groups of thyroid carcinoma patients in order to identify the suitable test for group comparisons. In case of homoscedasticity, unpaired t test was performed to compare values among different (e.g., the metastatic and the non-metastatic) groups. In case of heteroscedasticity, Mann-Whitney-U test was performed to compare values among the different groups. Group comparisons were performed for nodal negative vs. nodal positive patients and patients with restricted (thyroid gland and trachea) vs. advanced (trachea, esophagus, jugular vein) infiltration pattern. Since only one patient with distant metastatic disease was included, a sufficient group comparison between M0 and M1 patients could not be performed. p-Values ≤ 0.05 were considered as statistically significant.

5. Conclusions

This exploratory study revealed significant differences in ADC histogram skewness and kurtosis comparing nodal negative and nodal positive thyroid cancer. Significant correlations between different ADC parameters were identified with p53, Ki-67, and cell count, substantiating the potential of ADC as an important prognostic imaging biomarker. This information certainly has the potential to aid thyroid surgeons in identifying the optimal treatment strategy for patients with thyroid cancer. Further studies investigating a greater cohort of patients are necessary to confirm these findings.

Acknowledgments

We acknowledge funding by the German Research Foundation (DFG) and University Leipzig within the program of open access publishing.

Author Contributions

Stefan Schob and Alexey Surov conceived and designed the experiments and wrote the paper; Nikolaos Pazaitis and Anne Kathrin Höhn performed the immunohistopathological experiments; Diana Horvath-Rizea and Bhogal Pervinder analyzed the data; Nikita Garnov contributed the histogram analysis tool; Hans Jonas Meyer digitalized the immunohistological slides and performed image analysis; Karl-Titus Hoffmann performed MRI; Julia Dieckow wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MDPIMultidisciplinary Digital Publishing Institute
DOAJDirectory of Open Access Journals
TLAThree Letter Acronym
LDlinear Dichroism

References

  1. Katoh, H.; Yamashita, K.; Enomoto, T.; Watanabe, M. Classification and general considerations of thyroid cancer. Ann. Clin. Pathol. 2015, 3, 1–9. [Google Scholar]
  2. Dralle, H.; Machens, A.; Basa, J.; Fatourechi, V.; Franceschi, S.; Hay, I.D.; Nikiforov, Y.E.; Pacini, F.; Pasieka, J.L.; Sherman, S.I. Follicular cell-derived thyroid cancer. Nat. Rev. Dis. Prim. 2015, 1, 15077. [Google Scholar] [CrossRef] [PubMed]
  3. Paschke, R.; Lincke, T.; Müller, S.P.; Kreissl, M.C.; Dralle, H.; Fassnacht, M. The treatment of well-differentiated thyroid carcinoma. Dtsch. Arztebl. Int. 2015, 112, 452–458. [Google Scholar] [PubMed]
  4. Nixon, I.J.; Shaha, A.R. Management of regional nodes in thyroid cancer. Oral Oncol. 2013, 49, 671–675. [Google Scholar] [CrossRef] [PubMed]
  5. Shaha, A.R. Recurrent differentiated thyroid cancer. Endocr. Pract. 2012, 18, 600–603. [Google Scholar] [CrossRef] [PubMed]
  6. Ferrari, S.M.; Fallahi, P.; Politti, U.; Materazzi, G.; Baldini, E.; Ulisse, S.; Miccoli, P.; Antonelli, A. Molecular targeted therapies of aggressive thyroid cancer. Front Endocrinol. 2015, 6, 176. [Google Scholar] [CrossRef] [PubMed]
  7. Cabanillas, M.E.; Dadu, R.; Hu, M.I.; Lu, C.; Gunn, G.B.; Grubbs, E.G.; Lai, S.Y.; Williams, M.D. Thyroid gland malignancies. Hematol. Oncol. Clin. N. Am. 2015, 29, 1123–1143. [Google Scholar] [CrossRef] [PubMed]
  8. Wendler, J.; Kroiss, M.; Gast, K.; Kreissl, M.C.; Allelein, S.; Lichtenauer, U.; Blaser, R.; Spitzweg, C.; Fassnacht, M.; Schott, M.; et al. Clinical presentation, treatment and outcome of anaplastic thyroid carcinoma: Results of a multicenter study in Germany. Eur. J. Endocrinol. 2016, 175, 521–529. [Google Scholar] [CrossRef] [PubMed]
  9. Asimakopoulos, P.; Nixon, I.J.; Shaha, A.R. Differentiated and medullary thyroid cancer: Surgical management of cervical lymph nodes. Clin. Oncol. 2017, 29, 283–289. [Google Scholar] [CrossRef] [PubMed]
  10. Mizrachi, A.; Shaha, A.R. Lymph node dissection for differentiated thyroid cancer. Mol. Imaging Radionucl. Ther. 2016, 26, 10–15. [Google Scholar] [CrossRef] [PubMed]
  11. Schob, S.; Meyer, J.; Gawlitza, M.; Frydrychowicz, C.; Müller, W.; Preuss, M.; Bure, L.; Quäschling, U.; Hoffmann, K.-T.; Surov, A. Diffusion-weighted MRI reflects proliferative activity in primary CNS lymphoma. PLoS ONE 2016, 11, e0161386. [Google Scholar] [CrossRef] [PubMed]
  12. Surov, A.; Stumpp, P.; Meyer, H.J.; Gawlitza, M.; Höhn, A.-K.; Boehm, A.; Sabri, O.; Kahn, T.; Purz, S. Simultaneous 18F-FDG-PET/MRI: Associations between diffusion, glucose metabolism and histopathological parameters in patients with head and neck squamous cell carcinoma. Oral Oncol. 2016, 58, 14–20. [Google Scholar] [CrossRef] [PubMed]
  13. Khizer, A.T.; Raza, S.; Slehria, A.-U.-R. Diffusion-weighted MR imaging and ADC mapping in differentiating benign from malignant thyroid nodules. J. Coll. Physicians Surg. Pak. 2015, 25, 785–788. [Google Scholar] [PubMed]
  14. Lu, Y.; Moreira, A.L.; Hatzoglou, V.; Stambuk, H.E.; Gonen, M.; Mazaheri, Y.; Deasy, J.O.; Shaha, A.R.; Tuttle, R.M.; Shukla-Dave, A. Using diffusion-weighted MRI to predict aggressive histological features in papillary thyroid carcinoma: A novel tool for pre-operative risk stratification in thyroid cancer. Thyroid 2015, 25, 672–680. [Google Scholar] [CrossRef] [PubMed]
  15. Porter, D.A.; Heidemann, R.M. High resolution diffusion-weighted imaging using readout-segmented echo-planar imaging, parallel imaging and a two-dimensional navigator-based reacquisition. Magn. Reson. Med. 2009, 62, 468–475. [Google Scholar] [CrossRef] [PubMed]
  16. Schob, S.; Voigt, P.; Bure, L.; Meyer, H.J.; Wickenhauser, C.; Behrmann, C.; Höhn, A.; Kachel, P.; Dralle, H.; Hoffmann, K.-T.; Surov, A. Diffusion-weighted imaging using a readout-segmented, multishot EPI sequence at 3T distinguishes between morphologically differentiated and undifferentiated subtypes of thyroid carcinoma—A preliminary study. Transl. Oncol. 2016, 9, 403–410. [Google Scholar] [CrossRef] [PubMed]
  17. Just, N. Improving tumour heterogeneity MRI assessment with histograms. Br. J. Cancer 2014, 111, 2205–2213. [Google Scholar] [CrossRef] [PubMed]
  18. Hao, Y.; Pan, C.; Chen, W.; Li, T.; Zhu, W.; Qi, J. Differentiation between malignant and benign thyroid nodules and stratification of papillary thyroid cancer with aggressive histological features: Whole-lesion diffusion-weighted imaging histogram analysis. J. Magn. Reson. Imaging 2016, 44, 1546–1555. [Google Scholar] [CrossRef] [PubMed]
  19. Kierans, A.S.; Rusinek, H.; Lee, A.; Shaikh, M.B.; Triolo, M.; Huang, W.C.; Chandarana, H. Textural differences in apparent diffusion coefficient between low- and high-stage clear cell renal cell carcinoma. Am. J. Roentgenol. 2014, 203, W637–W644. [Google Scholar] [CrossRef] [PubMed]
  20. Liu, L.; Liu, Y.; Xu, L.; Li, Z.; Lv, H.; Dong, N.; Li, W.; Yang, Z.; Wang, Z.; Jin, E. Application of texture analysis based on apparent diffusion coefficient maps in discriminating different stages of rectal cancer. J. Magn. Reson. Imaging 2016. [Google Scholar] [CrossRef] [PubMed]
  21. Nowosielski, M.; Recheis, W.; Goebel, G.; Güler, O.; Tinkhauser, G.; Kostron, H.; Schocke, M.; Gotwald, T.; Stockhammer, G.; Hutterer, M. ADC histograms predict response to anti-angiogenic therapy in patients with recurrent high-grade glioma. Neuroradiology 2011, 53, 291–302. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, L.; Liu, M.; Bao, J.; Xia, Y.; Zhang, J.; Zhang, L.; Huang, X.; Wang, J. The correlation between apparent diffusion coefficient and tumor cellularity in patients: A meta-analysis. PLoS ONE 2013, 8, e79008. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, Z.; Sun, Y. Targeting p53 for novel anticancer therapy. Transl. Oncol. 2010, 3, 1–12. [Google Scholar] [CrossRef] [PubMed]
  24. Godballe, C.; Asschenfeldt, P.; Jørgensen, K.E.; Bastholt, L.; Clausen, P.P.; Hansen, T.P.; Hansen, O.; Bentzen, S.M. Prognostic factors in papillary and follicular thyroid carcinomas: P53 expression is a significant indicator of prognosis. Laryngoscope 1998, 108, 243–249. [Google Scholar] [CrossRef] [PubMed]
  25. Bachmann, K.; Pawliska, D.; Kaifi, J.; Schurr, P.; Zörb, J.; Mann, O.; Kahl, H.J.; Izbicki, J.R.; Strate, T. P53 is an independent prognostic factor for survival in thyroid cancer. Anticancer Res. 2006, 27, 3993–3997. [Google Scholar]
  26. Padhani, A.R.; Liu, G.; Mu-Koh, D.; Chenevert, T.L.; Thoeny, H.C.; Takahara, T.; Dzik-Jurasz, A.; Ross, B.D.; Van Cauteren, M.; Collins, D.; et al. Diffusion-weighted magnetic resonance imaging as a cancer biomarker: Consensus and recommendations. Neoplasia 2009, 11, 102–125. [Google Scholar] [CrossRef] [PubMed]
  27. Schlüter, C.; Duchrow, M.; Wohlenberg, C. The cell proliferation-associated antigen of antibody Ki-67: A very large, ubiquitous nuclear protein with numerous repeated elements, representing a new kind of cell cycle-maintaining proteins. J. Cell Biol. 1993, 123, 1–10. [Google Scholar] [CrossRef]
  28. Surov, A.; Caysa, H.; Wienke, A.; Spielmann, R.P.; Fiedler, E. Correlation between different ADC fractions, cell count, Ki-67, total nucleic areas and average nucleic areas in meningothelial meningiomas. Anticancer Res. 2015, 35, 6841–6846. [Google Scholar] [PubMed]
  29. Chen, L.; Zhang, J.; Chen, Y.; Wang, W.; Zhou, X.; Yan, X.; Wang, J. Relationship between apparent diffusion coefficient and tumour cellularity in lung cancer. PLoS ONE 2014, 9, e99865. [Google Scholar] [CrossRef] [PubMed]
  30. Chandarana, H.; Rosenkrantz, A.B.; Mussi, T.C.; Kim, S.; Ahmad, A.A.; Raj, S.D.; McMenamy, J.; Melamed, J.; Babb, J.S.; Kiefer, B.; et al. Histogram analysis of whole-lesion enhancement in differentiating clear cell from papillary subtype of renal cell cancer. Radiology 2012, 265, 790–798. [Google Scholar] [CrossRef] [PubMed]
  31. Surov, A.; Gottschling, S.; Mawrin, C.; Prell, J.; Spielmann, R.P.; Wienke, A.; Fiedler, E. Diffusion-weighted imaging in meningioma: Prediction of tumor grade and association with histopathological parameters. Transl. Oncol. 2015, 8, 517–523. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Imaging findings in a patient with follicular thyroid carcinoma. (A) Magnetic resonance imaging (T2w axial section) showing a massive inhomogenous enlargement of the right thyroid lobe; (BE) represent the apparent diffusion coefficient (ADC) maps of the tumor; (F) is the ADC histogram of the whole lesion. The calculated ADC parameters (×10−5 mm2·s−1) are as follows: ADCmin = 18.2; ADCmean = 113.3; ADCmax = 315.0, mode = 114.4, ADCmedian = 108.1, P10 = 58.2, P25 = 83.2, P75 = 138.7, and P90 = 176.6. Histogram based parameters are as follows: skewness = 0.59, kurtosis = 3.88, and entropy = 3.21. The z-axis in Figure 1F gives the voxel count.
Figure 1. Imaging findings in a patient with follicular thyroid carcinoma. (A) Magnetic resonance imaging (T2w axial section) showing a massive inhomogenous enlargement of the right thyroid lobe; (BE) represent the apparent diffusion coefficient (ADC) maps of the tumor; (F) is the ADC histogram of the whole lesion. The calculated ADC parameters (×10−5 mm2·s−1) are as follows: ADCmin = 18.2; ADCmean = 113.3; ADCmax = 315.0, mode = 114.4, ADCmedian = 108.1, P10 = 58.2, P25 = 83.2, P75 = 138.7, and P90 = 176.6. Histogram based parameters are as follows: skewness = 0.59, kurtosis = 3.88, and entropy = 3.21. The z-axis in Figure 1F gives the voxel count.
Ijms 18 00821 g001
Figure 2. Immunohistochemistry of follicular thyroid carcinoma. (A) Shows Ki-67 staining (cell count: 1407, Ki-67 immunoreactiviy: 11%) and (B) shows p53 staining (cell count: 1811, p53 immunoreactivity: 36%) of the tumor displayed in Figure 1.
Figure 2. Immunohistochemistry of follicular thyroid carcinoma. (A) Shows Ki-67 staining (cell count: 1407, Ki-67 immunoreactiviy: 11%) and (B) shows p53 staining (cell count: 1811, p53 immunoreactivity: 36%) of the tumor displayed in Figure 1.
Ijms 18 00821 g002
Figure 3. Graphic summary of the significant correlations between imaging and immunohistological findings. R2-values for the plots shown in Figure 3 are as follows; (A) ADCmean & p53: r2 = 0.438; (B) ADCmax & p53: r2 = 0.425; (C) ADCmax & Ki-67: r2 = 0.464; (D) ADCp75 & p53: r2 = 0.499; (E) ADCp90 & p53: r2 = 0.431; (F) ADCp90 & Ki-67: r2 = 0.360; (G) ADCmedian & p53: r2 = 0.440; (H) ADCmodus & p53: r2 = 0.377; (I) ADCkurtosis & cell count: r2 = 0.160.
Figure 3. Graphic summary of the significant correlations between imaging and immunohistological findings. R2-values for the plots shown in Figure 3 are as follows; (A) ADCmean & p53: r2 = 0.438; (B) ADCmax & p53: r2 = 0.425; (C) ADCmax & Ki-67: r2 = 0.464; (D) ADCp75 & p53: r2 = 0.499; (E) ADCp90 & p53: r2 = 0.431; (F) ADCp90 & Ki-67: r2 = 0.360; (G) ADCmedian & p53: r2 = 0.440; (H) ADCmodus & p53: r2 = 0.377; (I) ADCkurtosis & cell count: r2 = 0.160.
Ijms 18 00821 g003aIjms 18 00821 g003b
Figure 4. Graphically summarizes the differences in histogram parameters between nodal negative and nodal positive patients with thyroid carcinoma. (A) Shows significantly increased ADC histogram skewness in noda-positive compared to nodal-negative patients; (B) demonstrates significantly increased values of ADC histogram kurtosis in nodal-positive compared to nodal negative thyroid carcinomas.
Figure 4. Graphically summarizes the differences in histogram parameters between nodal negative and nodal positive patients with thyroid carcinoma. (A) Shows significantly increased ADC histogram skewness in noda-positive compared to nodal-negative patients; (B) demonstrates significantly increased values of ADC histogram kurtosis in nodal-positive compared to nodal negative thyroid carcinomas.
Ijms 18 00821 g004
Table 1. Diffusion weighted imaging and related histogram parameters of thyroid carcinoma based on n = 15 patients.
Table 1. Diffusion weighted imaging and related histogram parameters of thyroid carcinoma based on n = 15 patients.
DWI Related ParametersMedianRangeMinimum–Maximum
ADCmean, ×10−5 mm2·s−1124.309073–163
ADCmin, ×10−5 mm2·s−114.90530.2–53
ADCmax, ×10−5 mm2·s−1250.70179147–325
P10 ADC, ×10−5 mm2·s−172.108530–114
P25 ADC, ×10−5 mm2·s−191.908452–136
P75 ADC, ×10−5 mm2·s−1140.409993–192
P90 ADC, ×10−5 mm2·s−1172.8211697–213
Median ADC, ×10−5 mm2·s−1118.009471–165
Mode ADC, ×10−5 mm2·s−1101.408853–141
Kurtosis3.641.902.89–4.79
Skewness0.301.79−0.97–0.81
Entropy3.271.982.75–4.72
Table 2. Estimated immunohistopathological parameters of thyroid carcinoma (n = 15).
Table 2. Estimated immunohistopathological parameters of thyroid carcinoma (n = 15).
ParametersMedianRangeMinimum–Maximum
Cell count, n14071808439–2247
Ki 67, %32.0909–99
p53, %4.0940–94
Total nuclear area, µm271,735148,62014,649–163,269
Average nuclear area, µm253.06133–94
Table 3. Results of Spearman’s rank order correlation analysis between DWI and immunohistological parameters (n = 15).
Table 3. Results of Spearman’s rank order correlation analysis between DWI and immunohistological parameters (n = 15).
ADC Parameters and Histogram ValuesCell Countp53Ki-67Total Nuclear AreaAverage Nuclear Area
ADCmean, ×10−3 mm2·s−1r = 0.429r = 0.548r = −0.325r = 0.389r = 0.034
p = 0.111p = 0.034p = 0.237p = 0.152p = 0.904
ADCmin, ×10−3 mm2·s−1r = 0.256r = 0.244r = −0.241r = 0.163r = −0.208
p = 0.358p = 0.381p =0.386p = 0.562p = 0.456
ADCmax, ×10−3 mm2 s−1r = 0.372r = 0.645r = −0.646r = 0.461r = 0.155
p = 0.173p = 0.009p = 0.009p = 0.084p = 0.580
ADC p10, ×10−3 mm2·s−1r = 0.361r = 0.409r = 0.289r = 0.275r = −0.079
p = 0.187p =0.130p = 0.296p = 0.321p = 0.781
ADC p25, ×10−3 mm2·s−1r = 0.375r = 0.509r = 0.361r = 0.311r = −0.064
p = 0.168p = 0.053p = 0.187p = 0.260p = 0.820
ADC p75, ×10−3 mm2·s−1r = 0.450r = 0.537r = −0.343r = 0.411r = 0.055
p = 0.092p = 0.025p =0.211p =0.128p = 0.845
ADC p90, ×10−3 mm2·s−1r = 0.289r = 0.588r = −0.568r = 0.300r = 0.075
p = 0.296p = 0.021p = 0.027p = 0.277p = 0.790
Median ADC, ×10−3 mm2·s−1r = 0.414r = 0.556r = −0.314r = 0.361r = −0.020
p = 0.125p = 0.032p = 0.254p = 0.187p = 0.945
Mode ADC, ×10−3 mm2·s−1r = 0.496r = 0.534r = −0.357r = 0.432r = −0.149
p = 0.060p = 0.040p = 0.191p = 0.108p = 0.682
Kurtosisr = −0.571r = −0.262r = −0.314r = −0.411r = −0.182
p = 0.026p =0.346p = 0.254p = 0.128p = 0.516
Skewnessr = −0.229r = −0.004r = −0.389r = 0.011r = 0.186
p = 0.413p = 0.990p = 0.152p = 0.970p = 0.507
Entropyr = 0.243r = −0.240r = 0.289r = 0.225r = 0.316
p = 0.383p = 0.389p = 0.296p = 0.420p = 0.251
Table 4. Group comparison of ADC and histogram parameters of thyroid carcinomas with (N1/2, n = 10 patients) and without lymphatic metastatic dissemination (N0, n = 5 patients).
Table 4. Group comparison of ADC and histogram parameters of thyroid carcinomas with (N1/2, n = 10 patients) and without lymphatic metastatic dissemination (N0, n = 5 patients).
ADC Parameters and Histogram ValuesN0 Mean ± SDN1/2 Mean ± SDGroup Comparison: p-Values
ADCmean, ×10−5 mm2·s−1125.2534.1111.4125.000.513
ADCmin, ×10−5 mm2·s−128.2617.3014.0216.900.075
ADCmax, ×10−5 mm2·s−1238.4469.40259.4338.500.768
P10 ADC, ×10−5 mm2·s−182.1526.1769.1423.500.371
P25 ADC, ×10−5 mm2·s−1102.2530.0089.1923.300.440
P75 ADC, ×10−5 mm2·s−1147.2639.14131.7526.430.440
P90 ADC, ×10−5 mm2·s−1170.6944.15156.5528.500.440
Median ADC, ×10−5 mm2·s−1124.1434.86109.1925.500.513
Mode ADC, ×10−5 mm2·s−1112.3225.56101.3927.500.594
Kurtosis3.230.293.810.570.028
Skewness−0.120.640.410.210.031
Entropy3.560.663.50.710.768
Table 5. Demographic and pathological data of the investigated thyroid carcinoma patients.
Table 5. Demographic and pathological data of the investigated thyroid carcinoma patients.
CaseAgeGenderHistological SubtypeInfiltration PatternM StageN Stage
191femaleanaplastictrachea01
260femalepapillarytrachea01
373malepapillarytrachea, esophagus01
468femalepapillarytrachea, esophagus internal jugular vein00
573femalepapillarytrachea01
667femaleanaplasticTrachea internal jugular vein12
773femaleanaplastictrachea, esophagus00
841femalefolliculartrachea01
972femaleanaplasticnone01
1059femaleanaplastictrachea01
1183femalepapillarytrachea00
1277femalefolliculartrachea01
1352femaleanaplastictrachea00
1451femalefolliculartrachea00
1566femaleanaplastictrachea01

Share and Cite

MDPI and ACS Style

Schob, S.; Meyer, H.J.; Dieckow, J.; Pervinder, B.; Pazaitis, N.; Höhn, A.K.; Garnov, N.; Horvath-Rizea, D.; Hoffmann, K.-T.; Surov, A. Histogram Analysis of Diffusion Weighted Imaging at 3T is Useful for Prediction of Lymphatic Metastatic Spread, Proliferative Activity, and Cellularity in Thyroid Cancer. Int. J. Mol. Sci. 2017, 18, 821. https://doi.org/10.3390/ijms18040821

AMA Style

Schob S, Meyer HJ, Dieckow J, Pervinder B, Pazaitis N, Höhn AK, Garnov N, Horvath-Rizea D, Hoffmann K-T, Surov A. Histogram Analysis of Diffusion Weighted Imaging at 3T is Useful for Prediction of Lymphatic Metastatic Spread, Proliferative Activity, and Cellularity in Thyroid Cancer. International Journal of Molecular Sciences. 2017; 18(4):821. https://doi.org/10.3390/ijms18040821

Chicago/Turabian Style

Schob, Stefan, Hans Jonas Meyer, Julia Dieckow, Bhogal Pervinder, Nikolaos Pazaitis, Anne Kathrin Höhn, Nikita Garnov, Diana Horvath-Rizea, Karl-Titus Hoffmann, and Alexey Surov. 2017. "Histogram Analysis of Diffusion Weighted Imaging at 3T is Useful for Prediction of Lymphatic Metastatic Spread, Proliferative Activity, and Cellularity in Thyroid Cancer" International Journal of Molecular Sciences 18, no. 4: 821. https://doi.org/10.3390/ijms18040821

APA Style

Schob, S., Meyer, H. J., Dieckow, J., Pervinder, B., Pazaitis, N., Höhn, A. K., Garnov, N., Horvath-Rizea, D., Hoffmann, K. -T., & Surov, A. (2017). Histogram Analysis of Diffusion Weighted Imaging at 3T is Useful for Prediction of Lymphatic Metastatic Spread, Proliferative Activity, and Cellularity in Thyroid Cancer. International Journal of Molecular Sciences, 18(4), 821. https://doi.org/10.3390/ijms18040821

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop