- Department of Ultrasound, The First Affiliated Hospital of Shantou University Medical College, Shantou, China
Purpose: To assess the clinical utility of ultrasound grayscale ratio (UGSR) in distinguishing between benign and malignant thyroid nodules.
Methods: We conducted a retrospective analysis of patients diagnosed with thyroid nodules between January 2017 and December 2021. Malignancy and benignity were determined based on histopathology (biopsy or surgery) as the reference standard. Conventional ultrasonography (US) was performed to measure the maximum diameter of each nodule and assess positive features of the Chinese Thyroid Imaging Reporting and Data System (C-TIRADS). UGSR was calculated as the ratio of nodule grayscale value to surrounding thyroid parenchyma grayscale value, quantified using ImageJ software. Two experienced radiologists used ImageJ software for UGSR measurements. Logistic regression analysis examined the association between UGSR and thyroid malignancy. A receiver operating characteristic (ROC) curve analysis was conducted, and the area under the curve (AUC) was calculated to assess UGSR’s effectiveness in distinguishing between benign and malignant thyroid nodules. The UGSR cut-off value was established through ROC curve analysis.
Results: A total of 125 nodules (78 benign, 47 malignant) were evaluated. Significant differences were observed between benign and malignant nodules in maximum diameter (p = 0.042), UGSR (p < 0.001), and C-TIRADS features (p < 0.001), supporting UGSR’s diagnostic utility. However, no significant intergroup differences were observed in gender or age distribution (p>0.05 for both). Multivariable logistic regression analysis identified UGSR, irregular margins, taller-than-wide orientation, and microcalcifications as independent predictive factors for differentiating malignant from benign thyroid nodules (all p<0.05). The diagnostic performance evaluation demonstrated that UGSR achieved an AUC of 0.852 (95% CI: 0.792 - 0.912), with a sensitivity of 63.83% and specificity of 92.31%. UGSR showed significantly superior diagnostic accuracy compared to markedly hypoechogenicity (p<0.05).
Conclusion: UGSR demonstrated high specificity (92.31%) and reliability in differentiating malignant from benign thyroid nodules, suggesting its potential as a quantitative adjunct to ultrasound diagnosis, though sensitivity (63.83%) warrants combination with other features.
1 Introduction
Thyroid nodules are a widespread endocrine condition. While most thyroid nodules are benign, 5 - 10% are malignant (1). The high prevalence of this condition highlights the importance of effective screening and clinical management of thyroid nodules. Thyroid nodule detection rates vary by age and gender, with higher rates observed in women (2). Ultrasound (US) is the preferred method for assessing thyroid nodules as it offers crucial information on their size, shape, and morphological features, which aid in determining whether they are benign or malignant (3–6). Chinese Thyroid Imaging Reporting and Data System (C-TIRADS) stratifies nodules by malignancy risk using sonographic features (e.g., echogenicity, margins) (7).
The evaluation of thyroid nodules typically begins with US as the first-line imaging modality. On ultrasounds, the echogenic characteristics of thyroid nodules help differentiate benign from malignant lesions. Based on their echogenicity, a nodule may be classified as anechoic, markedly hypoechogenic, hypoechogenic, isoechogenic, or hyperechogenic. These classifications assist physicians in making initial diagnoses and determining the need for subsequent fine needle aspiration biopsies (FNA). Previous studies have shown that hypoechogenicity and markedly hypoechogenicity nodules are more frequently associated with malignancy, while isoechogenicity and hyperechogenicity nodules are typically benign (8, 9). However, these echogenic features lack specificity and often fail to meet clinical requirements. Although the link between thyroid nodule echogenicity and cancer risk has been established, it remains unclear in certain cases.
In clinical practice, evaluating thyroid nodules via ultrasonography often relies on subjective judgment, leading to inconsistent or ambiguous outcomes. Previous studies have reported significant interobserver variability in the interpretation of ultrasound images, which may contribute to unnecessary diagnostic or treatment interventions for thyroid nodules (10). To address this, the ultrasound grayscale ratio (UGSR) quantified by ImageJ software offers a potential solution by replacing visual interpretation with a quantitative measurement. This approach could provide greater objectivity and improved interobserver consistency, making it a reliable alternative for echogenicity assessment (11). Since different ultrasound machine parameters may affect raw grayscale values, we used UGSR instead of raw grayscale values for standardization.
This study aims to enhance diagnostic accuracy by quantitatively analyzing thyroid nodule echogenicity and investigating its association with benign and malignant pathology.
2 Materials and methods
2.1 Study population
We reviewed patients treated for thyroid nodules from January 2017 to December 2021. Malignancy and benignity were determined based on histopathology (biopsy or surgery) as the reference standard. Inclusion criteria: (1) thyroid ultrasound within a month before biopsy or surgery; (2) had a thyroid nodule biopsy or surgery at our hospital; (3) the ultrasound images and thyroid nodule pathology were clear. Exclusion Criteria: (1) ultrasound images inadequately displayed the thyroid nodule and its surrounding tissue; (2) nodules <5 mm or >40 mm; (3) cystic or calcified nodules with unclear borders; (4) significantly changed thyroid parenchyma echogenicity or Hashimoto’s thyroiditis; (5) previous thyroid biopsy, surgery, or neck radiation therapy. This cross-sectional study analyzed retrospective data from patients who underwent thyroid ultrasound and subsequent histopathology.
The protocol of this retrospective study has been reviewed and approved by the Ethics Review Committee of the First Affiliated Hospital of Shantou University Medical College (Ethics Review No.: B - 2024 - 239), and the committee waived the requirement for informed consent.
2.2 Instruments and methods
Our study used the Siemens Acuson S2000 ultrasound machine’s 4 – 9 MHz linear transducer. The patient lay supine with an extended head, exposing the neck to scan the nodule’s location, shape, borders, echogenicity, calcifications, and internal structures, and measure its maximum diameter. A radiologist acquired all images with 10 years of experience in thyroid imaging. Two experienced radiologists analyzed images of thyroid nodules, scoring features based on C-TIRADS criteria. Each positive feature—vertical position, solidity, very hypoechoicity, microcalcification, and margin blurring—received one point, while comet-tail artifacts resulted in a one-point deduction. The total score was calculated as follows: Category 1 for no nodule (0% malignancy rate), category 2 for a score of -1 (0% malignancy rate), category 3 for a score of 0 (<2% malignancy rate), category 4a for a score of 1 (2%-10% malignancy rate), category 4b for a score of 2 (>10%-50% malignancy rate), category 4c for scores of 3 to 4 (>50%-90% malignancy rate), and category 5 for a score of 5 (>90% malignancy rate) (7). The two radiologists were unaware of the pathology results. If two radiologists disagree on an image feature’s definition and diagnosis, the final decision will be based on their consensus.
Two experienced radiologists used ImageJ software for UGSR measurements. Two radiologists with 10 years of thyroid imaging experience performed all measurements, blinded to pathology results. ImageJ software measured the grayscale value of the entire nodule and surrounding thyroid tissue, ensuring consistent sampling depth (8, 12). The measurement area for the surrounding tissue matched the nodule’s size. Perform five measurements for each and average them. UGSR was calculated as the ratio of nodule grayscale value to surrounding thyroid parenchyma grayscale value, quantified using ImageJ software (Figure 1). ROIs were measured at the same depth in transverse/longitudinal planes. We sampled the entire nodule for nodules with heterogeneous echogenicity to avoid bias. A total of 125 nodules (78 benign, 47 malignant) were evaluated.

Figure 1. Process of extracting the grayscale ratio from thyroid nodules. Images (A, B) show the extracted gray value of the thyroid nodule, which is 7.808. Images (C, D) show the gray values extracted from the adjacent thyroid tissue (44.019). The ultrasound grayscale ratio (UGSR) is 7.808/44.019 = 17.74%.
Thirty randomly selected nodules were measured twice by the same radiologist to assess intra-observer reproducibility (ICC).
2.3 Statistical methods
All statistical analyses were performed using SPSS software (version 26.0; IBM Corp., Armonk, NY, USA) and MedCalc Statistical Software (version 18.0; MedCalc Software Ltd, Ostend, Belgium). Continuous variables were expressed as mean ± standard deviation (SD) for normally distributed data or median with interquartile range (IQR) for non-normally distributed data. Categorical variables were presented as frequencies and percentages. Comparative analyses were conducted using: Independent samples t-test or Mann-Whitney U test for continuous variables; Chi-square test or Fisher’s exact test for categorical variables; One-way ANOVA to compare UGSR distributions across different C-TIRADS categories, with post-hoc correlation analysis. For logistic regression modeling, we adopted a focused approach: Each model incorporated either UGSR or marked hypoechogenicity as the primary sonographic variable; Four additional covariates were consistently included: maximum diameter, margin characteristics, taller-than-wide orientation, and microcalcifications; Nodule malignancy status served as the dichotomous dependent variable. Diagnostic performance was evaluated through receiver operating characteristic (ROC) curve analysis, calculation of the area under the curve (AUC) with 95% confidence intervals (CI), and sensitivity and specificity determinations at optimal cut-off values. The UGSR cut-off value was established through ROC curve analysis. A two-tailed p-value <0.05 was considered statistically significant for all analyses.
3 Results
This retrospective study encompassed 118 patients, accounting for 125 nodules. Among these, 78 nodules were classified as benign, comprising 66 nodular goiters and 12 follicular adenomas. The remaining 47 nodules were categorized as malignant, including 46 papillary carcinomas and one follicular adenocarcinoma, as detailed in Table 1. In the cohort with benign nodules, the mean age was 46.79 ± 13.89 years, and the maximum nodule diameter was 20.05 mm (range: 14.10 – 26.30 mm). Conversely, in the cohort with malignant nodules, the mean age was 45.19 ± 13.73 years, with a maximum nodule diameter of 12.30 mm (range: 7.85 – 15.08 mm). Significant differences were observed between benign and malignant nodules in maximum diameter (p = 0.042), UGSR (p < 0.001), and C-TIRADS features (p < 0.001), supporting UGSR’s diagnostic utility. Conversely, the differences between the two groups concerning gender and age did not reach statistical significance (refer to Table 2). Despite a higher proportion of females in both benign (85.9%) and malignant (89.4%) groups, no significant gender difference was observed (p = 0.574), likely due to matched recruitment by nodule size/suspicion rather than population prevalence. The ICC for the measurement of thyroid nodule UGSR was 0.964 (p < 0.001), demonstrating a high level of reproducibility. Furthermore, the correlation between UGSR and the C-TIRADS was 0.540 (p < 0.001), as presented in Table 3.
In the logistic regression analysis, we constructed comprehensive models incorporating maximum diameter, margin characteristics, taller-than-wide orientation, microcalcifications, and one sonographic variable (UGSR or marked hypoechogenicity). Table 4 (Model 1 and Model 2) demonstrated that UGSR emerged as an independent predictive factor for thyroid malignancy (OR: 0.966, 95% CI: 0.936 - 0.997, p = 0.033). The diagnostic performance evaluation revealed: UGSR achieved an AUC of 0.852 (95% CI: 0.786 - 0.917) for malignancy discrimination; The optimal predictive probability cut-off was ≤45.027; Sensitivity: 63.83%; Specificity: 92.31%. Comparative analysis demonstrated UGSR’s statistically superior diagnostic accuracy over conventional marked hypoechogenicity (p<0.05 by DeLong’s test), with detailed performance metrics in Table 5 and corresponding ROC curves illustrated in Figure 2.

Table 5. Diagnostic performances of UGSR, margin, vertical orientation, microcalcification, and markedly hypoechoic.

Figure 2. Diagnostic performances of UGSR, margin, vertical orientation, microcalcification, and markedly hypoechoic.
4 Discussion
The US is the preferred method for thyroid nodule screening, and the TI-RADS system is increasingly used to standardize nodule grading based on sonographic features. Previous studies have demonstrated that positive C-TIRADS features effectively distinguish between benign and malignant thyroid nodules (7, 13–16). Nevertheless, the accuracy of radiologists in identifying these C-TIRADS features is significantly influenced by their experience and the resolution of the ultrasound equipment. Specifically, assessing nodule echogenicity is particularly vulnerable to subjective interpretation. Prior research has reported limited specificity for diagnosing malignant nodules using hypoechoic echogenicity characteristics when employing hypoechoic features (17, 18).
Additional biopsy procedures are necessary to evaluate further nodules that cannot be definitively classified as benign or malignant. Although FNA is a valuable diagnostic tool, it has several limitations (19), such as challenges related to the puncture procedure and potential inaccuracies due to sampling variability. Additionally, FNA has difficulties in differentiating between follicular adenomas and follicular carcinomas. Enhancing the diagnostic accuracy of thyroid nodules using conventional US and reducing unnecessary FNAs are of significant clinical importance. Therefore, quantifying nodule echogenicity could improve the efficacy of their diagnosis. Since the gain adjustment affects the image’s overall value, direct grayscale-based diagnosis is unreliable. Consequently, this study employed the grayscale ratio to represent nodule characteristics objectively. Since the calcified areas tend to overestimate the grayscale ratio and cystic components underestimate it, only solid nodules without significant calcification were included.
The findings demonstrated that the UGSR of benign nodules was significantly higher than malignant nodules (p < 0.001). Tissue echogenicity in ultrasound depends on the number of acoustic interfaces and the degree of acoustic impedance difference between adjacent media. Malignant thyroid nodules primarily consist of cancerous cells, microscopically resulting in a marked reduction of thyroid follicular components. Consequently, this structure does not easily establish an acoustic interface, leading to hypoechoic echogenicity (20). In contrast, thyroid adenomas are benign neoplasms composed of follicular epithelium. Histologically, they contain densely packed follicles in a gelatinous matrix, resulting in homogeneous iso- to hyperechoic patterns. Nodular goiter differs pathologically, featuring focal hyperplasia, small follicles, and areas of epithelial changes, often appearing as uneven or honeycomb-like iso- to hyperechoic lesions (21).
In nodules with heterogeneous echogenicity, manual selection of the sampling frame may introduce bias. Thus, this study used the overall gray value for analysis. The ICC for the UGSR measurement of thyroid nodules was 0.964 (p < 0.001), confirming high reproducibility and reliability. Our results suggest that UGSR measurements are consistent across different technical parameters, including gain, dynamic range, and transducer frequency. Furthermore, the grayscale ratio normalizes operator-dependent variability, providing an objective quantitative method.
Logistic regression analyses incorporated five sonographic parameters: maximum diameter, margin irregularity, taller-than-wide orientation, microcalcifications, and UGSR or marked hypoechogenicity. As shown in Table 4 (Models 1 and 2), UGSR was an independent predictor of thyroid malignancy (OR: 0.966; 95% CI: 0.936 - 0.997; p = 0.033). Diagnostic performance analysis showed: Superior discriminatory AUC for UGSR (0.852; 95% CI: 0.786 - 0.917); Optimal threshold: ≤45.027; Sensitivity: 63.83%; Specificity: 92.31%; Significantly better than marked hypoechogenicity (p<0.05 by DeLong’s test). Our UGSR results (AUC = 0.852) align with Yun et al.’s findings (AUC = 0.856–0.875), supporting its generalizability. However, their cohort achieved higher sensitivity (78.9%–86.9% vs. 63.8%), likely because they focused on 3 – 10 mm nodules (8).
Smaller malignant nodule size (12.3 vs. 20.1 mm, p = 0.042) reflects clinical practice: (1) Suspicious features (e.g., microcalcifications) prompt early intervention; (2) Large benign nodules are often resected for compressive symptoms, introducing selection bias.
The limitations of this study are as follows: (1) Single-center retrospective design, which may introduce selection bias and limit generalizability. Multicenter prospective studies are needed to validate UGSR. (2) Small sample size; future studies should expand cohorts for robustness. While our sample (n=125) suffices for pilot validation (power=0.8, effect size=0.5), multicenter studies with >500 nodules are needed for definitive cut-off establishment. (3) Potential variability across radiologists and ultrasound instruments requires further investigation. (4) Inclusion bias: Only biopsied/surgical nodules were analyzed, which tend to have higher malignancy rates and larger sizes than screening-detected nodules. (5) Reference standard bias: Biopsy/surgical histopathology may miss some benign lesions. (6) Technical variability: UGSR performance may differ with ultrasound systems; protocol standardization is needed. (7) Lack of subgroup analysis by nodule size due to sample size constraints.
5 Conclusion
UGSR demonstrated high specificity (92.31%) and reliability in differentiating malignant from benign thyroid nodules, suggesting its potential as a quantitative adjunct to ultrasound diagnosis, though sensitivity (63.83%) warrants combination with other features.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Ethical Review Committee at the First Affiliated Hospital of the Medical College of Shantou University Ethics review No. B - 2024 - 239). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin because Based on the characteristics of the retrospective study and the principle of data desensitization processing, all subjects' personal information was processed using encryption technology to exempt the informed consent process.
Author contributions
KL: Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. TL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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. Jiang L, Liu D, Long L, Chen J, Lan X, and Zhang J. Dual-source dual-energy computed tomography-derived quantitative parameters combined with machine learning for the differential diagnosis of benign and Malignant thyroid nodules. Quant Imaging Med Surg. (2022) 12:967–78. doi: 10.21037/qims-21-501
2. Li X, Du H, Luo J, Ding W, Lai B, He J, et al. Comparison of the clinical validity of droplet digital PCR to ARMS-PCR for BRAF V600E mutation detection in thyroid nodules. J Clin Lab Anal. (2020) 34:e23458. doi: 10.1002/jcla.23458
3. Min PH, Jihye L, Young KJ, Youngjean PV, Miribi R, Minah L, et al. Using ultrasonographic features to predict the outcomes of patients with small papillary thyroid carcinomas: a retrospective study implementing the 2015 ATA patterns and ACR TI-RADS categories. Ultrasonography. (2022) 41:298–306. doi: 10.14366/usg.21097
4. Zhang WB, Xu W, Fu WJ, He BL, Liu H, and Deng WF. Comparison of ACR TI-RADS, Kwak TI-RADS, ATA guidelines and KTA/KSThR guidelines in combination with SWE in the diagnosis of thyroid nodules. Clin Hemorheol Microcirc. (2021) 78:163–74. doi: 10.3233/CH-201021
5. Wang Y, Dong T, Nie F, Wang G, Liu T, and Niu Q. Contrast-enhanced ultrasound in the differential diagnosis and risk stratification of ACR TI-RADS category 4 and 5 thyroid nodules with non-hypovascular. Front Oncol. (2021) 11:662273. doi: 10.3389/fonc.2021.662273
6. Trimboli P, Castellana M, Piccardo A, Romanelli F, Grani G, Giovanella L, et al. The ultrasound risk stratification systems for thyroid nodule have been evaluated against papillary carcinoma. A meta-analysis. Rev Endocr Metab Disord. (2021) 22:453–60. doi: 10.1007/s11154-020-09592-3
7. Zhou J, Yin L, Wei X, Zhang S, Song Y, Luo B, et al. 2020 Chinese guidelines for ultrasound Malignancy risk stratification of thyroid nodules: the C-TIRADS. Endocrine. (2020) 70:256–79. doi: 10.1007/s12020-020-02441-y
8. Gong Y, Yao X, Yu L, Wei P, Han Z, Fang J, et al. Ultrasound grayscale ratio: a reliable parameter for differentiating between papillary thyroid microcarcinoma and micronodular goiter. BMC Endocr Disord. (2022) 22:75. doi: 10.1186/s12902-022-00994-9
9. Low G, Bara M, Du Y, Katlariwala P, Croutze R, Resch K, et al. Tips for improving consistency of thyroid nodule interpretation with ACR TI-RADS. J Ultrason. (2022) 22:e51–6. doi: 10.15557/JoU.2022.0009
10. Ma W, Li X, Zou L, Fan C, and Wu M. Symmetrical awareness network for cross-site ultrasound thyroid nodule segmentation. Front Public Health. (2023) 11:1055815. doi: 10.3389/fpubh.2023.1055815
11. Wu MH, Chen KY, Chen A, and Chen CN. Software-based analysis of the taller-than-wide feature of high-risk thyroid nodules. Ann Surg Oncol. (2021) 28:4347–57. doi: 10.1245/s10434-020-09463-w
12. Hu T, Zhou T, Zhang Y, Zhou L, Huang X, Cai Y, et al. The predictive value of the thyroid nodule benign and Malignant based on the ultrasound nodule-to-muscle grayscale ratio. J Clin Ultrasound. (2024) 52:51–8. doi: 10.1002/jcu.23601
13. Hu Y, Xu S, and Zhan W. Diagnostic performance of C-TIRADS in Malignancy risk stratification of thyroid nodules: A systematic review and meta-analysis. Front Endocrinol (Lausanne). (2022) 13:938961. doi: 10.3389/fendo.2022.938961
14. Renxu Li, Zhenwei L, Xiangyu W, and Luzeng C. Role of echogenic foci in ultrasonographic risk stratification of thyroid nodules: Echogenic focus scoring in the American College of Radiology Thyroid Imaging Reporting and Data System. Front Oncol. (2022) 12:929500. doi: 10.3389/fonc.2022.929500
15. Liu J, Luo T, Zhang H, Liu H, Gu Y, Chen X, et al. Markedly hypoechoic: a new definition improves the diagnostic performance of thyroid ultrasound. Eur Radiol. (2023) 33:7857–65. doi: 10.1007/s00330-023-09828-1
16. Qu C, Li HJ, Gao Q, Zhang JC, and Li WM. Alteration trend and overlap analysis of positive features in different-sized benign and Malignant thyroid nodules: based on chinese thyroid imaging reporting and data system. Int J Gen Med. (2024) 17:1887–95. doi: 10.2147/IJGM.S461076
17. Prieditis P, Radzina M, Mikijanska M, Liepa M, Stepanovs K, Grani G, et al. Non-marked hypoechogenic nodules: multicenter study on the thyroid Malignancy risk stratification and accuracy based on TIRADS systems comparison. Medicina (Kaunas). (2022) 58(2):257. doi: 10.3390/medicina58020257
18. Petersen M, Schenke SA, Firla J, Croner RS, and Kreissl MC. Shear wave elastography and thyroid imaging reporting and data system (TIRADS) for the risk stratification of thyroid nodules-results of a prospective study. Diagn (Basel). (2022) 12(1):109. doi: 10.3390/diagnostics12010109
19. Fabbri C, Fornelli A, Fuccio L, Giovanelli S, Tarantino I, Antonini F, et al. High diagnostic adequacy and accuracy of the new 20G procore needle for EUS-guided tissue acquisition: Results of a large multicentre retrospective study. Endosc Ultrasound. (2019) 8:261–8. doi: 10.4103/eus.eus_14_19
20. Su Z, Bao W, Yang G, Liu J, and Zhao B. SOX12 promotes thyroid cancer cell proliferation and invasion by regulating the expression of POU2F1 and POU3F1. Yonsei Med J. (2022) 63:591–600. doi: 10.3349/ymj.2022.63.6.591
Keywords: thyroid nodule, Chinese thyroid imaging reporting and data system, ultrasound grayscale ratio, malignant thyroid nodule, imageJ
Citation: Lian K and Lin T (2025) The predictive value of the ultrasound grayscale ratio for identifying malignant thyroid nodules. Front. Endocrinol. 16:1665113. doi: 10.3389/fendo.2025.1665113
Received: 13 July 2025; Accepted: 19 August 2025;
Published: 29 August 2025.
Edited by:
Erivelto Martinho Volpi, Hospital Alemão Oswaldo Cruz, BrazilReviewed by:
Maha Bondok, Alexandria University, EgyptXiyue Yang, The Eighth Affiliated Hospital of Guangxi Medical University, China
Copyright © 2025 Lian and Lin. 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: Teng Lin, dGVuZ2xpbmdkc3R1bUAxNjMuY29t