|
|
||||||||
1From the Statistical Genetics Section, Inherited Disease Research Branch, National Human Genome Research Institute, Baltimore, Maryland; and the 2Department of Ophthalmology and Visual Sciences, University of Wisconsin Medical School, Madison, Wisconsin.
| Abstract |
|---|
|
|
|---|
METHODS. Familial correlations were examined and segregation analysis was performed on the average refractive error measurements in the right and left eyes after adjustment for age, sex, and education. Analyses were based on 2138 individuals in 620 extended pedigrees with complete data on age, sex, education, and spherical equivalent.
RESULTS. Substantial positive correlation was found between siblings (0.33), parents and offspring (0.17), and cousins (0.10) and lower correlation among avuncular pairs (0.08) after adjustment for age, sex, and years of education. The results of this segregation analysis do not support the involvement of a single major locus throughout the entire range of refractive error. However, models allowing for familial correlation, attributable in part to polygenic effects, provided a better fit to the observed data than models without a polygenic component, suggesting that several genes of modest effect may influence refractive error, possibly in conjunction with environmental factors.
CONCLUSIONS. These results support the involvement of genetic factors in the etiology of refractive error and are consistent with reports of linkage to multiple regions of the genome.
Among the more distant family members, it has long been established that there is correlation in refractive error measurements.9 However, the only reported segregation analysis of refraction did not support the influence of a single major gene on refraction.10 These analyses did not take into account age, sex, and education effects. The age distribution of this study population was much broader (<10 to
70 years) compared with the Beaver Dam Eye Study (4384 years). In addition, Ashton10 transformed the data before analysis, unlike in the current study, in which the best-fitting transformation is estimated as part of the model fitting process of segregation analysis.
Environmental risk factors have also been associated with refractive error, myopia, or hyperopia. Previous studies, including studies conducted within the cohort used for these analyses, The Beaver Dam Eye Study, have shown that both age and sex are associated with refractive error.11 12 13 14 Education15 16 17 and near-work18 are both strongly associated with increasing severity of myopia. Therefore, these environmental factors must be taken into account when examining familial risk.
The goal of this study was to expand on our previous findings that refractive errors are moderately correlated among relative pairs in the following ways. First, we wanted to determine whether the observed strong familial correlations of refractive error measurements would remain after adjustment for age, education, and sex. Second, through the use of complex segregation analysis, we assessed whether these correlations between family members were due to shared environmental effects, genetic effects, or a combination of the two.
| Methods |
|---|
|
|
|---|
Study Population
The baseline exam of The Beaver Dam Eye Study, conducted between 1988 and 1990, involved 4926 participants (of the 5924 eligible individuals, who resided in the township of Beaver Dam).19 Follow-up examinations have been conducted every 5 years. However, in the present study, we used only data from the baseline examination. Recruitment methods and study procedures are described in detail elsewhere.20
During the baseline and subsequent examinations, eye examinations were performed, including automated refractive error measurements for all participants. Family relationship information was obtained from all participants at the baseline examination. During the first follow-up visit, conducted between 1993 and 1995, family relationships, including extended pedigree information were confirmed. Of the 5924 eligible individuals, 2783 had available information on familial relationships and could be classified into one of 602 pedigrees. Of these individuals, 2138 had complete age, sex, education, and refractive error data.
However, due to the limitations of the software used to analyze these data, several of the more complex pedigrees were split into smaller pedigrees yielding the 620 pedigrees used in these analyses. In this process, no individuals for whom refraction measurements were available were duplicated and only distant relationships (i.e., cousin pairs or more distant) were severed. The resultant pedigrees did not have any individuals included more than once.
Measurement of Refractive Error
Automated refractive error measurements were obtained from 96% of the eyes. When data were available from Early Treatment Diabetic Retinopathy Study (ETDRS) refraction, these measurements were used in the analysis (4% of eyes). When data from neither of these refractions were available, (<1% of eyes), refraction from the current prescription was used. Eyes without a lens, with an intraocular lens, or eyes with best corrected visual acuity of 20/200 or worse were excluded. Only individuals with data on both eyes were included in the analysis. Spherical equivalent (sphere power + [0.5 · cylinder power] measured in diopters) was calculated from the refraction measurements. The average of the spherical equivalent in the right and left eyes was used in these analyses.
Statistical Analysis
Familial correlation analysis was performed with FCOR, part of the Statistical Analysis for Genetic Epidemiology S.A.G.E. (S.A.G.E., ver. 4.5) statistical package. Correlations in refractive error measurements were calculated between the following relative pairs: parents and offspring, sibling, avuncular, and cousin. Equal weight was given to each pair of relatives.21
Commingling analysis and segregation analysis were performed with REGC version 2.1 and REGCHUNT.22 REGC is part of the S.A.G.E. version 2.1 statistical package.23 For both commingling and segregation analysis, Box-Cox transformation of the data was estimated as part of the analysis, denoted by parameters
1 and
2, to ensure data were on the proper scale.23 24 For all analyses
1, the power parameter, was freely estimated, whereas
2, the scale parameter, was fixed to 20.5 (to ensure all that adjusted trait values were non-negative before power transformation).
Commingling analysis (fitting mixtures of distributions) was used to determine whether there was a single normal distribution (described by mean and variance denoted µ and
, respectively) that provided an adequate description of the data or if a mixture of two or three normal distributions provided a significantly better description of the data. The proportion of the population in each of the distributions is denoted by
.
Segregation analysis involves fitting a series of genetic and nongenetic models, both with and without polygenic components, to determine whether there is evidence of a major gene or polygenic components that influence refractive error. If the genetic models, models in which these parameters are fixed to what is expected under Mendelian assumptions, describe the data as well as more general models, there is support for the involvement of a single major gene associated with refractive error. The specific parameters that comprise these models are detailed herein.
REGC uses the regressive models proposed by Bonney23 25 26 to perform segregation analysis of a continuous trait. Using these models, we tested for autosomal inheritance of a single biallelic major locus that influences refractive error as a quantitative trait by obtaining maximum-likelihood estimates for parameters designed to describe the distribution of refractive error in this population. The underling type of each individual is estimated. This represents an underlying discrete factor that influences refractive error.27 In the models that test for inheritance of a major gene, "type" represents a genotype, but in models that test for nongenetic factors, "type" (denoted AA, AB, and BB) is interpreted as levels of exposure to an unmeasured major environmental risk factor that is not correlated between family members. Three possible types are considered, which for Mendelian inheritance represent the two homozygotes AA and BB and the heterozygote AB. The mean (µ) and variance (
) of the refractive error phenotype for each of the types is estimated. However, given that these types must sum to 1, only two parameters are estimated, denoted qA and qB. When Hardy-Weinberg equilibrium is assumed, only a single parameter qA is estimated. In addition, transmission parameters (denoted by
), are estimated to test whether type is shared between parent and offspring in the proportions expected under Mendelian assumptions or whether this sharing shows patterns consistent with major random environmental exposures. The probability of a parents transmitting factor A (or A allele for genetic models) to an offspring is represented by these parameters. Under Mendelian expectations, transmission parameters are fixed to 1 for individuals of type AA, 0.5 for individuals of type AB, and 0 for individuals of type BB, denoted as
AA,
AB, and
BB, respectively. When these parameters are set equal to each other or equal to the type frequencies, an environmental model is represented.
These analyses were performed under the assumptions of a class-D model: Dependency between sets of siblings is equalthat is, it is not affected by birth order or other factors, but is not due to common parentage alone. In addition, familial correlations can be estimated within these analyses to account for other genes of small to modest effect (polygenes) or additional environmental factors that are shared among family members. Only correlation between members of nuclear families is estimated including: spousal (
fm), parent and offspring (
po), and sibling (
ss). In addition, because age, sex, and education are known to influence refractive error, they were incorporated in the analysis. Age and education effects were included as covariates (influencing mean values). The variance of refractive error was allowed to be different between men and women.
The most parsimonious model that adequately described these data was selected by using likelihood ratio tests and Akaikes information criterion A (AIC). Likelihood ratio tests were computed by 2 times the difference in lnLikelihood between the general model and a smaller model. This was then compared to a
2 distribution, in which the degrees of freedom were equal to the difference in the number of parameters estimated between the two models. In the situation when the parameters in the general model maximized at a boundary, a mixture of
2 distributions were used to compute probabilities.28 AIC allowed us to compare non-nested models by taking the 2 lnLikelihood of the model plus a correction of 2 times degrees of freedom of the model for estimation of additional parameters.29 The minimum AIC indicates the most parsimonious model.
No ascertainment correction was necessary for these data, given that they were obtained through a population-based survey.
| Results |
|---|
|
|
|---|
|
|
|
| Discussion |
|---|
|
|
|---|
Commingling analysis of these data indicated that fitting multiple distributions to the data provided a better fit to the data than did a single distribution. This suggests there is not a single distribution of refractive errors in the population, but several (at least three) that could be determined by either genetic factors, major environmental factors, or a combination of both.
The results of our segregation analysis indicated that neither a single-distribution model with polygenetic effects nor models in which there was a major environmental effect along with polygenes provided an adequate fit. However, the environmental models, which allowed for multiple underlying distributions, did provide a better fit to the data than did the single-distribution model. The models that incorporated a major gene effect with additional polygenetic effects provided a better fit to the data than did the environmental models. However, these Mendelian models still did not adequately fit the data and were rejected when compared with the general transmission models. This lack of fit of the Mendelian models may be due to a cohort effect, as is indicated in the familial correlation analysis, or to other environmental factors that we have not adjusted for in the analysis or to a very complex mix of several major loci, polygenes, and environmental factors that all influence the variability of refraction in this population. Because the models we tested only examined the possibility of a single major gene that controls refraction across the entire range of refractive values (i.e., from myopia to hyperopia) the lack of fit of the single-gene Mendelian models could indicate that multiple gene(s) influence this trait. In addition, models that included a polygenic component provided a better fit to the observed data than did models that did not, which supports the involvement of several genes of modest effect in the etiology of refractive error. These results are consistent with the twin studies by Hammond et al.7 8 and the segregation analysis by Ashton.10
Previous studies, including studies conducted within this cohort have demonstrated that refraction changes with age.11 12 31 In the Beaver Dam Eye Study population, the 10-year change in refraction was approximately +0.5 D in individuals aged 43 to 59, which contrasted with a 0.41-D 10-year change in individuals aged
70. Little change was observed for individuals aged 60 to 69 in the 10-year follow-up period. This change in refraction over a 10-year period did not seem to differ between myopes and hyperopes. Therefore, although individuals who are classified as myopic at younger ages may no longer meet the criteria for myopia at older ages, they will still (on average) have lower refractive error than an individual of the same age who was never myopic. Thus, given that we examined refractive error as a continuous trait and adjusted for the effect of age, we feel there is minimal misclassification in these data.
We may not have removed all of the effect of age, because we assumed the relationship with age to be linear, when from the longitudinal studies we know the relationship between age and refraction is not entirely linear. We attempted to include an age-squared term in the analysis, but inclusion of this covariate overparameterized the model. However, the nonlinearity is strongest in the oldest age groups, and only 6.7% of the study population was >80 years of age.
Although the Beaver Dam Eye Study was designed as a population-based study, given that this study was conducted in a small town in Wisconsin and that a high proportion of the participants were related within extended pedigrees, provides us the unique opportunity to study the genetics of refractive error in a population of families not ascertained based on refractive error measurements. The population of the entire Beaver Dam Eye Study was comparable to those in other U.S. cities of its size for income, occupation, sex distribution, and education attainment, as described in the 1980 census (BEKK, personal communication, 2003). The current analysis was limited to study participants who were also members of families. Participants who could be classified into families tended to be a bit older, were more likely to be male, had a lower level of education, and were more likely to have nuclear cataract than the entire Beaver Dam Eye Study cohort.30 No data were obtained from family members who did not reside in Beaver Dam.
Although no genes for refraction, myopia, and hyperopia have been identified, linkage has been reported to regions on chromosomes 12, 17, and 18 in studies of families with extremely high values of refraction (i.e., high myopia), on chromosome 22 in families in which there is a large degree of aggregation of moderate myopia, and on chromosomes 3, 4, 8, and 11 for refractive error as a quantitative trait in a cohort of dizygotic twins. In addition to the possibility that multiple genes act to influence overall refraction, these results may suggest that some of the genes that influence myopia are distinct from the genes that influence hyperopia. Linkage and association studies of the entire Beaver Dam Eye Study family resource, designed to localize the genes involved in refractive error, myopia, and hyperopia, are currently under way.
| Footnotes |
|---|
Submitted for publication July 7, 2004; revised October 18, 2004; accepted October 20, 2004.
Disclosure: A.P. Klein, None; P. Duggal, None; K.E. Lee, None; R. Klein, None; J.E. Bailey-Wilson, None; B.E.K. Klein, 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: Barbara E. K. Klein, Professor Ophthalmology and Visual Sciences, University of Wisconsin-Madison, Ocular Epidemiology, 610 North Walnut Street, Room 409, Madison, WI 53726; kleinb{at}epi.ophth.wisc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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. M. Ip, S. C. Huynh, D. Robaei, K. A. Rose, I. G. Morgan, W. Smith, A. Kifley, and P. Mitchell Ethnic Differences in the Impact of Parental Myopia: Findings from a Population-Based Study of 12-Year-Old Australian Children Invest. Ophthalmol. Vis. Sci., June 1, 2007; 48(6): 2520 - 2528. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Gupta, G. V. S. Murthy, N. Morrison, G. M. Price, M. Dherani, N. John, A. E. Fletcher, and U. Chakravarthy Prevalence of Early and Late Age-Related Macular Degeneration in a Rural Population in Northern India: The INDEYE Feasibility Study Invest. Ophthalmol. Vis. Sci., March 1, 2007; 48(3): 1007 - 1011. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. L. Young, R. Metlapally, and A. E. Shay Complex Trait Genetics of Refractive Error Arch Ophthalmol, January 1, 2007; 125(1): 38 - 48. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. P. Klein, P. Duggal, K. E. Lee, R. Klein, J. E. Bailey-Wilson, and B. E. K. Klein Confirmation of Linkage to Ocular Refraction on Chromosome 22q and Identification of a Novel Linkage Region on 1q Arch Ophthalmol, January 1, 2007; 125(1): 80 - 85. [Abstract] [Full Text] [PDF] |
||||
![]() |
S-M Saw, P-P Goh, A Cheng, A Shankar, D T H Tan, and L B Ellwein Ethnicity-specific prevalences of refractive errors vary in Asian children in neighbouring Malaysia and Singapore Br. J. Ophthalmol., October 1, 2006; 90(10): 1230 - 1235. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Klein and B. E. K. Klein The Epidemiology of Eye Disease: From Glycemia to Genetics The Friedenwald Lecture Invest. Ophthalmol. Vis. Sci., May 1, 2006; 47(5): 1747 - 1753. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |