- 1Section of Academic Ophthalmology, School of Life Course Sciences, Faculty of Life Sciences & Medicine, King's College London, London, United Kingdom
- 2Department of Ophthalmology and Epidemiology, University Medical Centre Groningen, University of Groningen, Groningen, Netherlands
- 3Department of Neuro-Ophthalmology, National Hospital for Neurology and Neurosurgery, London, United Kingdom
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.
Introduction
Pupil size is important not only because it controls the amount of light that is able to reach the retina but also because it affects the levels of chromatic and spherical aberration, thereby optimising visual perception (1). Pupil size is governed by the opposing actions of the iris sphincter and dilator muscles under the control of the parasympathetic and sympathetic nerves, respectively. In the healthy population there is significant variation in the resting size of the pupil and in the reflex constriction of the pupil to light. In any individual person pupil size is influenced by a number of factors including ambient light, retinal sensitivity, central cognitive processes (2, 3), alertness and emotional unrest (4, 5) with pupil measurements additionally being used by researchers as a biomarker of activity in the autonomic nervous system- both in health and in disease. Many studies have evaluated this variation in pupil size in different illumination levels, suggesting pupil size decreases in a linear manner for each luminance level (6). At a population level it has also been shown that pupils tend to be smaller with increasing age (7). This inverse correlation between age and resting pupil size has been confirmed pre-operatively in refractive surgery candidates across a range of different luminance levels (8, 9), but the effect seems to be most marked at low luminance suggesting that it is due to a progressive decline in sympathetic tone throughout adult life (7, 10).
Certain medical conditions such as diabetes and pseudoexfoliation can also decrease the effect of dilating drops on the size of the pupil (11). Different classes of glaucoma medication also influence pupil size (1). Although it was previously believed that pupil size may decrease significantly after cataract surgery, more recent literature suggests the reduction in size is temporary and returns to preoperative levels by 1 month (12).
In addition to this variability in resting pupil size, clinically there is wide variation in the dynamic pupillary light reflex (PLR) in healthy individuals. A more rigorous assessment of the PLR is now possible using pupillometric measurements using infrared video techniques. Measurements obtained include the amplitude of the light reflex, latency (13), velocity (peak or average) (14) and acceleration (15), all of which show wide variation in the healthy population. In part this may be due to the variation in resting pupil size, since a smaller pupil admits less of the stimulus under closed-loop conditions; an autonomic control system whereby a mechanism is regulated by feedback. However, even in studies where the light stimulus has been focussed to a size narrower than even the smallest pupil (i.e., “Maxwellian” optics which create open-loop conditions) and amplitude measurements are expressed as percentage change relative to the starting diameter, there is still a varied range in all PLR measurements across healthy individuals (16).
Over recent years, automated pupillometry has evolved into a robust and reliable measure of pupillary evaluation. This has removed the subjectivity of pupillary evaluation, providing data which is more accurate and clinically significant. Its use has been supported in the monitoring neurointensive care patients, for whom subtle changes can be a predictive factor of neurological deterioration (17). More recently a number of studies have highlighted its clinical relevance including in the evaluation of refractive surgery (18), sleep studies (19), monitoring the effects of certain medication on the central nervous system and even as a measure of emotional response in patients with psychiatric disorders (20). However, despite the significant clinical implication of understanding determinants of resting pupil size and reactivity to light, little is known about the role of genes. An improved understanding could enable clinicians to more reliably use pupil measurements to detect and monitor a wide range of diseases including glaucoma, other forms of optic neuropathy, and idiopathic intracranial hypertension (21). Often these diseases are bilateral and thus the fellow eye cannot be used as an internal control; distinguishing disease from health in these cases can be extremely challenging, instead pupil measurements must be compared with the normative database. However, the “normal ranges” for these pupil measurements are discouragingly wide, lowering the sensitivity for disease detection and the usefulness of pupil testing in clinical practise. The potentially high genetic determinate of pupil measurements will provide an improved belief on the stability of these results, their clinical implications and reduce the fear of environmental influence. Additionally, this improved understanding of pupil size and constriction may allow further appreciation of the mechanisms involved. This can assist in identifying individuals whose pupillary characteristics are different from those predicted by genetic risk scores thereby enabling earlier identification of abnormalities.
Twin studies provide us with the unique ability to examine the relative contribution of environmental and genetic factors to trait variation. Heritability is the quantification of the overall phenotypic variation that is attributable to genetic factors. Heritability estimates allow us to define boundaries for the ability of genetics to predict traits or disease and indirectly provides insights into the role of genes (22). Heritability calculations are performed by comparing the concordance of disease or covariance of quantitative traits in non-identical dizygotic (DZ) and identical monozygotic (MZ) twin pairs. Given that MZ twins share nearly 100% of their genetic data and DZ only half, any greater similarity within MZ twin pairs can be ascribed to this additional genetic sharing (11).
The purpose of this study was to undertake a classical twin study evaluating the relative contributions of genetic and environmental factors to resting pupil size under varying degrees of background illuminance and the PLR (the dynamic pupil response to a transient light stimulus).
Materials and Methods
Participants
Healthy twin volunteers were recruited from the TwinsUK registry based at St Thomas' Hospital, London (23). Twin participants attending for ocular examinations as part of ongoing studies were additionally recruited to undergo pupillometry between January 2010 and April 2011. Pupillometry was used to assess the resting pupil diameter in darkness and in increasing levels of ambient light, and also the dynamic pupillary response to a standardised transient light stimulus. Additional information was collected on patient demographics (age at examination, sex), iris colour (brown, blue or green), non-cycloplegic autorefraction (ARM-10 autorefractor, Takagi Seiko, Japan), intra-ocular pressure (non-contact tonometry using Visionix PT100), past ocular history, past medical history (including neurological disorders and diabetes) and a complete list of current medications including any topical eye drops. Twin zygosity was determined either by genotyping or through the “peas in the pod” questionnaire (PPQ). The PPQ questionnaire has previously been validated for its excellent accuracy as a proxy indicator of zygosity in the absence of genotyping information (24). Spherical equivalent was calculated in the standard manner for each eye as sphere + (cylinder/2).
Inclusion criteria included a minimum age of 18, the capacity to understand the study and consent to participate. Patients were excluded if they were a known diabetic given the well-accepted association with autonomic pupillary dysfunction. The study had appropriate local research ethics approval (EC04/015) and was conducted in accordance with the tenets of the Declaration of Helsinki.
Pupillometry and Protocols
All pupillometry measurements were taken using the Procyon P3000 pupillometer (Procyon Instruments Ltd.). This is a binocular dynamic pupilometer that allows for customised stimulus protocols and uses video cameras running at 25 Hz with a spatial resolution of 0.03 mm. This device and similar protocols have previously been described elsewhere (25). In brief, each study participant had the resting pupil diameter measured in increasing levels of ambient light: “dark,” “low,” “medium,” and “high” (labelled as resting states 1–4, respectively); the subject was equilibrated to each ambient light level for 30 s before measuring the resting pupil size (averaged over 3 s). The background illuminance levels used were the four standard light levels provided by the manufacturer and “factory calibrated” such that the illumination at the cornea was ~0.00 Lux (“dark”), 0.04 Lux (“low”), 0.40 Lux (“medium”), and 4.00 Lux (“high”). The “standard” light stimulus used to elicit the PLR = 4.00 Lux. A series of standard white light stimuli (duration 1.0 s, interstimulus interval 5 s) was then presented alternately to right and left eyes (total of three stimulus presentations to each eye) and the pupil responses averaged to measure various outcome parameters associated with the direct and consensual pupil light reflex (amplitude, latency, and maximum constriction velocity).
Analysis
Details of twin modelling methods have been previously described (26, 27). The methodology is based on the evaluation of the variance and covariance matrices between (MZ) and (DZ) twin pairs.
The approach utilises structural equation modelling to separate the observed phenotypic variance into additive (A) genetic effects, dominant genetic effects (D) or common environmental (C) effects, and unique environmental (E) effects, the latter also including measurement error. By subsequently dividing these individual components by the total variance you are able to determine the different standardised components of variance. To adjust for the effects of age on studied pupil traits, we performed linear regression models with pupil traits and age at visit, from which we obtained the residuals. Heritability models were then performed using these residuals. Structural equation modelling is used to estimate best-fit model that explains the variance and covariance of a trait, using the phenotype information from each twin pair and knowing their genetic relationship. Only three variables can be inferred to fit the data from two twins, so an ACDE model cannot be used; ACE or ADE models have to be used separately. An ADE model implies non-additive genetic effects, which may the case when the MZ correlation is more than twice the DZ correlation- otherwise an ACE model is chosen, implying additive genetic effects explain the greater MZ twin correlation. Univariant models were initially created using all specified parameters (ACE/ADE) and then using a stepwise approach, individual parameters were dropped, testing the deterioration in model fit accordingly. This creates three parsimonious models with increasing degrees of freedom (AE, CE, and E). Twice the difference in log likelihoods between the full and sub models follows a x2 distribution, with degrees of freedom equal to the difference in degrees of freedom between the full and subsequent models (likelihood ratio test). Individual sub models are subsequently compared with the full model with the best fitting model determined by identification of the model with the lowest Akaike information criterion (AIC). The AIC is a measure of the relative quality of statistical models for a given data set. Heritability is consequently estimated from the best-fitting model from the total contribution of genetic factors to trait variance.
Statistical analysis was performed using STATA (Stata Corp., College Station, TX). Twin modelling was completed using the OpenMx (http://openmx.psyc.virginia.edu) package in R (http://www.R-project.org), with results from the left eye reported.
Results
Three hundred eighty-five twins were considered for our study. Participants were excluded if they were unpaired or if zygosity was unknown (n = 7), diabetic (n = 30) or if incomplete measurements were obtained for the parameters investigated (n = 22). In total results from 326 twins were included [70 MZ (43%) and 93 DZ (57%) twin pairs] in heritability calculations. Zygosity was determined by genotyping in 94% of the population (n = 305) and through PPQ in 6% (n = 21). Mean age at time of examination was 63.7 years [Standard Deviation (±) 7.9 (range 37–81)] and all participants were female (100%). Only one twin pair was not of European ethnicity (99% European). Self-reported iris colour was available in 232 twins, of which 144 were blue (61.2%), 50 brown (22.4%) and 38 green (16.5%). Mean spherical equivalent was −0.25 D (±2.9) in the right eye and −0.25 D (±2.8) in the left eye. Mean IOP of both eyes was 15.1 mmHg (±3.3).
Regarding medical history, 1 (0.3%) had a neurological disorder and 10 (3.1%) were pseudophakic. Thirty-two patients took regular eye drops [indications included glaucoma (n = 5), dry eye disease (n = 25) and other (n = 2)]. Eighteen (5.5%) patients had other forms of eye surgery including squint surgery (n = 5), oculoplastic procedures (n = 6), retinal detachment surgery (n = 4) and undefined laser procedures (n = 3). There was no significant difference in the means of any of the pupil measurements between MZ and DZ twins. Full details can be seen in Table 1.
The left eye resting pupil size varied widely, from <3 mm to over 7 mm in darkness, with an approximately normal distribution (Figure 1); similarly, there was wide variation in the PLR measurements. However, in a linear regression model incorporating age, sex, spherical equivalent, intraocular pressure, eye colour and whilst clustering by unique family identifier, only age was significantly associated with resting pupil diameter (β = −0.02, p = 0.016), PLR amplitude (β = −0.01, p ≤ 0.001), constriction velocity (β = 0.03, p ≤ 0.001) and PLR amplitude (β = 0.98, p ≤ 0.001) (Table 2).
Table 2. Multivariable linear regression models examining the association of age, iris colour, intraocular pressure and spherical equivalent on pupil characteristics in the left eye, adjusted for family structure.
Within our cohort individuals with blue irides had the smallest resting pupil diameter (mean 5.21 mm, ±0.77) vs. 5.52 mm (±0.80) in green irides vs. 5.35 mm (±0.90) in brown irides, p trend = 0.468, no difference in mean age was seen between groups (p trend = 0.125). The PLR constriction velocity was slowest in blue irides (mean −3.96 mm/s, ±0.81) vs. −4.29 mm/s (±0.96) in green irides vs. −4.21 mm/s (±0.69) in brown irides, but again no statistical trend difference between iris colours was observed (p trend = 0.669) (Table 3).
Table 3. Mean left eye pupil measurements for resting pupil diameter, PLR amplitude, constriction velocity (CV) and pupil latency for brown, green and blue eyes.
Correlations between MZ-DZ pairs were calculated prior to the completion of twin modelling through calculating within-pair correlations by zygosity. Resting pupil size measurements demonstrated a stronger correlation within MZ twin pairs than DZ pairs at all four levels of illumination, and also at measurements of PLR latency and amplitude. The left eye resting state 1 (Dark) intra-pair correlation coefficient for MZ twins was 0.85, almost double that of DZ 0.44, suggesting a strong genetic effect on pupil size. Similar results were seen for dynamic parameters such as PLR amplitude, constriction velocity and pupil latency. Scatter plots depicting intra-pair correlation for MZ and DZ twins for the studied pupil parameters are shown in Figure 2.
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.
Estimation of Heritability
For the sample of complete twin pairs (n = 163) model fit statistics suggested the AE model provided the most parsimonious fit for all traits studied (Table 4). This would suggest the variance in resting pupil size, amplitude, constriction velocity and pupil latency are best explained by additive genetic and unique environment effects, whilst common family environment and dominant genetic effects could be dropped with no significant loss of fit.
Standardised parameter estimates for the best fitting models for the left eye can be seen in Table 5. This would suggest the left eye heritability of the resting pupil size in darkness (state 1) is 86% (95% CI 79–90%) with the remaining 14% (95% CI 10–21%) explained by unique environment. Similar results were seen for other resting states examined, with heritability estimated to be 76% (95% CI 65–83%) for resting state 2, 76% (95% CI 65–83%) for resting state 3 and 70% (95% CI 57–79%) for resting state 4, respectively.
Other dynamic variables investigated included PLR amplitude, constriction velocity and latency; the left eye heritability for these parameters was estimated to be between 69% (95% CI 54–79%), 64% (95% CI 48–75%), and 40% (95% CI 21–56%), respectively.
Discussion
We have evaluated the relative contributions of genetic and environmental factors to resting pupil size under varying degrees of illuminance, and also pupil reactivity to transient light stimuli. Our results indicate that resting pupil size in complete darkness is strongly heritable with additive genes explaining up to 86% of the variance. As ambient light luminance was increased the heritability estimates fell, as expected given that many ophthalmic and external factors may affect pupillary constriction. A comparably lower heritability estimate was seen for dynamic characteristics such as PLR amplitude, velocity and latency, with variation due to additive genes ranging between 40 and 69%, respectively. This suggests that environmental factors are only responsible for 14% of the variability in resting pupil size within the population, and for between 31 and 60% of the variability in the PLR. It is noteworthy that baseline pupil diameter and reactivity have an equivalent heritability to many other highly heritable ocular traits such as central corneal thickness (95%) (28), optic disc parameters (66–73%) (29) and refractive error (77%) (30). We additionally examined for any association between pupil characteristics and iris colour, intra-ocular pressure and spherical equivalent. It has been reported that blue pupils may have a larger resting pupil diameter and be more dynamic (31, 32) however recent studies suggest this may not be the case (33, 34). Although our sample size was limited, we found no significant association between eye colours. We did however confirm a significant association with age for all characteristics (pupil resting diameter β = −0.02, p = 0.016), pupil light reaction amplitude (β = −0.01, p ≤ 0.001), constriction velocity (β = 0.03, p ≤ 0.001) and pupil latency (β = 0.98, p ≤ 0.001).
To date only two studies have attempted to evaluate the heritability of pupil size. The Guangzhou Twin Eye Study investigated the distribution and heritability of iris thickness and pupil diameter using anterior segment OCT (ASOCT). Their population of predominately East Asian individuals were aged between 8 and 16 years. Pupil diameter was measured between the most central points of the iris border. They estimated that genetic influences account for 63% of pupil size in a dark room (<5 lux) (35). This heritability estimate appears comparable to our findings. Our findings of high heritability are also comparable to that identified for pupil size after pharmacological mydriasis, which has been reported to be between 78 and 80% (intraclass correlation coefficients 0.82 in MZ twins vs. 0.39 in DZ twins) (11). To our knowledge no previous study has quantified the heritability of dynamic PLR measurements.
Although a number of factors that impact the pupil size and its dynamics have been well-described (36–39), there may be other influences including mechanical/anatomical factors. For example, iris thickness has been postulated as influencing the amplitude of pupillary constriction at high intensities using photic stimulation with blue and red light, possibly due to a passive resistance (39). Such local anatomy is highly heritable and would support the need for further genetic evaluation (35). Although limited, family-based genetic linkage studies have previously identified a congenital microcoria locus (40, 41), associated with maldevelopment of the dilator pupillae muscle of the iris and often associated with juvenile-onset glaucoma. Rare genetic mutations causing abnormal development of the anterior segment and iris are well-known, with a number of hereditary disorders being easily recognised through characteristic pupillary changes, including aniridia, iris colobomas, ectopic, scalloped or peninsula pupils, polycoria, congenital mydriasis, persistent pupillary membrane remnants and congenital miosis (42). The role of non-coding polymorphisms in these genes influencing normal pupil size variation and reactivity is unknown.
Additionally, given the potential environmental influence on PLR, a number of already established factors must be appreciated. In healthy individuals this may include: (a) retinal sensitivity to light (43); (b) any extrinsic factor, including those of dietary or pharmacological aetiology which may affect the function of retinal ganglion cells, in particular, mitochondrial function (44); (c) central influences modulating the excitability of neurones in the olivary pretectal nuclei and the Edinger-Westphal nuclei (45), (d) factors that may influence transmission across the neuro-effector junction in the iris sphincter muscle such as exposure to drugs or chemicals that modulate the availability and function of muscarinic receptors and (e) anything that affects the iris and its mobility, including the size of the lens, iris perfusion, oxygenation and any previous anterior segment pathology, anterior segment surgery or laser treatment (46).
Heritability is a population-specific factor and although our twin population data has been shown to be generalizable to the singleton population [as twins have similar morbidity and mortality to the rest of the population (47)], it consists predominantly of British women of European ancestry of a relatively narrow age range of whom are predominately in their sixties. Other limitations include that of volunteer or recruitment bias which could lead to an overestimation of the heritability of a trait. Additionally, the simultaneous effects of both shared-environment and gene/environment interactions cannot be considered concurrently without including additional siblings in the design.
The lower heritability of dynamic PLR characteristics in comparison to resting pupil size may be in part limited by measurement error, environmental factors or unknown shared family factors. Pupil latency, for instance, was seen to have relatively low heritability, however, we believe this is reflective of inaccuracy in its measurement. The pupilometer used possesses a more accurate spatial resolution in comparison to its frequency resolution. Potentially explaining any imprecise results. Given these limitations, we believe our results warrant further exploration on a larger scale, to determine if these conclusions are applicable to the wider population.
Future work should include genome-wide association studies to identify genetic polymorphisms associated with pupil diameter variance. A greater understanding of the genetic factors and biological mechanisms underlying pupil size and reactivity might allow greater clinical utility of assessing normal pupil size and reactivity in health and disease. Assessing pupillary responses can also be a useful clinical tool in monitoring certain neurological conditions or identifying those that warrant further investigation. As an example, unilateral macular diseases rarely show a relative afferent pupillary defect (RAPD) of more than 0.5 log units, thus if a measured RAPD is greater than this, additional pathologies such as compressive optic neuropathy must be considered (21). Additionally, given the widespread disparity of normative measures in the general population, the clinical significance of a solitary pupillary measurement may hold limitations, but this could be improved by obtaining serial measurements from the same individual over a period of time, therefore offering a surrogate marker of disease progression.
In conclusion, this study provides a current, well-phenotyped normative dataset of pupillary size and reactivity in a largely European female population. Twin modelling demonstrates that genetic effects are a key determinant of both resting pupil size, particularly in dim light, and reactivity characteristics including latency, constriction velocity and amplitude. The knowledge that pupil size and reactivity is predominantly genetically determined should direct future research into understanding functional mechanisms influencing pupil size and ultimately improve our ability to use these measures clinically.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by King's College London - research ethics approval (EC04/015). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
FB and CH were involved in the conception and design of the study. FB and JV were involved in data collection. JV, CH, AA, FB, and KW carried out the statistical analysis and data interpretation, provided critical review of the manuscript, and contributed to the final write-up. AA and KW drafted the manuscript. All authors read and approved the final manuscript.
Funding
TwinsUK was funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. KW acknowledges funding from a Chadburn Academic Clinical Lectureship and AA from a NIHR Academic Clinical Fellowship. KW acknowledges funding from a Medical Research Council UK Clinical Research Training Fellowship (MRK0237211).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank all the participants and teams involved in the completion of this study.
References
1. Park JW, Kang BH, Kwon JW, Cho KJ. Analysis of various factors affecting pupil size in patients with glaucoma. BMC Ophthalmol. (2017) 17:1–8. doi: 10.1186/s12886-017-0564-6
2. Kucewicz MT, Dolezal J, Kremen V, Berry BM, Miller LR, Magee AL, et al. Pupil size reflects successful encoding and recall of memory in humans. Sci Rep. (2018) 8:4949. doi: 10.1038/s41598-018-23197-6
3. Aminihajibashi S, Hagen T, Foldal MD, Laeng B, Espeseth T. Individual differences in resting-state pupil size: evidence for association between working memory capacity and pupil size variability. Int J Psychophysiol. (2019) 140:1–7. doi: 10.1016/j.ijpsycho.2019.03.007
4. Nakamura K, Bissen-Miyajima H, Oki S, Onuma K. Pupil sizes in different Japanese age groups and the implications for intraocular lens choice. J Cataract Refract Surg. (2009) 35:134–8. doi: 10.1016/j.jcrs.2008.10.020
5. Kobashi H, Kamiya K, Ishikawa H, Goseki T, Shimizu K. Daytime variations in pupil size under photopic conditions. Optometry Vis Sci. (2012) 89:197–202. doi: 10.1097/OPX.0b013e31824048a9
6. Winn B, Whitaker D, Elliott DB, Phillips NJ. Factors affecting light-adapted pupil size in normal human subjects. Invest Ophthalmol Vis Sci. (1994) 35:1132–7.
7. Guillon M, Dumbleton K, Theodoratos P, Gobbe M, Wooley CB, Moody K. The effects of age, refractive status, and luminance on pupil size. Optom Vis Sci. (2016) 93:1093–100. doi: 10.1097/OPX.0000000000000893
8. Netto MV, Ambrósio R, Wilson SE. Pupil size in refractive surgery candidates. J Refract Surg. (2004) 20:337–42. doi: 10.3928/1081-597X-20040701-06
9. Linke SJ, Baviera J, Munzer G, Fricke OH, Richard G, Katz T. Mesopic pupil size in a refractive surgery population (13,959 eyes). Optom Vis Sci. (2012) 89:1156–64. doi: 10.1097/OPX.0b013e318263c165
10. Mathur A, Gehrmann J, Atchison DA. Influences of luminance and accommodation stimuli on pupil size and pupil center location. Investig Ophthalmol Visual Sci. (2014) 55:2166–72. doi: 10.1167/iovs.13-13492
11. Hammond CJ, Snieder H, Spector TD, Gilbert CE. Factors affecting pupil size after dilatation: the Twin Eye Study. Br J Ophthalmol. (2000) 84:1173–6. doi: 10.1136/bjo.84.10.1173
12. Hayashi K, Hayashi H. Pupil size before and after phacoemulsification in nondiabetic and diabetic patients. J Cataract Refract Surg. (2004) 30:2543–50. doi: 10.1016/j.jcrs.2004.04.045
13. Alexandridis E, Argyropoulos T, Krastel H. The latent period of the pupil light reflex in lesions of the optic nerve. Ophthalmologica. (1981) 182:211–7. doi: 10.1159/000309116
14. Patwari PP, Stewart TM, Rand CM, Carroll MS, Kuntz NL, Kenny AS, et al. Pupillometry in congenital central hypoventilation syndrome (CCHS): quantitative evidence of autonomic nervous system dysregulation. Pediatr Res. (2012) 71:280–5. doi: 10.1038/pr.2011.38
15. Giza E, Fotiou D, Bostantjopoulou S, Katsarou Z, Gerasimou G, Gotzamani-Psarrakou A, et al. Pupillometry and 123I-DaTSCAN imaging in Parkinson's disease: a comparison study. Int J Neurosci. (2012) 122:26–34. doi: 10.3109/00207454.2011.619285
16. Bremner FD. Pupillometric evaluation of the dynamics of the pupillary response to a brief light stimulus in healthy subjects. Investig Ophthalmol Visual Sci. (2012) 53:7343–7. doi: 10.1167/iovs.12-10881
17. Martínez-Ricarte F, Castro A, Poca M, Sahuquillo J, Expósito L, Arribas M, et al. Infrared pupillometry. Basic principles and their application in the non-invasive monitoring of neurocritical patients. Neurología. (2013) 28:41–51. doi: 10.1016/j.nrl.2010.07.028
18. Tandogan T, Son HS, Choi CY, Knorz MC, Auffarth GU, Khoramnia R. Laboratory evaluation of the influence of decentration and pupil size on the optical performance of a monofocal, bifocal, and trifocal intraocular lens. J Refract Surg. (2017) 33:808–12. doi: 10.3928/1081597X-20171004-02
19. Ostrin LA, Abbott KS, Queener HM. Attenuation of short wavelengths alters sleep and the IP RGC pupil response. Ophthalmic Physiol Optics. (2017) 37:440–50. doi: 10.1111/opo.12385
20. Wang C-A, Baird T, Huang J, Coutinho JD, Brien DC, Munoz DP. Arousal effects on pupil size, heart rate, and skin conductance in an emotional face task. Front Neurol. (2018) 9:1029. doi: 10.3389/fneur.2018.01029
21. Bremner FD. Pupil assessment in optic nerve disorders. Eye. (2004) 18:1175–81. doi: 10.1038/sj.eye.6701560
22. Mayhew AJ, Meyre D. Assessing the heritability of complex traits in humans: methodological challenges and opportunities. Curr Genom. (2017) 18:332–40. doi: 10.2174/1389202918666170307161450
23. Verdi S, Abbasian G, Bowyer RCE, Lachance G, Yarand D, Christofidou P, et al. TwinsUK: the UK adult twin registry update. Twin Res Hum Genet. (2019) 22:523–9. doi: 10.1017/thg.2019.65
24. Jarrar ZA, Ward KJ, Mangino M, Cherkas LF, Gill R, Gillham-Nasenya I, et al. Definitive zygosity scores in the peas in the pod questionnaire is a sensitive and accurate assessment of the zygosity of adult twins. Twin Res Hum Genet. (2018) 21:146–54. doi: 10.1017/thg.2018.9
25. Bremner F. Apraclonidine is better than cocaine for detection of horner syndrome. Front Neurol. (2019) 10:55. doi: 10.3389/fneur.2019.00055
26. Neale M, Cardon LR. Methodology for Genetic Studies of Twins and Families. Richmond: Springer Science & Business Media (2013).
27. Mahroo OA, Williams KM, Hossain IT, Yonova-Doing E, Kozareva D, Yusuf A, et al. Do twins share the same dress code? Quantifying relative genetic and environmental contributions to subjective perceptions of “the dress” in a classical twin study. J Vis. (2017) 17:29–29. doi: 10.1167/17.1.29
28. Toh TY, Liew SHM, Mackinnon JR, Hewitt AW, Poulsen JL, Spector TD, et al. Central corneal thickness is highly heritable: the twin eye studies. Investig Ophthalmol Visual Sci. (2005) 46:3718–22. doi: 10.1167/iovs.04-1497
29. Healey P, Carbonaro F, Taylor B, Spector TD, Mitchell P, Hammond CJ. The heritability of optic disc parameters: a classic twin study. Investig Ophthalmol Visual Sci. (2008) 49:77–80. doi: 10.1167/iovs.07-0962
30. Lopes MC, Andrew T, Carbonaro F, Spector TD, Hammond CJ. Estimating heritability and shared environmental effects for refractive error in twin and family studies. Investig Ophthalmol Visual Sci. (2009) 50:126–31. doi: 10.1167/iovs.08-2385
31. Said FS, Sawires WS. Age dependence of changes in pupil diameter in the dark. Optica Acta Int J Optics. (1972) 19:359–61. doi: 10.1080/713818582
32. Vaswani RS, Mudgil AV, Gleicher D. Correlation of pupil size to iris color in children. J Refract Surg. (2002) 18:189–189. doi: 10.3928/1081-597X-20020301-15
33. Bergamin O, Schoetzau A, Sugimoto K, Zulauf M. The influence of iris color on the pupillary light reflex. Graefes Arch Clin Exp Ophthalmol. (1998) 236:567–70. doi: 10.1007/s004170050122
34. Bradley JC, Bentley KC, Mughal AI, Bodhireddy H, Young RS, Brown SM. The effect of gender and iris color on the dark-adapted pupil diameter. J Ocul Pharmacol Ther. (2010) 26:335–40. doi: 10.1089/jop.2010.0061
35. He M, Wang D, Console JW, Zhang J, Zheng Y, Huang W. Distribution and heritability of iris thickness and pupil size in Chinese: the Guangzhou twin eye study. Investig Ophthalmol Visual Sci. (2009) 50:1593–7. doi: 10.1167/iovs.08-2735
36. Lowenstein O, Loewenfeld I, Davson H. The Pupil, the Eye. Editcd by Davson H. New York, NY:Academic (1969). p. 311–4.
37. Hess E. Pupillometrics. In: Greenfield NS, Stembach RA, editors. Handbook of Psychophysiology. New York, NY: Holt, Rinehart and Winston (1972).
38. Wyatt HJ, Musselman JF. Pupillary light reflex in humans: evidence for an unbalanced pathway from nasal retina, and for signal cancellation in brainstem. Vis Res. (1981) 21:513–25. doi: 10.1016/0042-6989(81)90097-3
39. Sharma S, Baskaran M, Rukmini AV, Nongpiur ME, Htoon H, Cheng C-Y, et al. Factors influencing the pupillary light reflex in healthy individuals. Graefes Arch Clin Exp Ophthalmol. (2016) 254:1353–9. doi: 10.1007/s00417-016-3311-4
40. Rouillac C, Roche O, Marchant D, Bachner L, Kobetz A, Toulemont PJ, et al. Mapping of a congenital microcoria locus to 13q31-q32. Am J Hum Genet. (1998) 62:1117–22. doi: 10.1086/301841
41. Ramprasad VL, Sripriya S, Ronnie G, Nancarrow D, Saxena S, Hemamalini A, et al. Genetic homogeneity for inherited congenital microcoria loci in an Asian Indian pedigree. Mol Vis. (2005) 11:934–40.
42. Morrison PJ. The iris–a window into the genetics of common and rare eye diseases. Ulster Med J. (2010) 79:3.
43. Roecklein K, Wong P, Ernecoff N, Miller M, Donofry S, Kamarck M, et al. The post illumination pupil response is reduced in seasonal affective disorder. Psychiatry Res. (2013) 210:150–8. doi: 10.1016/j.psychres.2013.05.023
44. Bremner FD, Shallo–Hoffmann J, Riordan–Eva P, Smith SE. Comparing pupil function with visual function in patients with Leber's hereditary optic neuropathy. Investig Ophthalmol Visual Sci. (1999) 40:2528–34.
45. Larsen RS, Waters J. Neuromodulatory correlates of pupil dilation. Front Neural Circuits. (2018) 12:21. doi: 10.3389/fncir.2018.00021
46. Yilmaz I, Perente I, Saracoglu B, Yazici AT, Taskapili M. Changes in pupil size following panretinal retinal photocoagulation: conventional laser vs pattern scan laser (PASCAL). Eye. (2016) 30:1359–64. doi: 10.1038/eye.2016.135
Keywords: pupil size, heritability, genetics, pupil, twin study
Citation: Ansari AS, Vehof J, Hammond CJ, Bremner FD and Williams KM (2021) Evidence That Pupil Size and Reactivity Are Determined More by Your Parents Than by Your Environment. Front. Neurol. 12:651755. doi: 10.3389/fneur.2021.651755
Received: 10 January 2021; Accepted: 29 March 2021;
Published: 29 April 2021.
Edited by:
Valerie Purvin, Midwest Eye Institute, United StatesReviewed by:
Juan R. Ordonana, University of Murcia, SpainMárta Janáky, University of Szeged, Hungary
Copyright © 2021 Ansari, Vehof, Hammond, Bremner and Williams. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katie M. Williams, katie.williams@kcl.ac.uk