Spectral characteristics of longitudinal corneal apex velocities and their relation to the cardiopulmonary system

The effect of changes in cardiovascular activity on corneal biomechanics and pulsation in rabbits

The aim was to assess the relationships between cardiovascular activity, corneal pulse characteristics, and corneal biomechanics in rabbits. Seventeen rabbits were randomly assigned to one of two anesthetic regimens to induce differences in arterial blood pressure and heart rate. Experimental protocol included measuring blood flow parameters in the ophthalmic artery by color Doppler imaging, corneal biomechanical parameters using a non-contact tonometer Corvis ST, and the corneal pulse (CP) signal using a non-contact ultrasonic technique. Statistically significantly lower mean values of normalized amplitudes of higher CP harmonics and changes in eight of the twelve corneal biomechanical parameters were observed in the rabbit group with lower arterial blood pressure and higher heart rate, intraocular pressure, and resistive index. The results of partial correlations showed that the CP signal energy and amplitude of its first harmonic correlate with the resistive index, diastolic and mean arterial pressures, whereas no statistically significant correlation was found between any of the CP parameters and intraocular pressure. Our pilot study indicates, for the first time, that non-contact and continuous measuring of corneal pulse allows indirectly assessing changes in cardiovascular activity when the confounding effect of intraocular pressure is eliminated.

Using the entropy of the corneal pulse signal to distinguish healthy eyes from eyes affected by primary open-angle glaucoma

OBJECTIVE

The purpose of this study was to evaluate whether the complexity of the corneal pulse (CP) signal can be used to differentiate patients with primary open-angle glaucoma (POAG) from healthy subjects.

APPROACH

The study sample consisted of 28 patients with POAG and a control, age-matched group of 30 subjects. After standard ophthalmic examination, the CP signal from a randomly selected eye of each participant was measured using non-contact ultrasonic micro-displacement measurement technology. After pre-processing, the complexity of the CP signal was estimated using refined composite multiscale fuzzy entropy (RCMFE) up to scale factor 50. The average RCMFE values were computed from three repeated measurements of the CP signals for each participant and each scale factor.

MAIN RESULTS

The complexity of the CP signal in glaucomatous eyes was higher than that observed in healthy ones. Also, RCMFE of the CP signal was found to differentiate (statistically significantly) between the two groups for scales in the range from 26 to 43. For these scales, the one for which the lowest p-value (t-test, p = 0.017) was obtained when comparing RCMFE between the two groups was selected as the optimal scale. Next, a receiver operating characteristic analysis for the optimal scale showed that the proposed approach of calculating the multiscale entropy of the CP signal has some potential to discriminate between patients with POAG and healthy controls (sensitivity, specificity and accuracy of 0.643, 0.700 and 0.672, respectively).

SIGNIFICANCE

In conclusion, RCMFE, as a complexity measure, may be considered an auxiliary indicator to support glaucoma diagnostics.

Corneal pulsation and biomechanics during induced ocular pulse. An ex-vivo pilot study

The purpose of this study was to ascertain the relationships between the amplitude of the corneal pulse (CP) signal and the parameters of corneal biomechanics during ex-vivo intraocular pressure (IOP) elevation experiments on porcine eyes with artificially induced ocular pulse cycles. Two experiments were carried out using porcine eyes. In the first one, a selected eye globe was subjected to three IOP levels (15, 30 and 45 mmHg), where changes in physical ocular pulse amplitude were controlled by infusion/withdrawal volumes (ΔV). In the second experiment, six eyes were subjected to IOP from 15 mmHg to 45 mmHg in steps of 5 mmHg with a constant ΔV, where corneal deformation parameters were measured using Corvis ST. In both experiments, at each IOP, the CP and IOP signals were acquired synchronically using a non-contact ultrasonic distance sensor and a pressure transmitter, respectively. Based on the amplitudes of the CP and IOP signals ocular pulse based corneal rigidity index (OPCRI) was calculated. Results indicate positive correlations between ΔV and the physical ocular pulse amplitude, and between ΔV and the corneal pulse amplitude (both p < 0.001). OPCRI was found to increase with elevated IOP. Furthermore, IOP statistically significantly differentiated changes in OPCRI, the amplitudes of CP and IOP signals and in most of the corneal deformation parameters (p < 0.05). The partial correlation analysis, with IOP as a control variable, revealed a significant correlation between the length of the flattened cornea during the first applanation (A1L) and the corneal pulse amplitude (p = 0.002), and between A1L and OPCRI (p = 0.003). In conclusion, this study proved that natural corneal pulsations, detected with a non-contact ultrasonic technique, reflect pressure-volume dynamics and can potentially be utilized to assess stiffness of the cornea. The proposed new rigidity index could be a simple approach to estimating corneal rigidity.

Changes in spectral parameters of corneal pulse following canaloplasty

PURPOSE

To ascertain whether changes in the spectral content of the corneal pulse (CP) signal, measured in vivo in primary open-angle glaucoma (POAG) patients, indirectly reflect changes in corneal biomechanics after canaloplasty.

METHODS

Fifteen eyes of 15 POAG patients who underwent canaloplasty combined with phacoemulsification were enrolled. Standard ophthalmic examinations were conducted before washout, pre-operatively, at days 1, 7, and 1, 3, 6, and 12 months after surgery. Non-contact measurements of the CP signal were performed at pre-washout, pre-operatively, and at 3, 6, and 12 months post-operatively. Then, amplitudes of the CP first five harmonics associated with the heart rate were estimated. Temporal changes of all considered parameters were tested at a Bonferroni-adjusted significance level set to 0.005.

RESULTS

A decrease in the amplitude of the first harmonic and an increase in the normalized amplitude of the third harmonic (A_CP3n) of the CP signal were noticed between the pre-washout and the pre-operative stages (p = 0.003 and p = 0.004, respectively). This corresponds to an increase in median intraocular pressure (IOP) values by 6.0 mmHg (p = 0.0045). After surgery, A_CP3n reached the highest value at 3 months post-operatively, compared with pre-washout level (p = 0.0045).

CONCLUSIONS

Alterations in corneoscleral stiffness caused by surgery are reflected in changes in the A_CP3n value. Hence, post-operative corneal biomechanics could be monitored indirectly by this supporting indicator that can be used to estimate the time at which measures of IOP are no longer biased by the changed cornea boundary conditions caused by canaloplasty.

CLINICAL TRIALS REGISTRATION

NCT02908633.

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.

The Bull's Eye Pattern of the Tear Film in Humans during Visual Fixation on En-Face Optical Coherence Tomography

The aim of the study was to define and characterize the optical behavior of the tear film during visual fixation in humans on en-face optical coherence tomography (OCT). We included 20 healthy participants, 60% female, aged from 25 to 42 years (33.05 ± 4.97 [mean ± SD]) and ten patients with severe dry eye, 50% female, aged from 26 to 42 years (33.7 ± 5.31). To perform high-resolution tear film imaging, participants were asked to gaze at the internal fixation point in the spectral-domain anterior segment OCT device, and meanwhile scanning session was executed at the following time-points after blinking: at the 2nd, 3rd, 4th, 5th, and 6th second. After one hour, OCT imaging was repeated (second session) by a different operator masked to the study to verify the reliability of results. During each measuring session, a pulse oximetry was used for continuously measuring the heart rate and oxygen saturation (SpO2%). A preliminary experiment was also performed to test the absence of geometric patterns from the anterior surface of a motionless artificial eye. OCT imaging showed a motionless, stable anterior surface of the artificial eye and in dry eye patients. Conversely, in the healthy participants of the study, a bull's eye pattern of the tear film was detected by OCT at the 2nd, 3rd, 4th, 5th, and 6th second after blinking, respectively, in 45%, 60%, 45%, 60%, and 40% of OCT scans during the first session, and in 35%, 65%, 65%, 60%, and 35% of cases in the second session. Overall, a total of 200 OCT scans were performed in normal human population. A significant correlation was found between the novel tear film pattern and heart rate during the first and the second session (p < 0.01) in healthy eyes. Conversely, no correlation was revealed with SpO2%. Intraclass correlation (ICC) analysis for OCT imaging of the tear film revealed a statistically significant reproducibility of the results (ICC = 0.838; p < 0.01), indicating the high level of reliability of the method, independently of heart rate and SpO2% variables. There exists a novel, geometric pattern of the tear film during visual fixation detectable by en-face OCT, which is mainly evident as heart rate increases. Its discovery implies in turn the presence of a specific vibration (or imperceptible motion) of the tear film that, at present, is not recognized and corrected by the OCT software (in image postprocessing) unlike other eyeball movements.

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.

Wavelet representation of the corneal pulse for detecting ocular dicrotism

Purpose

To develop a reliable and powerful method for detecting the ocular dicrotism from non-invasively acquired signals of corneal pulse without the knowledge of the underlying cardiopulmonary information present in signals of ocular blood pulse and the electrical heart activity.

Methods

Retrospective data from a study on glaucomatous and age-related changes in corneal pulsation [PLOS ONE 9(7),(2014):e102814] involving 261 subjects was used. Continuous wavelet representation of the signal derivative of the corneal pulse was considered with a complex Gaussian derivative function chosen as mother wavelet. Gray-level Co-occurrence Matrix has been applied to the image (heat-maps) of CWT to yield a set of parameters that can be used to devise the ocular dicrotic pulse detection schemes based on the Conditional Inference Tree and the Random Forest models. The detection scheme was first tested on synthetic signals resembling those of a dicrotic and a non-dicrotic ocular pulse before being used on all 261 real recordings.

Results

A detection scheme based on a single feature of the Continuous Wavelet Transform of the corneal pulse signal resulted in a low detection rate. Conglomeration of a set of features based on measures of texture (homogeneity, correlation, energy, and contrast) resulted in a high detection rate reaching 93%.

Conclusion

It is possible to reliably detect a dicrotic ocular pulse from the signals of corneal pulsation without the need of acquiring additional signals related to heart activity, which was the previous state-of-the-art. The proposed scheme can be applied to other non-stationary biomedical signals related to ocular dynamics.

Glaucomatous and age-related changes in corneal pulsation shape. The ocular dicrotism

PURPOSE

To ascertain whether the incidence of ocular dicrotic pulse (ODP) increases with age, it is more pronounced in glaucomatous than in normal eyes and whether it is related to cardiovascular activity.

METHODS

261 subjects aged 47 to 78 years were included in the study and classified into four groups: primary open angle glaucoma (POAG), primary angle-closure glaucoma (PACG), glaucoma suspects with glaucomatous optic disc appearance (GODA) and the controls (CG). Additionally, in each group, subjects with ODP were divided into two age subgroups around the median age. A non-contact ultrasonic method was used to measure corneal indentation pulse (CIP) synchronically with the acquisition of electrocardiography (ECG) and blood pulse signals. ODP was assessed from the acquired signals that were numerically processed in a custom written program.

RESULTS

ODP incidence was about 78%, 66%, 66% and 84% for CG, GODA, POAG, and PACG group, respectively. With advancing age, the ODP incidence increased for all subjects (Δ = 12%), the highest being for the PACG and POAG groups (Δ = 30%). GODA group did not show an age-related increase in the incidence of ODP.

CONCLUSIONS

The ocular dicrotism, measured with non-contact ultrasonic method, was found to be a common phenomenon in elderly subjects. The increased ODP incidence in PACG and POAG group may correspond to either higher stiffness of glaucoma eyes, biochemical abnormalities in eye tissues, changes in ocular hemodynamics, may reflect the effect of medications or be a combination of all those factors. The results of GODA group suggest different mechanisms governing their ocular pulse that makes them less susceptible to generating ODP and having decreased predisposition to glaucoma.

Age-related changes in corneal pulsation: ocular dicrotism

PURPOSE

To evaluate age-related changes in the shape of corneal indentation pulse (CIP) in relation to simultaneously registered blood pulsation (BPL) and electrocardiogram (ECG) signals in healthy subjects.

METHODS

Seventy-seven subjects were recruited for the study, including 36 young subjects (aged 23 to 32 years) and 41 older subjects (aged 44 to 72 years) with no reported ocular and cardiovascular abnormalities and free of any systemic diseases. Corneal indentation pulse was measured with a noninvasive ultrasonic distance sensor. Simultaneously, BPL and ECG were acquired with a pulse oximeter placed on the earlobe and Eindhoven triangle setup, respectively. Correlation analysis was applied to identify the interactions between the parameters of the averaged signals, estimated with dedicated signal processing algorithms, for all subject measurements and separately for both considered groups.

RESULTS

Twenty-nine out of 41 older subjects showed evidence of a double peak-shaped CIP waveform that could be interpreted as ocular pulse dicrotism. None of the young subjects exhibited this phenomenon. For the young group, significant correlations (p < 0.05) were found between the CIP crest time of a single peak-shaped CIP waveform and the other shape parameters of the CIP, BPL, and ECG signals. Similarly, for the older group of subjects, significant corresponding correlations were only found between the preliminary crest time of a double peak-shaped CIP waveform.

CONCLUSIONS

Using noninvasive corneal indentation pulse measurement revealed, for the first time, an ocular pulse dicrotism. This might be a natural sign of aging or an early indication of hemodynamic aspects of cardiovascular diseases.

Clinical utility of spectral analysis of intraocular pressure pulse wave

Background

To evaluate the clinical utility of spectral analysis of intraocular pressure pulse wave in healthy eyes of a control group (CG), patients having glaucomatous optic disc appearance or ocular hypertension, and patients with primary open angle glaucoma or primary angle closure glaucoma.

Methods

This is a prospective study that enrolled 296 patients from a single glaucoma clinic. Age matched CG consisted of 62 individuals. Subjects underwent comprehensive clinical diagnostic procedures including intraocular pressure (IOP) measurement with dynamic contour tonometry (DCT) and Goldmann applanation tonometry (GAT). DCT time series were analyzed with custom written software that included signal preprocessing, filtering and spectral analysis. An amplitude and energy content analysis, which takes into account non-stationarity of signals but also provides methodology that is independent of IOP and ocular pulse amplitude (OPA) levels, was applied. Spectral content up to the 6th harmonic of the pressure pulse wave was considered. Statistical analyses included descriptive statistics, normality test, and a multicomparison of medians for independent groups using Kruskal-Wallis test.

Results

GAT IOP showed statistical significance (Kruskal-Willis test p < 0.05) for three out of 10 considered multiple comparisons, DCT IOP and OPA showed statistically significant results in five and seven cases, respectively. Changes in heart rate and central corneal thickness between the groups were statistically significant in two cases. None of the above parameters showed statistically significant differences between CG and the suspects with glaucomatous optic disc appearance (GODA). On the other hand, spectral analysis showed statistically significant differences for that case.

Conclusions

Spectral analysis of the DCT signals was the only method showing statistically significant differences between healthy eyes and those of GODA suspects.

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.

Future directions in non-invasive measurements of tear film surface kinetics

PURPOSE

To review the methods for dynamic, non-invasive, and objective assessment of tear film surface quality and to outline their current state-of-the-art and their future potential.

METHODS

Among the methods available, high-speed videokeratoscopy, lateral shearing interferometry, and dynamic wavefront sensing are being considered.

RESULTS

The principles of operations, their advantages and disadvantages, and limitations of each method are being outlined. The possible future directions of each method are also discussed.

CONCLUSIONS

To gain a better understanding of tear film, its structure and function, it is essential to combine macroimaging technologies with those focusing on tear film microstructure. In this way, one can envisage a clinical device that could help, in future, early diagnosis of dry eye syndrome and development of better materials for contact lenses and eye lubricants.

Phase dependencies between longitudinal corneal apex displacement and cardiovascular signals: is the ocular pulse influenced by the electrical activity of the heart?

BACKGROUND

The aim was to establish phase relationships between the principal harmonic, related to the heart rate, of synchronically registered longitudinal corneal apex displacement (LCAD), blood pulsation (BP) and electrical heart activity signals in a group of healthy subjects.

METHODS

Longitudinal corneal apex displacement was non-invasively measured using an ultrasonic distance sensor. Synchronously, electrocardiographic (ECG) and blood pulsation signals were acquired. As all considered signals are non-stationary (that is, their spectral characteristics vary in time), a reliable and repeatable phase estimation method was sought. For this, a range of phase estimators were tested in the windowed regime of simulated non-stationary signals. Two robust estimators that showed minimum mean square error performance, were selected for further analysis of real signals registered for seven subjects participating in the study.

RESULTS

The windowed cross-correlation and the windowed minimum sum of squared error method achieved the best results among the estimators considered and their outputs were averaged to arrive at a robust phase estimator. Across the subjects, it was found that an increase in the time delay between the principal harmonic of BP and ECG signals, θ(BP,ECG), corresponds to a slight time delay increase between the corresponding harmonics of longitudinal corneal apex displacement and blood pulsation signals, θ(LCAD,BP) and a decrease in the time delay between those of longitudinal corneal apex displacement and ECG signals, θ(LCAD,ECG). Significant correlation (paired t-test, p < 0.05) were found between θ(BP,ECG) and θ(LCAD,BP) as well as between θ(BP,ECG) and θ(LCAD,ECG). There was no significant correlation found between θ(LCAD,BP) and θ(LCAD,ECG).

CONCLUSION

The results indicate that longitudinal corneal apex displacement and correspondingly the ocular pulse phenomenon have not only a vascular origin but could also be influenced by the electrical activity of the heart.

Ultrasonic in vivo measurement of ocular surface expansion

Our aim was to ascertain whether the ultrasonic measurement of longitudinal corneal apex displacements carried out in a proper headrest is a credible method of ocular pulse (OP) detection. To distinguish between longitudinal movements of the eye globe treated as a rigid body and ocular surface expansion caused by the variations of the eye-globe volume, two ultrasound distance sensors were applied to noninvasively measure displacements of cornea and sclera. The same sensors were used to examine the influence of the anterio-posterior movements of a fixed head on the registration of corneal apex pulsation. In both experiments, ECG signals were synchronically recorded. Time, spectral, and coherence analyses obtained for four healthy subjects showed that the ocular surface expansion due to pulsatile ocular blood flow (POBF) is the main component of longitudinal corneal displacement. Ocular surface pulsation is always affected by the head movement. However, there exist some unique properties of signals, which help to distinguish between head and eye movements. A rigid headrest and a bite bar are required to stabilize the head during OP measurement. Ultrasonic technique enables noninvasive and accurate in vivo measurement of corneal pulsation, which could be of interest for indirectly estimating intraocular pressure propagation and POBF component.

Comparison of high-speed videokeratoscopy and ultrasound distance sensing for measuring the longitudinal corneal apex movements

Two different methods to measure binocular longitudinal corneal apex movements were synchronously applied. High-speed videokeratoscopy at a sampling frequency of 15 Hz and a custom-designed ultrasound distance sensor at 100 Hz were used for the left and the right eye, respectively. Four healthy subjects participated in the study. Simultaneously, cardiac electric cycle (ECG) was registered for each subject at 100 Hz. Each measurement took 20 s. Subjects were asked to suppress blinking during the measurements. A rigid headrest and a bite-bar were used to minimize undesirable head movements. Time, frequency and time-frequency representations of the acquired signals were obtained to establish their temporal and spectral contents. Coherence analysis was used to estimate the correlation between the measured signals. The results showed close correlation between both corneal apex movements and the cardiopulmonary system. Unraveling these relationships could lead to better understanding of interactions between ocular biomechanics and vision. The advantages and disadvantages of the two methods in the context of measuring longitudinal movements of the corneal apex are outlined.

Mechanical interferometry of nanoscale motion and local mechanical properties of living zebrafish embryos

We present an interferometric imaging technique that permits local measurement of mechanical properties and nanomechanical motion in small living animals. Measurements of nanomechanical properties and spatially resolved pulsations of <60 nm were recorded for the developing eye of a living zebrafish (Danio rerio) embryo, an important model organism. We also used magnetic microreflectors to conduct contact nanomechanical indentation measurements of the stiffness of the embryonic eye.