Adduction induces large optic nerve head deformations in subjects with normal-tension glaucoma.

PURPOSE: To assess intraocular pressure (IOP)-induced and gaze-induced optic nerve head (ONH) strains in subjects with high-tension glaucoma (HTG) and normal-tension glaucoma (NTG). DESIGN: Clinic-based cross-sectional study. METHODS: The ONH from one eye of 228 subjects (114 subjects with HTG (pre-treatment IOP≥21 mm Hg) and 114 with NTG (pre-treatment IOP<21 mm Hg)) was imaged with optical coherence tomography (OCT) under the following conditions: (1) OCT primary gaze, (2) 20° adduction from OCT primary gaze, (3) 20° abduction from OCT primary gaze and (4) OCT primary gaze with acute IOP elevation (to approximately 33 mm Hg). We then performed digital volume correlation analysis to quantify IOP-induced and gaze-induced ONH tissue deformations and strains. RESULTS: Across all subjects, adduction generated high effective strain (4.4%±2.3%) in the LC tissue with no significant difference (p>0.05) with those induced by IOP elevation (4.5%±2.4%); while abduction generated significantly lower (p=0.01) effective strain (3.1%±1.9%). The lamina cribrosa (LC) of HTG subjects exhibited significantly higher effective strain than those of NTG subjects under IOP elevation (HTG: 4.6%±1.7% vs NTG: 4.1%±1.5%, p<0.05). Conversely, the LC of NTG subjects exhibited significantly higher effective strain than those of HTG subjects under adduction (NTG: 4.9%±1.9% vs HTG: 4.0%±1.4%, p<0.05). CONCLUSION: We found that NTG subjects experienced higher strains due to adduction than HTG subjects, while HTG subjects experienced higher strain due to IOP elevation than NTG subjects-and that these differences were most pronounced in the LC tissue.

The Structural Layers of the Porcine Iris Exhibit Inherently Different Biomechanical Properties.

Purpose: The purpose of this study was to isolate the structural components of the ex vivo porcine iris tissue and to determine their biomechanical properties. Methods: The porcine stroma and dilator tissues were separated, and their dimensions were assessed using optical coherence tomography (OCT). The stroma underwent flow test (n = 32) to evaluate for permeability using Darcy's Law (ΔP = 2000 Pa, A = 0.0391 mm2), and both tissues underwent stress relaxation experiments (ε = 0.5 with initial ramp of δε = 0.1) to evaluate for their viscoelastic behaviours (n = 28). Viscoelasticity was characterized by the parameters ß (half width of the Gaussian distribution), τm (mean relaxation time constant), E0 (instantaneous modulus), and E∞ (equilibrium modulus). Results: For the stroma, the hydraulic permeability was 9.49 ± 3.05 × 10-6 mm2/Pa · s, and the viscoelastic parameters were ß = 2.50 ± 1.40, and τm = 7.43 ± 4.96 s, with the 2 moduli calculated to be E0 = 14.14 ± 6.44 kPa and E∞ = 6.08 ± 2.74 kPa. For the dilator tissue, the viscoelastic parameters were ß = 2.06 ± 1.33 and τm = 1.28 ± 1.27 seconds, with the 2 moduli calculated to be E0 = 9.16 ± 3.03 kPa and E∞ = 5.54 ± 1.98 kPa. Conclusions: We have established a new protocol to evaluate the biomechanical properties of the structural layers of the iris. Overall, the stroma was permeable and exhibited smaller moduli than those of the dilator muscle. An improved characterization of iris biomechanics may form the basis to further our understanding of angle closure glaucoma.

Discriminating Between Papilledema and Optic Disc Drusen Using 3D Structural Analysis of the Optic Nerve Head.

BACKGROUND AND OBJECTIVES: The distinction of papilledema from other optic nerve head (ONH) lesions mimicking papilledema, such as optic disc drusen (ODD), can be difficult in clinical practice. We aimed the following: (1) to develop a deep learning algorithm to automatically identify major structures of the ONH in 3-dimensional (3D) optical coherence tomography (OCT) scans and (2) to exploit such information to robustly differentiate among ODD, papilledema, and healthy ONHs. METHODS: This was a cross-sectional comparative study of patients from 3 sites (Singapore, Denmark, and Australia) with confirmed ODD, those with papilledema due to raised intracranial pressure, and healthy controls. Raster scans of the ONH were acquired using OCT imaging and then processed to improve deep-tissue visibility. First, a deep learning algorithm was developed to identify major ONH tissues and ODD regions. The performance of our algorithm was assessed using the Dice coefficient. Second, a classification algorithm (random forest) was designed to perform 3-class classifications (1: ODD, 2: papilledema, and 3: healthy ONHs) strictly from their drusen and prelamina swelling scores (calculated from the segmentations). To assess performance, we reported the area under the receiver operating characteristic curve for each class. RESULTS: A total of 241 patients (256 imaged ONHs, including 105 ODD, 51 papilledema, and 100 healthy ONHs) were retrospectively included in this study. Using OCT images of the ONH, our segmentation algorithm was able to isolate neural and connective tissues and ODD regions/conglomerates whenever present. This was confirmed by an averaged Dice coefficient of 0.93 ± 0.03 on the test set, corresponding to good segmentation performance. Classification was achieved with high AUCs, that is, 0.99 ± 0.001 for the detection of ODD, 0.99 ± 0.005 for the detection of papilledema, and 0.98 ± 0.01 for the detection of healthy ONHs. DISCUSSION: Our artificial intelligence approach can discriminate ODD from papilledema, strictly using a single OCT scan of the ONH. Our classification performance was very good in the studied population, with the caveat that validation in a much larger population is warranted. Our approach may have the potential to establish OCT imaging as one of the mainstays of diagnostic imaging for ONH disorders in neuro-ophthalmology, in addition to fundus photography.

Differing Associations between Optic Nerve Head Strains and Visual Field Loss in Patients with Normal- and High-Tension Glaucoma.

PURPOSE: To study the associations between optic nerve head (ONH) strains under intraocular pressure (IOP) elevation with retinal sensitivity in patients with glaucoma. DESIGN: Clinic-based cross-sectional study. PARTICIPANTS: Two hundred twenty-nine patients with primary open-angle glaucoma (subdivided into 115 patients with high-tension glaucoma [HTG] and 114 patients with normal-tension glaucoma [NTG]). METHODS: For 1 eye of each patient, we imaged the ONH using spectral-domain OCT under the following conditions: (1) primary gaze and (2) primary gaze with acute IOP elevation (to approximately 35 mmHg) achieved through ophthalmodynamometry. A 3-dimensional strain-mapping algorithm was applied to quantify IOP-induced ONH tissue strain (i.e., deformation) in each ONH. Strains in the prelaminar tissue (PLT), the retina, the choroid, the sclera, and the lamina cribrosa (LC) were associated (using linear regression) with measures of retinal sensitivity from the 24-2 Humphrey visual field test (Carl Zeiss Meditec). This was performed globally, then locally according to a previously published regionalization scheme. MAIN OUTCOME MEASURES: Associations between ONH strains and values of retinal sensitivity from visual field testing. RESULTS: For patients with HTG, we found (1) significant negative linear associations between ONH strains and retinal sensitivity (P < 0.001; on average, a 1% increase in ONH strains corresponded to a decrease in retinal sensitivity of 1.1 decibels [dB]), (2) that high-strain regions colocalized with anatomically mapped regions of high visual field loss, and (3) that the strongest negative associations were observed in the superior region and in the PLT. In contrast, for patients with NTG, no significant associations between strains and retinal sensitivity were observed except in the superotemporal region of the LC. CONCLUSIONS: We found significant negative associations between IOP-induced ONH strains and retinal sensitivity in a relatively large glaucoma cohort. Specifically, patients with HTG who experienced higher ONH strains were more likely to exhibit lower retinal sensitivities. Interestingly, this trend in general was less pronounced in patients with NTG, which could suggest a distinct pathophysiologic relationship between the two glaucoma subtypes.

The three-dimensional structural configuration of the central retinal vessel trunk and branches as a glaucoma biomarker.

PURPOSE: To assess whether the 3-dimensional (3D) structural configuration of the central retinal vessel trunk and its branches (CRVT&B) could be used as a diagnostic marker for glaucoma. DESIGN: Retrospective, deep-learning approach diagnosis study. METHODS: We trained a deep learning network to automatically segment the CRVT&B from the B-scans of the optical coherence tomography (OCT) volume of the optic nerve head. Subsequently, 2 different approaches were used for glaucoma diagnosis using the structural configuration of the CRVT&B as extracted from the OCT volumes. In the first approach, we aimed to provide a diagnosis using only 3D convolutional neural networks and the 3D structure of the CRVT&B. For the second approach, we projected the 3D structure of the CRVT&B orthographically onto sagittal, frontal, and transverse planes to obtain 3 two-dimensional (2D) images, and then a 2D convolutional neural network was used for diagnosis. The segmentation accuracy was evaluated using the Dice coefficient, whereas the diagnostic accuracy was assessed using the area under the receiver operating characteristic curves (AUCs). The diagnostic performance of the CRVT&B was also compared with that of retinal nerve fiber layer (RNFL) thickness (calculated in the same cohorts). RESULTS: Our segmentation network was able to efficiently segment retinal blood vessels from OCT scans. On a test set, we achieved a Dice coefficient of 0.81 ± 0.07. The 3D and 2D diagnostic networks were able to differentiate glaucoma from nonglaucoma subjects with accuracies of 82.7% and 83.3%, respectively. The corresponding AUCs for the CRVT&B were 0.89 and 0.90, higher than those obtained with RNFL thickness alone (AUCs ranging from 0.74 to 0.80). CONCLUSIONS: Our work demonstrated that the diagnostic power of the CRVT&B is superior to that of a gold-standard glaucoma parameter, that is, RNFL thickness. Our work also suggested that the major retinal blood vessels form a "skeleton"-the configuration of which may be representative of major optic nerve head structural changes as typically observed with the development and progression of glaucoma.

Describing the Structural Phenotype of the Glaucomatous Optic Nerve Head Using Artificial Intelligence.

PURPOSE: To develop a novel deep-learning approach that can describe the structural phenotype of the glaucomatous optic nerve head (ONH) and can be used as a robust glaucoma diagnosis tool. DESIGN: Retrospective, deep-learning approach diagnosis study. METHOD: We trained a deep-learning network to segment 3 neural-tissue and 4 connective-tissue layers of the ONH. The segmented optical coherence tomography images were then processed by a customized autoencoder network with an additional parallel branch for binary classification. The encoder part of the autoencoder reduced the segmented optical coherence tomography images into a low-dimensional latent space (LS), whereas the decoder and the classification branches reconstructed the images and classified them as glaucoma or nonglaucoma, respectively. We performed principal component analysis on the latent parameters and identified the principal components (PCs). Subsequently, the magnitude of each PC was altered in steps and reported how it impacted the morphology of the ONH. RESULTS: The image reconstruction quality and diagnostic accuracy increased with the size of the LS. With 54 parameters in the LS, the diagnostic accuracy was 92.0 ± 2.3% with a sensitivity of 90.0 ± 2.4% (at 95% specificity), and the corresponding Dice coefficient for the reconstructed images was 0.86 ± 0.04. By changing the magnitudes of PC in steps, we were able to reveal how the morphology of the ONH changes as one transitions from a "nonglaucoma" to a "glaucoma" condition. CONCLUSIONS: Our network was able to identify novel biomarkers of the ONH for glaucoma diagnosis. Specifically, the structural features identified by our algorithm were found to be related to clinical observations of glaucoma.

Changes in Iris Stiffness and Permeability in Primary Angle Closure Glaucoma.

Purpose: To evaluate the biomechanical properties of the iris by evaluating iris movement during pupil constriction and to compare such properties between healthy and primary angle-closure glaucoma (PACG) subjects. Methods: A total of 140 subjects were recruited for this study. In a dark room, the anterior segments of one eye per subject were scanned using anterior segment optical coherence tomography imaging during induced pupil constriction with an external white light source of 1700 lux. Using a custom segmentation code, we automatically isolated the iris segments from the AS-OCT images, which were then discretized and transformed into a three-dimensional point cloud. For each iris, a finite element (FE) mesh was constructed from the point cloud, and an inverse FE simulation was performed to match the clinically observed iris constriction in the AS-OCT images. Through this optimization process, we were able to identify the elastic modulus and permeability of each iris. Results: For all 140 subjects (95 healthy and 45 PACG of Indian/Chinese ethnicity; age 60.2 ± 8.7 for PACG subjects and 57.7 ± 10.1 for healthy subjects), the simulated deformation pattern of the iris during pupil constriction matched well with OCT images. We found that the iris stiffness was higher in PACG than in healthy controls (24.5 ± 8.4 kPa vs. 17.1 ± 6.6 kPa with 40 kPa of active stress specified in the sphincter region; P < 0.001), whereas iris permeability was lower (0.41 ± 0.2 mm2/kPa s vs. 0.55 ± 0.2 mm2/kPa s; p = 0.142). Conclusions: This study suggests that the biomechanical properties of the iris in PACG are different from those in healthy controls. An improved understanding of the biomechanical behavior of the iris may have implications for the understanding and management of angle-closure glaucoma.

Towards label-free 3D segmentation of optical coherence tomography images of the optic nerve head using deep learning.

Recently proposed deep learning (DL) algorithms for the segmentation of optical coherence tomography (OCT) images to quantify the morphological changes to the optic nerve head (ONH) tissues during glaucoma have limited clinical adoption due to their device specific nature and the difficulty in preparing manual segmentations (training data). We propose a DL-based 3D segmentation framework that is easily translatable across OCT devices in a label-free manner (i.e. without the need to manually re-segment data for each device). Specifically, we developed 2 sets of DL networks: the 'enhancer' (enhance OCT image quality and harmonize image characteristics from 3 devices) and the 'ONH-Net' (3D segmentation of 6 ONH tissues). We found that only when the 'enhancer' was used to preprocess the OCT images, the 'ONH-Net' trained on any of the 3 devices successfully segmented ONH tissues from the other two unseen devices with high performance (Dice coefficients > 0.92). We demonstrate that is possible to automatically segment OCT images from new devices without ever needing manual segmentation data from them.

A viscoelastic framework for inflation testing of gastrointestinal tissue.

Gastrointestinal (GI) diseases are often associated with hypertrophy of the layers of the GI wall, along with dilatation and a denervation of smooth muscle cells which alters the biomechanical properties of the tissue. 'Balloon distension' is a specialised experimental protocol performed on hollow organs to investigate their biomechanical properties. A balloon is inserted and pressurized during this procedure and the change in external diameter is monitored as a function of the applied pressure. A hyperelastic framework is often used in this context to evaluate the stresses in the wall after inflation. However, this only gives an idea about the final equilibrium state of the tissue, whereas the intermediate states of deformations are overlooked. GI soft tissues are viscoelastic, thus, the stress values during inflation are loading rate dependent and are higher than the equilibrium state values. Therefore, it is necessary to consider the time- and rate-dependent material properties during a balloon distension test. The aim of this work was to develop a viscoelastic framework for interpreting balloon distension experiments under finite deformation. To demonstrate the efficacy of the framework, it was used to recreate experimental results from intestinal and colonic tissue segments. In all cases, the simulation results were well matched (R2>0.9) with the experimental data.

A finite element approach for gastrointestinal tissue mechanics.

The biomechanical properties of gastrointestinal (GI) tissue play a significant role in the normal functioning of the organ. GI soft tissues exhibit a highly nonlinear rate- and time-dependent stress-strain behaviour. In recent years, many constitutive relations have been proposed to characterize these properties. However, a constitutive relation is not sufficient to analyse the biomechanics at the organ level with complex loading and boundary conditions. Hence, for a refined mechanical analysis, a finite element (FE) implementation of the constitutive relation is needed. Here, we propose an FE implementation of a finite nonlinear hyperviscoelastic model suitable for soft biological tissues. The FE model has been validated at first by comparing its results with the analytical solutions of a standard linear solid, and then it has been used to recreate experimental observations performed on tissue strips obtained from different animals. We have also proposed a method, in this work, to construct a residually stressed FE model so that the consequences of residual stresses on GI mechanics can be examined. Our FE formulation was able to capture the nonlinear soft tissue properties and also demonstrated that the addition of residual stresses reduces stress concentrations and the stress gradient in the GI wall.

A finite nonlinear hyper-viscoelastic model for soft biological tissues.

Soft tissues exhibit highly nonlinear rate and time-dependent stress-strain behaviour. Strain and strain rate dependencies are often modelled using a hyperelastic model and a discrete (standard linear solid) or continuous spectrum (quasi-linear) viscoelastic model, respectively. However, these models are unable to properly capture the materials characteristics because hyperelastic models are unsuited for time-dependent events, whereas the common viscoelastic models are insufficient for the nonlinear and finite strain viscoelastic tissue responses. The convolution integral based models can demonstrate a finite viscoelastic response; however, their derivations are not consistent with the laws of thermodynamics. The aim of this work was to develop a three-dimensional finite hyper-viscoelastic model for soft tissues using a thermodynamically consistent approach. In addition, a nonlinear function, dependent on strain and strain rate, was adopted to capture the nonlinear variation of viscosity during a loading process. To demonstrate the efficacy and versatility of this approach, the model was used to recreate the experimental results performed on different types of soft tissues. In all the cases, the simulation results were well matched (R2⩾0.99) with the experimental data.