ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 16 July 2025

Sec. Biomechanics

Volume 13 - 2025 | https://doi.org/10.3389/fbioe.2025.1598546

Development and validation of the VAE-NT index: a novel biomechanical parameter for distinguishing subclinical corneal abnormalities

  • 1. Department of Ophthalmology, Huadong Hospital, Fudan Universiry, Shanghai, China

  • 2. National Engineering Research Center of Ophthalmology and Optometry, Eye Hospital, Wenzhou Medical University, Wenzhou, China

  • 3. National Clinical Research Center for Ocular Diseases, Eye Hospital, Wenzhou Medical University, Wenzhou, China

Article metrics

View details

1,4k

Views

300

Downloads

Abstract

Purpose:

The aim of this study is to develop an index for distinguishing between very asymmetric ectasia with normal topography (VAE-NT) eyes and normal eyes, with good performance in validity, reliability, and predictive values.

Methods:

In the training dataset, this single-center retrospective study involved 102 healthy eyes and 97 VAE-NT eyes. After propensity score matching (PSM), data from 53 healthy eyes and 53 VAE-NT eyes, including demographic and Corvis ST examination results, were collected. The area under the receiver operating characteristic curve (AUC), sensitivity, specificity, intraclass correlation coefficient (ICC), and positive and negative likelihood ratios were calculated for the dynamic corneal response (DCR) parameters of Corvis ST. The Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) model was used to objectively and comprehensively evaluate the Corvis ST DCRs, and logistic regression was used to determine the optimal combination of parameters that can accurately separate VAE-NT from normal corneas. In the validation dataset, 44 VAE-NT eyes and 49 normal eyes were involved. The validity, reliability, and predictive value of the index were further assessed using the validation dataset. The VAE-NT index was compared with the tomographic and biomechanical index (TBI) in both the training and validation datasets.

Results:

In the training dataset, the optimal parameter combination forming the VAE-NT index included the following DCRs: SP A1, SP HC, A1 Time, DA Ratio Max (2 mm), DA Ratio Max (1 mm), Integrated Radius, and stress–strain index version 2 (SSI2). The receiver operating characteristic (ROC) curve analysis showed an AUC value of 0.971, with a cut-off value of 0.425, an accuracy of 95.283%, a specificity of 94.340%, and a sensitivity of 96.230%. In the validation dataset, the AUC value of the VAE-NT index was 0.980. The sensitivity and specificity of the VAE-NT index were 93.180% and 95.920%, respectively. The positive and negative likelihood ratios of the VAE-NT index were 22.830 and 0.071, respectively. The ICC of the VAE-NT index was 0.835, and the accuracy was 94.624%. The VAE-NT index outperformed TBI in both the training and validation datasets.

Conclusion:

The VAE-NT index was developed, exhibiting high sensitivity, specificity, and AUC, along with favorable likelihood ratios and repeatability, suggesting that the VAE-NT index is a robust and reliable tool for distinguishing VAE-NT eyes from normal eyes. Further validation in broader populations and over longer follow-up periods is needed to support clinical translation.

Introduction

Forme fruste keratoconus (FFKC), recently termed very asymmetric ectasia with normal topography (VAE-NT), is a clinically significant condition characterized by normal topography and slit-lamp examination in one eye, whereas the fellow eye shows signs of keratoconus. This atypical manifestation implies an incomplete state of the disease, in which the cornea protrudes, causing irregular astigmatism and vision impairment (Henriquez et al., 2020). Early diagnosis of FFKC is crucial as it enables patients to proactively address the condition and prevent its progression into fully developed keratoconus (KC), thereby mitigating the risk of further vision loss (Rabinowitz, 1998).

The biomechanical properties of KC play a pivotal role in understanding FFKC. KC is a corneal disorder marked by alterations in the normal collagen fibril network, frequently exhibiting asymmetry between the two eyes. It has been proposed that the progression of keratoconus is driven by a biomechanical cycle of decompensation, involving corneal thinning, increased mechanical strain, and stress redistribution, initiated by a localized reduction in corneal material properties (Ruberti et al., 2011; Roberts and Dupps, 2014). In FFKC, subtle biomechanical abnormalities are often present despite the absence of overt clinical findings. These biomechanical alterations may precede morphological changes, potentially leading to progressive corneal protrusion, irregular astigmatism, and visual impairment if left undetected. To improve diagnostic precision, the term VAE-NT has been proposed to replace the previously used designation of FFKC (Ambrósio et al., 2023; Ambrósio et al., 2017; Hwang et al., 2018). Early identification of these subtle biomechanical alterations is critical for the diagnosis of VAE-NT, enabling timely intervention to prevent progression to clinically manifest keratoconus and preserve visual function.

Corneal biomechanics have gained significant attention in recent decades, with their importance recognized in several applications, including the measurement of intraocular pressure, evaluation of ectasia risk following refractive surgeries, and assessment of the effectiveness of corneal cross-linking (CXL) treatment (Roberts and Dupps, 2014; Herber et al., 2021; Ramm et al., 2019; Girard et al., 2015; Piñero and Alcón, 2015). Several in vivo methods have been developed to assess corneal biomechanics, among which the Corvis ST (Oculus, Wetzlar, Germany) is widely used. The Corvis ST uses an ultra-fast Scheimpflug camera that captures 140 frames over 32.11 ms, allowing detailed analysis of corneal deformation in response to an air puff stimulus (Tian et al., 2014). Analysis of the resulting deformation yields several dynamic corneal response (DCR) parameters, which correlate with corneal stiffness (Xian et al., 2023; Miao et al., 2023). These metrics include the stiffness parameter at first applanation (SP-A1) (Ambrósio et al., 2017; Zhang et al., 2021a), the deflection amplitude DA, and the ratio between the deflection amplitudes at apex and 2 mm away from apex (DA Ratio Max 2 mm) (Lu et al., 2022). Additional parameters, such as Ambrósio’s relational thickness to the horizontal profile (ARTh) (Tian et al., 2021a), the Corvis biomechanical index (CBI) (Wang et al., 2017), the stress–strain index (SSI) (Zhang et al., 2021b), and the Chinese CBI (cCBI) (Zhang et al., 2024), have all demonstrated clinical utility in the diagnosis of keratoconus (Ren et al., 2021).

The evaluation criteria for diagnostic indicators are primarily based on metrics derived from receiver operating characteristic (ROC) analysis, including the area under the ROC curve (AUC), sensitivity, specificity, and positive and negative likelihood ratios (Mandrekar, 2010). Although an ideal diagnostic indicator should exhibit high performance across all metrics, in practice, some indicators may not simultaneously achieve high AUC, sensitivity, and specificity. Therefore, clinicians often rely on their clinical experience to interpret these metrics, which can further complicate the diagnostic decision-making process. In addition, the intraclass correlation coefficient (ICC) is commonly used to assess the repeatability and stability of diagnostic indicators (Muller R, 1994). Therefore, this study aims to identify diagnostic indicators that demonstrate superior performance in terms of AUC, sensitivity, specificity, and ICC.

Multi-criteria decision analysis (MCDA) is a structured methodology developed to support decision-making processes involving multiple and often conflicting evaluation criteria (Talukder et al., 2018). As described by Keeney, MCDA provides a logical and systematic framework for evaluating options based on multiple criteria (Keeney, 1982). In this study, we used the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS), which is a widely used MCDA method. Hwang and Yoon (1981) originally proposed that TOPSIS evaluates alternatives based on their geometric proximity to an ideal solution and distance from a negative ideal solution. Due to its computational simplicity and robustness, TOPSIS has become one of the most widely adopted quantitative techniques in multi-criteria decision-making. Its applicability spans various domains and relies on a well-established mathematical foundation. It has been applied for more than three decades (Hwang and Yoon, 1981; Jahanshahloo et al., 2006), with extensive validation and documentation in the scientific literature (Huang et al., 2016; Yoon and Hwang, 1995). In the TOPSIS framework, the optimal alternative is defined as the one with the shortest distance to the positive ideal solution and the greatest distance from the negative ideal solution.

In this study, the validity, reliability, and predictive performance of Corvis ST DCR parameters for identifying VAE-NT were comprehensively evaluated using the TOPSIS approach. Based on sample size requirements and TOPSIS rankings, the seven highest-performing DCR parameters were selected and integrated to construct a novel composite biomechanical index for differentiating VAE-NT from normal corneas.

Methods

The steps followed in this study are illustrated in Figure 1.

FIGURE 1

Participants

All study participants were recruited at the Eye Hospital of Wenzhou Medical University. This single-center retrospective study initially enrolled 102 healthy eyes and 97 VAE-NT eyes. Each patient underwent a comprehensive eye examination, incorporating tests using the Pentacam and Corvis ST (Oculus Optikgeräte GmbH). The study complied with the ethical principles outlined in the Declaration of Helsinki and was approved by the Ethics Committee of the Eye Hospital of Wenzhou Medical University (ethics approval number: 2022-198-K-154).

In detecting KC, the following criteria were considered: (a) an irregular cornea, determined by distorted keratometry mires and distortion of the retinoscopic, ophthalmoscopic red reflex, or a combination of these and (b) the presence of at least one of the following biomicroscopic signs: Vogt’s striae, Fleischer’s ring of greater than 2-mm arc, or corneal scarring consistent with keratoconus (Rabinowitz, 1998; Arbelaez et al., 2012).

The criteria for very asymmetric ectasia (VAE) refer to the diagnosis of ectasia in one eye, according to previously established definitions, with the fellow eye considered clinically normal based on unremarkable corneal topography; VAE-NT eyes are the fellow eyes of these patients that have normal topography and a keratoconus percentage index (KISA%) score lower than 60, along with a paracentral inferior–superior (I–S value) asymmetry value at 6 mm (3-mm radii) of less than 1.45.

On the other hand, the inclusion criteria for healthy individuals included providing a signed informed consent form, qualification as a candidate for refractive surgery with the absence of topographic distortions (LASIK or SMILE), and having a Corvis ST assessment in the database. The exclusion criteria encompassed any prior ocular surgery or illness, myopia exceeding 10.00 diopters (D), concurrent or prior glaucoma, and the use of hypotonic treatments.

Corvis ST examinations

All DCR parameters were measured using the same Corvis ST device, and baseline values were recorded. In addition to the biomechanically corrected intraocular pressure (bIOP) and the central corneal thickness (CCT), the Corvis ST provided detailed information on corneal response to an air pulse. To eliminate inter-rater variations, the same technician, blinded to the study design, performed all the examinations. Only results with “OK” in the QS window indicating good image quality were included in the statistical analyses. The Corvis ST provided values of DCRs, including SP A1, DA Ratio Max (2 mm), ARTh, CBI, cCBI, and SSI version 2 (SSI2) (Zhang et al., 2021b). All the information regarding the included Corvis ST parameters is listed in Supplementary Table S1.

TOPSIS model

We implemented the TOPSIS model to create the VAE-NT index. This work used the TOPSIS method to build an evaluation system (Ramón-Canul et al., 2021).

In this process, positive and negative ideal solutions could be developed using Equation 1 and Equation 2, respectively:

The Euclidean distance from the ith evaluated object to the positive ideal solution is presented in Equation 3:

For each evaluation object, its Euclidean distance to the negative ideal solution is presented in Equation 4 as follows:

The closeness of the ith evaluation object to the ideal solution is described in Equation 5:

After normalization, the TOPSIS composite score could be obtained, as shown in Equation 6:where denotes the proximity of each evaluation object to the optimal solution. 0 ≤ ≤ 1; the closer was to 1, the better the evaluation object was.

Statistical analyses

The Shapiro–Wilk test assessed the normality of continuous variables. Descriptive statistics, including mean ± standard deviation (SD), were used for value description. ROC analysis evaluated the diagnostic efficacy of Corvis ST parameters and the new index for VAE-NT diagnosis. Reliability was assessed using the ICC for averaged measurements, using a two-way random-effects model with “single rater” type and “absolute agreement” ICC definition (Koo and Li, 2016). Propensity score matching (PSM) was used to match CCT between groups, reducing confounding bias. Binary logistic regression with backward stepwise inclusion determined the optimal combination of predictors from individual Corvis ST parameters for creating the VAE-NT index. Due to the requirements for the sample size based on research design and statistical methods, the top seven parameters ranked based on TOPSIS models and AUC values were considered. MedCalc software version 12.3.0.0 (Ostend, Belgium) was used for ROC analysis in both the TOPSIS and AUC groups. R was used for PSM matching. Other statistical analyses were conducted using SPSS (version 22.0, IBM, Inc.). A significance level of 0.05 was applied.

Results

Demographic data

Initially, the study included 102 healthy eyes and 97 VAE-NT eyes. After PSM, 53 healthy eyes and 53 VAE-NT eyes were selected for further analysis in the training dataset. The remaining eyes were included in the validation dataset. Table 1a provides information on participants’ age, CCT, bIOP, and gender in the training dataset. The data for CCT in both eyes were comparable (t = −0.019; P = 0.850).

TABLE 1

ParameterVAE-NT(53 eyes)Normal (53 eyes)Test statistic Zχ2t valueP
Age18.736 ± 5.36820.264 ± 3.187−2.0770.038
GenderF/M = 11/42F/M = 16/371.2420.265
CCT548.148 ± 27.664549.153 ± 26.904−0.1900.850
bIOP13.947 ± 2.07616.124 ± 1.672−5.947<0.001
ParameterVAE-NT (44 eyes)Normal (49 eyes)Test statistic Zχ2t valueP
Age20.955 ± 4.77920.837 ± 3.436−0.1050.917
GenderF/M = 15/29F/M = 9/402.9930.084
CCT509.898 ± 16.074572.553 ± 21.41815.811<0.001
bIOP14.376 ± 1.86015.894 ± 2.806−3.1170.002

(a) Demographic data of the training dataset. (b) Demographic data of the validation dataset.

Notes: CCT, central cornea thickness; bIOP, biomechanically corrected intraocular pressure.

Table 1b provides information on participants’ age, CCT, bIOP, and gender in the validation dataset. The data for CCT in both eyes were significantly different (t = 15.811; P < 0.001).

Assessment criteria included in the TOPSIS model

The Corvis DCR parameters were assessed in terms of validity, reliability, and predictive values. Validity included AUC, sensitivity, and specificity, while reliability was measured using ICC. The predictive value was evaluated through + LR (positive likelihood ratio) and −LR (negative likelihood ratio). Detailed criteria for evaluation are presented in Table 2.

TABLE 2

Assessment
aspect
Assessment indicatorDescription of the indicator
ValidityAUCArea under the receiver operating characteristic curve
SensitivityProportion of true positive tests out of all patients with a condition
SpecificityPercentage of true negatives out of all subjects who do not have a disease or condition
ReliabilityICCIntraclass correlation coefficient is a measure of the correlation between individuals clustered within the same context
Predictive value+LRThe probability that a positive test would be expected in a patient divided by the probability that a positive test would be expected in a patient without a disease
−LRThe probability of a patient testing negative who has a disease divided by the probability of a patient testing negative who does not have a disease

Assessment criteria included in the TOPSIS model.

Corvis ST parameters ranked based on the TOPSIS evaluation and AUC scores

After the TOPSIS model analysis, the Corvis ST DCRs were ranked based on the comprehensive evaluation score (Table 3a). Additionally, the ranking according to the AUC score is presented in Table 3b. The results of the ROC analysis and ICC scores for all included Corvis ST parameters are provided in Supplementary Table S2 and Supplementary Table S3, respectively.

TABLE 3

ParameterAUCSensitivitySpecificity+LR−LRICCUSNS
A1 Time [ms]0.87181.13092.45010.7500.2000.8260.9210.062
SP HC0.70143.40096.23011.5000.5900.7440.7800.052
SP A10.77262.26092.4508.2500.4100.6020.7180.048
DA Ratio Max (2 mm)0.66043.40094.3407.6700.6000.7960.6700.045
CBI0.71052.83092.4507.0000.5100.8280.6670.045
Integrated Radius [mm]0.65149.06086.7903.7100.5900.8000.5140.034
SSI20.73364.15077.3602.8300.4600.8120.5020.034
DA Ratio Max (1 mm)0.64450.94086.7903.8600.5700.5830.4910.033
A1 Velocity [m/s]0.60245.28086.7903.4300.6300.7370.4880.033
cCBI0.67456.60079.2502.7300.5500.8150.4870.033
HC Time [ms]0.72860.38081.1303.2000.4900.5430.4750.032
HC Deflection Area [mm2]0.64754.72073.5802.0700.6200.8220.4590.031
ARTh0.59156.60069.8101.8700.6200.8840.4580.031
Radius [mm]0.68952.83081.1302.8000.5800.6120.4580.031
PachySlope [µm]0.62564.15058.4901.5500.6100.8930.4510.030
Deformation Amp. Max [mm]0.65177.36047.1701.4600.4800.8290.4500.030
HC Deformation Amp. [mm]0.65177.36047.1701.4600.4800.8290.4500.030
Peak Dist. [mm]0.58792.45024.5301.2200.3100.7830.4410.030
Max Inverse Radius [mm^-1]0.67977.36058.4901.8600.3900.5670.4370.029
HC Deflection Length [mm]0.64645.28084.9103.0000.6400.3000.4100.027
A1 Deflection Velocity [m/s]0.51739.62079.2501.9100.7600.4550.3790.025
A1 Deformation Amp. [mm]0.55935.85086.7902.7100.7400.1410.3680.025
A1 Deflection Area [mm2]0.51430.19081.1301.6000.8600.4720.3650.024
dArc Length Max [mm]0.61460.38067.9201.8800.5800.1570.3610.024
Deflection Amp. Max [mm]0.61437.74083.0202.2200.7500.1360.3530.024
Deflection Amp Max [ms]0.54386.79028.3001.2100.4700.1360.3510.024
A1 dArc Length [mm]0.53420.75092.4502.7500.8600.0060.3490.023
HC Deflection Amp. [mm]0.63145.28075.4701.8500.7300.1660.3460.023
A1 Deflection Amp. [mm]0.62860.38066.0401.7800.6000.0140.3440.023
A1 Deflection Length [mm]0.59062.26060.3801.5700.6300.1040.3390.023
HC dArc Length [mm]0.57249.06071.7001.7300.7100.1440.3390.023
SSI0.50381.1305.6600.8603.3300.5780.3070.021
ParameterAUCSensitivitySpecificity+LR−LRICC
A1 Time [ms]0.87181.13092.45010.7500.2000.826
SP A10.77262.26092.4508.2500.4100.602
SSI20.73364.15077.3602.8300.4600.812
HC Time [ms]0.72860.38081.1303.2000.4900.543
CBI0.71052.83092.4507.0000.5100.828
SP HC0.70143.40096.23011.5000.5900.744
Radius [mm]0.68952.83081.1302.8000.5800.612
Max Inverse Radius [mm^-1]0.67977.36058.4901.8600.3900.567
cCBI0.67456.60079.2502.7300.5500.815
DA Ratio Max (2 mm)0.66043.40094.3407.6700.6000.796
Integrated Radius [mm]0.65149.06086.7903.7100.5900.800
Deformation Amp. Max [mm]0.65177.36047.1701.4600.4800.829
HC Deformation Amp. [mm]0.65177.36047.1701.4600.4800.829
HC Deflection Area [mm2]0.64754.72073.5802.0700.6200.822
HC Deflection Length [mm]0.64645.28084.9103.0000.6400.300
DA Ratio Max (1 mm)0.64450.94086.7903.8600.5700.583
HC Deflection Amp. [mm]0.63145.28075.4701.8500.7300.166
A1 Deflection Amp. [mm]0.62860.38066.0401.7800.6000.014
PachySlope [µm]0.62564.15058.4901.5500.6100.893
dArc Length Max [mm]0.61460.38067.9201.8800.5800.157
Deflection Amp. Max [mm]0.61437.74083.0202.2200.7500.136
A1 Velocity [m/s]0.60245.28086.7903.4300.6300.737
ARTh0.59156.60069.8101.8700.6200.884
A1 Deflection Length [mm]0.59062.26060.3801.5700.6300.104
Peak Dist. [mm]0.58792.45024.5301.2200.3100.783
HC dArc Length [mm]0.57249.06071.7001.7300.7100.144
A1 Deformation Amp. [mm]0.55935.85086.7902.7100.7400.141
Deflection Amp Max [ms]0.54386.79028.3001.2100.4700.136
A1 dArc Length [mm]0.53420.75092.4502.7500.8600.006
A1 Deflection Velocity [m/s]0.51739.62079.2501.9100.7600.455
A1 Deflection Area [mm2]0.51430.19081.1301.6000.8600.472
SSI0.50381.1305.6600.8603.3300.578

(a) Order of Corvis ST parameters based on the TOPSIS model comprehensive evaluation score. (b) Order of Corvis ST parameters based on the AUC score.

Notes: US, unnormalized score, NS, normalized score.

VAE-NT index formula

Backward stepwise logistic regression was used to analyze the top seven parameters based on the TOPSIS score, and the following formula was derived:where Beta = B0+B1 * A1 Time + B2 * SP A1 + B3 *SP HC + B4* DA Ratio Max (2 mm) +B5* DA Ratio Max (1 mm) +B6* SSI2+B7* Integrated Radius.

Moreover, B0 = 117.602, B1 = −19.943, B2 = −0.218, B3 = 1.594, B4 = −5.659, B5 = 52.669, B6 = −11.914, and B7 = −1.427. The results of logistic regression are shown in Table 4.

TABLE 4

ParameterβS.E.WalddfSigExp(β)
A1 Time [ms]−19.9435.21914.6021<0.0010.000
SP HC1.5940.46111.98210.0014.925
SP A1−0.2180.06013.0841<0.0010.804
DA Ratio Max (2 mm)−5.6592.7614.20110.0400.003
DA Ratio Max (1 mm)52.66921.0046.28810.0127.477E+22
SSI2−11.9145.8734.11510.0430.000
Integrated Radius [mm]−1.4270.7333.78710.0520.240
Constant117.60232.82212.8381<0.0011.185E+51

Variables in the equation based on TOPSIS-selected parameters.

The regression results of the top seven parameters are based on the AUC value, and the final equation only contains A1 Time, SP HC, SP A1, and HC Time. The details are shown in Supplementary Table S4.

Assessment of the VAE-NT detection index

The created VAE-NT index was then tested for validity, reliability, and predictive value in diagnosing VAE-NT from a normal cornea. In the training dataset, the AUC value of the VAE-NT index was 0.971, with a sensitivity of 96.230%, a specificity of 94.340%, and a cutoff value of 0.425. The ICC of VAE-NT was 0.777. Detailed results are presented in Table 5a. Additionally, the AUC value of the composite index based on the AUC value was 0.958, with a sensitivity of 90.570%, a specificity of 92.450%, and a cutoff value of 0.434. Supplementary Table S5 provides detailed results. The AUC value of the tomographic and biomechanical index (TBI) was 0.688, with a sensitivity of 43.400% and a specificity of 88.680%. Detailed results are presented in Table 5b. In the validation dataset, the AUC value of the VAE-NT index was 0.980. The sensitivity and specificity of the VAE-NT index were 93.180% and 95.920%, respectively. The positive and negative likelihood ratios of VAE-NT were 22.830 and 0.071, respectively. The ICC of VAE-NT was 0.835, and the accuracy was 94.624%. Detailed results are presented in Table 5a. For TBI, the receiver operating characteristic curve analysis showed an AUC value of 0.808, with a cutoff value of 0.309, an accuracy of 76.087%, a specificity of 77.550%, and a sensitivity of 72.730%. The positive and negative likelihood ratios were 3.240 and 0.350, respectively. The ICC was 0.881. Detailed results are presented in Table 5b. The ROC curves of VAE-NT and TBI in both the training and validation datasets are shown in Figures 2–5.

TABLE 5

DatasetOmnibus test of model coefficientsHosmer and Lemeshow testModel accuracy (%)AUCSensitivity (%)Specificity (%)+LR−LRICC
Training dataset<0.0010.27295.2830.97196.23094.34017.0000.0400.777
Validation dataset--94.6240.98093.18095.92022.8300.0710.835
DatasetAccuracy (%)AUCSensitivity (%)Specificity (%)+LR−LRICC
Training dataset66.0380.68843.40088.6803.8300.6400.751
Validation dataset76.0870.80872.73077.5503.2400.3500.881

(a) Logistic regression results and diagnostic effectiveness evaluation of the VAE-NT index. (b) Diagnostic effectiveness evaluation of TBI.

FIGURE 2

FIGURE 3

FIGURE 4

FIGURE 5

Discussion

In recent years, several composite diagnostic indices have been developed to diagnose keratoconus. Vinciguerra et al. (2016) introduced the CBI to differentiate between normal and keratoconus corneas. Ambrósio et al. (2017) developed the TBI, which exhibited superior diagnostic performance. Subsequent studies have consistently confirmed their high diagnostic accuracy for keratoconus; however, their performance in detecting FFKC has been comparatively limited (Tian et al., 2021b; Asroui et al., 2022; Wallace et al., 2023; Koh et al., 2019; Ferreira-Mendes et al., 2019; Kataria et al., 2019; Song et al., 2022). Although studies have shown that these indices can achieve an AUC value greater than 0.90, there remains a need to further improve sensitivity, specificity, likelihood ratios (positive and negative), and ICCs (Zhang et al., 2020). This situation requires clinicians to consider multiple diagnostic indices simultaneously when making clinical decisions, thereby increasing the complexity and burden of diagnosis. To address this issue, we developed a novel diagnostic index, VAE-NT. The VAE-NT index demonstrated robust diagnostic performance, achieving high AUC, sensitivity, specificity, favorable likelihood ratios (positive and negative), and ICC values in both the training and validation datasets. Furthermore, compared with the TBI, the VAE-NT index exhibited superior diagnostic efficiency and predictive performance in our dataset.

The primary distinction of the VAE-NT index from previous indices lies in its parameter selection methodology, which is based on the TOPSIS model. Unlike the CBI and TBI, which rely solely on AUC-based parameter selection, the VAE-NT index incorporates a TOPSIS-based multi-criteria evaluation approach (Ambrósio et al., 2017; Vinciguerra et al., 2016). In this study, the TOPSIS model was utilized to evaluate the Corvis ST parameters in terms of validity, predictive power, and reliability. Based on this evaluation, the parameters were ranked according to their TOPSIS-normalized scores. The top seven parameters, determined by the available sample size, were then used in a binary logistic regression analysis, and the model with the optimal diagnostic performance was retained. This study represents a novel integration of the MCDA technique and statistical modeling for the comprehensive evaluation of diagnostic parameters in keratoconus detection.

Corneal thickness is a well-recognized factor influencing the diagnosis of keratoconus. It can significantly affect the diagnostic performance of parameters used to differentiate VAE-NT or KC from normal corneas. It is important to note that corneal thickness was not matched during the development of the CBI and TBI (6, 34). In the present study, the PSM method was applied to balance the distribution of central corneal thickness and other potential confounders between the normal and VAE-NT groups (Kane et al., 2020). The matching process also ensured a consistent sample size across both groups. This approach improves the reliability of the subsequent ROC analysis by meeting sample size requirements and effectively controlling for the potential confounding effects of corneal thickness on diagnostic outcomes. In the validation dataset, which exhibited significant differences in central corneal thickness, the VAE-NT index maintained superior discriminatory ability between VAE-NT and normal eyes. Recently, Ambrósio et al. (2023) introduced the TBI version 2 (TBIV2), which showed high diagnostic accuracy for detecting VAE-NT, with an AUC value of 0.945 (DeLong test, P < 0.0001), a sensitivity of 84.4%, and a specificity of 90.1%.

In our dataset, TBIV2 demonstrated slightly lower diagnostic performance than the VAE-NT index, although it remains highly effective in identifying VAE-NT. Notably, the VAE-NT index is exclusively based on Corvis ST parameters, which reflect corneal biomechanical properties, whereas TBIV2 integrates both biomechanical and tomographic data. These findings highlight the important contribution of biomechanical information to VAE-NT diagnosis, aligning with the well-established role of collagen disruption in keratoconus pathogenesis (Rabinowitz, 1998; Roberts and Dupps, 2014). The components of the VAE-NT index—including SP A1, SSI2, A1 Time, SP HC, DA Ratio Max (2 mm), DA Ratio Max (1 mm), and Integrated Radius—were selected based on their demonstrated effectiveness in distinguishing VAE-NT from normal corneas (Miao et al., 2023; Tian et al., 2021b; Huo et al., 2023; Zhang et al., 2022). Moreover, SP HC, DA Ratio Max (2 mm), DA Ratio Max (1 mm), and Integrated Radius have also been reported to differ significantly between normal and keratoconus corneas (Song et al., 2023).

This study has several limitations. First, the single-center design and relatively moderate sample size (n = 106) may limit the generalizability of the findings. Future studies involving larger, multi-center cohorts are warranted to validate the results. Second, the retrospective nature of data collection may introduce selection bias, highlighting the need for well-designed prospective studies. Third, although the VAE-NT index demonstrated promising diagnostic performance, its validity, reliability, and predictive values require further evaluation using independent external datasets from other clinical centers.

Conclusion

In this study, the VAE-NT index was developed to distinguish VAE-NT from normal eyes. It demonstrated high sensitivity, specificity, AUC, favorable likelihood ratios, and good reliability, indicating strong diagnostic potential. The use of the TOPSIS model enabled a comprehensive evaluation of diagnostic indicators, facilitating the selection of features with superior overall diagnostic strength and providing clinicians with a more objective decision-making reference. Further validation in larger, more diverse populations and with longer follow-ups is necessary to support clinical implementation.

Statements

Data availability statement

The datasets presented in this article are not readily available because the dataset is not publicly available due to ethical/privacy restrictions. Requests to access the datasets should be directed to landy.yang@foxmail.com.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Eye Hospital, Wenzhou Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

LY: Data curation, Methodology, Validation, Visualization, Conceptualization, Writing – original draft, Formal analysis, Software, Investigation, Writing – review and editing. HX: Investigation, Formal analysis, Writing – original draft, Data curation. HJ: Writing – review and editing, Data curation. JZ: Funding acquisition, Project administration, Resources, Supervision, Writing – review and editing, Investigation. SC: Supervision, Investigation, Writing – review and editing, Conceptualization, Methodology, Resources, Visualization, Project administration, Validation.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by grants from Huadong Hospital affiliated with Fudan University, China (HDLC2022009). The funding body had no role in the design of the study, data collection and analysis, interpretation of data, or writing of the manuscript.

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.

The handling editor YS declared a shared parent affiliation with the author(s) LY, HJ, and JZ at the time of review.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2025.1598546/full#supplementary-material

References

  • 1

    AmbrósioR.LopesB. T.Faria-CorreiaF.SalomãoM. Q.BührenJ.RobertsC. J.et al (2017). Integration of scheimpflug-based corneal tomography and biomechanical assessments for enhancing ectasia detection. J. Refract Surg.33 (7), 434443. 10.3928/1081597x-20170426-02

  • 2

    AmbrósioR.MachadoA. P.LeãoE.LyraJ. M. G.SalomãoM. Q.EsporcatteL. G. P.et al (2023). Optimized artificial intelligence for enhanced ectasia detection using scheimpflug-based corneal tomography and biomechanical data. Am. J. Ophthalmol.251, 126142. 10.1016/j.ajo.2022.12.016

  • 3

    ArbelaezM. C.VersaciF.VestriG.BarboniP.SaviniG. (2012). Use of a support vector machine for keratoconus and subclinical keratoconus detection by topographic and tomographic data. Ophthalmology119 (11), 22312238. 10.1016/j.ophtha.2012.06.005

  • 4

    AsrouiL.DagherS. A.ElsheikhA.LopesB. T.RobertsC. J.AssouadM.et al (2022). Biomechanical evaluation of topographically and tomographically normal fellow eyes of patients with keratoconus. J. Refract Surg.38 (5), 318325. 10.3928/1081597x-20220225-01

  • 5

    Ferreira-MendesJ.LopesB. T.Faria-CorreiaF.SalomãoM. Q.Rodrigues-BarrosS.AmbrósioR. (2019). Enhanced ectasia detection using corneal tomography and biomechanics. Am. J. Ophthalmol.197, 716. 10.1016/j.ajo.2018.08.054

  • 6

    GirardM. J. A.DuppsW. J.BaskaranM.ScarcelliG.YunS. H.QuigleyH. A.et al (2015). Translating ocular biomechanics into clinical practice: current state and future prospects. Curr. Eye Res.40 (1), 118. 10.3109/02713683.2014.914543

  • 7

    HenriquezM. A.HadidM.IzquierdoL. (2020). A systematic review of subclinical keratoconus and forme fruste keratoconus. J. Refract Surg.36 (4), 270279. 10.3928/1081597x-20200212-03

  • 8

    HerberR.PillunatL. E.RaiskupF. (2021). Development of a classification system based on corneal biomechanical properties using artificial intelligence predicting keratoconus severity. Eye Vis.8 (1), 2111. 10.1186/s40662-021-00244-4

  • 9

    HuangK. W.HuangJ. H.TzengG. H. (2016). New hybrid multiple attribute decision-making model for improving competence sets: enhancing a company’s core competitiveness. Sustain8 (2), 175. 10.3390/su8020175

  • 10

    HuoY.ChenX.CaoH.LiJ.HouJ.WangY. (2023). Biomechanical properties analysis of forme fruste keratoconus and subclinical keratoconus. Graefe’s Arch. Clin. Exp. Ophthalmol.261 (5), 13111320. 10.1007/s00417-022-05916-y

  • 11

    HwangC.-L.YoonK. (1981). “Multiple attribute decision making methods and applications”. New York: Springer-Verlag Berlin Heidelb. 10.1007/978-3-642-48318-9

  • 12

    HwangE. S.Perez-StraziotaC. E.KimS. W.SanthiagoM. R.RandlemanJ. B. (2018). Distinguishing highly asymmetric keratoconus eyes using combined scheimpflug and spectral-domain OCT analysis. Ophthalmology125 (12), 18621871. 10.1016/j.ophtha.2018.06.020

  • 13

    JahanshahlooG. R.LotfiF. H.IzadikhahM. (2006). Extension of the TOPSIS method for decision-making problems with fuzzy data. Appl. Math. Comput.181 (2), 15441551. 10.1016/j.amc.2006.02.057

  • 14

    KaneL. T.FangT.GalettaM. S.GoyalD. K. C.NicholsonK. J.KeplerC. K.et al (2020). Propensity score matching A statistical method. 33(3):120122. 10.1097/bsd.0000000000000932

  • 15

    KatariaP.PadmanabhanP.GopalakrishnanA.PadmanabanV.MahadikS.AmbrósioR. (2019). Accuracy of scheimpflug-derived corneal biomechanical and tomographic indices for detecting subclinical and mild keratectasia in a South Asian population. J. Cataract. Refract Surg.45 (3), 328336. 10.1016/j.jcrs.2018.10.030

  • 16

    KeeneyR. L. (1982). Feature article—decision analysis: an overview. Oper. Res.30, 803838. 10.1287/opre.30.5.803

  • 17

    KohS.AmbrósioR.InoueR.MaedaN.MikiA.NishidaK. (2019). Detection of subclinical corneal ectasia using corneal tomographic and biomechanical assessments in a Japanese population. J. Refract Surg.35 (6), 383390. 10.3928/1081597x-20190417-01

  • 18

    KooT. K.LiM. Y. (2016). A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med.15 (2), 155163. 10.1016/j.jcm.2016.02.012

  • 19

    LuN. J.ElsheikhA.RozemaJ. J.HafeziN.AslanidesI. M.HillenM.et al (2022). Combining spectral-domain OCT and air-puff tonometry analysis to diagnose keratoconus. J. Refract Surg.38 (6), 374380. 10.3928/1081597x-20220414-02

  • 20

    MandrekarJ. N. (2010). Receiver operating characteristic curve in diagnostic test assessment. J. Thorac. Oncol.5 (9), 13151316. 10.1097/jto.0b013e3181ec173d

  • 21

    MiaoY.MaX.QuZ.EliasyA.WuB. W.XuH.et al (2023). Performance of corvis ST parameters including updated stress-strain index in differentiating between normal, forme-fruste, subclinical, and clinical keratoconic eyes. Am. J. Ophthalmol.258, 196207. 10.1016/j.ajo.2023.10.015

  • 22

    MullerR.ButtnerP. (1994). A critical discussion of intraclass correlation coefficients. Stat. Med.13, 24652476. 10.1002/sim.4780132310

  • 23

    PiñeroD. P.AlcónN. (2015). Corneal biomechanics: a review. Clin. Exp. Optom.98 (2), 107116. 10.1111/cxo.12230

  • 24

    RabinowitzY. S. (1998). Keratoconus. Surv. Ophthalmol.42 (4), 297319. 10.1016/s0039-6257(97)00119-7

  • 25

    RammL.HerberR.SpoerlE.RaiskupF.PillunatL. E.TeraiN. (2019). Intraocular pressure measurement using ocular response analyzer, dynamic contour tonometer, and scheimpflug analyzer corvis ST. J. Ophthalmol.2019, 38796513879659. 10.1155/2019/3879651

  • 26

    Ramón-CanulL. G.Margarito-CarrizalD. L.Limón-RiveraR.Morales-CarrreraU. A.Rodríguez-BuenfilI. M.Ramírez-SucreM. O.et al (2021). Technique for order of preference by similarity to ideal solution (TOPSIS) method for the generation of external preference mapping using rapid sensometric techniques. J. Sci. Food Agric.101 (8), 32983307. 10.1002/jsfa.10959

  • 27

    RenS.XuL.FanQ.GuY.YangK. (2021). Accuracy of new corvis ST parameters for detecting subclinical and clinical keratoconus eyes in a Chinese population. Sci. Rep.11 (1), 4962. 10.1038/s41598-021-84370-y

  • 28

    RobertsC. J.DuppsW. J.Jr (2014). Biomechanics of corneal ectasia and biomechanical treatments. J. Cataract. Refract Surg.40 (6), 991998. 10.1016/j.jcrs.2014.04.013

  • 29

    RubertiJ. W.Sinha RoyA.RobertsC. J. (2011). Corneal biomechanics and biomaterials. Annu. Rev. Biomed. Eng.13, 269295. 10.1146/annurev-bioeng-070909-105243

  • 30

    SongP.RenS.LiuY.LiP.ZengQ. (2022). Detection of subclinical keratoconus using a novel combined tomographic and biomechanical model based on an automated decision tree. Sci. Rep.12 (0123456789), 53165319. 10.1038/s41598-022-09160-6

  • 31

    SongY.FengY.QuM.MaQ.TianH.LiD.et al (2023). Analysis of the diagnostic accuracy of belin/ambrósio enhanced ectasia and corvis ST parameters for subclinical keratoconus. Int. Ophthalmol.43 (5), 14651475. 10.1007/s10792-022-02543-8

  • 32

    TalukderB.HipelK. W.vanLoonG. W. (2018). Using multi-criteria decision analysis for assessing sustainability of agricultural systems. Sustain Dev.26 (6), 781799. 10.1002/sd.1848

  • 33

    TianL.KoM. W.WangL. K.ZhangJ. Y.LiT. J.HuangY. F.et al (2014). Assessment of ocular biomechanics using dynamic ultra high-speed Scheimpflug imaging in keratoconic and normal eyes. J. Refract Surg.30 (11), 785791. 10.3928/1081597x-20140930-01

  • 34

    TianL.QinX.ZhangH.ZhangD.GuoL. L.ZhangH. X.et al (2021b). A potential screening index of corneal biomechanics in healthy subjects, forme fruste keratoconus patients and clinical keratoconus patients. Front. Bioeng. Biotechnol.9 (9(December), 766605766610. 10.3389/fbioe.2021.766605

  • 35

    TianL.ZhangD.GuoL.QinX.ZhangH.ZhangH.et al (2021a). Comparisons of corneal biomechanical and tomographic parameters among thin normal cornea, forme fruste keratoconus, and mild keratoconus. Eye Vis.8 (1), 4411. 10.1186/s40662-021-00266-y

  • 36

    VinciguerraR.AmbrósioR.ElsheikhA.RobertsC. J.LopesB.MorenghiE.et al (2016). Detection of keratoconus with a new biomechanical index. J. Refract Surg.32 (12), 803810. 10.3928/1081597x-20160629-01

  • 37

    WallaceH. B.VellaraH. R.GokulA.McGheeC. N. J.MeyerJ. J. (2023). Comparison of ectasia detection in early keratoconus using scheimpflug-based corneal tomography and biomechanical assessments. Cornea42 (12), 15281535. 10.1097/ico.0000000000003273

  • 38

    WangY. M.ChanT. C. Y.YuM.JhanjiV. (2017). Comparison of corneal dynamic and tomographic analysis in normal, forme fruste keratoconic, and keratoconic eyes. J. Refract Surg.33 (9), 632638. 10.3928/1081597x-20170621-09

  • 39

    XianY.ZhaoY.SunL.ZhangX.DingL.LiuZ.et al (2023). Comparison of bilateral differential characteristics of corneal biomechanics between keratoconus and normal eyes. Front. Bioeng. Biotechnol.11 (June), 116322310. 10.3389/fbioe.2023.1163223

  • 40

    YoonK. P.HwangC. L. (1995). Multiple attribute decision making: anintroduction. Berlin Heidelberg: Springer. 10.1007/978-3-642-48318-9

  • 41

    ZhangH.EliasyA.LopesB.AbassA.VinciguerraR.VinciguerraP.et al (2021b). Stress-strain index map: a new way to represent corneal material stiffness. Front. Bioeng. Biotechnol.9, 640434. 10.3389/fbioe.2021.640434

  • 42

    ZhangH.TianL.GuoL.QinX.ZhangD.LiL.et al (2021a). Comprehensive evaluation of corneas from normal, forme fruste keratoconus and clinical keratoconus patients using morphological and biomechanical properties. Int. Ophthalmol.41 (4), 12471259. 10.1007/s10792-020-01679-9

  • 43

    ZhangH.ZhangX.HuaL.LiL.TianL.ZhangX.et al (2022). An exploratory analysis of forme fruste keratoconus sensitivity diagnostic parameters. Int. Ophthalmol.42 (8), 24732481. 10.1007/s10792-022-02246-0

  • 44

    ZhangM.ZhangF.LiY.SongY.WangZ. (2020). Early diagnosis of keratoconus in Chinese myopic eyes by combining corvis ST with pentacam. Curr. Eye Res.45 (2), 118123. 10.1080/02713683.2019.1658787

  • 45

    ZhangP.YangL.MaoY.ZhangX.ChengJ.MiaoY.et al (2024). CorNet: autonomous feature learning in raw corvis ST data for keratoconus diagnosis via residual CNN approach. Comput. Biol. Med.172, 108286. 10.1016/j.compbiomed.2024.108286

Summary

Keywords

very asymmetric ectasia with normal tomography, dynamic corneal response parameters, technique for order of preference by similarity to ideal solution, keratoconus, corneal biomechanics

Citation

Yang L, Xu H, Jiang H, Zhu J and Chen S (2025) Development and validation of the VAE-NT index: a novel biomechanical parameter for distinguishing subclinical corneal abnormalities. Front. Bioeng. Biotechnol. 13:1598546. doi: 10.3389/fbioe.2025.1598546

Received

23 March 2025

Accepted

20 June 2025

Published

16 July 2025

Volume

13 - 2025

Edited by

Yang Shen, Fudan University, China

Reviewed by

Xuefeng Shi, Tianjin Eye Hospital, China

Zhipeng Yan, Third Hospital of Hebei Medical University, China

Updates

Copyright

*Correspondence: Jingyin Zhu, ; Shihao Chen,

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics