Evidence That Pupil Size and Reactivity Are Determined More by Your Parents Than by Your Environment

Abstract

Purpose: A classic twin study to evaluate the relative contributions of genetic and environmental factors to resting pupil size and reactivity. Methods: Pupillometry was performed on 326 female twins (mean age 64 years) from the TwinsUK Adult Twin Registry, assessing resting pupil diameter in darkness and increasing levels of ambient light, alongside dynamic pupillary characteristics. Maximum-likelihood structural equation models estimated the proportion of trait variance attributable to genetic factors. Results: Mean (SD) pupil diameter in darkness was 5.29 mm (0.81), decreasing to 3.24 mm (0.57) in bright light. Pupil light reaction (PLR) had a mean (SD) amplitude of 1.38 mm (0.27) and latency of 250.34 milliseconds (28.58). Pupil size and PLR were not associated with iris colour, intraocular pressure or refractive error, but were associated with age (diameter β = -0.02, p = 0.016, constriction amplitude β = -0.01, p < 0.001, velocity β = 0.03, p < 0.001, and latency β = 0.98, p < 0.001). In darkness the resting pupil size showed a MZ intraclass correlation coefficient of 0.85, almost double that of DZ (0.44), suggesting strong additive genetic effects, with the most parsimonious model estimating a heritability of 86% [95% confidence interval (CI) 79-90%] with 14% (95% CI 10-21%) explained by unique environmental factors. PLR amplitude, latency and constriction velocity had estimated heritabilities of 69% (95% CI 54-79%), 40% (95% CI 21-56%), and 64% (95% CI 48-75%), respectively. Conclusion: Genetic effects are key determinants of resting pupil size and reactivity. Future studies to identify these genetic factors could improve our understanding of variation in pupil size and pupillary reactions in health and disease.

Keywords: genetics; heritability; pupil; pupil size; twin study.

Figure 1

Histogram of left eye resting pupil diameter in darkness.

Figure 2

Measurements in twin pairs (twin 2 plotted against twin 1) for resting state 1 (Top), PLR latency, amplitude and constriction velocity (Bottom). Monozygotic (MZ) twin pairs are on the left and Dizygotic (DZ) on the right. Correlation coefficient (r) is given below for each parameter investigated.