- 1Department of Pediatrics, Kawasaki Disease Center, Kaohsiung Chang Gung Memorial Hospital, College of Medicine, Chang Gung University, Kaohsiung, Taiwan
- 2Department of Computer Science and Information Engineering, Tamkang University, New Taipei City, Taiwan
- 3Department of Information Management, Chang Gung University, Kaohsiung, Taiwan
- 4Department of E-Sport Technology Management, Cheng Shiu University, Kaohsiung, Taiwan
- 5Department of Medicine, College of Medicine, China Medical University, Taichung, Taiwan
- 6Department of Pediatric Cardiology, China Medical University Children's Hospital, China Medical University, Taichung, Taiwan
- 7Center of Structure and Congenital Heart Disease/Ultrasound and Department of Cardiology, Children's Hospital, China Medical University, Taichung, Taiwan
Introduction: Kawasaki disease (KD) may increase the risk of myocardial infarction or sudden death. In children, delayed KD diagnosis and treatment can increase coronary lesions (CLs) incidence by 25% and mortality by approximately 1%. This study focuses on the use of deep learning algorithm-based KD detection from cardiac ultrasound images.
Methods: Specifically, object detection for the identification of coronary artery dilatation and brightness of left and right coronary artery is proposed and different AI algorithms were compared. In infants and young children, a dilated coronary artery is only 1-2 mm in diameter than a normal one, and its ultrasound images demonstrate a large amount of noise background-this can be a considerable challenge for image recognition. This study proposes a framework, named Scaled-YOLOv4-HarDNet, integrating the recent Scaled-YOLOv4 but with the CSPDarkNet backbone replaced by the CSPHarDNet framework.
Results: The experimental result demonstrated that the mean average precision (mAP) of Scaled-YOLOv4-HarDNet was 72.63%, higher than that of Scaled YOLOv4 and YOLOv5 (70.05% and 69.79% respectively). In addition, it could detect small objects significantly better than Scaled-YOLOv4 and YOLOv5.
Conclusions: Scaled-YOLOv4-HarDNet may aid physicians in detecting KD and determining the treatment approach. Because relatively few artificial intelligence solutions about images for KD detection have been reported thus far, this paper is expected to make a substantial academic and clinical contribution.
1. Introduction
Kawasaki disease (KD), a systemic vasculitis predominantly affecting medium-sized arteries (1, 2). It may accelerate coronary arteriosclerosis and sudden death risks, and it is the leading cause of acquired heart disease in children of many developed countries (3). The global incidence of Kawasaki Disease (KD) has increased in the last 10–20 years, with an annual incidence rates per 100,000 children <5 years old around 100–300 in Japan and northern Asia regions (4). It is therefore a significant disease burden for children.
KD-associated coronary lesions (CLs) are the most serious cardiovascular sequelae and can lead to acute intra-coronary thrombosis and stenosis (5). These lethal complications often occur in KD patients with delayed use of intravaenous immnunosglssosbulin (IVIG) or refractory to initial IVIG (6). Abou Sherif et al. (7) demonstrated the effects of coronary artery aneurysms on human health. Therefore, timely assessment of patients with KD using imaging modalities is the key factor ensuring the effective lowering of long-term cardiovascular events in KD patients.
Among multi-modality imaging used to evaluate KD-associated CLs, echocardiography is the most common to identify patients with KD as it is non-invasive, easy-accessible, widely available, and cost-effective and provides a real-time qualitative assessment of the coronary system (8). At present time, echocardiography has to be assessed by cardiologist. The echocardiography contains large amount of information. It is therefore logic to employ the advanced bioinformatics technology to assess the echocardiography in patients with KD (9, 10).
Artificial intelligence (AI), deep learning (DL), is basically a number of neural network-related algorithms to identify patterns with purpose and has been widely applied in the medical field to assist in diagnosis. Its power comes from its ability to find these associations from large amounts of data and draw non-linear relationships between various predictors and an outcome of interest without background knowledge. There are numerous applications of deep learning in cardiac ultrasound image classification (11–14) described briefly in Table 1. In addition, some important object detection algorithms belonged to the category of DL were recently proposed (21–23) organized in Table 2. However, to the best of our knowledge, there are no relevant studies regarding applied AI for CLs assessment in patients with KD.
2. Method
There are five major procedures of this study depicted in Figure 1. First of all, this retrospective study included cardiac ultrasound images from patients with KD in Kaohsiung Chang Gung Hospital from June 1, 2000, to June 30, 2020 (IRB approving No. 202001238B0C502). The extracted data were patient age, body height, body weight, cardiovascular diameter of coronary arteries, and cardiac ultrasound images in the DICOM format. The inclusion criterion was the receipt of KD diagnosis (International Classification of Diseases, Ninth Revision, Clinical Modification code: 446.1, or International Classification of Diseases, Tenth Revision code: M303).
All the echocardiography was performed by board-certificated pediatric cardiologists. The echocardiographic images of coronary arteries were obtained according to the standard procedure as shown in the scientific statement of American Heart Association on Kawasaki Disease (10). The inter-observer variability was around 0.90–0.95 by laboratory quality check.
After anonymization of the obtained echo image data. As such, the privacy of the included patients was not infringed upon. We labeled region of interest on the images using LabelMe. The annotation results was validated by another experienced cardiologist. Considering the retrospective nature of this study, we did not consider additional behavioral factors of the included children. Thereafter, in-depth learning was performed using different models and evaluate which object detection algorithm will perform better.
2.1. Research definition
Coronary arteries can be divided into left coronary arteries (LCAs) and right coronary arteries (RCAs). One of the main symptoms of KD is CLs, which are detectable on cardiac ultrasound. Clinically, short axis ultrasound can be used to assess for abnormalities in LCAs and RCAs. Figure 2 presents ultrasound images from Kaohsiung Chang Gung Hospital; we took two features of CLs, brightness and dilatation (10), to be detected.
Figure 2. Two major echocardiographic coronary patterns in children with KD. (A) LCA brightness, (B) LCA dilatation, (C) RCA brightness, and (D) RCA dilatation.
We created object detection algorithms to learn the precise location of the two major symptoms. Image annotation was required to mark the location of the aneurysm as shown in Figure 3. Subsequent algorithms were designed to read the contents of the JSON materials, including file location, disease attribute, and corresponding location. In addition, the association between data amplification method and marking results, such as the image rotation-marking frame relationship, was considered so as to maximize the advantages of data amplification.
Figure 3. Object detection picture and annotation format in JSON. (A) Left coronary artery annotation and (B) JSON file.
2.2. Scaled-YOLOv4-HarDNet object detection algorithm
The main framework of Scaled-YOLOv4-HarDNet is depicted in Figure 4. Because the proposed algorithm adopts the main features of Scaled-YOLOv4 and HarDNet, we illustrate the original settings and also explain the key modification items of Scaled-YOLOv4. Scaled-YOLOv4 integrates various methods, such as the cross stage partial network (CSPNet), spatial pyramid pooling with CSP (SPPCSP), FPN, path aggregation network (PANet), mish activation function, label smoothing, and complete IOU (CIoU) loss function.
The use of CSPHarDNet as the backbone is one of the main features of the proposed algorithm, and the input of CSPHarDNet is CBM, which is formed through a combination of convolution layers and the batch normalization and Mish activation function. The input resolution of the first convolution layer is 416 × 416. The Mish function is a self-regularized nonmonotonic neural activation function, which increases the penetration of information in the neural network. Moreover, ZCRn comprises zero pool, CBM and CSPRn. CSPRn refers to the CSPNet framework, with n denoting the number of copies. CSPNet divides the function scheme into two parts: In the first part, gradient change is retained in the feature map from the beginning to the end to ensure the accuracy and reduce the amount of calculation and storage cost required. The second part involves the skip connection approach of ResNet. The first part is consistent with the feature map of the second part. The output resolutions of three output heads are 52 × 52, 26 × 26, and 13 × 13, respectively.
Scaled-YOLOv4 uses FPN and PANet methods in the Neck area, whereas YOLOv3 only uses FPN. FPN performs upsampling from smaller to larger resolutions and then connects with the larger ZCRn. The PANet framework applies bottom-up path enhancement to the previous local convolution layer by performing the upsampling operation. Its advantage includes shortening the information path between high- and low-resolution features.
Originating from YOLOv5 (22), Scaled-YOLOv4 adopts CSPNet in Neck to increase the training and detection speed further. The activation function in Neck can be replaced by the Mish function. Scaled-YOLOv4 is applicable in different environment based on depth requirements. For instance, the P5 architecture can be adopted in an edge computing environment, and the P6 architecture can be adopted in general environments. For larger GPU devices, such as the Tesla V100, the P7 architecture can be used for training and prediction. Thus, we compared our algorithm with YOLOv5 by using the P5 architecture.
Scaled-YOLOv4 uses a head area identical to that of YOLOv3 and outputs the same three dimensions with 76 × 76 (256 feature maps), 38 × 38 (512 feature maps), and 19 × 19 (1,024 feature maps) resolution when the input resolution is 608. When the resolution is reduced to 416 × 416, the output resolution of the three dimensions are 13 × 13 (1,024 feature maps), 26 × 26 (512 feature maps), and 52 × 52 (256 feature maps). The only change in the head area is the loss function, and the loss function for Scaled-YOLOv4 is CIoU, which measures the difference between the real and the prediction boxes.
Finally, in our proposed algorithm Scaled-YOLOv4-HarDNet, there are four variants which are the combination of the activation function and the number of layers. Firstly, we substitute the original activation function ReLU and to be Mish function. Secondly, because the number of layer in a deep learning algorithm also influences the performance, we select the number of layers in the CSPHarDNet includes 39 and 68. Instead of using the full name of Scaled-YOLOv4-HarDNet, the four variants are named HarDNet39-ReLU, HarDNet39-Mish, HarDNet68-ReLU, and HarDNet68-Mish in short. This study will conduct extensive experiments and compare the proposed algorithms with two existing benchmark algorithms.
3. Experimental results
In total, 1,283 images–of which 395 and 472 were on the LCA for brightness and dilatation, respectively; 174 and 242 were on the RCA for brightness and dilatation, respectively–were included in the validation dataset. The number of images varied due to the coronary artery symptoms of KD tending to occur in the LCA and RCA. After image arrangement, we use 70% of images to train our proposed algorithms, 20% of images are as the validation dataset during the training, and the rest of 10% of images as the test dataset. Two well-known benchmark algorithms, Scaled-YOLOv4 and YOLOv5, are applied into the comparison. The following training parameters: batch = 72 and epoch = 300 are used across the experiments. To carry out these experiments, we employed the PyTorch environment with Nvidia NGC container 21.08 on a Tesla V100 server of Taiwan Computing Cloud.
After each algorithm was trained, we used the rest of 10% of images to test the real performance of these algorithm. The test results are presented in Table 3. We list the performance of each algorithm on the four patterns of LCA and RCA via the mean average precision (mAP) metric in percentage. The best mAP value of each class is shown in bold. The overall column shows the average results of the four patterns for each compared algorithm. Among them, HarDNet68-Mish achieves the highest overall mAP score to be 72.63% and the HarDNet39-ReLU (71.08%) is the inferior version of our proposed method. However, HarDNet39-ReLU is still better than Scaled-YOLOv4 (70.05%) and YOLOv5 (69.79%). The models including HarDNet had significantly improved outcomes compared with the original Scaled YOLOv4; moreover, Mish function also led to a better performance than ReLU.
From the four image types, the average result of LCA dilatation is >83.32% which shows the LCA dilatation is easier to be detected. However, RCA dilatation is not as high as the one of LCA dilatation. Moreover, the pattern of LCA brightness has the lowest mAP in average which implies LCA brightness is the most difficult one to be detected.
The effectiveness of the aforementioned methods was compared using different Intersections over Union (IoU) as criteria; the comparison of AP50, APsmall, APmedium, ARsmall, and ARmedium is presented in Table 4. Here, AP represents the average precision or accuracy, AR the average recall, and AP50 the AP with IoU = 50% and “small” and “medium” represent the sizes of the bounding boxes; we did not include a large classification because no large-sized objects were considered. The best AP value of each class is also shown in bold. The result of the AP50 is quite similar to the one of mAP-50 because HarDNet68-Mish is best algorithm which outperform the Scaled-YOLOv4 and YOLOv5. When it comes to the metrics of APsmall or ARsmall, HarDNet68-Mish is also the best one compared with others. However, if we apply the metrics of APmedium and ARmedium, Scaled-YOLOv4 is better than the proposed algorithms.
4. Discussion
This pioneer, pilot study has demonstrated that several AI based algorithms showed encouraging result. The assessment of coronary morphology could be performed by AI is promising and showed that assessment of coronary morphology may be assisted by AI algorithm in the future and may thus greatly decrease the work-loading by cardiologist using conventional assessment of coronary morphology for echocardiography in patients with KD.
In the proposed Scaled-YOLOv4-HarDNet, we verified the 39 and 68 HarDNet layers and substituted the ReLU function with the Mish function. We the impact of the number of layers improves the solution quality which is expected. However, it does not guarantee to yield better prediction quality when it comes to 85 layers in our pilot experiments. Then, Mish function generally enhances the performance when it is compared to the ReLU function compared to its original setting. It is noticeable that the performance of HarDNet39-Mish is almost equal to the one of HarDNet68-ReLU. As a result, if the object detection model is designed to be installed on an edge device, we could consider the HarDNet39-Mish because the model size is smaller.
When we compared the proposed algorithm with Scaled-YOLOv4 and YOLOv5, even though HarDNet39-ReLU is the inferior version of the proposed method, it is better than the two benchmark algorithms. Then, HarDNet68-Mish yielded the optimal overall mAP (72.63%) across the four features compared with Scaled-YOLOv4 and YOLOv5 (70.05 and 69.79%, respectively).
The four patterns also have different difficulty levels. LCA dilatation has the better detection result from no matter which object detection algorithm is applied. However, the performance of the brightness pattern of LCA and RCA is poor. The primary reason is that the characteristic of dilatation shows that the color inside the vessel is black and white outside the vessel. However, the brightness pattern is the grayish white of the vessel wall compared to the grayish color in the background. As a result, the brightness pattern is harder to detect.
COCO metrics include average precision and average recall for bounding boxes of different sizes. Here, the proposed Scaled-YOLOv4-HarDNet, with 68 layers and the Mish function, significantly outperformed Scaled-YOLOv4 and YOLOv5 in terms of the average precision and average recall for small bounding boxes. However, for the medium-sized bounding box, Scaled-YOLOv4, but not YOLOv5, performed <0.66% better than did the proposed framework. Thus, our COCO metrics results indicated that proposed algorithm can detect small objects extremely well and medium-sized objects adequately. Finally, because we did not encounter any large objects in the medical images in this study, how the performance of our proposed algorithm for such objects compares with the performance of the two benchmark algorithms remains unknown.
According to these results, the strength of the proposed algorithms are to detect the two KD patterns on echocardiographic images, and also perform objective assistance for cardiologist. Finally, varies of echo equipment models, ultrasound resolution, case numbers, and different results from different cardiovascular specialist in different time may be the limitation of this study. Due to the medical images are available from the Philip EPIQ 7C and iE33 belonged to the studied hospital, it limits the capability of reading the echo images of other manufacturer due to the color and resolution, frame rate are different. On the other hand, when there are more cardiovascular specialists are involved in the research, it is beneficial to validate the labeling results in many rounds and yield results precisely. The results from this study still need further validation from different researches before a conclusion.
5. Conclusions
KD is the leading cause of acquired heart disease in children especial in Asia. Precisely diagnosis and treatment is the major point to prevent coronary lesions in KD. We proposed the Scaled-YOLOv4-HarDNet, which integrates HarDNet as the backbone rather than DarkNet. In the four variants of the proposed algorithm with different number of layers and activation functions, they are better than the two benchmark algorithms Scaled-YOLOv4 and YOLOv5. Moreover, the proposed algorithm is good at the smaller object size of coronary artery in children which is suitable to detect the lesion in LCA and RCA.
The results from this study will help clinician to identify coronary artery lesion objectively from the help of AI and further support the detection of KD.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
H-CK provides the ultrasound images for investigation and supervised the study. S-HC is the paper's guarantor. Y-HC contributed to conception and design of the study. Y-CL annotated the medical images. Y-CW wrote the code and ran the experiments. C-YC and T-DW performed the statistical analysis. I-HT performed most echocardiography imaging in the investigation process, provided insight into coronary arteritis, and critically revised the manuscript. K-SH supervised image labeling corrections. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This project was supported by the National Science and Technology Council with ID (NSTC 109-2221-E-032-042 and NSTC 109-2221-E-230-006) and Kaohsiung Chang Gung Memorial Hospital (Project Nos. CMRPD3M0011, CORPG8L0271, and CORPG8M0091).
Acknowledgments
We thank the two reviewers who gave constructive comments for this paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Darby JB, Jackson JM. Kawasaki disease and multisystem inflammatory syndrome in children: an overview and comparison. Am Family Phys. (2021) 104:244–252. Available online at: https://www.aafp.org/pubs/afp/issues/2021/0900/p244.html
3. Burns JC, Glodé MP. Kawasaki syndrome. Lancet. (2004) 364:533–44. doi: 10.1016/S0140-6736(04)16814-1
4. Ying-Hsien KMLSCHJHYMHL Huang, Ho-Chang K. Increased incidence of kawasaki disease in taiwan in recent years: a 15 years nationwide population-based cohort study. Front Pediatr. (2019) 7:121. doi: 10.3389/fped.2019.00121
5. Dallaire F, Dahdah N. New equations and a critical appraisal of coronary artery Z scores in healthy children. J Am Soc Echocardiogr. (2011) 24:60–74. doi: 10.1016/j.echo.2010.10.004
6. Suda K, Iemura M, Nishiono H, Teramachi Y, Koteda Y, Kishimoto S, et al. Long-term prognosis of patients with Kawasaki disease complicated by giant coronary aneurysms: a single-institution experience. Circulation. (2011) 123:1836–42. doi: 10.1161/CIRCULATIONAHA.110.978213
7. Abou Sherif S, Ozden Tok O, Taşköylü Ö, Goktekin O, Kilic ID. Coronary artery aneurysms: a review of the epidemiology, pathophysiology, diagnosis, and treatment. Front Cardiovasc Med. (2017) 4:24. doi: 10.3389/fcvm.2017.00024
8. Tai I, Hsieh KS. Imaging of coronary artery abnormalities in kawasaki disease with emphasis on multimodalities. In: Kawasaki Disease. Springer (2022). p. 63–81.
9. Fukazawa R, Kobayashi J, Ayusawa M, Hamada H, Miura M, Mitani Y, et al. JCS/JSCS 2020 guideline on diagnosis and management of cardiovascular sequelae in Kawasaki disease. Circ J. (2020) 84:1348–407. doi: 10.1253/circj.CJ-19-1094
10. McCrindle BW, Rowley AH, Newburger JW, Burns JC, Bolger AF, Gewitz M, et al. Diagnosis, treatment, and long-term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association. Circulation. (2017) 135:e927–99. doi: 10.1161/CIR.0000000000000484
11. Gao X, Li W, Loomes M, Wang L. A fused deep learning architecture for viewpoint classification of echocardiography. Information Fusion. (2017) 36:103–13. doi: 10.1016/j.inffus.2016.11.007
12. Chen H, Zheng Y., Park J. -H., Heng P. -A., Zhou S. K. (2016). “Iterative multi-domain regularized deep learning for anatomical structure detection and segmentation from ultrasound images,” in International Conference on Medical Image Computing and Computer-Assisted Intervention (Springer), 487–496.
13. Gungor D,. G., Rao B., Wolverton C., Guracar I. (2020). View Classification and Object Detection in Cardiac Ultrasound to Localize Valves via Deep Learning. Available online at: https://openreview.net/forum?id=eLblfUbJ05
14. Chen SH, Wang CW, Tai I, Weng KP, Chen YH, Hsieh KS. Modified YOLOv4-DenseNet algorithm for detection of ventricular septal defects in ultrasound images. Int J Interact Multimedia Artif Intell. (2021). 6:101–8. doi: 10.9781/ijimai.2021.06.001
15. Nascimento JC, Carneiro G. Multi-atlas segmentation using manifold learning with deep belief networks. In: Biomedical Imaging (ISBI), 2016 IEEE 13th International Symposium on. Prague: IEEE (2016). p. 867–71.
16. Bridge CP, Ioannou C, Noble JA. Automated annotation and quantitative description of ultrasound videos of the fetal heart. Med Image Anal. (2017) 36:147–61. doi: 10.1016/j.media.2016.11.006
17. Poudel R. K., Lamata P., Montana G. (2017). “Recurrent fully convolutional neural networks for multi-slice MRI cardiac segmentation,” in International Workshop on Reconstruction and Analysis of Moving Body Organs, International Workshop on Whole-Heart and Great Vessel Segmentation from 3D Cardiovascular MRI in Congenital Heart Disease (Springer), 83–94.
18. Avendi M, Kheradvar A, Jafarkhani H. A combined deep-learning and deformable-model approach to fully automatic segmentation of the left ventricle in cardiac MRI. Med Image Anal. (2016) 30:108–19. doi: 10.1016/j.media.2016.01.005
19. Wolterink JM, Leiner T, de Vos BD, van Hamersvelt RW, Viergever MA, Išgum I. Automatic coronary artery calcium scoring in cardiac CT angiography using paired convolutional neural networks. Med Image Anal. (2016) 34:123–36. doi: 10.1016/j.media.2016.04.004
20. Sirjani N, Moradi S, Oghli MG, Hosseinsabet A, Alizadehasl M Azin adn Yadollahi, Shiri I, et al. Automatic cardiac evaluations using a deep video object segmentation network. Insights Into Imaging. (2022) 69:13. doi: 10.1186/s13244-022-01212-9
21. Bochkovskiy A, Wang CY, Liao HYM. YOLOv4: optimal speed and accuracy of object detection. arXiv preprint arXiv:200410934. (2020). doi: 10.48550/arXiv.2004.10934
22. Jocher NKMTVR G,. YOLOv5. (2020). Available online at: https://github.com/ultralytics/yolov5
23. Wang CY, Bochkovskiy A, Liao HYM. Scaled-yolov4: scaling cross stage partial network. In: Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition. Nashville, TN: IEEE (2021). p. 13029–38.
24. Redmon J, Divvala S, Girshick R, Farhadi A. You only look once: unified, real-time object detection. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. Las Vegas, NV: IEEE (2016). p. 779–88.
25. Liu W, Anguelov D., Erhan D., Szegedy C., Reed S., Fu C. -Y., et al. (2016). “Ssd: Single shot multibox detector,” in European Conference on Computer Vision (Springer), 21–37.
26. Redmon J, Farhadi A. YOLO9000: better, faster, stronger. In: 2017. IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, HI: IEEE (2017). p. 6517–25.
27. Lin TY, Goyal P, Girshick R, He K, Dollár P. Focal loss for dense object detection. IEEE Transactions on Pattern Analysis and Machine Intelligence. Venice: IEEE (2018).
28. Redmon J, Farhadi A. Yolov3: an incremental improvement. arXiv preprint arXiv:180402767. (2018). doi: 10.48550/arXiv.1804.02767
Keywords: Kawasaki disease, echocardiography, deep learning, object detection, Scaled-YOLOv4, HarDNet, coronary dilatation and brightness
Citation: Kuo H-C, Chen S-H, Chen Y-H, Lin Y-C, Chang C-Y, Wu Y-C, Wang T-D, Chang L-S, Tai I-H and Hsieh K-S (2023) Detection of coronary lesions in Kawasaki disease by Scaled-YOLOv4 with HarDNet backbone. Front. Cardiovasc. Med. 9:1000374. doi: 10.3389/fcvm.2022.1000374
Received: 22 July 2022; Accepted: 30 December 2022;
Published: 20 January 2023.
Edited by:
Tommaso Gori, Johannes Gutenberg University Mainz, GermanyReviewed by:
Fuyong Jiao, Xi'an Jiaotong University, ChinaRakesh Kumar Pilania, Post Graduate Institute of Medical Education and Research (PGIMER), India
Copyright © 2023 Kuo, Chen, Chen, Lin, Chang, Wu, Wang, Chang, Tai and Hsieh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shih-Hsin Chen, c2hjaGVuJiN4MDAwNDA7bWFpbC50a3UuZWR1LnR3; Kai-Sheng Hsieh, a3Noc2llaCYjeDAwMDQwO2hvdG1haWwuY29t