|
|
||||||||
1From the MRC Centre of Epidemiology for Child Health, Institute of Child Health, and 2Institute of Ophthalmology, University College London, London, United Kingdom; 3Twin Research and Genetic Epidemiology Unit, Kings College London, London, United Kingdom; and 4Genome Research Center, The University of Hong Kong, Hong Kong, SAR, China.
| Abstract |
|---|
|
|
|---|
METHODS. Five hundred ninety-six persons from the 1958 British Birth Cohort, a nationally representative population, were randomly selected from the outer tertiles of the refractive error (RE) distribution and were genotyped using 25 tagSNPs across PAX6 and 3 tagSNPs across SOX2 and their putative control regions. This experiment had 80% power to exclude either gene contributing more than 10% of the variance of refractive error.
RESULTS. All SNPs were in Hardy-Weinberg equilibrium, and the genotyping failure rate was less than 5%. Accounting for multiple testing, no significant association (P < 0.05) was found between any of the SNPs or haplotypes and refractive error.
CONCLUSIONS. PAX6 and SOX2 are obvious candidates in RE genetic studies because of their biological roles and prior linkage studies. The present findings strongly suggest refractive error is not directly affected in this population by variants in either gene or by their known promoters/enhancers. The authors suggest that neither PAX6 nor SOX2 should be prioritized in the international search for genetic modifiers of refractive error. Their findings contribute to broader understanding of the pathophysiology of refractive error and highlight the critical role of replication in genetic research on complex disorders.
Refractive error (myopia, hyperopia, or astigmatism) is an archetypal complex quantitative trait; its development is considered to be influenced by multiple genes, interactions between genes, and interactions between genes and environmental factors. There is a particular and longstanding scientific interest in myopia—written accounts date from at least the time of Aristotle4 —that reflects its impact on the lives of affected persons and on the associated societal costs of lost productivity and of comorbidity from vision impairment.5 Recognizing this, the World Health Organization has identified the prevention and treatment of myopia as a priority within its current global initiative against avoidable blindness.6
The prevalence of myopia shows wide geographic variation. Mild to moderate (primary physiological or common) myopia affects approximately 25% of the population in Europe and North America7 8 9 but only approximately 5%9 in Africa and as many as 61% in Asia.10 The risk for sight loss increases with increasing severity of myopia. High myopia (pathologic or syndromic), which may be a distinct entity, affects approximately 2% to 3% of most populations.10 Recent marked increases in the prevalence and the severity of common myopia in many populations11 have been attributed to changing environmental influences. However, there has been a renewed interest in its genetic basis fueled by twin studies,12 13 14 indicating high heritability (60%–90%), and by familial aggregation studies.15 16 17 18
Most genetic studies to date have focused on rare forms of pathologic or syndromic myopia, which may have underlying causes different from those of common myopia. The first genomewide linkage scan to investigate common refractive error was reported in a sample of dizygotic twins drawn from a volunteer twin register.19 Of four linked loci reported, the two of the highest LOD scores were on 11p13 (LOD 6.1) and 3q26 (LOD 3.7). PAX6 underlies the highest point of the peak on 11p13. It plays an essential role in oculogenesis. Mutations in PAX6 cause familial and sporadic aniridia and a number of other severe ocular phenotypes.20 The purpose of its continued postnatal expression in the eye is unknown, but dosage of PAX6 in transgenic mice influences eye size.21 Excessive eye elongation is a key feature of myopia, and this has led to the hypothesis that common, more benign, polymorphisms in PAX6 may be associated with the development of refractive error, though no association has been found to date.19 Thus, the findings of the genomewide linkage scan described earlier have generated a great deal of recent interest in PAX6 as a candidate gene for myopia, and much research effort is being directed toward it.
The 3q26 locus contains the SOX2 gene, a member of the family of sex-determining region Y-box transcription factor genes. SOX2 also has a key role in eye development, and mutations in SOX2 are responsible for syndromic microphthalmia. Although mutations in SOX2 (as with those in PAX6) produce a severe phenotype, its important role in eye development supports the idea that because of common variation, this gene might be a good candidate for study in refractive error.
Thus, we have investigated common variants selected from HapMap across PAX6 and SOX2 and tested their association with refractive error in a well-characterized and nationally representative population.
| Methods |
|---|
|
|
|---|
Refractive error is quantified using the summary measure of spherical equivalent of refraction (SER) in diopters. Although subjects can be divided using thresholds of SER set arbitrarily into categories of myopia (negative SER), hypermetropia (positive SER), and emmetropia (no refractive error), greater power is obtained by treating refractive error as a continuous variable. From autorefraction readings, spherical equivalent of refraction (phenotype) was calculated in the conventional way (SER = S + C/2), where S is the sphere and C is the cylinder. Data on subjects whose discordance between the spherical equivalents in each eye was within the worst 5% of the sample were discarded. This difference was defined as
![]() |
Five hundred ninety-six persons were selected at random from the lowest and highest tertiles (myopic versus hypermetropic) of the mean SER distribution because comparison of extremely opposed phenotypes was considered more beneficial than analysis of the entire trait spectrum.24 25
Refractive error has a leptokurtotic distribution and is skewed toward the myopia end of the distribution (Fig. 1) . This means the hypermetropic tertile includes all those with hypermetropia and some with emmetropia, whereas the myopic tertile covers the full spectrum of myopia severity. Importantly, in analyzing RE as a quantitative trait, it is not the absolute value of SER but the ranking within the distribution that is the variable of interest.
|
The study presented here was part of a broader association study of refractive error, in which 1536 tag SNPs were chosen across 111 candidate genes. Of these, 25 SNPs were selected across a 530-kb region that included PAX6 and putative control regions. Three SNPs were selected across a 10.8-kb region centered on SOX2 and potential control regions. Correction for multiple testing was performed using the false discovery rate.29 This study had 80% power to exclude either PAX6 or SOX2, contributing to more than 10% of the variance of the refractive error in this cohort. Power calculations were performed using the Genetic Power Calculator30 based on an additive genetic variance of 0.1. In the cohort (N = 596) we analyzed and whose results we are reporting, the variance of refraction was 6.25 and the SD from the mean was 2.50. Using these data, we modeled the power to detect, at
= 0.05, any association between genotype and phenotype. The model can be characterized by two key parameters, regression slope and SD, of regression errors.
Results show that our panel was powered more than 80% to detect even minor effects, not only when, as generally expected, deviation from the predicted refraction value was no larger than the SD for the phenotype we observed (2.5) but also when it was twice as large. Our sample has more than 80% power to detect association at a liberal significance level, even in circumstances of small effects from the locus (changes in refraction for each additional copy of susceptibility allele are as small as 0.2–0.3 D). Power to confidently rule out association is more elusive to calculate using most of the standard statistical methods geared at assessing positive association. Power calculations reported assumed an aggressive significance level (
<.01) and did not provide sufficient information for lower contribution to the genetic variance.
| Results |
|---|
|
|
|---|
|
|
| Discussion |
|---|
|
|
|---|
We think it unlikely that potential errors in our study account for this unexpected finding. Unlike many of the early candidate gene studies that were underpowered or overinterpreted and whose findings were irreproducible,31 the design and size of our study are sufficiently robust to allow us to exclude these as explanations for our finding. Despite the undisputed practical importance of the conclusions of the HapMap project, it has been suggested that HapMap SNPs do not necessarily capture variants that have not been genotyped by the HapMap project.32 Therefore, it remains possible that other polymorphisms not typed by HapMap will not be in significant enough linkage disequilibrium (LD) with any of the tagSNPs chosen to be sufficiently captured. Most of the SNPs reported by dbSNP but not typed by the HapMap around PAX6 and SOX2 have low minor allele frequencies or are not polymorphic in Caucasians because both genes are highly conserved. The CD/CV hypothesis suggests that such rare variants are not likely to be of relevance to common traits. Given the extensive research literature on PAX6 and SOX2, it seems unlikely that the failure to identify an association between either of these genes and refractive error results from the existence of yet unidentified common variants. The theoretical alternative to the CD/CV hypothesis is the common disease, rare variant hypothesis (CD/RV),33 where multiple rare variants within genes are responsible for the development of complex traits. This hypothesis was not tested and is not ruled out by this study. There is little empiric evidence for rare variants in most complex traits to date, and we do not think this is a likely underlying paradigm in refractive error.
Our study serves to highlight the critical role of replication in the interpretation of genetic research. Unfortunately, good evidence indicates significant publication bias in favor of "positive" findings if they are novel and contrary to findings that either rule out important hypotheses ("negative" findings) or fail to replicate earlier work and thereby demonstrate prior false discovery. Notably, it is generally agreed34 35 that regardless of the relative quality of the science, it takes longer to achieve publication of negative results, and it is difficult to get them published in journals with higher impact factors. Indeed such findings are sometimes not even submitted for publication. Instead they are labeled "file drawer problem"36 in anticipation of the difficulties. Promotion of false gene discovery in this way is certain to mean that research efforts and funds are being directed to areas unlikely to be fruitful.
PAX6 and SOX2 were obvious candidates for refractive error, based on knowledge of biological roles and linkage data. We recognize that the lack of association in our study does not preclude PAX6 or SOX2 involvement in refractive error caused by mechanisms such as regulation of expression by environmental factors or undiscovered cis- or trans-acting elements. However, from the findings of the present study, we suggest that neither PAX6 nor SOX2 should be prioritized in the active international search for risk and phenotype modifiers of refractive error.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication February 22, 2007; revised April 23, 2007; accepted August 13, 2007.
Disclosure: C.L. Simpson, None; P. Hysi, None; S.S. Bhattacharya, None; C.J. Hammond, None; A. Webster, None; C.S. Peckham, None; P.C. Sham, None; J.S. Rahi, 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: Jugnoo S. Rahi, MRC Centre of Epidemiology for Child Health, Institute of Child Health, University College London, London WC1N 1EH, United Kingdom; j.rahi{at}ich.ucl.ac.uk.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Han, K. H. Leung, W. Y. Fung, J. Y. Y. Mak, Y. M. Li, M. K. H. Yap, and S. P. Yip Association of PAX6 Polymorphisms with High Myopia in Han Chinese Nuclear Families Invest. Ophthalmol. Vis. Sci., January 1, 2009; 50(1): 47 - 56. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |