|
|
||||||||
1From the Bloomberg School of Public Health, and the 3Dana Center for Preventive Ophthalmology, School of Medicine, Johns Hopkins University, Baltimore, Maryland; and the 2Inherited Disease Research Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland.
| Abstract |
|---|
|
|
|---|
METHODS. Seven hundred fifty-nine siblings (mean age, 73.4 years) in 241 families were recruited from the Salisbury Eye Evaluation (SEE) Study in eastern Maryland. Refractive error was determined by noncycloplegic subjective refraction (if presenting distance visual acuity was
20/40) or lensometry (if best corrected visual acuity was >20/40 with spectacles). Participants were considered plano (refractive error of zero) if uncorrected visual acuity was >20/40. Preoperative refraction from medical records was used for pseudophakic subjects. Heritability of refractive error was calculated with multivariate linear regression and was estimated as twice the residual between-sibling correlation after adjusting for age, gender, and race. Logistic regression models were used to estimate the odds ratio (OR) of myopia, given a myopic sibling relative to having a nonmyopic sibling.
RESULTS. The estimated heritability of refractive error was 61% (95% confidence interval [CI]: 34%88%) in this population. The age-, race-, and sex-adjusted ORs of myopia were 2.65 (95% CI: 1.674.19), 2.25 (95% CI: 1.313.87), 3.00 (95% CI: 1.565.79), and 2.98 (95% CI: 1.515.87) for myopia thresholds of 0.50, 1.00, 1.50, and 2.00 D, respectively. Neither race nor gender was significantly associated with an increased risk of myopia.
CONCLUSIONS. Refractive error and myopia are highly heritable in this elderly population.
There is also evidence for a recent increase in the prevalence of myopia in populations such as Eskimos, in whom myopia was rare only a few generations ago.4 Such apparent rapid changes in the prevalence of myopia within defined populations argue for strong environmental influences in the development of myopia. For example, much evidence exists that implicates near work as an important determinant of myopia.5 6 7 8 9 10
Nevertheless, there is also considerable evidence that myopia is under genetic control. Studies have consistently shown strong correlations between the refractive errors of first-degree relatives.11 12 13 Twin studies by Hammond et al.,14 Teikari et al.,15 16 and Lyhne et al.17 have suggested that the heritabilitythat is, the proportion of the population variance under genetic controlfor refraction may exceed 80%. It is therefore likely that the development of myopia is mediated by both environmental and genetic risk factors.
At least four chromosomal locations18 19 20 21 22 have been linked to high (that is, in excess of 5 D) myopia. No such loci had been identified for moderate myopia (15 D), which comprises some 80% of myopia among American adults,23 until a recent report by Stambolian et al.24 who found significant genetic linkage of myopia to a region on chromosome 22 in an Ashkenazi Jewish population. Moreover, in a genome-wide scan of 221 dizygotic twin pairs, Hammond et al.25 reported significant linkage signals to multiple genetic loci when analyzing refractive error as a quantitative trait. Their highest linkage peak was observed at chromosomal location 11p13 while evidence for linkage was also noted at regions 3q26, 8p23, and 4q12.
The relative paucity of evidence for linkage in lower levels of refractive error, compared with high myopia, is partly because high myopia may be more likely to present in an autosomal dominant fashion,22 whereas the genetics of lower refractive errors may be more complex and multifactorial. Nevertheless, Farbrother et al.26 argue that even high myopia phenotypes are unlikely to be inherited in an autosomal dominant fashion and that a multifactorial etiology is more likely. In view of the high prevalence (some 20 million American adults23 ) and the resultant significant impact on the health economy in the cost of spectacles, contact lenses, and refractive surgery, it certainly seems that moderate myopia is deserving of genetic study to elucidate its fundamental mechanisms, which may or may not be the same as those of high myopia.
Using linear regression models and generalized estimating equations (GEEs), we estimated the heritability of refractive error in a population-based sample of elderly sibships in eastern Maryland. We also used logistic regression models to calculate the odds of myopia in siblings of myopic individuals compared with siblings of nonmyopes, using four thresholds to define myopia.
| Subjects and Methods |
|---|
|
|
|---|
Methods
Participants underwent standard distance visual acuity testing using Early Treatment Diabetic Retinopathy Study (ETDRS) charts and protocols,27 wearing their habitual spectacle refractions when available. Participants in whom presenting visual acuity was 20/40 or worse in one or both eyes underwent noncycloplegic autorefractometry (Humphrey Autorefractometer model 595; Carl Zeiss Meditec, Dublin CA) followed by binocular subjective refraction. Refractive error was defined for each phakic eye of all subjects as follows: the spectacle lensometry reading for participants whose corrected visual acuities were better than 20/40, or the noncycloplegic subjective refraction yielding best visual acuity for all others. Participants in whom uncorrected distance visual acuity in both eyes was better than 20/40 did not undergo refraction and were assigned a refractive error of zero (plano). For the 154 bilaterally pseudophakic participants, preoperative records were requested from treating physicians, and refractive error was defined for each eye as the spherical equivalent for the last refraction before cataract surgery. In unilaterally pseudophakic subjects, refractive error was defined as described earlier for the sole phakic eye, and the pseudophakic eye was excluded from analysis.
Digital slit and retroillumination lens photographs were taken of all nonpseudophakic eyes of all participants by using a previously described protocol.23 These photographs were graded with the Wilmer Cataract Grading System24 by a team of trained graders.
In addition, participants underwent a slit lamp examination of the anterior segment as well as a dilated fundus examination by the study optometrist (HB). All eyes in which disease (such as corneal opacity) could have affected refractive error or precluded an accurate measurement of the refractive status were excluded from further analysis. The study was conducted in accordance with the tenets of the Declaration of Helsinki and was approved by the Johns Hopkins Medical Institutions Institutional Review Board.
Statistical Methods
Refractive error was analyzed as a continuous variable for heritability estimation. The continuous trait was defined as the spherical equivalent refractive error, averaged between eyes. Spherical equivalent refractive error for the only eye contributing data was used in individuals with monocular pseudophakia or other conditions precluding bilateral refraction. Age was defined as the age at examination or, for bilaterally pseudophakic participants, as the age at the time of last refraction before cataract surgery (averaged between the dates for the two eyes). In addition, analyses were performed using myopia as a binary outcome. For these binary trait analyses, several thresholds were used to define myopia. These included mean spherical equivalent refractive errors lower than or equal to 0.50, 1.00, 1.50, and 2.00 D. Analyses for various definitions of hyperopia were also performed and are reported elsewhere.28
Heritability is estimated from the degree of resemblance in a trait between siblings, or the proportion of additive genetic variance to total phenotypic variance. Mathematically, it is calculated as twice the phenotypic correlation between siblings or, in the case of multivariate regression, as twice the residual between-sibling correlation after adjustment for other variables. We estimated between-sibling correlations using multivariate linear regression models and extended GEEs, with the clustering variable being individual families.25 Only self-reported full siblings were included in the analysis. Main covariates included in the statistical models were age, gender, and race. In addition, we investigated the effect of height, weight, and body mass index (BMI) on refractive error. Cataract grade was not significantly associated with refractive error in models including age and thus was not included in the final model.
Possible modifying effects on refractive errors were assessed by including three first-order interaction terms between age, gender, and race in the regression models. However, interactions between age and gender and between age and race were minimal and are not reported. We also performed stratified analyses to obtain race-specific heritability estimates for refractive error.
Although even large departures from multivariate normality assumptions have been shown to yield unbiased parameter estimates in heritability studies,29 we also calculated the heritability of refractive error after applying a normalizing transformation, transformed refraction = ln(refraction + 12), as suggested by Blackie and Harris.30 The heritability estimate of the transformed data did not differ substantially from that obtained using the raw spherical equivalent refractive error (67% versus 62%). Hence, we report all analyses using the untransformed data.
For the binary trait analysis, we estimated the odds ratio (OR) of being myopic, given a myopic sibling relative to a nonmyopic sibling (i.e., recurrence OR) for our four definitions of myopia. This was accomplished by using logistic regression analysis and GEEs, as described by Liang and Beaty.31 Affection status was coded as a binary variable and covariates, including a first-order gender-race interaction term, were identical with those in the linear regression analysis. No ascertainment correction was used, because participants were recruited independent of refractive status.
Logistic models using GEEs are advantageous as they provide unbiased estimates of sibling correlations and ORs while allowing for the incorporation of covariates and accommodating various family sizes.26 Although the statistical properties of logistic regression lead to estimates of ORs, the parameter of interest for genetic linkage and association studies for binary traits is the sibling recurrence relative risk (
s) which is defined as the risk of being affected, given an affected sibling relative to the risk of being affected in the population (i.e., the population prevalence). ORs are biased estimates of relative risks, and the direction of the bias is always away from the null value of unity. For rare diseases the bias is negligible. However, for more common conditions, such as myopia, the difference can be substantial. We transformed our estimated sibling recurrence ORs into
s using
![]() | (1) |
The statistical analyses were performed in R version 1.7.132 using the GEE library, version 4.13-8. The GEE library was extended to accommodate logistic models and provide sibling recurrence ORs.
| Results |
|---|
|
|
|---|
|
|
The residual between-sibling correlation for refractive error in was 0.310 (95% CI: 0.1750.445). Thus, the estimated heritability for refractive error in this population was 62% (95% CI: 35.0%89%). When analyses were stratified by race, the age- and sex-adjusted estimated heritability of refractive error was 80% (95% CI: 22.9%100%) in African American and 50% (95% CI: 30.4%60.8%) in white participants. However, the heritability estimates for black and white participants did not differ significantly (P = 0.32).
Sibling recurrence ORs for myopia were estimated in logistic regression models that included age (in years), gender, and race variables as well as a gender-race interaction term. The sibling recurrence ORs for myopia ranged from 2.25 (95% CI: 1.373.87) to 3.00 (95% CI: 1.565.79) for the various myopia cutoffs (Table 3) . Estimated
s were 2.00, 1.90, 2.45, and 2.52 for myopia thresholds of 0.50, 1.00, 1.50, and 2.00 D, respectively (Table 3) .
|
| Discussion |
|---|
|
|
|---|
It is important to note that all heritability estimates are population-specific and can show large between-population differences. Specifically, heritabilities are estimated as the proportion of the additive genetic variance to the total (environmental and genetic) variance. Hence, populations in which environmental exposures predominate will have lower heritabilities relative to groups in which these factors are absent, given similar genetic backgrounds. Another possible source of bias in our study relates to the use of correlations between full siblings to estimate heritability. In full siblings, dominance genetic effects in addition to additive effects, contribute to the phenotypic covariance between individuals.34 Hence, heritabilities in studies of full siblings may be inflated relative to those estimated in monozygotic versus dizygotic twin comparisons or parent-offspring designs. We do not think that this bias is likely to have been substantial in this study, since our estimate of heritability is at the lower end of previously reported figures. In addition, substantial differences in environmental exposures may actually reduce the phenotypic resemblance between generations in parent-offspring designs, leading to reduced heritability estimates.
To our knowledge, the heritability of refractive error has never been estimated in African-American populations. The heritability of 80% for the black subcohort in this study was higher than the 50% estimated heritability for the white participants. This may be the result of a higher environmental variance contributing to the phenotype in the white participants or to differing genetic backgrounds between the races. Nevertheless, the CIs of the race-specific heritability estimates overlapped considerably and hence did not differ significantly and the difference in race-specific heritability estimates may have been the result of sampling variation.
The advantage of using linear regression to estimate residual sibling correlation is that covariates can readily be included to account for nongenetic determinants. The addition of extended GEEs (GEE2) allows for the calculation of regression coefficients while accounting for multiple within-family comparisons and families of various sizes.31 Our heritability estimates were age, gender, and race adjusted. Nevertheless, because environmental exposures, such as educational achievement and exposure to prolonged near-vision tasks, are likely to be correlated within sibships and were not completely accounted for in our statistical models, our estimates of heritability may have been artificially elevated in this population.
The mean age of our study population was >70 years. We can expect that nongenetic influences, especially those related to cumulative age effects, would have had a relatively large effect on refractive error in this age group, increasing the total phenotypic variance and weakening the estimated genetic effect. Furthermore, there may have been some phenotypic misclassification in our population, because participants with uncorrected visual acuities better than 20/40 were assigned a null refractive error. Any such misclassification would probably be nondifferential with respect to sibling refractive error and hence decrease the observed heritability estimate.
Physical stature has previously been found to be associated with refractive error.9 35 36 Our results suggest that height may be marginally associated with refractive error, although this association did not reach statistical significance (P = 0.09). Our age-, gender-, and race-adjusted estimate of heritability for height in this population was 92% (data not shown). Nevertheless, the addition of height in our regression models did not change our estimated heritability of refractive error. This suggests that the genetic and/or common environmental factors responsible for physical stature differ from those associated with refractive error.
The age-, gender-, and race-adjusted odds of myopia were, on average, 2.72 times higher in siblings of myopic individuals than in siblings of nonmyopic participants. Our estimates are somewhat lower than those reported by the Lee et al.,13 who obtained ORs of 2.82 to 4.25 (mean, 3.42) for myopia, defined as a spherical equivalent refraction less than 0.50 D, among 100 random samples of sibling pairs from the Beaver Dam population. As with our estimate of heritability, it is likely that the older age of our population compared that of the Beaver Dam cohort (mean age, 61.5 years) would have lead to relatively larger environmental influence and thus a lower estimated OR. Although we did not find a statistically significant relationship between cataract grades and ocular refraction it is possible that more subtle lens changes and differences in the anatomic relationship between ocular components causes age-related variations in ocular refraction.
Using regression models similar to ours, the Framingham Offspring Eye Study Group11 found that the recurrence OR of myopia was related to the age difference of siblings. The OR ranged from 2.50 for siblings whose ages differed by >10 years to 5.13 for siblings whose age differences were within 2 years. In our study, the mean between-sibling age difference was >4 years, and the ages of >20% of our sib pairs differed by 10 years or more. Hence, our recurrence ORs are consistent with Framingham Offspring Study results.
We have also shown that sibling recurrence ORs for myopia are largely invariant to the specific definition of refractive error used. This may reflect underlying quantitative trait loci for refractive error or an environmental dose-response relationship unaccounted for in our statistical models. Nevertheless, it is difficult to assess adequately the significance of this finding, because the tests for different myopia thresholds were correlated (i.e., the same population was used for all measures). Furthermore, we limited our analysis to low-cutoff thresholds for myopia of between 0.50 and 2.00 D, because we did not have a sufficient sample size to use more extreme definitions of the myopic phenotype.
We also provide estimates of
s for myopia.
s is more easily interpreted than the sibling recurrence OR and is the parameter of interest when calculating the power of genetic linkage and association studies for binary traits. The age-, gender-, and race-adjusted risk of myopia was 1.90 to 2.52 higher in siblings of myopic participants than in the general population. Hence, the estimated prevalence of myopia of at least 0.50 D among siblings of myopic individuals was 40%, as opposed to 20% among all study participants. Similarly, at a threshold of 2.00 D, the estimated prevalence of myopia among siblings of myopes was 23%, or 2.55 times the population prevalence of 9%.
It may be that high myopia, which can be transmitted through autosomal dominant genetic mechanisms and for which several chromosomal loci have been identified in genetic linkage studies,18 19 20 21 22 37 clusters in families to a greater extent than the more common, moderate myopia. In a reanalysis of Goldschmidts 1968 population-based sample of Danish teenagers,38 Guggenheim et al.39 calculated
s for high myopia (i.e., refractive error
6 D) to be 20. However, this estimate was based on only 39 high myopes in a population pool of 9243 children. In a study of 296 randomly selected high myopes from British optometric practices, Farbrother et al.26 estimated
s for high myopia to be 4.9 (95% CI: 2.87.6). Unfortunately, they inferred the presence of high myopia among siblings through a questionnaire to be an age of onset of spectacle wear of 9.1 years or younger, which limits the interpretation of their findings.
Our results confirm reports of previous studies that nonpathologic myopia is substantially determined by heredity.4 12 13 14 15 16 However, the genetic mechanism responsible for the development of low myopia is likely to be complex, significantly complicating the search for susceptibility loci. For example, Ashton,40 using segregation analysis, rejected a Mendelian inheritance model for refraction in nuclear families of European and Japanese ancestries.
The current study should be understood in the context of its limitations. Less than complete levels of participation at all levels (in the parent SEE study, or in this nested genetic study) could in theory have led to selection bias with regard to refractive error. However, even if such bias occurred, it would not have biased our estimates of myopia heritability unless myopes with affected siblings were more or less likely to participate than myopes without such sibs, which does not seem likely. Numerous studies have implicated environmental factors,5 6 7 8 10 12 41 42 43 such as prolonged exposure to near work, in refractive error development. In our study, data were not available to assess these effects and their influence on the heritability of refractive error. Efforts should be made to evaluate more fully these risk factors in future genetic epidemiologic studies of refractive error.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication June 23, 2004; revised December 16, 2004; accepted January 4, 2005.
Disclosure: R. Wojciechowski, None; N. Congdon, None; H. Bowie, None; B. Munoz, None; D. Gilbert, None; S.K. West, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Robert Wojciechowski, Inherited Disease Research Branch, National Human Genome Research Institute (NIH), 333 Cassell Drive, Suite 1200, Baltimore, MD 21224; rwojciec{at}jhsph.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. R Rudnicka, C. G Owen, M. Richards, M. E. Wadsworth, and D. P Strachan Effect of breastfeeding and sociodemographic factors on visual outcome in childhood and adolescence Am. J. Clinical Nutrition, May 1, 2008; 87(5): 1392 - 1399. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schache, A. J. Richardson, K. K. Pertile, M. Dirani, K. Scurrah, and P. N. Baird Genetic Mapping of Myopia Susceptibility Loci Invest. Ophthalmol. Vis. Sci., November 1, 2007; 48(11): 4924 - 4929. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Fotouhi, A. Etemadi, H. Hashemi, H. Zeraati, J. E Bailey-Wilson, and K. Mohammad Familial aggregation of myopia in the Tehran eye study: estimation of the sibling and parent offspring recurrence risk ratios Br. J. Ophthalmol., November 1, 2007; 91(11): 1440 - 1444. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Peet, M.-F. Cotch, R. Wojciechowski, J. E. Bailey-Wilson, and D. Stambolian Heritability and Familial Aggregation of Refractive Error in the Old Order Amish Invest. Ophthalmol. Vis. Sci., September 1, 2007; 48(9): 4002 - 4006. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A Guggenheim, R. Pong-Wong, C. S Haley, G. Gazzard, and S. M. Saw Correlations in refractive errors between siblings in the Singapore Cohort Study of Risk factors for Myopia Br. J. Ophthalmol., June 1, 2007; 91(6): 781 - 784. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Tarczy-Hornoch, M. Ying-Lai, R. Varma, and the Los Angeles Latino Eye Study Group Myopic Refractive Error in Adult Latinos: The Los Angeles Latino Eye Study Invest. Ophthalmol. Vis. Sci., May 1, 2006; 47(5): 1845 - 1852. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |