Simultaneous microstructural and mechanical characterization of human corneas at increasing pressure

Measuring intraocular pressure with OCT: the first approach

The variability of corneal OCT speckle statistics is indirectly related to changes in corneal microstructure, which may be induced by intraocular pressure (IOP). A new approach is considered, which attempts to estimate IOP based on corneal speckle statistics in OCT images. An area (A) under trajectories of contrast ratio with respect to stromal depth was calculated. The proposed method was evaluated on OCT images from the ex-vivo study on porcine eyeballs and in-vivo study on human corneas. A statistically significant multivariate linear regression model was obtained from the ex-vivo study: IOP = 0.70 · A - 6.11, in which IOP was precisely controlled in the anterior chamber. The ex-vivo study showed good correlation between A and IOP (R = 0.628, at the least) whereas the in-vivo study showed poor correlation between A and clinical air-puff tonometry based estimates of IOP (R = 0.351, at the most), indicating substantial differences between the two studies. The results of the ex-vivo study show the potential for OCT speckle statistics to be utilized for measuring IOP using static corneal imaging that does not require corneal deformation. Nevertheless, further work is needed to validate this approach in living human corneas.

Unveiling the lamellar structure of the human cornea over its full thickness using polarization-resolved SHG microscopy

A key property of the human cornea is to maintain its curvature and consequently its refraction capability despite daily changes in intraocular pressure. This is closely related to the multiscale structure of the corneal stroma, which consists of 1-3 µm-thick stacked lamellae made of thin collagen fibrils. Nevertheless, the distribution, size, and orientation of these lamellae along the depth of the cornea are poorly characterized up to now. In this study, we use second harmonic generation (SHG) microscopy to visualize the collagen distribution over the full depth of 10 intact and unstained human corneas (500-600 µm thick). We take advantage of the small coherence length in epi-detection to axially resolve the lamellae while maintaining the corneal physiological curvature. Moreover, as raw epi-detected SHG images are spatially homogenous because of the sub-wavelength size of stromal collagen fibrils, we use a polarimetric approach to measure the collagen orientation in every voxel. After a careful validation of this approach, we show that the collagen lamellae (i) are mostly oriented along the inferior-superior axis in the anterior stroma and along the nasal-temporal axis in the posterior stroma, with a gradual shift in between and (ii) exhibit more disorder in the anterior stroma. These results represent the first quantitative characterization of the lamellar structure of the human cornea continuously along its entire thickness with micrometric resolution. It also shows the unique potential of P-SHG microscopy for imaging of collagen distribution in thick dense tissues.

Noninvasive quantitative assessment of collagen degradation in parchments by polarization-resolved SHG microscopy

Nondestructive and noninvasive investigation techniques are highly sought-after to establish the degradation state of historical parchments, which is up to now assessed by thermal techniques that are invasive and destructive. We show that advanced nonlinear optical (NLO) microscopy enables quantitative in situ mapping of parchment degradation at the micrometer scale. We introduce two parameters that are sensitive to different degradation stages: the ratio of two-photon excited fluorescence to second harmonic generation (SHG) signals probes severe degradation, while the anisotropy parameter extracted from polarization-resolved SHG measurements is sensitive to early degradation. This approach is first validated by comparing NLO quantitative parameters to thermal measurements on artificially altered contemporary parchments. We then analyze invaluable parchments from the Middle Ages and show that we can map their conservation state and assess the impact of a restoration process. NLO quantitative microscopy should therefore help to identify parchments most at risk and optimize restoration methods.

Quantitative structural imaging of keratoconic corneas using polarization-resolved SHG microscopy

The human cornea is mainly composed of collagen fibrils aligned together within stacked lamellae. This lamellar structure can be affected in pathologies such as keratoconus, which is characterized by progressive corneal thinning and local steepening. In this study, we use polarization-resolved second harmonic generation (P-SHG) microscopy to characterize 8 control and 6 keratoconic human corneas. Automated processing of P-SHG images of transverse sections provides the collagen orientation in every pixel with sub-micrometer resolution. Series of P-SHG images recorded in the most anterior region of the stroma evidence sutural lamellae inclined at 22° ± 5° to the corneal surface, but show no significant difference between control and keratoconic corneas. In contrast, series of P-SHG images acquired along the full thickness of the stroma show a loss of order in the lamellar structure of keratoconic corneas, in agreement with their defective mechanical properties. This structural difference is analyzed quantitatively by computing the entropy and the orientation index of the collagen orientation distribution and significant differences are obtained along the full thickness of the stroma. This study shows that P-SHG is an effective tool for automatic quantitative analysis of structural defects of human corneas and should be applied to other collagen-rich tissues.

Corneal Optical Coherence Tomography Speckle in Crosslinked and Untreated Rabbit Eyes in Response to Elevated Intraocular Pressure

Purpose

To ascertain the influence of intraocular pressure (IOP) on corneal optical coherence tomography (OCT) speckle in untreated and ultraviolet A–riboflavin induced corneal collagen crosslinked rabbit eyes.

Methods

Left corneas of eight rabbits were de-epithelialized and crosslinked by applying riboflavin and 30-minute ultraviolet A light exposure. After enucleation (6 months after treatment), each eyeball (treated and untreated) was mounted in a measurement setup, in which IOP was increased from 15 to 45 mm Hg in steps of 5 mm Hg. At each IOP value, single B-scans of the central cornea were acquired three times with the spectral-domain OCT Copernicus-HR. Then, three regions of interest, including the anterior, posterior, and entire corneal stroma, were automatically extracted. Five different probability distributions were used as a model for the corneal speckle and the one with the best goodness of fit was chosen for further analysis.

Results

The generalized gamma distribution achieved the best goodness of fit and its scale (a) and shape (v) parameters statistically significantly changed with increasing IOP in the three regions of analysis (two-way repeated measures analysis of variance, all P < 0.05). The statistically significant difference between untreated and crosslinked eyes was observed for the shape parameters of the posterior and entire corneal stroma.

Conclusions

Corneal OCT speckle is influenced by IOP and shows to be significantly different in untreated and crosslinked eyes. Corneal OCT speckle analysis has the potential to be indirectly used for assessing changes in corneal stroma in ex vivo and in vivo studies.

Translational Relevance

Investigation of corneal OCT speckle statistics can offer additional diagnostic biomarkers related to changes in the corneal stroma after ocular surgeries.

Comparative analysis of the corneal birefringence pattern in healthy children and adults

PURPOSE

To undertake a comparative analysis of the corneal shape, thickness and isochromatics in the eyes of children and adults in order to determine the extent of similarities and differences between the cohorts.

METHODS

The study involved 24 children (aged 8 years) and 37 young White adults (aged between 22-24 years) with no apparent or known health or ocular conditions. Measurements were made of corneal radius of curvature, both central (CCT) and paracentral (PCT) corneal thickness and intraocular pressure (IOP). Images of the isochromatics were captured using a slit lamp and a circular polarizer. The geometry of fringe I and II of the isochromatics was analysed.

RESULTS

Statistically significant differences were found between CCT and PCT in nasal and temporal regions for both the children and adult cohorts. The same trends were observed in the radii of the cornea. Statistically significant differences between side lengths and angles of isochromatic fringes were found. No differences in asymmetry of shape for fringe I between adults and children were detected; greater symmetry was seen in fringe II in children than for adults.

CONCLUSIONS

The asymmetry in corneal shape and curvature contributes to the shape of the isochromatic fringes. This is likely linked to the orientation and parameters of the collagen fibres and to the muscles' forces, and be relevant for surgical procedures such as corneal transplantation.

The effect of intraocular pressure elevation and related ocular biometry changes on corneal OCT speckle distribution in porcine eyes

The aim of this study was to evaluate the influence of increase in intraocular pressure (IOP) and cooccurring changes in ocular biometry parameters on the corneal optical coherence tomography (OCT) speckle distribution in ex-vivo experiments on porcine intact eyes. Twenty-three eyeballs were used in the inflation test where IOP in the anterior chamber was precisely set from 10 mmHg to 40 mmHg in steps of 5 mmHg and where eye biometry was utilized (IOL Master 700). To assess the influence of the duration of the experiment on the OCT speckle statistics, the second experiment was performed with 10 eyeballs at the constant IOP of 15 mmHg. Based on the OCT scans of central cornea (Copernicus REVO), spatial maps of the scale parameter (a) and the shape parameter (v) of the gamma distribution speckle model were estimated. The means of both parameters for each spatial map were computed within the 2 mm of the central stroma. Both distributional parameters statistically significantly varied with IOP and time (one way repeated measures ANOVA, all p-values < 0.001). The a parameter revealed a faster statistically significant increase in IOP up to 25 mmHg, regardless of time. Central corneal thickness (CCT), the anterior chamber depth, and the mean equivalent spherical power varied significantly with IOP, whereas CCT and axial length changed statistically significantly with time. Statistically significant correlation was found between CCT and the a parameter, after removing IOP as a confounding factor (r = −0.576, p < 0.001). The parameters of the gamma distribution can be used not only for identifying IOP induced changes in the optical scattering within the corneal stroma, but also in corneal geometry. The approach of corneal speckle analysis could be potentially utilized for an indirect and noninvasive assessment of some properties of corneal stroma.

Theoretical Analysis of Wave-Front Aberrations Induced from Conventional Laser Refractive Surgery in a Biomechanical Finite Element Model

Purpose

To examine the biomechanical effects-induced wave-front aberrations after conventional laser refractive surgery.

Methods

A finite element model of the human eye was established to simulate conventional laser refractive surgery with corrected refraction from –1 to –15 diopters (D). The deformation of the anterior and posterior corneal surfaces was obtained under the intraocular pressure (IOP). Then, the surface displacement was converted to wave-front aberrations.

Results

Following conventional refractive surgery, significant deformation of the anterior and posterior corneal surfaces occurred because of the corneal biomechanical effects, resulting in increased residual wave-front aberrations. Deformation of the anterior surface resulted in a hyperopic shift, which was significantly increased with the increasing refractive correction. The residual high-order aberrations consisted of spherical aberration, vertical coma, and y-trefoil. Spherical aberration was significantly positively correlated to enhanced refraction correction. The effect of posterior corneal surface on induced wave-front aberration was less than the anterior corneal surface. The IOP slightly affects the postoperative defocus, coma, and spherical aberration. When treatment decentration occurred during the procedure, the hyperopic shift decreased as the eccentricity increased. Treatment decentration had a significant impact on the spherical aberration and the coma. In addition, the ocular tissue elasticity played a key role in hyperopic shift, whereas it had little effect on the other aberrations.

Conclusions

Among the many factors that affect high-order aberrations after conventional laser refractive surgery, the alterations in corneal morphology caused by biomechanical effects must be considered, as they can lead to an increase in postoperative residual wave-front aberrations.

Second harmonic generation imaging reveals asymmetry in the rotational helicity of collagen lamellae in chicken corneas

High tensile strength and optical clarity are unique properties of the cornea. These features are dictated by the three-dimensional architecture of corneal lamellae. Therefore, understanding the microscopic details of the cornea's structural organization may contribute to the development of artificial cornea for the treatment of corneal diseases. In this study, the combination of forward second harmonic generation (SHG) microcopy and fast Fourier-transform based image analysis was used to characterize the depth-dependent superstructure of chicken corneal stroma. Our results show that from the surface, adjacent lamellae of anterior chicken cornea lamella rotate in a counterclockwise direction, and the same rotational helicity is observed in left and right corneas. Furthermore, the overall average rotational pitch of lamellae is 0.92 ± 0.11 degree/µm which persists for 176 ± 14 µm in the anterior stroma. As depth further increased, the rate of lamellar rotation decreases. Upon reaching posterior stroma, lamellar orientation remains constant. Throughout the stroma, collagen lamellae in chicken rotate a total of 169 ± 21 degrees. The lack of lamellar rotation in posterior stroma suggests that packing efficiency cannot be used to explain the helicity of depth-dependent rotation of anterior stroma. In addition, although the right cornea has a higher rotational pitch (0.95 ± 11 vs 0.90 ± 10 degrees/µm) and thinner anterior stroma (173 ± 13 vs 179 ± 14 µm) than the left cornea, the two effects cancel each other out and result in similar total angular rotation of anterior stroma (161 ± 23 and 165 degrees ± 21). Finally, our observation of a total angular rotation of 169 ± 21 degrees shows that within experimental error, chicken cornea lamellae rotate around 180 degrees or half of a complete turn. Additional studies are needed to arrive at an explanation of chicken superstructure in three dimensions.

Super-shear evanescent waves for non-contact elastography of soft tissues

We describe surface wave propagation in soft elastic media at speeds exceeding the bulk shear wave speed. By linking these waves to the elastodynamic Green's function, we derive a simple relationship to quantify the elasticity of a soft medium from the speed of this supershear evanescent wave (SEW). We experimentally probe SEW propagation in tissue-mimicking phantoms, human cornea ex vivo, and skin in vivo using a high-speed optical coherence elastography system. Measurements confirm the predicted relationship between SEW and bulk shear wave speeds, agreeing well with both theoretical and numerical models. These results suggest that SEW measurements may be a robust method to quantify elasticity in soft media, particularly in complex, bounded materials where dispersive Rayleigh-Lamb modes complicate measurements.

Implementation of artifact-free circular dichroism SHG imaging of collagen

Second harmonic generation (SHG) enables in situ imaging of fibrillar collagen architecture in connective tissues. Recently, Circular Dichroism SHG (CD-SHG) microscopy has been implemented to take advantage of collagen chirality to improve 3D visualization. It measures the normalized difference in the SHG signal obtained upon excitation by left versus right circular polarizations. However, CD-SHG signal is not well characterized yet, and quite different CD-SHG values are reported in the literature. Here, we identify two major artifacts that may occur in CD-SHG experiments and we demonstrate that thorough optimization and calibration of the experimental setup are required for CD-SHG imaging. Notably it requires a careful calibration of the incident circular polarizations and a perfect mechanical stabilization of the microscope stage. Finally, we successfully record CD-SHG images in human cornea sections and confirm that this technique efficiently reveals collagen fibrils oriented out of the focal plane.

Monitoring dynamic collagen reorganization during skin stretching with fast polarization-resolved second harmonic generation imaging

The mechanical properties of biological tissues are strongly correlated to the specific distribution of their collagen fibers. Monitoring the dynamic reorganization of the collagen network during mechanical stretching is however a technical challenge, because it requires mapping orientation of collagen fibers in a thick and deforming sample. In this work, a fast polarization-resolved second harmonic generation microscope is implemented to map collagen orientation during mechanical assays. This system is based on line-to-line switching of polarization using an electro-optical modulator and works in epi-detection geometry. After proper calibration, it successfully highlights the collagen dynamic alignment along the traction direction in ex vivo murine skin dermis. This microstructure reorganization is quantified by the entropy of the collagen orientation distribution as a function of the stretch ratio. It exhibits a linear behavior, whose slope is measured with a good accuracy. This approach can be generalized to probe a variety of dynamic processes in thick tissues.

Analysis of X-ray scattering microstructure data for implementation in numerical simulations of ocular biomechanical behaviour

This study aimed to analyse microstructure data on the density and orientation of collagen fibrils in whole eye globes and to propose an effective method for the preparation of data for use in numerical simulations of the eye’s biomechanical performance. Wide-angle X-ray scattering was applied to seven healthy ex-vivo human eyes. Each eye was dissected into an anterior and a posterior cup, and radial incisions were used to flatten the tissue before microstructure characterisation. A method was developed to use the microstructure data obtained for the dissected tissue to build realistic 3D maps of fibril density and orientation covering the whole eye globe. At the central cornea, 61.5±2.3% of fibrils were aligned within 45° sectors surrounding the two orthogonal directions. In contrast, more than one-third of the total fibril content was concentrated along the circumferential direction at the limbus (37.0±2.4%) and around the optic nerve head (34.8±2.1%). The insertion locations of the four recti muscles exhibited a preference in the meridional direction near the equator (38.6±3.9%). There was also a significant difference in fibril density between the limbus and other regions (ratio = 1.91±0.45, p <0.01 at the central cornea and ratio = 0.80±0.21, p <0.01 at the posterior pole). Characterisation of collagen fibril density and orientation across the whole ocular surface has been possible but required the use of a technique that involved tissue dissection and hence caused tissue damage. The method presented in this paper aimed to minimise the effect of dissection on the quality of obtained data and was successful in identifying fibril distribution trends that were compatible with earlier studies, which concentrated on localised areas of the ocular globe.

Combination of Traction Assays and Multiphoton Imaging to Quantify Skin Biomechanics

An important issue in tissue biomechanics is to decipher the relationship between the mechanical behavior at macroscopic scale and the organization of the collagen fiber network at microscopic scale. Here, we present a protocol to combine traction assays with multiphoton microscopy in ex vivo murine skin. This multiscale approach provides simultaneously the stress/stretch response of a skin biopsy and the collagen reorganization in the dermis by use of second harmonic generation (SHG) signals and appropriate image processing.

Corneal biomechanics after laser refractive surgery: Unmasking differences between techniques

The hypothesis that small-incision lenticule extraction provides better preservation of corneal biomechanics than previous laser refractive techniques has led to a growth in the interest in clinical and experimental research in this field. This hypothesis is based on the fact that corneal layers with greater stiffness are preserved with this new technique. However, this hypothesis is controversial because clinical research has shown a great disparity in the outcomes. In this review, we performed an in-depth analysis of the factors that might affect corneal biomechanics in laser refractive surgery procedures from a macrostructural to a microstructural viewpoint. New advances in algorithms with current devices or the introduction of new devices might help unmask the possible advantages of small-incision lenticule extraction in corneal biomechanics.

How aging impacts skin biomechanics: a multiscale study in mice

Skin aging is a complex process that strongly affects the mechanical behavior of skin. This study aims at deciphering the relationship between age-related changes in dermis mechanical behavior and the underlying changes in dermis microstructure. To that end, we use multiphoton microscopy to monitor the reorganization of dermal collagen during mechanical traction assays in ex vivo skin from young and old mice. The simultaneous variations of a full set of mechanical and microstructural parameters are analyzed in the framework of a multiscale mechanical interpretation. They show consistent results for wild-type mice as well as for genetically-modified mice with modified collagen V synthesis. We mainly observe an increase of the tangent modulus and a lengthening of the heel region in old murine skin from all strains, which is attributed to two different origins that may act together: (i) increased cross-linking of collagen fibers and (ii) loss of water due to proteoglycans deterioration, which impedes inner sliding within these fibers. In contrast, the microstructure reorganization upon stretching shows no age-related difference, which can be attributed to opposite effects of the decrease of collagen content and of the increase of collagen cross-linking in old mice.

Stromal striae: a new insight into corneal physiology and mechanics

We uncover the significance of a previously unappreciated structural feature in corneal stroma, important to its biomechanics. Vogt striae are a known clinical indicator of keratoconus, and consist of dark, vertical lines crossing the corneal depth. However we detected stromal striae in most corneas, not only keratoconus. We observed striae with multiple imaging modalities in 82% of 118 human corneas, with pathology-specific differences. Striae generally depart from anchor points at Descemet’s membrane in the posterior stroma obliquely in a V-shape, whereas in keratoconus, striae depart vertically from posterior toward anterior stroma. Optical coherence tomography shear wave elastography showed discontinuity of rigidity, and second harmonic generation and scanning electron microscopies showed undulation of lamellae at striae locations. Striae visibility decreased beyond physiological pressure and increased beyond physiological hydration. Immunohistology revealed striae to predominantly contain collagen VI, lumican and keratocan. The role of these regions of collagen VI linking sets of lamellae may be to absorb increases in intraocular pressure and external shocks.

Recent advances in studying single bacteria and biofilm mechanics

Bacterial biofilms correspond to surface-associated bacterial communities embedded in hydrogel-like matrix, in which high cell density, reduced diffusion and physico-chemical heterogeneity play a protective role and induce novel behaviors. In this review, we present recent advances on the understanding of how bacterial mechanical properties, from single cell to high-cell density community, determine biofilm tri-dimensional growth and eventual dispersion and we attempt to draw a parallel between these properties and the mechanical properties of other well-studied hydrogels and living systems.

Third-harmonic generation Stokes-Mueller polarimetric microscopy

An experimental implementation of the nonlinear Stokes-Mueller polarimetric (NSMP) microscopy in third-harmonic generation modality is presented. The technique is able to extract all eight 2D-accessible χ⁽³⁾ components for any sample from 64 polarization measurements, and can be applied to noninvasive ultrastructural characterization. The polarization signature of an isotropic glass coverslip is presented, and carotenoid crystallites in the root of orange carrot (Daucus carota) are investigated, showing complex χ⁽³⁾ components with a significant chiral contribution.