Study of association of PAX6 polymorphisms with susceptibility to high myopia in a Japanese population
Authors Kanemaki N, Meguro A, Yamane T, Takeuchi M, Okada E, Iijima Y, Mizuki N
Received 26 August 2015
Accepted for publication 23 September 2015
Published 27 October 2015 Volume 2015:9 Pages 2005—2011
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 2
Editor who approved publication: Dr Scott Fraser
Nobuyuki Kanemaki,1 Akira Meguro,2 Takahiro Yamane,2 Masaki Takeuchi,2,3 Eiichi Okada,4 Yasuhito Iijima,5 Nobuhisa Mizuki2
1Veterinary Teaching Hospital, Azabu University, Sagamihara, 2Department of Ophthalmology and Visual Science, Yokohama City University Graduate School of Medicine, Yokohama, Kanagawa, Japan; 3Inflammatory Disease Section, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA; 4Okada Eye Clinic, 5Aoto Eye Clinic, Yokohama, Kanagawa, Japan
Purpose: Many studies have investigated the relationship of paired box 6 (PAX6) gene polymorphisms with the risk of high myopia, but the results across studies remain inconsistent and ambiguous. In the present work, we investigated whether PAX6 polymorphisms are associated with high myopia in a Japanese population.
Methods: A total of 1,585 Japanese patients with high myopia (spherical equivalent [SE] <-9.00 diopters [D]) and 1,011 Japanese healthy controls (SE≥-1.00 D) were recruited. To compare genotype frequencies between cases and controls, we genotyped five single nucleotide polymorphisms in the PAX6 gene that are reportedly associated with high/extreme myopia: rs662702, rs3026393, rs644242, rs3026390, and rs667773.
Results: For rs662702, rs644242, and rs667773, odds ratios (ORs) for their risk alleles tended to increase with the progression of SE and axial length in the additive and recessive models. Of these, rs644242 had the highest OR (2.56) in patients with SE<-15 D in both eyes in the recessive model. On the other hand, for rs3026393 and rs3026390, the ORs for their risk alleles tended to increase according to the progression of SE and axial length in the dominant model. Of the two, rs3026393 had the highest OR (2.32) in patients with SE<-15 D in both eyes in the dominant model. However, no significant associations were identified in this study.
Conclusion: We found that these PAX6 single nucleotide polymorphisms were associated with an increased risk of extreme myopia. Although the results, which are in agreement with some previous studies, did not reach statistical significance, PAX6 single nucleotide polymorphisms may be important risk factors for the development of extreme myopia. Further genetic studies with larger sample sizes and taking into account the degree of myopia are needed to clarify the contribution of PAX6 variants in myopia development.
Keywords: high myopia, PAX6, association study, polymorphism
Myopia is a refractive visual disorder that is considered an important public health problem, especially in Asian populations. High myopia is generally defined by a spherical equivalent (SE) refractive error <−6 diopters (D) or an axial length (AL) >26 mm that can cause blindness associated with an increased risk of various ocular diseases, including retinal detachment, glaucoma, cataracts, and macular degeneration.1 Because of its high prevalence compared to the global average in Asian countries such as the People’s Republic of China, Singapore, and Japan, high myopia imposes enormous economic and social burdens in these regions.2–5
Although the cause of myopia is unclear, family-based studies and twin studies have shown that genetic factors contribute to its development and progression.6–13 Therefore, the identification of genetic factors is important for establishing preventive strategies for myopia. To date, 18 genetic loci for myopia (MYP1 to MYP18) have been reported by familial linkage studies.14 In the past several years, genome-wide association studies have led to the identification of many important genetic loci/genes involved in myopia, refractive error, and elongation of AL.15–18 Recently, environmental factors such as education have been reported to influence the association between genetic factors and myopia.19–21
The paired box 6 (PAX6) gene, located on chromosome 11p13, is a highly conserved family of transcription factors containing the paired and homeobox DNA-binding domains. PAX6 regulates the expression of various molecules, such as cell adhesion molecules, cell–cell signaling molecules, hormones, and structural eye proteins, and the gene plays a critical role in the development of the central nervous system and the eye.22 Mutations in the PAX6 gene result in various ocular diseases including aniridia, foveal hypoplasia, presenile cataract, and aniridia-related keratopathy.23 In 2004, a genome-wide linkage analysis in UK dizygotic twins showed the highest peak at the 11p13 locus and strong linkage for refractive error in the vicinity of PAX6 at this locus.24 Subsequently, the results of some genetic surveys suggested an association of single nucleotide polymorphisms (SNPs) in PAX6 with high and/or extreme myopia in Taiwanese, Chinese, and Japanese populations,25–28 suggesting that PAX6 may contribute to the development of high/extreme myopia through genetic polymorphisms. However, other studies have reported a lack of association of PAX6 SNPs with refractive error, myopia, or high myopia in British, Chinese, and mixed populations of some ethnic groups (including mainly Caucasians).29–32
Further genetic studies thus are needed to clarify the contribution of the PAX6 region to myopia development. The aim of the current study was to investigate whether SNPs in the PAX6 region are associated with high/extreme myopia in our Japanese population.
We recruited 1,585 unrelated Japanese patients with high myopia (SE<−9.00 D in at least one eye) and 1,011 unrelated healthy Japanese controls (SE≥−1.00 D in both eyes) at Yokohama City University, Okada Eye Clinic, and Aoto Eye Clinic in Japan. Participants were unrelated to each other, had similar social backgrounds, and resided in the same urban area. All participants were diagnosed by comprehensive ophthalmologic tests, including AL, fundus examination, spherical power, and corneal curvature (autorefractor; NIDEK [Gamagori, Japan]; ARK-730A, ARK-700A, TOPCON [Tokyo, Japan]; KP-8100P, Biometer/Pachymeter AL-2000; Tomey Corporation, Nagoya, Japan). Those with high myopia had no known genetic diseases associated with myopia and/or high myopia, including glaucoma, keratoconus, and Marfan syndrome. We excluded individuals younger than 20 years, which may include potential myopia patients, from the control cohort. Informed consent was obtained from all participants. The study methodology adhered to the tenets of the Declaration of Helsinki and was approved by the relevant ethics committees in Yokohama City University, Okada Eye Clinic, and Aoto Eye Clinic.
DNA and PAX6 genotyping
The QIAamp DNA Blood Maxi Kit (Qiagen, Hilden, Germany) was used to collect peripheral blood lymphocytes and extract genomic DNA from peripheral blood cells. Procedures were performed under standardized conditions to prevent variation in DNA quality.
We evaluated five SNPs in the PAX6 gene regions: rs662702, rs3026393, rs644242, rs3026390, and rs667773. All five have previously been reported to be associated with high or extreme myopia.25–28 Genotyping of all SNPs was performed using the TaqMan 5′ exonuclease assay with primer–probe sets supplied by Applied Biosystems (Foster City, CA, USA). Polymerase chain reaction (PCR) was performed using a reaction mixture with a total volume of 10 μL containing 1× TaqMan Universal PCR Master Mix (Applied Biosystems), 24 nM of each primer–probe set, and 3 ng genomic DNA. The PCR conditions were as follows: 95°C for 10 minutes, followed by 40 cycles of denaturation at 92°C for 15 seconds and annealing/extension at 60°C for 1 minute. The probe fluorescence signal was detected using the StepOnePlus Real-Time PCR System (Applied Biosystems) following the manufacturer’s instructions.
Genetic data were analyzed using SNP and Variation Suite v8 (Golden Helix, Inc., Bozeman, MT, USA, http://www.goldenhelix.com). Multiple inheritance models were applied in analysis of genotype data to assess each risk allele; additive (R/R vs R/nR vs nR/nR), dominant (R/R + R/nR vs nR/nR), and recessive (R/R vs R/nR + nR/nR) models (assuming that R is the risk allele and nR is the non-risk allele) were assessed. All P-values and odds ratios (ORs) were adjusted for age and sex. Linkage disequilibrium (LD) was assessed using Haploview 4.1 (Broad Institute, Cambridge, MA, USA).33
The clinical characteristics of the study populations are shown in Table 1. Patient age ranged from 12 to 78 years (mean 38.8±12.2 years), and 44.2% of patients were male. The average SEs were −10.94±2.06 D (range −4.50 to −22.75 D) in the right eye and −10.82±2.07 D (range −4.50 to −24.50 D) in the left eye. The average AL was 27.55±1.22 mm (range 23.92–33.85 mm) for the right eye and 27.51±1.24 mm (range 23.99–34.74 mm) for the left eye. For the controls, the age ranged from 20 to 87 years (mean 58.2±12.3 years), and 41.6% of them were male. The average SEs were 0.48±0.64 D (range −1.00 to 3.50 D) in the right eye and 0.48±0.63 D (range −1.00 to 3.00 D) in the left eye. Average AL values were 23.22±0.80 mm (range 18.76–25.96 mm) and 23.20±0.79 mm (range 18.99–26.05 mm) for the right and left eyes, respectively.
The genotype frequencies of all five SNPs in PAX6 were in Hardy–Weinberg equilibrium in controls (P>0.05). The magnitude of LD among three SNPs (rs662702, rs644242, and rs667773) was extremely high (r2≥0.98) while the remaining two SNPs (rs3026393 and rs3026390) were in complete LD (r2=1.00). These two LD blocks were not linked to each other (r2≤0.10).
Table 2 shows genotypic association results for the SNPs after stratification according to SE in high myopia cases; that is, <−9.0 D in at least one eye, <−9.0 D in both eyes, <−11.0 D in both eyes, <−13.0 D in both eyes, and <−15.0 D in both eyes. For rs662702, rs644242, and rs667773, ORs for their risk alleles tended to increase according to the progression of SE in the additive and recessive models, and marginal associations were found in patients with SE<−15.0 D in both eyes for rs644242 and rs667773 (rs644242: P=0.069 and OR=2.51; rs667773: P=0.071 and OR=2.50 for the additive model; rs644242: P=0.093 and OR=2.56; rs667773: P=0.095 and OR=2.55 for the recessive model).On the other hand, for rs3026393 and rs3026390, ORs for their risk alleles tended to increase with the progression of SE in the dominant model, and marginal associations were found in patients with SE<−15.0 D in both eyes (rs3026393: P=0.068 and OR=2.32; rs3026390: P=0.075 and OR=2.28).
Table 3 shows genotypic association results after stratification according to AL in high myopia cases, as follows: >26, >28, >29, and >30 mm in both eyes. For rs662702, rs644242, and rs667773, ORs for their risk alleles tended to increase according to the progression of AL in the additive and recessive models, and the highest ORs were observed in patients with AL>30 mm in both eyes in the recessive model (OR=1.87–1.90). On the other hand, for rs3026393 and rs3026390, ORs for their risk alleles tended to increase according to the progression of AL in the dominant model, and the highest ORs were observed in patients with AL>30 mm in both eyes (OR=1.79 and 1.78, respectively). However, no significant associations were found in this genotypic analysis.
The aim of the present study was to assess whether polymorphisms in the PAX6 gene region affect the development of high/extreme myopia in our Japanese population. Patients with SE<−9.00 D in at least one eye were selected for this study because higher degree of myopia suggests that genetic factors may show stronger involvement. We genotyped five polymorphisms in the region and found that the relationship between PAX6 SNPs and myopia became stronger with the progression of myopia, although it did not reach statistical significance. Our findings suggest that PAX6 SNPs play an important role in the risk of extreme myopia in the Japanese population.
Tsai et al25 found a significantly higher frequency of the CC genotype of rs667773 in extreme myopia (SE<−10.0 D in both eyes) (OR=5.265) but not high myopia (SE<−6.0 D in both eyes) compared to controls in Chinese Taiwanese. Liang et al27 reported that the frequencies of the GG genotype of rs662702 and rs644242, which are in extremely high LD with rs667773, significantly increased with the progression of myopia in Chinese Taiwanese and that the GG genotypes of rs662702 and rs644242 had the highest ORs (2.1 and 2.3, respectively) in extreme myopia (SE≤−11.0 D in both eyes). Additionally, Miyake et al28 reported that rs644242 showed a significant association with high (AL≥26 mm in both eyes) and extreme (AL≥28 mm in both eyes) myopia in a Japanese population and that the G allele was protective for both. In the present study, in a Japanese population, we found that the CC genotype of rs667773 and the GG genotype of rs644242 (recessive model) had the highest ORs (2.55 and 2.56, respectively) in extreme myopia with SE<−15.0 D in both eyes and that the GG genotype of rs662702 (recessive model) had the highest OR (1.90) in extreme myopia with AL≥30 mm in both eyes. These findings are in line with the previous studies. Considering the above, the risk allele(s) of rs667773, rs644242, and/or rs662702 may contribute to the development of extreme myopia under the recessive model.
On the other hand, other studies have not found a relationship between rs667773/rs662702 and myopia in several populations.26,29–32 However, these authors, with the exception of Dai et al,32 did not perform an analysis targeting extreme myopia, and moreover, almost all of them, including Dai et al,32 used limited sample sizes, with the potential to yield false results.26,29–32 Therefore, it is assumed that these previous studies had not been performed under appropriate conditions to assess associations between rs667773/rs644242/rs662702 and extreme myopia.
Han et al26 demonstrated that the T allele of rs3026393 and the A allele of rs3026390 are associated with a significantly increased risk of high myopia (<−6.0 D) in dominant and recessive models in Han Chinese. In the present study, we found that the ORs of the G allele of rs3026393 and the G allele of rs3026390 increased according to the progression of myopia in the dominant model. Our result diverges from that of Han et al,26 in that the disease risk-related allele differs between the present and previous studies for rs3026393 and rs3026390. Moreover, other studies found no relationship between the genotypes and myopia in Han Chinese, Chinese Taiwanese, and British populations.27,29,31,32 None of these previous studies26,27,29,31,32 involved analyses targeting extreme myopia with SE<−13.0 D or <−15.0 D, which had higher ORs (OR>1.50) at rs3026393 and rs3026390 in the present study (the study by Han et al26 used the lowest SE cutoff point [≤−10.0 D] among the previous studies). Because the present study found that ORs for rs3026393 and rs3026390 were proportionate to the progression of myopia in the dominant model and that lower ORs (≤1.23) were observed in myopia patients with SE<−11.0 D, there is a possibility that these earlier studies26,27,29,31,32 failed to detect the association between rs3026393/rs3026390 and extreme myopia.
In conclusion, we found that all of the PAX6 SNPs tested in this study were associated with an increased risk of extreme myopia. However, the results did not reach statistical significance; this is probably due to small sample sizes after stratification by SE or AL. Further genetic studies with larger sample sizes and consideration of degree of myopia are needed to elucidate the possible contributions of the PAX6 region in myopia development.
This study was supported by JSPS KAKENHI grant number 23590382. We thank all the participants for their participation in this study and all the medical staff involved in sample collection and diagnosis.
The authors report no conflicts of interest in this work.
Saw SM, Gazzard G, Shih-Yen EC, Chua WH. Myopia and associated pathological complications. Ophthalmic Physiol Opt. 2005;25(5):381–391.
He M, Zeng J, Liu Y, Xu J, Pokharel GP, Ellwein LB. Refractive error and visual impairment in urban children in Southern China. Invest Ophthalmol Vis Sci. 2004;45(3):793–799.
Qian YS, Chu RY, He JC, et al. Incidence of myopia in high school students with and without red-green color vision deficiency. Invest Ophthalmol Vis Sci. 2009;50(4):1598–1605.
Seet B, Wong TY, Tan DT, et al. Myopia in Singapore: taking a public health approach. Br J Ophthalmol. 2001;85(5):521–526.
Matsumura H, Hirai H. Prevalence of myopia and refractive changes in students from 3 to 17 years of age. Surv Ophthalmol. 1999;44:S109–S115.
Wu MM, Edwards MH. The effect of having myopic parents: an analysis of myopia in three generations. Optom Vis Sci. 1999;76(6):387–392.
Guggenheim JA, Kirov G, Hodson SA. The heritability of high myopia: a reanalysis of Goldschmidt’s data. J Med Genet. 2000;37(3): 227–231.
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(12):3633–3640.
Farbrother JE, Kirov G, Owen MJ, Guggenheim JA. Family aggregation of high myopia: estimation of the sibling recurrence risk ratio. Invest Ophthalmol Vis Sci. 2004;45(9):2873–2878.
Liang CL, Yen E, Su JY, et al. Impact of family history of high myopia on level and onset of myopia. Invest Ophthalmol Vis Sci. 2004;45(10):3446–3452.
Klein AP, Duggal P, Lee KE, Klein R, Bailey-Wilson JE, Klein BE. Support for polygenic influences on ocular refractive error. Invest Ophthalmol Vis Sci. 2005;46(2):442–446.
Sorsby A, Fraser GR. Statistical note on the components of ocular refraction in twins. J Med Genet. 1964;1(1):47–49.
Lyhne N, Sjølie AK, Kyvik KO, Green A. The importance of genes and environment for ocular refraction and its determiners: a population based study among 20–45 year old twins. Br J Ophthalmol. 2001;85(12):1470–1476.
Yu L, Li ZK, Gao JR, Liu JR, Xu CT. Epidemiology, genetics and treatments for myopia. Int J Ophthalmol. 2011;4(6):658–669.
Hysi PG, Wojciechowski R, Rahi JS, Hammond CJ. Genome-wide association studies of refractive error and myopia, lessons learned, and implications for the future. Invest Ophthalmol Vis Sci. 2014;55(5):3344–3351.
Fan Q, Barathi VA, Cheng CY, et al. Genetic variants on chromosome 1q41 influence ocular axial length and high myopia. PLoS Genet. 2012;8(6):e1002753.
Cheng CY, Schache M, Ikram MK, et al. Nine loci for ocular axial length identified through genome-wide association studies, including shared loci with refractive error. Am J Hum Genet. 2013;93(2):264–277.
Miyake M, Yamashiro K, Tabara Y, et al. Identification of myopia-associated WNT7B polymorphisms provides insights into the mechanism underlying the development of myopia. Nat Commun. 2015;6:6689.
Fan Q, Wojciechowski R, Kamran Ikram M, et al. Education influences the association between genetic variants and refractive error: a meta-analysis of five Singapore studies. Hum Mol Genet. 2014;23(2):546–554.
Verhoeven VJ, Buitendijk GH, Consortium for Refractive Error and Myopia (CREAM), et al. Education influences the role of genetics in myopia. Eur J Epidemiol. 2013;28(12):973–980.
Goldschmidt E, Jacobsen N. Genetic and environmental effects on myopia development and progression. Eye (Lond). 2014;28(2):126–133.
Simpson TI, Price DJ. Pax6; a pleiotropic player in development. Bioessays. 2002;24(11):1041–1051.
Tsonis PA, Fuentes EJ. Focus on molecules: Pax-6, the eye master. Exp Eye Res. 2006;83(2):233–234.
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(2):294–304.
Tsai YY, Chiang CC, Lin HJ, Lin JM, Wan L, Tsai FJ. A PAX6 gene polymorphism is associated with genetic predisposition to extreme myopia. Eye (Lond). 2008;22(4):576–581.
Han W, Leung KH, Fung WY, et al. Association of PAX6 polymorphisms with high myopia in Han Chinese nuclear families. Invest Ophthalmol Vis Sci. 2009;50(1):47–56.
Liang CL, Hsi E, Chen KC, Pan YR, Wang YS, Juo SH. A functional polymorphism at 3′UTR of the PAX6 gene may confer risk for extreme myopia in the Chinese. Invest Ophthalmol Vis Sci. 2011;52(6):3500–3505.
Miyake M, Yamashiro K, Nakanishi H, et al. Association of paired box 6 with high myopia in Japanese. Mol Vis. 2012;18:2726–2735.
Simpson CL, Hysi P, Bhattacharya SS, et al. The roles of PAX6 and SOX2 in myopia: lessons from the 1958 British Birth Cohort. Invest Ophthalmol Vis Sci. 2007;48(10):4421–4425.
Mutti DO, Cooper ME, O’Brien S, et al. Candidate gene and locus analysis of myopia. Mol Vis. 2007;13:1012–1019.
Jiang B, Yap MK, Leung KH, et al. PAX6 haplotypes are associated with high myopia in Han Chinese. PLoS One. 2011;6(5):e19587.
Dai L, Li Y, Du CY, et al. Ten SNPs of PAX6, Lumican, and MYOC genes are not associated with high myopia in Han Chinese. Ophthalmic Genet. 2012;33(3):171–178.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–265.
This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.Download Article [PDF]