Microstructure visualization of conventional outflow pathway and finite element modeling analysis of trabecular meshwork
© The Author(s) 2016
Published: 28 December 2016
The intraocular pressure (IOP) is maintained through a dynamic equilibrium between the production and drainage of aqueous humor. Elevation of intraocular pressure is mainly caused by the blocking of aqueous humor outflow pathway. Therefore, it is particularly important to study the structure of drainage pathway and the effect of ocular hypertension at the process of aqueous humor outflow.
Conventional drainage pathway of aqueous humor, including trabecular meshwork (TM), Schlemm’s canal (SC) and aqueous vein, were imaged by using trans-scleral imaging method with lateral resolution of 2 μm. For quantitative assessment, the morphological parameters of the TM were measured with different IOP levels via a combination of measurements and simulations.
Images of the TM and the adjacent tissues were obtained. The porosity of TM with normal intraocular pressure varies from 0.63 to 0.74 as the depth increases, while in high IOP it is changed from 0.44 to 0.59. The diameter of aqueous vein varies from 32 to 43 μm, and is smaller than that of SC, which varies from 48 to 64.67 μm.
Our research provides a non-contact method to visualize the microstructure of tissue for clinical examination associated with the blocking of the outflow pathway of aqueous humor in humans. The three-dimensional (3D) microstructures of limbus and the results of finite element modeling analysis of the TM model will serve for the future evaluation of new glaucoma surgical techniques.
KeywordsConventional outflow pathway Trans-scleral imaging method 3D reconstruction Porosity Finite element modeling
Glaucoma, which can result in the damage of optic nerve and defect of visual field, is thought to be the highest irreversible blinding eye disease in the world [1, 2]. Elevation of intraocular pressure (IOP) mainly caused by the blocking of aqueous outflow system is a high risk factor for glaucoma . Aqueous humor leaves the eye mainly from conventional outflow pathway which contains the trabecular meshwork (TM), Schlemm’s canal (SC) and collector channels (CC) before finally draining into aqueous vein. The aqueous humor flow through the TM is driven by a pressure gradient across the tissue. Thus, the pressure difference between inside (intraocular pressure, IOP) and outside (episcleral venous pressure, EVP) provides the force to make aqueous humor flow out of the TM . During primary open angle glaucoma, the greatest resistance to aqueous humor flow is provided by the morphological changes in TM [5–7]. To clarify the production of outflow resistance, it is necessary to obtain the microstructure of drainage pathway.
Consequently, there is a high demand for an effective method to reveal the detailed structure of trabecular meshwork and adjacent tissues. There are several techniques used to inspect ocular tissues. As one of the frequently-used clinical diagnostic equipment, optical coherence tomography (OCT) has a resolution of 15–25 μm , while the light scattering properties of sclera causes imaging blur. Another common technique is ultrasound biomicroscopy (UBM) . The limitations of these two technologies lead to the low resolution of image. It remains unknown about three-dimensional (3D) microstructure of aqueous humor outflow pathway in physiological state. Besides that, because of the very tiny size of tissues in the aqueous outflow system, the detailed 3D information of trabecular meshwork tissues is difficult to be detected by them. Two-photon microscopy (2PM) techniques, including two-photon autofluorescence (2PAF) and second harmonic generation (SHG), which have been widely regarded as a useful tool for 3D imaging, have advantages on depth of field and lateral resolution [10–12]. Meanwhile its high penetration rate can be used for deep tissue imaging with long wavelength excitation. And the long wavelengths of the infrared lasers used in 2PM have considerably less chance to thermally damage . Zoumi et al. showed that 2PM can be used to excite endogenous biological chromophores such as elastin/collagen and melanin . TM with component of collagen fibers and elastin could be imaged by 2PM, which is an effective method for imaging of the TM in its native state and can obviate histology distortion by the cause of fixation and histological processing. As TM locates in deep tissue, 2PM is an available imaging method to obtain image in situ.
In this paper, images of TM and its surrounding tissues were obtained simultaneously. Using these images to measure the size of SC and porosity of TM were particularly close to morphology. Comparing to the structure of 3D reconstruction using histology section, using optical tomography images to reconstruct the TM region of intact eye can get a structure which is based on physiological state. Microstructure under high IOP may be considered in clinical diagnosis with our 3D reconstruction results. 3D structures of TM in different intraocular pressure were reconstructed respectively, which can be used to get the deformation field of the TM after finite element modeling analysis. Moreover, analysis and comparison are made between the results of the morphologic changes of TM from the experiment and numerical simulation. The findings in this report would make sense for a better understanding and wider applications of diagnosis and therapy of glaucoma.
Rat eye tissue
Six SD rats with no ocular phenotype were provided by the animal department (Capital medical University, Beijing, China, IACUC: AEEI-2013-x-123). The weights of the rats were between 300 and 350 g. After euthanasia, we divided them into two groups. Group A (n = 3): eyes were enucleated and held in phosphate-buffered saline (PBS; 8 g/L sodium chloride, 0.2 g/L potassium phosphate monobasic, 2.16 g/L sodium phosphate dibasic heptahydrate, pH 7.4) for intact eye imaging. Group B (n = 3): eyes were enucleated and preserved in 4% paraformaldehyde for 12 h at 4 °C for fixation, and then washed in PBS. After that, anterior segments consisting of cornea, anterior sclera, and trabecular meshwork with iris were prepared by cutting along the equator of the globe and placed in a culture dish with the TM side up for broken eye imaging.
Two-photon microscopy imaging
The entire stack of images from the 2PM was used to reconstruct the TM structures. The format of the images of TM is transformed into software MIMICS10.0 ((Materialise Inc., Leuven, Belgium). In this study, the images are processed using the tools of threshold segmentation, region glowing and edit. Then, the three remodeling model of TM region can be obtained.
The TM and sclera materials were assumed to be linearly elastic and isotropic. The simulated pressures were set to the pressure drop between IOP and EVP in the rat eye. The material properties of TM model were a Poisson ratio of 0.48 and Young’s modulus of 162 Pa . The simulated pressures were set from 10 to 160 Pa with an interval of 10 Pa. The direction of pressure is perpendicular to the plane in parallel to the inner wall of SC. Pressure was applied on the TM and the sclera edge was set as the constraint for the analysis of pressure drop across the TM tissue.
Trabecular meshwork imaging and porosity
Schlemm’s canal imaging
3D reconstruction of TM and Schlemm’s canal
Morphological analysis of TM
Although analysis based on experiment has contributed to the understanding of the morphological structure of the TM, the detailed structural information of TM during the progression of intraocular pressure is not well-known. For the requirement of analysis, 3D finite element method (FEM) analysis on TM model was carried out using ANASYS 15. In the FEM analysis, if the TM model was combined with the sclera, TM biomechanical properties with more constraints could be acquired, resulting in a more accurate model. On this basis, 3D TM models included a 3D sclera model in the rat eye. The entity model of TM was processed by meshing to fit in FEM. As the final assembled model, shown in Fig. 7c, is the 3D remodeling model of TM. As can be seen in Fig. 7d, is the TM model of intact eye after the meshing, it was tested with the 193,566 elements and 575,076 nodes.
Aqueous vein imaging and 3D reconstruction
In our study, 2PM was used to image the trabecular drainage channels of intact enucleated rat eyes. Due to the intrinsic properties of drainage tissue, autofluorescence from collagen fibers can be used to image the microstructure of the pathway. Compared to external fluorescent label, autofluorescence can avoid tissue damage caused by fluorochrome. Meanwhile, 2PM is suited to tissue imaging under the sclera, and has been successfully used to image living skin at a depth of 350 μm by visualizing the 2PAF of the skin’s extracellular matrix and melanin . Most importantly of all, we were able to successfully get the high resolution images under high scattering tissue of eye. Our imaging had a lateral resolution of 2 μm, and a depth resolution of 1 μm. It is found that the high resolution imaging can achieve resolution of tissue section, and can avoid tissue distortions caused by infusion of fixatives.
Aqueous humor drains into aqueous vein through TM region firstly. In this progress, TM and SC are both significant to fluid drainage. Obtaining the microstructures of these tissues is useful to investigate the production mechanism of outflow resistance in the whole draining tunnel. Imaging of TM in human cadaver eyes and obtaining the location of SC and TM in intact enucleated mouse eyes using 2PM has been reported respectively [20, 21].In this research work, we obtained the images of TM, the surrounding SC and aqueous veins totally. It also shows that rat and mouse eye tissue as a suitable animal model is useful for tussue imaging in the outflow pathway of aqueous humor.
Moreover, we calculated the porosity of TM in the intact eye and found that porosity became larger as depth increases. The porosity of TM in the broken eye was smaller compared to the intact eye. It can be seen that IOP plays an important role in the structure of TM. Force balance between IOP and episcleral venous pressure can maintain the physiological state of aqueous outflow pathway. Breaking the force balance can lead to the collapse of tissues in the trabecular drainage channels and then affect the IOP in return.
Some researches reconstructed limbus using histologic section with manual add CC to obtain the distribution of pressure and flow velocity in the model of eye and results offered were comparable to physiology references . We also performed 3D reconstructions of serial images of limbus tissue to obtain the 3D structures of aqueous humor outflow pathway. From stereo image matching, we found that 3D TM model with some pores of various sizes was consistent with the form obtained from histologic sections. Which means that our 3D model based on ex vivo rat eye approaches to the physiology state, and can be used to simulate the physiological flow in TM region. Such observations of TM are interesting and important with pressure dependent changes in the trabecular meshwork.
There are several problems to be solved in the next step. Depth of penetration decreased image resolution due to scattering and tissue in homogeneities . Previous research used gonioscopic lens to reduce highly scattering of tissue and obtain the structures of TM in the porcine eye . The porcine eye has a series of CCs instead of a single large circumferential SC . Compared with porcine eye, rat eye is smaller as well as the size of SC and aqueous veins [26, 27]. The resolution of image is the limitation of imaging quality in the drainage channels. Some researchers indicated that filtration patterns in human TM appeared segmental over both macro-scale and micro-scale dimensions . To get a comprehensive description of aqueous humor outflow pathway, we will image the circumferential TM/SC/aqueous vein patterns instead of local regions. In order to provide sufficient evidence for diagnosis and therapy of glaucoma, accurately assessing structural changes of TM is very significant. To investigate the relation between the structural changes of limbus and the intraocular pressure, we performed 3D reconstructions of the TM region. In addition, static analysis of the deformation of a TM model was performed using the FEM. Moreover, the simulated results were combined with real changes obtained in our experiment. Early detection and clinical intervention in the prevention of glaucoma disease will greatly benefit from the images with high resolution and the results of deformation analysis. However, there are some limitations of the FEM analysis in this study. In the FEM analysis, if the TM model was combined with Schlemm’s canal and aqueous vein, TM biomechanical properties could be acquired, resulting in a more accurate model. To further modeling research, the model will be refined and the method will be improved, which may bring the hardness to interpret in terms of pathology. Potentially, it will become a valuable method for finding the resistance to drainage of aqueous humor and comprehensive assessment of treatment efficacy.
Trans-scleral imaging method was used to obtain morphology of trabecular drainage channels, including TM/SC/aqueous vein, in the intact eye at different depths without staining and fixing. We also reconstructed the tissues of the drainage pathway to get 3D morphology and structure. Furthermore, we simulated and analyzed the structural changes of TM with different IOP levels by the FEM analysis. Combining the experiment results and simulations, our research provides a reliable method to visualize aqueous humor outflow pathway and a potential technique for clinical examination associated with the blocking of aqueous outflow pathway in humans.
JZ was responsible for the computational modeling, the statistical analysis and overall investigation. LR and XM were responsible for the data analysis part. JZ, QX and WZ were responsible for the experiment part. ZL was responsible for the design. All authors have made substantial contributions to conception and design, or acquisition of data, or analysis and interpretation of data and have been involved in drafting the manuscript or revising it critically for important intellectual content. Each author has participated sufficiently in the work to take public responsibility for appropriate portions of the content. All authors read and approved the final manuscript.
This work was supported by the National Natural Science Foundation of China (Nos. 31070840, 10802053, 81201123), the Natural Science Foundation of Beijing (No. 3122010), the Natural Science Foundation of Guangdong Province (No. S2012040006576), the Shenzhen Science and Technology Innovation Committee Grant (No. KQCX20120816155352228), and the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (No. PHR201110506). We thank Dr. Qingjun Lu in Tongren Hospital of Beijing for valuable discussion.
The authors declare that they have no other competing interest.
About this supplement
This article has been published as part of BioMedical Engineering OnLine Volume 15 Supplement 2, 2016. Computational and experimental methods for biological research: cardiovascular diseases and beyond. The full contents of the supplement are available online http://biomedical-engineering-online.biomedcentral.com/articles/supplements/volume-15-supplement-2.
Availability of data and materials
Data is accessible upon request.
Ethics approval and consent to participate
All the animal experiments reported in our study were carried out in accordance with the National Institute of Health Guide for the Care and Use of Laboratory Animals and approved by the Institutional Animal Care and Use Committee of China.
Publication of this article was funded by the National Natural Science Foundation of China (No. 31570952).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Weinreb RN, Khaw PT. Primary open-angle glaucoma. Lancet. 2004;363:1711–20.View ArticleGoogle Scholar
- Tham YC, Li X, Wong TY, Quigley HA, Aung T, Cheng CY. Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis. Ophthalmology. 2014;121(11):2081–90.View ArticleGoogle Scholar
- Llobet A, Gasull X, Gual A. Understanding trabecular meshwork physiology: a key to the control of intraocular pressure. News Physiol Sci. 2003;18:205–9.Google Scholar
- Mäepea O, Bill A. The pressures in the episcleral veins, Schlemm’s canal and the trabecular meshwork in monkeys: effects of changes in intraocular pressure. Exp Eye Res. 1989;49(4):645–63.View ArticleGoogle Scholar
- Kass MA, Heuer DK, Higginbotham EJ, Johnson CA, Keltner JL, Miller JP, Parrish RK 2nd, Wilson MR, Gordon MO. The Ocular Hypertension Treatment Study: a randomized trial determines that topical ocular hypotensive medication delays or prevents the onset of primary open-angle glaucoma. Arch Ophthalmol. 2002;120(6):701–13.View ArticleGoogle Scholar
- Barany EH. In vitro studies on the resistance to flow through the angle of the anterior chamber. Acta Soc Med Ups. 1954;59(3–4):260–76.Google Scholar
- Johnson M. What controls aqueous humour outflow resistance? Exp Eye Res. 2006;82(4):545–57.View ArticleGoogle Scholar
- Leung CK, Li H, Weinreb RN, et al. Anterior chamber angle measurement with anterior segment optical coherence tomography: a comparison between slit lamp OCT and Visante OCT. Invest Ophthalmol Vis Sci. 2008;49:3469–74.View ArticleGoogle Scholar
- Knspik DA, Starkoski B, Pavlin CJ, Foster FS. A 100–200 MHz ultrasound biomicroscope. IEEE Trans Ultrason Ferroelectr Freq Control. 2000;47:1540–9.View ArticleGoogle Scholar
- Kobat D, Horton NG, Xu C. In vivo two-photon microscopy to 1.6-mm depth in mouse cortex. J Biomed Opt. 2011;16(10):106014.View ArticleGoogle Scholar
- Morishige N, Wahlert AJ, Kenney MC, Brown DJ, Kawamoto K, Chikama T, Nishida T, Jester JV. Second-harmonic imaging microscopy of normal human and keratoconus cornea. Invest Ophthalmol Vis Sci. 2007;48:1087–94.View ArticleGoogle Scholar
- Zhang Xianzeng, Liu Nenrong, Mak PU, et al. Three-dimensional segmentation and quantitative measurement of the aqueous outflow system of intact mouse eyes based on spectral two-photon microscopy techniques. Invest Ophthalmol Vis Sci. 2016;57:3159–67.View ArticleGoogle Scholar
- Masihzadeh O, Ammar DA, Kahook MY, Lei TC. Coherent anti-stokes Raman scattering (CARS) microscopy: a novel technique for imaging the retina. Invest Ophthalmol Vis Sci. 2013;54(5):3094–101.View ArticleGoogle Scholar
- Zoumi A, Yeh A, Tromberg BJ. Imaging cells and extracellular matrix in vivo by using second-harmonic generation and two-photon excited fluorescence. Proc Natl Acad Sci USA. 2002;99(17):11014–9.View ArticleGoogle Scholar
- Li D, Zheng W, Qu JY. Two-photon autofluorescence microscopy of multicolor excitation. Opt Lett. 2009;34(2):202–4.View ArticleGoogle Scholar
- Huang Jianyong, Camras Lucinda J, Yuan Fan. Mechanical analysis of rat trabecular meshwork. Soft Matter. 2015;11(14):2857–65.View ArticleGoogle Scholar
- Chu ER, Gonzalez JM Jr, Tan JCH. Tissue-based imaging model of human trabecular meshwork. J Ocul Pharmacol Ther. 2014;30(2):1–10.Google Scholar
- Tamm ER. The trabecular meshwork outflow pathways: structural and functional aspects. Exp Eye Res. 2009;88:648–55.View ArticleGoogle Scholar
- Rajadhyaksha M, Gonzalez S, Zavislan JM, Anderson RR, Webb RH. In vivo confocal scanning laser microscopy of human skin II: advances in instrumentation and comparison with histology. J Invest Dermatol. 1999;113:293–303.View ArticleGoogle Scholar
- Ammar David A, Lei Tim C, Gibson Emily A, Kahook Malik Y. Two-photon imaging of the trabecular meshwork. Mol Vis. 2010;16:935–44.Google Scholar
- Johnson AW, Ammar DA. Two-photon imaging of the mouse eye. Invest Ophthalmol Vis Sci. 2011;52(7):4098–105.View ArticleGoogle Scholar
- Villamarina Adan, Roya Sylvain, Vardoulis Reda Hasballa Orestis, Reymond Philippe, Stergiopulos Nikolaos. 3D simulation of the aqueous flow in the human eye. Med Eng Phys. 2012;34:1462–70.View ArticleGoogle Scholar
- Masihzadeh O, Lei TC, Ammar DA, Kahook MY, Gibson EA. A multiphoton microscope platform for imaging the mouse eye. Mol Vis. 2012;18:1840–8.Google Scholar
- Masihzadeh Omid, Ammar David A, Kahook Malik Y, Gibson Emily A, Lei Tim C. Direct trabecular meshwork imaging in porcine eyes through multiphoton gonioscopy. J Biomed Opt. 2013;18(3):036009.View ArticleGoogle Scholar
- McMenamin PG, Steptoe RJ. Normal anatomy of the aqueous humour outflow system in the domestic pig eye. J Anat. 1991;178:65–77.Google Scholar
- Hann Cheryl R, Bentley Michael D, Vercnocke Andrew, Ritman Erik L, Fautsch Michael P. Imaging the aqueous humor outflow pathway in human eyes by three-dimensional micro-computed tomography (3D micro-CT). Exp Eye Res. 2011;92(2):104–11.View ArticleGoogle Scholar
- Irshad FA, Mayfield MS, Zurakowski D, Ayyala RS. Variation in Schlemm’s canal diameter and location by ultrasound biomicroscopy. Ophthalmology. 2010;117(5):916–20.View ArticleGoogle Scholar
- Chang JYH, Folz SL, Laryea SN, Overby DR. Multi-scale analysis of segmental outflow patterns in human trabecular meshwork with changing intraocular pressure. J Ocul Pharmacol Ther. 2014;30(2–3):213–23.View ArticleGoogle Scholar