Quantification of collagen fiber structure using second harmonic generation imaging and two-dimensional discrete Fourier transform analysis: Application to the human optic nerve head

Collagen organization and structure in FBLN5-/- mice using label-free microscopy: implications for pelvic organ prolapse

Pelvic organ prolapse (POP) is a gynecological disorder described by the descent of superior pelvic organs into or out of the vagina as a consequence of disrupted muscles and tissue. A thorough understanding of the etiology of POP is limited by the availability of clinically relevant samples, restricting longitudinal POP studies on soft-tissue biomechanics and structure to POP-induced models such as fibulin-5 knockout (FBLN5-/- ) mice. Despite being a principal constituent in the extracellular matrix, little is known about structural perturbations to collagen networks in the FBLN5-/- mouse cervix. We identify significantly different collagen network populations in normal and prolapsed cervical cross-sections using two label-free, nonlinear microscopy techniques. Collagen in the prolapsed mouse cervix tends to be more isotropic, and displays reduced alignment persistence via 2-D Fourier transform analysis of images acquired using second harmonic generation microscopy. Furthermore, coherent Raman hyperspectral imaging revealed elevated disorder in the secondary structure of collagen in prolapsed tissues. Our results underscore the need for in situ multimodal monitoring of collagen organization to improve POP predictive capabilities.

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.

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.

Collagen organization and structure in FLBN5-/- mice using label-free microscopy: implications for pelvic organ prolapse

Pelvic organ prolapse (POP) is a gynecological disorder described by the descent of superior pelvic organs into or out of the vagina as a consequence of disrupted muscles and tissue. A thorough understanding of the etiology of POP is limited by the availability of clinically relevant samples, restricting longitudinal POP studies on soft-tissue biomechanics and structure to POP-induced models such as fibulin-5 knockout (FBLN5-/-) mice. Despite being a principal constituent in the extracellular matrix, little is known about structural perturbations to collagen networks in the FBLN5-/- mouse cervix. We identify significantly different collagen network populations in normal and prolapsed cervical cross-sections using two label-free, nonlinear microscopy techniques. Collagen in the prolapsed mouse cervix tends to be more isotropic, and displays reduced alignment persistence via 2-D Fourier Transform analysis of images acquired using second harmonic generation microscopy. Furthermore, coherent Raman hyperspectral imaging revealed elevated disorder in the secondary structure of collagen in prolapsed tissues. Our results underscore the need for in situ multimodal monitoring of collagen organization to improve POP predictive capabilities.

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.

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.

In situ measurement of the stiffness increase in the posterior sclera after UV-riboflavin crosslinking by optical coherence elastography

Scleral crosslinking may provide a way to prevent or treat myopia by stiffening scleral tissues. The ability to measure the stiffness of scleral tissues in situ pre and post scleral crosslinking would be useful but has not been established. Here, we tested the feasibility of optical coherence elastography (OCE) to measure shear modulus of scleral tissues and evaluate the impact of crosslinking on different posterior scleral regions using ex vivo porcine eyes as a model. From measured elastic wave speeds at 6 - 16 kHz, we obtained out-of-plane shear modulus value of 0.71 ± 0.12 MPa (n = 20) for normal porcine scleral tissues. After riboflavin-assisted UV crosslinking, the shear modulus increased to 1.50 ± 0.39 MPa (n = 20). This 2-fold change was consistent with the increase of static Young's modulus from 5.5 ± 1.1 MPa to 9.3 ± 1.9 MPa after crosslinking, which we measured using conventional uniaxial extensometry on tissue stripes. OCE revealed regional stiffness differences across the temporal, nasal, and deeper posterior sclera. Our results show the potential of OCE as a noninvasive tool to evaluate the effect of scleral crosslinking.

Remodelling of the fibre-aggregate structure of collagen gels by cancer-associated fibroblasts: A time-resolved grey-tone image analysis based on stochastic modelling

Introduction

Solid tumors consist of tumor cells associated with stromal and immune cells, secreted factors and extracellular matrix (ECM), which together constitute the tumor microenvironment. Among stromal cells, activated fibroblasts, known as cancer-associated fibroblasts (CAFs) are of particular interest. CAFs secrete a plethora of ECM components including collagen and modulate the architecture of the ECM, thereby influencing cancer cell migration. The characterization of the collagen fibre network and its space and time-dependent microstructural modifications is key to investigating the interactions between cells and the ECM. Developing image analysis tools for that purpose is still a challenge because the structural complexity of the collagen network calls for specific statistical descriptors. Moreover, the low signal-to-noise ratio of imaging techniques available for time-resolved studies rules out standard methods based on image segmentation.

Methods

In this work, we develop a novel approach based on the stochastic modelling of the gel structure and on grey-tone image analysis. The method is then used to study the remodelling of a collagen matrix by migrating breast cancer-derived CAFs in a three-dimensional spheroid model of cellular invasion imaged by time-lapse confocal microscopy.

Results

The structure of the collagen at the scale of a few microns consists in regions with high fibre density separated by depleted regions, which can be thought of as aggregates and pores. The approach developped captures this two-scale structure with a clipped Gaussian field model to describe the aggregates-and-pores large-scale structure, and a homogeneous Boolean model to describe the small-scale fibre network within the aggregates. The model parameters are identified by fitting the grey-tone histograms and correlation functions of the images. The method applies to unprocessed grey-tone images, and it can therefore be used with low magnification, noisy time-lapse reflectance images. When applied to the CAF spheroid time-resolved images, the method reveals different matrix densification mechanisms for the matrix in direct contact or far from the cells.

Conclusion

We developed a novel and multidisciplinary image analysis approach to investigate the remodelling of fibrillar collagen in a 3D spheroid model of cellular invasion. The specificity of the method is that it applies to the unprocessed grey-tone images, and it can therefore be used with noisy time-lapse reflectance images of non-fluorescent collagen. When applied to the CAF spheroid time-resolved images, the method reveals different matrix densification mechanisms for the matrix in direct contact or far from the cells.

Integrated Cells and Collagen Fibers Spatial Image Analysis

Modern technologies designed for tissue structure visualization like brightfield microscopy, fluorescent microscopy, mass cytometry imaging (MCI) and mass spectrometry imaging (MSI) provide large amounts of quantitative and spatial information about cells and tissue structures like vessels, bronchioles etc. Many published reports have demonstrated that the structural features of cells and extracellular matrix (ECM) and their interactions strongly predict disease development and progression. Computational image analysis methods in combination with spatial analysis and machine learning can reveal novel structural patterns in normal and diseased tissue. Here, we have developed a Python package designed for integrated analysis of cells and ECM in a spatially dependent manner. The package performs segmentation, labeling and feature analysis of ECM fibers, combines this information with pre-generated single-cell based datasets and realizes cell-cell and cell-fiber spatial analysis. To demonstrate performance and compatibility of our computational tool, we integrated it with a pipeline designed for cell segmentation, classification, and feature analysis in the KNIME analytical platform. For validation, we used a set of mouse mammary gland tumors and human lung adenocarcinoma tissue samples stained for multiple cellular markers and collagen as the main ECM protein. The developed package provides sufficient performance and precision to be used as a novel method to investigate cell-ECM relationships in the tissue, as well as detect structural patterns correlated with specific disease outcomes.

Study of collagen remodeling in structural disorders of the temporomandibular joint using second-harmonic generation microscopy

Structural disorder of the temporomandibular joint (TMJ) is a progressive disease with poor prognosis due to its physiological three-dimensional anatomical position and the complicated relationship among muscles, ligaments, and cartilage. The lack of detection methods for changes in the collagen structure of the TMJ disc makes the diagnosis untimely and unclear. This work aimed to explore the feasibility of using a promising detection technique, second-harmonic generation (SHG), to characterize collagen fibers in a TMJ disc with structural disorders. The TMJ discs with structural disorder were observed using SHG microscopy, and assessment of collagen orientation was conducted by analyzing digitized images. The SHG images were also compared with the scanning electron microscopy images and microscopic images acquired after hematoxylin and eosin and Masson's trichrome staining. The SHG imaging showed that the collagen fibers in diseased TMJ were distributed in a disorderly manner, and pixel intensities in diseased TMJ discs were significantly different from those acquired in healthy TMJs. Moreover, the three dimensions of collagen fibers and dynamic images acquired by SHG nonlinear optical microscopy showed the structural disorder of the collagen fibers in a diseased TMJ. In summary, SHG imaging could provide three-dimensional and quantitative data, with dynamic and critical pathological information for clinical diagnosis, showing its potential value in the diagnosis and evaluation of structural disorders of the TMJ disc.

Mechanosensitive channel inhibition attenuates TGFβ2-induced actin cytoskeletal remodeling and reactivity in mouse optic nerve head astrocytes

Astrocytes within the optic nerve head undergo actin cytoskeletal rearrangement early in glaucoma, which coincides with astrocyte reactivity and extracellular matrix (ECM) deposition. Elevated transforming growth factor beta 2 (TGFβ2) levels within astrocytes have been described in glaucoma, and TGFβ signaling induces actin cytoskeletal remodeling and ECM deposition in many tissues. A key mechanism by which astrocytes sense and respond to external stimuli is via mechanosensitive ion channels. Here, we tested the hypothesis that inhibition of mechanosensitive channels will attenuate TGFβ2-mediated optic nerve head astrocyte actin cytoskeletal remodeling, reactivity, and ECM deposition. Primary optic nerve head astrocytes were isolated from C57BL/6J mice and cell purity was confirmed by immunostaining. Astrocytes were treated with vehicle control, TGFβ2 (5 ng/ml), GsMTx4 (a mechanosensitive channel inhibitor; 500 nM), or TGFβ2 (5 ng/ml) + GsMTx4 (500 nM) for 48 h. FITC-phalloidin staining was used to assess the formation of f-actin stress fibers and to quantify the presence of crosslinked actin networks (CLANs). Cell reactivity was determined by immunostaining and immunoblotting for GFAP. Levels of fibronectin and collagen IV deposition were also quantified. Primary optic nerve head astrocytes were positive for the astrocyte marker GFAP and negative for markers for microglia (F4/80) and oligodendrocytes (OSP1). Significantly increased %CLAN-positive cells were observed after 48-h treatment with TGFβ2 vs. control in a dose-dependent manner. Co-treatment with GsMTx4 significantly decreased %CLAN-positive cells vs. TGFβ2 treatment and the presence of f-actin stress fibers. TGFβ2 treatment significantly increased GFAP, fibronectin, and collagen IV levels, and GsMTx4 co-treatment ameliorated GFAP immunoreactivity. Our data suggest inhibition of mechanosensitive channel activity as a potential therapeutic strategy to modulate actin cytoskeletal remodeling within the optic nerve head in glaucoma.

Characterization of Collagen I Fiber Thickness, Density, and Orientation in the Human Skin In Vivo Using Second-Harmonic Generation Imaging

The assessment of dermal alterations is necessary to monitor skin aging, cancer, and other skin diseases and alterations. The gold standard of morphologic diagnostics is still histopathology. Here, we proposed parameters to distinguish morphologically different collagen I structures in the extracellular matrix and to characterize varying collagen I structures in the skin with similar SAAID (SHG-to-AF Aging Index of Dermis, SHG—second-harmonic generation; AF—autofluorescence) values. Test datasets for the papillary and reticular extracellular matrix from images in 24 female subjects, 36 to 50 years of age, were generated. Parameters for SAAID, edge detection, and fast Fourier transformation directionality were determined. Additionally, textural analyses based on the grey level co-occurrence matrix (GLCM) were conducted. At first, changes in the GLCM parameters were determined in the native greyscale images and, furthermore, in the Hilbert-transformed images. Our results demonstrate a robust set of parameters for noninvasive in vivo classification for morphologically different collagen I structures in the skin, with similar and different SAAID values. We anticipate our method to enable an automated prevention and monitoring system with an age- and gender-specific algorithm.

Decorin-An Antagonist of TGF-β in Astrocytes of the Optic Nerve

During the pathogenesis of glaucoma, optic nerve (ON) axons become continuously damaged at the optic nerve head (ONH). This often is associated with reactive astrocytes and increased transforming growth factor (TGF-β) 2 levels. In this study we tested the hypothesis if the presence or absence of decorin (DCN), a small leucine-rich proteoglycan and a natural inhibitor of several members of the TGF family, would affect the expression of the TGF-βs and connective tissue growth factor (CTGF/CCN2) in human ONH astrocytes and murine ON astrocytes. We found that DCN is present in the mouse ON and is expressed by human ONH and murine ON astrocytes. DCN expression and synthesis was significantly reduced after 24 h treatment with 3 nM CTGF/CCN2, while treatment with 4 pM TGF-β2 only reduced expression of DCN significantly. Conversely, DCN treatment significantly reduced the expression of TGF-β1, TGF-β2 and CTGF/CCN2 vis-a-vis untreated controls. Furthermore, DCN treatment significantly reduced expression of fibronectin (FN) and collagen IV (COL IV). Notably, combined treatment with DCN and triciribine, a small molecule inhibitor of protein kinase B (AKT), attenuated effects of DCN on CTGF/CCN2, TGF-β1, and TGF-β2 mRNA expression. We conclude (1) that DCN is an important regulator of TGF-β and CTGF/CCN2 expression in astrocytes of the ON and ONH, (2) that DCN thereby regulates the expression of extracellular matrix (ECM) components and (3) that DCN executes its negative regulatory effects on TGF-β and CTGF/CCN2 via the pAKT/AKT signaling pathway in ON astrocytes.

A Quantitative Analysis of Bone Lamellarity and Bone Collagen Linearity Induced by Distinct Dosing and Frequencies of Teriparatide Administration in Ovariectomized Rats and Monkeys

The lamellar structure of bone, which endows biomechanical rigidity to support the host organism, is observed in mammals, including humans. It is therefore essential to develop a quantitative analysis to evaluate the lamellarity of bone, which would especially be useful for the pharmacological evaluation of anti-osteoporotic drugs. This study applied a current system for the semi-automatic recognition of fluorescence signals to the analysis of un-decalcified bone sections from rat and monkey specimens treated with teriparatide (TPTD). Our analyses on bone formation pattern and collagen topology indicated that TPTD augmented bone lamellarity and bone collagen linearity, which were possibly associated with the recovery of collagen cross-linking, thus endowing bone rigidity.

Substrate-dependent control of the chiral orientation of mesenchymal stem cells: image-based quantitative profiling

Stem-cell (SC) chirality or left-right (LR) asymmetry is an essential attribute, observed during tissue regeneration. The ability to control the LR orientation of cells by biophysical manipulation is a promising approach for recapitulating their inherent function. Despite remarkable progress in tissue engineering, the development of LR chirality in SCs has been largely unexplored. Here, we demonstrate the role of substrate stiffness on the LR asymmetry of cultured mesenchymal stem cells (MSCs). We found that MSCs acquired higher asymmetricity when cultured on stiffer PCL/collagen matrices. To confirm cellular asymmetry, different parameters such as the aspect ratio, orientation angle and intensity of polarized proteins (Par) were investigated. The results showed a significant (p < 0.01) difference in the average orientation angle, the cellular aspect ratio, and the expression of actin and Par proteins in MSCs cultured on matrices with different stiffnesses. Furthermore, a Gaussian support-vector machine was applied to classify cells cultured on both (2% and 10% PCL/Collagen) matrices, with a resulting accuracy of 96.2%. To the best of our knowledge, this study is the first that interrelates and quantifies MSC asymmetricity with matrix properties using a simple 2D model.

Spatial frequency metrics for analysis of microscopic images of musculoskeletal tissues

Purpose: Imaging-based metrics for analysis of biological tissues are powerful tools that can extract information such as shape, size, periodicity, and many other features to assess the requested qualities of a tissue. Muscular and osseous tissues consist of periodic structures that are directly related to their function, and so analysis of these patterns likely reflects tissue health and regeneration.Methods: A method for assessment of periodic structures is by analyzing them in the spatial frequency domain using the Fourier transform. In this paper, we present two filters which we developed in the spatial frequency domain for the purpose of analyzing musculoskeletal structures. These filters provide information about 1) the angular orientation of the tissues and 2) their periodicity. We explore periodic structural patterns in the mitochondrial network of skeletal muscles that are reflective of muscle metabolism and myogenesis; and patterns of collagen fibers in the bone that are reflective of the organization and health of bone extracellular matrix.Results: We present an analysis of mouse skeletal muscle in healthy and injured muscles. We used a transgenic mouse that ubiquitously expresses fluorescent protein in their mitochondria and performed 2-photon microscopy to image the structures. To acquire the collagen structure of the bone we used non-linear SHG microscopy of mouse flat bone. We analyze and compare juvenile versus adult mice, which have different structural patterns.Conclusions: Our results indicate that these metrics can quantify musculoskeletal tissues during development and regeneration.

The Effects of Glaucoma on the Pressure-Induced Strain Response of the Human Lamina Cribrosa

Purpose

To measure the ex vivo pressure-induced strain response of the human optic nerve head and analyze for variations with glaucoma diagnosis and optic nerve axon damage.

Methods

The posterior sclera of 16 eyes from 8 diagnosed glaucoma donors and 10 eyes from 6 donors with no history of glaucoma were inflation tested between 5 and 45 mm Hg. The optic nerve from each donor was examined for degree of axon loss. The posterior volume of the lamina cribrosa (LC) was imaged with second harmonic generation and analyzed using volume correlation to calculate LC strains between 5 and 10 and 5 and 45 mm Hg.

Results

Eye length and LC area were larger in eyes diagnosed with glaucoma (P≤ 0.03). Nasal-temporal E_XX and circumferential E_θθ strains were lower in the LC of diagnosed glaucoma eyes at 10 mm Hg (P≤ 0.05) and 45 mm Hg (P≤ 0.07). E_XX was smaller in the LC of glaucoma eyes with <25% axon loss compared with undamaged normal eyes (P = 0.01, 45 mm Hg). In general, the strains were larger in the peripheral than central LC. The ratio of the maximum principal strain E_max in the peripheral to central LC was larger in glaucoma eyes with >25% axon loss than in glaucoma eyes with milder damage (P = 0.004, 10 mm Hg).

Conclusions

The stiffness of the LC pressure-strain response was greater in diagnosed glaucoma eyes and varied with glaucomatous axon damage. Lower LC strains in glaucoma eyes with milder damage may represent baseline biomechanical behavior that contributes to axon loss, whereas greater LC strain and altered radial LC strain variation in glaucoma eyes with more severe damage may be caused by glaucoma-related remodeling.

Astrocyte responses to experimental glaucoma in mouse optic nerve head

PURPOSE

To delineate responses of optic nerve head astrocytes to sustained intraocular pressure (IOP) elevation in mice.

METHODS

We elevated IOP for 1 day to 6 weeks by intracameral microbead injection in 4 strains of mice. Astrocyte alterations were studied by transmission electron microscopy (TEM) including immunogold molecular localization, and by laser scanning microscopy (LSM) with immunofluorescence for integrin β1, α-dystroglycan, and glial fibrillary acidic protein (GFAP). Astrocyte proliferation and apoptosis were quantified by Ki67 and TUNEL labeling, respectively.

RESULTS

Astrocytes in normal optic nerve head expressed integrin β1 and α-dystroglycan by LSM and TEM immunogold labeling at electron dense junctional complexes that were found only on cell membrane zones bordering their basement membranes (BM) at the peripapillary sclera (PPS) and optic nerve head capillaries. At 1-3 days after IOP elevation, abnormal extracellular spaces appeared between astrocytes near PPS, and axonal vesical and mitochondrial accumulation indicated axonal transport blockade. By 1 week, abnormal spaces increased, new collagen formation occurred, and astrocytes separated from their BM, leaving cell membrane fragments. Electron dense junctional complexes separated or were absent at the BM. Astrocyte proliferation was modest during the first week, while only occasional apoptotic astrocytes were observed by TEM and TUNEL.

CONCLUSIONS

Astrocytes normally exhibit junctions with their BM which are disrupted by extended IOP elevation. Responses include reorientation of cell processes, new collagen formation, and cell proliferation.

Biomechanical properties of the rat sclera obtained with inverse finite element modeling

It is widely accepted that biomechanics plays an important role in glaucoma pathophysiology, but the mechanisms involved are largely unknown. Rats are a common animal model of glaucoma, and finite element models are being developed to provide much-needed insight into the biomechanical environment of the posterior rat eye. However, material properties of rat ocular tissues, including the sclera, are currently unknown. Since the sclera plays a major role in posterior ocular biomechanics, our goal was to use inverse finite element modeling to extract rat scleral material properties. We first used digital image correlation to measure scleral surface displacement during whole-globe inflation testing. We modeled the sclera as a nonlinear material with embedded collagen fibers and then fit modeled displacements to experimental data using a differential evolution algorithm. Subject-specific models were constructed in which 3 parameters described the stiffness of the ground substance and collagen fibers in the posterior eye, and 16 parameters defined the primary orientation and alignment of fibers within eight scleral sub-regions. We successfully extracted scleral material properties for eight rat eyes. Model displacements recreated general patterns of the experimental displacements but did not always match local patterns. The fiber directions and fiber concentration parameters were highly variable, but on average, fibers were aligned circumferentially and were more aligned in the peripapillary sclera than in the peripheral sclera. The material properties determined here will be used to inform future finite element models of the rat posterior eye with the goal of elucidating the role of biomechanics in glaucoma pathophysiology.

Fibrillar Collagen Quantification With Curvelet Transform Based Computational Methods

Quantification of fibrillar collagen organization has given new insight into the possible role of collagen topology in many diseases and has also identified candidate image-based bio-markers in breast cancer and pancreatic cancer. We have been developing collagen quantification tools based on the curvelet transform (CT) algorithm and have demonstrated this to be a powerful multiscale image representation method due to its unique features in collagen image denoising and fiber edge enhancement. In this paper, we present our CT-based collagen quantification software platform with a focus on new features and also giving a detailed description of curvelet-based fiber representation. These new features include C++-based code optimization for fast individual fiber tracking, Java-based synthetic fiber generator module for method validation, automatic tumor boundary generation for fiber relative quantification, parallel computing for large-scale batch mode processing, region-of-interest analysis for user-specified quantification, and pre- and post-processing modules for individual fiber visualization. We present a validation of the tracking of individual fibers and fiber orientations by using synthesized fibers generated by the synthetic fiber generator. In addition, we provide a comparison of the fiber orientation calculation on pancreatic tissue images between our tool and three other quantitative approaches. Lastly, we demonstrate the use of our software tool for the automatic tumor boundary creation and the relative alignment quantification of collagen fibers in human breast cancer pathology images, as well as the alignment quantification of in vivo mouse xenograft breast cancer images.

Factors affecting optic nerve head biomechanics in a rat model of glaucoma

Glaucoma is the leading cause of irreversible blindness and is characterized by the death of retinal ganglion cells, which carry vision information from the retina to the brain. Although it is well accepted that biomechanics is an important part of the glaucomatous disease process, the mechanisms by which biomechanical insult, usually due to elevated intraocular pressure (IOP), leads to retinal ganglion cell death are not understood. Rat models of glaucoma afford an opportunity for learning more about these mechanisms, but the biomechanics of the rat optic nerve head (ONH), a primary region of damage in glaucoma, are only just beginning to be characterized. In a previous study, we built finite-element models with individual-specific rat ONH geometries. Here, we developed a parametrized model of the rat ONH and used it to perform a sensitivity study to determine the influence that six geometric parameters and 13 tissue material properties have on rat optic nerve biomechanical strains due to IOP elevation. Strain magnitudes and patterns in the parametrized model generally matched those from individual-specific models, suggesting that the parametrized model sufficiently approximated rat ONH anatomy. Similar to previous studies in human eyes, we found that scleral properties were highly influential: the six parameters with highest influence on optic nerve strains were optic nerve stiffness, IOP, scleral thickness, the degree of alignment of scleral collagen fibres, scleral ground substance stiffness and the scleral collagen fibre uncrimping coefficient. We conclude that a parametrized modelling strategy is an efficient approach that allows insight into rat ONH biomechanics. Further, scleral properties are important influences on rat ONH biomechanics, and additional efforts should be made to better characterize rat scleral collagen fibre organization.

Characterizing the Collagen Network Structure and Pressure-Induced Strains of the Human Lamina Cribrosa

Purpose

The purpose of this study was to measure the 2D collagen network structure of the human lamina cribrosa (LC), analyze for the correlations with age, region, and LC size, as well as the correlations with pressure-induced strains.

Methods

The posterior scleral cups of 10 enucleated human eyes with no known ocular disease were subjected to ex vivo inflation testing from 5 to 45 mm Hg. The optic nerve head was imaged by using second harmonic generation imaging (SHG) to identify the LC collagen structure at both pressures. Displacements and strains were calculated by using digital volume correlation of the SHG volumes. Nine structural features were measured by using a custom Matlab image analysis program, including the pore area fraction, node density, and beam connectivity, tortuosity, and anisotropy.

Results

All strain measures increased significantly with higher pore area fraction, and all but the radial-circumferential shear strain (E_rθ) decreased with higher node density. The maximum principal strain (E_max) and maximum shear strain (Γ_max) also increased with larger beam aspect ratio and tortuosity, respectively, and decreased with higher connectivity. The peripheral regions had lower node density and connectivity, and higher pore area fraction, tortuosity, and strains (except for E_rθ) than the central regions. The peripheral nasal region had the lowest E_max, Γ_max, radial strain, and pore area fraction.

Conclusions

Features of LC beam network microstructure that are indicative of greater collagen density and connectivity are associated with lower pressure-induced LC strain, potentially contributing to resistance to glaucomatous damage.

Scleral structure and biomechanics

As the eye's main load-bearing connective tissue, the sclera is centrally important to vision. In addition to cooperatively maintaining refractive status with the cornea, the sclera must also provide stable mechanical support to vulnerable internal ocular structures such as the retina and optic nerve head. Moreover, it must achieve this under complex, dynamic loading conditions imposed by eye movements and fluid pressures. Recent years have seen significant advances in our knowledge of scleral biomechanics, its modulation with ageing and disease, and their relationship to the hierarchical structure of the collagen-rich scleral extracellular matrix (ECM) and its resident cells. This review focuses on notable recent structural and biomechanical studies, setting their findings in the context of the wider scleral literature. It reviews recent progress in the development of scattering and bioimaging methods to resolve scleral ECM structure at multiple scales. In vivo and ex vivo experimental methods to characterise scleral biomechanics are explored, along with computational techniques that combine structural and biomechanical data to simulate ocular behaviour and extract tissue material properties. Studies into alterations of scleral structure and biomechanics in myopia and glaucoma are presented, and their results reconciled with associated findings on changes in the ageing eye. Finally, new developments in scleral surgery and emerging minimally invasive therapies are highlighted that could offer new hope in the fight against escalating scleral-related vision disorder worldwide.