- Department of Ultrasonography, Shanxi Bethune Hospital, Shanxi Academy of Medical Sciences, Tongji Shanxi Hospital, The Third Hospital of Shanxi Medical University, Taiyuan, China
Objective: To analyze the ultrasonic characteristics of false-negative and false-positive results of shear wave elastography (SWE) in the diagnosis of thyroid nodules to clarify the influence of nodular characteristics on SWE and to guide the clinical application of SWE.
Methods: A total of 435 thyroid nodules from 343 patients with the diagnosis confirmed by surgical pathology were analyzed. Preoperative ultrasonography and SWE were conducted. The conventional ultrasound characteristics of thyroid nodules and the maximum Young’s modulus were recorded. The false negativity and false positivity of SWE for the diagnosis of thyroid nodules were calculated. The ultrasonic characteristics of thyroid nodules with SWE false results were analyzed, and logistic regression analysis was adopted to determine the ultrasonic characteristics associated with SWE false results of thyroid nodules.
Results: Among 323 malignant nodules, the SWE false negativity was 27.2% (88/323). The false positivity of SWE in 112 benign nodules was 19.6% (22/112). Regression analysis showed that an increase in the nodule volume increased the risk of SWE false-positive results (odds ratio [OR] 3.286; 95% confidence interval [CI]: 1.572–6.871; P = 0.002) and decreased the risk of false-negative results (OR 0.238; 95% CI: 0.115–0.493; P < 0.001). Nodules with coarse calcification had an increased risk of SWE false-positive results compared with those without calcification (OR 5.303; 95% CI: 1.098–25.619; P = 0.038). However, nodules with scattered hyperechoic foci had a reduced risk of SWE false-negative results (OR 0.515; 95% CI: 0.280–0.951; P = 0.034).
Conclusion: Nodular size and calcification were correlated with SWE false results, and the clinical application of SWE should be combined with conventional ultrasound features. Fine needle aspiration or a puncture biopsy should be conducted if necessary.
Introduction
Thyroid nodules are a common clinical problem that 19%–68% of the population suffer from, and thyroid cancers account for 7%–15% among these (1). Ultrasound is the first choice to determine the nature of thyroid nodules. Morphological characteristics are an important basis for differentiating between benign and malignant thyroid nodules. In addition, changes in tissue texture are closely correlated with the pathological development of thyroid nodules (2).
Shear wave elastography (SWE) is a new method to qualitatively and quantitatively evaluate tissue hardness (3) and is suggested to be used as a supplement to conventional ultrasound to determine the nature of thyroid nodules, rather than replacing conventional ultrasound as a separate diagnostic tool. Accumulating studies have shown that SWE has a certain value in distinguishing benign and malignant thyroid nodules, but no widely applicable diagnostic criteria have been established, which is mainly due to the fact that the results of SWE will be affected by many factors (4, 5).
The histological characteristics of nodules may affect the results of SWE (6), and the sonographic features can, to some extent, reflect the characteristics of nodules. A comparative analysis of the nodular features of misdiagnosed and correctly diagnosed nodules by SWE and the identification of which features are likely to produce incorrect results will help to determine the accuracy of SWE results. Therefore, in order to guide the clinical application of SWE, this study investigated the impact of nodular features on SWE through the analysis of SWE in the diagnosis of thyroid nodules in the false negative and false positive results of the relevant sonographic features.
Materials and Methods
General Data
From June 2016 to July 2018, 343 patients with 435 thyroid nodules were enrolled in the study. All the patients underwent thyroid surgery and had pathological results. There were 78 males and 265 females aged from 24 to 77 years, with an average age of 47.9 ± 10.5 years. The maximum diameter of the nodules ranged from 3.0 mm to 57.0 mm, with an average of 11.0 ± 7.7 mm.
The inclusion criteria were as follows: (1) patients aged >18 years old without treatment or puncture before the ultrasound examination; (2) patients who underwent two-dimensional ultrasound and SWE ultrasound examination before surgery.
The exclusion criteria were as follows: (1) patients with cystic changes in the nodules >25%; (2) patients with loss of elastic image information because of coarse calcification in or around the nodules; (3) patients with nodules in the isthmus or adjacent to the trachea cartilage or the arteria carotis communis, since the lateral displacement of the pulsing may cause artifacts in the hardness of the nodules; (4) patients with benign and malignant nodules in the same thyroid lobe.
Research Method
A SuperSonic Imagine Aixplorer® diasonograph made in France was used for thyroid conventional ultrasound and SWE with the linear probe set at a frequency of 4–15 MHz or 2–10 MHz. All images were collected by the same sonographer, who had more than 10 years of experience.
Conventional ultrasound examination: The patients were maintained in a supine position with full exposure of the neck. The probe was placed lightly on the neck without pressure. The patients could not move or swallow during the examination. The nodule features observed were (1) size; (2) composition: mixture, solid; (3) echo: hyperechoic or isoechoic, hypoechoic, extremely hypoechoic; (4) shape: an anteroposterior/transverse diameter (AP/T) ratio >1, AP/T ratio <1; (5) edge: clear, lobed, or irregular extraglandular invasion; (6) calcification: microcalcification, extensive calcification, no calcification. All nodules were scored according to the American College of Radiology Thyroid Imaging Reporting and Data System (ACR–TIRADS) (7).
Shear wave elastography: After a conventional ultrasound examination, the target nodule was identified on the longitudinal section of the thyroid gland, and the instrument was switched to the SWE mode (with a scale plate of 0–100 kPa). A sampling frame with the inclusion of the entire lesion and the surrounding normal thyroid tissue was placed on the nodule, and the patient was asked to hold their breath and stand for 3–5 s. Once stabilized, the image was frozen and an elastic image was obtained. On the image, blue represented the relatively soft tissue and red represented the relatively hard tissue. The circular region of interest (Q–Box™) was placed at the hardest part of the solid nodules with a diameter of 2–4 mm to avoid visible calcification. The maximum elastic value (Emax), minimum elastic value (Emin), and the average elastic value (Emean) of each nodule were automatically generated. Each nodule was measured three times, and the average value was calculated.
Analysis Indicators
The characteristics of SWE in the diagnosis of thyroid nodules were compared between the group with false positivity and the group with false negativity, and the related factors leading to false-negative and false-positive results were analyzed.
Statistical Analysis
The SPSS® Statistics 22.0 statistical software was used for statistical analysis, and the pathological results were taken as a reference standard to construct the receiver operating characteristic (ROC) curves of Emax, Emean, and Emin. The area under the curve (AUC) of each parameter was calculated. The parameter with the highest AUC value was taken as the SWE diagnostic parameter, and the optimal cut-off point was taken as the diagnostic criterion. In pathologically malignant nodules, those with an SWE diagnostic parameter greater than the diagnostic criterion were classified as the true-positive group, and those with an SWE diagnostic parameter less than the diagnostic criterion were classified as the false-negative group. In pathologically benign nodules, those with an SWE diagnostic parameter less than the diagnostic criterion were classified as the true-negative group, and those with an SWE diagnostic parameter greater than the diagnostic criterion were classified as the false-positive group.
The measurement data that satisfied the normal distribution were expressed as means ± standard deviations, and the differences between groups were compared by the independent sample t-test. The measurement data that did not satisfy the normal distribution were represented by medians and inter-quartile ranges and analyzed by the nonparametric rank-sum test. The enumeration data were analyzed by the chi-squared test or Fisher’s exact test. Multivariate analysis of the SWE false results was conducted by logistic regression analysis; P < 0.05 was considered statistically significant.
Results
Pathologic Results
A total of 343 patients had 435 nodules, including 112 (25.7%) benign nodules. There were 103 cases with nodular goiter, 6 cases with adenoma, 1 case with subacute thyroiditis, and 2 cases with Hashimoto’s thyroiditis. There were 323 (74.3%) malignant nodules, including 318 papillary carcinomas, 3 follicular adenocarcinomas, and 2 medullary carcinomas.
Diagnostic Performance of the SWE Parameters
The SWE parameters Emax, Emean, and Emin of the malignant nodules were significantly higher than those of the benign ones (P < 0.001 in all) (Table 1). Among the quantitative parameters of SWE, when Emax was at the optimal cut-off point of 36.15 kPa, the AUC reached 0.818 (Table 2), with a sensitivity of 72.8%, a specificity of 80.4%, a positive predictive value of 91.4%, a negative predictive value of 50.6%, and 74.7% accuracy.
Table 1 Comparison of the SWE quantitative parameters between the benign and malignant thyroid nodules.
Univariate Analysis of the Factors Influencing the Diagnostic Results of SWE
When “Emax > 36.15 kPa” was used as the standard for a positive result for SWE, 19.6% (22/112) of the benign nodules had an Emax greater than 36.15 kPa and were classified in the false-positive group (Figure 1), while 80.4% (90/112) of the benign nodules had an Emax less than 36.15 kPa and were classified in the true-negative group (Table 3). The sonographic findings of the nodule size, edge, and calcification were correlated with the SWE false-positive results (P < 0.05). In the ACR–TIRADS classification, the proportion of nodules with an ACR–TIRADS grade of 1–3 in the benign nodule false-positive group (22.7%, 5/22) was significantly lower than that in the true-negative group (60.0%, 54/90) (P = 0.021).
Figure 1 Nodular goiter (shear wave elastography false-positive nodule) by two-dimensional ultrasound and elastic ultrasound. (A) A solid isoechoic nodule in the middle of the right lobe of the thyroid, 18.1 × 9.6 mm, with an ill-defined margin and coarse calcification. The ACR–TIRADS score was 4 and ACR–TIRADS category was 4. (B) The Emax value of the nodule on SWE imaging was 72.6 kPa.
When “Emax > 36.15 kPa” was used as the standard for a positive result for SWE, 27.2% (88/323) of the malignant nodules had an Emax less than 36.15 kPa and were classified in the false-negative group (Figure 2), while 72.8% (235/323) of the malignant nodules had an Emax greater than 36.15 kPa and were classified in the true-positive group (Table 3). Nodule size, edge, AP/T ratio, and calcification were correlated with the SWE false-negative results (P < 0.05). In the ACR–TIRADS classification, there was no significant difference in the proportion of nodules with an ACR–TIRADS grade of 4–5 between the false-negative group (96.6%, 85/88) and the true-positive group (99.2%, 233/235) (P = 0.053).
Figure 2 Papillary thyroid carcinoma (SWE false-negative nodule) by two-dimensional ultrasound and elastic ultrasound. (A) A solid, very hypoechoic nodule in the middle-left lobe of the thyroid, 6.0 × 3.0 mm, with an ill-defined margin and taller than wide in shape. The ACR-TIRADS score was 8 and ACR–TIRADS category was 5. (B) The Emax value of the nodule on SWE imaging was 20.4 kPa.
Multivariate Analysis of the Factors Influencing the Diagnostic Results of SWE
Multiple logistic regression analysis of the factors correlated with the false-positive results in SWE: Multiple logistic regression equations were established by incorporating the nodule size, edge, and calcification. The results suggested that the risk of SWE false-positive results increased with an increase in the nodule volume (OR 3.286; 95% CI: 1.572–6.871; P = 0.002). Nodules with coarse calcification had an increased risk of false-positive SWE results compared with those without nodular calcification (OR 5.303; 95% CI: 1.098–25.619; P = 0.038) (Table 4).
Multiple logistic regression analysis of the factors correlated with the false-negative results in SWE: Multiple logistic regression equations were established by incorporating the nodule size, edge, AP/T ratio, and calcification. The results suggested that the risk of SWE false-negative results decreased with an increase in the nodule volume (OR 0.238; 95% CI: 0.115–0.493; P < 0.001). Conversely, the smaller the nodule volume, the greater the risk of SWE false-negative results. Compared with the nodules without calcification, nodules with scattered hyperechoic foci had a reduced risk of SWE false-negative results (OR 0.51; 95% CI: 0.280–0.951; P = 0.034) (Table 5).
Discussion
Many studies have shown that SWE is valuable in diagnosing the nature of thyroid nodules by evaluating the inherent characteristics of thyroid nodules. However, there are still some problems in the diagnosis of SWE, such as the large overlap of the elastic values of some benign and malignant nodules, and the wide range of cutoff values for the diagnosis of thyroid malignant nodules, ranging from 26 to 85kPa (8–10), which may be related to the great heterogeneity between different studies or between nodules. Swan et al. (6) reported that the diagnostic efficacy of SWE at the individual level is not satisfactory, so it seems difficult to establish a clear and widely applicable diagnostic criteria in the short term. However, the clinical application of SWE is for individual patients, so it is particularly important to interpret the reliability of SWE results. Therefore, by comparing the sonographic characteristics of false negative and false positive results of SWE in the diagnosis of thyroid nodules, and analyzing the influence of ultrasound characteristics of thyroid nodules on SWE, this study helps to judge whether the results of SWE are reliable. In this study, the ROC curve was used to determine the best diagnostic index Emax. When 36.15kPa, the cutoff value of Emax, was used as the diagnostic criteria, the false positive rate of benign nodules was 19.6% (22cm/112) and the false negative rate of malignant nodules was 27.2% (88/323).
Shang et al. (11) showed that the internal composition of nodules changes with the increase of nodular volume, and the hardness increases. The results of the present study show that the size of nodules impacted the SWE results. The greater the nodule volume, the higher the risk of SWE false-positive results, and the lower the risk of false-negative results. Therefore, too large or too small nodules will affect the diagnostic efficacy of SWE, which is consistent with the research results of Wang et al. (12) Elastography has a low sensitivity for nodules smaller than 10 mm and low specificity for nodules larger than 20 mm. The components of benign and malignant thyroid nodules change during their development. When malignant nodules are small, the fibrous components are relatively small with low hardness, but as they increase in volume, fiber composition increases and calcification may occur. As benign nodules increase in volume, extensive fibrosis and hyaline changes occur, both of which will cause an increase in the nodular hardness (13, 14).
In addition, calcification has been thought to be an important factor in increasing the hardness value. Among malignant nodules in the present study, the proportion of nodules with scattered hyperechoic foci in the true-positive group (51.1%, 120/235) was significantly higher than that in the false-negative group (25.0%, 22/88). Multivariate regression analysis showed that compared with the non-calcified nodules, nodules with scattered hyperechoic foci had a lower risk of SWE false-negative results. In the benign nodules, the percentages of coarse calcification (18.2%, 4/22) and scattered hyperechoic foci (27.3%, 7/22) in the false-positive group were significantly higher than those in the true-negative group (5.6%, 5/90; and 7.8%, 7/90, respectively). Multivariate regression analysis showed that coarse calcification was a risk factor for the SWE false-positive results. This is consistent with the results of Gregory et al. (15), which meant that cluster microcalcification and coarse calcification might lead to a significant increase in the hardness value. The scattered hyperechoic foci in the benign nodules was not a risk factor for the SWE false-positive results, possibly because part of the scattered hyperechoic foci in the benign nodules was not real calcification but concentrated colloid or fibrous tissue (16). The effect on the SWE value was relatively small; therefore, attention should be paid to identification in clinical application.
The results of multivariate analysis in the present study showed that the AP/T ratio and edge had no significant correlation with the SWE false results. The AP/T ratio is the ratio of the anterior and posterior diameter to the transverse diameter of nodules, reflecting the growth mode of nodules, while SWE mainly reflects the internal characteristics of nodules (13). An AP/T ratio >1 is considered an important sign of thyroid malignant nodules in conventional ultrasound, especially for nodules smaller than 10 mm, with a high sensitivity (81.4%) and specificity (96.8%) for the diagnosis of malignant nodules (17). Therefore, for nodules smaller than 10 mm, the combination of two-dimensional ultrasonic characteristics may be more helpful for judging the nature of nodules. Yoo et al. (18) found that the combination of the Emax value of nodules and the aspect ratio information of nodules is not only valuable in distinguishing benign from malignant nodules, but also helpful to predict the pathological types of thyroid nodules. However, two-dimensional ultrasound is of limited value when large, benign nodules with coarse calcification present a high risk of a false-positive result in SWE. In the present study, the proportion of nodules with an ACR–TIRADS grade of 4–5 in the benign nodules’ SWE false-positive group (77.3%, 17/22) was significantly higher than that in the true-negative group (40.0%, 36/90). Therefore, for nodules with an ACR–TIRADS score of Class 4 or above, fine needle aspiration or biopsy is necessary for further diagnosis.
Limitations
This study only analyzed the correlation between the sonographic features of nodules and the results of SWE diagnosis, without considering other related factors, such as nodular location, depth and surrounding tissue. Moreover, this study is a retrospective study with a high proportion of malignant nodules and a small number of benign nodules, and there is a certain selection bias. Finally, the malignant nodules in this study are mainly papillary carcinoma, and other pathological properties such as follicular carcinoma and medullary carcinoma need to be further studied.
Conclusion
when the volume of malignant nodules is less than 10mm, SWE is predisposed to produce false negative results, and combining two-dimensional ultrasound features is more helpful to improve the diagnostic accuracy. For nodules with clustered microcalcification, the risk of false negative is reduced. For benign nodules with large volume and coarse calcification, SWE has a higher risk of producing false positive results, and further diagnosis by FNA or fine needle biopsy is needed if necessary. Therefore, this study has a certain guiding significance in defining the scope and value of clinical application of SWE.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by The Third Hospital of Shanxi Medical University, Shanxi Bethune Hospital, Shanxi Academy of Medical Sciences. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conception and design of the research: C-sK, J-pX; Acquisition of data: J-wM and Y-xZ; Analysis and interpretation of the data: X-yK and H-zL; Statistical analysis: X-yK, J-wM, and F-cY; Obtaining financing: J-pX. Writing of the manuscript: J-pX. Critical revision of the manuscript for intellectual content: C-sK. All authors read and approved the final draft.
Funding
Key research and Development Projects in Shanxi Province (grant number 201803D31143).
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. Haugen BR, Alexander EK, Bible KC, Doherty GM, Mandel SJ, Nikiforov YE, et al. American Thyroid Association Management Guidelines for Adult Patients With Thyroid Nodules and Differentiated Thyroid Cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid (2016) 26(1):1–133. doi: 10.1089/thy.2015.0020
2. Bhatia KS, Tong CS, Cho CC, Yuen EH, Lee YY, Ahuja AT. Shear Wave Elastography of Thyroid Nodules in Routine Clinical Practice:Preliminary Observations and Utility for Detecting Malignancy. Eur Radiol (2012) 22(11):2397–406. doi: 10.1007/s00330-012-2495-1
3. Fu C, Cui KF, Qing SC. Clinical Application of Real-Time Shear Wave Elastography in Diagnosis of Benign and Malignant Thyroid Solid Nodules. Chin J Ultrason (2012) 21(1):49–51.
4. Szczepanek-Parulska E, Woliński K, Stangierski A, Gurgul E, Ruchała M. Biochemical and Ultrasonographic Parameters Influencing Thyroid Nodules Elasticity. Endocrine (2014) 47(2):519–27. doi: 10.1007/s12020-014-0197-y
5. Kim HG, Moon HJ, Kwak JY, Kim EK. Diagnostic Accuracy of the Uhrasonographic Features for Subcentimeter Thyroid Nodules Suggested by the Revised American Thyroid Association Guidelines. Thyroid (2013) 23(12):1583–9. doi: 10.1089/thy.2012.0586
6. Swan KZ, Nielsen VE, Bonnema SJ. Evaluation of Thyroid Nodules by Shear Wave Elastography: A Review of Current Knowledge. J Endocrinol Invest (2021) 44(10):2043–56. doi: 10.1007/s40618-021-01570-z
7. Tessler FN, Middleton WD, Grant EG, Hoang JK, Berland LL, Teefey SA, et al. ACR Thyroid Imaging, Reporting and Data System (TI-RADS): White Paper of the ACR TI-RADS Committee. J Am Coll Radiol (2017) 14(5):587–95. doi: 10.1016/j.jacr.2017.01.046
8. Cosgrove D, Barr R, Bojunga J, Cantisani V, Chammas MC, Dighe M, et al. WFUMB Guidelines and Recommendations on the Clinical Use of Ultrasound Elastography: Part 4. Thyroid. Ultrasound Med Biol (2017) 43(1):4–26. doi: 10.1016/j.ultrasmedbio.2016.06.022
9. Dobruch-Sobczak K, Zalewska EB, Guminska A, Słapa RZ, Mlosek K, Wareluk P, et al. Diagnostic Performance of Shear Wave Elastography Parameters Alone and in Combination With Conventional B-Mode Ultrasound Parameters for the Characterization of Thyroid Nodules: A Prospective, Dual-Center Study. Ultrasound Med Biol (2016) 42(12):2803–11. doi: 10.1016/j.ultrasmedbio.2016.07.010
10. Swan KZ, Nielsen VE, Bibby BM, Bonnema SJ. Is the Reproducibility of Shear Wave Elastography of Thyroid Nodules High Enough for Clinical Use? A Methodological Study. Clin Endocrinol (2017) 86(4):606–13. doi: 10.1111/cen.13295
11. Shang H, Wu B, Liu Z, Liu Y, Cheng W. The Effectiveness of Shear Wave Elastography in the Diagnosis of PTMC. Technol Health Care (2020) 28(2):221–6. doi: 10.3233/THC-191895
12. Wang F, Chang C, Chen M, Gao Y, Chen YL, Zhou SC, et al. Does Lesion Size Affect the Value of Shear Wave Elastography for Differentiating Between Benign and Malignant Thyroid Nodules? J Ultrasound Med (2018) 37(3):601–9. doi: 10.1002/jum.14367
13. Fukuhara T, Matsuda E, Endo Y, Takenobu M, Izawa S, Fujiwara K, et al. Correlation Between Quantitative Shear Wave Elastography and Pathologic Structures of Thyroid Lesions. Ultrasound Med Biol (2015) 41(9):2326–32. doi: 10.1016/j.ultrasmedbio.2015.05.001
14. Ha EJ, Baek JH, Lee JH, Song DE, Kim JK, Shong YK, et al. Sonographically Suspicious Thyroid Nodules With Initially Benign Cytologic Results: The Role of a Core Needle Biopsy. Thyroid (2013) 23(6):703–8. doi: 10.1089/thy.2012.0426
15. Gregory A, Bayat M, Denis M, Mehrmohammadi M, Fatemi M, Alizad A, et al. An Experimental Phantom Study on the Effect of Calcifications on Ultrasound Shear Wave Elastography. Annu Int Conf IEEE Eng Med Biol Soc (2015) 2015:3843–6. doi: 10.1109/EMBC.2015.7319232
16. Malhi H, Beland MD, Cen SY, Allgood E, Daley K, Martin SE, et al. Echogenic Foci in Thyroid Nodules: Significance of Posterior Acoustic Artifacts. AJR Am J Roentgenol (2014) 203(6):1310–6. doi: 10.2214/AJR.13.11934
17. Ren J, Liu B, Zhang LL, Li HY, Zhang F, Li S, et al. A Taller-Than-Wide Shape Is a Good Predictor of Papillary Thyroid Carcinoma in Small Solid Nodules. J Ultrasound Med (2015) 34(1):19–26. doi: 10.7863/ultra.34.1.19
Keywords: thyroid, nodule characteristics, shear wave elastography, influence, ultrasonography
Citation: Xue J-p, Kang X-y, Miao J-w, Zhang Y-x, Li H-z, Yao F-c and Kang C-s (2022) Analysis of the Influence of Thyroid Nodule Characteristics on the Results of Shear Wave Elastography. Front. Endocrinol. 13:858565. doi: 10.3389/fendo.2022.858565
Received: 20 January 2022; Accepted: 09 May 2022;
Published: 10 June 2022.
Edited by:
Bernadette Biondi, University of Naples Federico II, ItalyReviewed by:
Serena Ippolito, ASL Napoli 1 centro, ItalyErivelto Martinho Volpi, Centro de referencia no ensino do diagnóstico por imagem (CETRUS), Brazil
Sara Garberoglio, Humanitas Cellini Clinic, Italy
Copyright © 2022 Xue, Kang, Miao, Zhang, Li, Yao and Kang. 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: Xiao-yan Kang, kcsdr_09sea@163.com