- 1Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
- 2Department of Epidemiology and Biostatistics, School of Public Health, Peking University, Beijing, China
- 3Medical Informatics Center, Peking University, Beijing, China
- 4X Lab, the Second Academy of CASIC, Beijing, China
- 5School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen, China
Orthokeratology has been widely used to control myopia, but the mechanism is still unknown. To further investigate the underlying mechanism of corneal reshaping using orthokeratology lenses via the finite element method, numerical models with different corneal curvatures, corneal thicknesses, and myopia reduction degrees had been developed and validated to simulate the corneal response and quantify the changes in maximum stress in the central and peripheral corneal areas during orthokeratology. The influence of the factors on corneal response had been analyzed by using median quantile regression. A partial eta squared value in analysis of variance models was established to compare the effect size of these factors. The results showed central and peripheral corneal stress responses changed significantly with increased myopia reduction, corneal curvature, and corneal thickness. The target myopia reduction had the greatest effect on the central corneal stress value (partial eta square = 0.9382), followed by corneal curvature (partial eta square = 0.5650) and corneal thickness (partial eta square = 0.1975). The corneal curvature had the greatest effect on the peripheral corneal stress value (partial eta square = 0.5220), followed by myopia reduction (partial eta square = 0.2375) and corneal thickness (partial eta square = 0.1972). In summary, the biomechanical response of the cornea varies significantly with the change in corneal conditions and lens designs. Therefore, the orthokeratology lens design and the lens fitting process should be taken into consideration in clinical practice, especially for patients with high myopia and steep corneas.
Introduction
As one of the most effective methods of myopia control, orthokeratology (ortho-k) has been widely used around the world in the past two decades (Wolffsohn et al., 2016; Cho and Tan, 2019; Jonas et al., 2021). The ortho-k lens, a specially designed rigid contact lens, is different from ordinary optical contact lenses because it can change the corneal shape to achieve the purpose of clear vision during the daytime when the lens is removed after overnight wear. The strong effect on myopia control of ortho-k lens for children and adolescents has been proved and reported by numerous clinical researchers. The investigations revealed a slight variance in myopia control effect with an unweighted average 43% (Cho and Cheung, 2012; Hiraoka et al., 2012; Cheung and Cho, 2016; Cho and Cheung, 2017).
Currently, the clinician fits the ortho-k lens based on the patient’s examination results, which include refractive target degree, corneal curvature, corneal e value, and corneal diameter, under the guidance of the manufacturer and clinical experience. Although the effectiveness in reducing myopia and controlling myopia progression of ortho-k lens has been proven by statistical data, the results vary for the individual subjects (Cho et al., 2005). Furthermore, the mechanism by which it controls the development of myopia, the influences, and the weight of these influencing factors are still unknown, while the clinical application of the ortho-k lens has been reported extensively.
Previous researchers had mentioned that the ortho-k lens plays its role mainly via the mechanical interaction between the lens and the cornea or the eyeball (Swarbrick, 2004; Wan et al., 2020). Numerous studies (Kobayashi et al., 2008; Lau et al., 2019; Mohidin et al., 2020; Queiros et al., 2020) had attempted to reveal the mechanical effect of wearing ortho-k lenses on the cornea by comparing the change of refractive states, corneal topographies, and corneal thicknesses before and after ortho-k treatment, but few researchers went into details of the analysis due to the complexity of measurement and calculation.
In the current study, finite element modeling was used to analyze the corneal stress responses produced by changing the ortho-k lens parameters and fitting the lenses on various corneas. The purpose of our research was to establish a deeper understanding of the mechanical response of the cornea during ortho-k treatment. The interaction between the lenses and the cornea had been simulated, and the stress responses of the cornea with different designs of ortho-k lens had been analyzed quantitatively. By investigating the corneal stress response, our research strives to reveal the mechanism of ortho-k and contribute to the evaluation of the lens performance and corneal response during the treatment.
Materials and Methods
Study Design
The finite element methods were used to analyze the ortho-k process. First, the corneal models were built based on the common corneal thicknesses and corneal curvatures. Second, the corresponding ortho-k lenses were selected in accordance with the manufacturer’s fitting instructions, and the lens parameters were extracted for simulation. Then, the interactions during the wearing process were numerically modeled, and corneal mechanical responses were investigated under these specific conditions. According to the post-fitting changes in the corneal topography, the central area was defined as the central 4 mm area of the corneal apex, outside of which was defined as the peripheral area. The central and peripheral corneal stress under the ortho-k lens was quantified, and the impact of a refractive target, corneal curvature, and corneal thickness on the corneal stress were analyzed.
Corneal Geometry
Based on the common corneal curvatures and corneal thicknesses, cornea models were developed, with the curvature of 40.25 D, 41.25 D, 42.25 D, 43.25 D, and 44.25 D, respectively, and thicknesses of 500 μm, 525 μm, 550 μm, and 575 μm, respectively, as shown in Table 1. The corneal radius was calculated by the following formula:
where
Ortho-k Lens Geometry
Four-zone ortho-k lenses (α ortho-k; Alpha Corp., Nagoya, Japan) made of highly oxygen permeable lens material (Boston EM, nominal Dk of 104 × 10–11 (cm2/sec) (mL O2/mL
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743745/fbioe-09-743745-HTML/image_m/fbioe-09-743745-g001.jpg)
FIGURE 1. (A) OCT B-Scan image of ortho-k lens from the Optimec is830. (B) Three-dimensional model of ortho-k lens.
Finite Element Model
Histologically, the microanatomical structure of the cornea can be divided into five layers from the anterior to the posterior: epithelium, Bowman’s membrane, stroma, Descemet’s membrane, and endothelium. In our modeling of the cornea, three layers of different elastic modulus were assigned to represent the real corneal structure. The corneal epithelium and Bowman’s membrane were regarded as one layer, as was Descemet’s membrane and endothelium. The layer thickness of the epithelial and Bowman’s membrane was presumed to be 50–57.5 μm, 10 percent of the total corneal thickness, as in routine clinical findings (Ehlers et al., 2010). The cornea was assumed to be a nonhomogeneous elastic material (Pinsky et al., 2005; Elsheikh et al., 2009; Last et al., 2009), and all material properties are shown in Table 3 (Wollensak et al., 2003; Pinsky et al., 2005; Elsheikh et al., 2008; Elsheikh et al., 2009; Last et al., 2009). According to the previous studies (Elsheikh et al., 2008; Elsheikh et al., 2009), the stiffness of both the outermost and innermost layers were assumed to be 10 percent as stiff as the corneal stroma. The finite element mesh of structured elements was generated. A mesh dependency study was conducted such that the relative error in two consecutive mesh refinements was <0.5% for computation. The mesh and boundary conditions are shown in Figure 2. Clinical studies have proven that principally only the thickness of the corneal epithelium and the curvature of the anterior cornea are changed when wearing an ortho-k lens, without bending of the entire cornea or affecting the depth of the anterior chamber (Alharbi and Swarbrick, 2003; Kim et al., 2018; Swarbrick et al., 2020). Therefore, in our model, the sclera was constrained from moving in all directions. To best match with the model of the whole eye, the orientation of corneal edge support was set to 23 (Elsheikh and Wang, 2007). The intraocular pressure (IOP) was set to 13 mmHg and homogeneously applied on the inner corneal surface, while the eyelid pressure was set to 9 mmHg and applied on the outer surface of the lens (Shaw et al., 2010). Eighty static finite element simulations were performed using the Abaqus software (version 6.14 Dassault Systemes, Simulia Corp., RI, the United States) for revealing the stress and displacement status of ortho-k lenses’ wearing.
Statistical Analysis
The normality of residuals and homogeneity of variance were tested by using the Shapiro–Wilk test and Levene’s test, respectively. Because data were found to violate normality, median quantile regression was conducted to explore factors influencing central corneal stress and peripheral corneal stress. Quantile regression is an approach that models the distribution of the responses and performs better than conventional means of regression when the normality assumption is violated. The central corneal stress value and the peripheral corneal stress value were, respectively, used as dependent variables (continuous variables), and the degree of myopia, corneal curvature, and corneal thickness were included as independent variables (categorical variables) in the models. Furthermore, to estimate the effect size of the independent variables on dependent variables, we used the analysis of variance (ANOVA) method and reported partial eta square as the measure of effect size. Partial eta squared (
Since residuals in the stress value model are non-normally distributed and residuals in the curvature change do not have constant variance, the Box–Cox transformation was performed on the stress value and the curvature change to meet the ANOVA assumptions. The Box–Cox transformation can obtain a normal distribution of the transformed data (after transformation) and a constant variance. The one-parameter Box–Cox transformation is defined as
and the two-parameter Box–Cox transformations as
Moreover, the first transformation holds for
Results
Corneal Biomechanical Responses Under ortho-k With FEA
In the current study, the central corneal stress and the peripheral corneal stress were regarded as continuous variables, while myopia, corneal curvature, and corneal thickness were regarded as categorical variables. The maximum von Mises stress in the central corneal area and the peripheral area are shown in Table 4 and Table 5. It can be found that the maximum von Mises stress of the corneal central area showed an upward trend with the increase of myopia reduction degree, corneal curvature, and corneal thickness. The average of central cornea stress was 0.019178 MPa (SD
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743745/fbioe-09-743745-HTML/image_m/fbioe-09-743745-g003.jpg)
FIGURE 3. Distribution of the von Mises stress in the cornea under varied conditions. (A) Sagittal face: −2.00 D/40.25 D/500 μm; (B) Sagittal face: −5.00 D/44.25 D/575 μm; (C) En face: −2.00 D/40.25 D/500 μm; (D) En face: −5.00 D/44.25 D/575 μm.
The trend of the maximum von Mises stress in the peripheral cornea was opposite to that in the corneal center. As the cornea thickness increased, the maximum von Mises stress in the peripheral area decreased gradually under the same curvature. The change of the maximum von Mises stress under various degrees was not large, and the maximum von Mises stress ranged from 0.00985 to 0.01733 MPa, with an average stress of the 80 models being 0.014451 MPa (SD = 0.011908). For −3.00 D, −4.00 D, and −5.00 D lens wearing, the maximum von Mises stress occurred in the central cornea, while for −2.00 D, the maximum von Mises stress was at the peripheral cornea. The representative maximum von Mises stress appearing in the central and peripheral areas are shown in Figure 4 and Figure 5.
Quantile Regression Results of Dependent Variables and Independent Variables
The quantile regression results of central corneal stress and independent variables are shown in Table 6. Compared with myopia of −2.00 D, the central corneal stress value increased when the myopia was −3.00 D (β = 0.0036, p < 0.001), −4.00 D (β = 0.0074, p < 0.01), and −5.00 D (β = 0.0111, p < 0.01). Compared with corneal curvature of 40.25 D, the central corneal stress value increased when the corneal curvature was 41.25 D (β = 0.0007, p < 0.01), 42.25 D (β = 0.0014, p < 0.01), 43.25 D (β = 0.0019, p < 0.01), and 44.25 D (β = 0.0031, p < 0.01). Compared with a corneal thickness of 500 μm, the central corneal stress value increased when the corneal thickness was 525 μm (β = 0.0004, p = 0.0298), 550 μm (β = 0.0007, p < 0.01), and 575 μm (β = 0.0011, p < 0.01). All models were controlled for the other covariates.
The quantile regression results of corneal peripheral stress and independent variables are shown in Table 7. In patients with myopia of −3.00 D (β = 0.0007, p < 0.01) and −4.00 D (β = 0.0009, p < 0.01), the corneal peripheral stress increased compared with patients with myopia of −2.00 D. In patients with corneal curvature of 42.25 D (β = 0.0007, p = 0.0324), 43.25 D (β = 0.0015, p < 0.01), and 44.25 D (β = 0.0040, p < 0.01), the peripheral corneal stress value increased compared with patients with corneal curvature of 40.25 D. In patients with corneal thickness of 525 μm (β = −0.0005, p < 0.01), 550 μm (β = −0.0009, p < 0.01), and 575 μm (β = −0.0013, p < 0.01), the peripheral corneal stress value decreased compared with patients with corneal thickness of 500 μm. All models were controlled for the other covariates.
Comparison of the Effect of Each Independent Variable on the Dependent Variable
The larger the partial eta value, the greater the effect of the independent variable on the dependent variable. The results of the variance (ANOVA) method are as follows: The degree of myopia had the greatest effect on the central corneal stress value (partial eta square
Discussion
The rapidly growing myopic prevalence has sparked intense research interest in the methods of myopia control (Cho and Cheung, 2012; Gonzalez-Meijome et al., 2016; Bremond-Gignac, 2020; Bullimore and Johnson, 2020). Ortho-k, a most promising treatment approach to slow down myopia progression, is gaining popularity among myopic patients due to its effectiveness and convenience. By using specially designed gas-permeable contact lenses to reshape the cornea, ortho-k has an effect in temporarily correcting myopia and controlling myopia progression (Cho et al., 2008; Hiraoka et al., 2012). However, the corneal integrity can be compromised if an improperly designed lens is fitted (Li S.-M. et al., 2016; Cheng et al., 2016; Bullimore et al., 2021). Therefore, there is an urgent need to understand the biomechanical change during ortho-k treatment and to identify the population most suitable for this treatment.
Currently all the existing theoretical hypotheses that the ortho-k lens controls myopia progression, such as the retinal peripheral defocus theory (Gardner et al., 2015; Smith et al., 2020), the aberration theory (Hiraoka et al., 2015; Zhang X. et al., 2020; Lau et al., 2020), and the accommodation theory (Batres et al., 2020; Song et al., 2021), are based on the change of the corneal shape or reduction of the refractive error after ortho-k lens wear. The biomechanical interaction between the lens and the cornea or the entire eyeball is vitally important in the course of ortho-k treatment. Different from the regular rigid contact lenses, the ortho-k lens has a special design, which is reversed against the corneal surface to reshape the cornea. The difference in the lens parameters will create different pressure to the cornea. Under the action of this force, the corneal morphology changes. Then, the lens can achieve its effects of reducing refractive error and controlling the myopia progression. However, the force of the interaction between the lens and cornea, or the pressure under the ortho-k lens on the cornea is difficult to measure or calculate. To the best of the authors’ knowledge, the research in this field remains very limited. Therefore, mathematical modeling quantifying the pressure and force that the cornea receives can play a key role to reveal the mechanism.
A finite element analysis (FEA), a commonly used numerical method, is an effective way to study and analyze the mechanical behavior and functional mechanism of the cornea, which has been extensively implemented to study corneal mechanical problems (Kling et al., 2014; Li ZD. et al., 2016; Fang et al., 2020; Meng et al., 2020). Not only has the FEA been widely recognized for its potential in simulating and analyzing the effects of corneal surgery, it has also been used in corneal diseases and ocular trauma because of its predictive feature. The focus of our current work was to use the FEA to analyze the biomechanical response of the cornea under the ortho-k lens. Based on the finite element models, this article compared and analyzed the von Mises stress in different areas of the corneas of various thicknesses and curvatures within the wearing course of a variety of ortho-k lens designs.
The whole-eye model with both the cornea and the sclera would provide a better representation of the actual state. However, it was extremely time-consuming. Some studies showed that if the boundary conditions were defined properly along the corneal edge, the cornea-only models can represent the whole-eye model very well. We set the orientation of corneal edge support to 23 in FEA according to Elsheikh and Wang (2007). The final results of our study were fair and showed consistent trends with clinical studies. In our model, the selection ranges of corneal curvature and corneal thickness were based on common clinical values. Although dividing the cornea into five layers according to the physiological structure can better show its mechanical response, due to the lack of relevant research data, it is impossible to accurately define the material properties of each layer and the interaction between the layers. Many articles showed that the corneal changes after wearing the ortho-k lens were mainly confined to the epithelial layer (Kim et al., 2018; Mohidin et al., 2020). Furthermore, the other layers of the cornea are relatively much thinner than the stromal layer, which may cause little effect on the simulation results. Therefore, the hypothesis of the three-layer corneal model is reasonable. Recently, many studies performed the corneal simulation using the hyper-elastic material model (Elsheikh and Wang, 2007; Bao et al., 2018). However, the behavior of a nonlinear material model is very similar to that of the elastic material model under 13 mmHg (1.73 KPa) IOP (Elsheikh and Wang, 2007). Nonlinear elastic materials can deform linearly for very small deformation (Wollensak et al., 2003). In order to improve the simulation’s efficiency, the corneal material was assumed to be elastic in this study.
According to the above results, the change of the targeted myopia reduction degree had the greatest impact on corneal stress. In the ortho-k lens design, the reduction degree mainly affects the base curve radius (BCR), a means of controlling the refractive change, which makes the cornea flatter if the target reduction is increased and base curve becomes flatter. From a biomechanical point of view, increasing the target degree means applying more compressive pressure on the corneal center. We analyzed the corneal responses with the change of targeted reduction degree. As the refractive target increased, both the central and peripheral corneal stress responses increased significantly, especially in the corneal center. Although the FDA or other health authorities approve the use of ortho-k treatment in myopia up to −6.00 D in some lens designs, the safest and most successful outcome will be achieved for those patients with lower baseline levels of myopia (Mika et al., 2007; Liu and Xie, 2016; Singh et al., 2020; Hu et al., 2021). Therefore, for patients with high myopia or for doctors who used to increase the target degrees, more attention must be paid to the possible side effects of the ortho-k lens wear, such as damage to the integrity of the cornea, since those lenses are designed with increased compressive pressure.
The corneal curvature is a very important parameter in ortho-k lens design. The calculation of BCR and the design of the lens alignment curve are based on the corneal curvature. The range of the corneal curvatures for most patients is between 40 and 46 D, and the most common is between 42 and 43 D. If the cornea is too flat or too steep, it can cause problems during lens fitting. When the cornea is too flat, the range of corneal deformation is limited, and adequate treatment may not be achieved, especially for high myopia. However, if the corneal is too steep, the incident of corneal staining is higher than that of the common corneal curvature (Cheng et al., 2016). In our study, it was found that the cornea had a larger stress response when the corneal curvature became steeper for both the central and peripheral areas. Therefore, practitioners need to be more cautious in treating patients with steep corneas and high myopia for safety concerns.
From our research results, it can be seen that corneal thickness also affects corneal stress in the central and peripheral areas. However, in the normal range, the corneal thickness is the factor that has the least influence on the corneal stress response among all the factors. After wearing ortho-k lenses, the corneal thickness has a certain amount of thinning (Wan et al., 2021). However, a large number of studies have found that the changes of corneal thickness are mainly confined to the corneal epithelium (Kim et al., 2018; Zhang J. et al., 2020; Zhou et al., 2020a; b). Therefore, it is relatively safe to wear an ortho-k lens within the range of normal corneal thickness.
There are some limitations in the current study. First, our model was mainly designed to analyze the positive pressure of the lens on the cornea without analyzing the tear layer force. Factors such as squeeze film force, or tension force, which are related to the tear viscosity and the tear thickness under the lens, are important in helping to alter the corneal shape and ideally centralize the lens position. In the future research, we will establish the tear flow field and analyze the tear layer under the lens. Second, this research assumed the cornea to be an ideal sphere. However, there have been some controversy about the changes of corneal shape during ortho-k, and the initial corneal eccentricity may affect lens selection. Creating the definite model of the eyeball according to the actual shape of the human cornea could be the foundation for the individual predictive model. Third, in this study, we were trying to reveal the biomechanical response of the cornea in the ortho-k treatment under the condition of various corneal shapes and thicknesses, and a refractive change with the uniform setting criteria of the pressures, including the IOP and eyelid tension. However, the lid tension and IOP may have potential effects on ortho-k. Therefore, we will explore the influence of these factors on the corneal mechanical response during ortho-k treatment in the future.
Conclusion
This article investigated the interaction between the cornea and the ortho-k lens through finite element methods to clearly reveal the mechanical response of the cornea under ortho-k. The current work provides ortho-k lens practitioners with an objective, reliable, and noncontact approach for quantitative assessment of risk and effectiveness of ortho-k for specific patients. This research also provides insights into the design of future ortho-k lenses. In ortho-k treatment, both the refractive target and the corneal curvature have a significant impact on corneal stress. Any change in corneal thickness within the normal range has minimal impact on stress. In clinical practice, more attention should be paid to patients with high myopia and high corneal curvature.
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 authors.
Author Contributions
JW, DZ, and QR contributed to the conception of the study. JW performed lens measurement and wrote the manuscript. WF, DZ, JW, and XL contributed to the finite element analysis. JW and HX performed data analyses. QR and DZ reviewed the manuscript, and they contributed equally to this study. All authors approved it for publication.
Funding
This work was supported by the National Key Research and Development Program (2016YFB1101503).
Conflict of Interest
XL was employed by the company the Second Academy Of CASIC.
The remaining 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.
Acknowledgments
We thank all participants of the study at the Peking University.
References
Alharbi, A., and Swarbrick, H. A. (2003). The Effects of Overnight Orthokeratology Lens Wear on Corneal Thickness. Invest. Ophthalmol. Vis. Sci. 44 (6), 2518–2523. doi:10.1167/iovs.02-0680
Bao, F., Wang, J., Cao, S., Liao, N., Shu, B., Zhao, Y., et al. (2018). Development and Clinical Verification of Numerical Simulation for Laser In Situ Keratomileusis. J. Mech. Behav. Biomed. Mater. 83, 126–134. doi:10.1016/j.jmbbm.2018.04.016
Batres, L., Peruzzo, S., Serramito, M., and Carracedo, G. (2020). Accommodation Response and Spherical Aberration during Orthokeratology. Graefes Arch. Clin. Exp. Ophthalmol. 258 (1), 117–127. doi:10.1007/s00417-019-04504-x
Bremond-Gignac, D. (2020). Myopie de l'enfant. Med. Sci. (Paris) 36 (8-9), 763–768. doi:10.1051/medsci/2020131
Bullimore, M. A., and Johnson, L. A. (2020). Overnight Orthokeratology. Contact Lens and Anterior Eye 43 (4), 322–332. doi:10.1016/j.clae.2020.03.018
Bullimore, M. A., Mirsayafov, D. S., Khurai, A. R., Kononov, L. B., Asatrian, S. P., Shmakov, A. N., et al. (2021). Pediatric Microbial Keratitis with Overnight Orthokeratology in Russia. Eye Contact Lens 47, 420–425. doi:10.1097/ICL.0000000000000801
Cheng, H.-C., Liang, J.-B., Lin, W.-P., and Wu, R. (2016). Effectiveness and Safety of Overnight Orthokeratology with Boston XO2 High-Permeability Lens Material: A 24 Week Follow-Up Study. Contact Lens and Anterior Eye 39 (1), 67–71. doi:10.1016/j.clae.2015.07.002
Cheung, S.-W., and Cho, P. (2016). Long-term Effect of Orthokeratology on the Anterior Segment Length. Contact Lens and Anterior Eye 39 (4), 262–265. doi:10.1016/j.clae.2016.02.003
Cho, P., and Cheung, S.-W. (2017). Protective Role of Orthokeratology in Reducing Risk of Rapid Axial Elongation: A Reanalysis of Data from the ROMIO and TO-SEE Studies. Invest. Ophthalmol. Vis. Sci. 58 (3), 1411–1416. doi:10.1167/iovs.16-20594
Cho, P., and Cheung, S.-W. (2012). Retardation of Myopia in Orthokeratology (ROMIO) Study: a 2-year Randomized Clinical Trial. Invest. Ophthalmol. Vis. Sci. 53 (11), 7077–7085. doi:10.1167/iovs.12-10565
Cho, P., Cheung, S. W., and Edwards, M. (2005). The Longitudinal Orthokeratology Research in Children (LORIC) in Hong Kong: a Pilot Study on Refractive Changes and Myopic Control. Curr. Eye Res. 30 (1), 71–80. doi:10.1080/02713680590907256
Cho, P., Cheung, S. W., Mountford, J., and White, P. (2008). Good Clinical Practice in Orthokeratology. Contact Lens and Anterior Eye 31 (1), 17–28. doi:10.1016/j.clae.2007.07.003
Cho, P., and Tan, Q. (2019). Myopia and Orthokeratology for Myopia Control. Clin. Exp. Optom. 102 (4), 364–377. doi:10.1111/cxo.12839
Coldrick, B. J., Richards, C., Sugden, K., Wolffsohn, J. S., and Drew, T. E. (2016). Developments in Contact Lens Measurement: A Comparative Study of Industry Standard Geometric Inspection and Optical Coherence Tomography. Contact Lens and Anterior Eye 39 (4), 270–276. doi:10.1016/j.clae.2016.01.002
Ehlers, N., Heegaard, S., Hjortdal, J., Ivarsen, A., Nielsen, K., and Prause, J. U. (2010). Morphological Evaluation of normal Human Corneal Epithelium. Acta Ophthalmol. 88 (8), 858–861. doi:10.1111/j.1755-3768.2009.01610.x
Elsheikh, A., Alhasso, D., and Rama, P. (2008). Assessment of the Epithelium's Contribution to Corneal Biomechanics. Exp. Eye Res. 86 (2), 445–451. doi:10.1016/j.exer.2007.12.002
Elsheikh, A., Ross, S., Alhasso, D., and Rama, P. (2009). Numerical Study of the Effect of Corneal Layered Structure on Ocular Biomechanics. Curr. Eye Res. 34 (1), 26–35. doi:10.1080/02713680802535263
Elsheikh, A., and Wang, D. (2007). Numerical Modelling of Corneal Biomechanical Behaviour. Comp. Methods Biomech. Biomed. Eng. 10 (2), 85–95. doi:10.1080/10255840600976013
Fang, L., Wang, Y., Yang, R., Deng, S., Deng, J., and Wan, L. (2020). Effects of the LASIK Flap Thickness on Corneal Biomechanical Behavior: a Finite Element Analysis. BMC Ophthalmol. 20 (1), 67. doi:10.1186/s12886-020-01338-8
Gardner, D. J., Walline, J. J., and Mutti, D. O. (2015). Choroidal Thickness and Peripheral Myopic Defocus during Orthokeratology. Optom. Vis. Sci. 92 (5), 579–588. doi:10.1097/OPX.0000000000000573
González-Méijome, J. M., Carracedo, G., Lopes-Ferreira, D., Faria-Ribeiro, M. A., Peixoto-de-Matos, S. C., and Queirós, A. (2016). Stabilization in Early Adult-Onset Myopia with Corneal Refractive Therapy. Contact Lens and Anterior Eye 39 (1), 72–77. doi:10.1016/j.clae.2015.06.009
Hiraoka, T., Kakita, T., Okamoto, F., and Oshika, T. (2015). Influence of Ocular Wavefront Aberrations on Axial Length Elongation in Myopic Children Treated with Overnight Orthokeratology. Ophthalmology 122 (1), 93–100. doi:10.1016/j.ophtha.2014.07.042
Hiraoka, T., Kakita, T., Okamoto, F., Takahashi, H., and Oshika, T. (2012). Long-term Effect of Overnight Orthokeratology on Axial Length Elongation in Childhood Myopia: a 5-year Follow-Up Study. Invest. Ophthalmol. Vis. Sci. 53 (7), 3913–3919. doi:10.1167/iovs.11-8453
Hu, P., Zhao, Y., Chen, D., and Ni, H. (2021). The Safety of Orthokeratology in Myopic Children and Analysis of Related Factors. Contact Lens and Anterior Eye 44 (1), 89–93. doi:10.1016/j.clae.2020.08.011
Jonas, J. B., Ang, M., Cho, P., Guggenheim, J. A., He, M. G., Jong, M., et al. (2021). IMI Prevention of Myopia and its Progression. Invest. Ophthalmol. Vis. Sci. 62 (5), 6. doi:10.1167/iovs.62.5.6
Kim, W. K., Kim, B. J., Ryu, I.-H., Kim, J. K., and Kim, S. W. (2018). Corneal Epithelial and Stromal Thickness Changes in Myopic Orthokeratology and Their Relationship with Refractive Change. PLoS One 13 (9), e0203652. doi:10.1371/journal.pone.0203652
Kling, S., Bekesi, N., Dorronsoro, C., Pascual, D., and Marcos, S. (2014). Corneal Viscoelastic Properties from Finite-Element Analysis of In Vivo Air-Puff Deformation. PLoS One 9 (8), e104904. doi:10.1371/journal.pone.0104904
Kobayashi, Y., Yanai, R., Chikamoto, N., Chikama, T.-i., Ueda, K., and Nishida, T. (2008). Reversibility of Effects of Orthokeratology on Visual Acuity, Refractive Error, Corneal Topography, and Contrast Sensitivity. Eye Contact Lens 34 (4), 224–228. doi:10.1097/ICL.0b013e318165d501
Last, J. A., Liliensiek, S. J., Nealey, P. F., and Murphy, C. J. (2009). Determining the Mechanical Properties of Human Corneal Basement Membranes with Atomic Force Microscopy. J. Struct. Biol. 167 (1), 19–24. doi:10.1016/j.jsb.2009.03.012
Lau, J. K., Vincent, S. J., Cheung, S.-W., and Cho, P. (2020). Higher-order Aberrations and Axial Elongation in Myopic Children Treated with Orthokeratology. Invest. Ophthalmol. Vis. Sci. 61 (2), 22. doi:10.1167/iovs.61.2.22
Lau, J. K., Wan, K., Cheung, S.-W., Vincent, S. J., and Cho, P. (2019). Weekly Changes in Axial Length and Choroidal Thickness in Children during and Following Orthokeratology Treatment with Different Compression Factors. Trans. Vis. Sci. Tech. 8 (4), 9. doi:10.1167/tvst.8.4.9
Li, S.-M., Kang, M.-T., Wu, S.-S., Liu, L.-R., Li, H., Chen, Z., et al. (2016a). Efficacy, Safety and Acceptability of Orthokeratology on Slowing Axial Elongation in Myopic Children by Meta-Analysis. Curr. Eye Res. 41 (5), 600–608. doi:10.3109/02713683.2015.1050743
Li, Z. D., Bao, F. J., Wang, Q. M., and Huang, J. H. (2016b). [Finite Element Analysis of Astigmatic Keratotomy Based on Corneal Biomechanical Properties]. Zhonghua Yan Ke Za Zhi 52 (9), 674–680. doi:10.3760/cma.j.issn.0412-4081.2016.09.008
Liu, Y. M., and Xie, P. (2016). The Safety of Orthokeratology-A Systematic Review. Eye Contact Lens 42 (1), 35–42. doi:10.1097/ICL.0000000000000219
Meng, Q., Wang, X., Chen, W., Li, X., and He, R. (2020). [Finite Element Analysis of Determining Corneal Biomechanical Properties In Vivo Based on Corvis ST]. Sheng Wu Yi Xue Gong Cheng Xue Za Zhi 37 (4), 608–613. doi:10.7507/1001-5515.201911081
Mika, R., Morgan, B., Cron, M., Lotoczky, J., and Pole, J. (2007). Safety and Efficacy of Overnight Orthokeratology in Myopic Children. Optom. - J. Am. Optometric Assoc. 78 (5), 225–231. doi:10.1016/j.optm.2006.12.013
Mohidin, N., Mat Yacob, A., and Norazman, F. N. (2020). Corneal Thickness and Morphology after Orthokeratology of Six-Month Lens Wear Among Young Malay Adults. Med. J. Malaysia 75 (5), 538–542.
Pinsky, P. M., van der Heide, D., and Chernyak, D. (2005). Computational Modeling of Mechanical Anisotropy in the Cornea and Sclera. J. Cataract Refract Surg. 31 (1), 136–145. doi:10.1016/j.jcrs.2004.10.048
Queirós, A., Lopes‐ferreira, D., Yeoh, B., Issacs, S., Amorim‐de‐sousa, A., Villa‐collar, C., et al. (2020). Refractive, Biometric and Corneal Topographic Parameter Changes during 12-months of Orthokeratology. Clin. Exp. Optom. 103 (4), 454–462. doi:10.1111/cxo.12976
Shaw, A. J., Collins, M. J., Davis, B. A., and Carney, L. G. (2010). Eyelid Pressure and Contact with the Ocular Surface. Invest. Ophthalmol. Vis. Sci. 51 (4), 1911–1917. doi:10.1167/iovs.09-4090
Singh, K., Bhattacharyya, M., Goel, A., Arora, R., Gotmare, N., and Aggarwal, H. (2020). Orthokeratology in Moderate Myopia: A Study of Predictability and Safety. Jovr 15 (2), 210–217. doi:10.18502/jovr.v15i2.6739
Smith, E. L., Arumugam, B., Hung, L.-F., She, Z., Beach, K., and Sankaridurg, P. (2020). Eccentricity-dependent Effects of Simultaneous Competing Defocus on Emmetropization in Infant Rhesus Monkeys. Vis. Res. 177, 32–40. doi:10.1016/j.visres.2020.08.003
Song, Y., Zhu, S., Yang, B., Wang, X., Ma, W., Dong, G., et al. (2021). Accommodation and Binocular Vision Changes after Wearing Orthokeratology Lens in 8- to 14-Year-Old Myopic Children. Graefes Arch. Clin. Exp. Ophthalmol. 259, 2035–2045. doi:10.1007/s00417-021-05106-2
Swarbrick, H. A., Kang, P., and Peguda, R. (2020). Corneal Total and Epithelial Thickness Measured by Sonogage Ultrasound Pachometry and High-Resolution Optical Coherence Tomography. Optom. Vis. Sci. 97 (5), 346–350. doi:10.1097/OPX.0000000000001508
Swarbrick, H. A. (2004). Orthokeratology (Corneal Refractive Therapy). Eye Contact Lens: Sci. Clin. Pract. 30 (4), 181–185. discussion 205-186. doi:10.1097/01.icl.0000140221.41806.6e
Wan, K., Lau, J. K.-k., Cheung, S. W., and Cho, P. (2020). Refractive and Corneal Responses of Young Myopic Children to Short-Term Orthokeratology Treatment with Different Compression Factors. Contact Lens and Anterior Eye 43 (1), 65–72. doi:10.1016/j.clae.2019.10.134
Wan, K., Yau, H. T., Cheung, S. W., and Cho, P. (2021). Corneal Thickness Changes in Myopic Children during and after Short‐term Orthokeratology Lens Wear. Ophthalmic Physiol. Opt. 41, 757–767. doi:10.1111/opo.12824
Wolffsohn, J. S., Calossi, A., Cho, P., Gifford, K., Jones, L., Li, M., et al. (2016). Global Trends in Myopia Management Attitudes and Strategies in Clinical Practice. Contact Lens and Anterior Eye 39 (2), 106–116. doi:10.1016/j.clae.2016.02.005
Wollensak, G., Spoerl, E., and Seiler, T. (2003). Stress-strain Measurements of Human and Porcine Corneas after Riboflavin-Ultraviolet-A-Induced Cross-Linking. J. Cataract Refract Surg. 29 (9), 1780–1785. doi:10.1016/s0886-3350(03)00407-3
Zhang, J., Li, J., Li, X., Li, F., and Wang, T. (2020a). Redistribution of the Corneal Epithelium after Overnight Wear of Orthokeratology Contact Lenses for Myopia Reduction. Contact Lens and Anterior Eye 43 (3), 232–237. doi:10.1016/j.clae.2020.02.015
Zhang, X., Ma, J.-H., Xi, X., and Guan, L. (2020b). Characteristics of Corneal High-Order Aberrations in Adolescents with Mild to Moderate Myopia. BMC Ophthalmol. 20 (1), 465. doi:10.1186/s12886-020-01727-z
Zhou, J., Xue, F., Zhou, X., Naidu, R. K., and Qian, Y. (2020a). Correction to: Thickness Profiles of the Corneal Epithelium along the Steep and Flat Meridians of Astigmatic Corneas after Orthokeratology. BMC Ophthalmol. 20 (1), 338. doi:10.1186/s12886-020-01607-6
Keywords: orthokeratology, myopia control, ocular biomechanics, finite element simulation, corneal response
Citation: Wu J, Fang W, Xu H, Liu X, Zhao D and Rong Q (2021) The Biomechanical Response of the Cornea in Orthokeratology. Front. Bioeng. Biotechnol. 9:743745. doi: 10.3389/fbioe.2021.743745
Received: 19 July 2021; Accepted: 06 September 2021;
Published: 11 October 2021.
Edited by:
Weiyi Chen, Taiyuan University of Technology, ChinaReviewed by:
Ashkan Eliasy, University of Liverpool, United KingdomXingtao Zhou, Fudan University, China
Copyright © 2021 Wu, Fang, Xu, Liu, Zhao and Rong. 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: Qiguo Rong, qrong@pku.edu.cn; Dongliang Zhao, zhaodl@pku.edu.cn