Molecular Vision 2011; 17:2028-2039 <http://www.molvis.org/molvis/v17/a221>
Received 24 March 2011 | Accepted 18 July 2011 | Published 22 July 2011

Identification of novel suggestive loci for high-grade myopia in Polish families

Malgorzata Rydzanicz,1 Swapan K. Nath,2 Celi Sun,2 Monika Podfigurna-Musielak,3 Agata Frajdenberg,4,5,6 Malgorzata Mrugacz,7 Daniel Winters,8 Uppala Ratnamala,9 Uppala Radhakrishna,9 Bassem A. Bejjani,10 Marzena Gajecka1,11

1Institute of Human Genetics, Polish Academy of Sciences, Poznan, Poland; 2Arthritis and Immunology Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK; 3Department of Ophthalmology, Leszno Hospital, Leszno, Poland; 4Department of Ophthalmology, Marcinkowski University of Medical Sciences, Poznan, Poland; 5Namsos Hospital, Department of Ophthalmology, Namsos, Norway; 6University Hospital in Linköping, Department of Ophthalmology, Linköping, Sweden; 7Department of Pediatric Ophthalmology, Medical University of Bialystok, Bialystok, Poland; 8School of Molecular Biosciences, Washington State University, Spokane, WA; 9Department of Surgery-Transplant, University of Nebraska Medical Center, Omaha, NE; 10Signature Genomic Laboratories, LLC, Spokane, WA; 11Basic Medical Sciences Program, WWAMI, Spokane, WA

Correspondence to: Marzena Gajecka, Ph.D., Institute of Human Genetics, Polish Academy of Sciences, Strzeszynska 32, Poznan, 60-479, Poland; Phone: (061) 657-9160; FAX: (061) 823-3235; email: gamar@man.poznan.pl

Abstract

Purpose: Myopia is the most common human eye disorder with complex genetic and environmental causes. To date, several myopia loci have been identified in families of different geographic origin. However, no causative gene(s) have yet been identified. The aim of this study was the characterization of Polish families with high-grade myopia, including genetic analysis.

Methods: Forty-two multiplex Polish families with non–syndromic high-grade myopia participated in the study. All family members underwent detailed ophthalmic examination and high-grade myopia was defined as ≤-6.0 diopters (D) based on the spherical refractive error. A genome-wide single nucleotide polymorphism (SNP)-based high-density linkage scan was performed using Affymetrix Human SNP Array 6.0 on a selected family (HM-32) with multiple affected individuals.

Results: Nonparametric linkage analysis identified three novel loci in family HM-32 at chromosome 7p22.1–7p21.1 ([NPL] 8.26; p=0.006), chromosome 7p12.3–7p11.2 ([NPL] 8.23; p=0.006), and chromosome 12p12.3–12p12.1 ([NPL] 8.02; p=0.006), respectively. The effect of linkage disequilibrium on linkage due to dense SNP map was addressed by systematically pruning SNPs from the linkage panel.

Conclusions: Haplotype analysis with informative crossovers in affected individuals defined a 12.2; 10.9; and 9.5 Mb genomic regions for high-grade myopia spanned between SNP markers rs11977885/rs10950639, rs11770622/rs9719399, and rs4763417/rs10842388 on chromosomes 7p22.1–7p21.1, 7p12.3–7p11.2, and 12p12.3–12p12.1, respectively.

Introduction

Myopia, also known as shortsightedness, is the most common eye disorder worldwide. In myopic subjects, the image of distant objects falls in front of the retina, either as the eye is too long (axial myopia), the cornea is too convex or the index of refraction of the lens is too high (refractive myopia) [1]. The myopic eye is generally vulnerable and persons with ≤-6.0 diopters (D) are more liable to a wide range of ocular pathologies. The development of high-grade myopia involves anterior-posterior enlargement of the eye, scleral thinning, changes in the diameter of scleral collagen fibrils, and frequent detachment of the retina resulting from stress related with axial elongation [2].

The estimated prevalence of high grade myopia is ~2.5 to 9.6% in the elderly world population [3,4]. However, its highest prevalence rates are in Asians, in whom almost 50 to 80% of the adult populations are myopic [5-7]. Recent population-based studies suggest that the prevalence is increasing, specifically in Asian populations [8,9]. The frequency of myopia in the Polish population is unknown, and there is a paucity of data about the epidemiology of high-myopia in Poland. Until the present study, no analysis has yet been made on familial high-grade myopia. However, in Poland the main cause of blindness and ~12% childhood visual impairment is due to high-grade myopia [10].

Myopia may be of diverse etiology, including environmental and genetic factors [11-17]. However, high-grade myopia is highly heritable and genetic predispositions seem to be a dominant factor of its development and progression [18,19]. Families with autosomal dominant, autosomal recessive, and X-linked inheritance of high-grade myopia have been described, though the majority of the reports deal with the autosomal dominant form [20-28]. Although there were no obvious phenotypic differences between the affected subjects of families used in the previously published linkage analyses, the data were inconsistent, suggesting genetic heterogeneity among populations of different geographic origin [20,23]. Several studies have demonstrated a significant genetic component in the familial aggregation of high-grade and/or moderate myopia. Up to 17 loci on different chromosomal regions have been identified [20,22,24-37]. Additionally, increased incidence of single nucleotide polymorphism (SNP)s association between the insulin-like growth factor 1 (IGF-1) gene and high-grade myopia [38], and polymorphisms in the promoter regions of matrix metalloproteinase (MMP) genes were reported [39]. Moreover, two recent independent genome-wide association studies conducted on large cohorts of refractive error patients identified loci at chromosome 15q14 and 15q25 and suggest that the genetic variance in refractive errors is mostly determined by multiple variants with a low to moderate penetrance [40,41].

In this study we present clinical characteristics of forty two Polish families with non-syndromic high-grade myopia and the results of a high-density SNP-based linkage analysis for one selected large high-grade myopia family with multiple affected and normal individuals. Our findings provided evidence of suggestive linkage at three distinct novel loci on chromosome 7p22.1–7p21.3, 7p12.3–7p11.2, and 12p12.3–12p12.1 in the analyzed family.

Methods

Recruitment and clinical evaluation of high-grade myopia families

The study population consisted of 42 multiplex high-grade myopia families from Poland, who were ascertained at three independent Polish institutions: 1) Department of Ophthalmology, Marcinkowski University of Medical Sciences, in Poznan, 2) Department of Pediatric Ophthalmology, University of Medical Sciences in Bialystok, and 3) Department of Ophthalmology, Hospital, Leszno. A constant clinical evaluation procedure was applied at all clinical sites. Informed consent was obtained from all study subjects after the possible consequences of participating in the study were explained, in accordance with the Declaration of Helsinki.

All study subjects underwent a detailed ophthalmic evaluation using computer-assisted equipment included: a visual acuity testing, best-corrected visual acuity testing, a slit lamp evaluation, intraocular pressure examination, fundoscopy, axial length determination, keratometry and refractometry. Biometric axial length (including anterior-chamber depth, lens thickness, and total axial length) was measured using ultrasonography (A, OPTOPOL, Desmin F/H, version 2.06.21). In children ≤15 years old, the refractive error was measured with an autorefractor after cycloplegia. A complete questionnaire was filled for each subject with clinical and family history.

To minimize misclassification, clear diagnostic criteria were established for all high-grade myopia study subjects including spherical refractive error analysis. The subjects were classified into three groups, 1) Affected individuals with high-grade myopia, 2) Individuals with an unknown status and 3) Unaffected persons. All affected individuals showed: 1) bilateral axial high-grade myopia, in excess of or equal to −6.0 D (≤-6.0 D) for at least one eye and in excess of or equal to −5.0 D (≤-5.0 D) for the second eye; 2) a history of onset of myopia at age ≤15 years, and 3) individual with affected status while high-grade myopia was identified in multiple members of their family in different generations. Individuals who were classified as unknown were: 1) all children ≤15 years unless they fulfill criteria for affected status as specified above, or 2) individuals who have myopia with −6.0 D < X ≤ −4.0 D, or 3) individuals, with a refractive error of ≤-6.0 D for one eye and a refractive error >-5.0 D for the second eye, or 4) individuals with late age of onset (>15 years). All the remaining were treated as unaffected as neither of them were classified as affected nor unknown for the analysis.

For all 42 Polish HM families we have performed the analysis using microsatellite markers to exclude or confirm linkage with known high myopia loci (data not shown). In all families previously suggested candidate loci/genes for high myopia were excluded (data not shown).

Statistical analysis in clinical evaluation

Differences in ophthalmic parameters obtained for respective groups, as well as comparison of age were analyzed by the Kruskal–Wallis test [42]. Gender distribution was calculated by χ2 test. All analyzed features were compared among groups according to the scheme: affected versus unaffected, affected versus unknown and unaffected versus unknown. The differences between examined groups were considered significant if the value of probability (p) did not exceed 0.05. Axial length in affected individuals helped reveal whether a patient had corrective surgery in the past. Affected subjects who underwent corrective surgery were not included in analysis of mean refraction value for high-grade myopes versus non-highly myopic subjects.

Genome-wide genotyping in family HM-32

The family HM-32 was chosen for genome-wide genotyping analysis. The selected pedigree was the largest, multigenerational, representative family with many available family members, including patients with high-grade myopia, as well as unaffected relatives.

A genome-wide SNP-based high-density linkage scan was performed using the Affymetrix Human SNP array 6.0 (Affymetrix Inc., Santa Clara, CA) which features 1.8 million genetic markers, including 906,600 SNPs and 946,000 probes for the detection of copy number variation. The assay was performed with 500 ng of genomic DNA, and more than 99% of the SNPs were determined unequivocally for each sample. Scanned images were processed with gene microarray software (Affymetrix) and the data was analyzed (GDAS ver. 2, software; GeneChip Data Analysis; Affymetrix). PEDCHECK [43] was used to identify Mendelian inconsistencies, and the MERLIN [44] program was used to detect double recombination events over short genetic distances that were probably due to genotyping errors. After the quality control (QC) of the raw genotype data which deleted SNPs with missing genotypes and all SNPs at which all individuals have the ‘BB’ or ‘AA’ genotypes, there were a total 550,441 SNPs left for analysis. The genotypes and the markers generated from the Affymetrix 6.0 platform were so dense that linkage disequilibrium between many markers will result in severe biases in linkage calculations. Since we analyzed one family with 16 samples genotypes available, we used the allele frequencies from HapMap CEU population for the analysis. To prevent bias from linkage disequilibrium (LD) in linkage calculations, we first created genotype subset with 4,417 SNPs for the genome wide linkage scan by selecting one SNP from every 100,000 bases in the QC genotype data set. While selecting SNPs, the minor allele frequency (MAF) >1% at each SNP was also a criterion for selection. To maximize the heterozygosity, we always selected the SNPs with high MAF. In a case of identified candidate interval(s) further analysis was performed, selecting one SNP from every two SNPs from candidate interval(s) for the analysis to eliminate the bias from LD. Moreover, linkage analysis was also done by another way to account for LD effect by using genetic distances, selecting one SNP per 0.5 cM, 1 cM, 1.5 cM, and 2 cM for the linkage analysis, respectively.

SNP genotype data were imported into the linkage-analysis programs GENEHUNTER [45] and MERLIN [44]. In the initial genome scan, evidence of linkage was assessed with a nonparametric, penetrance-independent, affected-only, and allele-sharing analysis (Z-mean from MERLIN and nonparametric linkage (NPL) from GENEHUNTER). With MERLIN, one can convert this into a nonparametric logarithm (base 10) of odds (LOD*) score by maximizing the likelihood with respect to a scalar parameter, δ, that measures the amount of excess sharing of identical-by-descent alleles among affected relatives (with δ=0) corresponding to the null hypothesis of no linkage [46]. We used the Sall scoring function and the exponential allele-sharing model to generate the relevant linkage statistic. When significant evidence of linkage was found by exceeding the predetermined threshold (p<0.01), two-point as well as multipoint LOD scores maximized over various plausible genetic model parameters were calculated. For the parametric linkage analysis the best model was estimated as an autosomal dominant mode of inheritance with reduced penetrance (0.6) and phenocopy rate (0.01) and a disease allele frequency of 0.0001. In addition, for the parametric linkage analysis an affected only analysis was performed under an autosomal dominant mode of inheritance allowing for phenocopies. Genetic map distances were derived from the Rutgers combined linkage-physical map of the human genome [47], either directly or by interpolation. Haplotypes were reconstructed using both GENEHUNTER and SIMWALK2 programs [48,49].

Results

Clinical and demographic characteristics of studied families

The forty-two large Polish pedigrees enrolled in the study, had families with five-generation (n=10), four-generations (n=7), three-generations (n=27), and two-generations (n=7) with an average number of individuals in each generation per family of 3.5 (Figure 1, Figure 2, and Appendix 1). The mean family size was 8.3 individuals (range 3–45), with the average affected individuals per family of 3.0 (range 2–10), unaffected 4.1 (range 1–28) and unknown status 0.8 (range 1–7). The specified information was based on individuals who underwent an ophthalmic examination.

A complete eye examination was performed for 331 participated individuals. In accordance with our classification criteria, 128 individuals were considered as affected, 171 as unaffected, and 32 with unknown status. The characteristics and details of the ophthalmic examinations in particular study groups are given in Table 1. Additionally, Appendix 2 shows detailed clinical findings in family HM-32. Spherical refractive error alone (without cylindrical refractive error) was enough to classify 128 affected individuals in this category.

There are several individuals in this study with unknown disease status, which partly due to the involvement of children <15 years. In 50% of cases the unknown disease status is due to the inclusion of individuals with an average spherical refractive error (SPH) ranges between −6.0 D < X ≤ −4.0 D. For example, in families 10, 14, 39 and 75, individuals 10-III: 3, 14-II: 2, 39-III: 4 and 75-I: 4 presented SPH as follow: −4.75/-4.75 D, −4.5/-5.25 D, −4.25/-4.25 D and −4.5/-6.0 D, respectively.

Based on medical records and/or self-reports the average age of onset in myopic subjects was ~8 years (range 2–15). Affected females had slightly earlier onset than affected males (7.78 years, range 2–14 versus 8.88 years, range 2–15, respectively), however the difference was not statistically significant (p=0.077). Some of the affecteds were found with various associated anomalies including glaucoma (11.7%; n=15), cataract (4.7%; n=6), retinal detachment (RD; 5.5%; n=7), and RD in both eyes (n=2). In both unknown and unaffected individuals no other anomalies were identified except one normal individual with glaucoma (0.6%).

We have found statistically significant differences in spherical refractive error, spherical equivalent refractive error (SE) and axial length (AL) between the studied groups (p<0.001; Appendix 3). In affected subjects the average SPH was −9.32 D (±3.89), where the median value was −7.87 D, compared with −2.82 D (±2.02, median −4.00 D) for unknown and 0.00 D (±1.23, median 0.00 D) for unaffected individuals. The average spherical equivalent refractive error of the affected was −9.69 D (±4.00, median −8.25 D), where for unknown and unaffected respectively, −2.91 D (±2.00, median −3.87 D) and −0.19 D (±1.22, median −0.12 D). The mean AL (27.27±2.03, median 26.95 mm) of high-grade myopic eyes was significantly higher than observed for unknown (24.31±1.14, median 24.18 mm) and unaffected individuals (23.36±0.83, median 23.37 mm; Table 1 and Appendix 3). There was no difference in keratometry among the studied groups.

Genome-wide genotyping in family HM-32

Genome-wide (SNP)-based high density linkage analysis in a large high-grade myopia family (HM-32) has shown evidence for susceptibility loci on chromosomes 7p22.1–7p21.1 ([NPL] 8.26; p=0.006); 7p12.3–7p11.2 ([NPL] 8.23; p=0.006) and 12p12.3–12p12.1 ([NPL] 8.02; p=0.006: Table 2 and Figure 3) suspected to be implicated in high-grade myopia in the present family. These results were also supported by multipoint parametric linkage analysis with maximum multipoint LOD scores of 1.33, 2.12, and 1.81 for the respective chromosomal regions. Haplotypes were constructed using over 30 informative SNP markers on three variant linked chromosomal regions. Segregation analysis identified a risk haplotypes that was shared by all affected individuals at three loci and was not found in any of the unrelated unaffected spouses. Critical recombination events that occurred in two affected patients M-32–31 and M-32–15 allowed us to define a 21.4 cM disease interval delimited by SNPs rs11977885 and rs10950639 on chromosome 7p22.1–7p21.1 (Figure 4A). Informative recombination events in affected individuals M-32–04, M-32–15, M-32–19, and M-32–20, confined the myopia candidate locus on chromosome 7p12.3–7p11.2 to a region of 9.3 cM between SNPs rs11770622 and rs9719399 (Figure 4B). Similar haplotype analysis and critical recombination events across the affected family members M-32–15 and M-32–20 on chromosome 12p12.3–12p12.1 narrowed the genomic region to 12.6 cM (Figure 4C). The area is bordered by proximal marker rs4763417 and distal marker rs10842388 (Appendix 4A-C). We also tried another way to prune out the LD by using genetic distances, selecting one SNP per 0.5 cM, 1 cM, 1.5 cM and 2 cM for the linkage analysis. The results were similar compared to what we described earlier except for the second interval 68 cM- 80 cM on chromosome 7p13–7p11.2, which shows a change from 0.5 cM to 1 cM, but still remains significant P value (Appendix 5).

The SNP-based copy number analyses did not reveal any variation in the linked genomic regions.

Discussion

Myopia is the most common ocular disorder in all populations and the incidence is increasing in all parts of the world [3]. After the first genetic linkage studies of X-linked syndromic form of myopia, many loci were mapped in highly-selected families that aggregated severe forms of myopia in diverse populations of various geographic origins and some attempted to replicate the data [20,22,24,50-53]. However, to date no disease-causing mutation(s) in any gene(s) have been identified. Most of the myopia loci were identified using family-based linkage analysis with microsatellite markers and the only available recent genome-wide association studies (GWAS) using SNP arrays, in a large cohort of refractive error patients identified two different loci [40,41].

The present study provides further evidence for genetic heterogeneity and indicates that more than one locus may contribute to high myopia. It also suggests that the high-grade myopia in Polish families is not allelic to any of the previously described candidate loci (personal communication, M.G.). In addition, exclusion of newly identified candidate loci in other Polish families indicates possible genetic heterogeneity within Polish population signifying that genome-wide linkage analysis in these families may reveal novel locus/loci for high myopia. Naiglin et al. [21] also reported genetic heterogeneity in families with high myopia. Earlier linkage for high myopia was reported on chromosome 7p and 12p; however these loci do not overlap with the genomic regions identified in the present family HM-32 [51,54,55]. Since we are drawing the linkage inferences from one large family with a high-density SNP data, we took proper care for accounting the false positive due to high LD, and our results were very consistent. At the same time, we have chosen the SNPs in our initial linkage panel in a way (MAF>1%) that increased the marker heterozygosity, hence, increased the linkage information content that improves the likelihood of detecting a recombinant event.

The 12.2 Mb candidate interval on chromosome 7p22.1–7p21.1 contains 61 known transcripts. These include genes involved in regulation of cell proliferation, growth and extension: β-actin (ACTB [OMIM 102630]), fascin homolog 1, actin-bundling protein (FSCN1 [OMIM 602689]), ras-related C3 botulinum toxin substrate 1 (RAC1 [OMIM 602048]), as well as in gene expression: zinc finger protein 12 (ZNF12 [OMIM 194536]). Further, genes considered in this locus are: islet cell autoantigen 1 (ICA1 [OMIM 147625]) and collagen, type XXVIII, alpha 1 (COL28A1 [OMIM 609996]).

The second interesting locus identified at chromosome 7p12.3–7p11.2 that maps to a 10.9 Mb region comprises 30 known transcripts, including growth factor receptor-bound protein 10 (GRB10 [OMIM 601523]) and epidermal growth factor receptor (EGFR [OMIM 131550]). GRB10 encodes a growth factor receptor-binding protein that interacts with insulin receptors and insulin-like growth-factor receptors [55]. Based on an animal model, it has been established, that GRB10 acts as an inhibitor of intracellular signaling pathways regulating growth and metabolism. Gene knockouts in mice results in disproportionate overgrowth of the embryo and placenta [56]. EGFR and its ligands are cell signaling molecules involved in diverse cellular functions, including cell proliferation, differentiation, motility, and survival, and in tissue development [57]. EGFR is the prototypical tyrosine kinase receptor localized to basal and differentiated epithelia in the cornea and is a key regulator for maintaining a healthy cornea and promoting regrowth of the wounded cornea [58,59]. Furthermore, Domínguez et al. [60] reported various EGFR functions in Drosophila eye development.

The locus on chromosome 12p12.3–12p12.1 maps within 9.5 Mb region that has 32 known transcripts which are involved in cell signaling and proliferation: Rho GDP dissociation inhibitor (GDI) beta (ARHGDIB [OMIM 602843]), RAS-like, estrogen-regulated, growth inhibitor (RERG [OMIM 612664]), and epidermal growth factor receptor pathway substrate 8 (EPS8 [OMIM 600206]), phosphoinositide-3-kinase, class 2, gamma polypeptide (PIK3C2G [OMIM 609001]).

In the present linkage analysis, all HM-32 family members who carried the three disease-related haplotypes were found with high-grade myopia, indicating that more than one locus contributes to the high myopia phenotype in this pedigree. It is also possible that one of these linked loci is a major dominant determinant and that the others are modifier genomic variants. Therefore, we hypothesize that the high-grade myopia phenotype in this family could be due to multifactorial inheritance; however, it is difficult to prove this hypothesis until we identify the pathologic mutations. Further research is needed to understand role of multifactorial inheritance and how high-grade myopia can be prevented and/or treated.

Appendix

Acknowledgments

We thank the members of all families for their collaboration. The authors thank Kristen A. Bailey for assistance with data checking and pedigree analyses. Supported by Ministry of Education and Science, Poland, Grant 2 PO5A 095 29.

References

  1. Llorente L, Barbero S, Cano D, Dorronsoro C, Marcos S. Myopic versus hyperopic eyes: axial length, corneal shape and optical aberrations. J Vis. 2004; 4:288-98. [PMID: 15134476]
  2. Burton TC. The influence of refractive error and lattice degeneration on the incidence of retinal detachment. Trans Am Ophthalmol Soc. 1989; 87:143-55. [PMID: 2562517]
  3. Kempen JH, Mitchell P, Lee KE, Tielsch JM, Broman AT, Taylor HR, Ikram MK, Congdon NG, O'Colmain BJ. The prevalence of refractive errors among adults in the United States, Western Europe, and Australia. Arch Ophthalmol. 2004; 122:495-505. [PMID: 15078666]
  4. Wong TY, Foster PJ, Hee J, Ng TP, Tielsch JM, Chew SJ, Johnson GJ, Seah SK. Prevalence and risk factors for refractive errors in adult Chinese in Singapore. Invest Ophthalmol Vis Sci. 2000; 41:2486-94. [PMID: 10937558]
  5. Hosaka A. Population studies–myopia experience in Japan. Acta Ophthalmol Suppl. 1988; 185:37-40. [PMID: 2853537]
  6. Chandran S. Comparative study of refractive errors in West Malaysia. Br J Ophthalmol. 1972; 56:492-5. [PMID: 5069190]
  7. Wu HM, Seet B, Yap EP, Saw SM, Lim TH, Chia KS. Does education explain ethnic differences in myopia prevalence? A population-based study of young adult males in Singapore. Optom Vis Sci. 2001; 78:234-9. [PMID: 11349931]
  8. Kleinstein RN, Jones LA, Hullett S, Kwon S, Lee RJ, Friedman NE, Manny RE, Mutti DO, Yu JA, Zadnik K. Refractive error and ethnicity in children. Arch Ophthalmol. 2003; 121:1141-7. [PMID: 12912692]
  9. Lin LL, Shih YF, Tsai CB, Chen CJ, Lee LA, Hung PT, Hou PK. Epidemiologic study of ocular refraction among schoolchildren in Taiwan in 1995. Optom Vis Sci. 1999; 76:275-81. [PMID: 10375241]
  10. Seroczynska M, Prost ME, Medrun J, Lukasiak E, Oleksiak E. The causes of childhood blindness and visual impairment in Poland. Klin Oczna. 2001; 103:117-20. [PMID: 11873409]
  11. Mutti DO, Mitchell GL, Moeschberger ML, Jones LA, Zadnik K. Parental myopia, near work, school achievement, and children's refractive error. Invest Ophthalmol Vis Sci. 2002; 43:3633-40. [PMID: 12454029]
  12. Jones LA, Sinnott LT, Mutti DO, Mitchell GL, Moeschberger ML, Zadnik K. Parental history of myopia, sports and outdoor activities, and future myopia. Invest Ophthalmol Vis Sci. 2007; 48:3524-32. [PMID: 17652719]
  13. Cordain L, Eaton SB, Brand MJ, Lindeberg S, Jensen C. An evolutionary analysis of the aetiology and pathogenesis of juvenile-onset myopia. Acta Ophthalmol Scand. 2002; 80:125-35. [PMID: 11952477]
  14. Czepita D, Gosławski W, Mojsa A, Muszyńska-Lachota I. Role of light emitted by incandescent or fluorescent lamps in the development of myopia and astigmatism. Med Sci Monit. 2004; 10:CR168-71. [PMID: 15039648]
  15. Baird PN, Schäche M, Dirani M. The GEnes in Myopia (GEM) study in understanding the aetiology of refractive errors. Prog Retin Eye Res. 2010; 29:520-42. [PMID: 20576483]
  16. Dirani M, Chamberlain M, Garoufalis P, Chen C, Guymer RH, Baird PN. Refractive errors in twin studies. Twin Res Hum Genet. 2006; 9:566-72. [PMID: 16899164]
  17. Dirani M, Chamberlain M, Shekar SN, Islam AF, Garoufalis P, Chen CY, Guymer RH, Baird PN. Heritability of refractive error and ocular biometrics: the Genes in Myopia (GEM) twin study. Invest Ophthalmol Vis Sci. 2006; 47:4756-61. [PMID: 17065484]
  18. Hammond CJ, Snieder H, Gilbert CE, Spector TD. Genes and environment in refractive error: the twin eye study. Invest Ophthalmol Vis Sci. 2001; 42:1232-6. [PMID: 11328732]
  19. Liang CL, Yen E, Su JY, Liu C, Chang TY, Park N, Wu MJ, Lee S, Flynn JT, Juo SH. Impact of family history of high myopia on level and onset of myopia. Invest Ophthalmol Vis Sci. 2004; 45:3446-52. [PMID: 15452048]
  20. Young TL, Ronan SM, Drahozal LA, Wildenberg SC, Alvear AB, Oetting WS, Atwood LD, Wilkin DJ, King RA. Evidence that a locus for familial high myopia maps to chromosome 18p. Am J Hum Genet. 1998; 63:109-19. [PMID: 9634508]
  21. Naiglin L, Clayton J, Gazagne C, Dallongeville F, Malecaze F, Calvas P. Familial high myopia: evidence of an autosomal dominant mode of inheritance and genetic heterogeneity. Ann Genet. 1999; 42:140-6. [PMID: 10526656]
  22. Young TL, Ronan SM, Alvear AB, Wildenberg SC, Oetting WS, Atwood LD, Wilkin DJ, King RA. A second locus for familial high myopia maps to chromosome 12q. Am J Hum Genet. 1998; 63:1419-24. [PMID: 9792869]
  23. Naiglin L, Gazagne C, Dallongeville F, Thalamas C, Idder A, Rascol O, Malecaze F, Calvas P. A genome wide scan for familial high myopia suggests a novel locus on chromosome 7q36. J Med Genet. 2002; 39:118-24. [PMID: 11836361]
  24. Paluru P, Ronan SM, Heon E, Devoto M, Wildenberg SC, Scavello G, Holleschau A, Makitie O, Cole WG, King RA, Young TL. New locus for autosomal dominant high myopia maps to the long arm of chromosome 17. Invest Ophthalmol Vis Sci. 2003; 44:1830-6. [PMID: 12714612]
  25. Paluru PC, Nallasamy S, Devoto M, Rappaport EF, Young TL. Identification of a novel locus on 2q for autosomal dominant high-grade myopia. Invest Ophthalmol Vis Sci. 2005; 46:2300-7. [PMID: 15980214]
  26. Zhang Q, Guo X, Xiao X, Jia X, Li S, Hejtmancik JF. A new locus for autosomal dominant high myopia maps to 4q22-q27 between D4S1578 and D4S1612. Mol Vis. 2005; 11:554-60. [PMID: 16052171]
  27. Edwards M, Lewis WH. Autosomal recessive inheritance of myopia in Hong Kong Chinese infants. Ophthalmic Physiol Opt. 1991; 11:227-31. [PMID: 1766686]
  28. Schwartz M, Haim M, Skarsholm D. X-linked myopia: Bornholm eye disease. Linkage to DNA markers on the distal part of Xq. Clin Genet. 1990; 38:281-6. [PMID: 1980096]
  29. Yang Z, Xiao X, Li S, Zhang Q. Clinical and linkage study on a consanguineous Chinese family with autosomal recessive high myopia. Mol Vis. 2009; 15:312-8. [PMID: 19204786]
  30. Farbrother JE, Kirov G, Owen MJ, Pong-Wong R, Haley CS, Guggenheim JA. Linkage analysis of the genetic loci for high myopia on 18p, 12q, and 17q in 51 U.K. families. Invest Ophthalmol Vis Sci. 2004; 45:2879-85. [PMID: 15326098]
  31. Zhang Q, Guo X, Xiao X, Jia X, Li S, Hejtmancik JF. Novel locus for X linked recessive high myopia maps to Xq23-q25 but outside MYP1. J Med Genet. 2006; 43:e20 [PMID: 16648373]
  32. Stambolian D, Ibay G, Reider L, Dana D, Moy C, Schlifka M, Holmes T, Ciner E, Bailey-Wilson JE. Genomewide linkage scan for myopia susceptibility loci among Ashkenazi Jewish families shows evidence of linkage on chromosome 22q12. Am J Hum Genet. 2004; 75:448-59. [PMID: 15273935]
  33. Hammond CJ, Andrew T, Mak YT, Spector TD. A susceptibility locus for myopia in the normal population is linked to the PAX6 gene region on chromosome 11: a genomewide scan of dizygotic twins. Am J Hum Genet. 2004; 75:294-304. [PMID: 15307048]
  34. Wojciechowski R, Moy C, Ciner E, Ibay G, Reider L, Bailey-Wilson JE, Stambolian D. Genomewide scan in Ashkenazi Jewish families demonstrates evidence of linkage of ocular refraction to a QTL on chromosome 1p36. Hum Genet. 2006; 119:389-99. [PMID: 16501916]
  35. Nallasamy S, Paluru PC, Devoto M, Wasserman NF, Zhou J, Young TL. Genetic linkage study of high-grade myopia in a Hutterite population from South Dakota. Mol Vis. 2007; 13:229-36. [PMID: 17327828]
  36. Lam CY, Tam PO, Fan DS, Fan BJ, Wang DY, Lee CW, Pang CP, Lam DS. A genome-wide scan maps a novel high myopia locus to 5p15. Invest Ophthalmol Vis Sci. 2008; 49:3768-78. [PMID: 18421076]
  37. Ciner E, Wojciechowski R, Ibay G, Bailey-Wilson JE, Stambolian D. Genomewide scan of ocular refraction in African-American families shows significant linkage to chromosome 7p15. Genet Epidemiol. 2008; 32:454-63. [PMID: 18293391]
  38. Metlapally R, Ki CS, Li YJ, Tran-Viet KN, Abbott D, Malecaze F, Calvas P, Mackey DA, Rosenberg T, Paget S, Guggenheim JA, Young TL. Genetic association of insulin-like growth factor-1 polymorphisms with high-grade myopia in an international family cohort. Invest Ophthalmol Vis Sci. 2010; 51:4476-9. [PMID: 20435602]
  39. Nakanishi H, Hayashi H, Yamada R, Yamashiro K, Nakata I, Shimada N, Ohno-Matsui K, Mochizuki M, Ozaki M, Yoshitake S, Kuriyama S, Saito M, Iida T, Matsuo K, Matsuda F, Yoshimura N. Single-nucleotide polymorphisms in the promoter region of matrix metalloproteinase-1, −2, and −3 in Japanese with high myopia. Invest Ophthalmol Vis Sci. 2010; 51:4432-6. [PMID: 20435584]
  40. Hysi PG, Young TL, Mackey DA, Andrew T, Fernandez-Medarde A, Solouki AM, Hewitt AW, Macgregor S, Vingerling JR, Li YJ, Ikram MK, Fai LY, Sham PC, Manyes L, Porteros A, Lopes MC, Carbonaro F, Fahy SJ, Martin NG, van Duijn CM, Spector TD, Rahi JS, Santos E, Klaver CC, Hammond CJ. A genome-wide association study for myopia and refractive error identifies a susceptibility locus at 15q25. Nat Genet. 2010; 42:902-5. [PMID: 20835236]
  41. Solouki AM, Verhoeven VJ, van Duijn CM, Verkerk AJ, Ikram MK, Hysi PG, Despriet DD, van Koolwijk LM, Ho L, Ramdas WD, Czudowska M, Kuijpers RW, Amin N, Struchalin M, Aulchenko YS, van Rij G, Riemslag FC, Young TL, Mackey DA, Spector TD, Gorgels TG, Willemse-Assink JJ, Isaacs A, Kramer R, Swagemakers SM, Bergen AA, van Oosterhout AA, Oostra BA, Rivadeneira F, Uitterlinden AG, Hofman A, de Jong PT, Hammond CJ, Vingerling JR, Klaver CC. A genome-wide association study identifies a susceptibility locus for refractive errors and myopia at 15q14. Nat Genet. 2010; 42:897-901. [PMID: 20835239]
  42. Kruskal WH, Wallis WA, Brown DL. Use of ranks in one-criterion variance analysis. J Am Stat Assoc. 1952; 47:583-621.
  43. O'Connell JR, Weeks DE. PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998; 63:259-66. [PMID: 9634505]
  44. Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin–rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002; 30:97-101. [PMID: 11731797]
  45. Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES. Parametric and nonparametric linkage analysis: a unified multipoint approach. Am J Hum Genet. 1996; 58:1347-63. [PMID: 8651312]
  46. Kong A, Cox NJ. Allele-sharing models: LOD scores and accurate linkage tests. Am J Hum Genet. 1997; 61:1179-88. [PMID: 9345087]
  47. Matise TC, Chen F, Chen W, De La Vega FM, Hansen M, He C, Hyland FC, Kennedy GC, Kong X, Murray SS, Ziegle JS, Stewart WC, Buyske S. A second-generation combined linkage physical map of the human genome. Genome Res. 2007; 17:1783-6. [PMID: 17989245]
  48. Weeks DE, Sobel E, O'Connell JR, Lange K. Computer programs for multilocus haplotyping of general pedigrees. Am J Hum Genet. 1995; 56:1506-7. [PMID: 7762577]
  49. Sobel E, Lange K. Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet. 1996; 58:1323-37. [PMID: 8651310]
  50. Young TL. Molecular genetics of human myopia: an update. Optom Vis Sci. 2009; 86:E8-22. [PMID: 19104467]
  51. Hornbeak DM, Young TL. Myopia genetics: a review of current research and emerging trends. Curr Opin Ophthalmol. 2009; 20:356-62. [PMID: 19587595]
  52. Li YJ, Guggenheim JA, Bulusu A, Metlapally R, Abbott D, Malecaze F, Calvas P, Rosenberg T, Paget S, Creer RC, Kirov G, Owen MJ, Zhao B, White T, Mackey DA, Young TL. An international collaborative family-based whole-genome linkage scan for high-grade myopia. Invest Ophthalmol Vis Sci. 2009; 50:3116-27. [PMID: 19324860]
  53. Paget S, Julia S, Vitezica ZG, Soler V, Malecaze F, Calvas P. Linkage analysis of high myopia susceptibility locus in 26 families. Mol Vis. 2008; 14:2566-74. [PMID: 19122830]
  54. Wojciechowski R, Stambolian D, Ciner E, Ibay G, Holmes TN, Bailey-Wilson JE. Genomewide linkage scans for ocular refraction and meta-analysis of four populations in the Myopia Family Study. Invest Ophthalmol Vis Sci. 2009; 50:2024-32. [PMID: 19151385]
  55. Morrione A. Grb10 adapter protein as regulator of insulin-like growth factor receptor signaling. J Cell Physiol. 2003; 197:307-11. [PMID: 14566960]
  56. Holt LJ, Siddle K. Grb10 and Grb14: enigmatic regulators of insulin action–and more? Biochem J. 2005; 388:393-406. [PMID: 15901248]
  57. Wang K, Yamamoto H, Chin JR, Werb Z, Vu TH. Epidermal growth factor receptor-deficient mice have delayed primary endochondral ossification because of defective osteoclast recruitment. J Biol Chem. 2004; 279:53848-56. [PMID: 15456762]
  58. Zieske JD, Takahashi H, Hutcheon AE, Dalbone AC. Activation of epidermal growth factor receptor during corneal epithelial migration. Invest Ophthalmol Vis Sci. 2000; 41:1346-55. [PMID: 10798649]
  59. Yu FS, Yin J, Xu K, Huang J. Growth factors and corneal epithelial wound healing. Brain Res Bull. 2010; 81:229-35. [PMID: 19733636]
  60. Domínguez M, Wasserman JD, Freeman M. Multiple functions of the EGF receptor in Drosophila eye development. Curr Biol. 1998; 8:1039-48. [PMID: 9768358]