**Purpose**:
Previous studies have shown that iris mechanical properties may play a role in the pathophysiology of primary angle-closure glaucoma (PACG). Such studies, however, were not conducted in vivo and as such were limited in application and scope, especially for the development of diagnostic methods or new treatment options. The purpose of this study was to quantify in vivo iris mechanical properties both in patients with a history of angle-closure glaucoma and in healthy volunteers.

**Methods**:
We acquired optical coherence tomography scans of anterior segments under standard and dim light conditions. Using a combination of finite element simulation and an inverse fitting algorithm, we quantified the stiffness of the iris.

**Results**:
The irides in the eyes of patients with a history of PACG were significantly stiffer when compared with healthy control irides, a result consistent with ex vivo studies. This result was independent of the compressibility assumption (incompressible: 0.97 ± 0.14 vs. 2.72 ± 0.71, *P* = 0.02; compressible: 0.89 ± 0.13 vs. 2.57 ± 0.69, *P* = 0.02) when comparing the normalized elastic modulus of the iris between patients with PACG and healthy controls.

**Conclusions**:
Our noninvasive, in vivo quantification is free of numerous ethical issues and potential limitations involved with ex vivo examinations. If further studies confirm that the iris stiffness is an omnipresent PACG risk factor and a mechanistic role between increased iris stiffness and angle-closure glaucoma does exist, treatment methods such as lowering the iris stiffness can be developed.

^{1}By 2040, 111.8 million people will be affected by glaucoma.

^{1}In the United States, glaucoma accounts for 9% to 12% (i.e., approximately 120,000) of all cases of blindness.

^{2,3}Two major types of glaucoma, that is, primary open angle glaucoma (POAG) and primary angle-closure glaucoma (PACG), are categorized based on the angle between the iris and the cornea, known as the anterior chamber angle (ACA; Fig. 1A). In PACG, the iris is abnormally bowed toward the anterior, thereby blocking the aqueous humor outflow pathway (Fig. 1B). Although not as prevalent among Caucasian men, PACG primarily affects women and non-Caucasian ethno-racial groups, especially Eskimos and Asians.

^{2,4–6}In fact, because of its high prevalence in East Asia, PACG affects at least as many people as does POAG.

^{7}

^{8,9}The higher pressure in the posterior chamber generates a net force that pushes the iris periphery toward the anterior, leading to a significant narrowing or complete closure of the ACA and consequently blockage of the outflow pathway.

^{10–12}Laser peripheral iridotomy (LPI) remains the most common method for the treatment of pupillary block.

^{13}During this procedure, a laser is used to make holes through the iris to reduce the pressure difference between the anterior and posterior chambers and to restore the normal iris configuration. The LPI procedure, however, may not necessarily widen the ACA in all cases.

^{14,15}In addition, although in many cases the angle opens immediately following the LPI procedure,

^{16}the long-term outcomes may not be as promising, and further surgical procedures may still be necessary to reduce the intraocular pressure.

^{17}Recent studies have identified the success rate of LPI even as low as 24%.

^{18,19}The existence of occludable ACAs following LPI indicates that other mechanisms in addition to pupillary block must be involved in the pathophysiology of PACG.

^{20–27}(Narayanaswamy AK, et al.

*IOVS*2015;56:ARVO E-Abstract 6139). Briefly,

- Histological characterizations of excised iris tissues have shown a higher density of collagen I, a major component of the iris extracellular matrix, in patients suffering from acute angle closure.
^{24} - Real-time quantitative polymerase chain reaction examinations of iris samples obtained from PACG and POAG patients have shown higher level of collagen type 1 alpha 1 chain (
*COL1A1*) mRNA expression in PACG patients.^{25}It is noteworthy that the amount of*COL1A1*expressed in tissues directly scales with the mechanical stiffness across all types of tissues.^{28} - Genomewide studies of more than 10,000 PACG patients in 24 countries have shown an association between the locus of the collagen type XI alpha 1 chain (
*COL11A1*) gene and PACG.^{26,27} - In vivo imaging of the iris cross-sectional area have shown that, in PACG patients, the iris volume does not change significantly after dilation.
^{22,23}Such a phenomenon indicates that although the iris is spongy and compressible in healthy and POAG eyes, in PACG and angle closure suspects, the iris exhibits an incompressible behavior.^{29} - Our in silico models developed based on typical anterior segment geometric parameters have shown that following dilatation, the ACA narrows much more if the iris is incompressible and less spongy.
^{20} - Ex vivo mechanical testing on the iris samples have shown a larger stiffness value in the irides of patients suffering from PACG when compared with those with POAG (Narayanaswamy AK, et al.
*IOVS*2015;56:ARVO E-Abstract 6139).

^{30}Finally, noninvasive in vivo procedures are more suitable for the development of potential new diagnostic methods and/or treatment strategies. As such, we have used clinical anterior segment optical coherence tomography (AS-OCT) images and computer simulation to estimate the mechanical properties of the irides in vivo. With this noninvasive method, we quantified the in vivo iris mechanical properties in patients who had occludable ACA after undergoing LPI and compared them with those of the healthy volunteers.

*n*= 8) was randomly selected from a pool of 30 patients who had undergone LPI for angle closure, but their ACA remained occluded following the procedure. Patients who exhibited abnormal pupillary light reflex response (i.e., their pupils did not dilate) were not included in the study. The second group (

*n*= 8) was randomly selected from a group of healthy volunteers using the same exclusion criteria. The demographic and clinical factors of the enrolled population are shown in Table 1. For each group, images were taken under standard light conditions and following dim light–induced dilation (Fig. 2). For the standard light conditions, a standard light intensity of 350 lux was maintained. No light was targeted toward the eye, and standard room illumination was used in all cases to stimulate the eye. For the dim light conditions, the 350-lux light source was turned off. For all scans, the images were aligned and rotated correctly in such a way that the corneal axis is in the center axis. The tenets of the Declaration of Helsinki were followed for this study.

^{31}:

- AS-OCT images were imported into Solidworks (Dassault Systèmes, Velizy-Villacoublay, France) and appropriately scaled. The boundary of the iris was then manually tracked and a solid two-dimensional model of the iris was created (Fig. 3A).
- Because Abaqus is limited in that it can generate only up to 8-node quadrilateral element mesh, the output of the Abaqus mesh was subsequently imported into an internally developed program in C, in which an extra node was added to each element to generate a 9-node biquadratic finite element mesh.

^{32}and was based on OCT scans. The OCT scans were chosen randomly from either of the nasal, temporal, inferior, or superior configurations, and only one side of the iris was modeled. The iris was modeled as a neo-Hookean solid material governed by the following stress balance equation:

*σ*represents the Cauchy stress tensor. The Cauchy stress tensor was defined by the neo-Hookean (

*σ*

_{NH}) and active sphincter (

*σ*

_{S}) stress tensors:

*G*is the shear modulus;

*ν*is the Poisson's ratio;

*I*is the identity tensor;

*J*is the determinant of

*F*, where

*F*is the deformation gradient tensor; and

*B*is the left Cauchy–Green deformation tensor. The tensors

*F*and

*B*are defined as

*x*is the current position of a material point and

*X*is its resting position. The sphincter stress

*σ*

_{S}, which was applied in the sphincter region, is defined by the scalar active muscle contraction stress (

*σ*

_{Act}) and the unit vector representing the direction of the sphincter muscle (

*e*

*).*

_{θ}*σ*

_{Act}, based on the estimated values in our previous models.

^{32}

^{8,9,20,31,32}The Newton–Raphson iteration and the direct linear solver MUltifrontal Massively Parallel sparse direct Solver (MUMPS; University of Bordeaux, Bordeaux, Nouvelle-Aquitaine, France) were employed in solving the nonlinear algebraic equations.

^{33}To examine the effect of the mesh size on the problem outcomes, a mesh convergence study was performed similar to the one employed in our previous study.

^{34}

^{35}that was combined with our finite element model of the iris relaxation during pupil constriction was used to calculate the iris shear modulus

*G*. We have used a similar approach in our previous study for the estimation of

*G*in the trabecular meshwork.

^{31}A basic overview of the inverse approach is shown in Figure 4. The objective function was defined as the absolute difference between distance from the root to the tip of the iris measured in the experiments (

*d*) and the simulated distance between the root and tip of the iris (

_{exp}*d*).

_{sim}*G*for each set of the iris images, two cases using different values for the Poisson's ratio

*ν*were considered: (1) the iris was considered as a compressible material with a Poisson's ratio of

*ν*= 0.35 and (2) the iris was considered as a nearly incompressible material; hence, a Poisson's ratio of

*ν*= 0.49 was used.

*G*was chosen between 1 and 91 kPa. All simulations were performed using an HP Intel Xeon machine (Intel, Santa Clara, CA, USA) at the Ohio Supercomputing Center (Columbus, OH, USA).

^{36}The maximum number of generations (i.e., the iteration levels for the inverse algorithm) was set to 50 to allow for the solution to converge. At the end of 50 iterations, convergence on a solution for a particular objective function was achieved when the individual values of

*G*in the population set came satisfyingly close to each other, meaning that the maximum difference between each population and the best population was below a threshold value at the end of 50 iterations; otherwise, the population was defined not to have converged to a unique solution.

^{37}

*σ*

_{Act}) to obtain a normalized elastic modulus value (

*E*′).

*σ*

_{Act}has any influence on the calculated normalized elastic modules, the value of

*σ*

_{Act}was perturbed in two cases.

*t*-test was performed to compute the statistical significance (

*P*< 0.05); (2) no significant difference was found to exist between the mean normalized elastic moduli of the compressible/nearly incompressible healthy and compressible/nearly incompressible post-LPI irides, and a one-tailed Student's

*t*-test was performed to compute the statistical significance (

*P*< 0.05).

*d*) in normal patients under standard light and following dim light dilation was found to be 4.48 ± 0.15 mm (average ± standard error) and 4.00 ± 0.14 mm, respectively. The average iris length in PACG patients under standard light and following dim light dilation was found to be 4.28 ± 0.09 mm and 4.09 ± 0.07 mm, respectively. No significant difference was found in the iris length of normal and PACG patients during both standard light and following dim light-induced dilation (

*P*> 0.05, using the Student's

*t*-test).

*P*> 0.05, Fig. 6). However, when comparing the healthy and post-LPI cases for different possible combinations of compressible and incompressible cases, we found out that the post-LPI group had a significantly larger normalized elastic modulus than the healthy control group (

*P*< 0.05, Fig. 7).

*G*obtained from the inverse method did not change beyond a total number of 1161 elements (Table 4). A global size of approximately 0.05 mm generated the 1161 elements. The same global size was used for generating the mesh for all samples. As for the validation of the simulated outcomes, perturbation of the model for typical simulated cases showed that only one absolute minimum error was obtained from the inverse modeling approach (Fig. 8). Regarding the effect of the magnitude of the muscle stresses on the calculated stiffness, we have shown that irrespective of the magnitude of the muscle stress, the normalized elastic modulus remains unchanged. We have shown this concept for two different cases of a compressible healthy iris (Fig. 9).

^{24}However, in this study, such a difference was not observed between the chronic angle closure eyes and the healthy control group eyes. Conversely, the amount of

*COL1A1*mRNA expression in PACG iris samples has been shown to be larger when compared with POAG patients.

^{25}Although many of these recent quantitative analyses may suffer from limitations such as small sample size, it has been clearly demonstrated that, from soft tissues such as fat and tissues in the brain to extremely stiff tissues such as those found in bones, the amount of

*COL1A1*mRNA expression directly scales with the mechanical stiffness.

^{28}As such, the results of limited studies on the potential link between the iris stroma extracellular matrix collagen content and angle closure agrees with the outcome of our in vivo stiffness calculation.

^{38}Their models, however, did not include the active deformation of the iris. In addition, due to simplification in the form of constitutive tissue model, less complexity in the tissue behavior was considered. For example, the iris compressibility was not taken into consideration in those models.

^{32}where we observed a pronounced bowing of the iris toward the posterior, a phenomenon that was not present in the AS-OCT images of the human iris after switching from standard light conditions to a dim light level. Other limitations of our study pertained to the identification of the sphincter region in the irides and the magnitude of the sphincter stress. Although a sphincter region was chosen in proportion to the iris size based on previously quantified histological measurements,

^{39}more accurate patient-specific quantifications of the sphincter region could offer more accurate stiffness measurements. We also made the assumption that the magnitude of the sphincter stress was similar in all simulated cases and defined the normalized elastic modulus (normalized stiffness) as a parameter for comparison. Although the significant difference between the normalized values of the stiffness are most likely due to the extracellular collagen content of the stromal layer of the irides, theoretically, same stromal stiffness and different sphincter stresses could produce similar outcomes. Calculating the in vivo patient-specific sphincter muscular stress remains a challenge.

^{40,41}In addition, a sphincter muscle, which consists of circumferentially aligned smooth muscle fibers, is located close to the pupillary edge of the iris. As such, the isotropic assumption of this model can be improved upon by incorporating more advanced anisotropic material models.

^{42}Therefore, an additional nonparametric, one-tailed Mann-Whitney

*U*test was also perfomed to evaluate the significant difference between the healthy control iris and the post-LPI occludable iris. The outcomes were similar to that of the Student's

*t*-test with the exception of one case, where no significant difference was found between the nearly incompressible healthy control iris and the compressible post-LPI occludable iris. This particular case could be considered as a nonrealistic comparison, as our previous study has shown that the normal healthy iris behaves in a compressible manner, whereas a post-LPI occludable iris displays incompressible behavior (Wojcik M, et al.

*IOVS*2017;58:ARVO E-Abstract 2087). Nonethess, some level of caution is required in interpretating the data presented in this study.

^{43,44}However, for our study, this color dependency can be disregarded because we used the same source of illumination for examining the eyes of all individuals who participated in our study.

^{34}; in the additional study, we fixed the Poisson's ratio and only optimized the shear modulus. We found that using the standalone chord length as the objective function was sufficient to accurately calculate the shear modulus of 9 kPa, which was initially used for the generation of numerically simulated displacement data.

*G*. Perturbation of our results showed that the error was confidently minimized at the value obtained from the model. We have shown that simultaneous optimization of more than one parameter required more detailed scans of the anterior segment cross-section and more complex objective functions.

^{34}With more accurate deformation measurements, for example, using digital image correlation, one would expect to improve the current methods.

^{45}Digital image correlation has been used in ocular tissues alone, along with ultrasound, OCT, or microscopy for measurement of scleral and corneal strain.

^{46–48}

**A.D. Pant**, None;

**P. Gogte**, None;

**V. Pathak-Ray**, None;

**S.K. Dorairaj**, None;

**R. Amini**, None

*. 2014; 121: 2081–2090.*

*Ophthalmology**. 2006; 90: 262–267.*

*Br J Ophthalmol**. 2004; 122: 477–485.*

*Arch Ophthalmol**. 1981; 65: 46–49.*

*Br J Ophthalmol**. 1988; 66: 556–564.*

*Acta Ophthalmol**. 1989; 25: 115–119.*

*Zhonghua Yan Ke Za Zhi**. 1996; 80: 389–393.*

*Br J Ophthalmol**. 2012; 45: 2293–2296.*

*J Biomech**. 2010; 132: 071010.*

*J Biomech Eng**. 2002; 17: 69–78.*

*Semin Ophthalmol**. 2003; 14: 70–73.*

*Curr Opin Ophthalmol**. 1968; 52: 19.*

*Br J Ophthalmol**. 1994; 101: 1749–1758.*

*Ophthalmology**. 2007; 114: 494–500.*

*Ophthalmology**. 2007; 114: 1513–1519.*

*Ophthalmology**. 2003; 110: 630–638.*

*Ophthalmology**. 2000; 107: 2300–2304.*

*Ophthalmology**. 2014; 55: 6795–6801.*

*Invest Ophthalmol Vis Sci**. 2013; 28: 4–8.*

*Semin Ophthalmol**. 2013; 54: 2977–2984.*

*Invest Ophthalmol Vis Sci**. 2001; 123: 540–547.*

*J Biomech Eng**. 2009; 18: 173–179.*

*J Glaucoma**. 2010; 117: 3–10.*

*Ophthalmology**. 2008; 17: 386–392.*

*J Glaucoma**. 2016; 44: 684–692.*

*Clin Exp Ophthalmol**. 2012; 44: 1142–1146.*

*Nat Genet**. 2016; 48: 556–562.*

*Nat Genet**. 2013; 341: 1240104.*

*Science**. 2010; 117: 1–2.*

*Ophthalmology**. 2001; 81: 69.*

*MMVR**. 2017; 1: 100–111.*

*J Model Ophthalmol**. 2012; 53: 1188–1194.*

*Invest Ophthalmol Vis Sci**. 2000; 184: 501–520.*

*Comput Methods Appl Mech Eng**. 2018; 140: 074502.*

*J Biomech Eng**. 1997; 11: 341–359.*

*J Global Optim**Ohio Supercomputer Center*. Available at: http://osc.edu/ark:/19495/f5s1ph73.

*. Vol 143. Springer, Berlin, Heidelberg; 2008: 111–138.*

*Advances in Differential Evolution. Studies in Computational Intelligence**. 2014; 13: 66.*

*Biomed Eng Online**. Philadelphia: Saunders; 1971.*

*Histology of the Human Eye: An Atlas and Textbook**. Philadelphia, PA: Saunders Ltd.; 2015.*

*The Eye E-Book: Basic Sciences in Practice**. 2011; 93: 475–481.*

*Exp Eye Res**. 2006; 36: 171–176.*

*Turk J Med Sci**. 2009; 116: 1564–1573.*

*Ophthalmology**. 2011; 2: 10.*

*Front Neurol**. 2009; 20: 062001.*

*Meas Sci Technol**. 2014; 40: 1041–1047.*

*J Cataract Refract Surg**. 2016; 138: 021015.*

*J Biomech Eng**. 2017; 53: 123–139.*

*Acta Biomater*