Abstract

Purpose. To determine three-year change of the corneal biomechanical parameter stress-strain index (SSI) in schoolchildren aged 7– 9 years and their correlation with refractive error and axial length (AL). Methods. This is a prospective cohort study. Data of the AL, refractive error, and corneal biomechanical parameter SSI were collected at baseline and a 3-year follow-up for 217 schoolchildren. SSI, AL, and refractive error were measured via corneal visualization Scheimpflug technology (Corvis ST), IOLMaster biometry, and cycloplegic refraction. Three years of changes in SSI and its association with refractive error and AL were analyzed. Participants were divided into persistent nonmyopia (PNM), newly developed myopia (NDM), and persistent myopia (PM). The three-year difference in SSI among the three groups was analyzed. Results. After three years of follow-up, the corneal biomechanical parameter SSI decreased in all participants (). There was a negative correlation between the change in SSI and the change in AL (r = −0.205, ) and a positive correlation between the change in refractive error (r = 0.183, ). After three years of follow-up, there was a decrease in the SSI for the NDM, PM, and PNM participants, with a median change of −0.05 for PNM and −0.13 and −0.09 for the NDM and PM, respectively. There was a significant decrease in corneal biomechanical properties for NDM patients compared with PNM patients (). Conclusion. In 7- to 9-year-old schoolchildren, SSI decreased after three years of the longitudinal study, and the change in SSI was correlated with the change in AL and refractive error. There was a rapid decrease in corneal biomechanical properties among newly developed myopic patients.

1. Introduction

Myopia is a global public health problem and its progression can lead to visual impairment, such as myopic macular degeneration and retinal detachment. Previous studies have shown that progressive myopia reduced the elasticity, stress, and load-bearing capacity of the sclera [1, 2]. Human scleral biomechanics is measured mainly in vitro, which can only reflect the state of myopia after stabilization but not during development [1]. Nguyen et al. studied human and donor eyes and found that both corneal and scleral material properties contribute to the observed corneal response to air-puff-induced deformation [3]. It has been hypothesized that the biomechanical properties of the cornea may reflect those of the sclera since the corneal and the sclera tissue of the eye are structurally connected [4].

Prior studies have demonstrated that the biomechanical characteristics of the cornea were influenced by both central corneal thickness (CCT) and intraocular pressure (IOP) [57]. It is therefore difficult to obtain accurate corneal biomechanical measurements that are unaffected by CCT and IOP [8]. Eliasy et al. developed a novel biomechanical parameter, known as the material stiffness parameter stress-strain index (SSI), for the estimation of the in vivo elastic properties of the cornea [8]. There exists a significant correlation between the SSI value and the stiffness of the cornea. Specifically, when the SSI value increases, the corneal stiffness also tends to increase, and conversely, when the SSI value decreases, the corneal stiffness tends to decrease [8]. SSI is also a marker for material stiffness and corneal biomechanical parameters, which do not correlate with the CCT or IOP [8].

Recent cross-sectional studies found that eyes with higher myopia have lower corneal stiffness [9, 10] and may be less stable and more susceptible to stress [1113]. Similarly, in a cross-sectional study conducted by Liu et al., it was observed that the SSI value exhibited a sequential decrease when categorized according to the degree of myopia, namely, low, moderate, high, and ultra-high [14]. The difference in the rate of change of spherical equivalent (SE) refractive error may be related to the change in corneal stiffness and needs to be investigated in a longitudinal study. Results have been mixed from a limited number of longitudinal studies in this area. Shah et al. found no significant association between changes in corneal biomechanical properties and axial length (AL) or SE over two years [15]. In contrast, Momeni-Moghaddam et al. reported a decrease in most of the corneal biomechanical parameters after a four-year follow-up [16]. It is noteworthy that they utilized the ocular response analyzer, which could not differentiate corneal behavior from the effects of IOP and CCT [17]. A longitudinal study is warranted to investigate the causality relationship between corneal biomechanics and axial length using corneal biomechanical parameters that are independent of IOP and CCT.

Therefore, the current study that included participants aged 7–9 years was designed to explore the changes in SSI in a three-year longitudinal cohort and investigate the potential role of the SSI in the development of myopia.

2. Methods

2.1. Study Design

Data for this study were collected as part of an ongoing prospective study that was approved by the institutional review board of Sun Yat-sen University School of Public Health. The study was conducted in accordance with the Declaration of Helsinki. Parents were informed about the study’s purpose and content before obtaining their written consent.

In November and December 2018, 12 primary schools in Guangzhou City, China, randomly selected Grade 2 and Grade 3 students to participate in the study. After three years of follow-up, the study was completed in November and December 2022.

This study excluded participants with strabismus, amblyopia, and nystagmus. Participants who used orthokeratology were also excluded. This study measured AL using a laser interferometer (IOLMaster 500; Zeiss). Cycloplegic autorefraction was performed by an autorefractor (KR 8800; Topcon Corp.). A cyclopentolate solution of 1% was administered three times over five minutes to induce cycloplegia. Following complete cycloplegia, measurements of the refractive index were taken. We also performed cycloplegic subjective refraction on children with vision worse than 20/20 in either eye.

A child with persistent nonmyopia (PNM) was defined as having a SE value greater than −0.50D at both visits. Newly developed myopia (NDM) was defined as a child who was nonmyopic (SE > −0.50D) at baseline but had myopia (SE ≤ −0.50D) at the 3-year follow-up visit. A child who was myopic (SE ≤ −0.50D) at both visits was defined as having persistent myopia (PM) [18].

2.2. Corvis ST and SSI

A visual dynamic intraocular pressure analyzer, Corvis ST (software version 1.6b2507), was used to measure corneal biomechanical properties. A quality score of “OK” was used for the analysis of all measurements. A repeat examination was conducted if the quality was not rated “OK.” Subjects whose multiple attempts did not meet an “OK” assessment were excluded.

Software used least squares analysis and numerical modeling to estimate SSI. Parameters such as CCT, bIOP, and SP-HC were used as inputs [8].

The equation of SSI reported in the Eliasy et al.’s study was SSI = f(a1 + a2C1 + a3C2 + a4C + a5C1C'2 + a6C+ a7C + a8CC2 + a9C1C + C) + ln (SP-HC)) [8]. The C1 and C2 are obtained by dividing the values of CCT and bIOP by 545 and 20, respectively. The natural logarithm of the difference between SP and HC, denoted as Ln (SP-HC), is calculated. The constants a1 to a9 are determined through a fitting process, where the equation is adjusted to match the given numerical input and output values [8].

2.3. Data Analysis

Statistics analysis was performed by IBM SPSS (IBM SPSS Statistics version 25.0, IBM SPSS, USA). The Shapiro–Wilk test was performed to check data normality. Right eyes were the only ones included for Spearman analysis since the right and left eyes showed high correlations (r, 0.91 for SE and 0.88 for AL). We analyzed baseline and follow-up differences using Kruskal–Wallis H tests. The Wilcoxon test was used to compare the difference between baseline and follow-up. Spearman correlation was used to analyze the correlation among SSI, refraction error, and AL. A value lower than 0.01 was considered statistically significant.

3. Results

3.1. General Characteristics

217 schoolchildren with a median age of 8.1 years at baseline, of which 52.07% were girls, were available for analysis. Boys had a greater AL than girls at baseline and the follow-up visit (). No significant differences were observed in SSI and SE between genders (Table 1).

After a three-year follow-up period, the SSI decreased significantly (), the eye elongated significantly, and myopia shifted significantly (Table 1). There was no significant difference between the 7.0–8.0 age groups and the 8.1–9.0 age groups in these ocular parameters (Table 2).

3.2. Relationship between Changes in the SSI and Other Ocular Components

There was a positive correlation between changes in the SSI and SE (r = 0.183, ) and a negative correlation between changes in the SSI and AL (r = −0.205, ) (Table 3).

3.3. SSI Changes over Time among Refractive Status Groups

The 217 participants were divided into 92 PNMs, 97 NDMs, and 28 PMs. At baseline, PMs showed lower SSI values than PNMs; NDMs and PNMs showed no significant differences (Table 4). After three years of follow-up, all participants decreased in SSI, with PNMs losing 0.05, NDM 0.13, and PMs 0.09; the decrease of SSI in NDM was greater than in PNM. SE and AL changed more significantly in the NDM and PM than in the PNM (Table 4).

4. Discussion

This study found that after three years of follow-up, the corneal biomechanical parameter stress-strain index decreased as the degree of myopia increased. The change in SSI correlated with the change in SE and AL. Moreover, there was no significant difference between the PNM and NDM groups in the baseline. After three years of follow-up, however, a significant difference was observed between the PNM and NDM groups, suggesting that an increase in axial length may lead to a decrease in corneal biomechanical properties. In addition, patients with newly developed myopia showed greater attenuation in corneal stiffness than those with persistent nonmyopia.

Our study found that SSI decreased after three years of follow-up. This decrease can be attributed to the significant impact of axial elongation. In addition, the change in SSI was negatively correlated with the change in axial length, which is typically associated with myopia development. Momeni-Moghaddam et al. also found a similar difference with the ocular response analyzer (ORA) [16]. However, a two-year change in AL/SE and CH/CRF was not significantly associated with Shah et al.’s study [15]. The difference between our study and Shahe et al.’s study may be because the shift of AL and SE in our study was more prominent (0.82 mm and −2.0 DS) than in theirs (0.22 mm and −0.7 DS).

One of the major findings from this three-year longitudinal study was that the change of SSI in newly developed myopia is larger than in persistent nonmyopia. The dramatic decrease of the SSI in the NDM indicates that the corneal biomechanical properties also underwent dramatic changes with myopia increase.

Many cross-sectional studies have found that the SSI value was lower in eyes with higher spherical error [9, 10, 19]. In our longitudinal study, SSI in the NDM group showed lower values than PNM after three years of follow-up, indicating that SSI changes are secondary to the increase in axial length. The findings suggest that as children progress from emmetropia to myopia, there is a discernible decrease in SSI. However, it remains to be explored whether SSI can reliably indicate children at risk of developing myopia. NDM-related factors, such as reduced outdoor activity and heightened educational pressure, as well as excessive smartphone use, have been identified as influencing myopia development [20]. Notably, greater outdoor activity appears to be associated with a reduced risk of myopia [20]. In addition, higher degrees of myopic refractive error and longer axial lengths correlate with increased screen time on smartphones and computers [21]. These lifestyle factors likely contribute to the emergence of NDM. While this study did not incorporate underlying lifestyle factors that differentiate between PNM and NDM, investigating the causal relationships between myopia, lifestyle, and corneal biomechanics promises to advance our understanding of myopia progression. Future research endeavors should further explore the intricate interplay of lifestyle factors, myopia progression, and corneal biomechanics. The changes in corneal collagenous may be physiological factors of corneal biomechanical decrease in newly developed myopia. Scleras are fibrous collagenous shells that are continuous with the cornea anteriorly, creating essentially closed chambers in the anterior, equatorial, and posterior eye. Scleral thinning occurs very rapidly in response to the onset of myopia development [2]. Thus, corneal biomechanics—inseparable from whole-eye biomechanics—may undergo a decrease in corneal biomechanical properties.

Age-associated biomechanical changes in the cornea have been reported in a few studies in vivo. In our study, there was no difference in the biomechanical parameters of SSI in relation to age, similar to Matalia et al.’s study, where corneal deformation parameters were unaffected by age in vivo [22]. However, glycation-induced cross-linking of collagen molecules has been shown in other studies to cause corneal stiffening as people age [23, 24]. According to Kamiya et al.’s study, corneal biomechanical parameters decline significantly with age [25]. Matalia et al.’s study also showed that corneal biomechanical parameters stiffer with increasing age [22, 26]. Discrepancies in patient ages may explain this discrepancy in the results between us. In Kamiya et al. and Elsheikh et al.’s studies, the participants were with a wider range, aged 19–89 years old and 50–95 years old, respectively. In our study, a narrower age (7–9 years) of children was included; the two age groups in our study may be too close to produce a statistical difference. Furthermore, a study revealed that SSI exhibited a notable level of stability before age 35, followed by an increase as individuals increased in age [27]. The alteration in corneal biomechanics observed with advancing age may be attributed to modifications in the molecular composition of the cornea, which show an increased trend of collagen glycation and collagen molecules in individuals aged ≥40 years [23, 24, 28].

There is a recognized limitation in our study. We only included participants of Grade 2 and Grade 3, with a limited age and refraction range. Further research should be conducted involving a wider range of ages and refractions to gain a better understanding of the topic.

In summary, by conducting a study with a 3-year follow-up, we observed the role of the corneal biomechanical parameter stress-strain index during myopia development in schoolchildren and found that stress-strain index decreases and is associated with an increase in myopia. Moreover, a rapid attenuation of SSI was observed in newly developed myopic patients, suggesting its potential significance in myopia onset. Future research should consider including different age groups to provide a more comprehensive understanding. In addition, investigating the intricate relationship among lifestyle factors, myopia progression, and corneal biomechanics warrants further exploration in our upcoming studies.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

This study was approved by the Ethics Committee of the School of Public Health, Sun Yat-sen University (No. 004).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (82070994).