In vivo biomechanical mapping of normal and keratoconus corneas

The Pathophysiology of Keratoconus

PURPOSE

Keratoconus is a progressive disease characterized by changes in corneal shape, resulting in loss of visual function. There remains a lack of comprehensive understanding regarding its underlying pathophysiology. This review aims to bridge this gap by exploring structural failures and inflammatory processes involved in the etiology and progression of keratoconus.

METHODS

A literature review was conducted using PubMed and Google Scholar databases, screening for articles published in English using the keyword combinations of "keratoconus" with "pathophysiology," "pathology," "metabolism," "inflammatory," "oxidative stress," "cytokines," "enzymes," "collagen," and "cornea." Articles published between January 1, 1970, and June 1, 2023, were queried and reviewed, with greater emphasis placed on more recent data. Fifty-six relevant studies were examined to develop a thorough review of the pathophysiological mechanisms at play in keratoconus.

RESULTS

Biomechanical structural failures in the cornea seem to be the primary militating factors in keratoconus etiology and progression. These include disruptions in the arrangement in the collagen lamellae, a decrease in collagen levels, a decrease in natural collagen crosslinking, and changes in lysosomal enzyme activity. Immunologic changes have also been identified in keratoconus, challenging the traditional view of the condition as noninflammatory. Elevated levels of proinflammatory cytokines like IL-1b, IL-6, IL-17, and TNF-α have been observed, along with increased apoptosis of keratocytes. Increased oxidative stress leads to the activation of collagenase and gelatinase enzymes.

CONCLUSIONS

Keratoconus is a complex condition influenced by both structural defects and inflammatory processes. Understanding these mechanisms can inform clinical management and potentially lead to more effective treatments.

CorNet: Autonomous feature learning in raw Corvis ST data for keratoconus diagnosis via residual CNN approach

PURPOSE

To ascertain whether the integration of raw Corvis ST data with an end-to-end CNN can enhance the diagnosis of keratoconus (KC).

METHOD

The Corvis ST is a non-contact device for in vivo measurement of corneal biomechanics. The CorNet was trained and validated on a dataset consisting of 1786 Corvis ST raw data from 1112 normal eyes and 674 KC eyes. Each raw data consists of the anterior and posterior corneal surface elevation during air-puff induced dynamic deformation. The architecture of CorNet utilizes four ResNet-inspired convolutional structures that employ 1 × 1 convolution in identity mapping. Gradient-weighted Class Activation Mapping (Grad-CAM) was adopted to visualize the attention allocation to diagnostic areas. Discriminative performance was assessed using metrics including the AUC of ROC curve, sensitivity, specificity, precision, accuracy, and F1 score.

RESULTS

CorNet demonstrated outstanding performance in distinguishing KC from normal eyes, achieving an AUC of 0.971 (sensitivity: 92.49%, specificity: 91.54%) in the validation set, outperforming the best existing Corvis ST parameters, namely the Corvis Biomechanical Index (CBI) with an AUC of 0.947, and its updated version for Chinese populations (cCBI) with an AUC of 0.963. Though the ROC curve analysis showed no significant difference between CorNet and cCBI (p = 0.295), it indicated a notable difference between CorNet and CBI (p = 0.011). The Grad-CAM visualizations highlighted the significance of corneal deformation data during the loading phase rather than the unloading phase for KC diagnosis.

CONCLUSION

This study proposed an end-to-end CNN approach utilizing raw biomechanical data by Corvis ST for KC detection, showing effectiveness comparable to or surpassing existing parameters provided by Corvis ST. The CorNet, autonomously learning comprehensive temporal and spatial features, demonstrated a promising performance for advancing KC diagnosis in ophthalmology.

Subclinical Keratoconus Detection and Characterization Using Motion-Tracking Brillouin Microscopy

PURPOSE

To characterize focal biomechanical alterations in subclinical keratoconus (SKC) using motion-tracking (MT) Brillouin microscopy and evaluate the ability of MT Brillouin metrics to differentiate eyes with SKC from normal control eyes.

DESIGN

Prospective cross-sectional study.

PARTICIPANTS

Thirty eyes from 30 patients were evaluated, including 15 eyes from 15 bilaterally normal patients and 15 eyes with SKC from 15 patients.

METHODS

All patients underwent Scheimpflug tomography and MT Brillouin microscopy using a custom-built device. Mean and minimum MT Brillouin values within the anterior plateau region and anterior 150 μm were generated. Scheimpflug metrics evaluated included inferior-superior (IS) value, maximum keratometry (Kmax), thinnest corneal thickness, asymmetry indices, Belin/Ambrosio display total deviation, and Ambrosio relational thickness. Receiver operating characteristic (ROC) curves were generated for all Scheimpflug and MT Brillouin metrics evaluated to determine the area under the ROC curve (AUC), sensitivity, and specificity for each variable.

MAIN OUTCOME MEASURES

Discriminative performance based on AUC, sensitivity, and specificity.

RESULTS

No significant differences were found between groups for age, sex, manifest refraction spherical equivalent, corrected distance visual acuity, Kmax, or KISA% index. Among Scheimpflug metrics, significant differences were found between groups for thinnest corneal thickness (556 μm vs. 522 μm; P < 0.001), IS value (0.29 diopter [D] vs. 1.05 D; P < 0.001), index of vertical asymmetry (IVA; 0.10 vs. 0.19; P < 0.001), and keratoconus index (1.01 vs. 1.05; P < 0.001), and no significant differences were found for any other Scheimpflug metric. Among MT Brillouin metrics, clear differences were found between control eyes and eyes with SKC for mean plateau (5.71 GHz vs. 5.68 GHz; P < 0.0001), minimum plateau (5.69 GHz vs. 5.65 GHz; P < 0.0001), mean anterior 150 μm (5.72 GHz vs. 5.68 GHz; P < 0.0001), and minimum anterior 150 μm (5.70 GHz vs. 5.66 GHz; P < 0.001). All MT Brillouin plateau and anterior 150 μm mean and minimum metrics fully differentiated groups (AUC, 1.0 for each), whereas the best performing Scheimpflug metrics were keratoconus index (AUC, 0.91), IS value (AUC, 0.89), and IVA (AUC, 0.88).

CONCLUSIONS

Motion-tracking Brillouin microscopy metrics effectively characterize focal corneal biomechanical alterations in eyes with SKC and clearly differentiated these eyes from control eyes, including eyes that were not differentiated accurately using Scheimpflug metrics.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Corneal Stress Distribution Evolves from Thickness-Driven in Normal Corneas to Curvature-Driven with Progression in Keratoconus

Purpose

To introduce the novel parameter of Corneal Contribution to Stress (CCS) and compare stress distribution patterns between keratoconus (KCN) and normal corneas.

Design

Prospective, observational, cross-sectional study.

Participants

The study included 66 eyes of 40 subjects diagnosed with KCN and 155 left eyes from 155 normal control (NRL) subjects.

Methods

Tomography was obtained to calculate the newly proposed CCS, defined according to the hoop stress formula without intraocular pressure, R/2t, where R is the radius of curvature and t is the thickness. CCS maps were calculated from pachymetry and tangential curvature maps. Custom software identified the 2-mm-diameter zones of greatest curvature (Cspot-max), thinnest pachymetry (Pach-min), greatest stress (CCSmax), and lowest stress (CCSmin). Stress difference (CCSdiff) was calculated as CCSmax - CCSmin. Distances between Cspot-max vs. Pach-min, vs. CCSmax, and vs. CCSmin, as well as between Pach-min vs. CCSmax and vs. CCSmin, were calculated. t tests were performed between cohorts, and paired t tests were performed within cohorts. Univariate linear regression analyses were performed between parameters and distances. The significance threshold was P < 0.05.

Main Outcome Measures

Corneal stress parameters, corneal features of maximum curvature, minimum thickness, and distances between corneal stress parameters and corneal features.

Results

CCSmax was significantly closer to Pach-min (0.79 ± 0.92) and Cspot-max (2.04 ± 0.85) than CCSmin (3.17 ± 0.38, 2.73 ± 1.53, respectively) in NRL, P < 0.0001, whereas CCSmin was significantly closer to Cspot-max (1.35 ± 1.43) than CCSmax (2.52 ± 0.72) in KCN, P < 0.0001. Cspot-max (severity) was significantly related to CCSdiff in KCN (P < 0.0001; R2 = 0.5882) with a weak relationship in NRL (P < 0.0080, R2 = 0.0451). Cspot-max was significantly related to the distance from Pach-min to CCSmax (P < 0.0001; R2 = 0.3737) without significance in NRL (P = 0.8011).

Conclusions

Corneal stress is driven by thickness in NRL, with greatest stress at thinnest pachymetry and greatest curvature. However, maximum stress moves away from thinnest pachymetry with progression in KCN, and minimum stress is associated with maximum curvature. Severity in KCN is significantly related to greater difference between maximum and minimum stress, consistent with the biomechanical cycle of decompensation.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

In Vivo Evaluation of Corneal Biomechanics Following Cross-Linking Surgeries Using Optical Coherence Elastography in a Rabbit Model of Keratoconus

Purpose

Validation of the feasibility of novel acoustic radiation force optical coherence elastography (ARF-OCE) for the evaluation of biomechanical enhancement of the in vivo model of keratoconus by clinical cross-linking (CXL) surgery.

Methods

Twelve in vivo rabbit corneas were randomly divided into two groups. Both groups were treated with collagenase type II, and a keratoconus model was obtained. Then, the two groups were treated with CXL procedures with different irradiation energy of 15 J and 30 J (CXL-15 J and CXL-30 J, respectively). An ARF-OCE probe with an ultrasmall ultrasound transducer was used to detect the biomechanical properties of cornea. An antisymmetric Lamb wave model was combined with the frequency dispersion relationship to achieve depth-resolved elastography.

Results

Compared with the phase velocity of the Lamb wave in healthy corneas (approximately 3.96 ± 0.27 m/s), the phase velocity of the Lamb wave was lower in the keratoconus region (P < 0.05), with an average value of 3.12 ± 0.12 m/s. Moreover, the corneal stiffness increased after CXL treatment (P < 0.05), and the average phase velocity of the Lamb wave was 4.3 ± 0.19 m/s and 4.54 ± 0.13 m/s after CXL-15 J and CXL-30 J treatment.

Conclusions

The Young's moduli of the keratoconus regions were significantly lower than the healthy corneas. Moreover, the Young's modulus of the keratoconus regions was significantly higher after CXL-30 J treatment than after CXL-15 J treatment. We demonstrated that the ARF-OCE technique has great potential in screening keratoconus and guiding clinical CXL treatment.

Translational Relevance

This work accelerates the clinical translation of OCE systems using ultrasmall ultrasound transducers and is used to guide CXL procedures.

Performance of Corvis ST Parameters Including Updated Stress-Strain Index in Differentiating Between Normal, Forme-Fruste, Subclinical, and Clinical Keratoconic Eyes

PURPOSE

This study seeks to evaluate the ability of the updated stress strain index (SSIv2) and other Corvis ST biomechanical parameters in distinguishing between keratoconus at different disease stages and normal eyes.

DESIGN

Diagnostic accuracy analysis to distinguish disease stages.

METHODS

1084 eyes were included and divided into groups of normal (199 eyes), forme fruste keratoconus (FFKC, 194 eyes), subclinical keratoconus (SKC, 113 eyes), mild clinical keratoconus (CKC-Ⅰ, 175 eyes), moderate clinical keratoconus (CKC-Ⅱ, 204 eyes), and severe clinical keratoconus (CKC-Ⅲ, 199 eyes). Each eye was subjected to a Corvis ST examination to determine the central corneal thickness (CCT), biomechanically corrected intraocular pressure (bIOP), SSIv2 (updated stress-strain index), and other 8 Corvis parameters including the stress-strain index (SSIv1), stiffness parameter at first applanation (SP-A1), first applanation time (A1T), Ambrósio relational thickness to the horizontal profile (ARTh), integrated inverse radius (IIR), maximum deformation amplitude (DAM), ratio between deformation amplitude at the apex and at 2 mm nasal and temporal (DARatio2), and Corvis biomechanical index (CBI). The sensitivity and specificity of these parameters in diagnosing keratoconus were analyzed through receiver operating characteristic curves.

RESULTS

Before and after correction for CCT and bIOP, SSIv2 and ARTh were significantly higher and IIR and CBI were significantly lower in the normal group than in the FFKC group, SKC group and the 3 CKC groups (all P < .05). There were also significant correlations between the values of SSIv2, ARTh, IIR, CBI, and the CKC severity (all P < .05). AUC of SSIv2 was significantly higher than all other Corvis parameters in distinguishing normal eyes from FFKC, followed by IIR, ARTh and CBI.

CONCLUSION

Corvis ST's updated stress-strain index, SSIv2, demonstrated superior performance in differentiating between normal and keratoconic corneas, and between corneas with different keratoconus stages. Similar, but less pronounced, performance was demonstrated by the IIR, ARTh and CBI.

Keratoconus cone location influences ocular biomechanical parameters measured by the Ocular Response Analyzer

BACKGROUND

Keratoconus is characterized by asymmetry in the biomechanical properties of the cornea, with focal weakness in the area of cone formation. We tested the hypothesis that centrally-measured biomechanical parameters differ between corneas with peripheral cones and corneas with central cones.

METHODS

Fifty participants with keratoconus were prospectively recruited. The mean ± standard deviation age was 38 ± 13 years. Axial and tangential corneal topography were analyzed in both eyes, if eligible. Cones in the central 3 mm of the cornea were considered central, and cones outside the central 3 mm were considered peripheral. Each eye was then measured with the Ocular Response Analyzer (ORA) tonometer. T-tests compared differences in ORA-generated waveform parameters between cohorts.

RESULTS

Seventy-eight eyes were analyzed. According to the axial topography maps, 37 eyes had central cones and 41 eyes had peripheral cones. According to the tangential topography maps, 53 eyes had central cones, and 25 eyes had peripheral cones. For the axial-topography algorithm, wave score (WS) was significantly higher in peripheral cones than central cones (inter-cohort difference = 1.27 ± 1.87). Peripheral cones had a significantly higher area of first peak, p1area (1047 ± 1346), area of second peak, p2area (1130 ± 1478), height of first peak, h1 (102 ± 147), and height of second peak, h2 (102 ± 127), than central cones. Corneal hysteresis (CH), width of the first peak, w1, and width of the second peak, w2, did not significantly differ between cohorts. There were similar results for the tangential-topography algorithm, with a significant difference between the cohorts for p1area (855 ± 1389), p2area (860 ± 1531), h1 (81.7 ± 151), and h2 (92.1 ± 131).

CONCLUSIONS

Cone location affects the biomechanical response parameters measured under central loading of the cornea. The ORA delivers its air puff to the central cornea, so the fact that h1 and h2 and that p1area and p2area were smaller in the central cone cohort than in the peripheral cone cohort suggests that corneas with central cones are softer or more compliant centrally than corneas with peripheral cones, which is consistent with the location of the pathology. This result is evidence that corneal weakening in keratoconus is focal in nature and is consistent with localized disruption of lamellar orientation.

Motion-Tracking Brillouin Microscopy Evaluation of Normal, Keratoconic, and Post-Laser Vision Correction Corneas

PURPOSE

To characterize focal biomechanical differences between normal, keratoconic, and post-laser vision correction (LVC) corneas using motion-tracking Brillouin microscopy.

DESIGN

Prospective cross-sectional study.

METHODS

Thirty eyes from 30 patients (10 normal controls [Controls], 10 post-LVC, and 10 stage I or II keratoconus [KC]) had Scheimpflug and motion-tracking Brillouin microscopy imaging using a custom-built device. Mean, maximum (max) and minimum (min) Brillouin shift, spatial standard deviation, and max-min values were compared. Min values were correlated with local Brillouin values at multiple Scheimpflug imaging locations.

RESULTS

Mean (P < .0003), min (P < .00001), spatial standard deviation (P < .01), and max-min (P < .001) were significantly different between the groups. In post hoc pairwise comparisons, the best differentiators for group comparisons were mean (P = .0004) and min (P = .000002) for Controls vs KC, min (P = .0022) and max-min (P = .002) for Controls vs LVC, and mean (P = .0037) and min (P = .0043) for LVC vs KC. Min (area under the receiver operating characteristic = 1.0) and mean (area under the receiver operating characteristic =  0.96) performed well in differentiating Control and KC eyes. Min values correlated best with Brillouin shift values at the thinnest corneal point (r2 = 0.871, P = .001) and maximum keratometry value identified in the tangential curvature map (r2 = 0.840, P = .002).

CONCLUSIONS

Motion-tracking Brillouin microscopy effectively characterized focal corneal biomechanical alterations in LVC and KC and clearly differentiated these groups from Controls. Primary motion-tracking Brillouin metrics performed well in differentiating groups as compared with basic Scheimpflug metrics, in contrast to previous Brillouin studies, and identified focal changes after LVC where prior Brillouin studies did not.

Corneal Biomechanical Measures for Glaucoma: A Clinical Approach

Over the last two decades, there has been growing interest in assessing corneal biomechanics in different diseases, such as keratoconus, glaucoma, and corneal disorders. Given the interaction and structural continuity between the cornea and sclera, evaluating corneal biomechanics may give us further insights into the pathogenesis, diagnosis, progression, and management of glaucoma. Therefore, some authorities have recommended baseline evaluations of corneal biomechanics in all glaucoma and glaucoma suspects patients. Currently, two devices (Ocular Response Analyzer and Corneal Visualization Schiempflug Technology) are commercially available for evaluating corneal biomechanics; however, each device reports different parameters, and there is a weak to moderate agreement between the reported parameters. Studies are further limited by the inclusion of glaucoma subjects taking topical prostaglandin analogues, which may alter corneal biomechanics and contribute to contradicting results, lack of proper stratification of patients, and misinterpretation of the results based on factors that are confounded by intraocular pressure changes. This review aims to summarize the recent evidence on corneal biomechanics in glaucoma patients and insights for future studies to address the current limitations of the literature studying corneal biomechanics.

Localized Refractive Changes Induced by Symmetric and Progressive Asymmetric Intracorneal Ring Segments Assessed with a 3D Finite-Element Model

To build a representative 3D finite element model (FEM) for intracorneal ring segment (ICRS) implantation and to investigate localized optical changes induced by different ICRS geometries, a hyperelastic shell FEM was developed to compare the effect of symmetric and progressive asymmetric ICRS designs in a generic healthy and asymmetric keratoconic (KC) cornea. The resulting deformed geometry was assessed in terms of average curvature via a biconic fit, sagittal curvature (K), and optical aberrations via Zernike polynomials. The sagittal curvature map showed a locally restricted flattening interior to the ring (Kmax -11 to -25 dpt) and, in the KC cornea, an additional local steepening on the opposite half of the cornea (Kmax up to +1.9 dpt). Considering the optical aberrations present in the model of the KC cornea, the progressive ICRS corrected vertical coma (-3.42 vs. -3.13 µm); horizontal coma (-0.67 vs. 0.36 µm); and defocus (2.90 vs. 2.75 µm), oblique trefoil (-0.54 vs. -0.08 µm), and oblique secondary astigmatism (0.48 vs. -0.09 µm) aberrations stronger than the symmetric ICRS. Customized ICRS designs inspired by the underlying KC phenotype have the potential to achieve more tailored refractive corrections, particularly in asymmetric keratoconus patterns.

Sequence variants contributing to dysregulated inflammatory responses across keratoconic cone surface in adolescent patients with keratoconus

Background

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Methods

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

Results

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

Conclusion

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

The Impaired Wound Healing Process Is a Major Factor in Remodeling of the Corneal Epithelium in Adult and Adolescent Patients With Keratoconus

Purpose

Keratoconus (KTCN) is the most common corneal ectasia, characterized by pathological cone formation. Here, to provide an insight into the remodeling of the corneal epithelium (CE) during the course of the disease, we evaluated topographic regions of the CE of adult and adolescent patients with KTCN.

Methods

The CE samples from 17 adult and 6 adolescent patients with KTCN, and 5 control CE samples were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated toward RNA sequencing and MALDI-TOF/TOF Tandem Mass Spectrometry. Data from transcriptomic and proteomic investigations were consolidated with the morphological and clinical findings.

Results

The critical elements of the wound healing process, epithelial-mesenchymal transition, cell-cell communications, and cell-extracellular matrix interactions were altered in the particular corneal topographic regions. Abnormalities in pathways of neutrophils degranulation, extracellular matrix processing, apical junctions, IL, and IFN signaling were revealed to cooperatively disorganize the epithelial healing. Deregulation of the epithelial healing, G2M checkpoints, apoptosis, and DNA repair pathways in the middle CE topographic region in KTCN explains the presence of morphological changes in the corresponding doughnut pattern (a thin cone center surrounded by a thickened annulus). Despite similar morphological characteristics of CE samples in adolescents and adults with KTCN, their transcriptomic features were different. Values of the posterior corneal elevation differentiated adults with KTCN from adolescents with KTCN and correlated with the expression of TCHP, SPATA13, CNOT3, WNK1, TGFB2, and KRT12 genes.

Conclusions

Identified molecular, morphological, and clinical features indicate the effect of impaired wound healing on corneal remodeling in KTCN CE.

In Vivo Corneal Stiffness Mapping by the Stress-Strain Index Maps and Brillouin Microscopy

The study of corneal stiffness in vivo has numerous clinical applications such as the measurement of intraocular pressure, the preoperative screening for iatrogenic ectasia after laser vision correction surgery and the diagnosis and treatment of corneal ectatic diseases such as keratoconus. The localised aspect of the microstructure deterioration in keratoconus leading to local biomechanical softening, corneal bulging, irregular astigmatism and ultimately loss of vision boosted the need to map the corneal stiffness to identify the regional biomechanical failure. Currently, two methods to map the corneal stiffness in vivo are integrated into devices that are either already commercially available or about to be commercialised: the stress-strain index (SSI) maps and the Brillouin Microscopy (BM). The former method produces 2D map of stiffness across the corneal surface, developed through numerical simulations using the corneal shape, its microstructure content, and the deformation behaviour under air-puff excitation. It estimates the whole stress-strain behaviour, making it possible to obtain the material tangent modulus under different intraocular pressure levels. On the other hand, BM produces a 3D map of the corneal longitudinal modulus across the corneal surface and thickness. It uses a low-power near-infrared laser beam and through a spectral analysis of the returned signal, it assesses the mechanical compressibility of the tissue as measured by the longitudinal modulus. In this paper, these two techniques are reviewed, and their advantages and limitations discussed.

Keratoconus: A Biomechanical Perspective

PURPOSE

The relevance of corneal biomechanics and the importance of including it in the clinical assessment of corneal ectasias are being increasingly recognized. The connection between corneal ultrastructure, biomechanical properties, and optical function is exemplified by a condition like keratoconus. Biomechanical instability is seen as the underlying basis for the secondary morphological changes in the cornea. Asymmetric biomechanical weakening is believed to drive progressive corneal steepening and thinning. Biomechanical strengthening is the principle of collagen crosslinking that has been shown to effectively arrest progression of the keratoconus. Corneal biomechanics has therefore ignited the interest of researchers and clinicians alike and has given us new insights into the cause and course of the disease. This article is an overview of the extensive work published, predominantly in the last two decades, on the biomechanical aspect of keratoconus.

METHODS

Published articles on corneal biomechanics in the specific context of keratoconus were reviewed, based on an electronic search using PubMed, Elsevier, and Science Direct. The search terms used included "Corneal Biomechanics," "Mechanical properties of the cornea," "Corneal ultrastructure," "Corneal Collagen," and "Keratoconus". Articles pertaining to refractive surgery, keratoplasty, collagen crosslinking, or intrastromal rings were excluded.

RESULTS

The electronic search revealed more than 500 articles, from which 80 were chosen for this article.

CONCLUSIONS

The structural and organizational pattern of the corneal stroma determines its mechanical properties and are responsible for the maintenance of the normal shape and function of the cornea. Changes in the ultrastructure are responsible for the biomechanical instability that leads to corneal ectasia. As non-invasive methods for evaluating corneal biomechanics in vivo evolve, our ability to diagnose subclinical keratoconus will improve, allowing identification of patients at risk to develop ectasia and to allow early treatment to arrest progression of the disease.

Clinical Ocular Biomechanics: Where Are We after 20 Years of Progress?

Purpose: Ocular biomechanics is an assessment of the response of the structures of the eye to forces that may lead to disease development and progression, or influence the response to surgical intervention. The goals of this review are (1) to introduce basic biomechanical principles and terminology, (2) to provide perspective on the progress made in the clinical study and assessment of ocular biomechanics, and (3) to highlight critical studies conducted in keratoconus, laser refractive surgery, and glaucoma in order to aid interpretation of biomechanical parameters in the laboratory and in the clinic.Methods: A literature review was first conducted of basic biomechanical studies related to ocular tissue. The subsequent review of ocular biomechanical studies was limited to those focusing on keratoconus, laser refractive surgery, or glaucoma using the only two commercially available devices that allow rapid assessment of biomechanical response in the clinic.Results: Foundational studies on ocular biomechanics used a combination of computer modeling and destructive forces on ex-vivo tissues. The knowledge gained from these studies could not be directly translated to clinical research and practice until the introduction of non-contact tonometers that quantified the deformation response of the cornea to an air puff, which represents a non-destructive, clinically appropriate load. The corneal response includes a contribution from the sclera which may limit corneal deformation. Two commercial devices are available, the Ocular Response Analyzer which produces viscoelastic parameters with a customized load for each eye, and the Corvis ST which produces elastic parameters with a consistent load for every eye. Neither device produces the classic biomechanical properties reported in basic studies, but rather biomechanical deformation response parameters which require careful interpretation.Conclusions: Research using clinical tools has enriched our understanding of how ocular disease alters ocular biomechanics, as well as how ocular biomechanics may influence the pathophysiology of ocular disease and response to surgical intervention.

Motion-tracking Brillouin microscopy for in-vivo corneal biomechanics mapping

Corneal biomechanics play a critical role in maintaining corneal shape and thereby directly influence visual acuity. However, direct corneal biomechanical measurement in-vivo with sufficient accuracy and a high spatial resolution remains an open need. Here, we developed a three-dimensional (3D) motion-tracking Brillouin microscope for in-vivo corneal biomechanics mapping. The axial tracking utilized optical coherence tomography, which provided a tracking accuracy better than 3 µm. Meanwhile, 10 µm lateral tracking was achieved by tracking pupils with digital image processing. The 3D tracking enabled reconstruction of depth-dependent Brillouin distribution with a high spatial resolution. This superior technical performance enabled the capture of high-quality mechanical mapping in vivo even while the subject was breathing normally. Importantly, we improved Brillouin spectral measurements to achieve relative accuracy better than 0.07% verified by rubidium absorption frequencies, with 0.12% stability over 2000 seconds. These specifications finally yield the Brillouin measurement sensitivity that is required to detect ophthalmology-relevant corneal biomechanical properties.

Determining the center of a keratoconus: Comparison of different tomographic parameters and impact of disease severity

Purpose

There exists remarkable variation in definitions for the location of the center of a keratoconus. The objective of this study was to analyze deviations between locations obtained by different tomographic maps for that purpose. Furthermore, it was investigated whether these deviations are influenced by disease severity.

Methods

In 162 eyes with keratoconus, corneal tomographic maps derived by Scheimpflug technology were retrospectively analyzed to determine the cone location with 5 different methods: maximum axial curvature of the front surface (Kmax), maximum tangential curvature of the front surface (tKmax), minimum pachymetry (Pachymin), maximum elevation of the front surface (ELEF), and maximum elevation of the back surface (ELEB). Distances between the locations were calculated and tested for a correlation with keratoconus severity and distance between cone and corneal vertex.

Results

Cone locations derived from the curvature maps (Kmax, tKmax) showed the lowest agreement with the locations determined by pachymetry or elevation maps. The largest distances were found between Kmax and Pachymin [Median and Interquartile range: 1.19 mm (0.87, 1.60)], Kmax and ELEB [1.12 mm (0.79, 1.41)], and Kmax and ELEF [0.97 mm (0.64, 1.27)]. Low distances (<0.5 mm) were calculated between ELEB and ELEF, and ELEB and Pachymin. All of the calculated distances between the locations showed a significant negative correlation with keratoconus severity and most of them increased significantly with a more peripheral position of the cone (p < 0.05).

Conclusions

There was low consistency between different methods for describing the location of a keratoconus. Curvature-based determinations of the cone center (Kmax, tKmax) showed the highest deviations and should not be used for that purpose. However, the discrepancies between different cone location methods diminished with increasing disease severity and more central position of the cone.

High-frequency ultrasound detects biomechanical weakening in keratoconus with lower stiffness at higher grade

In vivo biomechanical characterization of the cornea remains a challenge. We have developed a high-frequency ultrasound elastography method, the ocular pulse elastography (OPE), to measure corneal axial displacement (CAD) induced by the ocular pulse. Here we compared CAD and a stiffness index derived from CAD between keratoconus patients and normal controls. We also explored the trend of these parameters with keratoconus grade. Twenty normal subjects and twenty keratoconus patients were recruited in this study. Corneal topography, tomography, intraocular pressure (IOP) and ocular pulse amplitude (OPA) were obtained in each measured eye. The cornea’s heartbeat-induced cyclic axial displacements were measured by high-frequency (50 MHz) ultrasound. A corneal stiffness index (CSI) was derived from CAD normalized against OPA. CAD and CSI were compared between normal and keratoconus groups, and across keratoconus grades. Keratoconus corneas had significantly greater CAD and lower CSI than normal controls (p’s<0.01). Both parameters correlated strongly with grade, in which CAD increased significantly (p = 0.002) and CSI decreased significantly (p = 0.011) with grade. These results suggested a biomechanical weakening in keratoconus which worsens at higher disease severity. This study also demonstrated the ability of high-frequency ultrasound elastography to provide a safe, quick, and accurate evaluation of the cornea’s biomechanical condition in vivo. The OPE-measured biomechanical metrics, when integrated with existing diagnostic criteria, may aid the decision-making in the early and definitive diagnosis and staging of keratoconus.

In vivo Assessment of Localised Corneal Biomechanical Deterioration With Keratoconus Progression

Purpose: To evaluate the regional corneal biomechanical deterioration with keratoconus (KC) progression as measured by the Stress-Strain Index (SSI) maps.

Methods: The preoperative examinations of 29 progressive KC cases that were submitted to corneal cross-linking (CXL) were evaluated. The examinations included the tomography and the SSI measured by the Pentacam HR and the Corvis ST (Oculus, Wetzlar, Germany), respectively. The results were recorded twice, the latter of which was at the last visit before the CXL procedure. The patient-specific SSI maps were built, using data at each examination, based on finite element modelling and employing inverse analysis to represent the regional variation of biomechanical stiffness across the cornea.

Results: All cases presented significant shape progression (above the 95% CI of repeatability) in anterior and posterior curvatures and minimum thickness. The overall corneal stiffness as measured by the SSI within the central 8 mm-diameter area underwent slight but significant reductions from the first to the last examination (−0.02 ± 0.02, range: −0.09 to 0, p < 0.001). In all 29 cases, the reduction in stiffness was localised and concentred in the area inside the keratoconus cone. The SSI values inside the cone were significantly lower in the last examination (by 0.15 ± 0.09, range: −0.42 to −0.01, p < 0.001), while the SSI outside the cone presented minimal, non-significant variations (0 ± 0.01, range: −0.04 to 0.01, p = 0.999).

Conclusion: It has been observed through the SSI maps that the regional deterioration in stiffness was concerted inside the area of pathology, while only mild non-significant alterations were observed outside the area of pathology.

Torsional wave elastography to assess the mechanical properties of the cornea

Corneal mechanical changes are believed to occur before any visible structural alterations observed during routine clinical evaluation. This study proposed developing an elastography technique based on torsional waves (TWE) adapted to the specificities of the cornea. By measuring the displacements in the propagation plane perpendicular to the axis of the emitter, the effect of guided waves in plate-like media was proven negligible. Ex vivo experiments were carried out on porcine corneal samples considering a group of control and one group of alkali burn treatment (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {NH}_\text {4}$$\end{document}NH4OH) that modified the mechanical properties. Phase speed was recovered as a function of intraocular pressure (IOP), and a Kelvin-Voigt rheological model was fitted to the dispersion curves to estimate viscoelastic parameters. A comparison with uniaxial tensile testing with thin-walled assumptions was also performed. Both shear elasticity and viscosity correlated positively with IOP, being the elasticity lower and the viscosity higher for the treated group. The viscoelastic parameters ranged from 21.33 to 63.17 kPa, and from 2.82 to 5.30 Pa s, for shear elasticity and viscosity, respectively. As far as the authors know, no other investigations have studied this mechanical plane under low strain ratios, typical of dynamic elastography in corneal tissue. TWE reflected mechanical properties changes after treatment, showing a high potential for clinical diagnosis due to its rapid performance time and paving the way for future in vivo studies.

Biomechanical Evaluation of Topographically and Tomographically Normal Fellow Eyes of Patients With Keratoconus

PURPOSE

To determine the effectiveness of parameters and indices based on biomechanical measures at discriminating fellow eyes with topographically and tomographically normal corneas in patients with keratoconus from normal control corneas.

METHODS

The study included 47 keratoconus suspect eyes, defined as the topographically and tomographically normal fellow eyes of patients with frank keratoconus in the other eye. Eyes were imaged using the Pentacam HR and Corvis ST (both Oculus Optikgeräte GmbH). Fellow eyes were then categorized as topographically/tomographically normal fellow eyes (TNF) and topographically/tomographically borderline fellow eyes (TBF). The ability of each of the Corvis Biomechanical Index (CBI), Tomographic and Biomechanical Index (TBI), stiffness parameter at applanation 1 (SP-A1), and stress-strain index (SSI) at discriminating between normal controls and keratoconus suspects was assessed.

RESULTS

The TBI had the best discriminative ability with the greatest area under the receiver operating characteristic (AUROC) curve value of 0.946 for normal controls versus TBF eyes, and 0.824 for normal controls versus TNF eyes. Compared to the TBI AUROC curves, SP-A1 and CBI had AUROC curve values of 0.833 (P = .09) and 0.822 (P = .01) for normal controls versus TBF eyes, respectively, and AUROC curve values of 0.822 (P = .96) and 0.550 (P = .0002) for normal controls versus TNF eyes, respectively. The TBI had the best positive predictive value for TNF and TBF eyes, followed by CBI and SP-A1.

CONCLUSIONS

The TBI and the purely biomechanical parameter SP-A1 were of moderate utility in distinguishing between normal and keratoconus suspect eyes. In the absence of topographic/tomographic evidence of keratectasia, an independently abnormal biomechanical parameter may suggest an increased risk of ectasia. [J Refract Surg. 2022;38(5):318-325.].

Localization-assisted stimulated Brillouin scattering spectroscopy

Brillouin spectroscopy has emerged as a promising modality to noninvasively probe the mechanical properties of biologically relevant materials. Stimulated Brillouin scattering (SBS) has the potential to improve measurement speed and resolution by exploiting a resonant amplification of the scattered signal, yet current SBS spectrometers have not provided significant improvements due to fundamental and practical limitations of illumination and detection parameters. To overcome this challenge, here we derive a signal localization theory for the Brillouin spectral domain and accordingly design an SBS spectrometer with much improved performances compared to state-of-the-art systems. We present experimental and simulated data validating our theory, which result in a tenfold improvement in acquisition speed, or an order of magnitude improved spectral precision, for SBS spectral measurements when properly optimizing the SBS photon detection architecture.

Determine Corneal Biomechanical Parameters by Finite Element Simulation and Parametric Analysis Based on ORA Measurements

Purpose: The Ocular Response Analyzer (ORA) is one of the most commonly used devices to measure corneal biomechanics in vivo. Until now, the relationship between the output parameters and corneal typical biomechanical parameters was not clear. Hence, we defined the output parameters of ORA as ORA output parameters. This study aims to propose a method to determine corneal biomechanical parameters based on ORA measurements by finite element simulation and parametric analysis.

Methods: Finite element analysis was used to simulate the mechanics process of ORA measurements with different intraocular pressure (IOP), corneal geometrical parameters and corneal biomechanical parameters. A simplified geometrical optics model was built to simulate the optical process of the measurements to extract ORA output parameters. After that, 70% of the simulated data was used to establish the quantitative relationship between corneal biomechanical parameters and ORA output parameters by parametric analysis and 30% of the simulated data was used to validate the established model. Besides, ten normal subjects were included to evaluate the normal range of corneal biomechanical parameters calculated from ORA.

Results: The quantitative relationship between corneal biomechanical parameters and ORA output parameters is established by combining parametric analysis with finite element simulation. The elastic modulus (E) and relaxation limit (G_∞) of the ten normal subjects were 0.65 ± 0.07 MPa and 0.26 ± 0.15, respectively.

Conclusions: A method was proposed to determine corneal biomechanical parameters based on the results of ORA measurements. The magnitude of the corneal biomechanical parameters calculated according to our method was reasonable.

Multimodal Heartbeat and Compression Optical Coherence Elastography for Mapping Corneal Biomechanics

The biomechanical properties of the cornea have a profound influence on the health, structural integrity, and function of the eye. Understanding these properties may be critical for diagnosis and identifying disease pathogenesis. This work demonstrates how two different elastography techniques can be combined for a multimodal approach to measuring corneal biomechanical properties. Heartbeat optical coherence elastography (Hb-OCE) and compression OCE were performed simultaneously to measure the stiffness of the cornea in an in vivo rabbit model. Measurements were further performed after collagen crosslinking to demonstrate how the combined technique can be used to measure changes in corneal stiffness and map mechanical contrast. The results of this work further suggest that measurements from Hb-OCE and compression OCE are comparable, meaning that Hb-OCE and compression OCE may be used interchangeably despite distinct differences in both techniques.

Measuring mechanical anisotropy of the cornea with Brillouin microscopy

Load-bearing tissues are typically fortified by networks of protein fibers, often with preferential orientations. This fiber structure imparts the tissues with direction-dependent mechanical properties optimized to support specific external loads. To accurately model and predict tissues’ mechanical response, it is essential to characterize the anisotropy on a microstructural scale. Previously, it has been difficult to measure the mechanical properties of intact tissues noninvasively. Here, we use Brillouin optical microscopy to visualize and quantify the anisotropic mechanical properties of corneal tissues at different length scales. We derive the stiffness tensor for a lamellar network of collagen fibrils and use angle-resolved Brillouin measurements to determine the longitudinal stiffness coefficients (longitudinal moduli) describing the ex vivo porcine cornea as a transverse isotropic material. Lastly, we observe significant mechanical anisotropy of the human cornea in vivo, highlighting the potential for clinical applications of off-axis Brillouin microscopy.

Here, Brillouin optical microscopy noninvasively visualizes microscale anisotropy of the porcine cornea owing to its lamellar fiber structure and quantifies the longitudinal moduli of the bulk tissue. Anisotropy is also detected in angle-resolved measurement of the human cornea in vivo.

Universal architecture of corneal segmental tomography biomarkers for artificial intelligence-driven diagnosis of early keratoconus

AIMS

To develop a comprehensive three-dimensional analyses of segmental tomography (placido and optical coherence tomography) using artificial intelligence (AI).

METHODS

Preoperative imaging data (MS-39, CSO, Italy) of refractive surgery patients with stable outcomes and diagnosed with asymmetric or bilateral keratoconus (KC) were used. The curvature, wavefront aberrations and thickness distributions were analysed with Zernike polynomials (ZP) and a random forest (RF) AI model. For training and cross-validation, there were groups of healthy (n=527), very asymmetric ectasia (VAE; n=144) and KC (n=454). The VAE eyes were the fellow eyes of KC patients but no further manual segregation of these eyes into subclinical or forme-fruste was performed.

RESULTS

The AI achieved an excellent area under the curve (0.994), accuracy (95.6%), recall (98.5%) and precision (92.7%) for the healthy eyes. For the KC eyes, the same were 0.997, 99.1%, 98.7% and 99.1%, respectively. For the VAE eyes, the same were 0.976, 95.5%, 71.5% and 91.2%, respectively. Interestingly, the AI reclassified 36 (subclinical) of the VAE eyes as healthy though these eyes were distinct from healthy eyes. Most of the remaining VAE (n=104; forme fruste) eyes retained their classification, and were distinct from both KC and healthy eyes. Further, the posterior surface features were not among the highest ranked variables by the AI model.

CONCLUSIONS

A universal architecture of combining segmental tomography with ZP and AI was developed. It achieved an excellent classification of healthy and KC eyes. The AI efficiently classified the VAE eyes as 'subclinical' and 'forme-fruste'.

Dissecting Organismal Morphogenesis by Bridging Genetics and Biophysics

Multicellular organisms develop complex shapes from much simpler, single-celled zygotes through a process commonly called morphogenesis. Morphogenesis involves an interplay between several factors, ranging from the gene regulatory networks determining cell fate and differentiation to the mechanical processes underlying cell and tissue shape changes. Thus, the study of morphogenesis has historically been based on multidisciplinary approaches at the interface of biology with physics and mathematics. Recent technological advances have further improved our ability to study morphogenesis by bridging the gap between the genetic and biophysical factors through the development of new tools for visualizing, analyzing, and perturbing these factors and their biochemical intermediaries. Here, we review how a combination of genetic, microscopic, biophysical, and biochemical approaches has aided our attempts to understand morphogenesis and discuss potential approaches that may be beneficial to such an inquiry in the future. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Comparative Analysis of Safety and Efficacy of Topography-Guided Customized Cross-linking and Standard Cross-linking in the Treatment of Progressive Keratoconus

PURPOSE

To compare the safety and efficacy of topography-guided customized corneal cross-linking (PiXL) with standard cross-linking (CXL) for the treatment of progressive keratoconus.

METHODS

In a prospective interventional analysis, the eyes of patients with progressive keratoconus underwent standard cross-linking (homogenous 9-mm ultraviolet-A irradiation of 9 mW/cm2 delivering a total fluence of 5.4J/cm2) versus topography-guided customized cross-linking (30 mW/cm2 pulsed irradiance with a total fluence ranging from 5.4 to 15 J/cm2, in concentric circles centered on the posterior float maximum). The following parameters were analyzed at the preoperative, 1-month, 6-month, and 1-year postoperative visits: corrected spectacle distance visual acuity, manifest refraction including sphere, cylinder and mean refractive spherical equivalent, corneal tomography, higher order aberration profile, and endothelial cell count. Anterior segment optical coherence tomography evaluation was performed at the 1-month postoperative visit to assess the depth of the demarcation line.

RESULTS

Sixty-four eyes of 45 patients (32 eyes in each group) were included. There was a significant reduction in maximum keratometry and IS asymmetry in the PiXL group at both the 6-month and 1-year postoperative visits (P = 0.001 and 0.06). Corrected spectacle distance visual acuity improved significantly in the PiXL (0.05 ± 0.08 logarithm of the minimum angle of resolution, P = 0.02) versus the standard CXL (0.01 ± 0.025 logarithm of the minimum angle of resolution, P = 0.26) group. A greater depth of a stromal demarcation line was observed in the customized CXL group (P = 0.02). No significant complications were noted in either cohort.

CONCLUSIONS

PiXL demonstrated similar safety with significantly greater keratometry flattening and corneal regularization vis-à-vis standard CXL. This resulted in significant improvement of spectacle corrected visual acuity for eyes with mild-to-moderate keratoconus.

Regional variation of corneal stromal deformation measured by high-frequency ultrasound elastography

The cornea's mechanical response to intraocular pressure elevations may alter in ectatic diseases such as keratoconus. Regional variations of mechanical deformation in normal and keratoconus eyes during intraocular pressure elevation have not been well-characterized. We applied a high-frequency ultrasound elastography technique to characterize the regional deformation of normal and keratoconus human corneas through the full thickness of corneal stroma. A cross-section centered at the corneal apex in 11 normal and 2 keratoconus human donor eyes was imaged with high-frequency ultrasound during whole globe inflation from 5 to 30 mmHg. An ultrasound speckle tracking algorithm was used to compute local tissue displacements. Radial, tangential, and shear strains were mapped across the imaged cross-section. Strains in the central (1 mm surrounding apex) and paracentral (1 to 4 mm from apex) regions were analyzed in both normal and keratoconus eyes. Additional regional analysis was performed in the eye with severe keratoconus presenting significant thinning and scarring. Our results showed that in normal corneas, the central region had significantly smaller tangential stretch than the paracentral region, and that within the central region, the magnitudes of radial and shear strains were significantly larger than that of tangential strain. The eye with mild keratoconus had similar shear strain but substantially larger radial strains than normal corneas, while the eye with severe keratoconus had similar overall strains as in normal eyes but marked regional heterogeneity and large strains in the cone region. These findings suggested regional variation of mechanical responses to intraocular pressure elevation in both normal and keratoconus corneas, and keratoconus appeared to be associated with mechanical weakening in the cone region, especially in resisting radial compression. Comprehensive characterization of radial, tangential, and shear strains through corneal stroma may provide new insights to understand the biomechanical alterations in keratoconus.

Brillouin Spectroscopy: From Biomedical Research to New Generation Pathology Diagnosis

Brillouin spectroscopy has recently gained considerable interest within the biomedical field as an innovative tool to study mechanical properties in biology. The Brillouin effect is based on the inelastic scattering of photons caused by their interaction with thermodynamically driven acoustic modes or phonons and it is highly dependent on the material's elasticity. Therefore, Brillouin is a contactless, label-free optic approach to elastic and viscoelastic analysis that has enabled unprecedented analysis of ex vivo and in vivo mechanical behavior of several tissues with a micrometric resolution, paving the way to a promising future in clinical diagnosis. Here, we comprehensively review the different studies of this fast-moving field that have been performed up to date to provide a quick guide of the current literature. In addition, we offer a general view of Brillouin's biomedical potential to encourage its further development to reach its implementation as a feasible, cost-effective pathology diagnostic tool.

Update on corneal collagen crosslinking for ectasia

PURPOSE OF REVIEW

Corneal collagen crosslinking (CXL) is a minimally invasive treatment that can stabilize corneal ectatic disorders including keratoconus, pellucid marginal degeneration, or postrefractive surgery ectasia. The benefits of CXL have been well documented. New research is focused on modifying current treatment protocols with the goals of maximizing corneal stability while also shortening overall procedure time.

RECENT FINDINGS

Accelerated CXL protocols have the goal of delivering the same ultraviolet A intensity as conventional protocols, but over a shorter time period. Accelerated protocols have shown success to date, but there are concerns for long-term corneal stability. Pulsed protocols may increase the long-term efficacy of the accelerated designs. In addition, transepithelial crosslinking protocols have been designed with the goal of reducing postoperative pain and lower the risk of infectious complications of epithelial-off conventional protocols.

SUMMARY

Newer CXL protocols attempt to make the procedure safer and more effective. Current research is promising, but long-term studies are essential to understand how the new protocols may affect corneal stability.

Corneal biomechanical parameters in keratoconus eyes with abnormal elevation on the back corneal surface only versus both back and front surfaces

Corneal biomechanical parameters were compared in 100 keratoconus eyes with abnormal elevation on the back corneal surface only (group 1), versus both the back and front surfaces (group 2). Scheimpflug tomography with Pentacam HR, corneal biomechanical assessments using Corvis ST and Ocular Response Analyzer (ORA) and corneal epithelium thickness maps using anterior segment optical coherence tomography were assessed. There were no significant differences in the IOP measured using Corvis ST and ORA, age or sex between the two groups. Statistically significant differences were found in all corneal shape parameters and all new parameters of Corvis ST: corneal stiffness parameter at first applanation (SP-A1), integrated inverse radius (IR) and deformation amplitude ratio (DAR)) between groups (p < 0.001). The classic parameters of ORA including corneal hysteresis (CH) and corneal resistance factor (CRF) were about 1.00 mmHg higher in group 1 (p < 0.001). In conclusion, keratoconus eyes with abnormal elevation limited to the back corneal surface have lower grade, stiffer corneal biomechanical parameters and less asymmetric shape. This is consistent with progressive biomechanical weakening from the first detectable back surface elevation to manifestation on the front surface as the severity overwhelms the ability of the epithelium to compensate.

Analysis of OPD-Scan and Pentacam Parameters for Early Keratoconus Detection

BACKGROUND

PURPOSE: To evaluate Pentacam and OPD-Scan parameters in the early detection of keratoconus.

DESIGN

Retrospective case-control study.

METHODS

Case group included 50 clinically unaffected fellow eyes diagnosed with asymmetric keratoconus showing subtle qualitative changes at the 0.5-D sensitivity OPD-Scan scale, as well as normal anterior and back elevation difference map at Belin/Ambrósio enhanced ectasia display (BAD) at the Pentacam. Control group included 172 normal eyes that underwent Lasik surgery and presented no complications throughout the 2-year follow-up period. OPD-Scan and Pentacam parameters were compared, calculating sensitivity, specificity, and area under the receiver operating characteristic curve (AUC). A multivariate analysis was performed using Pentacam or OPD-Scan variables, and a model using variables of both devices.

RESULTS

Pentacam variables with AUC ≥0.8 were keratoconus index (0.85), index of height decentration (0.81), and overall deviation at BAD (0.8). OPD-Scan variables with AUC ≥0.8 were keratoconus prediction index (0.83), surface asymmetry index (0.83), and total of higher-order trefoil aberration (0.8). In the multivariate analysis, the AUC was 0.85 in the case of OPD-Scan whereas it was 0.89 in the case of Pentacam. When combining all variables from the 2 devices, the AUC was 0.93, with a sensitivity of 82% and a specificity of 94%.

CONCLUSIONS

Several parameters of OPD-Scan and Pentacam can be useful to differentiate cases from normal control eyes, demonstrating even better results when combining parameters of both devices. Anterior corneal indexes were the most important parameters to discriminate both groups.

Corneal Crosslinking for Progressive Keratoconus and Corneal Ectasia: Summary of US Multicenter and Subgroup Clinical Trials

Purpose

The multicenter studies reviewed were designed to evaluate the safety and efficacy of corneal collagen crosslinking (CXL) for the treatment of progressive keratoconus and corneal ectasia after laser refractive surgery. The results of these studies led to approval by the United States Food and Drug Agency for both conditions in 2016. This paper reviews these studies, as well as single-center substudies investigating other aspects of crosslinking outcomes.

Methods

As part of prospective, randomized, controlled clinical trials, the treatment group received standard CXL, and the sham control group received only riboflavin ophthalmic solution. The primary efficacy criterion was maximum keratometry (K_max) 1 year after CXL. Secondary outcomes were corrected distance visual acuity (CDVA) and uncorrected distance visual acuity (UDVA). Safety and adverse events were analyzed. In single-center substudies, corneal topography, ocular aberrations, corneal haze measurements, corneal thickness, corneal biomechanics, subjective visual function, and outcomes predictors were also investigated. This paper presents a general review of the design and outcomes of crosslinking in these studies.

Results

In the crosslinking treatment group, K_max flattened by 1.6 diopters (D) and 0.7 D in eyes with keratoconus and ectasia, respectively. In both studies, there was continued progression in the control group. The CDVA improved by an average of 5.7 logMAR letters (LL) in the keratoconus treatment group and by 5.0 LL in the ectasia group. In both studies, corneal haze was the most frequently reported crosslinking-related adverse finding. This was most prominent at 1 month and generally returned to baseline between 3 and 12 months. In general, corneal topography, ocular aberrations, and subjective visual function improved after crosslinking.

Conclusions

In the US multicenter trials, CXL was shown to be safe and effective in stabilizing K_max, CDVA, and UDVA in eyes with progressive keratoconus or corneal ectasia.

Translational Relevance

Corneal crosslinking was originally developed in the laboratory at the University of Dresden in the late 1990s. The combination of ultraviolet-A light and riboflavin was found to be the most effective of a number of different modalities tested to increase the biomechanical strength of the cornea. The clinical study design for the US multicenter clinical trials of crosslinking demonstrated the safety and effectiveness of this technique for treatment of progressive keratoconus and corneal ectasia, bringing this important advancement to patients in the United States.

Label-free histological imaging of tissues using Brillouin light scattering contrast

Brillouin light scattering offers a unique label-free approach to measure biomechanical properties non-invasively. While this technique is used in biomechanical analysis of cells and tissues, its potential for visualizing structural features of tissues based on the biomechanical contrast has not been much exploited. Here, we present high-resolution Brillouin microscopy images of four basic tissue types: muscular, connective, epithelial, and nervous tissues. The Brillouin contrast distinguishes between muscle fiber cells and endomysium in skeletal muscle and reveals chondrocytes along with spatially varying stiffness of the extracellular matrix in articular cartilage. The hydration-sensitive contrast can visualize the stratum corneum, epidermis, and dermis in the skin epithelium. In brain tissues, the Brillouin images show the mechanical heterogeneity across the cortex and deeper regions. This work demonstrates the versatility of using the Brillouin shift as histological contrast for examining intact tissue substructures via longitudinal modulus without the need for laborious tissue processing steps.

Advances in Imaging Technology of Anterior Segment of the Eye

Advances in imaging technology and computer science have allowed the development of newer assessment of the anterior segment, including Corvis ST, Brillouin microscopy, ultrahigh-resolution optical coherence tomography, and artificial intelligence. They enable accurate and precise assessment of structural and biomechanical alterations associated with anterior segment disorders. This review will focus on these 4 new techniques, and a brief overview of these modalities will be introduced. The authors will also discuss the recent advances in research regarding these techniques and potential application of these techniques in clinical practice. Many studies on these modalities have reported promising results, indicating the potential for more detailed comprehensive understanding of the anterior segment tissues.

Corneal biomechanics: Measurement and structural correlations

The characterization of corneal biomechanical properties has important implications for the management of ocular disease and prediction of surgical responses. Corneal refractive surgery outcomes, progression or stabilization of ectatic disease, and intraocular pressure determination are just examples of the many key clinical problems that depend highly upon corneal biomechanical characteristics. However, to date there is no gold standard measurement technique. Since the advent of a 1-dimensional (1D) air-puff based technique for measuring the corneal surface response in 2005, advances in clinical imaging technology have yielded increasingly sophisticated approaches to characterizing the biomechanical properties of the cornea. Novel analyses of 1D responses are expanding the clinical utility of commercially-available air-puff-based instruments, and other imaging modalities-including optical coherence elastography (OCE), Brillouin microscopy and phase-decorrelation ocular coherence tomography (PhD-OCT)-offer new opportunities for probing local biomechanical behavior in 3-dimensional space and drawing new inferences about the relationships between corneal structure, mechanical behavior, and corneal refractive function. These advances are likely to drive greater clinical adoption of in vivo biomechanical analysis and to support more personalized medical and surgical decision-making.

Mapping mechanical properties of biological materials via an add-on Brillouin module to confocal microscopes

Several techniques have been developed over the past few decades to assess the mechanical properties of biological samples, which has fueled a rapid growth in the fields of biophysics, bioengineering, and mechanobiology. In this context, Brillouin optical spectroscopy has long been known as an intriguing modality for noncontact material characterization. However, limited by speed and sample damage, it had not translated into a viable imaging modality for biomedically relevant materials. Recently, based on a novel spectroscopy strategy that substantially improves the speed of Brillouin measurement, confocal Brillouin microscopy has emerged as a unique complementary tool to traditional methods as it allows noncontact, nonperturbative, label-free measurements of material mechanical properties. The feasibility and potential of this innovative technique at both the cell and tissue level have been extensively demonstrated over the past decade. As Brillouin technology is rapidly recognized, a standard approach for building and operating Brillouin microscopes is required to facilitate the widespread adoption of this technology. In this protocol, we aim to establish a robust approach for instrumentation, and data acquisition and analysis. By carefully following this protocol, we expect that a Brillouin instrument can be built in 5-9 days by a person with basic optics knowledge and alignment experience; the data acquisition as well as postprocessing can be accomplished within 2-8 h.

Heartbeat-Induced Corneal Axial Displacement and Strain Measured by High Frequency Ultrasound Elastography in Human Volunteers

Purpose

The purpose of this study was to establish in vivo data acquisition and processing protocols for repeatable measurements of heartbeat-induced corneal displacements and strains in human eyes, using a high-frequency ultrasound elastography method, termed ocular pulse elastography (OPE).

Methods

Twenty-four volunteers with no known ocular diseases were recruited for this study. Intraocular pressure (IOP) and ocular pulse amplitude (OPA) were measured using a PASCAL Dynamic Contour Tonometer (DCT). An in vivo OPE protocol was developed to measure heartbeat-induced corneal displacements. Videos of the central 5.7 mm of the cornea were acquired using a 50-MHz ultrasound probe at 128 frames per second. The radiofrequency data of 1000 frames were analyzed using an ultrasound speckle tracking algorithm to calculate corneal displacements and quantify spectral and temporal characteristics. The intrasession and intersession repeatability of OPE- and DCT-measured parameters were also analyzed.

Results

The in vivo OPE protocol and setup were successful in tracking heartbeat-induced corneal motion using high-frequency ultrasound. Corneal axial displacements showed a strong cardiac rhythm, with good intrasession and intersession repeatability, and high interocular symmetry. Corneal strain was calculated in two eyes of two subjects, showing substantially different responses.

Conclusions

We demonstrated the feasibility of high-frequency ultrasound elastography for noninvasive in vivo measurement of the cornea's biomechanical responses to the intrinsic ocular pulse. The high intrasession and intersession repeatability suggested a robust implementation of this technique to the in vivo setting.

Translational Relevance

OPE may offer a useful tool for clinical biomechanical evaluation of the cornea by quantifying its response to the intrinsic pulsation.

Optimization of surgical parameters based on patient-specific models: Application to arcuate keratotomy

PURPOSE

To determine surgical parameters for arcuate keratotomy by simulating the intervention with a patient-specific model.

SETTING

University Eye Clinic Salzburg, Paracelsus Medical University, Austria, and Institute for Surgical Technology and Biomechanics, University of Bern, Switzerland.

DESIGN

Computational modeling study.

METHODS

A new approach to plan arcuate keratotomy based on personalized finite element simulations was developed. Using this numeric tool, an optimization algorithm was implemented to determine the incision parameters that best met the surgeon's requirements while preserving the orientation of the astigmatism. Virtual surgeries were performed on patients to compare the performance of the simulation-based approach with results based on the Lindstrom and Donnenfeld nomograms and with intrastromal interventions.

RESULTS

Retrospective data on 28 patients showed that personalized simulation reproduced the surgically induced change in astigmatism (Pearson correlation = 0.8). Patient-specific simulation was used to examine strategies for arcuate interventions on 621 corneal topographies. The Lindstrom nomogram resulted in low postoperative astigmatism (mean 0.03 diopter [D] ± 0.3 [SD]) but frequent overcorrections (20%). The Donnenfeld nomogram and intrastromal incisions resulted in a small amount of overcorrection (1.5%) but a wider spread in astigmatism (mean 0.63 ± 0.35 D and 0.48 ± 0.50 D, respectively). In contrast, the new numeric parameter optimization approach led to postoperative astigmatism values (mean 0.40 ± 0.08 D, 0.20 ± 0.08 D, and 0.04 ± 0.13 D) that closely matched the target astigmatism (0.40 D, 0.20 D, and 0.00 D), respectively, while keeping the number of overcorrections low (<1.5%).

CONCLUSION

Using numeric modeling to optimize surgical parameters for arcuate keratotomy led to more reliable postoperative astigmatism, limiting the risk for overcorrection.

Characterization of retinal biomechanical properties using Brillouin microscopy

SIGNIFICANCE

The retina is critical for vision, and several diseases may alter its biomechanical properties. However, assessing the biomechanical properties of the retina nondestructively is a challenge due to its fragile nature and location within the eye globe. Advancements in Brillouin spectroscopy have provided the means for nondestructive investigations of retina biomechanical properties.

AIM

We assessed the biomechanical properties of mouse retinas using Brillouin microscopy noninvasively and showed the potential of Brillouin microscopy to differentiate the type and layers of retinas based on stiffness.

APPROACH

We used Brillouin microscopy to quantify stiffness of fresh and paraformaldehyde (PFA)-fixed retinas. As further proof-of-concept, we demonstrated a change in the stiffness of a retina with N-methyl-D-aspartate (NMDA)-induced damage, compared to an undamaged sample.

RESULTS

We found that the retina layers with higher cell body density had higher Brillouin modulus compared to less cell-dense layers. We have also demonstrated that PFA-fixed retina samples were stiffer compared with fresh samples. Further, NMDA-induced neurotoxicity leads to retinal ganglion cell (RGC) death and reactive gliosis, increasing the stiffness of the RGC layer.

CONCLUSION

Brillouin microscopy can be used to characterize the stiffness distribution of the layers of the retina and can be used to differentiate tissue at different conditions based on biomechanical properties.

Characterization of non-linear mechanical behavior of the cornea

The objective of this study was to evaluate which hyperelastic model could best describe the non-linear mechanical behavior of the cornea, in order to characterize the capability of the non-linear model parameters to discriminate structural changes in a damaged cornea. Porcine corneas were used, establishing two different groups: control (non-treated) and NaOH-treated (damaged) corneas (n = 8). NaOH causes a chemical burn to the corneal tissue, simulating a disease associated to structural damage of the stromal layer. Quasi-static uniaxial tensile tests were performed in nasal-temporal direction immediately after preparing corneal strips from the two groups. Three non-linear hyperelastic models (i.e. Hamilton-Zabolotskaya model, Ogden model and Mooney-Rivlin model) were fitted to the stress–strain curves obtained in the tensile tests and statistically compared. The corneas from the two groups showed a non-linear mechanical behavior that was best described by the Hamilton-Zabolotskaya model, obtaining the highest coefficient of determination (R² > 0.95). Moreover, Hamilton-Zabolotskaya model showed the highest discriminative capability of the non-linear model parameter (Parameter A) for the tissue structural changes between the two sample groups (p = 0.0005). The present work determines the best hyperelastic model with the highest discriminative capability in description of the non-linear mechanical behavior of the cornea.

Detecting Mechanical Anisotropy of the Cornea Using Brillouin Microscopy

Purpose

The purpose of this study was to detect the mechanical anisotropy of the cornea using Brillouin microscopy along different perturbation directions.

Methods

Brillouin frequency shift of both whole globes (n = 10) and cornea punches (n = 10) were measured at different angles to the incident laser, thereby probing corneal longitudinal modulus of elasticity along different directions. Frequency shift of virgin (n = 26) versus cross-linked corneas (n = 15) over a large range of hydration conditions were compared in order to differentiate the contributions to Brillouin shift due to hydration from those due to stromal tissue.

Results

We detected mechanical anisotropy of corneas, with an average frequency shift increase of 53 MHz and 96 MHz when the instrument probed from 0° to 15° and 30° along the direction of the stromal fibers. Brillouin microscopy did not lose sensitivity to mechanical anisotropy up to 96% water content. We experimentally measured and theoretically modeled how mechanical changes independent of hydration affect frequency shift as a result of corneal cross-linking by isolating an approximately 100 MHz increase in frequency shift following a cross-linking procedure purely due to changes of stromal tissue mechanics.

Conclusions

Brillouin microscopy is sensitive to mechanical anisotropy of the stroma even in highly hydrated corneas. The agreement between model and experimental data suggested a quantitative relationship between Brillouin frequency shift, hydration state of the cornea, and stromal tissue stiffness.

Translational Relevance

The protocol and model validated throughout this study offer a path for comprehensive measurements of corneal mechanics within the clinic; allowing for improved evaluation of the long-term mechanical efficacy of cross-linking procedures.

Depth-Dependent Corneal Biomechanical Properties in Normal and Keratoconic Subjects by Optical Coherence Elastography

Purpose

Compare depth-resolved biomechanical properties in normal and keratoconic corneas in live human subjects using optical coherence elastography (OCE).

Methods

In a prospective series of normal and keratoconus (KC) eyes, a corneal perturbation was applied by a custom swept-source OCE system using a transparent flat lens coupled to force transducers. Cross-correlation was applied to track frame-by-frame OCT speckle displacement. Regional displacements for the anterior and posterior stroma were plotted in force versus displacement (k) graphs. A spatial biomechanical property ratio (k_a/k_p ) was defined by dividing the maximum total displacement by the maximum force for the anterior (k_a ) and posterior cornea (k_p) and was compared between normal and KC groups with the Mann-Whitney U test. Area under the receiver operating characteristics curve (AUROC) for differentiating normal and KC eyes was calculated for k_a/k_p , k_max, and thinnest point of corneal thickness (TPCT).

Results

Thirty-six eyes were analyzed (21 eyes of 12 normal subjects and 15 KC eyes of 12 subjects). The k_a/k_p for the normal group was 1.135 ± 0.07 (mean ± standard deviation) and 1.02 ± 0.08 for the KC group (P < 0.001), indicating a relative deficit in anterior stromal stiffness in KC eyes. AUROC was 0.91 for k_a /k_p , 0.95 for k_max, and 1 for TPCT.

Conclusions

Significant differences in depth-dependent corneal biomechanical properties were observed between normal and KC subjects.

Translational Relevance

OCE was applied for the first time to human KC subjects and revealed alterations in the normal anterior-to-posterior stromal stiffness gradient, a novel and clinically accessible disease biomarker.

Clinical Corneal Optical Coherence Elastography Measurement Precision: Effect of Heartbeat and Respiration

Purpose

Normal physiological movements (e.g., respiration and heartbeat) induce eye motions during clinical measurements of human corneal biomechanical properties using optical coherence elastography (OCE). We quantified the effects of respiratory and cardiac-induced eye motions on clinical corneal OCE measurement precision and repeatability.

Methods

Corneal OCE was performed using low-force, micro-air-pulse tissue stimulation and high-resolution phase-sensitive optical coherence tomography (OCT) imaging. Axial surface displacements of the corneal apex were measured (M-mode) at a 70-kHz sampling rate and three different stimulation pressures (20-60 Pa). Simultaneously, the axial corneal position was tracked with structural OCT imaging, while the heartrate and respiration were monitored over a 90 second period.

Results

Respiratory- and cardiac-induced eye motions have distinctly lower frequency (0.1-1 Hz) and much greater amplitude (up to ± 50 µm movements) than air-pulse-induced corneal tissue deformations (∼250 Hz, <1 µm). The corneal displacements induced during OCE measurements in vivo were -0.41 ± 0.06 µm (n = 22 measurements, coefficient of variation [CV]: 14.6%) and -0.44 ± 0.07 µm (n = 50 measurements, CV: 15.9%), respectively, from two human subjects at 40 Pa stimulation pressure. Observed variation in corneal tissue displacements were not associated with tissue stimulation magnitude, or the amplitude of physiologically induced axial eye motion.

Conclusions

The microsecond timescale and submicron tissue displacements observed during corneal OCE measurements are separable from normal involuntary physiological movements, such as the oculocardiac pulse and respiratory movements.

Translational Relevance

This work advances innovations in biomedical imaging and engineering for clinical diagnostic applications for soft-tissue biomechanical testing.

SNR enhancement in brillouin microspectroscopy using spectrum reconstruction

Brillouin spectroscopy can suffer from low signal-to-noise ratios (SNRs). Such low SNRs can render common data analysis protocols unreliable, especially for SNRs below ∼10. In this work we exploit two denoising algorithms, namely maximum entropy reconstruction (MER) and wavelet analysis (WA), to improve the accuracy and precision in determination of Brillouin shifts and linewidth. Algorithm performance is quantified using Monte-Carlo simulations and benchmarked against the Cramér-Rao lower bound. Superior estimation results are demonstrated even at low SNRs (≥ 1). Denoising is furthermore applied to experimental Brillouin spectra of distilled water at room temperature, allowing the speed of sound in water to be extracted. Experimental and theoretical values were found to be consistent to within ±1% at unity SNR.

Optical Coherence Elastography-Based Corneal Strain Imaging During Low-Amplitude Intraocular Pressure Modulation

Purpose: Optical coherence elastography (OCE) is a promising technique for high-resolution strain imaging in ocular tissues. A major strain-inducing factor in the eye is intraocular pressure (IOP), with diurnal physiological fluctuations reaching up to 5 mmHg. We study herein low-amplitude IOP modulation to assess local corneal strain patterns.

Methods:Ex vivo porcine eye globes were adjusted to an initial IOP of 15 mmHg and subsequently 25 mmHg. Corneal strain was induced by two subsequent pressure cycles, in which IOP was first increased and then decreased, each by a total of 5 mmHg. Two-dimensional optical coherence tomography (2D-OCT) B-scans were recorded after each loading step. Axial strain maps were obtained from magnitude and phase changes and supra-pixel displacements from cross-correlation. The strain detection sensitivity was evaluated in an isotropic material.

Results: Deformations arising from a single 1-mmHg step could be resolved. The largest strain amplitudes (5.11·10⁻³) were observed in the posterior stroma at a low initial IOP. Strain amplitude was 1.34 times higher at 15 mmHg than at 25 mmHg (p = 0.003). Upon IOP increase, the anterior cornea was compressed, whereas the posterior cornea showed axial expansion. Both morphological images and strain maps were sensitive to postmortem time. Strains that are larger than 2.44·10⁻⁵ could be reliably measured.

Conclusions: Low-amplitude IOP modulation, similar to diurnal physiological changes, induced measurable deformations in corneal tissue. Axial strain maps permit a localized comparison of the corneal biomechanical response. Small-strain OCE can likely be extended to other domains.