Quantization of collagen organization in the stroma with a new order coefficient

Optimizing near-infrared polariscopic imaging for the living human eye

Hardware architectures and image interpretation can be simplified by partial polarimetry. Mueller matrix (MM) polarimetry allows the investigation of partial polarimeter designs for a given scientific task. In this work, we use MM measurements to solve for a fixed polarization illumination and analyzer state that maximize polariscopic image contrast of the human eye. The eye MM image acquisition takes place over 15 seconds which motivates the development of a partial polarimeter that has snapshot operation. Within the eye, the birefringent cornea produces spatially-varying patterns of retardance exceeding half of a wave with a fast-axis varying from linear, to circular, and elliptical states in between. Our closed-form polariscopic pairs are a general solution that maximizes contrast between two non-depolarizing pure retarder MMs. For these MMs, there is a family of polariscopic pairs that maximize contrast. This range of solutions creates an opportunity to use the distance from optimal as a criteria to adjust polarimetric hardware architecture. We demonstrate our optimization approach by performing both Mueller and polariscopic imaging of an in vivo human eye at 947 nm using a dual-rotating-retarder polarimeter. Polariscopic images are simulated from Mueller measurements of 19 other human subjects to test the robustness of this optimal solution.

A review of the epithelial and stromal effects of corneal collagen crosslinking

Induced corneal collagen crosslinking and mechanical stiffening via ultraviolet-A photoactivation of riboflavin (UVA CXL) is now a common treatment for corneal ectasia and Keratoconus. Some effects of the procedure such as induced mechanical stiffening, corneal flattening, and cellular toxicity are well-known, but others remain more controversial. Authors report a variety of contradictory effects, and provide evidence based on individual results and observations. A full understanding of the effects of and mechanisms behind this procedure are essential to predicting its outcome. A growing interest in modifications to the standard UVA CXL protocol, such as transepithelial or accelerated UVA CXL, makes analyzing the literature as a whole more urgent. This review presents an analysis of both the agreed-upon and contradictory results reported and the various methods used to obtain them.

Three-dimensional analysis of load-dependent changes in the orientation of dermal collagen fibers in human skin: A pilot study

The availability of quantitative structural data on the orientation of collagen fibers is of crucial importance for understanding the behavior of connective tissues. These fibers can be visualized using a variety of imaging techniques, including second harmonic generation (SHG) microscopy. However, characterization of the collagen network requires the accurate extraction of parameters from imaging data. To this end, several automated processes have been developed to identify the preferred orientation of collagen fibers. Common methods include fast Fourier transforms and curvelet transforms, but these tools are mostly used to infer a single preferred orientation. The purpose of this pilot study was to develop an easy procedure for comprehensively comparing multiple human skin samples with the goal of analyzing load-dependent changes via SHG microscopy. We created a 3D model based upon 2D image stacks that provide fiber orientation data perpendicular and parallel to the plane of the epidermis. The SHG images were analyzed by CurveAlign to obtain angle histogram plots containing information about the multiple fiber orientations in each single image. Subsequently, contour plots of the angle histogram intensities were created to provide a useful visual plotting method to clearly show the anomalies in the angle histograms in all samples. Our results provided additional details on how the collagen network carries a load. In fact, analysis of SHG images indicated that increased stretch was accompanied by an increase in the alignment of fibers in the loading direction. Moreover, these images demonstrated that more than one type of preferred orientation is present. In particular, the 3D network of fibers appears to have two preferred orientations in the planes both perpendicular and parallel to the plane of the epidermis.

Estimation of scleral mechanical properties from air-puff optical coherence tomography

We introduce a method to estimate the biomechanical properties of the porcine sclera in intact eye globes ex vivo, using optical coherence tomography that is coupled with an air-puff excitation source, and inverse optimization techniques based on finite element modeling. Air-puff induced tissue deformation was determined at seven different locations on the ocular globe, and the maximum apex deformation, the deformation velocity, and the arc-length during deformation were quantified. In the sclera, the experimental maximum deformation amplitude and the corresponding arc length were dependent on the location of air-puff excitation. The normalized temporal deformation profile of the sclera was distinct from that in the cornea, but similar in all tested scleral locations, suggesting that this profile is independent of variations in scleral thickness. Inverse optimization techniques showed that the estimated scleral elastic modulus ranged from 1.84 ± 0.30 MPa (equatorial inferior) to 6.04 ± 2.11 MPa (equatorial temporal). The use of scleral air-puff imaging holds promise for non-invasively investigating the structural changes in the sclera associated with myopia and glaucoma, and for monitoring potential modulation of scleral stiffness in disease or treatment.

Quantitative Classification of 3D Collagen Fiber Organization From Volumetric Images

Collagen fibers in biological tissues have a complex 3D organization containing rich information linked to tissue mechanical properties and are affected by mutations that lead to diseases. Quantitative assessment of this 3D collagen fiber organization could help to develop reliable biomechanical models and understand tissue structure-function relationships, which impact diagnosis and treatment of diseases or injuries. While there are advanced techniques for imaging collagen fibers, published methods for quantifying 3D collagen fiber organization have been sparse and give limited structural information which cannot distinguish a wide range of 3D organizations. In this article, we demonstrate an algorithm for quantitative classification of 3D collagen fiber organization. The algorithm first simulates five groups, or classifications, of fiber organization: unidirectional, crimped, disordered, two-fiber family, and helical. These five groups are widespread in natural tissues and are known to affect the tissue's mechanical properties. We use quantitative metrics based on features such as preferred 3D fiber orientation and spherical variance to differentiate each classification in a repeatable manner. We validate our algorithm by applying it to second-harmonic generation images of collagen fibers in tendon and cervix tissue that has been sectioned in specified orientations, and we find strong agreement between classification from simulated data and the physical fiber organization. Our approach provides insight for interpreting 3D fiber organization directly from volumetric images. This algorithm could be applied to other fiber-like structures that are not necessarily made of collagen.

Three-dimensional collagen fiber mapping and tractography of human uterine tissue using OCT

Automatic quantification and visualization of 3-D collagen fiber architecture using Optical Coherence Tomography (OCT) has previously relied on polarization information and/or prior knowledge of tissue-specific fiber architecture. This study explores image processing, enhancement, segmentation, and detection algorithms to map 3-D collagen fiber architecture from OCT images alone. 3-D fiber mapping, histogram analysis, and 3-D tractography revealed fiber groupings and macro-organization previously unseen in uterine tissue samples. We applied our method on centimeter-scale mosaic OCT volumes of uterine tissue blocks from pregnant and non-pregnant specimens revealing a complex, patient-specific network of fibrous collagen and myocyte bundles.

Corneal Collagen Ordering After In Vivo Rose Bengal and Riboflavin Cross-Linking

Purpose

Photoactivated cornea collagen cross-linking (CXL) increases corneal stiffness by initiating formation of covalent bonds between stromal proteins. Because CXL depends on diffusion to distribute the photoinitiator, a gradient of CXL efficiency with depth is expected that may affect the degree of stromal collagen organization. We used second harmonic generation (SHG) microscopy to investigate the differences in stromal collagen organization in rabbit eyes after corneal CXL in vivo as a function of depth and time after surgery.

Methods

Rabbit corneas were treated in vivo with either riboflavin/UV radiation (UVX) or Rose Bengal/green light (RGX) and evaluated 1 and 2 months after CXL. Collagen fibers were imaged with a custom-built SHG scanning microscope through the central cornea (350 µm depth, 225 × 225 µm en face images). The order coefficient (OC), a metric for collagen organization, and total SHG signal were computed for each depth and compared between treatments.

Results

OC values of CXL-treated corneas were larger than untreated corneas by 27% and 20% after 1 month and 38% and 33% after 2 months for the RGX and UVX, respectively. RGX OC values were larger than UVX OC values by 3% and 5% at 1 and 2 months. The SHG signal was higher in CXL corneas than untreated corneas, both at 1 and 2 months after surgery, by 18% and 26% and 1% and 10% for RGX and UVX, respectively.

Conclusions

Increased OC corresponded with increased collagen fiber organization in CXL corneas. Changes in collagen organization parallel reported temporal changes in cornea stiffness after CXL and also, surprisingly, are detected deeper in the stroma than the regions stiffened by collagen cross-links.

Correlations between corneal biomechanics and specular microscopy in patient with cataract

This study aimed to analyze the connection between corneal biomechanics (corneal hysteresis, CH) and endothelial cell density of cornea (mean endothelial cell density, MCD) in patients diagnosed with cataract. This retrospective, observational study was performed in the Ophthalmology Clinic of the University Emergency Hospital in Bucharest. Of 60 patients (120 eyes) with cataract, who were included in this study, we analyzed the CH values obtained using with the Ocular Response Analyzer (ORA) and the MCD values obtained using the specular microscopy. The study groups comprised both men and women with ages ranging from 45 to 63 years. Patients were divided into three study groups according to CH values. In each batch, the CH values obtained with the Ocular Response Analyzer (ORA) were correlated with age, gender and MCD, then the subgroups were compared. All the data gathered showed no correlation to be statistically significant regarding the biomechanical properties of the cornea and the corneal endothelial cell density in patients with cataract.

Quantitative Discrimination of Healthy and Diseased Corneas With Second Harmonic Generation Microscopy

Purpose

To analyze the spatial organization of pathological corneas with second harmonic generation (SHG) imaging and to provide a proof of concept to objectively distinguish these from the healthy corneas.

Methods

A custom-built SHG microscope was used to image the anterior stroma of ex vivo corneas, both control and affected by some representative pathologies. The structure tensor (ST) was employed as a metric to explore and quantify the alterations in the spatial distribution of the collagen lamellae.

Results

The collagen arrangement differed between healthy and pathological samples. The former showed a regular distribution and a low structural dispersion (SD < 40°) within the stroma with a well-defined dominant orientation. This regular arrangement drastically turns into a disorganized pattern in pathological corneas (SD > 40°).

Conclusions

The combination of SHG imaging and the ST allows obtaining quantitative information to differentiate the stromal collagen organization in healthy and diseased corneas. This approach represents a feasible and powerful technique with potential applications in clinical corneal diagnoses.

Translational Relevance

The ST applied to SHG microscopy images of the corneal stroma provides an experimental objective score to differentiate control from pathological or damaged corneas. Future implementations of this technique in clinical environments might might be a promising tool in Ophthalmology, not only to diagnose and monitor corneal diseases, but also to follow-up surgical outcome.

High-speed collagen fiber modeling and orientation quantification for optical coherence tomography imaging

Quantifying collagen fiber architecture has clinical and scientific relevance across a variety of tissue types and adds functionality to otherwise largely qualitative imaging modalities. Optical coherence tomography (OCT) is uniquely suited for this task due to its ability to capture the collagen microstructure over larger fields of view than traditional microscopy. Existing image processing techniques for quantifying fiber architecture, while accurate and effective, are very slow for processing large datasets and tend to lack structural specificity. We describe here a computationally efficient method for quantifying and visualizing collagen fiber organization. The algorithm is demonstrated on swine atria, bovine anterior cruciate ligament, and human cervical tissue samples. Additionally, we show an improved performance for images with crimped fiber textures and low signal to noise when compared to similar methods.

Quantitative Analysis of the Corneal Collagen Distribution after In Vivo Cross-Linking with Second Harmonic Microscopy

Corneal cross-linking (CXL) is a surgical procedure able to modify corneal biomechanics and stabilize keratoconus progression. Although it is known that CXL produces changes in corneal collagen distribution, these are still a topic of discussion. Here we quantitatively compare the corneal stroma architecture between two animal models four weeks after in vivo conventional CXL treatment, with second harmonic generation (SHG) imaging microscopy and the structure tensor (ST). The healing stage and the stroma recovery were also analyzed by means of histological sections. Results show that the CXL effects depend on the initial arrangement of the corneal collagen. While the treatment increases the order in corneas with a low level of initial organization, corneas presenting a fairly regular pattern are hardly affected. Histological samples showed active keratocytes in anterior and middle stroma, what means that the recovery is still in progress. The combination of SHG imaging and the ST is able to objectively discriminate the changes suffered by the collagen arrangement after the CXL treatment, whose effectiveness depends on the initial organization of the collagen fibers within the corneal stroma.