Subjects were 45 young adults, 33 (73.3%) of whom were female, with an average age (± SD) of 24.1 ± 1.8 years (range, 21.4–29.7 years). Refractive error was measured using the open-view Grand Seiko WR-5100K (Grand Seiko Co., Ltd., Hiroshima, Japan; distributed by AIT Industries, Bensenville, IL, USA). Ten readings were taken on each eye without cycloplegia but with the use of a Badal lens and track to relax accommodation and provide a clear target while subjects viewed the smaller lines on an acuity card at the far point of the eye. Grand Seiko autorefraction has been reported to be comparable to subjective refraction in adults.
49 Myopia was defined as having a spherical equivalent (SEQ) of –0.50 diopter (D) or more myopia in both eyes and an improvement in logMAR acuity by at least two lines with minus lens correction (
n = 28, 62.2%). All subjects classified as myopic habitually used corrections for their refractive error. Subjects less myopic/more hyperopic than –0.50 D SEQ in both eyes with distance acuity of 20/20 or better without correction were classified as non-myopic (
n = 16, 35.6%). One subject was unclassified because of having a low myopic, primarily spherical, refractive error but better than 20/20 distance acuity without correction. All subjects were included in the analysis of refractive error as a continuous variable, defined as the average SEQ of both eyes. No subject had significant anisometropia (≥1 D). Only five subjects had astigmatism ≥1 D, and none was a case of simple or mixed astigmatism. The average refractive error (± SD) was –1.73 ± 2.02 D (range, –6.33 to +1.70 D). Most of the sample was white (
n = 39, 87%); two subjects identifying as black, two as Hispanic, and two as Asian. The limited ethnic diversity did not allow for statistical analysis of this factor.
Because pupillary responses may be affected by prior activity,
47 environmental light exposure was measured as a covariate using the Daysimeter (Lighting Research Center, Rensselaer Polytechnic Institute, Troy, NY, USA). The device recorded average photopic lux values every minute along with a date and time stamp. It also recorded spectral data on three channels (blue, green, and red) in addition to lux. Spectral channel data were not analyzed due to the extremely high level of correlation between their output and lux. Subjects were instructed to wear the device mounted on an elastic strap on their upper arm outside their outermost layer of clothing during their normal waking hours for 7 days. Subjects were told to have the device facing outward from the torso while making no special attempt to have the device at any particular angle to the sun. When engaged in water sports such as boating or swimming, subjects were to have the device in their immediate environment but not in the water. The total lux-minutes of exposure were found for several time intervals prior to pupil testing: 1, 3, 12, 24, 72, and 120 hours. The shorter intervals included various amounts of night, morning, noon, or afternoon time, as testing was done at the convenience of the subject at any time of day during business hours. Season and time of day of testing were also recorded, as these have been shown to affect the melanopsin-mediated pupillary response in adults.
50–52 Birth month was noted because of a possible association with myopia prevalence.
53 Subjects were also asked to estimate the percent of time outdoors that they wore sunglasses. Exposure to environmental light was quantified as log
10 lux-minute values to create normal distributions for all intervals (Kolmogorov–Smirnov test for normality on the log-transformed data,
P < 0.20) except for 24 hours (excessive skew and kurtosis, Kolmogorov–Smirnov test for normality,
P = 0.015). Time outdoors was counted as any minute where illuminance exceeded 1000 lux.
54,55 A log transform of time outdoors was also used to create normal distributions.
Pupil testing was done using the RAPDx pupillometer (Konan Medical USA, Irvine, CA, USA). Pupillometry began within 10 minutes of arrival for testing after subjects’ ordinary daily activity. Subjects were positioned in front of the pupillometer without refractive correction after 5 minutes of dark adaptation (0.01 lux) before each of three separate trials. The length of time chosen for dark adaptation was somewhat arbitrary but was within the range of the 2 to 10 minutes used in other studies.
41,45,46 The light stimuli were presented in a 26.5° circular field to both eyes for 5 seconds interleaved with 5 seconds of dark (0.1 Hz). The light stimuli in the three trials were (1) pulses alternating between red and blue, (2) red only, and (3) blue only. The alternating presentation lasted for 2 minutes (six presentations of red interleaved with six of blue) and each of the single-color conditions lasted for 1 minute (six presentations of either red or blue). The order of presentation was the same for all subjects: the alternating red/blue stimulus was applied first, followed by the red-only stimulus and finally the blue-only stimulus (
Fig. 1A). The dark periods between the light stimuli had an overall radiance of 2.65 × 10
−3 W/sr/m
2, with a peak at 448 nm, due to minimal residual lighting from the liquid-crystal display screen within the instrument. At this peak wavelength, the irradiance was 2.9 log units lower than it was during exposure to the blue light stimulus. The spectral distributions and the irradiances of the red and blue stimuli used in these experiments were measured using a spectroradiometer (PR-670; Photo Research, Inc., Chatsworth, CA, USA) (
Fig. 1B). The peak intensity for the long-wavelength stimulus (red, for simplicity) was at 608 nm, and the peak for the short-wavelength stimulus (blue) was at 448 nm. This peak differs from the 480-nm maximum sensitivity of ipRGCs, but only on the order of 0.1 or less in log relative sensitivity.
56
Blink artifacts were either removed from the data by the RAPDx or identified as spikes and removed manually. The RAPDx records at 40 Hz but subject data were binned into averages at 0.25-second intervals for the purposes of simplifying data handling. Using the equation below, the diameter of the pupil was normalized at each 0.25-second interval separately for each eye relative to the smallest pupil diameter (100% constriction = 1.0) and largest pupil diameter (0% constriction = 0) that occurred for that eye across all three stimulus conditions, as in previous work.
48 Normalized pupillary responses were then averaged for the two eyes at each 0.25-second interval. Pupillary responses are reported as normalized pupillary constriction and are dimensionless unless otherwise marked in millimeters (mm).
\begin{eqnarray*}
&&Normalized\,pupillary\,constriction\\
&& = \frac{{Maximum\,diameter - Pupil\,diameter}}{{Maximum\,diameter - Minimum\,diameter}}\end{eqnarray*}
Four features of pupillary responses were selected for analysis. The ∆Blue variable was the average difference at each corresponding time point throughout testing between normalized pupillary responses during the presentation of blue in the single-color condition and those in the alternating-color condition. These differences were calculated across all time points, both when the stimulus was on during constriction and when it was off during redilation. This comparison between alternating and single-color presentation was chosen in order to highlight the adaptive change in pupil size between the earlier alternating and the later single-color condition. The ∆Red variable was the average difference at each corresponding time point throughout the testing between normalized pupillary responses during the presentation of red in the single-color condition and those in the alternating-color condition. Post-illumination pupillary redilation rates were calculated as the average of the six coefficients β for the exponential decay function eβt during the final 3 seconds (t) of each of the six pupillary redilations following presentation of blue only (ExpBlue). Only the data from the final 3 seconds were used because of the rapid rate of change during the first 2 seconds of redilation. The decay coefficient for redilation following presentation of red only was calculated in a similar manner (ExpRed).
Statistical analysis was performed using SPSS Statistics 21 (IBM, Armonk, NY, USA). Non-myopic and myopic subject characteristics were compared using independent t-tests or Fisher's exact test. The t-test P values were not adjusted for multiple testing. The six time periods for light exposure and time outdoors were compared between myopic and non-myopic subjects in separate repeated-measures ANOVA. The four pupillary response outcomes (∆Blue, ∆Red, ExpBlue, and ExpRed) were also analyzed using repeated-measures ANOVA with color (red or blue) and outcome type (Δ or Exp) as repeated factors. Myopic or non-myopic was a between-subject factor. Bivariate correlations were examined among the pupillary outcomes, light exposure, and SEQ using SPSS. Significant linear relationships were fit using the orthogonal regression procedure in JMP 10 (SAS Institute, Cary, NC, USA). General linear models were then used to examine multivariate associations among SEQ, environmental light exposure, age, and sex, including all two-way interactions. P < 0.05 was considered significant.