Capturing sclera anisotropy using direct collagen fiber models. Linking microstructure to macroscopic mechanical properties.

Because of the crucial role of collagen fibers on soft tissue mechanics, there is great interest in techniques to incorporate them in computational models. Recently we introduced a direct fiber modeling approach for sclera based on representing the long-interwoven fibers. Our method differs from the conventional continuum approach to modeling sclera that homogenizes the fibers and describes them as statistical distributions for each element. At large scale our method captured gross collagen fiber bundle architecture from histology and experimental intraocular pressure-induced deformations. At small scale, a direct fiber model of a sclera sample reproduced equi-biaxial experimental behavior from the literature. In this study our goal was a much more challenging task for the direct fiber modeling: to capture specimen-specific 3D fiber architecture and anisotropic mechanics of four sclera samples tested under equibiaxial and four non-equibiaxial loadings. Samples of sclera from three eyes were isolated and tested in five biaxial loadings following an approach previously reported. Using microstructural architecture from polarized light microscopy we then created specimen-specific direct fiber models. Model fiber orientations agreed well with the histological information (adjusted R2's>0.89). Through an inverse-fitting process we determined model characteristics, including specimen-specific fiber mechanical properties to match equibiaxial loading. Interestingly, the equibiaxial properties also reproduced all the non-equibiaxial behaviors. These results indicate that the direct fiber modeling method naturally accounted for tissue anisotropy within its fiber structure. Direct fiber modeling is therefore a promising approach to understand how macroscopic behavior arises from microstructure.

Comparing continuum and direct fiber models of soft tissues: An ocular biomechanics example reveals that continuum models may artificially disrupt the strains at both the tissue and fiber levels.

Collagen fibers are the main load-bearing component of soft tissues but difficult to incorporate into models. Whilst simplified homogenization models suffice for some applications, a thorough mechanistic understanding requires accurate prediction of fiber behavior, including both detailed fiber-level strains and long-distance transmission. Our goal was to compare the performance of a continuum model of the optic nerve head (ONH) built using conventional techniques with a fiber model we recently introduced which explicitly incorporates the complex 3D organization and interaction of collagen fiber bundles [1]. To ensure a fair comparison, we constructed the continuum model with identical geometrical, structural, and boundary specifications as for the fiber model. We found that: 1) although both models accurately matched the intraocular pressure (IOP)-induced globally averaged displacement responses observed in experiments, they diverged significantly in their ability to replicate specific 3D tissue-level strain patterns. Notably, the fiber model faithfully replicated the experimentally observed depth-dependent variability of radial strain, the ring-like pattern of meridional strain, and the radial pattern of circumferential strain, whereas the continuum model failed to do so; 2) the continuum model disrupted the strain transmission along each fiber, a feature captured well by the fiber model. These results demonstrate limitations of the conventional continuum models that rely on homogenization and affine deformation assumptions, which render them incapable of capturing some complex tissue-level and fiber-level deformations. Our results show that the strengths of explicit fiber modeling help capture intricate ONH biomechanics. They potentially also help modeling other fibrous tissues. STATEMENT OF SIGNIFICANCE: Understanding the mechanics of fibrous tissues is crucial for advancing knowledge of various diseases. This study uses the ONH as a test case to compare conventional continuum models with fiber models that explicitly account for the complex fiber structure. We found that the fiber model captured better the biomechanical behaviors at both the tissue level and the fiber level. The insights gained from this study demonstrate the significant potential of fiber models to advance our understanding of not only glaucoma pathophysiology but also other conditions involving fibrous soft tissues. This can contribute to the development of therapeutic strategies across a wide range of applications.

Morphological comparison of astrocytes in the lamina cribrosa and glial lamina.

Purpose: Although the mechanisms underlying glaucomatous neurodegeneration are not yet well understood, cellular and small animal models suggest that LC astrocytes undergo early morphologic and functional changes, indicating their role as early responders to glaucomatous stress. These models, however, lack the LC found in larger animals and humans, leaving the in situ morphology of LC astrocytes and their role in glaucoma initiation underexplored. In this work, we aimed to characterize the morphology of LC astrocytes in situ and determine differences and similarities with astrocytes in the mouse glial lamina (GL), the analogous structure in a prominent glaucoma model. Methods: Astrocytes in the LCs of twenty-two eyes from goats, sheep, and pigs were stochastically labeled via Multicolor DiOlistics and imaged in situ using confocal microscopy. 3D models of DiOlistically-labeled LC astrocytes and hGFAPpr-GFP mouse GL astrocytes were constructed to quantify morphological features related to astrocyte functions. LC and GL astrocyte cross-pore contacts, branching complexity, branch tortuosity, and cell and branch span were compared. Results: LC astrocytes displayed distinct spatial relationships with collagen, greater branching complexity, and higher branch tortuosity compared to GL astrocytes. Despite substantial differences in their anatomical environments, LC and GL astrocytes had similar cell and branch spans. Conclusions: Astrocyte morphology in the LC was characterized through Multicolor DiOlistic labeling. LC and GL astrocytes have both distinct and shared morphological features. Further research is needed to understand the potentially unique roles of LC astrocytes in glaucoma initiation and progression.

Impact of anatomic variability and other vascular factors on lamina cribrosa hypoxia.

Insufficient oxygenation in the lamina cribrosa (LC) may contribute to axonal damage and glaucomatous vision loss. To understand the range of susceptibilities to glaucoma, we aimed to identify key factors influencing LC oxygenation and examine if these factors vary with anatomical differences between eyes. We reconstructed 3D, eye-specific LC vessel networks from histological sections of four healthy monkey eyes. For each network, we generated 125 models varying vessel radius, oxygen consumption rate, and arteriole perfusion pressure. Using hemodynamic and oxygen supply modeling, we predicted blood flow distribution and tissue oxygenation in the LC. ANOVA assessed the significance of each parameter. Our results showed that vessel radius had the greatest influence on LC oxygenation, followed by anatomical variations. Arteriole perfusion pressure and oxygen consumption rate were the third and fourth most influential factors, respectively. The LC regions are well perfused under baseline conditions. These findings highlight the importance of vessel radius and anatomical variation in LC oxygenation, providing insights into LC physiology and pathology. Pathologies affecting vessel radius may increase the risk of LC hypoxia, and anatomical variations could influence susceptibility. Conversely, increased oxygen consumption rates had minimal effects, suggesting that higher metabolic demands, such as those needed to maintain intracellular transport despite elevated intraocular pressure, have limited impact on LC oxygenation.

Comparing continuum and direct fiber models of soft tissues. An ocular biomechanics example reveals that continuum models may artificially disrupt the strains at both the tissue and fiber levels.

Collagen fibers are the main load-bearing component of soft tissues but difficult to incorporate into models. Whilst simplified homogenization models suffice for some applications, a thorough mechanistic understanding requires accurate prediction of fiber behavior, including both detailed fiber-level strains and long-distance transmission. Our goal was to compare the performance of a continuum model of the optic nerve head (ONH) built using conventional techniques with a fiber model we recently introduced which explicitly incorporates the complex 3D organization and interaction of collagen fiber bundles [1]. To ensure a fair comparison, we constructed the continuum model with identical geometrical, structural, and boundary specifications as for the fiber model. We found that: 1) although both models accurately matched the intraocular pressure (IOP)-induced globally averaged displacement responses observed in experiments, they diverged significantly in their ability to replicate specific 3D tissue-level strain patterns. Notably, the fiber model faithfully replicated the experimentally observed depth-dependent variability of radial strain, the ring-like pattern of meridional strain, and the radial pattern of circumferential strain, whereas the continuum model failed to do so; 2) the continuum model disrupted the strain transmission along each fiber, a feature captured well by the fiber model. These results demonstrate limitations of the conventional continuum models that rely on homogenization and affine deformation assumptions, which render them incapable of capturing some complex tissue-level and fiber-level deformations. Our results show that the strengths of explicit fiber modeling help capture intricate ONH biomechanics. They potentially also help modeling other fibrous tissues.

Impact of elevated IOP on lamina cribrosa oxygenation; A combined experimental-computational study on monkeys.

Purpose: Our goal is to evaluate how lamina cribrosa (LC) oxygenation is affected by the tissue distortions resulting from elevated IOP. Design: Experimental study on monkeys. Subjects: Four healthy monkey eyes with OCT scans with IOP of 10 to 50 mmHg, and then with histological sections of LC. Methods: Since in-vivo LC oxygenation measurement is not yet possible, we used 3D eye-specific numerical models of the LC vasculature which we subjected to experimentally-derived tissue deformations. We reconstructed 3D models of the LC vessel networks of 4 healthy monkey eyes from histological sections. We also obtained in-vivo IOP-induced tissue deformations from a healthy monkey using OCT images and digital volume correlation analysis techniques. The extent that LC vessels distort under a given OCT-derived tissue strain remains unknown. We biomechanics-based mapping techniques: cross-sectional and isotropic. The hemodynamics and oxygenations of the four vessel networks were simulated for deformations at several IOPs up to 60mmHg. The results were used to determine the effects of IOP on LC oxygen supply, assorting the extent of tissue mild and severe hypoxia. Main Outcome Measures: IOP-induced deformation, vasculature structure, blood supply, and oxygen supply for LC region. Result: IOP-induced deformations reduced LC oxygenation significantly. More than 20% of LC tissue suffered from mild hypoxia when IOP reached 30 mmHg. Extreme IOP(>50mmHg) led to large severe hypoxia regions (>30%) in the isotropic mapping cases. Conclusion: Our models predicted that moderately elevated IOP can lead to mild hypoxia in a substantial part of the LC, which, if sustained chronically, may contribute to neural tissue damage. For extreme IOP elevations, severe hypoxia was predicted, which would potentially cause more immediate damage. Our findings suggest that despite the remarkable LC vascular robustness, IOP-induced distortions can potentially contribute to glaucomatous neuropathy.

The Robust Lamina Cribrosa Vasculature: Perfusion and Oxygenation Under Elevated Intraocular Pressure.

Purpose: Elevated intraocular pressure (IOP) is thought to cause lamina cribrosa (LC) blood vessel distortions and potentially collapse, adversely affecting LC hemodynamics, reducing oxygenation, and triggering, or contributing to, glaucomatous neuropathy. We assessed the robustness of LC perfusion and oxygenation to vessel collapses. Methods: From histology, we reconstructed three-dimensional eye-specific LC vessel networks of two healthy monkey eyes. We used numerical simulations to estimate LC perfusion and from this the oxygenation. We then evaluated the effects of collapsing a fraction of LC vessels (0%-36%). The collapsed vessels were selected through three scenarios: stochastic (collapse randomly), systematic (collapse strictly by the magnitude of local experimentally determined IOP-induced compression), and mixed (a combination of stochastic and systematic). Results: LC blood flow decreased linearly as vessels collapsed-faster for stochastic and mixed scenarios and slower for the systematic one. LC regions suffering severe hypoxia (oxygen <8 mm Hg) increased proportionally to the collapsed vessels in the systematic scenario. For the stochastic and mixed scenarios, severe hypoxia did not occur until 15% of vessels collapsed. Some LC regions had higher perfusion and oxygenation as vessels collapsed elsewhere. Some severely hypoxic regions maintained normal blood flow. Results were equivalent for both networks and patterns of experimental IOP-induced compression. Conclusions: LC blood flow was sensitive to distributed vessel collapses (stochastic and mixed) and moderately vulnerable to clustered collapses (systematic). Conversely, LC oxygenation was robust to distributed vessel collapses and sensitive to clustered collapses. Locally normal flow does not imply adequate oxygenation. The actual nature of IOP-induced vessel collapse remains unknown.

Lamina Cribrosa Insertions Into the Sclera Are Sparser, Narrower, and More Slanted in the Anterior Lamina.

Purpose: The lamina cribrosa (LC) depends on the sclera for support. The support must be provided through the LC insertions. Although a continuous insertion over the whole LC periphery is often assumed, LC insertions are actually discrete locations where LC collagenous beams meet the sclera. We hypothesized that LC insertions vary in number, size, and shape by quadrant and depth. Methods: Coronal cryosections through the full LCs from six healthy monkey eyes were imaged using instant polarized light microscopy. The images were registered into a stack, on which we manually marked LC insertion outlines, nothing their position in-depth and quadrant (inferior, superior, nasal, or temporal). From the marks, we determined the insertion number, width, angle to the canal wall (90 degrees = perpendicular), and insertion ratio (fraction of LC periphery represented by insertions). Using linear mixed effect models, we determined if the insertion characteristics were associated with depth or quadrant. Results: Insertions in the anterior LC were sparser, narrower, and more slanted than those in deeper LC (P values < 0.001). There were more insertions spanning a larger ratio of the canal wall in the middle LC than in the anterior and posterior (P values < 0.001). In the nasal quadrant, the insertion angles were significantly smaller (P < 0.001). Conclusions: LC insertions vary substantially and significantly over the canal. The sparser, narrower, and more slanted insertions of the anterior-most LC may not provide the robust support afforded by insertions of the middle and posterior LC. These variations may contribute to the progressive deepening of the LC and regional susceptibility to glaucoma.

Direct measurements of collagen fiber recruitment in the posterior pole of the eye.

Collagen is the main load-bearing component of the peripapillary sclera (PPS) and lamina cribrosa (LC) in the eye. Whilst it has been shown that uncrimping and recruitment of the PPS and LC collagen fibers underlies the macro-scale nonlinear stiffening of both tissues with increased intraocular pressure (IOP), the uncrimping and recruitment as a function of local stretch have not been directly measured. This knowledge is crucial to understanding their functions in bearing loads and maintaining tissue integrity. In this project we measured local stretch-induced collagen fiber bundle uncrimping and recruitment curves of the PPS and LC. Thin coronal samples of PPS and LC of sheep eyes were mounted and stretched biaxially quasi-statically using a custom system. At each step, we imaged the PPS and LC with instant polarized light microscopy and quantified pixel-level (1.5 µm/pixel) collagen fiber orientations. We used digital image correlation to measure the local stretch and quantified collagen crimp by the circular standard deviation of fiber orientations, or waviness. Local stretch-recruitment curves of PPS and LC approximated sigmoid functions. PPS recruited more fibers than the LC at the low levels of stretch. At 10% stretch the curves crossed with 75% bundles recruited. The PPS had higher uncrimping rate and waviness remaining after recruitment than the LC: 0.9º vs. 0.6º and 3.1º vs. 2.7º. Altogether our findings support describing fiber recruitment of both PPS and LC with sigmoid curves, with the PPS recruiting faster and at lower stretch than the LC, consistent with a stiffer tissue. STATEMENT OF SIGNIFICANCE: Peripapillary sclera (PPS) and lamina cribrosa (LC) collagen recruitment behaviors are central to the nonlinear mechanical behavior of the posterior pole of the eye. How PPS and LC collagen fibers recruit under stretch is crucial to develop constitutive models of the tissues but remains unclear. We used image-based stretch testing to characterize PPS and LC collagen fiber bundle recruitment under local stretch. We found that fiber-level stretch-recruitment curves of PPS and LC approximated sigmoid functions. PPS recruited more fibers at a low stretch, but at 10% bundle stretch the two curves crossed with 75% bundles recruited. We also found that PPS and LC fibers had different uncrimping rates and non-zero waviness's when recruited.

Fibrous finite element modeling of the optic nerve head region.

The optic nerve head (ONH) region at the posterior pole of the eye is supported by a fibrous structure of collagen fiber bundles. Discerning how the fibrous structure determines the region biomechanics is crucial to understand normal physiology, and the roles of biomechanics on vision loss. The fiber bundles within the ONH structure exhibit complex three-dimensional (3D) organization and continuity across the various tissue components. Computational models of the ONH, however, usually represent collagen fibers in a homogenized fashion without accounting for their continuity across tissues, fibers interacting with each other and other fiber-specific effects in a fibrous structure. We present a fibrous finite element (FFE) model of the ONH that incorporates discrete collagen fiber bundles and their histology-based 3D organization to study ONH biomechanics as a fibrous structure. The FFE model was constructed using polarized light microscopy data of porcine ONH cryosections, representing individual fiber bundles in the sclera, dura and pia maters with beam elements and canal tissues as continuum structures. The FFE model mimics the histological in-plane orientation and width distributions of collagen bundles as well as their continuity across different tissues. Modeling the fiber bundles as linear materials, the FFE model predicts the nonlinear ONH response observed in an inflation experiment from the literature. The model also captures important microstructural mechanisms including fiber interactions and long-range strain transmission among bundles that have not been considered before. The FFE model presented here advances our understanding of the role of fibrous collagen structure in the ONH biomechanics. STATEMENT OF SIGNIFICANCE: The microstructure and mechanics of the optic nerve head (ONH) are central to ocular physiology. Histologically, the ONH region exhibits a complex continuous fibrous structure of collagen bundles. Understanding the role of the fibrous collagen structure on ONH biomechanics requires high-fidelity computational models previously unavailable. We present a computational model of the ONH that incorporates histology-based fibrous collagen structure derived from polarized light microscopy images. The model predictions agree with experiments in the literature, and provide insight into important microstructural mechanisms of fibrous tissue biomechanics, such as long-range strain transmission along fiber bundles. Our model can be used to study the microstructural basis of biomechanical damage and the effects of collagen remodeling in glaucoma.

2D or not 2D? Mapping the in-depth inclination of the collagen fibers of the corneoscleral shell.

The collagen fibers of the corneoscleral shell play a central role in the eye mechanical behavior. Although it is well-known that these fibers form a complex three-dimensional interwoven structure, biomechanical and microstructural studies often assume that the fibers are aligned in-plane with the tissues. This is convenient as it removes the out-of-plane components and allows focusing on the 2D maps of in-plane fiber organization that are often quite complex. The simplification, however, risks missing potentially important aspects of the tissue architecture and mechanics. In the cornea, for instance, fibers with high in-depth inclination have been shown to be mechanically important. Outside the cornea, the in-depth fiber orientations have not been characterized, preventing a deeper understanding of their potential roles. Our goal was to characterize in-depth collagen fiber organization over the whole corneoscleral shell. Seven sheep whole-globe axial sections from eyes fixed at an IOP of 50 mmHg were imaged using polarized light microscopy to measure collagen fiber orientations and density. In-depth fiber orientation distributions and anisotropy (degree of fiber alignment) accounting for fiber density were quantified over the whole sclera and in 15 regions: central cornea, peripheral cornea, limbus, anterior equator, equator, posterior equator, posterior sclera and peripapillary sclera on both nasal and temporal sides. Orientation distributions were fitted using a combination of a uniform distribution and a sum of π-periodic von Mises distributions, each with three parameters: primary orientation µ, fiber concentration factor k, and weighting factor a. To study the features of fibers that are not in-plane, i.e., fiber inclination, we quantified the percentage of inclined fibers and the range of inclination angles (half width at half maximum of inclination angle distribution). Our measurements showed that the fibers were not uniformly in-plane but exhibited instead a wide range of in-depth orientations, with fibers significantly more aligned in-plane in the anterior parts of the globe. We found that fitting the orientation distributions required between one and three π-periodic von Mises distributions with different primary orientations and fiber concentration factors. Regions of the posterior globe, particularly on the temporal side, had a larger percentage of inclined fibers and a larger range of inclination angles than anterior and equatorial regions. Variations of orientation distributions and anisotropies may imply varying out-of-plane tissue mechanical properties around the eye globe. Out-of-plane fibers could indicate fiber interweaving, not necessarily long, inclined fibers. Effects of small-scale fiber undulations, or crimp, were minimized by using tissues from eyes at high IOPs. These fiber features also play a role in tissue stiffness and stability and are therefore also important experimental information.

Posterior scleral birefringence measured by triple-input polarization-sensitive imaging as a biomarker of myopia progression.

In myopic eyes, pathological remodelling of collagen in the posterior sclera has mostly been observed ex vivo. Here we report the development of triple-input polarization-sensitive optical coherence tomography (OCT) for measuring posterior scleral birefringence. In guinea pigs and humans, the technique offers superior imaging sensitivities and accuracies than dual-input polarization-sensitive OCT. In 8-week-long studies with young guinea pigs, scleral birefringence was positively correlated with spherical equivalent refractive errors and predicted the onset of myopia. In a cross-sectional study involving adult individuals, scleral birefringence was associated with myopia status and negatively correlated with refractive errors. Triple-input polarization-sensitive OCT may help establish posterior scleral birefringence as a non-invasive biomarker for assessing the progression of myopia.

Direct measurements of collagen fiber recruitment in the posterior pole of the eye.

Collagen is the main load-bearing component of the peripapillary sclera (PPS) and lamina cribrosa (LC) in the eye. Whilst it has been shown that uncrimping and recruitment of the PPS and LC collagen fibers underlies the macro-scale nonlinear stiffening of both tissues with increased intraocular pressure (IOP), the uncrimping and recruitment as a function of local stretch have not been directly measured. This knowledge is crucial for the development of constitutive models associating micro and macro scales. In this project we measured local stretch-induced collagen fiber bundle uncrimping and recruitment curves of the PPS and LC. Thin coronal samples of PPS and LC of sheep eyes were mounted and stretched biaxially quasi-statically using a custom system. At each step, we imaged the PPS and LC with instant polarized light microscopy and quantified pixel-level (1.5 µm/pixel) collagen fiber orientations. We used digital image correlation to measure the local stretch and quantified collagen crimp by the circular standard deviation of fiber orientations, or waviness. Local stretch-recruitment curves of PPS and LC approximated sigmoid functions. PPS recruited more fibers than the LC at the low levels of stretch. At 10% stretch the curves crossed with 75% bundles recruited. The PPS had higher uncrimping rate and waviness remaining after recruitment than the LC: 0.9° vs. 0.6° and 3.1° vs. 2.7°. Altogether our findings support describing fiber recruitment of both PPS and LC with sigmoid curves, with the PPS recruiting faster and at lower stretch than the LC, consistent with a stiffer tissue.

Instant polarized light microscopy pi (IPOLπ) for quantitative imaging of collagen architecture and dynamics in ocular tissues.

Collagen architecture determines the biomechanical environment in the eye, and thus characterizing collagen fiber organization and biomechanics is essential to fully understand eye physiology and pathology. We recently introduced instant polarized light microscopy (IPOL) that encodes optically information about fiber orientation and retardance through a color snapshot. Although IPOL allows imaging collagen at the full acquisition speed of the camera, with excellent spatial and angular resolutions, a limitation is that the orientation-encoding color is cyclic every 90 degrees (π/2 radians). In consequence, two orthogonal fibers have the same color and therefore the same orientation when quantified by color-angle mapping. In this study, we demonstrate IPOLπ, a new variation of IPOL, in which the orientation-encoding color is cyclic every 180 degrees (π radians). Herein we present the fundamentals of IPOLπ, including a framework based on a Mueller-matrix formalism to characterize how fiber orientation and retardance determine the color. The improved quantitative capability of IPOLπ enables further study of essential biomechanical properties of collagen in ocular tissues, such as fiber anisotropy and crimp. We present a series of experimental calibrations and quantitative procedures to visualize and quantify ocular collagen orientation and microstructure in the optic nerve head, a region in the back of the eye. There are four important strengths of IPOLπ compared to IPOL. First, IPOLπ can distinguish the orientations of orthogonal collagen fibers via colors, whereas IPOL cannot. Second, IPOLπ requires a lower exposure time than IPOL, thus allowing faster imaging speed. Third, IPOLπ allows visualizing non-birefringent tissues and backgrounds from tissue absorption, whereas both appear dark in IPOL images. Fourth, IPOLπ is cheaper and less sensitive to imperfectly collimated light than IPOL. Altogether, the high spatial, angular, and temporal resolutions of IPOLπ enable a deeper insight into ocular biomechanics and eye physiology and pathology.

A direct fiber approach to model sclera collagen architecture and biomechanics.

Sclera collagen fiber microstructure and mechanical behavior are central to eye physiology and pathology. They are also complex, and are therefore often studied using modeling. Most models of sclera, however, have been built within a conventional continuum framework. In this framework, collagen fibers are incorporated as statistical distributions of fiber characteristics such as the orientation of a family of fibers. The conventional continuum approach, while proven successful for describing the macroscale behavior of the sclera, does not account for the sclera fibers are long, interwoven and interact with one another. Hence, by not considering these potentially crucial characteristics, the conventional approach has only a limited ability to capture and describe sclera structure and mechanics at smaller, fiber-level, scales. Recent advances in the tools for characterizing sclera microarchitecture and mechanics bring to the forefront the need to develop more advanced modeling techniques that can incorporate and take advantage of the newly available highly detailed information. Our goal was to create a new computational modeling approach that can represent the sclera fibrous microstructure more accurately than with the conventional continuum approach, while still capturing its macroscale behavior. In this manuscript we introduce the new modeling approach, that we call direct fiber modeling, in which the collagen architecture is built explicitly by long, continuous, interwoven fibers. The fibers are embedded in a continuum matrix representing the non-fibrous tissue components. We demonstrate the approach by doing direct fiber modeling of a rectangular patch of posterior sclera. The model integrated fiber orientations obtained by polarized light microscopy from coronal and sagittal cryosections of pig and sheep. The fibers were modeled using a Mooney-Rivlin model, and the matrix using a Neo-Hookean model. The fiber parameters were determined by inversely matching experimental equi-biaxial tensile data from the literature. After reconstruction, the direct fiber model orientations agreed well with the microscopy data both in the coronal plane (adjusted R2 = 0.8234) and in the sagittal plane (adjusted R2 = 0.8495) of the sclera. With the estimated fiber properties (C10 = 5746.9 MPa; C01 = -5002.6 MPa, matrix shear modulus 200 kPa), the model's stress-strain curves simultaneously fit the experimental data in radial and circumferential directions (adjusted R2's 0.9971 and 0.9508, respectively). The estimated fiber elastic modulus at 2.16% strain was 5.45 GPa, in reasonable agreement with the literature. During stretch, the model exhibited stresses and strains at sub-fiber level, with interactions among individual fibers which are not accounted for by the conventional continuum methods. Our results demonstrate that direct fiber models can simultaneously describe the macroscale mechanics and microarchitecture of the sclera, and therefore that the approach can provide unique insight into tissue behavior questions inaccessible with continuum approaches.

Predatory bacteria can reduce Pseudomonas aeruginosa induced corneal perforation and proliferation in a rabbit keratitis model.

PURPOSE: Pseudomonas aeruginosa keratitis is a severe ocular infection that can lead to perforation of the cornea. In this study we evaluated the role of bacterial quorum sensing in generating corneal perforation and bacterial proliferation and tested whether co-injection of the predatory bacteria Bdellovibrio bacteriovorus could alter the clinical outcome. P. aeruginosa with lasR mutations were observed among keratitis isolates from a study collecting samples from India, so an isogenic lasR mutant strain of P. aeruginosa was included. METHODS: Rabbit corneas were intracorneally infected with P. aeruginosa strain PA14 or an isogenic ΔlasR mutant and co-injected with PBS or B. bacteriovorus. After 24 h, eyes were evaluated for clinical signs of infection. Samples were analyzed by scanning electron microscopy, optical coherence tomography, sectioned for histology, and corneas were homogenized for CFU enumeration and for inflammatory cytokines. RESULTS: We observed that 54% of corneas infected by wild-type PA14 presented with a corneal perforation (n = 24), whereas only 4% of PA14 infected corneas that were co-infected with B. bacteriovorus perforate (n = 25). Wild-type P. aeruginosa proliferation was reduced 7-fold in the predatory bacteria treated eyes. The ΔlasR mutant was less able to proliferate compared to the wild-type, but was largely unaffected by B. bacteriovorus. CONCLUSION: These studies indicate a role for bacterial quorum sensing in the ability of P. aeruginosa to proliferate and cause perforation of the rabbit cornea. Additionally, this study suggests that predatory bacteria can reduce the virulence of P. aeruginosa in an ocular prophylaxis model.

Individual astrocyte morphology in the collagenous lamina cribrosa revealed by multicolor DiOlistic labeling.

Astrocytes in the lamina region of the optic nerve head play vital roles in supporting retinal ganglion cell axon health. In glaucoma, these astrocytes are implicated as early responders to stressors, undergoing characteristic changes in cell function as well as cell morphology. Much of what is currently known about individual lamina astrocyte morphology has been learned from rodent models which lack a defining feature of the human optic nerve head, the collagenous lamina cribrosa (LC). Current methods available for evaluation of collagenous LC astrocyte morphology have significant shortcomings. We aimed to evaluate Multicolor DiOlistic labeling (MuDi) as an approach to reveal individual astrocyte morphologies across the collagenous LC. Gold microcarriers were coated with all combinations of three fluorescent cell membrane dyes, DiI, DiD, and DiO, for a total of seven dye combinations. Microcarriers were delivered to 150 µm-thick coronal vibratome slices through the LC of pig, sheep, goat, and monkey eyes via MuDi. Labeled tissues were imaged with confocal and second harmonic generation microscopy to visualize dyed cells and LC collagenous beams, respectively. GFAP labeling of DiOlistically-labeled cells with astrocyte morphologies was used to investigate cell identity. 3D models of astrocytes were created from confocal image stacks for quantification of morphological features. DiOlistic labeling revealed fine details of LC astrocyte morphologies including somas, primary branches, higher-order branches, and end-feet. Labeled cells with astrocyte morphologies were GFAP+. Astrocytes were visible across seven distinct color channels, allowing high labeling density while still distinguishing individual cells from their neighbors. MuDi was capable of revealing tens to hundreds of collagenous LC astrocytes, in situ, with a single application. 3D astrocyte models allowed automated quantification of morphological features including branch number, length, thickness, hierarchy, and straightness as well as Sholl analysis. MuDi labeling provides an opportunity to investigate morphologies of collagenous LC astrocytes, providing both qualitative and quantitative detail, in healthy tissues. This approach may open doors for research of glaucoma, where astrocyte morphological alterations are thought to coincide with key functional changes related to disease progression.

Predatory Bacteria can Reduce Pseudomonas aeruginosa Induced Corneal Perforation and Proliferation in a Rabbit Keratitis Model.

Purpose: Pseudomonas aeruginosa keratitis is a severe ocular infection that can lead to perforation of the cornea. In this study we evaluated the role of bacterial quorum sensing in generating corneal perforation and bacterial proliferation and tested whether co-injection of the predatory bacteria Bdellovibrio bacteriovorus could alter the clinical outcome. P. aeruginosa with lasR mutations were observed among keratitis isolates from a study collecting samples from India, so an isogenic lasR mutant strain of P. aeruginosa was included. Methods: Rabbit corneas were intracorneally infected with P. aeruginosa strain PA14 or an isogenic Δ lasR mutant and co-injected with PBS or B. bacteriovorus . After 24 h, eyes were evaluated for clinical signs of infection. Samples were analyzed by scanning electron microscopy, optical coherence tomography, sectioned for histology, and corneas were homogenized for CFU enumeration and for inflammatory cytokines. Results: We observed that 54% of corneas infected by wild-type PA14 presented with a corneal perforation (n=24), whereas only 4% of PA14 infected corneas that were co-infected with B. bacteriovorus perforate (n=25). Wild-type P. aeruginosa proliferation was reduced 7-fold in the predatory bacteria treated eyes. The Δ lasR mutant was less able to proliferate compared to the wild-type, but was largely unaffected by B. bacteriovorus . Conclusion: These studies indicate a role for bacterial quorum sensing in the ability of P. aeruginosa to proliferate and cause perforation of the rabbit cornea. Additionally, this study suggests that predatory bacteria can reduce the virulence of P. aeruginosa in an ocular prophylaxis model.

Who bears the load? IOP-induced collagen fiber recruitment over the corneoscleral shell.

Collagen is the main load-bearing component of cornea and sclera. When stretched, both of these tissues exhibit a behavior known as collagen fiber recruitment. In recruitment, as the tissues stretch the constitutive collagen fibers lose their natural waviness, progressively straightening. Recruited, straight, fibers bear substantially more mechanical load than non-recruited, wavy, fibers. As such, the process of recruitment underlies the well-established nonlinear macroscopic behavior of the corneoscleral shell. Recruitment has an interesting implication: when recruitment is incomplete, only a fraction of the collagen fibers is actually contributing to bear the loads, with the rest remaining "in reserve". In other words, at a given intraocular pressure (IOP), it is possible that not all the collagen fibers of the cornea and sclera are actually contributing to bear the loads. To the best of our knowledge, the fraction of corneoscleral shell fibers recruited and contributing to bear the load of IOP has not been reported. Our goal was to obtain regionally-resolved estimates of the fraction of corneoscleral collagen fibers recruited and in reserve. We developed a fiber-based microstructural constitutive model that could account for collagen fiber undulations or crimp via their tortuosity. We used experimentally-measured collagen fiber crimp tortuosity distributions in human eyes to derive region-specific nonlinear hyperelastic mechanical properties. We then built a three-dimensional axisymmetric model of the globe, assigning region-specific mechanical properties and regional anisotropy. The model was used to simulate the IOP-induced shell deformation. The model-predicted tissue stretch was then used to quantify collagen recruitment within each shell region. The calculations showed that, at low IOPs, collagen fibers in the posterior equator were recruited the fastest, such that at a physiologic IOP of 15 mmHg, over 90% of fibers were recruited, compared with only a third in the cornea and the peripapillary sclera. The differences in recruitment between regions, in turn, mean that at a physiologic IOP the posterior equator had a fiber reserve of only 10%, whereas the cornea and peripapillary sclera had two thirds. At an elevated IOP of 50 mmHg, collagen fibers in the limbus and the anterior/posterior equator were almost fully recruited, compared with 90% in the cornea and the posterior sclera, and 70% in the peripapillary sclera and the equator. That even at such an elevated IOP not all the fibers were recruited suggests that there are likely other conditions that challenge the corneoscleral tissues even more than IOP. The fraction of fibers recruited may have other potential implications. For example, fibers that are not bearing loads may be more susceptible to enzymatic digestion or remodeling. Similarly, it may be possible to control tissue stiffness through the fraction of recruited fibers without the need to add or remove collagen.

Instant polarized light microscopy pi (IPOLπ) for quantitative imaging of collagen architecture and dynamics in ocular tissues.

Collagen architecture determines the biomechanical environment in the eye, and thus characterizing collagen fiber organization and biomechanics is essential to fully understand eye physiology and pathology. We recently introduced instant polarized light microscopy (IPOL) that encodes optically information about fiber orientation and retardance through a color snapshot. Although IPOL allows imaging collagen at the full acquisition speed of the camera, with excellent spatial and angular resolutions, a limitation is that the orientation-encoding color is cyclic every 90 degrees (π/2 radians). In consequence, two orthogonal fibers have the same color and therefore the same orientation when quantified by color-angle mapping. In this study, we demonstrate IPOLπ, a new variation of IPOL, in which the orientation-encoding color is cyclic every 180 degrees (π radians). Herein we present the fundamentals of IPOLπ, including a framework based on a Mueller-matrix formalism to characterize how fiber orientation and retardance determine the color. The improved quantitative capability of IPOLπ enables further study of essential biomechanical properties of collagen in ocular tissues, such as fiber anisotropy and crimp. We present a series of experimental calibrations and quantitative procedures to visualize and quantify ocular collagen orientation and microstructure in the optic nerve head, a region in the back of the eye. There are four important strengths of IPOLπ compared to IPOL. First, IPOLπ can distinguish the orientations of orthogonal collagen fibers via colors, whereas IPOL cannot. Second, IPOLπ requires a lower exposure time than IPOL, thus allowing faster imaging speed. Third, IPOLπ allows visualizing non-birefringent tissues and backgrounds from tissue absorption, whereas both appear dark in IPOL images. Fourth, IPOLπ is cheaper and less sensitive to imperfectly collimated light than IPOL. Altogether, the high spatial, angular, and temporal resolutions of IPOLπ enable a deeper insight into ocular biomechanics and eye physiology and pathology.