Development of a novel instrument to measure the pulsatile movement of ocular tissues

Synchronous functional magnetic resonance eye imaging, video ophthalmoscopy, and eye surface imaging reveal the human brain and eye pulsation mechanisms

The eye possesses a paravascular solute transport pathway that is driven by physiological pulsations, resembling the brain glymphatic pathway. We developed synchronous multimodal imaging tools aimed at measuring the driving pulsations of the human eye, using an eye-tracking functional eye camera (FEC) compatible with magnetic resonance imaging (MRI) for measuring eye surface pulsations. Special optics enabled integration of the FEC with MRI-compatible video ophthalmoscopy (MRcVO) for simultaneous retinal imaging along with functional eye MRI imaging (fMREye) of the BOLD (blood oxygen level dependent) contrast. Upon optimizing the fMREye parameters, we measured the power of the physiological (vasomotor, respiratory, and cardiac) eye and brain pulsations by fast Fourier transform (FFT) power analysis. The human eye pulsated in all three physiological pulse bands, most prominently in the respiratory band. The FFT power means of physiological pulsation for two adjacent slices was significantly higher than in one-slice scans (RESP1 vs. RESP2; df = 5, p = 0.045). FEC and MRcVO confirmed the respiratory pulsations at the eye surface and retina. We conclude that in addition to the known cardiovascular pulsation, the human eye also has respiratory and vasomotor pulsation mechanisms, which are now amenable to study using non-invasive multimodal imaging of eye fluidics.

Noncontact Acoustic Micro-Tapping Optical Coherence Elastography for Quantification of Corneal Anisotropic Elasticity: In Vivo Rabbit Study

Purpose

The purpose of this study was to demonstrate accurate measurement of corneal elastic moduli in vivo with noncontact and noninvasive optical coherence elastography.

Methods

Elastic properties (in-plane Young's modulus, E, and both in-plane, μ, and out-of-plane, G, shear moduli) of rabbit cornea were quantified in vivo using noncontact dynamic acoustic micro-tapping optical coherence elastography (AµT-OCE). The intraocular pressure (IOP)-dependence of measured mechanical properties was explored in extracted whole globes following in vivo measurement. A nearly incompressible transverse isotropic (NITI) model was used to reconstruct moduli from AµT-OCE data. Independently, cornea elastic moduli were also measured ex vivo with traditional, destructive mechanical tests (tensile extensometry and shear rheometry).

Results

Our study demonstrates strong anisotropy of corneal elasticity in rabbits. The in-plane Young's modulus, computed as E = 3μ, was in the range of 20 MPa to 44 MPa, whereas the out-of-plane shear modulus was in the range of 34 kPa to 261 kPa. Both pressure-dependent ex vivo OCE and destructive mechanical tests performed on the same samples within an hour of euthanasia strongly support the results of AµT-OCE measurements.

Conclusions

Noncontact AµT-OCE can noninvasively quantify cornea anisotropic elastic properties in vivo.

Translational Relevance

As optical coherence tomography (OCT) is broadly accepted in ophthalmology, these results suggest the potential for rapid translation of AµT-OCE into clinical practice. In addition, AµT-OCE can likely improve diagnostic criteria of ectatic corneal diseases, leading to early diagnosis, reduced complications, customized surgical treatment, and personalized biomechanical models of the eye.

Retinal blood flow reversal quantitatively monitored in out-of-plane vessels with laser Doppler holography

Laser Doppler holography is a planar blood flow imaging technique recently introduced in ophthalmology to image human retinal and choroidal blood flow non-invasively. Here we present a digital method based on the Doppler spectrum asymmetry that reveals the local direction of blood flow with respect to the optical axis in out-of-plane vessels. This directional information is overlaid on standard grayscale blood flow images to depict flow moving towards the camera in red and flow moving away from the camera in blue, as in ultrasound color Doppler imaging. We show that thanks to the strong contribution of backscattering to the Doppler spectrum in out-of-plane vessels, the local axial direction of blood flow can be revealed with a high temporal resolution, which enables us to evidence pathological blood flow reversals. We also demonstrate the use of optical Doppler spectrograms to quantitatively monitor retinal blood flow reversals.

Laser Doppler holography of the anterior segment for blood flow imaging, eye tracking, and transparency assessment

Laser Doppler holography (LDH) is a full-field blood flow imaging technique able to reveal human retinal and choroidal blood flow with high temporal resolution. We here report on using LDH in the anterior segment of the eye without making changes to the instrument. Blood flow in the bulbar conjunctiva and episclera as well as in corneal neovascularization can be effectively imaged. We additionally demonstrate simultaneous holographic imaging of the anterior and posterior segments by simply adapting the numerical propagation distance to the plane of interest. We used this feature to track the movements of the retina and pupil with high temporal resolution. Finally, we show that the light backscattered by the retina can be used for retro-illumination of the anterior segment. Hence digital holography can reveal opacities caused by absorption or diffusion in the cornea and eye lens.

Blind Source Separation of Retinal Pulsatile Patterns in Optic Nerve Head Video-Recordings

Dynamic optical imaging of retinal hemodynamics is a rapidly evolving technique in vision and eye-disease research. Video-recording, which may be readily accessible and affordable, captures several distinct functional phenomena such as the spontaneous venous pulsations (SVP) of central vein or local arterial blood supply etc. These phenomena display specific dynamic patterns that have been detected using manual or semi-automated methods. We propose a pioneering concept in retina video-imaging using blind source separation (BSS) serving as an automated localizer of distinct areas with temporally synchronized hemodynamics. The feasibility of BSS techniques (such as spatial principal component analysis and spatial independent component analysis) and K-means based post-processing method were successfully tested on the monocular and binocular video-ophthalmoscopic (VO) recordings of optic nerve head (ONH) in healthy subjects. BSSs automatically detected three spatially distinct reproducible areas, i.e. SVP, optic cup pulsations (OCP) that included areas of larger vessels in the nasal part of ONH, and "other" pulsations (OP). The K-means post-processing reduced a spike noise from the patterns' dynamics while high linear dependence between the non-filtered and post-processed signals was preserved. Although the dynamics of all patterns were heart rate related, the morphology analysis demonstrated significant phase shifts between SVP and OCP, and between SVP and OP. In addition, we detected low frequency oscillations that may represent respiratory-induced effects in time-courses of the VO recordings.

Spatio-temporal filtering in laser Doppler holography for retinal blood flow imaging

Laser Doppler holography (LDH) is a full-field interferometric imaging technique recently applied in ophthalmology to measure blood flow, a parameter of high clinical interest. From the temporal fluctuations of digital holograms acquired at ultrafast frame rates, LDH reveals retinal and choroidal blood flow with a few milliseconds of temporal resolution. However, LDH experiences difficulties to detect slower blood flow as it requires to work with low Doppler frequency shifts which are corrupted by eye motion. We here demonstrate the use of a spatio-temporal decomposition adapted from Doppler ultrasound that provides a basis appropriate to the discrimination of blood flow from eye motion. A singular value decomposition (SVD) can be used as a simple, robust, and efficient way to separate the Doppler fluctuations of blood flow from those of strong spatial coherence such as eye motion. We show that the SVD outperforms the conventional Fourier based filter to reveal slower blood flow, and dramatically improves the ability of LDH to reveal vessels of smaller size or with a pathologically reduced blood flow.

Full-field displacement measurement of corneoscleral shells by combining multi-camera speckle interferometry with 3D shape reconstruction

Changes in the biomechanical properties of the connective tissue of the eye occur with age and underlie the development of several ocular diseases, such as glaucoma, myopia, and keratoconus. The biomechanical dynamics of ocular connective tissue are measured by ex vivo inflation testing, in which intraocular pressure (IOP) is varied and optical methods are used to produce maps of corneal and scleral displacement. Current optical methods are limited by acquisition rate, occlusions, poor spatial resolution, and insufficient 3D mapping. We developed an interferometric optical method integrates four-camera electronic speckle pattern interferometry (ESPI) and a novel three-dimensional (3D) shape reconstruction process to measure shape and full-field mechanical deformations of corneal and scleral shells during ex vivo inflation testing. Each camera provides accurate measurements of the laser beam phase related to deformations of the specimen surface; a multi-view stereovision method generates the shape of the specimen and a functional form that links every pixel of a given camera to 3D points on the specimen's visible surface. In this way, dynamic deformations of the specimen are localized, with quantification of the time-dependent 3D displacements of the specimen at nanometric accuracy. The ESPI-3D system is suitable for analyzing scleral deformation and morphological changes caused by time-varying IOP.

Relationship Between the Parameters of Corneal and Fundus Pulse Signals Acquired With a Combined Ultrasound and Laser Interferometry Technique

Purpose

To estimate the relationship between the characteristics of the corneal pulse (CP) signal and those of the fundus pulse (FP) signal measured with a combined noncontact ultrasonic and laser interferometry technique in healthy subjects.

Methods

Twenty-two healthy subjects participated in experiments that included measurements of intraocular pressure, ocular pulse amplitude, ocular biometry, blood pressure, and heart rate. Additionally, simultaneous recordings of CP and FP signals were acquired with a noncontact ultrasonic device combined with laser interferometry. Subsequently, ocular perfusion pressure (OPP) and the time and spectral parameters of CP and FP signals were computed. A system model was proposed to relate the FP signal to the CP signal.

Results

The system model revealed that the eye globe transfers information between signals of the posterior and anterior eye, relatively amplifying higher spectral harmonics. The amplitude of the second CP harmonic is predicted by FP_RMS and OPP (R² = 0.468, P = 0.002). Partial correlation analysis showed that the CP signal parameters are statistically significantly correlated with those of the FP signal and OPP, after correcting for age and sex.

Conclusions

The eye globe can be viewed as a high pass filter, in which the CP characteristic changes in relation to the fundus pulsation. The FP signal and OPP have an impact on the variations of the CP signal morphology.

Translational Relevance

Investigation of differences between the characteristics of the anterior and posterior tissue movements is a promising method for evaluating the role of circulatory and biomechanical components in the pathophysiology of ocular diseases.

Ocular fundus pulsations within the posterior rat eye: Chorioscleral motion and response to elevated intraocular pressure

A multi-functional optical coherence tomography (OCT) approach is presented to determine ocular fundus pulsations as an axial displacement between the retina and the chorioscleral complex in the albino rat eye. By combining optical coherence elastography and OCT angiography (OCTA), we measure subtle deformations in the nanometer range within the eye and simultaneously map retinal and choroidal perfusion. The conventional OCT reflectivity contrast serves as a backbone to segment the retina and to define several slabs which are subsequently used for quantitative ocular pulsation measurements as well as for a qualitative exploration of the multi-functional OCT image data. The proposed concept is applied in healthy albino rats as well as in rats under acute elevation of the intraocular pressure (IOP). The evaluation of this experiment revealed an increased pulsatility and deformation between the retinal and chorioscleral complex while increasing the IOP level from 15 mmHg to 65 mmHg. At IOP levels exceeding 65 mmHg, the pulsatility decreased significantly and retinal as well as choroidal perfusion vanished in OCTA. Furthermore, the evaluation of the multi-parametric experiment revealed a spatial correlation between fundus pulsatility and choroidal blood flow. This indicates that the assessed pulsatility may be a valuable parameter describing the choroidal perfusion.

In-vivo corneal pulsation in relation to in-vivo intraocular pressure and corneal biomechanics assessed in-vitro. An animal pilot study

The aim was to ascertain whether the characteristics of the corneal pulse (CP) measured in-vivo in a rabbit eye change after short-term artificial increase of intraocular pressure (IOP) and whether they correlate with corneal biomechanics assessed in-vitro. Eight New Zealand white rabbits were included in this study and were anesthetized. In-vivo experiments included simultaneous measurements of the CP signal, registered with a non-contact method, IOP, intra-arterial blood pressure, and blood pulse (BPL), at the baseline and short-term elevated IOP. Afterwards, thickness of post-mortem corneas was determined and then uniaxial tensile tests were conducted leading to estimates of their Young's modulus (E). At the baseline IOP, backward stepwise regression analyses were performed in which successively the ocular biomechanical, biometric and cardiovascular predictors were separately taken into account. Results of the analysis revealed that the 3rd CP harmonic can be statistically significantly predicted by E and central corneal thickness (Models: R² = 0.662, p < 0.005 and R² = 0.832, p < 0.001 for the signal amplitude and power, respectively). The 1st CP harmonic can be statistically significantly predicted by the amplitude and power of the 1st BPL harmonic (Models: R² = 0.534, p = 0.015 and R² = 0.509, p < 0.018, respectively). For elevated IOP, non-parametric analysis indicated significant differences for the power of the 1st CP harmonic (Kruskal-Wallis test; p = 0.031) and for the mean, systolic and diastolic blood pressures (p = 0.025, p = 0.019, p = 0.033, respectively). In conclusion, for the first time, the association between parameters of the CP signal in-vivo and corneal biomechanics in-vitro was confirmed. In particular, spectral analysis revealed that higher amplitude and power of the 3rd CP harmonic indicates higher corneal stiffness, while the 1st CP harmonic correlates positively with the corresponding harmonic of the BPL signal.

Optical Coherence Tomography as a Tool for Ocular Dynamics Estimation

Purpose. The aim of the study is to demonstrate that the ocular dynamics of the anterior chamber of the eye can be estimated quantitatively by means of optical coherence tomography (OCT). Methods. A commercial high speed, high resolution optical coherence tomographer was used. The sequences of tomographic images of the iridocorneal angle of three subjects were captured and each image from the sequence was processed in MATLAB environment in order to detect and identify the contours of the cornea and iris. The data on pulsatile displacements of the cornea and iris and the changes of the depth of the gap between them were retrieved from the sequences. Finally, the spectral analysis of the changes of these parameters was performed. Results. The results of the temporal and spectral analysis manifest the ocular microfluctuation that might be associated with breathing (manifested by 0.25 Hz peak in the power spectra), heart rate (1–1.5 Hz peak), and ocular hemodynamics (3.75–4.5 Hz peak). Conclusions. This paper shows that the optical coherence tomography can be used as a tool for noninvasive estimation of the ocular dynamics of the anterior segment of the eye, but its usability in diagnostics of the ocular hemodynamics needs further investigations.

Visualization of the ocular pulse in the anterior chamber of the mouse eye in vivo using phase-sensitive optical coherence tomography

We report on a phase-based method for accurately measuring the ocular pulse in the anterior chamber in vivo. Using phase-sensitive optical coherence tomography with optimized scanning protocols and equations for compensating bulk motion and environmental vibrations, a high sensitivity of 0.9 μm/s minimal velocity is demonstrated at a wide detection band of 0 to 380 Hz. The pulsatile relative motion between cornea and crystalline lens in rodents is visualized and quantified. The relative motion is most likely caused by respiration (1.6 Hz) and heartbeat (6.6 Hz). The velocity amplitude of the relative motion is 10.3 ± 2.4 μm/s. The displacement amplitudes at the respiratory and cardiac frequencies are 202.5 ± 64.9 and 179.9 ± 49.4 nm, respectively. The potential applications the measurement technique can be found in the evaluation of intraocular pressure and the measurement of biomechanical properties of the ocular tissue, which are important in several ocular diseases.

The collagen matrix of the human trabecular meshwork is an extension of the novel pre-Descemet's layer (Dua's layer)

BACKGROUND

The trabecular meshwork (TM) located at the angle of the anterior chamber of the eye contributes to aqueous drainage. A novel layer in the posterior part of the human cornea has recently been reported (the pre-Descemet's layer (Dua's layer (PDL)). We examined the peripheral part of this layer in relation to the origin of the TM.

METHODS

The PDL and TM of 19 human donor eyes and one exenterated sample were studied. Samples were examined by light and electron microscopy (EM) for tissue architecture and by immunohistology for four matricellular proteins, five collagen types and CD34.

RESULTS

EM revealed that beams of collagen emerged from the periphery of PDL on the anterior surface of the Descemet's membrane and divided and subdivided to continue as the beams of the TM. Long-spacing collagen was seen in the PDL and TM. Trabecular cells (CD34-ve) associated with basement membrane were seen in the peripheral part of the PDL and corresponded to the start of the separation of the collagen lamellae of PDL. Collagen VI was present continuously in PDL and extended into the TM. Matricellular proteins were seen predominantly in the TM with only laminin extending into the periphery of PDL.

CONCLUSIONS

This study provides an insight into the origins of the collagen core of the TM as an extension of the PDL of the cornea. This finding adds to the knowledge base of the TM and cornea and has the potential to impact future research into the TM and glaucoma.

Variability of Heidelberg Retina Tomograph parameters during exercise

PURPOSE

The purpose of this study was to determine whether exercise affects the stereometric parameters representing optic nerve head (ONH) topography.

METHODS

ONH topography, intraocular pressure (IOP) and blood pressure of 30 healthy volunteers were monitored before, during and after dynamic exercise raising systolic blood pressure by a minimum of 30 mmHg. Change in the stereometric ONH parameters was calculated.

RESULTS

IOP decreased and blood pressure increased during exercise, resulting in an increase in mean ocular perfusion pressure. Exercise was associated with an increase in variance in 17 of the 18 stereometric ONH parameters. The increase in variance was statistically significant in eight parameters, including rim area, cup/disc area ratio and cup shape measure. There was no statistically significant change in image quality. The absolute change from baseline in rim area, cup area, cup/disc area ratio, rim/disc area ratio and rim volume showed a statistically significant (p < 0.05) correlation with change in mean ocular perfusion pressure.

CONCLUSIONS

Exercise increases variability in stereometric ONH parameters. To avoid increased variance in the stereometric parameters, ONH imaging should be performed after allowing sufficient time to rest.

Analysis of pulsatile retinal movements by spectral-domain low-coherence interferometry: influence of age and glaucoma on the pulse wave

Recent studies have shown that ocular hemodynamics and eye tissue biomechanical properties play an important role in the pathophysiology of glaucoma. Nevertheless, better, non-invasive methods to assess these characteristics in vivo are essential for a thorough understanding of degenerative mechanisms. Here, we propose to measure ocular tissue movements induced by cardiac pulsations and study the ocular pulse waveform as an indicator of tissue compliance. Using a novel, low-cost and non-invasive device based on spectral-domain low coherence interferometry (SD-LCI), we demonstrate the potential of this technique to differentiate ocular hemodynamic and biomechanical properties. We measured the axial movement of the retina driven by the pulsatile ocular blood flow in 11 young healthy individuals, 12 older healthy individuals and 15 older treated glaucoma patients using our custom-made SD-OCT apparatus. The cardiac pulse was simultaneously measured through the use of an oximeter to allow comparison. Spectral components up to the second harmonic were obtained and analyzed. For the different cohorts, we computed a few parameters that characterize the three groups of individuals by analyzing the movement of the retinal tissue at two locations, using this simple, low-cost interferometric device. Our pilot study indicates that spectral analysis of the fundus pulsation has potential for the study of ocular biomechanical and vascular properties, as well as for the study of ocular disease.

Smoothness assessment of corneal stromal surfaces

PURPOSE

To assess the accuracy of the scanning electron microscopy (SEM) and present alternative approaches to quantify surface roughness based on numerical analysis.

SETTING

Department of Ophthalmology, Maisonneuve-Rosemont Hospital, University of Montreal, Montreal, Quebec, Canada.

DESIGN

Experimental study.

METHODS

Lamellar stromal cuts were performed on human corneas using a femtosecond laser or a microkeratome. The photodisrupted stromal surfaces were processed for SEM, and images were acquired at ×1000 magnification. First, images were evaluated by independent observers. Second, images were analyzed based on first-order and second-order statistics of gray-level intensities. Third, 3-dimensional (3-D) surface reconstructions were generated from pairs of SEM images acquired at 2 angles.

RESULTS

Results show that traditional assessment of roughness based on evaluating SEM images by independent observers can be replaced by computer-image texture analysis; an algorithm was developed to avoid subjective and time-consuming observations. The 3-D reconstructions allowed additional characterization of surface properties that was not possible with SEM images alone. Significant fluctuations in surface height were lost, although they could be retrieved using 3-D reconstructions.

CONCLUSIONS

Image texture analysis allowed objective and repeatable assessment of stromal surface roughness; however, full assessments of surface-height fluctuations required 3-D reconstruction. These complementary methodologies offer a more comprehensive assessment of corneal surface roughness in clinical applications.

Terahertz detection using spectral domain interferometry

In this work, we present a novel method based on spectral domain interferometry for the electro-optic (EO) sampling of terahertz (THz) electric fields. This technique allows the use of thick crystals without the drawback of the over-rotation that may occur with intense THz sources, allowing longer temporal scans and thus, better spectral resolution. Using this technique, a phase difference of approximately 8898π can be measured, which is 18,000 times larger than the phase difference that could be measured using EO sampling.

Phase-sensitive optical coherence tomography characterization of pulse-induced trabecular meshwork displacement in ex vivo nonhuman primate eyes

Glaucoma is a blinding disease for which intraocular pressure (IOP) is the only treatable risk factor. The mean IOP is regulated through the aqueous outflow system, which contains the trabecular meshwork (TM). Considerable evidence indicates that trabecular tissue movement regulates the aqueous outflow and becomes abnormal during glaucoma; however, such motion has thus far escaped detection. The purpose of this study is to describe anovel use of a phase-sensitive optical coherence tomography (PhS-OCT) method to assess pulse-dependent TM movement. For this study, we used enucleated monkey eyes, each mounted in an anterior segment holder. A perfusion system was used to control the mean IOP as well as to provide IOP sinusoidal transients (amplitude 3 mmHg, frequency 1 pulse/second) in all experiments. Measurements were carried out at seven graded mean IOPs (5, 8, 10, 20, 30, 40, and 50 mm Hg). We demonstrate that PhS-OCT is sensitive enough to image/visualize TM movement synchronous with the pulse-induced IOP transients, providing quantitative measurements of dynamic parameters such as velocity, displacement, and strain rate that are important for assessing the biomechanical compliance of the TM. We find that the largest TM displacement is in the area closest to Schlemm's canal (SC) endothelium. While maintaining constant ocular pulse amplitude, an increase of mean IOP results in a decrease of TM displacement and mean size of the SC. These results demonstrate that the PhS-OCT is a useful imaging technique capable of assessing functional properties necessary to maintain IOP in a healthy range, offering a new diagnostic alternative for glaucoma.

Analyzing speckle contrast for HiLo microscopy optimization

HiLo microscopy is a recently developed technique that provides both optical sectioning and fast imaging with a simple implementation and at a very low cost. The methodology combines widefield and speckled illumination images to obtain one optically sectioned image. Hence, the characteristics of such speckle illumination ultimately determine the quality of HiLo images and the overall performance of the method. In this work, we study how speckle contrast influence local variations of fluorescence intensity and brightness profiles of thick samples. We present this article as a guide to adjust the parameters of the system for optimizing the capabilities of this novel technology.