Goran Benčić, Veronique Vitart, Igor Rudan, Zoran Vatavuk, Harry Campbell, Ozren Polasek, Jelena Škunca Herman, Ivana Kolcic, Kajo Bućan, Alan F. Wright, Susan Campbell, Caroline Hayward, Jennifer E. Huffman, and Lina Zgaga
Studying quantitative endophenotypes was advocated to help unravel the genetic architecture of common diseases.1,2 Successes met by this approach include mapping of genes modulating QT elongation measured by ECG and cardiac arrhythmia risk,3 IgE levels and asthma risk,4 serum uric acid level and gout risk,5 and lipid levels and coronary heart disease risk.6 Ocular conditions, in particular the most common one, refractive error, lend themselves very well to this approach. Myopia and hypermetropia can be viewed largely as defects in the eye growth processes that normally adjust AL of the eye to the optical power of the cornea and lens. The values of the separate refractive components (axial length [AL], power of the cornea, and power of the lens), which if uncoordinated lead to refractive errors, have long been recognized as being normally distributed in general population surveys, whereas the distribution of refraction itself has a greater density around emmetropic values.7 Researchers in several large studies of unselected individuals, predominantly twins, have investigated to what extent genetic variation contributes to ocular quantitative components, and results have generally supported a substantial polygenic contribution. These include reports on AL, anterior chamber depth (ACD), corneal curvature (CC), and spherical equivalent refraction (SER) in a Sardinian isolate (n = 741; mean age, 41 years)8; in the Australian GEM twin study (n = 1224; mean age, 52 years)9; and in a Danish twin cohort (n = 114; age range, 20–45 years)10 together with lens thickness (LT), and analysis of refraction alone in a UK female twin cohort (n = 506; mean age, 62.4 years)11 and in the Beaver Dam population study (n = 2138; age range, 43–84 years).12 For corneal thickness (CT) there is, to our knowledge, only one previous report of heritability, 95% in a European sample of UK and Australian twins (n = 256; mean age, 38 years).13 This trait is now a recognized risk factor for progression from ocular hypertension to primary open-angle glaucoma,14 as well as a determinant of corneal refractive power. There has been a call for caution regarding the high heritabilities reported for refraction and AL from twin studies, ranging from 75% to 94%,9,11 in view of the much lower heritabilities, 18%8 to 34%,12 obtained from parent–offspring correlations.15 Heritability estimates, in both twin and family studies depend on different assumptions and are likely to be divergent for traits strongly influenced by environmental cues, such as myopia.15,16 Cross-population studies in isolated populations offer the advantage of accessing large complex pedigrees where heritabilities can be drawn simultaneously from the comparison of multiple pairs of relatives. They also benefit from a more stable and uniform diet, climate, and living conditions. However, oculometric traits were analyzed in only a few of those studies.8 The resemblance between distant relatives is less likely to be biased by nongenetic factors, but their genetic covariance is typically small and not well estimated if based on pedigree knowledge only (due to the stochasticity of segregation and recombination). In the present study, we measured the heritability of six ocular biometric traits in two isolated Croatian populations based on realized co-ancestry coefficients drawn from molecular marker information, allowing better estimates of true sharing and thus of heritability. Because of the ethnic variations in ocular morphology, we compared our results only with data derived from populations of European descent. Although numerous studies have shown that ocular biometry can be affected by the amount of near work,17,18 it has also been hypothesized that a diet rich in processed foods plays a role in the increase in juvenile-onset myopia.19 The extent to which body stature and level of education contribute to the values of the traits analyzed was therefore examined. To a large extent, these covariates, including educational achievement,20 are themselves known to have a strong genetic component. Although heritabilities are population specific in principle, in practice they are very similar across populations for morphometric traits and are usually high.21 Estimates of trait heritability in our study should thus inform on the contribution of genetic variants underlying these traits in the studied populations as well as others and guide the choice of covariates to take into account for follow-up gene-mapping studies.