- 1Shanxi Eye Hospital Affiliated to Shanxi Medical University, Taiyuan, Shanxi, China
- 2Ruiz Department of Ophthalmology and Visual Science, McGovern Medical School at The University of Texas Health Science Center at Houston, Houston, TX, United States
Purpose: To compare the astigmatism prediction accuracy of total keratometry (TK) from the IOLMaster 700 and total corneal power (TCP) from Anterion based on swept-source optical coherence tomography (SS-OCT) technology in toric intraocular lens (toric IOL) calculation.
Design: A retrospective observational study.
Methods: Total corneal astigmatism (TCA) were obtained using IOLMaster 700 and Anterion. Z CALC 2.0 was used to calculate the expected postoperative refractive astigmatism in conjunction with TCA. Prediction errors (PE) in refractive outcomes was analyzed 1 month postoperatively using the vector analysis by the Holladay method, including the mean vector PE magnitude, percentage of cases with vector PE in certain intervals, and the centroid PE.
Results: A total of 56 eyes from 56 patients were enrolled in the study with an insertion of an AT TORBI 709 toric IOL. The difference in mean vector PE of postoperative refractive astigmatism between TK and TCP was not statistically significant (0.48D versus 0.46D, P = 0.281). TK and TCP yielded 27.3 and 40.0% of eyes with vector PE ≤ 0.25D, and 58.2 and 63.6% with vector PE ≤ 0.5D (both P > 0.05), respectively. TK and TCP resulted in similar ATR centroid PE of 0.10D@35° ± 0.60D and 0.15D@22° ± 0.57D, respectively, and there were no significant differences between x-PE component and y-PE component.
Conclusion: IOLMaster 700 and Anterion provided comparable astigmatic predictability in toric IOL implantation using total keratometry and Z CALC 2.0.
1 Introduction
Corneal astigmatism is one of the critical factors leading to poor visual acuity after cataract surgery. Although toric intraocular lens (toric IOL) implantation can dramatically improve the refractive outcome of patients with pre-existing corneal astigmatism (1, 2), unexpected residual astigmatism remains a concern, which can be caused by a variety of factors, such as inaccurate preoperative corneal astigmatism measurement (2), surgical-induced corneal astigmatism (SIA), misalignment of toric IOL (3), and inappropriate way to calculate toric IOL diopter. Among the above factors, posterior corneal astigmatism (PCA) has the greatest influence on astigmatic prediction error (4). Much literature has discussed the contribution of PCA to total corneal astigmatism (TCA) (5, 6). Studies have shown that ignoring PCA may cause estimation errors in TCA, leading to overcorrection of with the rule (WTR) astigmatism and undercorrection of the against the rule (ATR) astigmatism (4, 5, 7). Therefore, when calculating the toric IOL power, considering PCA can improve the accuracy of refractive results (8, 9).
At present, with the development of examination technology, there are several biometers that are now available to acquire ocular anterior segment parameters and measure astigmatism on the front and back surface of the cornea, such as Pentacam (Optikgerate GmbH, Wetzlar, Germany), IOLMaster 700 (Carl Zeiss Meditec, Jena, Germany), and Anterion (Heidelberg Engineering, Heidelberg, Germany), of which the latter two were based on swept-source optical coherence tomography (SS-OCT) technology. Using the Alcon online toric calculator, Choi et al. (10) compared the accuracy of TCA from the IOLMaster 700 and Pentacam when implanting of Acrysof IQ Panoptix toric TFNT IOLs (Alcon Laboratories, Inc., Fort Worth, TX). It was found that favorable refractive outcomes were provided by TCA from the two biometers. Wang and Koch (11) compared the accuracy of astigmatic prediction by Barrett toric calculator using the PCA measured by the IOLMaster 700 and the predicted PCA, respectively, and the results showed that astigmatic prediction error (PE) was significantly improved with the measured PCA. However, to our knowledge, no literature has explored the performance of TCP from Anterion in astigmatic prediction. Therefore, this study aimed to compare the PEs in postoperative refractive astigmatism based on TK with the IOLMaster 700 and TCP with the Anterion.
2 Patients and methods
2.1 Study population
This retrospective study was reviewed and approved by the institutional review board of Shanxi Eye Hospital, Affiliated to Shanxi Medical University (No. SXYYLL-20210107), which conformed to the Declaration of Helsinki principles. A list of patients who underwent toric IOL implantation at Shanxi Eye Hospital was retrospectively obtained from Jan 2021 to July 2022. The inclusion criteria were senile cataract patients with reliable preoperative measurements and regular corneal astigmatism between 0.75D and 5D. The exclusion criteria were cases with systemic diseases and ophthalmologic diseases that affect postoperative refractive and visual acuity, intraoperative or postoperative complications, keratoconus and contact lens wearing history within 2 weeks before preoperative and postoperative measurement. Patients with IOL misalignment of > 5° by the toriCAM Application (12) and postoperative corrected distance visual acuity (CDVA) of < 20/40 were excluded.
2.2 Preoperative examination
A thorough ophthalmic examination was performed preoperatively, including uncorrected visual acuity and best-corrected visual acuity, slit-lamp examination, subjective and objective refraction, tonometry, and adequate fundoscopy after dilatation. In this study, two anterior segment SS-OCT biometers (IOLMaster 700 and Anterion) were used. Reliable measurements (quality check “green”) were obtained. Pentacam checking was performed in order to determine the regularity of astigmatism.
2.3 IOL power calculation
For postoperative target refraction, 6 eyes were targeted for myopia between −2.0D and −3.0D, and 50 eyes were aimed to emmetropia. The spherical equivalent (SE) and the cylinder power of IOL were determined using Z CALC 2.0 (Carl Zeiss Meditec AG) based on the TK mode (13), in combination with TK, anterior chamber depth and axial length of the IOLMaster 700. It was assumed that the SIA is 0.3D@120° for all cases according to our previous study.
2.4 Surgical technique
Before the surgery, all patients were manually marked the meridian of the incision and alignment axis of the intraocular lens in a sitting position under natural pupil by one experienced surgeon (JWL.), who performed all surgeries under topical anesthesia. Phacoemulsification was carried out with a 2.2 mm clear corneal primary incision, followed by an insertion of an AT TORBI 709 toric IOL (Carl Zeiss Meditec AG, Jena, Germany). After IOL implantation, the reference marker of IOL was aligned to the corneal preoperative marks.
2.5 Postoperative examination and calculations
Postoperative evaluation was performed at 1 month after surgery, including visual acuity and manifest refraction. Predicted residual astigmatism with the implanted toric IOL was calculated using two sets of values: (1) TK from the IOLMaster 700, and (2) TCP from the Anterion. PEs were calculated as the difference between the actual refractive astigmatism postoperatively and the predicted residual astigmatism (14). Vector analysis was carried out in all calculations according to the Holladay method (15). The mean vector PE magnitude, percentage of eyes within certain intervals of vector PE, as well as the centroid PE were calculated (16).
According to the steep meridian of TK obtained by the IOLMaster 700, eyes were divided into WTR group (the steep meridian at 67.5° to 112.5°) and ATR group (the steep meridian at 0° to 22.5° or 157.5° to 180°), and then subgroup analysis was performed in the two groups. Oblique astigmatism (the steep meridian at 22.5° to 67.5° or 112.5° to 157.5°) was excluded due to small sample size.
2.6 Statistical analysis
With a significance level of 5% and a power of 80%, sample size calculation was performed to detect a prediction error of 0.2 D and standard deviation of PE of approximately 0.5 D showed in a preliminary study. The calculation showed that 52 eyes were required.
The main outcomes were vector PE magnitudes, percentage of eyes with vector PEs in certain intervals, and centroid PE of TK and TCP. The R software1 and related functions outlined in the study by Holladay et al. (17) were used for statistical analysis. Univariate normality analysis was performed using Shapiro–Wilk test. The difference in total variance of preoperative corneal and actual postoperative refractive astigmatism for each method were compared. The vector PEs for TK method and TCP method were compared using a bootstrap-t method. McNemar test was carried out to compare the percentage of eyes with vector PEs at various intervals between the two methods. P-values were corrected by Holm-Bonferroni method. A bivariate analysis of x-PE component and y-PE component of the centroid value was conducted according to the bootstrap-t method. P < 0.05 was considered statistically significant.
3 Results
Demographic data of the study population are shown in Table 1. A total of 56 eyes from 56 patients were enrolled in the study.
3.1 Double-angle plots and total variances
Double-angle plots and cumulative histograms of preoperative corneal astigmatism and postoperative refractive astigmatism at corneal plane for TK and TCP were presented in Figure 1. As shown in the figure, there was no statistical difference in the percentages of eyes at various intervals of corneal astigmatism between TK and TCP (± 0.25D to ± 2.0D, with 0.25D interval) before surgery. There were no significant differences in preoperative bivariate variances between TK and TCP (5.1223 versus 5.5019, P = 0.081). The postoperative refractive astigmatism values were significantly decreased compared to both preoperative TK and TCP astigmatism (both P < 0.001 using Holm correction). Statistical significance was found between preoperative and postoperative bivariate variances (5.1223&5.5019 versus 0.3899, both P < 0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1276186/fmed-10-1276186-HTML/image_m/fmed-10-1276186-g001.jpg)
Figure 1. Double-angle plots (A,B,C) and a cumulative histogram (D,E) of the preoperative corneal and actual postoperative refractive astigmatism for IOLMaster 700 and Anterion. (A) TK from IOLMaster 700; (B) TCP from Anterion; (C) postoperative refractive astigmatism.
3.2 Mean vector prediction error
The mean vector magnitude PEs were 0.48D and 0.46D for TK and TCP, respectively, with no statistically significant difference (P = 0.281) (Figure 2). Similarly, no statistical differences were found for percentages of eyes within different intervals of vector PEs (all P > 0.05). Compared with TK, the percentages of eyes with TCP increased by 12.7 and 5.4% within the intervals of vector PEs ≤ 0.25D and ≤ 0.5D, respectively (Table 2); however, these improvements were not statistically significant.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1276186/fmed-10-1276186-HTML/image_m/fmed-10-1276186-g002.jpg)
Figure 2. Double-angle plots of astigmatic prediction error with 95% confidence convex polygon using total keratometry (TK) from IOLMaster and total corneal power (TCP) from Anterion for all eyes (top, A: IOLMaster; B: Anterion), WTR eyes (middle, C: IOLMaster; D: Anterion) and ATR eyes (bottom, E: IOLMaster; F: Anterion). ATR, against-the-rule astigmatism; WTR, with-the-rule astigmatism; prediction error = postoperative refractive astigmatism (spectacle plane) - predicted residual astigmatism (spectacle plane).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1276186/fmed-10-1276186-HTML/image_m/fmed-10-1276186-t002.jpg)
Table 2. The percentage of eyes with vector PE magnitudes below certain intervals and the P-value for comparison (Holm corrected) between TK and TCP.
Subgroup analysis showed that the vector PEs for TK and TCP were 0.52D and 0.51D in the WTR eyes, and 0.46D and 0.43D in the ATR eyes, respectively. There were no significant differences between TK and TCP in WTR eyes (P = 0.593) and ATR eyes (P = 0.245).
3.3 Centroid prediction error
Figure 2 and Table 3 presented double-angle plots of prediction errors for TK and TCP, with a brownish-yellow convex polygon representing 95% confidence boundary of the data. The centroid values were 0.10D@35° ± 0.60D for TK and 0.15D@22° ± 0.57D for TCP. There were no significant differences in x-PE and y-PE components between the two methods.
When subgroup analysis was performed, the WTR eyes had centroid PE of 0.17D@64° ± 0.60D for TK, aligned along the vertical axis, namely, the WTR prediction error; and 0.15D@31° ± 0.61D for TCP. In contrast, the ATR eyes had centroid value of 0.13D@17° ± 0.59D for TK and 0.15D@17° ± 0.55D for TCP, both were ATR prediction errors.
4 Discussion
Accurate astigmatism correction requires accurate measurement and calculation (18). At present, there are several instruments to measure total corneal astigmatism, which improves the accuracy of toric IOL calculation (10, 11). Both IOLMaster 700 and Anterion are based on swept-source optical coherence tomography. Previous studies have shown good repeatability and agreement of biometric measurements between them (19). Different from previous research, we evaluated the accuracy of TK from IOLMaster 700 and TCP from Anterion in astigmatism prediction based on Z CALC total keratometry mode. Strengths of the current study include the evaluation of TCA from two SS-OCT biometers (IOLMaster 700 and Anterion) in astigmatic prediction accuracy instead of PCA, the comparison of bivariate variances of two SS-OCT biometers, the inclusion of one eye per patient, and the use of confidence convex polygons instead of confidence ellipses on Double-angle plots.
No significant difference was found in the bivariate variances of preoperative corneal astigmatism between TK and TCP (5.1223 versus 5.5019, P = 0.081). The bivariate variances between preoperative corneal astigmatism and postoperative refractive astigmatism were compared for each method, and the difference reached statistical significance (5.1223 & 5.5019 versus 0.3899, both P < 0.001). Compared to the preoperative corneal astigmatism measured with TK and TCP, postoperative refractive astigmatism was significantly smaller, and the percentages of eyes within certain amounts of astigmatism magnitudes were significantly higher after surgery. The above outcomes indicated that each method significantly improved the postoperative refractive astigmatism.
The mean vector PEs were comparable between TK and TCP (0.48D versus 0.46D, P = 0.281). Compared with TK, the percentages of eyes with TCP increased by 12.7 and 5.4% with vector PEs of ≤ 0.25D and ≤ 0.5D, respectively. Although there was no statistical difference between the two methods, it did indicate that, if toric IOL calculation was performed using TCP from Anterion, 1 in 10 eyes may obtain more accurate refractive astigmatism results. The reasons for the difference could be as follows: the IOLMaster 700 uses a combination of SS-OCT and telecentric technology to directly measure the corneal curvature at a total of 18 points in the central 1.5, 2.5, and 3.5 mm rings (20). In comparison, the Anterion, based on SS-OCT technology, rapidly captures the corneal curvature of the central 3 mm ring of the cornea. A total of 65 radial B-scan images are obtained in a single scan, 256 A-scans per B-scan, which can theoretically calculate TCP more accurately (21).
Our results were compatible with previously published findings (10). Choi et al. (10) compared the astigmatic prediction accuracy of TCA from IOLMaster 700 and Pentacam when implanting Acrysof IQ Panoptix toric TFNT IOLs. In their study, the PE was ≤ 0.50 D for 45.5% (66/145) eyes when calculating based on TK, in contrast, 42.8% (62/145) eyes were ≤ 0.50 D when TCRP4 was used. For TK and TCRP4, the MAE were 0.60D and 0.61D, respectively (10). Our results are slightly higher than those of Choi et al. (10) possibly for the following three reasons: (1) The toric IOL calculators adopted in these studies were different: Z CALC 2.0 was used in current study, while Choi et al. used Alcon online calculator; (2) Different toric IOLs were implanted: The implanted IOL in this study was ZEISS 709M toric IOL, and Choi et al. inserted Acrysof IQ Panoptix toric TFNT IOL. It was showed that different designs of toric IOLs may result in different ELPs, which may affect the postoperative results (22). (3) Different biometers were utilized in the two studies.
In the current study, the x-PE component and y-PE component for TK and TCP were not normally distributed, and 95% confidence convex polygons were performed. No significant differences were found in x-PE component and y-PE component between TK and TCP. In subgroup analysis, the centroid PEs of TK and TCP in the WTR group were 0.17D@64° ± 0.60D and 0.15D@31° ± 0.61D, respectively, and those in the ATR group were 0.13D@17° ± 0.59D and 0.15D@17° ± 0.55D, respectively. This indicated that the IOLMaster 700 may have overestimated PCA in the WTR eyes and underestimated PCA in the ATR eyes, and the Anterion may have underestimated PCA in all eyes. Since the sample size in each group was not large enough and the SD values were large, subsequent studies need to include more eyes to draw a conclusion about whether these differences were in fact significant.
The limitation of our study is that the total corneal curvature measured preoperatively was used in the toric IOL calculation, and the SIA may have an impact on the analysis. The SIA used in the study is the average personalized SIA previously calculated for the same surgeon. Clinical practice can be simulated by using preoperative parameters and a personalized SIA. In addition, the sample size of the study was relatively small, only 20 WTR eyes were included, and oblique astigmatism was not included. Therefore, a larger population study is needed to evaluate the role of TCA in toric IOL calculation. Furthermore, further studies are needed to be extrapolated to other similar devices.
5 Conclusion
This study showed that, compared with TK from the IOLMaster 700, TCP from the Anterion achieved comparable accuracy in residual astigmatism prediction and resulted in 12.7 and 5.4% more eyes with vector PEs of ≤ 0.25D and ≤ 0.5D, respectively. This demonstrates that Anterion, as a new SS-OCT biometer, may provide a new choice for astigmatism prediction.
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 Shanxi Eye Hospital Institutional Review Board (No. SXYYLL-20210107). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
FC: Conceptualization, Data curation, Funding acquisition, Writing – original draft, Writing – review & editing. JM: Conceptualization, Software, Writing – review & editing, Methodology. XW: Conceptualization, Funding acquisition, Writing – original draft, Writing – review & editing. JiL: Formal analysis, Investigation, Validation, Writing – review & editing, Resources. LJ: Formal analysis, Validation, Visualization, Writing – review & editing, Investigation. WW: Data curation, Validation, Writing – review & editing, Formal analysis, Investigation. JuL: Conceptualization, Project administration, Supervision, Writing – original draft, Writing – review & editing, Methodology.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (81971697) and Scientific Research Foundation of Shanxi Eye Hospital (C201902). The funder had no role in the study, data analysis, manuscript writing, or publication.
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.
Footnotes
References
1. Goggin M, Moore S, Esterman A. Toric intraocular lens outcome using the manufacturer’s prediction of corneal plane equivalent intraocular lens cylinder power. Arch Ophthalmol. (2011) 129:1004–8. doi: 10.1001/archophthalmol.2011.178
2. Hirnschall N, Hoffmann PC, Draschl P, Maedel S, Findl O. Evaluation of factors influencing the remaining astigmatism after toric intraocular lens implantation. J Refract Surg. (2014) 30:394–400. doi: 10.3928/1081597X-20140429-01
3. Potvin R, Kramer BA, Hardten DR, Berdahl JP. Toric intraocular lens orientation and residual refractive astigmatism: an analysis. Clin Ophthalmol. (2016) 10:1829–36. doi: 10.2147/OPTH.S114118
4. Savini G, Naeser K. An analysis of the factors influencing the residual refractive astigmatism after cataract surgery with toric intraocular lenses. Invest Ophthalmol Vis Sci. (2015) 56:827–35. doi: 10.1167/iovs.14-15903
5. Koch DD, Jenkins RB, Weikert MP, Yeu E, Wang L. Correcting astigmatism with toric intraocular lenses: effect of posterior corneal astigmatism. J Cataract Refract Surg. (2013) 39:1803–9. doi: 10.1016/j.jcrs.2013.06.027
6. Ho JD, Tsai CY, Liou SW. Accuracy of corneal astigmatism estimation by neglecting the posterior corneal surface measurement. Am J Ophthalmol. (2009) 147(5):788–95, 795.e1-2. doi: 10.1016/j.ajo.2008.12.020
7. Tonn B, Klaproth OK, Kohnen T. Anterior surface-based keratometry compared with Scheimpflug tomography-based total corneal astigmatism. Invest Ophthalmol Vis Sci. (2014) 56:291–8. doi: 10.1167/iovs.14-15659
8. Koch DD, Ali SF, Weikert MP, Shirayama M, Jenkins R, Wang L. Contribution of posterior corneal astigmatism to total corneal astigmatism. Cataract Refract Surg. (2012) 38:2080–7. doi: 10.1016/j.jcrs.2012.08.036
9. Reitblat O, Levy A, Kleinmann G, Abulafia A, Assia EI. Effect of posterior corneal astigmatism on power calculation and alignment of toric intraocular lenses: comparison of methodologies. J Cataract Refract Surg. (2016) 42:217–25. doi: 10.1016/j.jcrs.2015.11.036
10. Choi A, Kwon H, Jeon S. Accuracy of total corneal power calculation for multifocal toric intraocular lens implantation: swept-source OCT-based biometer vs scheimpflug tomographer. J Refract Surg. (2021) 37:686–92. doi: 10.3928/1081597X-20210610-01
11. Wang L, Koch DD. Comparison of accuracy of a toric calculator with predicted versus measured posterior corneal astigmatism. J Cataract Refract Surg. (2023) 49:29–33. doi: 10.1097/j.jcrs.0000000000001025
12. Pallas A, Yeo TK, Trevenen M, Barrett G. Evaluation of the accuracy of two marking methods and the novel toriCAM application for toric intraocular lens alignment. J Refract Surg. (2018) 34:150–5. doi: 10.3928/1081597X-20180115-03
13. Lesieur G. Microincision cataract surgery with implantation of a bitoric intraocular lens using an enhanced program for intraocular lens power calculation. Eur J Ophthalmol. (2020) 30:1308–13. doi: 10.1177/1120672119877767
14. Gundersen KG, Potvin R. Clinical outcomes with toric intraocular lenses planned using an optical low coherence reflectometry ocular biometer with a new toric calculator. Clin Ophthalmol. (2016) 10:2141–7. doi: 10.2147/OPTH.S120414
15. Holladay JT, Moran JR, Kezirian GM. Analysis of aggregate surgically induced refractive change, prediction error, and intraocular astigmatism. J Cataract Refract Surg. (2001) 27:61–79. doi: 10.1016/s0886-3350(00)00796-3
16. Abulafia A, Koch DD, Holladay JT, Wang L, Hill W. Pursuing perfection in intraocular lens calculations: IV. Rethinking astigmatism analysis for intraocular lens-based surgery: suggested terminology, analysis, and standards for outcome reports. J Cataract Refract Surg. (2018) 44:1169–74. doi: 10.1016/j.jcrs.2018.07.027
17. Holladay JT, Wilcox RR, Koch DD, Wang L. Astigmatism analysis and reporting of surgically induced astigmatism and prediction error. J Cataract Refract Surg. (2022) 48:799–812. doi: 10.1097/j.jcrs.0000000000000871
18. Browne AW, Osher RH. Optimizing precision in toric lens selection by combining keratometry techniques. J Refract Surg. (2014) 30:67–72. doi: 10.3928/1081597X-20131217-07
19. Cheng SM, Zhang JS, Shao X, Wu ZT, Li TT, Wang P, et al. Repeatability of a new swept-source optical coherence tomographer and agreement with other three optical biometers. Graefes Arch Clin Exp Ophthalmol. (2022) 260:2271–81. doi: 10.1007/s00417-022-05579-9
20. Hoffer KJ, Hoffmann PC, Savini G. Comparison of a new optical biometer using swept-source optical coherence tomography and a biometer using optical low-coherence reflectometry. J Cataract Refract Surg. (2016) 42:1165–72. doi: 10.1016/j.jcrs.2016.07.013
21. Kim KY, Choi GS, Kang MS, Kim US. Comparison study of the axial length measured using the new swept-source optical coherence tomography ANTERION and the partial coherence interferometry IOL Master. PLoS One. (2020) 15:e0244590. doi: 10.1371/journal.pone.0244590
Keywords: total keratometry, total corneal power, prediction error, IOLMaster 700, Anterion
Citation: Chai F, Ma JX, Wang X, Liu J, Jiang L, Wu W and Li J (2024) Total keratometry for toric intraocular lens calculation: comparison from two swept-source optical coherence tomography biometers. Front. Med. 10:1276186. doi: 10.3389/fmed.2023.1276186
Received: 11 August 2023; Accepted: 19 December 2023;
Published: 08 January 2024.
Edited by:
Georgios A. Kounis, GNEMS, GreeceReviewed by:
Francesco D’Oria, Azienda Ospedaliero Universitaria Consorziale Policlinico di Bari, ItalyPablo De Gracia, University of Detroit Mercy, United States
Copyright © 2024 Chai, Ma, Wang, Liu, Jiang, Wu and Li. 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: Junhong Li, junhongli2022@163.com