Imaging video plethysmography shows reduced signal amplitude in glaucoma patients in the area of the microvascular tissue of the optic nerve head

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.

A retinal imaging system for combined measurement of optic nerve head vascular pulsation and stimulated vasodilation in humans

Vascular pulsation at the optic nerve head (ONH) reflects vessel properties. Reduction in the stimulated retinal vasodilatory capacity has been reported in diabetes, but its relation with vascular pulsation is unknown. Here we report a new retinal imaging system for correlative assessment of ONH vascular pulsation and stimulated retinal vasodilation. Retinal reflectance images were acquired before and during light flicker stimulation to quantify arterial and venous vasodilation (DAR, DVR) in subjects with and without diabetic retinopathy (N = 25). ONH vascular pulsation amplitude and frequency (PA, PF), were quantified by curve fitting of periodic intensity waveforms acquired in retinal vasculature (RV) and ONH tissue (ONHT) regions. The relationships between pulsation metrics, heart rate (HR), intraocular pressure (IOP), and vasodilatory responses were evaluated. Pulsation metrics were not significantly different between regions (p ≥ 0.70). In RV, inter-image variabilities of PA and PF were 10% and 6%, whereas inter-observer variabilities were 7% and 2% respectively. In both regions, PF was correlated with HR (p ≤ 0.001). PA was associated with DAR in both regions (p ≤ 0.03), but only with DVR in RV (p ≤ 0.05). Overall, ONH vascular pulsation was associated with stimulated retinal vasodilation, suggesting diabetes may have concomitant effects on retinal vasculature compliance and neurovascular coupling.

A multi-color video-ophthalmoscopes allows to measure the spectral distribution of light absorption of blood in the human retina

Based on our previously developed mono-color video-ophthalmoscope a multi-color video-ophthalmoscope was developed. Using narrow band transmission filters, this instrument allows to measure the pulsatile cardiac cycle induced blood volume changes in the human retina for any wavelength in the sensitivity range of the used CMOS-camera. In this key experiment, video sequences (8 s, 25 fps, 200 frames) of the optic nerve head (ONH) were acquire for seven wavelengths between 475 nm and 677 nm one after the other. After image registration of all frames of each video sequence (to compensate for eye movements) and trend correction (to compensate for slow intensity changes), the amplitude of the cardiac cycle induced light intensity changes (pulsatile absorption amplitude PAA) can be calculated for all seven wavelengths. The results confirmed that the spectral distribution of PAA (λ) follows the distribution of the light absorption of blood. The measured values correspond to the absorption of a thin blood layer of about 0.5 μm thickness.

Comparative analysis of retinal photoplethysmographic spatial maps and thickness of retinal nerve fiber layer

The paper presents a comparative study of the pulsatile attenuation amplitude (PAA) within the optic nerve head (ONH) at four different areas calculated from retinal video sequences and its relevance to the retinal nerve fiber layer thickness (RNFL) changes in normal subjects and patients with different stages of glaucoma. The proposed methodology utilizes processing of retinal video sequences acquired by a novel video ophthalmoscope. The PAA parameter measures the amplitude of heartbeat-modulated light attenuation in retinal tissue. Correlation analysis between PAA and RNFL is performed in vessel-free locations of the peripapillary region with the proposed evaluating patterns: 360° circular area, temporal semi-circle, nasal semi-circle. For comparison, the full ONH area is also included. Various positions and sizes of evaluating patterns in peripapillary region were tested which resulted in different outputs of correlation analysis. The results show significant correlation between PAA and RNFL thickness calculated in proposed areas. The highest correlation coefficient Rtemp = 0.557 (p<0.001) reflects the highest PAA-RNFL correspondence in the temporal semi-circular area, compared to the lowest value in the nasal semi-circular area (Rnasal = 0.332, p<0.001). Furthermore, the results indicate the most relevant approach to calculate PAA from the acquired video sequences is using a thin annulus near the ONH center. Finally, the paper shows the proposed photoplethysmographic principle based on innovative video ophthalmoscope can be used to analyze changes in retinal perfusion in peripapillary area and can be potentially used to assess progression of the RNFL deterioration.

Multispectral retinal video-ophthalmoscope with fiber optic illumination

Multispectral imaging is used in various applications including astronomy, industry and agriculture. In retinal imaging, the single-shot multispectral image stack is typically acquired and analyzed. This multispectral analysis can provide information on various structural or metabolic properties. This paper describes the multispectral improvement of a video-ophthalmoscope, which can acquire retinal video sequences of the optic nerve head and peripapillary area using various spectral light illumination. The description of the multispectral video imaging is provided and several applications are described. These applications include multispectral retinal photoplethysmography, visualization of spontaneous vein pulsation and multispectral RGB image generation.

Heart rate and age modulate retinal pulsatile patterns

Theoretical models of retinal hemodynamics showed the modulation of retinal pulsatile patterns (RPPs) by heart rate (HR), yet in-vivo validation and scientific merit of this biological process is lacking. Such evidence is critical for result interpretation, study design, and (patho-)physiological modeling of human biology spanning applications in various medical specialties. In retinal hemodynamic video-recordings, we characterize the morphology of RPPs and assess the impact of modulation by HR or other variables. Principal component analysis isolated two RPPs, i.e., spontaneous venous pulsation (SVP) and optic cup pulsation (OCP). Heart rate modulated SVP and OCP morphology (p_FDR < 0.05); age modulated SVP morphology (p_FDR < 0.05). In addition, age and HR demonstrated the effect on between-group differences. This knowledge greatly affects future study designs, analyses of between-group differences in RPPs, and biophysical models investigating relationships between RPPs, intracranial, intraocular pressures, and cardiovascular physiology.

Video-recordings of retinal pulsatile patterns in human eye are modulated by heart rate and age. This invivo evidence fundamentally impacts modeling of retinal function and clinical design in wide range of medical specialties.

Photoplethysmographic analysis of retinal videodata based on the Fourier domain approach

Assessment of retinal blood flow inside the optic nerve head (ONH) and the peripapillary area is an important task in retinal imaging. For this purpose, an experimental binocular video ophthalmoscope that acquires precisely synchronized video sequences of the optic nerve head and peripapillary area from both eyes has been previously developed. It enables to compare specific characteristics of both eyes and efficiently detect the eye asymmetry. In this paper, we describe a novel methodology for the analysis of acquired video data using a photoplethysmographic approach. We describe and calculate the pulsatile attenuation amplitude (PAA) spatial map, which quantifies the maximum relative change of blood volume during a cardiac cycle using a frequency domain approach. We also describe in detail the origin of PAA maps from the fundamental (the first) and the second harmonic component of the pulsatile signal, and we compare the results obtained by time-based and frequency-based approaches. In several cases, we show the advantages and possibilities of this device and the appropriate image analysis approach - fast measurement and comparison of blood flow characteristics of both eyes at a glance, the robustness of this approach, and the possibility of easy detection of asymmetry.

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.