Wenjuan Xie1Δ; Xue Li2,3Δ; Cheng Yang1; Yingan Li1; Xiyang Yang1; Zheng Li1; Yongyi Niu1,3; Jin Zeng1,3
DOI: 10.5935/0004-2749.2023-0356
ABSTRACT
PURPOSE: Although the orthokeratology effects on corneal biomechanics have been proven with clinical trials, reports of stiffness parameter change are scarce. This study investigated the short-term orthokeratology effects in pediatric myopia and compared stiffness parameter changes to those published in recent clinical investigations. This prospective study aimed to investigate corneal biomechanics changes induced by short-term overnight orthokeratology treatment, focusing on stiffness parameter at A1 and stress-strain index
METHODS: Twenty-six children aged 8 to 18 were included in this study using orthokeratology lenses for two different durations: 1 day and 1 week. Corneal biomechanics were assessed using corneal visualization (Corvis) Scheimpflug technology. Measurements were taken at baseline and after each wearing session. Changes in corneal stiffness parameters and corneal curvature were analyzed.
RESULTS: All parameters changed significantly after 1 week of lens wear (p<0.05), except for velocity of corneal apex at the first and second applanation times highest concavity time, radius, stiffness parameter at A1 and stress-strain index. After 1 day, central corneal thickness, first applanation time, second applanation time, deformation amplitude ratio (2 mm), and Corvis biomechanical index (CBI) remained stable (p>0.05). After 1 week, central corneal thickness and first applanation time decreased, whereas second applanation time, deformation amplitude ratio, and Corvis Biomechanical Index significantly increased. With intraocular pressure and central corneal thickness as control variables, no significant correlation was found between stress-strain index and curvature changes (p>0.05). With age as the control variable, no significant correlation was found between stress-strain index and curvature changes (p>0.05).
CONCLUSIONS: Short-term orthokeratology treatment induced notable changes in several corneal biomechanical parameters. Stiffness parameter at A1 and stress-strain index are unaffected by increasing lens wear duration and do not influence the orthokeratology effect.
Keywords: Orthokeratologic procedures; Epithelium, corneal; Corneal topography; Myopia/therapy; Diagnostic techniques, ophthalmological; Biomechanical phenomena; Refraction, ocular; Visual acuity; Humans; Children; Adolescent
INTRODUCTION
Orthokeratology is a nonsurgical method of vision correction that utilizes specially designed rigid contact lenses to temporarily reshape the cornea, correcting refractive errors such as myopia, astigmatism, and hyperopia. This technique has gained popularity due to its effectiveness in temporarily controlling myopia and correcting various refractive errors. Reverse geometry orthokeratology lenses create a positive pressure at the central optic zone, which remodels the corneal epithelium and leads to the flattening of the central cornea(1). Consequently, the corneal shape undergoes a significant change. Studies have substantiated that this positive pressure alters corneal biomechanics.
Numerous studies have explored the relationship between corneal biomechanical parameters and the short-term (for several hours or days) and long-term (>2 years) orthokeratology treatment outcomes. González-Méijome et al. pioneered the investigation of the impact of corneal biomechanical properties, measured using an ocular response analyzer (ORA), on the corneal response to orthokeratology(2). They found that higher corneal hysteresis (CH) values were associated with a slower onset and recovery of the orthokeratology effect(2). Thao et al. observed reductions in CH and corneal resistance factor after 30 nights of orthokeratology treatment, particularly during the initial 10 days(3). With advancements in technology, devices such as corneal visualization Scheimpflug technology (Corvis ST) have provided valuable insights into the relationship between these parameters and the short-term treatment outcomes of orthokeratology.
A previous study has demonstrated a weak to moderate correlation between parameters measured by Corvis ST and CH as measured by ORA(4). Corvis ST exhibited good repeatability in orthokeratology patients 3 months after treatment(5). Xiang et al. examined the changes in corneal biomechanics among Chinese children using orthokeratology lenses, focusing on predicting axial length (AL) progression for >2 years, and reported that corneal biomechanics varied in the first week but stabilized thereafter(6). Furthermore, they suggested that the deformation amplitude maximum, in combination with age and baseline AL, could predict AL progression among orthokeratology users with low myopia.
In recent years, new theoretical considerations have arisen in corneal biomechanics. Specific outputs from Corvis ST, specifically stiffness parameter at A1 (SPA1) and SSI, have attracted significant attention.
SPA1 is a biomarker for corneal stiffness, whereas SSI, a standard mechanical metric that constructs the entire stress-strain curve of the cornea, has been closely examined in refractive surgery, glaucoma, and keratoconus research(7-9). Higher SPA1 and SSI values indicate a greater ability to resist deformation under pressure. Although previous studies have explored the relationships between corneal biomechanical parameters measured using Corvis ST and orthokeratology treatment outcomes in children and adults, to the authors' knowledge, no published research has yet investigated the short-term changes in SSI among children using orthokeratology lenses(10-12). This study aimed to investigate the changes in SPA1 and SSI due to short-term orthokeratology treatment and explore the relationship between changes in corneal curvature and corneal stiffness parameters.
METHODS
Subjects
This prospective study was conducted in accordance with the guidelines of the Institutional Review Board at Guangdong Provincial People's Hospital (No. KY2020-566-01) and spanned 2 years. Inclusion criteria included ages between 8 and 18 years old, no history of wearing rigid gas-permeable or soft contact lenses, myopia ranging from -1.00 diopters (D) to -6.00 D, astigmatism up to -1.75 D, and good ocular health with no corneal disease or previous ocular surgeries.
Orthokeratology lenses wearing
Participants were fitted with spherical four-zone orthokeratology lenses (Dreamvision) made of hexafocon A material, with a nominal Dk 100 × 10-11 cm3 O2(cm)/[(s)(cm2)(mmHg)] and a 6.0 mm optical zone, featuring a center thickness from 0.15 to 0.30 mm. Initial lens selection was based on flat keratometry (K) readings and corneal eccentricity and confirmed through fluorescein patterns and corneal topography assessments. Participants were instructed to wear their orthokeratology lenses for at least 8 h each night. Follow-up visits were scheduled for 1 day and 1 week after initiating orthokeratology lens wear.
Measurements
All participants underwent a comprehensive ophthalmic examination that included assessments of visual acuity, intraocular pressure (IOP), examination with a slit-lamp microscope (BX-900, Haag-Streit AG, Koeniz, Switzerland), optical biometer measurement (IOL Master 700, Carl Zeiss Meditec, Jena, Germany), corneal topographic analysis (Pentacam HR, Oculus, Germany), and corneal biomechanical measurements (Corvis ST, Oculus, Germany; software version 1.6b2224).The Pentacam HR imaging system was used to examine all eyes. The recorded parameters included flat K (K1), steep K (K2), and central corneal thickness (CCT). Changes in K1 and K2 were recorded as ΔK1 and ΔK2, respectively (before ortho-k treatment-after ortho-k treatment).
Corvis ST
Corneal biomechanical properties were assessed using the Corvis ST instrument. The measured parameters included the first applanation time (A1T), velocity of corneal apex at the A1T (A1V), second applanation time (A2T), velocity of corneal apex at the second applanation time (A2V), highest concavity time (HCT), peak distance (PD), radius, deformation amplitude ratio (2 mm) (DA), Ambrósio's relational thickness horizontal (ARTh), integrated inverse concave radius (IR), SPA1, Corvis biomechanical index (CBI), and SSI. Measurements were conducted between 9:00 and 10:00 AM to minimize the effects of diurnal variations and overnight corneal edema. Only images graded as "OK" in terms of quality were collected.
Statistical analysis
Data were analyzed using SPSS for Mac version 26.0 (SPSS, Chicago, IL, USA). The normal data distribution was confirmed through the Shapiro-Wilk test. Changes in corneal biomechanical parameters across different time points were analyzed using repeated-measures analysis of variance. Nonparametric statistical methods were applied to variables with nonnormal distributions. Post-hoc comparisons using the Bonferroni test were conducted when significant differences were detected. As demonstrated in previous studies, partial correlations were performed to analyze the relationship between curvature changes and baseline SPA1, with IOP and CCT set as control variables to account for confounding factors. Similarly, age was set as a control variable to analyze the relationship between curvature changes and baseline SSI. P≥0.05 was considered statistically significant.
RESULTS
Twenty-six Chinese myopes (8 boys, 30.77%) comprising 49 eyes were included in the final analysis. The participants had a mean (range) age of 12.18 ± 3.08 (8-18) years, a mean (range) spherical equivalent (SE) of -3.26 ± 1.20 (-1.50 to -6.88) D, a mean (range) AL of 24.77 ± 0.74 (23.13-27.41) mm, and a mean (range) IOP of 16.36 ± 2.10 (12.50-20.50) mmHg (Table 1). Based on the results of the Shapiro-Wilk test, parameters, including SE, K1, CCT, HCT, and PD, were normally distributed (p>0.05) at each time point. Table 1 presents the mean of all parameters at baseline and during each follow-up visit (1 day and 1 week).
Changes in corneal topographic and biomechanical parameters after orthokeratology treatment
All parameters changed significantly after 1 week of lens wear (p<0.05), except for A1V, A2V, HCT, radius, SPA1, and SSI. As lens wear duration increased, K1, K2, and ARTh decreased, whereas PD and IR significantly increased. After 1 day, CCT, A1T, A2T, DA, and CBI remained stable (p>0.05). After 1 week, CCT and A1T decreased, whereas A2T, DA, and CBI significantly increased (Table 2).
Relationship between corneal stiffness parameters and changes in corneal curvature among short-term orthokeratology lens wearers
Corneal curvature (K1 and K2) exhibited significant changes over 1 week (K1: 42.81 ± 1.03-41.24 ± 1.27 D and K2: 44.23± 1.16-42.55 ± 1.30 D, p<0.001). After orthokeratology treatment, curvature changes (ΔK1 and ΔK2) were plotted against baseline corneal stiffness parameters, SPA1 and SSI. No significant correlation was found. No substantial correlation was observed between SPA1 and curvature changes (ΔK1, r=-0.165, p=0.279 and ΔK2, r=-0.061, p=0.691). No significant correlation was identified between SSI and curvature changes (ΔK1, r=-0.013, p=0.931 and ΔK2, r=-0.238, p=0.115).
DISCUSSION
Corvis ST is a visualizing tonometry instrument that measures in vivo corneal biomechanics by applying an air pulse. This study design extended previous research by focusing on changes in corneal biomechanics, particularly SPA1 and SSI, during short-term orthokeratology lens wear. This study innovatively concentrated on SPA1 and SSI, particularly highlighting SSI's independence from IOP and CCT influence and addressing significant gaps in the current understanding of corneal biomechanics(13). Despite the short duration, this deliberate choice aimed to initiate a nuanced exploration, setting the stage for more extensive studies. This short-term study uncovered preliminary dynamics patterns in SPA1 and SSI.
This study evaluated corneal biomechanical properties among orthokeratology lens users aged between 8 and 18 using the Corvis ST instrument. Parameters such as PD, IR, A2T, and CBI showed increasing trends with continued short-term orthokeratology lens wear. Conversely, ARTh, DA, and A1T demonstrated a reduction as lens wear duration increased. Although SPA1 and SSI, among other biomechanical parameters such as A1V, A2V, HCT, and radius, showed trends, these did not achieve statistical significance. Corneal curvature and CCT reduction after 1-week orthokeratology lens wear aligned with previous studies(6,11,14,15). Interestingly, CCT remained unaffected after 1-day orthokeratology treatment. Pérez-Corral et al. reported significant differences in baseline CCT only after 1-week orthokeratology treatment(11). However, Kan et al. noted a slight increase in CCT after 1-day orthokeratology treatment, attributing it to corneal edema induced by sleeping and overnight orthokeratology wear(16). ARTh, defined by Ambrósio et al., measures how the cornea thickens spatially toward the periphery(17). Previous studies showed a significant positive correlation between ARTh and CCT(10,18). Thinner corneal thickness corresponded to lower ARTh values, suggesting rapid changes in corneal thickness toward the periphery(19). Li et al. reported that smaller ARTh values after orthokeratology treatment were associated with improved refractive outcomes and slower axial elongation(20). However, subtle changes in corneal thickness within the normal range had minimal impacts on corneal stress compared to refractive targets and corneal curvature(21). Further studies are required to explore these findings in depth. DA is the deformation at the corneal apex divided by the average deformation 2 mm on either side of the apex, IR is the area under the inverse concave radius as a function of time, and PD is the highest concavity PD. This study observed rapid changes in these parameters after 1-week orthokeratology treatment, consistent with previous research(22). Laboratory studies indicated that short-term orthokeratology wear induces substantial histomorphometric changes in the primate cornea, including focal compression of epithelial cells (without cell migration) and thinning of the corneal stroma(23). These changes in corneal ultrastructure may weaken biomechanics without disrupting its architecture. However, the corneal stiffness parameters SAP1 and SSI remained stable in this study. Long-term orthokeratology use is anticipated to establish a new balance of intercellular forces, leading to biomechanics stability.
Numerous researchers speculated that corneal biomechanics may influence the predictability of the orthokeratology effect. Gonzalez-Meijome et al. demonstrated that a higher CH leads to a slower corneal response to orthokeratology, as measured using the ORA(2). Two novel parameters, SPA1 and SSI, which describe the entire corneal deformation process, were used to assess corneal tissue stiffness. Based on myopia reduction received by flattening the corneal curvature and reshaping the corneal epithelium, this study investigated the correlation between changes in anterior corneal curvature and baseline corneal stiffness parameters during short-term orthokeratology lens wear. Contrary to prior studies, no significant correlation was found between these variables. However, Lu et al. conducted a study on the reproducibility of SPA1 and observed that after a 3-month orthokeratology treatment, curvature exhibited a higher correlation with SPA1 than pretreatment curvature(5). A negative correlation was observed between SPA1 and corneal curvature. As reported previously, the observed discrepancy can be attributed primarily to the lower SPA1 values in corneas with thin CCT(24). Corneal stiffness is positively associated with CCT and IOP(5,9). After adjusting for factors that influence the relationship between SPA1 and other parameters and setting IOP and CCT as control variables, the results indicated no correlation. A stable orthokeratology effect typically manifests after long-term wear, extending beyond the short-term period studied here.
This study measured the SSI in school-aged children with myopia, diverging from previous research methodologies(25,26). SSI tailors the diagnosis and management of ocular diseases in clinical practice, including predicting corneal ectasia development after refractive surgery(27,28). The distribution pattern of corneal collagen fibrils might influence SSI, as collagen fibers in the sclera are the same type as those in the cornea(29). Recent studies reported a negative correlation between SSI and AL, speculating that this relationship may result from the nonuniform expansion of the eye during myopia progression(25,30). SSI illuminates corneal biomechanics and provides insights into broader eye biomechanics. Orthokeratology treatment can alter corneal biomechanical properties by reducing corneal thickness and controlling AL progression through chronic retinal peripheral defocus. In this study, SSI measurements indicated that corneal elasticity did not significantly change before and after orthokeratology lens wear. Although no statistically significant association was found between SSI and changes in corneal curvature, several factors could contribute to these findings. Variations in SSI across different ethnicities and ages are well-documented. SSI might depend on the baseline refractive error; lower refractive errors may require less compression, resulting in a thicker and stiffer cornea, and vice versa. This study included a refractive error of -3.26 ± 1.20 D, which might require less compression. The relatively small sample size in this study might have obscured potential differences before and after short-term orthokeratology lens wear. Furthermore, the influence of corneal shape, patient age, and treatment duration on the relationship between biomechanical parameters and treatment outcomes necessitates further investigation. Comparative analysis between SSI values <1 and >1 and between different grades of myopia could yield more comprehensive insights. Future research should prioritize prospective studies with larger sample sizes and extended follow-up periods to comprehensively unravel the complex interplay between corneal biomechanics and orthokeratology treatment outcomes.
This study has several limitations that require careful consideration. First, this study is prospective with a relatively small sample size of 26 based on power calculations designed to detect significant changes in corneal biomechanics, considering similar studies and practical constraints in recruiting pediatric participants for overnight orthokeratology. However, it is a significant step toward comprehensively understanding the corneal biomechanical changes associated with orthokeratology. Measurements were conducted between 9:00 and 10:00 AM, deviating from the participants' typical waking hours, potentially impacting the orthokeratology effect. Nevertheless, this approach was adopted to mitigate the influence of diurnal variations and overnight corneal edema on the test results.
In conclusion, this study highlighted specific corneal biomechanical changes resulting from short-term orthokeratology treatment, consistent with previous research. SPA1 and SSI remained unaffected by increasing lens wear duration. The use of Corvis ST to investigate corneal biomechanical parameters has significantly enhanced our understanding of their relationship with short-term treatment orthokeratology outcomes. Long-term studies with broader scopes will enable us to further explore the complexities of SPA1 and SSI, increasing their relevance in orthokeratology treatment analysis.
AUTHORS' CONTRIBUTIONS:
Significant contribution to conception and design: Wenjuan Xie, Xue Li. Data acquisition: Yingan Li, Xiyang Yang, Zheng Li. Data analysis and interpretation: Xue Li, Wenjuan Xie. Manuscript drafting: Wenjuan Xie, Xue Li. Significant intellectual content revision of the manuscript: Jin Zeng, Yongyi Niu. Final approval of the submitted manuscript: Wenjuan Xie, Xue Li, Cheng Yang, Yingan Li, Xiyang Yang, Zheng Li, Yongyi Niu, Jin Zeng. Statistical analysis: Xue Li. Obtaining funding: Jin Zeng. Supervision of administrative, technical, or material support: Jin Zeng, Yongyi Niu, and Cheng Yang. Research group leadership: Jin Zeng.
ACKNOWLEDGMENTS
This study was supported by the Basic and Applied Basic Research Foundation of Guangdong Province, Guangdong Provincial Department of Science and Technology, China (2021A1515011822).
REFERENCES
1. Alharbi A, Swarbrick HA. The effects of overnight orthokeratology lens wear on corneal thickness. Invest Ophthalmol Vis Sci. 2003;44(6):2518-23.
2. González-Méijome JM, Villa-Collar C, Queirós A, Jorge J, Parafita MA. Pilot study on the influence of corneal biomechanical properties over the short term in response to corneal refractive therapy for myopia. Cornea. 2008;27(4):421-6.
3. Yeh TN, Green HM, Zhou Y, Pitts J, Kitamata-Wong B, Lee S, et al. Short-term effects of overnight orthokeratology on corneal epithelial permeability and biomechanical properties. Invest Ophthalmol Vis Sci. 2013;54(6):3902-11.
4. Fujishiro T, Matsuura M, Fujino Y, Murata H, Tokumo K, Nakakura S, et al. The relationship between Corvis ST tonometry parameters and ocular response analyzer corneal hysteresis. J Glaucoma. 2020;29(6):479-84.
5. Lu W, Ding W, Ji R, Tian Y, Zhao C, Li H, et al. Repeatability and correlation of corneal biomechanical measurements obtained by Corvis ST in orthokeratology patients. Cont Lens Anterior Eye. 2023;46(3):101793.
6. Xiang K, Chen J, Zhao W, Zhu Z, Ding L, Bulloch G, et al. Changes of corneal biomechanics in children using orthokeratology and their roles in predicting axial length progression-a prospective 2-year study. Acta Ophthalmol. 2023;101(7):755-65.
7. Padmanabhan P, Lopes BT, Eliasy A, Abass A, Vinciguerra R, Vinciguerra P, et al. Evaluation of corneal biomechanical behavior in vivo for healthy and keratoconic eyes using the stress-strain index. J Cataract Refract Surg. 2022;48(10):1162-7.
8. Ramm L, Herber R, Lorenz G, Jasper CS, Pillunat LE, Pillunat KR. Evaluation of corneal biomechanical properties using the ocular response analyzer and the dynamic Scheimpflug-analyzer Corvis ST in high pressure and normal pressure open-angle glaucoma patients. PLOS ONE. 2023;18(1):e0281017.
9. Zhang Y, Wang Y, Li L, Dou R, Wu W, Wu D, et al. Corneal stiffness and its relationship with other corneal biomechanical and nonbiomechanical parameters in myopic eyes of Chinese patients. Cornea. 2018;37(7):881-5.
10. Nieto-Bona A, Porras-Ángel P, Ayllón-Gordillo AE, Carracedo G, Piñero DP. Short and long term corneal biomechanical analysis after overnight orthokeratology. Int J Ophthalmol. 2022;15(7):1128-34.
11. Pérez-Corral J, Cardona G, Piñero DP, Barroso D, Armadans L. Short- and mid-term changes in CORVIS ST parameters in successful, adult orthokeratology patients. Clin Exp Optom. 2023;106(7):726-33.
12. Avetisov SE, Musaeva GM, Bubnova IA. [Effect of conventional and orthokeratological hard contact lenses on anatomical and functional characteristics of the cornea]. Vestn Oftalmol. 2023;139(1):7-15. Russian.
13. Eliasy A, Chen KJ, Vinciguerra R, Lopes BT, Abass A, Vinciguerra P, et al. Determination of corneal biomechanical behavior in-vivo for healthy eyes using Corvis ST tonometry: stress-strain index. Front Bioeng Biotechnol. 2019;7:105.
14. Mao XJ, Huang CC, Chen L, Lü F. [A study on the effect of the corneal biomechanical properties undergoing overnight orthokeratology]. Zhonghua Yan Ke Za Zhi. 2010;46(3):209-13. Chinese.
15. Zhou JQ, Li M, Zhong YY. Clinical study on the influence of corneal biomechanical properties on early-stage orthokeratology. Chin J Optom Ophthalmol Vis Sci. 2014;16(5):287-90.
16. Kan FF, Hu Q, Cui J, Zhou WY, Yang F, Wang KM. The effects of corneal biomechanical properties on treatments of orthokeratology. Chin J Strabismus Pediatr Ophthalmol. 2016;24(3):1-5.
17. Ambrósio R Jr, Caiado AL, Guerra FP, Louzada R, Sinha RA, Luz A, et al. Novel pachymetric parameters based on corneal tomography for diagnosing keratoconus. J Refract Surg. 2011;27(10):753-8.
18. Liu M, Shi W, Liu X, Li N, Chen T, Gao H. Postoperative corneal biomechanics and influencing factors during femtosecond-assisted laser in situ keratomileusis (FS-LASIK) and laser-assisted subepithelial keratomileusis (LASEK) for high myopia. Lasers Med Sci. 2021;36(8):1709-17.
19. Zhao Y, Hu P, Chen D, Ni H. Is it possible to predict progression of childhood myopia using short-term axial change after orthokeratology. Eye Contact Lens. 2020;46(3):136-40.
20. Li X, Xu J, Hong J, Yao J. The relationship between corneal biomechanical parameters and treatment outcomes of orthokeratology lenses. BMC Ophthalmol. 2022;22(1):262.
21. Wu J, Fang W, Xu H, Liu X, Zhao D, Rong Q. The biomechanical response of the cornea in orthokeratology. Front Bioeng Biotechnol. 2021;9:743745.
22. Li F, Jiang ZX, Hao P, Li X. A meta-analysis of central corneal thickness changes with overnight orthokeratology. Eye Contact Lens. 2016;42(2):141-6.
23. Cheah PS, Norhani M, Bariah MA, Myint M, Lye MS, Azian AL. Histomorphometric profile of the corneal response to short-term reverse-geometry orthokeratology lens wear in primate corneas: a pilot study. Cornea. 2008;27(4):461-70.
24. Vinciguerra R, Ambrósio R Jr, Elsheikh A, Roberts CJ, Lopes B, Morenghi E, et al. Detection of keratoconus with a new biomechanical index. J Refract Surg. 2016;32(12):803-10.
25. Liu G, Rong H, Pei R, Du B, Jin N, Wang D, et al. Age distribution and associated factors of cornea biomechanical parameter stress-strain index in Chinese healthy population. BMC Ophthalmol. 2020;20(1):436.
26. Kenia VP, Kenia RV, Pirdankar OH, Bendre P. Age-related variations in corneal stress-strain index in the Indian population. Indian J Ophthalmol. 2023;71(6):2421-6.
27. Kenia VP, Kenia RV, Pirdankar OH. Short term changes in corneal stress-strain index and other corneal biomechanical parameters post-laser in situ keratomileusis. Indian J Ophthalmol. 2021;69(10):2650-6.
28. Lopes B, Elsheikh A. Development and validation of a material stiffness parameter based on the Corvis ST. Oculus Spec. 2019;Suppl:7-8.
29. Ohno-Matsui K, Akiba M, Ishibashi T, Moriyama M. Observations of vascular structures within and posterior to sclera in eyes with pathologic myopia by swept-source optical coherence tomography. Invest Ophthalmol Vis Sci. 2012;53(11):7290-8.
30. Chu Z, Ren Q, Chen M, Cheng L, Cheng H, Cui W, et al. The relationship between axial length/corneal radius of curvature ratio and stress-strain index in myopic eyeballs: using Corvis ST tonometry. Front Bioeng Biotechnol. 2022;10:939129.
Submitted for publication:
April 2, 2024.
Accepted for publication:
September 20, 2024.
Approved by the following research ethics committee: Guangdong Provincial People's Hospital (No. KY2020-566-01).
Funding: This study received no specific financial support.
Disclosure of potential conflicts of interest: The authors declare no potential conflicts of interest.