Next Article in Journal
‘I Didn’t Even Recognise Myself’: Survivors’ Experiences of Altered Appearance and Body Image Distress during and after Treatment for Head and Neck Cancer
Previous Article in Journal
Chromosomal Translocations in NK-Cell Lymphomas Originate from Inter-Chromosomal Contacts of Active rDNA Clusters Possessing Hot Spots of DSBs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deep Learning Fast Screening Approach on Cytological Whole Slides for Thyroid Cancer Diagnosis

1
Department of Pathology, Tri-Service General Hospital, Taipei 11490, Taiwan
2
Institute of Pathology and Parasitology, National Defense Medical Center, Taipei 11490, Taiwan
3
Graduate Institute of Biomedical Engineering, National Taiwan University of Science and Technology, Taipei 106335, Taiwan
4
Graduate Institute of Applied Science and Technology, National Taiwan University of Science and Technology, Taipei 106335, Taiwan
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(15), 3891; https://doi.org/10.3390/cancers13153891
Submission received: 25 June 2021 / Revised: 29 July 2021 / Accepted: 30 July 2021 / Published: 2 August 2021
(This article belongs to the Section Methods and Technologies Development)

Abstract

:

Simple Summary

Papillary thyroid carcinoma is the most common type of thyroid cancer and could be cured if diagnosed and treated early. In clinical practice, the primary method for determining diagnosis of papillary thyroid carcinoma is manual visual inspection of cytopathology slides, which is difficult, time consuming and subjective with a high inter-observer variability and sometimes causes suboptimal patient management due to false-positive and false-negative results. This study presents a fast, fully automatic and efficient deep learning framework for fast screening of cytological slides for thyroid cancer diagnosis. We confirmed the robustness and effectiveness of the proposed method based on evaluation results from two different types of slides: thyroid fine needle aspiration smears and ThinPrep slides.

Abstract

Thyroid cancer is the most common cancer in the endocrine system, and papillary thyroid carcinoma (PTC) is the most prevalent type of thyroid cancer, accounting for 70 to 80% of all thyroid cancer cases. In clinical practice, visual inspection of cytopathological slides is an essential initial method used by the pathologist to diagnose PTC. Manual visual assessment of the whole slide images is difficult, time consuming, and subjective, with a high inter-observer variability, which can sometimes lead to suboptimal patient management due to false-positive and false-negative. In this study, we present a fully automatic, efficient, and fast deep learning framework for fast screening of papanicolaou-stained thyroid fine needle aspiration (FNA) and ThinPrep (TP) cytological slides. To the authors’ best of knowledge, this work is the first study to build an automated deep learning framework for identification of PTC from both FNA and TP slides. The proposed deep learning framework is evaluated on a dataset of 131 WSIs, and the results show that the proposed method achieves an accuracy of 99%, precision of 85%, recall of 94% and F1-score of 87% in segmentation of PTC in FNA slides and an accuracy of 99%, precision of 97%, recall of 98%, F1-score of 98%, and Jaccard-Index of 96% in TP slides. In addition, the proposed method significantly outperforms the two state-of-the-art deep learning methods, i.e., U-Net and SegNet, in terms of accuracy, recall, F1-score, and Jaccard-Index ( p < 0.001 ). Furthermore, for run-time analysis, the proposed fast screening method takes 0.4 min to process a WSI and is 7.8 times faster than U-Net and 9.1 times faster than SegNet, respectively.

1. Introduction

Thyroid cancer is the most prevalent cancer in the endocrine system and accounts for the majority of head and neck cancer cases [1]. Thyroid cancer has been on the rise globally for the past two decades, including in the United States, despite a decline in the incidence of certain other cancer forms [1]. Thyroid cancer is three times more prevalent in women than in men. The types of thyroid cancer are papillary carcinoma, follicular carcinoma, Hürthle (oncocytic) cell carcinoma, poorly differentiated carcinoma, medullary carcinoma and anaplastic (undifferentiated) carcinoma [2]. Papillary thyroid carcinoma (PTC) is the most common type of thyroid carcinoma, which accounts for 70% to 80% of all thyroid malignancies. The prognosis of PTC is better than other types of thyroid carcinoma [3]. Thyroid fine needle aspiration (FNA) is an important, safe tool for diagnosing PTC with an accuracy of approximately 94%, and a high degree of sensitivity, specificity [4,5]. Thyroid FNA is applied to distinguish benign from neoplastic or malignant thyroid nodules [6]. The wide use of thyroid FNA has greatly reduced the unnecessary thyroid surgical intervention and thus increased the percent of malignant nodules among all nodules surgically removed. The Bethesda System for Reporting Thyroid Cytopathology (TBSRTC) [7] is the universally accepted reporting system for thyroid FNA diagnosis. Although the cytologic feature of PTC is well documented, including enlarged overlapping nuclei, irregular nuclear contours, intranuclear pseudoinclusions, nuclear grooving, and fine, pale chromatin [8,9], traditionally a time-consuming cytologic analysis is performed by an experienced pathologist who manually examines the glass slides under a light microscope. The most common stain for cytological preparations is the Papanicolaou stain. May-Grünwald Giemsa Stain is one of the common Romanwsky stains used in cytology. Digital pathology has emerged as a potential new standard of care, in which glass slides are digitized into whole slide images (WSIs) using digital slide scanners. With over 100 million pixels in a typical WSI, pathologists find it difficult to manually identify all the information in histopathological images [10]. Thus, an automated diagnosis methods based on artificial intelligence are developed to overcome the constraints of manual and complex diagnosis process [11,12]. In recent years, deep learning has emerged as a potential approach for the automated analysis of medical images. Automating the diagnostic process helps the pathologists to make correct diagnosis in a short period of time. Deep learning has been commonly used to identify diseases such as retinal disease [13], skin cancer [14], colorectal polyp [15], cardiac arrhythmia [16], neurological problems [17], psychiatric problems [18], acute intracranial hemorrhage [19], and autism [20]. Deep learning has also shown the ability to help pathologists diagnose, classify, and segment cancer. For example, Courtiol et al. [21] trained a deep convolutional neural network (MesoNet) to automatically and accurately estimate the overall survival in mesothelioma patients from diagnostic unannotated histopathology images. Yamamoto et al. [22] trained a deep learning based framework that can derive explainable features from diagnostic unannotated histopathology images and anticipate predictions more accurately than humans. Zhang et al. [23] proposed a deep learning-based framework for automating the human-like diagnostic reasoning process, which would include second opinions and thereby encourage clinic consensus. Sanyal et al. [24] trained a convolutional neural network to classify PTC and non PTC on microphotographs from thyroid fine needle aspiration cytology (FNAC). In comparison, Sanyal et al.’s method [24] obtains the diagnostic accuracy of 85.06% on microphotographs of size 512 × 512 pixels from thyroid FNAC while the proposed method achieves an accuracy of 99% on gigapixels WSI of papanicolaou-stained thyroid FNA and ThinPrep (TP) cytological slides for detection and segmentation of PTC. Furthermore, Sanyal et al.’s method [24] can only operate on FNA slides while the proposed method performs consistently well on both thyroid FNA and TP slides. Ke et al. [25] trained a deep convolutional neural network (Faster R-CNN [26]) for detection of PTC from ultrasonic images. To the best of the authors’ knowledge, this is the first study to build an automated deep learning framework for detection and segmentation of PTC from papanicolaou-stained thyroid FNA and TP cytological slides. Figure 1 presents the proposed framework structure and the dataset information. Figure 1a presents the workflow of the system from collection of data to analysis of outcome. Figure 1(ai) shows the thyroid smears are obtained through FNA and TP; in Figure 1(aii), slides of thyroid FNA and TP are prepared with papanicolaou’s staining; in Figure 1(aiii), stained slides are digitalized at 20× objective magnification using Leica AT Turbo scanner; in Figure 1(aiv), digitized whole slide gigapixel images are distributed into a separate training (21%) set and a testing (79%) set; in Figure 1(av), WSIs are processed with fast background filtering of the proposed system; in Figure 1(avi), cytological samples of PTC of individual WSIs are rapidly identified by the proposed deep learning model in seconds. Figure 1b presents the distribution thyroid FNA and TP cytological slides for training and testing. Figure 1c presents the distribution of the number of tiles per WSI. Figure 1d presents the size distribution of the WSIs w.r.t. the width and height. In evaluation, as this is the first study on automatic segmentation of PTC in papanicolaou-stained thyroid FNA and TP cytological slides, we compare the proposed method with the two state-of-the-art deep learning models, including U-Net [27] and SegNet [28].

2. Materials and Methods

2.1. The Dataset

De-identified and digitized 131 WSIs, including 120 PTC cytologic slides (smear, papanicolaou-stained, n = 120) and 11 PTC cytologic slides (TP, papanicolaou-stained, n = 11) were obtained from the Department of Pathology, Tri-Service General Hospital, Taipei, Taiwan. All papillary thyroid carcinoma smears were cytologically diagnosed and histologically confirmed by the two expert pathologists. The well-preserved thyroid FNAs, which were done within the last two years, are selected. Ethical approvals have been obtained from the research ethics committee of the Tri-Service General Hospital (TSGHIRB No.1-107-05-171 and No.B202005070), and the data were de-identified and used for a retrospective study without impacting patient care. All the stained slides were scanned using Leica AT Turbo (Leica, Germany), at 20× objective magnification. The average slide dimensions are 77,338 × 37,285 pixels with physical size 51.13 × 23.21 mm 2 . The ground truth annotations were produced by two expert pathologists. The training model uses a total of 28 papanicolaou-stained WSIs (21%), including 25 thyroid FNA and 3 TP cytologic slides. The remaining 103 papanicolaou-stained WISs (79%), including 95 thyroid FNA and 8 TP cytologic slides, are used as a separate testing set for evaluation. The detailed information on the distribution of data could be found in Figure 1b.

2.2. Methods

In this work, we propose a fast and efficient deep learning based framework for segmentation of PTC in papanicolaou-stained thyroid FNA and TP cytological slides. Figure 2a presents the workflow of the proposed framework. Initially, each WSI is formatted into hierarchical tile-based data structure and assessed by the proposed deep learning model to produce the segmentation results of PTC in papanicolaou-stained thyroid FNA and TP WSIs. Figure 2b shows the detailed architecture of the proposed deep learning model.

2.2.1. Whole Slide Image Processing

Figure 2a presents the framework for segmentation of PTC from papanicolaou-stained thyroid FNA and TP WSIs. Initially, each WSI is formatted into hierarchical data structure and processed by fast background filtering to efficiently discard all the background and reduce the amount of computation per slide and the proposed deep learning model is used to produce the segmentation results of PTC on papanicolaou-stained thyroid FNA and TP WSIs. The details of the proposed WSI processing framework is described in Section S1.1 of the Supplementary Methods and Evaluation Metrics.

2.2.2. Proposed Convolution Network Architecture

The proposed deep learning network is built using VGG16 model as a backbone and adapted from a fully convolutional network framework [29], which has been widely employed in the field of pathology such as neuropathology [30], histopathology [31], and microscopy [32]. The proposed deep learning network consists of a padding layer, six convolutional layers, five max-pooling layers, two dropout layers, one deconvolutional layer, and a cropping layer (see Section S1.2 in Supplementary Methods and Evaluation Metrics for details). The detailed architecture of the proposed deep learning network is shown in Table 1 and Figure 2c.

2.2.3. Implementation details

The proposed method uses the VGG16 model as the backbone for training, with the network optimized using stochastic gradient descent (SGD) optimization and the cross entropy function as a loss function. Furthermore, the network training parameters of the proposed method, including the learning rate, dropout ratio, and weight decay, are set as 1 × 10 10 , 0.5, and 0.0005, respectively. The benchmark methods (U-Net and SegNet) are implemented using the keras implementation [33]. For training, the benchmark methods (U-Net and SegNet) are initialized using a pre-trained VGG16 model, and the network is optimized using Ada-delta optimization with the cross entropy function as a loss function. Furthermore, the network training parameters of the benchmark methods, including the learning rate, dropout ratio, and the weight decay are set to 0.0001, 0.2, and 0.0002, respectively. The proposed method and the benchmark methods (U-Net and SegNet) uses the same framework to process a WSI, which is described in Section 2.2.

3. Results

3.1. Evaluation Metrics

The quantitative evaluation is produced using five measurements, i.e., accuracy, precision, recall, F1-score, and Jaccard-Index. The evaluation metrics are described in Section S2 of Supplementary Methods and Evaluation Metrics.

3.2. Quantitative Evaluation with Statistical Analysis

The aim of this study is to develop a deep learning framework that can automatically detect PTC from both papanicolaou-stained thyroid FNA and TP cytological slides. For quantitative evaluation, we compared the performance of the proposed method with the state-of-the-art deep learning models, including U-Net and SegNet. Table 2 shows the quantitative evaluation results for segmentation of PTC from papanicolaou-stained thyroid FNA and TP cytological slides. The experimental results show that overall the proposed method achieves the highest accuracy 99%, precision 86%, recall 94%, F1-score 88% and Jaccard 82%, and outperforms the two benchmark approaches. For TP slides, the proposed method obtains even better results with accuracy 99%, precision 97%, recall 98%, F1-score 98% and Jaccard 96%. Figure 3 presents the box plots of the quantitative evaluation results, showing that (a) the proposed method works constantly well overall and significantly outperforms the benchmark methods in terms of accuracy, recall, F1-score, and Jaccard-Index ( p < 0.001 ) and (b) the type of cytological slides, i.e., FNA or TP, does not affect the performance of the proposed model, which consistently performs well for both kinds of data while the benchmark approaches tend to perform better in TP than FNA w.r.t. accuracy and precision ( p < 0.05 ).
For statistical analysis, the quantitative scores were analyzed with the Fisher’s Least Significant Difference (LSD) test using SPSS software (see Table 3). In comparison, the proposed method significantly outperforms the benchmark methods (U-Net and SegNet) in terms of accuracy, recall, F1-score, and Jaccard-Index, based on LSD test ( p < 0.001 ). The experimental results demonstrate the high accuracy, efficiency, and reliability of the proposed method on papanicolaou-stained thyroid FNA and TP cytological slides. Figure 4 compares the qualitative segmentation results of the proposed method and two benchmark methods (U-Net and SegNet) in FNA and TP WSIs, showing that the proposed method is able to segment PTC from papanicolaou-stained thyroid FNA and TP slides consistent with the reference standard while the state-of-the-art benchmark methods (U-Net and SegNet) are unable to detect the PTC in some cases. Figure 5 further shows the annotations produced by the expert pathologists with typical PTC features, including papillary like structure, elongated nucleus, pale nucleus, pseudoinclusions in nuclear cytoplasm, nucleoli and longitudinal grooves in FNA and TP slides and the results by the proposed method with typical PTC features as well.

3.3. Run Time Analysis

Due to enormous size of WSI, the computing time of WSI analysis is crucial for practical clinical use. We examined the computing time using various hardware configurations (see Table 4). Table 4 compares the computational efficiency of the proposed method with the benchmark methods (U-Net and SegNet), showing that the proposed method takes 0.4 minute to process a WSI using four GeForce GTX 1080 Ti GPUs and 1.7 minute using a single GeForce GTX 1080 Ti GPU, whereas the U-Net model takes 13.2 minutes and the SegNet model takes 15.4 minutes. In addition, even with a single low-cost GPU, the proposed method outperforms the benchmark approaches with less computing time and is 7.8 times faster than U-Net and 9.1 times faster than SegNet. Overall, the proposed method is shown to be capable of detecting PTC reliably in both FNA and TP WSIs and rapidly in seconds, making it suitable for practical clinical use.

4. Discussion

In this study, we present a fully automatic and efficient deep learning framework for segmentation of PTC from both papanicolaou-stained thyroid FNA and TP cytological slides. PTC is the most common form of the thyroid cancer with best prognosis and most patients can be cured if treated appropriately and early enough. Thyroid FNA, in addition to pathological examination, is considered the most effective approach for the clinical diagnosis of PTC due to its diagnostic safety, minimal invasiveness and high accuracy. Manual pathological diagnosis is sometimes difficult, time-consuming, and laborious task. In cytopathological diagnosis, pathologists have to conduct a thorough inspection of all information on the glass slides under a light microscope. In recent years, digital pathology in which glass slides are converted into WSIs, has emerged as a potential new standard of care, allowing pathological images to be examined using computer-based algorithms. A typical WSI comprises of more than 100 million pixels, which makes it difficult for pathologists to manually conduct a thorough inspection of all information on cytopathological and histopathological slides. Algorithm-assisted pathologists on WSI diagnosis revealed higher accuracy than either the algorithm or the pathologist alone on review of lymph nodes for metastatic breast cancer, especially improved the sensitivity of detection for micrometastases (91% vs. 83%, p = 0.02) [34]. Pathologists easily over-calculate the percent of tumor cells, and the use of AI-based analysis increases the accuracy in applying tumor cell count to genetic analysis [35]. Although artificial intelligence has the potential to provide advantages in accuracy, precision, and efficiency through the automation of digital pathology. Artificial intelligence-related applications are also facing challenges, including regulatory roadblocks, quality of data, interpretability, algorithm validation, reimbursement, and clinical adoption [36]. There are compelling reasons to believe that digital pathology in addition to artificial intelligence for cytological PTC diagnosis is a viable answer to this problem because it aids in the production of more accurate diagnoses, shortens examination times, reduces the pathologists’ workload, and lowers examination cost. For diagnosis of WSIs, many current algorithms are not well adapted to clinical applications due to the high computational cost of employing computational techniques. For practical clinical usage, we develop a fast, efficient, and fully automatic deep learning framework for fast screening of both FNA and TP slides. The experimental results show that the proposed deep learning framework has been demonstrated to be effective, as the proposed method achieves the accuracy and recall of over 90%. Furthermore, the proposed method achieves significantly superior performance than the state-of-the art deep learning models, including U-Net and SegNet using Fisher’s LSD test ( p < 0.001 ) . The results demonstrate that the proposed method is able to segment PTC with high accuracy, precision, and sensitivity, comparable to the referenced standard produced by pathologists in seconds laying the groundwork for the use of computational decision support systems in clinical practice. The including 131 cytologic slides are all definitive papillary thyroid carcinoma under TBSRTC criteria. The limitation of our study is that in practice, cytopathologists have to face quantitatively insufficient specimens for a definitive diagnosis, for instance, the diagnosis of “suspicious for PTC”. In future work, the proposed framework could be extended to the detection and segmentation of “suspicious for PTC” slides and different types of carcinoma, such as follicular carcinoma, Hürthle (oncocytic) cell carcinoma, medullary carcinoma, poorly differentiated carcinoma and anaplastic carcinoma.

5. Conclusions

In this work, we introduce a deep learning-based framework for automatic detection and segmentation of PTC in papanicolaou-stained thyroid FNA and TP cytological slides. We evaluated the proposed framework on a dataset of 131 WSIs, including 120 PTC cytologic slides (smear, papanicolaou-stained, n = 120) and 11 PTC cytologic slides (TP, papanicolaou-stained, n = 11), and the experimental results show that the proposed method achieves high accuracy, precision, recall, F1-score, and Jaccard-Index. In addition, we compared the proposed method with the state-of-the-art deep learning models, including U-Net and SegNet. Based on Fisher’s LSD test, the proposed method significantly outperforms the two benchmark methods in terms of accuracy, recall, F1-score, and Jaccard-Index ( p < 0.001 ) .

Supplementary Materials

The following are available at https://www.mdpi.com/article/10.3390/cancers13153891/s1, Supplementary Methods and Evaluation Metrics.

Author Contributions

C.-W.W. and T.-K.C. conceived the idea of this work. C.-W.W. designed the methodology and the software of this work. M.-A.K., Y.-C.L. and D.-Z.H. carried out the validation of the methodology of this work. M.-A.K. performed the formal analysis of this work. Y.-C.L. performed the investigation. J.-J.W. and Y.-J.L. participated in annotation or curation of the dataset. C.-W.W., M.-A.K., J.-J.W. and Y.-J.L. prepared and wrote the manuscript. C.-W.W. reviewed and revised the manuscript. M.-A.K. prepared the visualization of the manuscript. C.-W.W. supervised this work. C.-W.W. and T.-K.C. administered this work and also acquired funding for this work. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the Ministry of Science and Technology of Taiwan, under a grant (MOST109-2221-E-011-018-MY3) and Tri-Service General Hospital, Taipei, Taiwan (TSGH-C-108086, TSGH-D-109094 and TSGH-D-110036).

Institutional Review Board Statement

In this study, to enable the development of recognizing-PTC DL algorithms, de-identified, digitized 131 WSIs, including 119 PTC cytologic slides (smear, papanicolaou-stained, n = 120) and 11 PTC cytologic slides (TP, papanicolaou-stained, n = 11) were obtained from the Department of Pathology, Tri-Service General Hospital, National Defense Medical Center, Taipei, Taiwan. Ethical approvals have been obtained from the research ethics committee of the Tri-Service General Hospital (TSGHIRB No.1-107-05-171 and No.B202005070).

Informed Consent Statement

Patient consent was formally waived by the approving review board, and the data were de-identified and used for a retrospective study without impacting patient care.

Data Availability Statement

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

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 145–164. [Google Scholar] [CrossRef] [Green Version]
  2. Bai, Y.; Kakudo, K.; Jung, C.K. Updates in the pathologic classification of thyroid neoplasms: A review of the world health organization classification. Endocrinol. Metab. 2020, 35, 696. [Google Scholar] [CrossRef] [PubMed]
  3. Limaiem, F.; Rehman, A.; Mazzoni, T. Cancer, Papillary Thyroid Carcinoma (PTC); StatPearls: Treasure Island, FL, USA, 2020. [Google Scholar]
  4. Kini, S.; Miller, J.; Hamburger, J.; Smith, M. Cytopathology of papillary carcinoma of the thyroid by fine needle aspiration. Acta Cytol. 1980, 24, 511–521. [Google Scholar]
  5. Akhtar, M.; Ali, M.A.; Huq, M.; Bakry, M. Fxine-needle aspiration biopsy of papillary thyroid carcinoma: Cytologic, histologic, and ultrastructural correlations. Diagn. Cytopathol. 1991, 7, 373–379. [Google Scholar] [CrossRef]
  6. Sidawy, M.K.; Vecchio, D.M.D.; Knoll, S.M. Fine-needle aspiration of thyroid nodules: Correlation between cytology and histology and evaluation of discrepant cases. Cancer Cytopathol. Interdiscip. Int. J. Am. Cancer Soc. 1997, 81, 253–259. [Google Scholar] [CrossRef]
  7. Cibas, E.S.; Ali, S.Z. The 2017 Bethesda system for reporting thyroid cytopathology. Thyroid 2017, 27, 1341–1346. [Google Scholar] [CrossRef]
  8. Kini, S. Guides to clinical aspiration biopsy: Thyroid. In Thyroid Adequacy Criteria and Thinprep; Igaku-Shoin: New York, NY, USA, 1996. [Google Scholar]
  9. Kumar, S.; Singh, N.; Siddaraju, N. “Cellular swirls” and similar structures on fine needle aspiration cytology as diagnostic clues to papillary thyroid carcinoma: A report of 4 cases. Acta Cytol. 2010, 54, 939–942. [Google Scholar]
  10. Mukhopadhyay, S.; Feldman, M.D.; Abels, E.; Ashfaq, R.; Beltaifa, S.; Cacciabeve, N.G.; Cathro, H.P.; Cheng, L.; Cooper, K.; Dickey, G.E.; et al. Whole slide imaging versus microscopy for primary diagnosis in surgical pathology: A multicenter blinded randomized noninferiority study of 1992 cases (pivotal study). Am. J. Surg. Pathol. 2018, 42, 39. [Google Scholar] [CrossRef]
  11. Yu, K.H.; Beam, A.L.; Kohane, I.S. Artificial intelligence in healthcare. Nat. Biomed. Eng. 2018, 2, 719–731. [Google Scholar] [CrossRef]
  12. Syrykh, C.; Abreu, A.; Amara, N.; Siegfried, A.; Maisongrosse, V.; Frenois, F.X.; Martin, L.; Rossi, C.; Laurent, C.; Brousset, P. Accurate diagnosis of lymphoma on whole-slide histopathology images using deep learning. NPJ Digit. Med. 2020, 3, 1–8. [Google Scholar] [CrossRef]
  13. De Fauw, J.; Ledsam, J.R.; Romera-Paredes, B.; Nikolov, S.; Tomasev, N.; Blackwell, S.; Askham, H.; Glorot, X.; O’Donoghue, B.; Visentin, D.; et al. Clinically applicable deep learning for diagnosis and referral in retinal disease. Nat. Med. 2018, 24, 1342–1350. [Google Scholar] [CrossRef]
  14. Esteva, A.; Kuprel, B.; Novoa, R.A.; Ko, J.; Swetter, S.M.; Blau, H.M.; Thrun, S. Correction: Corrigendum: Dermatologist-level classification of skin cancer with deep neural networks. Nature 2017, 546, 686. [Google Scholar] [CrossRef]
  15. Chen, P.J.; Lin, M.C.; Lai, M.J.; Lin, J.C.; Lu, H.H.S.; Tseng, V.S. Accurate classification of diminutive colorectal polyps using computer-aided analysis. Gastroenterology 2018, 154, 568–575. [Google Scholar] [CrossRef]
  16. Hannun, A.Y.; Rajpurkar, P.; Haghpanahi, M.; Tison, G.H.; Bourn, C.; Turakhia, M.P.; Ng, A.Y. Cardiologist-level arrhythmia detection and classification in ambulatory electrocardiograms using a deep neural network. Nat. Med. 2019, 25, 65–69. [Google Scholar] [CrossRef] [PubMed]
  17. Titano, J.J.; Badgeley, M.; Schefflein, J.; Pain, M.; Su, A.; Cai, M.; Swinburne, N.; Zech, J.; Kim, J.; Bederson, J.; et al. Automated deep-neural-network surveillance of cranial images for acute neurologic events. Nat. Med. 2018, 24, 1337–1341. [Google Scholar] [CrossRef]
  18. Durstewitz, D.; Koppe, G.; Meyer-Lindenberg, A. Deep neural networks in psychiatry. Mol. Psychiatry 2019, 24, 1583–1598. [Google Scholar] [CrossRef]
  19. Lee, H.; Yune, S.; Mansouri, M.; Kim, M.; Tajmir, S.H.; Guerrier, C.E.; Ebert, S.A.; Pomerantz, S.R.; Romero, J.M.; Kamalian, S.; et al. An explainable deep-learning algorithm for the detection of acute intracranial haemorrhage from small datasets. Nat. Biomed. Eng. 2019, 3, 173–182. [Google Scholar] [CrossRef]
  20. Hazlett, H.C.; Gu, H.; Munsell, B.C.; Kim, S.H.; Styner, M.; Wolff, J.J.; Elison, J.T.; Swanson, M.R.; Zhu, H.; Botteron, K.N.; et al. Early brain development in infants at high risk for autism spectrum disorder. Nature 2017, 542, 348–351. [Google Scholar] [CrossRef]
  21. Courtiol, P.; Maussion, C.; Moarii, M.; Pronier, E.; Pilcer, S.; Sefta, M.; Manceron, P.; Toldo, S.; Zaslavskiy, M.; Le Stang, N.; et al. Deep learning-based classification of mesothelioma improves prediction of patient outcome. Nat. Med. 2019, 25, 1519–1525. [Google Scholar] [CrossRef]
  22. Yamamoto, Y.; Tsuzuki, T.; Akatsuka, J.; Ueki, M.; Morikawa, H.; Numata, Y.; Takahara, T.; Tsuyuki, T.; Tsutsumi, K.; Nakazawa, R.; et al. Automated acquisition of explainable knowledge from unannotated histopathology images. Nat. Commun. 2019, 10, 1–9. [Google Scholar] [CrossRef] [Green Version]
  23. Zhang, Z.; Chen, P.; McGough, M.; Xing, F.; Wang, C.; Bui, M.; Xie, Y.; Sapkota, M.; Cui, L.; Dhillon, J.; et al. Pathologist-level interpretable whole-slide cancer diagnosis with deep learning. Nat. Mach. Intell. 2019, 1, 236–245. [Google Scholar] [CrossRef]
  24. Sanyal, P.; Mukherjee, T.; Barui, S.; Das, A.; Gangopadhyay, P. Artificial intelligence in cytopathology: A neural network to identify papillary carcinoma on thyroid fine-needle aspiration cytology smears. J. Pathol. Inform. 2018, 9, 43. [Google Scholar] [PubMed]
  25. Ke, W.; Wang, Y.; Wan, P.; Liu, W.; Li, H. An Ultrasonic Image Recognition Method for Papillary Thyroid Carcinoma Based on Depth Convolution Neural Network. In International Conference on Neural Information Processing; Springer: Berlin/Heidelberg, Germany, 2017; pp. 82–91. [Google Scholar]
  26. Ren, S.; He, K.; Girshick, R.; Sun, J. Faster R-CNN: Towards Real-Time Object Detection with Region Proposal Networks. IEEE Trans. Pattern Anal. Mach. Intell. 2017, 39, 1137–1149. [Google Scholar] [CrossRef] [Green Version]
  27. Falk, T.; Mai, D.; Bensch, R.; Çiçek, Ö.; Abdulkadir, A.; Marrakchi, Y.; Böhm, A.; Deubner, J.; Jäckel, Z.; Seiwald, K.; et al. U-Net: Deep learning for cell counting, detection, and morphometry. Nat. Methods 2019, 16, 67–70. [Google Scholar] [CrossRef]
  28. Badrinarayanan, V.; Kendall, A.; Cipolla, R. Segnet: A deep convolutional encoder-decoder architecture for image segmentation. IEEE Trans. Pattern Anal. Mach. Intell. 2017, 39, 2481–2495. [Google Scholar] [CrossRef]
  29. Shelhamer, E.; Long, J.; Darrell, T. Fully convolutional networks for semantic segmentation. IEEE Trans. Pattern Anal. Mach. Intell. 2017, 39, 640–651. [Google Scholar] [CrossRef]
  30. Signaevsky, M.; Prastawa, M.; Farrell, K.; Tabish, N.; Baldwin, E.; Han, N.; Iida, M.A.; Koll, J.; Bryce, C.; Purohit, D.; et al. Artificial intelligence in neuropathology: Deep learning-based assessment of tauopathy. Lab. Investig. 2019, 99, 1019. [Google Scholar] [CrossRef] [PubMed]
  31. Naylor, P.; Laé, M.; Reyal, F.; Walter, T. Segmentation of nuclei in histopathology images by deep regression of the distance map. IEEE Trans. Med. Imaging 2018, 38, 448–459. [Google Scholar] [CrossRef]
  32. Zhu, R.; Sui, D.; Qin, H.; Hao, A. An extended type cell detection and counting method based on FCN. In Proceedings of the 2017 IEEE 17th International Conference on Bioinformatics and Bioengineering (BIBE), Washington, DC, USA, 23–25 October 2017; pp. 51–56. [Google Scholar]
  33. Gupta, D.; Jhunjhunu Wala, R.; Juston, M.; MC, J. Image Segmentation Keras: Implementation of Segnet, FCN, UNet, PSPNet and Other Models in Keras. Available online: https://github.com/divamgupta/image-segmentation-keras (accessed on 15 October 2020).
  34. Steiner, D.F.; MacDonald, R.; Liu, Y.; Truszkowski, P.; Hipp, J.D.; Gammage, C.; Thng, F.; Peng, L.; Stumpe, M.C. Impact of deep learning assistance on the histopathologic review of lymph nodes for metastatic breast cancer. Am. J. Surg. Pathol. 2018, 42, 1636. [Google Scholar] [CrossRef]
  35. Sakamoto, T.; Furukawa, T.; Lami, K.; Pham, H.H.N.; Uegami, W.; Kuroda, K.; Kawai, M.; Sakanashi, H.; Cooper, L.A.D.; Bychkov, A.; et al. A narrative review of digital pathology and artificial intelligence: Focusing on lung cancer. Transl. Lung Cancer Res. 2020, 9, 2255. [Google Scholar] [CrossRef] [PubMed]
  36. Bera, K.; Schalper, K.A.; Rimm, D.L.; Velcheti, V.; Madabhushi, A. Artificial intelligence in digital pathology—New tools for diagnosis and precision oncology. Nat. Rev. Clin. Oncol. 2019, 16, 703–715. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of the proposed framework structure and the dataset information. (a) The proposed framework structure. (i) Collection of thyroid smear samples of patients through FNA and TP; (ii) preparation of thyroid FNA and TP slides with papanicolaou’s staining; (iii) digitalization of cytological slides at 20× objective magnification using Leica AT Turbo scanner; (iv) distribution of digitized gigapixel WSIs into separate training set (21%) and testing set (79%); (v) processing of WSIs with fast background filtering; (vi) identification of PTC tissues of individual WSIs using the proposed deep learning model in seconds. (b) Distribution of thyroid FNA and TP cytological slides for training and testing. (c) Distribution of tile numbers per WSI. (d) Size distribution of the WSIs with width and height as blue and red, respectively.
Figure 1. Illustration of the proposed framework structure and the dataset information. (a) The proposed framework structure. (i) Collection of thyroid smear samples of patients through FNA and TP; (ii) preparation of thyroid FNA and TP slides with papanicolaou’s staining; (iii) digitalization of cytological slides at 20× objective magnification using Leica AT Turbo scanner; (iv) distribution of digitized gigapixel WSIs into separate training set (21%) and testing set (79%); (v) processing of WSIs with fast background filtering; (vi) identification of PTC tissues of individual WSIs using the proposed deep learning model in seconds. (b) Distribution of thyroid FNA and TP cytological slides for training and testing. (c) Distribution of tile numbers per WSI. (d) Size distribution of the WSIs with width and height as blue and red, respectively.
Cancers 13 03891 g001
Figure 2. The proposed deep learning framework for segmentation of PTC from papanicolaou-stained thyroid FNA and TP WSIs. (a) The workflow of the proposed deep learning framework. Initially, each WSI is formatted into hierarchical data structure and processed by fast background filtering and the proposed deep learning model is used to produce the segmentation results of PTC on papanicolaou-stained thyroid FNA and TP WSIs. (b) The detailed architecture of the proposed deep learning model.
Figure 2. The proposed deep learning framework for segmentation of PTC from papanicolaou-stained thyroid FNA and TP WSIs. (a) The workflow of the proposed deep learning framework. Initially, each WSI is formatted into hierarchical data structure and processed by fast background filtering and the proposed deep learning model is used to produce the segmentation results of PTC on papanicolaou-stained thyroid FNA and TP WSIs. (b) The detailed architecture of the proposed deep learning model.
Cancers 13 03891 g002
Figure 3. (a) The box plot of quantitative evaluation results of the proposed method and the benchmark methods for PTC segmentation. The results of LSD tests shows that the proposed method significantly outperforms the benchmark methods in terms of accuracy, recall, F1-score, and Jaccard-Index ( p < 0.001 ). (b) The box plot of quantitative evaluation results of the proposed method and the benchmark methods for PTC segmentation from papanicolaou-stained thyroid FNA and TP cytological slides. The results shows that the proposed deep learning framework performs consistently well for the segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs and the cell types of thyroid FNA and TP cytological slides do not affect the judgment of the proposed deep learning model while the benchmark methods perform inconsistent for the segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs. The box plot of quantitative evaluation results of PTC segmentation where the outliers > 1.5 × interquartile range are marked with a dot and the outlier > 3 × interquartile are marked with an asterisk.
Figure 3. (a) The box plot of quantitative evaluation results of the proposed method and the benchmark methods for PTC segmentation. The results of LSD tests shows that the proposed method significantly outperforms the benchmark methods in terms of accuracy, recall, F1-score, and Jaccard-Index ( p < 0.001 ). (b) The box plot of quantitative evaluation results of the proposed method and the benchmark methods for PTC segmentation from papanicolaou-stained thyroid FNA and TP cytological slides. The results shows that the proposed deep learning framework performs consistently well for the segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs and the cell types of thyroid FNA and TP cytological slides do not affect the judgment of the proposed deep learning model while the benchmark methods perform inconsistent for the segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs. The box plot of quantitative evaluation results of PTC segmentation where the outliers > 1.5 × interquartile range are marked with a dot and the outlier > 3 × interquartile are marked with an asterisk.
Cancers 13 03891 g003
Figure 4. Qualitative segmentation results of the proposed method and two benchmark methods (U-Net and SegNet) for segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs.
Figure 4. Qualitative segmentation results of the proposed method and two benchmark methods (U-Net and SegNet) for segmentation of PTC on papanicolaou-stained thyroid FNA and TP WSIs.
Cancers 13 03891 g004
Figure 5. Typical PTC image features on manual annotations by the expert pathologists and automatic segmentation results by the proposed method. (a) Annotations by the expert pathologists on FNA and TP WSIs with typical PTC features, including papillary structure, enlarged nucleus, pale nucleus, pseudoinclusions in nuclear cytoplasm, nucleoli and longitudinal grooves and irregular cell membranes. (b) Automatic segmentation results by the proposed method with typical PTC features.
Figure 5. Typical PTC image features on manual annotations by the expert pathologists and automatic segmentation results by the proposed method. (a) Annotations by the expert pathologists on FNA and TP WSIs with typical PTC features, including papillary structure, enlarged nucleus, pale nucleus, pseudoinclusions in nuclear cytoplasm, nucleoli and longitudinal grooves and irregular cell membranes. (b) Automatic segmentation results by the proposed method with typical PTC features.
Cancers 13 03891 g005
Table 1. The architecture of the proposed deep learning model.
Table 1. The architecture of the proposed deep learning model.
LayerFeatures (Train)Features (Inference)Kernel SizeStride
Input512 × 512 × 3512 × 512 × 3--
Padding712 × 712 × 3712 × 712 × 3--
Conv1_1 + relu1_1710 × 710 × 64710 × 710 × 643 × 31
Conv1_2 + relu1_2710 × 710 × 64710 × 710 × 643 × 31
Pool1355 × 355 × 64355 × 355 × 642 × 22
Conv2_1 + relu2_1355 × 355 × 128355 × 355 × 1283 × 31
Conv2_2 + relu2_2355 × 355 × 128355 × 355 × 1283 × 31
Pool2178 × 178 × 128178 × 178 × 1282 × 22
Conv3_1 + relu3_1178 × 178 × 256178 × 178 × 2563 × 31
Conv3_2 + relu3_2178 × 178 × 256178 × 178 × 2563 × 31
Conv3_3 + relu3_3178 × 178 × 256178 × 178 × 2563 × 31
Pool389 × 89 × 25689 × 89 × 2562 × 22
Conv4_1 + relu4_189 × 89 × 51289 × 89 × 5123 × 31
Conv4_2 + relu4_289 × 89 × 51289 × 89 × 5123 × 31
Conv4_3 + relu4_389 × 89 × 51289 × 89 × 5123 × 31
Pool445 × 45 × 51245 × 45 × 5122 × 22
Conv5_1 + relu5_145 × 45 × 51245 × 45 × 5123 × 31
Conv5_2 + relu5_245 × 45 × 51245 × 45 × 5123 × 31
Conv5_3 + relu5_345 × 45 × 51245 × 45 × 5123 × 31
Pool523 × 23 × 51223 × 23 × 5122 × 22
Conv6 + relu6 + Drop617 × 17 × 409617 × 17 × 40967 × 71
Conv7 + relu7 + Drop717 × 17 × 409617 × 17 × 40961 × 11
Conv817 × 17 × 317 × 17 × 31 × 11
Deconv9576 × 576 × 3576 × 576 × 364 × 6432
Cropping512 × 512 × 3512 × 512 × 3--
Output Class Map512 × 512 × 1512 × 512 × 1--
Table 2. Quantitative evaluation for segmentation of PTC in thyroid FNA and TP slides.
Table 2. Quantitative evaluation for segmentation of PTC in thyroid FNA and TP slides.
Proposed MethodU-Net [27]SegNet [28]
AllFNATPAllFNATPAllFNATP
Accuracy0.990.990.990.920.920.980.920.920.97
Precision0.860.850.970.740.730.870.810.800.98
Recall0.940.940.980.610.610.660.560.560.55
F1-score0.880.870.980.640.630.740.620.610.67
Jaccard-Index0.820.800.960.490.480.600.480.470.54
Table 3. Statistical analysis: multiple comparisons for segmentation of PTC.
Table 3. Statistical analysis: multiple comparisons for segmentation of PTC.
Measurement(I) Method(J) MethodMean Diff. (I–J)Std. ErrorSig.95% C.I.
Lo. BoundUp. Bound
AccuracyProposed MethodU-Net0.0784732 *0.0068 < 0.001 0.06510.0918
SegNet0.0761947 *0.0068 < 0.001 0.06290.0895
PrecisionProposed MethodU-Net0.1187516 *0.0289 < 0.001 0.06200.1755
SegNet0.04520.02890.1180−0.01160.1020
RecallProposed MethodU-Net0.3336451 *0.0271 < 0.001 0.28030.3870
SegNet0.3856975 *0.0271 < 0.001 0.33230.4391
F1-scoreProposed MethodU-Net0.2392282 *0.0266 < 0.001 0.18680.2916
SegNet0.2578479 *0.0266 < 0.001 0.20550.3102
Jaccard IndexProposed MethodU-Net0.3238887 *0.0290 < 0.001 0.26680.3809
SegNet0.3391651 *0.0290 < 0.001 0.28210.3962
* The proposed method is significantly better than the benchmark methods (U-net [27] and SegNet [28]) (p < 0.001).
Table 4. Comparison on hardware and computing time per WSI.
Table 4. Comparison on hardware and computing time per WSI.
MethodCPURAMGPUTime (min) *
Proposed MethodIntel Xeon Gold 6134 CPU @ 3.20GHz × 16128 GB4 × GeForce GTX 1080 Ti0.4
Proposed MethodIntel Xeon CPU E5-2650 v2 @ 2.60GHz × 1632 GB1 × GeForce GTX 1080 Ti1.7
U-Net [27]Intel Xeon CPU E5-2650 v2 @ 2.60GHz × 1632 GB1 × GeForce GTX 1080 Ti13.2
SegNet [28]Intel Xeon CPU E5-2650 v2 @ 2.60GHz × 1632 GB1 × GeForce GTX 1080 Ti15.4
* The size of WSI in this evaluation is 4,069,926,912 pixels (91,632 × 44,416 pixels).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, Y.-J.; Chao, T.-K.; Khalil, M.-A.; Lee, Y.-C.; Hong, D.-Z.; Wu, J.-J.; Wang, C.-W. Deep Learning Fast Screening Approach on Cytological Whole Slides for Thyroid Cancer Diagnosis. Cancers 2021, 13, 3891. https://doi.org/10.3390/cancers13153891

AMA Style

Lin Y-J, Chao T-K, Khalil M-A, Lee Y-C, Hong D-Z, Wu J-J, Wang C-W. Deep Learning Fast Screening Approach on Cytological Whole Slides for Thyroid Cancer Diagnosis. Cancers. 2021; 13(15):3891. https://doi.org/10.3390/cancers13153891

Chicago/Turabian Style

Lin, Yi-Jia, Tai-Kuang Chao, Muhammad-Adil Khalil, Yu-Ching Lee, Ding-Zhi Hong, Jia-Jhen Wu, and Ching-Wei Wang. 2021. "Deep Learning Fast Screening Approach on Cytological Whole Slides for Thyroid Cancer Diagnosis" Cancers 13, no. 15: 3891. https://doi.org/10.3390/cancers13153891

APA Style

Lin, Y. -J., Chao, T. -K., Khalil, M. -A., Lee, Y. -C., Hong, D. -Z., Wu, J. -J., & Wang, C. -W. (2021). Deep Learning Fast Screening Approach on Cytological Whole Slides for Thyroid Cancer Diagnosis. Cancers, 13(15), 3891. https://doi.org/10.3390/cancers13153891

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