Changes in Retinal Optical Coherence Tomography Angiography Indexes Over 24 Hours.

Purpose: To evaluate changes in the retinal microvasculature of young adults over 24 hours using optical coherence tomography angiography (OCT-A).
Methods: Participants (n = 44, mean age 23.2 ± 4.1 years, 24 myopes and 20 nonmyopes) with normal ophthalmological findings were recruited. Two macular OCT-A and OCT scans, systemic blood pressure, intraocular pressure (IOP), and biometry measurements were taken every four hours over 24 hours. Superficial and deep retinal layer en face images were analyzed to extract magnification-corrected vascular indexes using image analysis including foveal avascular zone metrics, vessel density, and perfusion density for the foveal, parafoveal, and perifoveal zones.
Results: Significant diurnal variations (P < 0.001) were observed in the vessel and perfusion density in the three superficial retinal layer regions, with acrophase between 4:30 PM and 8:30 PM. Only foveal and parafoveal regions of the deep retinal layer exhibited significant diurnal variations with acrophase between 9 AM and 3 PM. Myopes and nonmyopes had different acrophases but not amplitudes in the parafoveal perfusion density of superficial retinal layer (P = 0.039). Significant correlations were observed between diurnal amplitudes or acrophases of superficial retinal layer indexes and systemic pulse pressure, IOP, axial length and retinal thickness. Conclusions: This study shows, for the first time, that significant diurnal variation exists in OCT-A indexes of macular superficial and deep retinal layer over 24 hours and were related to variations in various ocular and systemic measurements. Myopes and nonmyopes showed differences in the timing but not in amplitude of the superficial retinal layer parafoveal perfusion density variations but not in deep retinal layer.

Diurnal variation in retinal microvasculature have been relatively unexplored, despite evidence that ocular diurnal variation plays an important role in the regulation of eye growth and pathogenesis of ocular diseases such as glaucoma., Optical coherence tomography angiography (OCT-A), a relatively new technique, allows rapid noninvasive imaging of retinal microvasculature at different depths such as the superficial and deep retinal layer and can provide indexes such as vessel and perfusion density that show excellent repeatability.– The superficial retinal layer is the innermost region and extends from the internal limiting membrane to the inner plexiform layer whereas the deep retinal layer extends from the inner nuclear layer to the outer plexiform layer.

Past studies have explored daytime variations in the OCT-A indexes of superficial and deep retinal layers among healthy participants between 7 AM and 8 PM (Penteado R, et al. IOVS 2018;59:ARVO E-Abstract 2856; Rommel F, et al. IOVS 2018;59:ARVO E-Abstract 2860).– Most of these studies did not detect significant variations during the day, except for one study that reported a significant increase (P = 0.048) in parafoveal perfusion density between 12 PM and 8 PM using linear spline model analysis (Penteado R, et al. IOVS 2018;59:ARVO E-Abstract 2856). Other ocular measurements also showed significant diurnal variations. Retinal thickness, axial length, and intraocular pressure (IOP) demonstrate a reduction at night whereas choroidal thickness and mean ocular perfusion pressure (MOPP) demonstrate peaks during the night or early morning.– Additionally, significant fluctuations exists in systemic blood pressure, vascular resistance, and heart rate during nighttime., These ocular and systemic measurements also show significant associations with OCT-A indexes and contribute to the regulation of ocular blood flow.– Therefore it is expected that retinal OCT-A indexes may also demonstrate significant variations over 24 hours. However, to the best of our knowledge, no previous studies have evaluated OCT-A measurements over a complete 24-hour period. In addition, OCT-A indexes are also affected by axial length and refractive error.– Myopes had lower superficial and deep vessel/perfusion density and higher vascular resistance index than emmetropes in a number of previous studies.,,,– However, the role of these factors in diurnal variation of OCT-A indexes over 24 hours is not known.

This study therefore assessed the magnitude and pattern of variation in the microvasculature of superficial and deep retinal layers by capturing a series of macular angiography scans using OCT-A every four hours over 24 hours in young healthy participants. Findings were also compared between myopic and nonmyopic participants. Diurnal variations in axial length, IOP, retinal thickness and systemic blood pressure were also assessed to investigate their relationship with OCT-A indexes. This improved knowledge about diurnal variation will help to establish the expected physiological changes over a 24-hour period in young healthy adults and to assist in the interpretation of clinical changes in OCT-A indexes.

Methods

Participants

This was a prospective observational study conducted among 44 participants including 19 females and 25 males (mean age 23.2 ± 4.1 years, range 18–35 years). Older participants were not included to control for age-related changes in retinal microvasculature.– Participants with best-corrected vision of 6/7.5 or better, noncycloplegic refraction between ± 6.00 DS and astigmatism ≤2.00 DC, normal ophthalmological findings, normal sleep/wake cycle, and good sleep quality (confirmed with the Pittsburgh Sleep Quality Index [PSQI] questionnaire, PSQI score ≤7) were included. Participants were screened a week ± 2 days before their diurnal measurement session through a detailed medical and ocular history, visual acuity assessment, refraction, slit-lamp examination, IOP measurement, and OCT fundus examination. Noncycloplegic autorefraction (Tomey RC-800; Tomey Corporation, Nayoga, Japan) was done followed by careful subjective refraction with blur back and duochrome technique ensuring relaxed accommodation and aiming for maximum plus to achieve best visual acuity. Participants were asked to refrain from caffeine, alcohol, or tobacco use, vigorous physical activity and abstain from contact lens wear (i.e., to wear spectacles if required) on the experiment day. The Queensland University of Technology Human Research Ethics Committee approved this study and written informed consent was obtained after providing clear explanation of the testing procedures. All participants were treated in accordance with the tenets of the declaration of Helsinki. Participants were classified based upon the subjective, noncycloplegic spherical equivalent refraction (SER) of their right eye as being myopic (SER ≤ −0.50D) or nonmyopic (SER > −0.50D) and included 24 myopes and 20 nonmyopes. Table 1 summarizes the participants demographic information.

Table 1.

Demographics of the Participants

	Total	Myopes	Non-Myopes
Participants (Female/Male)	N = 44 (19/25)	N = 24 (8/16)	N = 20 (11/9)
Age (y), mean ± SD (range)	23.2 ± 4.1 (19.9–35)	23.5 ± 4.3 (19.9–35)	23.0 ± 4.0 (19.9–34.4)
SER (Diopter)	−1.26 ± 1.92 (−7.00 to 1.25)	−2.41 ± 1.95 (−0.50 to −7.00)	0.11 ± 0.38 (−0.38 to +1.25)
Astigmatism (Diopter)	0.40 ± 0.48 (0.00 to 2.00)	0.52 ± 0.56 (0.00 to 2.00)	0.26 ± 0.32 (0.00 to 0.75)

SER, spherical equivalent refraction; SD, standard deviation.

Age (P = 0.65) and gender distribution (P = 0.22) was not significantly different between the two refractive groups.

Data Collection Protocol

For eligible participants, measurements were collected at QUT Optometry Clinic over an approximate 24-hour period including seven sessions starting around 9 AM and ending around 9 AM the following day. Before each measurement session (excluding 1 AM and 5 AM), 10 minutes wash-out period was provided (watching television at 4 meters) with their habitual spectacles for relaxing accommodation and standardizing the conditions to limit the influence of previous activities upon the measurements. Participants carried out their normal daily activities and were permitted to leave the clinic between 9 AM and 11 PM. They were provided with bedding and remained in the clinic overnight with lights off from 11 PM to 7 AM and were encouraged to sleep. For the measurements at night (1 AM and 5 AM), the room illumination was lowered (∼2 Lux) to reduce the impact of measurements on participants’ normal circadian rhythm. The following systemic and ocular measurements (only on the right eye) in this order: OCT-A and OCT scans, blood pressure, IOP, and biometry measurements and required 15 to 20 minutes to complete.

OCT-A and OCT Scanning Protocol and Image Analysis

Two 3 mm × 3 mm and 6 mm × 6 mm macular angiography enhanced depth imaging scans each were captured through nondilated pupils in dim illumination using Zeiss OCT-A (AngioPlex software, version 11.0;Cirrus HD-OCT 5000; Carl Zeiss Meditec Inc, Dublin, CA, USA) with FastTrac Retinal Tracking. To ensure image quality, the OCT-A images were recaptured if signal strength was <9 or if they had severe motion artifacts, were not centered on the fovea or exhibited scan tilt > 5°., A larger scan size effectively reduces the image resolution, decreasing the ability to pick up fine capillaries in an image.– So, the foveal and parafoveal vascular indexes were quantified from 3 mm scans for detailed assessment and perifoveal indexes from 6 mm scans. The superficial retinal layer (internal limiting membrane to inner plexiform layer) and deep retinal layer (inner nuclear layer to outer plexiform layer) images were exported for extraction of the vascular indexes.

The details of the image analysis and extraction of indexes have been provided in the supplementary section 1. Using a custom image analysis MATLAB program (MathWorks, Natick, MA, USA), the superficial and deep retinal layer images from the seven sessions were registered with each other to ensure that OCT-A indexes were obtained from the same retinal region within the participants. A magnification correction factor determined using a previously described schematic eye method incorporating individual biometry measures (axial length, corneal power, anterior chamber depth and refractive error) was then used to rescale the foveal avascular zone (FAZ) and extract its measurements. Then, the superficial and deep retinal layer images were binarized and skeletonized from which the OCT-A indexes were derived (See Supplementary Section 1). A modified Early Treatment Diabetic Retinopathy Study grid adjusted for magnification which ensured that the indexes between subjects were derived from the same sized retinal regions. The grid was then centered on the images to extract OCT-A indexes from the foveal (central 1 mm), parafoveal (2.5 mm) and perifoveal (5 mm) regions of the superficial and deep retinal layer images (Fig. 1).

Figure 1.

Example of the modified early treatment diabetic retinopathy study (ETDRS) grid for the 6 mm × 6 mm (left) and 3 mm × 3 mm (right) superficial retinal layer (SRL) and deep retinal layer (DRL) en face image. Foveal avascular zone: highlighted in green in the center; central/foveal zone (1 mm): inside yellow circle; inner/parafoveal zone (2.5 mm annulus) with quadrants: region between yellow and green circles (left); and outer/perifoveal zone (5 mm annulus) with quadrants: region between green and blue circles (right).

The OCT-A indexes extracted from the superficial and deep retinal layers were the FAZ area, FAZ perimeter and the perfusion density and vessel density. The average of the values from the two images captured at each session were taken for analysis. The FAZ area (mm2) and FAZ perimeter (mm) were extracted from the superficial and deep retinal layer images of the 3 mm scan as FAZ metrics are not significantly affected by scan size., Perfusion density is calculated in percentage as the ratio of image area occupied by the vasculature to the total measured area in the binarized OCT-A image. Perfusion density takes both vessel length and diameter into consideration. Vessel density is calculated in mm/mm2 as the ratio of the length occupied by the blood vessels to the total measured area in the skeletonized OCT-A image. Vessel density indicates the presence of vessels by quantifying the vessel length density irrespective of the vessel diameters. Large vessels and small capillaries contribute equally to the vessel density quantification. As vessel density has a linear component in it, an additional transverse magnification correction (Equation 1) is required to appropriately adjust the extracted vessel density values for magnification. These corrected values were then used for further analysis.

A macular raster enhanced depth imaging–OCT scan was captured twice and analyzed using previously described automated methods in MATLAB to segment different retinal layers (Fig. 2). Magnification corrected total retinal, superficial and deep retinal layer thicknesses of the foveal, parafoveal and perifoveal regions were used for analysis (See Supplementary Section 1).

Figure 2.

Example of the segmented boundaries and the analysis zone of a typical OCT image. Segmentation boundaries included internal limiting membrane (red), inner nuclear layer/inner plexiform layer boundary (blue), outer plexiform layer/outer plexiform layer boundary (yellow), and outer boundary of retinal pigment epithelium (green). Total retinal thickness is the axial distance from the retinal pigment epithelium to the inner limiting membrane. Superficial retinal layer thickness is the axial distance from the inner nuclear layer/inner plexiform layer boundary to the inner limiting membrane and deep retinal layer thickness is the axial distance from the inner nuclear layer/inner plexiform layer boundary to the outer plexiform layer/outer plexiform layer boundary. Data were binned into the fovea (1 mm), parafovea (2.5 mm), and perifovea (5 mm).

Ocular Biometry, IOP and Blood Pressure Measurements

Ocular biometrics including axial length, corneal curvature, and anterior chamber depth were measured using IOL master (version 5.4; Carl Zeiss, Jena, Germany). A set of five measurements were captured and averaged. Two blood pressure readings including systolic and diastolic blood pressure were taken using an automated sphygmomanometer (version 42NOB; Welch Allyn, Skaneateles Falls, New York, USA) and were averaged for further analysis. IOP was measured using a noncontact tonometer (Canon TX-10 noncontact tonometer; Canon Inc., Tokyo, Japan). An average of three readings of IOP were taken for analysis. The pulse pressure, mean arterial pressure (MAP) and MOPP were calculated for each session from the following equations.,

Statistical Analysis

Statistical analysis was performed using SPSS (IBM Corp, Armonk, NY, USA), Microsoft Excel, and SigmaPlot (Systat Software, San Jose, CA) software. Repeated measures analysis of variance (ANOVA) was used to investigate the effects of time of day and refractive group for each of the measured parameters along with Bonferroni adjusted pairwise comparison. However, if the amplitudes and acrophase do not show any large within-group variations, a repeated measures ANOVA may not be sensitive enough to detect the diurnal changes and circular statistics can be more appropriate. Using circular statistical tests provides a sensitive test for detecting the occurrence of acrophase and its variation between the groups. Therefore Fourier analysis was conducted to evaluate the acrophase (peak time) and diurnal amplitude for each subject independently.,,

The Rayleigh test was used to assess the dispersion of acrophases across the day. Rayleigh's r values ranges between 0 and 1 with higher values representing more significant clustering of the acrophases indicating significant diurnal variations. The Rayleigh test has been reported to possess very good control of type 1 error rate., Mann-Whitney U and Watson's U2 test were performed to investigate refractive group differences in diurnal amplitudes and acrophase respectively. Spearman's Rho correlation (rρ) and circular correlation (rcc) assessed the relationship between amplitude and acrophases of different measurements respectively. The measurements were regarded as in-phase if their acrophases occurred at the same time (±2 hours) and antiphase if their acrophases occurred 12 hours apart (±2 hours).

Results

Significant diurnal variations were found in a number of the superficial and deep retinal layer OCT-A indexes. Tables 2 and 3 show the daily mean, diurnal amplitudes, diurnal acrophases, Rayleigh test results and repeated measures ANOVA outcomes for the OCT-A indexes of the superficial and deep retinal layers for all participants.

Table 2.

Daily Mean, Amplitude, and Acrophase of Diurnal Variation and Its Dispersion, and Time of Day and Refractive Group Comparisons for the OCT-A Indexes in the Foveal, Parafoveal, and Perifoveal Zones of the Superficial and Deep Retinal Layer En Face Images for All Participants

	Daily Mean Mean ± SEM	Diurnal Amplitude Mean ± SEM* (%)	Acrophase, Peak Time ± CD (Hours)	Rayleigh's r, P Value	Time of Day (df = 6, 258)†	Time by Refractive Error (df = 6, 252)†
SRL-FAZ area (mm2)	0.26 ± 0.01	0.02 ± 0.01 (7.69%)	10:16 AM ± 4.66	0.11, >0.09	0.24	0.31
SRL-FAZ perimeter (mm)	2.11 ± 0.05	0.13 ± 0.01 (6.16%)	9:48 AM ± 5.06	0.25, 0.05	0.66	0.66
DRL-FAZ area (mm2)	0.68 ± 0.02	0.08 ± 0.01 (11.76%)	3:47 PM ± 4.61	0.42, <0.001‡	0.002‡	0.57
DRL-FAZ perimeter (mm)	3.20 ± 0.06	0.22 ± 0.02 (6.88%)	4:14 PM ± 4.07	0.43, <0.001‡	0.02‡	0.60
SRL-perfusion density (%)
Fovea	28.13 ± 0.99	1.68 ± 0.13 (6.01%)	5:29 PM ± 4.68	0.26, 0.05‡	0.33	0.38
Parafovea	48.62 ± 0.26	1.55 ± 0.14 (3.19%)	6:13 PM ± 3.98	0.46, <0.001‡	0.001‡	0.04‡
Perifovea	34.98 ± 0.11	0.65 ± 0.06 (1.85%)	6:55 PM ± 3.90	0.48, <0.001‡	0.004‡	0.73
SRL-vessel density (mm/mm2)
Fovea	8.15 ± 0.28	0.61 ± 0.05 (7.59%)	4:50 PM ± 4.55	0.48, <0.001‡	0.045‡	0.42
Parafovea	14.35 ± 0.15	0.60 ± 0.06 (4.45%)	5:53 PM ± 4.13	0.41, <0.001‡	0.005‡	0.82
Perifovea	14.15 ± 0.01	0.46 ± 0.04 (3.29%)	7:13 PM ± 4.57	0.28, 0.04‡	0.20	0.64
DRL-perfusion density (%)
Fovea	20.92 ± 0.95	3.88 ± 0.43 (18.55%)	2:58 PM ± 4.68	0.26, 0.05‡	0.005‡	0.65
Parafovea	50.05 ± 0.55	1.92 ± 0.15 (3.83%)	11:17 PM ± 4.54	0.29, 0.03‡	0.08	0.67
Perifovea	38.72 ± 0.48	0.54 ± 0.05 (1.39%)	9:27 PM ± 5.17	0.08, >0.9	0.91	0.76
DRL-vessel density (mm/mm2)
Fovea	7.16 ± 0.29	1.14 ± 0.12 (16.24%)	2:28 PM ± 4.64	0.27, 0.045‡	0.002‡	0.53
Parafovea	13.74 ± 0.19	0.67 ± 0.05 (4.92%)	9:29 PM ± 4.49	0.31, 0.01‡	0.09	0.72
Perifovea	16.73 ± 0.21	0.37 ± 0.04 (2.24%)	10:24 PM ± 5.13	0.10, 0.67	0.54	0.76

SEM, standard error of mean; CD, circular deviation; SRL, superficial retinal layer; DRL, deep retinal layer.

Diurnal amplitude expressed as the percentage of the daily mean value.

P values from repeated measures ANOVA for time of day, and time of day by refractive error.

P < 0.05.

Table 3.

Amplitude of Diurnal Variation, Acrophase of Diurnal Variation and its Dispersion, and Time of Day and Refractive Group Comparisons for Vessel and Perfusion Density in the Parafoveal and Perifoveal Zone Quadrants of the Superficial and Deep Retinal Layer En Face Images for all Participants

	Diurnal Amplitude Mean ± SEM	Acrophase, Peak Time ± CD (Hours)	Rayleigh's r, P Value	Time of Day (df = 6, 258)*	Time by Refractive Error (df = 6, 252)*
SRL-perfusion density (%)
Parafovea-superior	2.21 ± 0.20	7:27 PM ± 4.59	0.28, 0.04†	0.13	0.51
Parafovea-nasal	1.65 ± 0.16	6:55 PM ± 4.31	0.36, 0.001†	0.01†	0.04†
Parafovea-inferior	2.09 ± 0.19	4:23 PM ± 3.89	0.48, <0.001†	0.001†	0.93
Parafovea-temporal	1.68 ± 0.0.13	6:16 PM ± 4.16	0.41, <0.001†	0.003†	0.36
Perifovea-superior	0.98 ± 0.10	8:22 PM ± 4.30	0.37, 0.001†	0.02†	0.43
Perifovea-nasal	0.90 ± 0.07	5:52 PM ± 4.31	0.36, 0.001†	0.05†	0.38
Perifovea-inferior	0.91 ± 0.7	6:19 PM ± 4.49	0.31, 0.01†	0.24	0.96
Perifovea-temporal	0.83 ± 0.09	6:33 PM ± 4.09	0.43, <0.001†	0.08	0.55
SRL-vessel density (mm/mm2)
parafovea-superior	0.80 ± 0.08	6:57 PM ± 4.69	0.26, 0.05†	0.18	0.91
parafovea-nasal	0.65 ± 0.06	5:48 PM ± 4.36	0.35, 0.002†	0.05†	0.19
parafovea-inferior	0.73 ± 0.06	5:02 PM ± 3.79	0.26, 0.03†	0.001†	0.95
parafovea-temporal	0.65 ± 0.06	6:26 PM ± 4.32	0.36, 0.001†	0.01†	0.65
perifovea-superior	0.65 ± 0.06	8:26 PM ± 4.83	0.20, 0.14	0.10	0.34
perifovea-nasal	0.63 ± 0.05	6:40 PM ± 4.69	0.25, 0.06	0.47	0.65
perifovea-inferior	0.61 ± 0.05	6:56 PM ± 4.69	0.24, 0.06	0.56	0.61
perifovea-temporal	0.57 ± 0.06	8:43 PM ± 4.75	0.23, 0.09	0.75	0.35
DRL-perfusion density (%)
Parafovea-superior	3.08 ± 0.28	7:49 PM ± 4.51	0.30, 0.01†	0.07	0.12
Parafovea-nasal	3.07 ± 0.29	11:41 PM ± 4.15	0.41, <0.001†	0.03†	0.82
Parafovea-inferior	3.70 ± 0.34	2:55 AM ± 5.20	0.07, <0.9	0.71	0.68
Parafovea-temporal	2.84 ± 0.24	1:46 AM ± 4.99	0.15, 0.36	0.13	0.51
Perifovea-superior	0.71 ± 0.06	2:02 AM ± 4.87	0.19, 0.20	0.24	0.72
Perifovea-nasal	0.84 ± 0.08	6:29 PM ± 5.15	0.09, >0.90	0.86	0.89
Perifovea-inferior	0.76 ± 0.08	7:59 AM ± 5.12	0.10, 0.67	0.69	0.18
Perifovea-temporal	0.71 ± 0.07	6:29 PM ± 5.11	0.11, 0.61	0.33	0.79
DRL-vessel density (mm/mm2)
Parafovea-superior	0.99 ± 0.09	7:49 PM ± 4.46	0.32, 0.005†	0.04†	0.41
Parafovea-nasal	0.99 ± 0.08	10:18 PM ± 4.38	0.34, 0.004†	0.08	0.79
Parafovea-inferior	0.98 ± 0.09	7:44 PM ± 5.10	0.11, 0.61	0.43	0.77
Parafovea-temporal	0.89 ± 0.08	10:29 PM ± 5.09	0.11, 0.61	0.10	0.44
Perifovea-superior	0.49 ± 0.04	11:23 PM ± 4.65	0.25, 0.06	0.07	0.89
Perifovea-nasal	0.57 ± 0.06	4:34 PM ± 4.84	0.20, 0.14	0.60	0.37
Perifovea-inferior	0.51 ± 0.05	5:38 AM ± 5.24	0.06, >0.9	0.76	0.27
Perifovea-temporal	0.48 ± 0.05	2:10 AM ± 5.38	0.01, >0.9	0.34	0.53

SEM, Standard Error of Mean; CD, Circular Deviation; df, degree of freedom; SRL, superficial retinal layer; DRL, deep retinal layer.

P values from repeated measures ANOVA for time of day and time of day by refractive error.

P < 0.05.

Diurnal Variation in the Superficial and Deep Retinal Layer OCT-A Indexes

Significant variations over 24 hours were observed in superficial retinal layer foveal and parafoveal vessel density (P = 0.045 and P = 0.005, Amplitude:0.62 mm/mm2 ± 0.32 and 0.61 mm/mm2 ± 0.39) and parafoveal and perifoveal perfusion density (P = 0.001 and P = 0.037, Amplitude: 1.55% ± 0.92 and 0.65% ± 0.42), respectively. Post hoc analysis with a Bonferroni adjustment revealed that foveal vessel and perfusion density at 5 AM were significantly different from 1 PM (P = 0.048 and P = 0.006) and 9 PM (P = 0.034 and P = 0.002) and parafoveal perfusion density at 5 AM was significantly different from 1 PM (P = 0.04), 5 PM (P = 0.01) and 9 PM (P = 0.004). The variation in perifoveal perfusion density was no longer significant after Bonferroni adjustment.

The Rayleigh test revealed significant acrophase clustering for vessel and perfusion density in all zones indicating significant diurnal variations across the macular regions. Peak superficial retinal layer vessel and perfusion density were observed between 4:30 PM and 8:30 PM (Table 2) for all the zones. The FAZ area and perimeter at superficial retinal layer did not demonstrate significant variations (all P > 0.05).

Repeated measures ANOVA showed significant diurnal variations in vessel and perfusion density for the superficial retinal layer parafoveal quadrants except superior whereas only perifoveal superior and nasal perfusion density showed significant results (Table 3). However, significant clustering of acrophase in the perfusion density of parafoveal and perifoveal quadrants and vessel density of the parafoveal quadrants were observed (Rayleigh's r, P < 0.05 for all).

Only foveal vessel and perfusion density of the deep retinal layer exhibited significant diurnal variations between the sessions (P = 0.005 and P = 0.002) which was no longer significant after Bonferroni adjustment. However, both deep retinal layer foveal and parafoveal vessel and perfusion density showed clustering of acrophase (Table 2) indicating significant diurnal variations. The deep retinal layer parafoveal vessel density and perfusion density peaks were observed between 9 PM and 11:30 PM and earlier in the mid-afternoon for the foveal zone. The FAZ area and perimeter at deep retinal layer showed significant diurnal variation (time of day effect, P = 0.002 and 0.019 and Rayleigh's r = 0.42 and 0.43, P < 0.001) with acrophase at 3:47 PM and 4:14 PM. Post hoc analysis revealed significant differences between FAZ measurements at 1 PM and 1 AM and between measurements at 9 PM and 9 AM (P < 0.05 for all). Only the perfusion density of the superior parafoveal zone demonstrated significant diurnal variations (P = 0.026), whereas the perifoveal quadrants did not exhibit significant variations.

Time by Refractive Groups Interaction in the Superficial and Deep Retinal Layer OCT-A Indexes

Myopes and nonmyopes showed different pattern of variation in the superficial retinal layer parafoveal perfusion density i.e., a significant time of day by refractive group interaction (P = 0.04). Further analysis showed a significant difference (P = 0.02) in the acrophase of superficial retinal layer parafoveal perfusion density between myopes (5:10 PM ± 3.31 hours) and nonmyopes (8:57 PM ± 4.36 hours). Figure 3 demonstrates the mean changes in superficial retinal layer parafoveal vessel and perfusion density changes over 24 hours for the two refractive groups. Figure 4 depicts polar diagrams for parafoveal vessel and perfusion density of superficial retinal layer plotting acrophase against amplitude for individual myopes and nonmyopes. A significant time of day by refractive group interaction (P = 0.04) was found in the parafoveal nasal perfusion density with a significant difference in acrophase of superficial retinal layer nasal perfusion density (P = 0.02) between myopes (4:54 PM ± 3.26 hours) and nonmyopes (8:35 PM ± 4.50 hours).

Figure 3.

Perfusion density (top) and vessel density (bottom) of the superficial retinal layer (SRL) at parafovea in myopes and nonmyopes plotted across 24 hours. The error bars represent one standard error of the mean. Significant changes over 24 hours were found (P < 0.05) for both parafoveal superficial retinal layer vessel and perfusion density. Both indexes were lowest at 5 AM for the two refractive groups. A significant time by refractive group interaction (P = 0.04) was observed only for parafoveal perfusion density.

Figure 4.

Diurnal amplitude of change (distances from center) plotted against acrophase time (around the perimeter, 24-hour clock) for myopes and nonmyopes for superficial retinal layer (SRL) vessel density (right) and perfusion density (left). Mean acrophase for myopes and nonmyopes are indicated by the time in bold. The acrophase of superficial retinal layer parafoveal perfusion density in myopes was significantly (P = 0.02) earlier than nonmyopes whereas the vessel density acrophase difference did not reach statistical significance.

None of the indexes for either of the layers showed any significant difference (all P > 0.05) in diurnal amplitude between the two refractive groups for any of the considered zones. The daily mean vessel density (mm/mm2) in the parafoveal and perifoveal zones of the superficial retinal layer (Myopes: 13.32 ± 0.92 and 13.86 ± 0.80; Non-myopes: 14.08 ± 0.76 and 14.94 ± 0.92) and deep retinal layer (Myopes: 13.37 ± 1.39 and 14.18 ± 0.98; Non-myopes: 16.12 ± 0.76 and 17.46 ± 1.52) were significantly greater in the nonmyopes (P < 0.05 for all). Similar results were also noted in the daily mean vessel density of all the perifoveal and parafoveal quadrants (P < 0.05 for all). The remaining indexes in each of the layers and zones did not show any significant refractive group effects (P > 0.05 for all).

Diurnal Rhythm in Ocular and Systemic Measurements

Axial length, total retinal, superficial, and deep retinal layer thicknesses, IOP, MAP, pulse pressure and MOPP all demonstrated significant diurnal variations (Table 4). The total retinal and superficial retinal layer thickness exhibited significant time of day effects and time by refractive group interactions at the parafovea and perifovea respectively whereas deep retinal layer thickness demonstrated time of day effects only at parafovea and perifovea. A significant effect of refractive group was only noted in axial length, due to the well-documented longer axial length associated with myopia.

Table 4.

Amplitude of Diurnal Variation, Acrophase of Diurnal Variation and Its Dispersion, and Time of Day and Refractive Group Comparisons for Ocular and Systemic Measurements for All Participants

	Diurnal Amplitude Mean ± SEM	Acrophase, Peak Time ± CD (Hours)	Rayleigh's r, P Value	Time of Day (df = 6, 258)*	Time by Refractive Error (df = 6, 252) *
Axial Length (mm)	0.04 ± 0.01	2:14 PM ± 3.17	0.65, <0.001†	<0.001†	0.32
Total Retinal Thickness (µm)
Fovea	5.89 ± 0.36	3:44 PM ± 2.08	0.85, <0.001†	<0.001†	0.84
Parafovea	3.39 ± 0.33	1:14 PM ± 4.40	0.34, 0.003†	0.001†	0.69
Perifovea	3.57 ± 0.60	5:02 AM ± 4.58	0.28, 0.04†	0.02†	0.04†
SRL Thickness (µm)
Fovea	4.02 ± 0.31	6:24 AM ± 5.19	0.08, 0.900	0.77	0.50
Parafovea	2.46 ± 0.21	4:05 AM ± 4.19	0.40, <0.001†	<0.001†	<0.001†
Perifovea	2.87 ± 0.30	2:50 AM ± 4.86	0.19, 0.20	0.12	0.34
DRL Thickness (µm)
Fovea	4.95 ± 0.47	11.28 AM ± 5.07	0.12, 0.65	0.44	0.66
Parafovea	4.93 ± 0.58	12.31 PM ± 4.59	0.28, 0.04†	0.04†	0.47
Perifovea	2.59 ± 0.36	5.40 AM ± 3.82	0.50, <0.001†	0.05†	0.53
Pulse Pressure (mm Hg)	7.04 ± 0.75	5:57 PM ± 4.10	0.42, <0.001†	0.05†	0.78
Mean arterial pressure(mm Hg)	7.68 ± 0.63	11:50 PM ± 4.25	0.38, <0.001†	0.004†	0.60
IOP (mm Hg)	3.81 ± 0.29	12:13 PM ± 2.98	0.69, <0.001†	<0.001†	0.74
MOPP (mm Hg)	6.77 ± 0.54	11:59 PM ± 3.15	0.66, <0.001†	<0.001†	0.86

SEM, standard error of mean; CD, circular deviation; df, degree of freedom; SRL, superficial retinal layer; DRL, deep retinal layer.

P values from repeated measures ANOVA for time of day, and time of day by refractive error.

P < 0.05.

Relationships of the Superficial and Deep Retinal Layer OCT-A Indices Diurnal Rhythm

The diurnal amplitude of the superficial retinal layer perfusion density was negatively correlated with the parafoveal retinal thickness (rρ = −0.30, P = 0.046) and MOPP diurnal amplitudes (rρ = −0.44, P = 0.003). The diurnal amplitudes of superficial retinal layer foveal vessel density and perfusion density were weakly correlated with the diurnal amplitudes of pulse pressure (rρ = −0.39, P = 0.0093 and rρ = −0.36, P = 0.02) and IOP (rρ = −0.37, P = 0.02 and rρ = −0.39, P = 0.010). The deep retinal layer perfusion density diurnal amplitude showed significant correlation with parafoveal deep retinal layer thickness and total retinal thickness amplitudes (rρ = 0.40, P = 0.007 and rρ = −0.45, P = 0.002). The deep retinal layer parafoveal vessel density and perfusion density and axial length amplitudes were also correlated (rρ = 0.33, P = 0.03 and rρ = 0.49, P = 0.001).

The acrophases of superficial retinal layer perfusion density and thickness at parafovea were also negatively correlated (rcc = −0.34, P = 0.01) (Fig. 5). The acrophases of superficial retinal layer perfusion and vessel density at parafovea were significantly correlated with the acrophases of IOP (rcc = 0.40, P = 0.001 and rcc = 0.42, P < 0.001) and pulse pressure (rcc = 0.31, P = 0.04 for both) (Fig. 6). The IOP and FAZ metrics of deep retinal layer acrophases were also correlated (rcc = 0.45, P < 0.001 and rcc = 0.42, P < 0.001).

Figure 5.

Mean (±1 standard error of the mean) change in perfusion density (black) and thickness (gray) of superficial retinal layer (SRL) in the parafoveal region across 24 hours demonstrating their variations in approximately antiphase with each other and their acrophases were significantly correlated (rcc = −0.34, P = 0.01).

Figure 6.

Mean (±1 standard error of the mean) change in perfusion density (black) of superficial retinal layer (SRL) in the parafoveal region and pulse pressure (gray) across 24 hours demonstrating their variations in phase with each other and their acrophases were significantly correlated (rcc = 0.31, P = 0.04).

Discussion

This is the first study to evaluate the pattern and magnitude of diurnal variation in the retinal microvasculature over a full 24-hour period using OCT-A with comparisons between myopes and nonmyopes. Perfusion and vessel density demonstrated significant diurnal variations in the foveal, parafoveal, and perifoveal zones of the superficial and deep retinal layers, with their peaks between afternoon and late evening and troughs during the night and early morning hours except for perifoveal deep retinal layer vessel and perfusion density, which did not show significant change. The diurnal pattern or amplitude of the OCT-A indexes were also found to be correlated with changes in retinal thickness, axial length, IOP and pulse pressure. Furthermore, the findings showed significant differences in the timing of the diurnal rhythms of superficial retinal layer indexes at the parafovea between myopic and nonmyopic healthy adults but the diurnal rhythms in the foveal and perifoveal zones or the deep retinal layer indexes were not different between refractive groups.

A change in retinal blood flow over 24 hours could potentially cause a corresponding change in the detectability of fine capillaries by the OCT-A imaging method. This will eventually result in variation in the OCT-A indexes. Studies in the past have evaluated OCT-A variations using different protocols with measurements captured between 7 AM and 8 PM in the superficial and deep retinal layers among healthy adults (Penteado R, et al. IOVS 2018;59:ARVO E-Abstract 2856; Rommel F, et al. IOVS 2018;59:ARVO E-Abstract 2860).– We hypothesized that ocular blood flow may be different at night compared to day, so the present study recruited young healthy participants and captured measurements over 24 hours encompassing seven measurement points. This allowed a more comprehensive understanding of the normal diurnal variation in the macular microvasculature.

The diurnal variation of the superficial retinal layer perfusion and vessel density in the foveal and parafoveal regions showed similar trends, with measurements highest in late afternoon and lowest during sleeping hours. Most previous studies found insignificant daily fluctuations in superficial retinal layer OCT-A indexes between 7 AM to 8 PM (Rommel F, et al. IOVS 2018;59:ARVO E-Abstract 2860).– In the present study, most changes occurred during night between 9 PM to 5 AM with very small changes taking place between 9 AM and 9 PM. Only one previous study found a small but significant increase in the parafoveal region perfusion density between 12 PM and 8 PM, however the layer was not specified (Penteado R, et al. IOVS 2018;59:ARVO E-Abstract 2856). A small number of studies evaluated variations in deep retinal layer OCT-A indexes between 8 AM to 8 PM among healthy adults and reported them to be stable throughout the day.,,, Blood flow in the deep retinal layer might be expected to increase at night, since it is nearest to the photoreceptors layer and supplies for its increased metabolic demands in the dark. Nesper et al. found a significant transient increase in superficial retinal layer and decrease in deep retinal layer perfusion density during the transition from dark to light. However, the current study observed lowest deep indexes in the foveal region during the sleeping hours (2 AM and 3 AM) and in the parafoveal region during late morning (9 AM and 11:30 AM). The deep retinal layer is supplied by the vertical anastomoses from the superficial vascular plexus whereas the superficial retinal layer is directly supplied by the central retina artery. The capillary flow between the superficial and deep network happens through the anatomical capillary interconnections. This interconnection between these two capillary networks could explain the different physiological variations taking place in these layers.,,

In the current study, the diurnal amplitudes of superficial and deep retinal layer perfusion density and of total retinal thickness at the parafovea were negatively correlated, but not their acrophase. This is in contrast with the finding reported by Milani et al. where the parafoveal superficial retinal layer perfusion density and retinal thickness diurnal changes during the day showed a weak positive correlation. However, that study included healthy participants and those with glaucoma and ocular hypertension, and only two time points (morning and evening) were assessed. Though the reason is not clear, it is likely that changes in the outer retinal layers occurring over 24 hours may have contributed towards the relationship. In addition, the vessel and perfusion density variations (acrophase for superficial retinal layer: 5:29 PM and 4:50 PM; acrophase for deep retinal layer: 2:58 PM and 2:28 PM) were in-phase with total retinal thickness in the foveal region (acrophase: 3:44 PM). Retinal perfusion and thickness have also been more highly correlated in the foveal region than parafoveal region in the previous cross-sectional studies.,–

This study also, reports for the first time, significant diurnal variation in the parafoveal superficial and deep retinal layers thickness over 24 hours and evaluates their relationship with OCT-A indexes variations. The photoreceptor's metabolic demands are known to affect retinal vasculature under normal physiological conditions. At night, usually the photoreceptor outer segment + RPE layer is thickest whereas photoreceptor inner segment is thinnest when evaluated over 24 hours corresponding to the changes observed in superficial and deep retinal layer thickness respectively in the present study. We speculate that these similarities between the variations of different retinal layers might suggest a possible interaction between the physiological variations of retinal thickness and perfusion. The presence of diurnal changes in the retinal perfusion could also be related to the metabolic demands of the inner retinal cells or the presence of autonomous inner retinal clocks that regulate dopamine release and influence retinal perfusion.–

The superficial and deep retinal layer perfusion and vessel density variations (acrophase for superficial retinal layer: 6:13 PM and 5:53 PM; acrophase for deep retinal layer: 11:17 PM and 9:29 PM) were approximately antiphase to the superficial retinal layer thickness (acrophase: 4:05 AM) and deep retinal layer thickness (acrophase: 12:31 PM) variations in parafovea, respectively (Fig. 5). The mechanism behind this antiphase relationship is unclear, although we speculate that the thickness of the outer and inner segments of the photoreceptors at night might inversely affect the perfusion of the inner retinal layer's perfusion. Macular perfusion has been reported to significantly correlate with superficial retinal layer thickness.

Systemic assessment of the blood flow dynamics in the nocturnal and morning hours has been strongly emphasized in the literature owing to its association with abnormal cardiovascular events. It has been reported that retinal perfusion/blood flow is sensitive toward vascular factors such as MAP, pulse pressure and MOPP., Both MAP and pulse pressure showed significant diurnal variation with lowest readings around midday and early morning making pulse pressure approximately in-phase with superficial retinal layer perfusion and vessel density (Fig. 6). In addition, their diurnal amplitudes and acrophases were weakly correlated. This could be theoretically explained by high vascular resistance during morning, in turn increasing diastolic blood pressure and reducing pulse pressure and blood flow during the early hours., The current study also involved postural changes from supine to standing followed by sitting and no adaptation time was provided to adapt to these postural changes before the night OCT-A scans were taken. It has been recommended to wait about one to five minutes for measuring the maximum effect of posture change.– So, the effect of supine posture on the current study measurements of ocular blood flow at night cannot be completely ruled out.–

Axial length changes can result in subtle microstructural retinal changes especially in high myopia.,, Previous studies have shown positive correlation between deep retinal layer perfusion density and axial length. The diurnal amplitudes of the deep parafoveal vessel and perfusion density and axial length were positively correlated, and their variations were in-phase with each other (peaks between 2 PM and 3 PM) in the present study. In addition, myopes also showed significantly lower mean vessel density across the day in the parafoveal and perifoveal zones of both the layers than nonmyopes. Diurnal fluctuations in IOP and MOPP and their acrophases were also significantly correlated to that of superficial retinal layer vessel and perfusion density in this study considering that MOPP plays an integral role in regulating blood flow in the retina and is also affected by posture change. Baek et al. also found that the diurnal change in macular OCT-A indexes were associated with MOPP changes.

Myopes and nonmyopes showed significantly different patterns of variation over 24 hours for superficial retinal layer perfusion density of the parafovea with the acrophase occurring 4-hours earlier in myopes. This may arise due to the differences between myopes and nonmyopes in terms of ocular blood flow. For instance, a number of studies– showed increased vascular resistance index and decreased pulsatile ocular blood flow with increasing myopia or axial length. Considering that vascular resistance index is already high among myopes compared to nonmyopes, it further increases during the morning resulting in more myopes demonstrating their trough in OCT-A indexes in the morning. However, the present study did not find any difference in the variation pattern in axial length and systemic parameters between myopes and nonmyopes. It may be that local retinal vascular dynamics varied between the refractive groups. Though the retinal blood flow is autoregulated, Kur et al. reported local regulation of retinal blood flow in response to flickering light.

Though the diurnal changes observed are of relatively small magnitude, the current result provides normative data for OCT-A diurnal amplitude of vascular indexes over 24 hours. When expressed in terms of percentage of the daily mean, the diurnal amplitude ranged between 1.85% to 18.55% in different zones of the two layers. These results can serve as a baseline of expected change when analyzing the diurnal variation over 24 hours in conditions such as glaucoma and hypertension. The diurnal changes of OCT-A indexes were more marked in the parafoveal region and were significantly related to the ocular and systemic parameters, similar to the findings in patients with glaucoma. Though Bonferroni adjusted comparisons were conducted to reduce the chance of type 1 errors within each of the individual ANOVA comparing seven sessions, the use of multiple repeated measures ANOVA analysis for different Early Treatment Diabetic Retinopathy Study regions and layers might also increase the chances of type 1 errors. Replication of 24-hour diurnal variation studies in retinal microvasculature is warranted for more definitive conclusions. Our results might also be affected by the “first night effect” causing altered sleep due to participants sleeping in a new environment. It should be noted that our results are confined to young healthy adults and might not apply to children due to their ongoing refractive development or to older adults due to age-related vascular changes., Further studies are required to profile complete diurnal variation over 24 hours across age groups and ocular vascular diseases.

In conclusion, this study for the first time demonstrates that significant diurnal variation over 24 hours exists in OCT-A indexes of macular retinal microvasculature including superficial and deep retinal layer in healthy young adults. Myopes and nonmyopes showed variation only in the diurnal pattern of the parafoveal superficial retinal layer perfusion density but not in diurnal amplitude. The diurnal variation in the deep retinal layer indexes did not demonstrate any difference between the two refractive groups. The diurnal pattern and amplitude of the OCT-A indexes were found to be correlated with ocular and systemic measurements.