Molecular Vision 2010; 16:855-861 <http://www.molvis.org/molvis/v16/a95>
Received 20 March 2010 | Accepted 6 May 2010 | Published 13 May 2010

Association analysis of retinoic acid receptor beta (RARβ) gene with high myopia in Chinese subjects

Yang Ding,1,2 Xiaoyan Chen,1,2 Dongsheng Yan,1,2 Anquan Xue,1 Fan Lu,1,2 Jia Qu,1,2 Xiangtian Zhou1,2

The first two authors contributed equally to this work

1School of Optometry and Ophthalmology and Eye Hospital, Wenzhou Medical College. Wenzhou, Zhejiang, China; 2State Key Laboratory Cultivation Base and Key Laboratory of Vision Science, Ministry of Health P.R.China and Zhejiang Provincial Key Laboratory of Ophthalmology and Optometry. Wenzhou, Zhejiang, China

Correspondence to: Xiangtian Zhou, M.D., School of Ophthalmology and Optometry, Wenzhou Medical College, 270 Xueyuan Road, Wenzhou, Zhejiang 325003, P. R. China; Phone: 011-86-577-88824116; FAX: 011-86-577-88824115; email: zxt-dr@wz.zj.cn

Abstract

Purpose: High myopia or pathological myopia is a common refractive error. Individuals with high myopia are subject to increased risk of serious eye complications. Accumulating evidence has demonstrated the role for heritability in ocular growth and in the development of high myopia. Retinoic acid and retinoic acid receptors play important roles in ocular development and in experimentally induced myopia. The purpose of this study was to determine if high myopia is associated with single nucleotide polymorphism (SNP) variants in the retinoic acid receptor beta (RARβ) gene in Chinese subjects.

Methods: DNA samples were purified from venous lymphocytes of 175 unrelated Chinese patients with high myopia (less than −8.00 diopters) and 101 Chinese control subjects without high myopia (±1.00 diopters). Direct nucleotide sequence analysis in the RARβ gene was performed, and the detected variations were further confirmed by reverse sequencing. Allelic frequencies of all detected SNPs were assessed for Hardy–Weinberg equilibrium.

Results: Five variations in RARβ were detected in Chinese subjects with high myopia, including 32574G>A, 32629G>A, 32645C>T, 32647T>G, and 151973C>T, of which only 32647T>G (NCBI notes as rs58244688 and rs2067964) had already been reported. The majority of SNP genotypes were heterozygous. While 32647T>G, 32629G>A, and 32645C>T were located in introns and 32574G>A and 151973C>T were located in coding regions, none of the SNPs affected the amino acid sequence. In the present study, no evidence of association was found between variations in the nucleotide sequence of RARβ and high myopia.

Conclusions: Five SNP variants in RARβ were detected in Chinese subjects with high myopia, none of them were associated significantly with high myopia. Further studies are needed to identify which genes are responsible for high myopia.

Introduction

High myopia is a major cause of visual impairment worldwide, especially in certain young Asian populations such as Chinese, Japanese, and Singaporean [1]. The development of high myopia is associated with excessive axial elongation, scleral thinning, and enlargement of the vitreous chamber depth, all of which are assumed to be the result of scleral remodeling [2,3]. People with high myopia are more likely to have other pathologic changes including primary open angle glaucoma, retinal detachment, macular degeneration, and subretinal hemorrhage [4,5]. High myopia is a complex eye disease involving genetic and environmental factors. While the exact mechanism underlying this abnormal ocular development is still unclear, there is genomic and clinical evidence in different ethnic populations that genetics plays a role in its development [6-10]. Though several mutation genes have been mapped in autosomal dominant or X-recessive inherited high myopia, no specific gene has yet identified [11-17].

All-trans retinoic acid (RA), dopamine, fibroblast growth factor, and transforming growth factor beta are implicated in the control of scleral remodeling and the development of myopia [18-22]. Different from the other candidate etiologies, changes in RA concentrations show sign of defocus selectivity, consistent with the morphological changes of the eye [18]. After visual deprivation in the chick, expression of RA increases in the retina, and the retinoic acid receptor beta (RARβ) is upregulated in the sclera [18,20,23]. Feeding RA to guinea pigs and chicks causes a rapid increase in ocular elongation [24,25]. RA also decreases scleral glycosaminoglycan synthesis while increasing vitreous chamber length during monocular deprivation in marmosets [26]. RA can elicit certain biologic responses by binding to and activating specific retinoic acid receptors (RARs) and the retinoid X receptors (RXRs) [27,28]. In our previous study, we found that six retinoic acid receptor subtypes (RARα, RARβ, RARγ, RXRα, RXRβ, and RXRγ) are present in cultured human scleral fibroblasts, and the proliferation of scleral fibroblasts is inhibited by retinoic acid [29]. We hypothesized that RARs are good candidate genes for myopia, especially high myopia. However a recent study found no association between RARα and either myopia or hypermetropia [30]. Therefore, the purpose of this study was to determine if there are significant associations between high myopia and single nucleotide polymorphism (SNP) variants in the RARβ exons of Chinese subjects.

Methods

Subjects

A total of 175 unrelated Chinese subjects with high myopia, having a spherical power of less than −8.00 diopters (D) in both eyes, were recruited for this study. The age range of the high myopia subjects examined was 7–68 years, and the average age was 34.2 years (Table 1). The average spherical power and the average cylindrical power were mentioned in Table 1. The subjects were from the Optometry Eye Hospital of Wenzhou Medical College, Wenzhou, China and were matched with 101 ethnically and socially similar, unrelated control subjects with refractive errors within ±1.00 D in each eye. None of the subjects had any known ocular, genetic, or systemic connective tissue disorders associated with myopia. The study protocol had the approval of the Human Subjects Ethics Committees of Wenzhou Medical College and Eye Hospital, Wenzhou, China. This study also adhered to the tenets of the Declaration of Helsinki with subsequent revisions. Written consent was obtained from every subject after being fully informed of the purpose and procedures of the study. Every participant received a complete ocular examination including visual acuity (Topcon RM-8800; Topcon Corp., Tokyo, Japan), slit-lamp evaluation of the anterior segment (Topcon SL-1E Slit Lamp; Topcon Corp.), dilated fundus examination (Heine Omega 180 Binocular Indirect Ophthalmoscope; Heine Optotechnik, Herrsching, Germany) and axial-length measurements (Zeiss IOL Master; Carl Zeiss Meditec, Jena, Germany).

DNA extraction and amplification

Genomic DNA for polymerase chain reaction (PCR) was extracted from 2 to 5 ml of peripheral venous blood from all participants. DNA was purified from lymphocytes according to the manufacturer’s instructions using a kit (BBI, Toronto, ON, Canada). Briefly, blood cells were lysed using lysis buffer and digested with proteinase K at 55 °C for 10 min. Then ethanol was added into the mixture and transferred into a spin column. The flow through was discarded after centrifugation. Genomic DNA was eluted from the column. The purity of genomic DNA was determined by the measurement of absorbance ratio A260/280. The accession number for the RARβ sequence used for the construction of our primers was NT_022517. PCR, using primers designed to amplify RARβ (Table 2), was performed in a 50 μl volume with 31.25 μl double distilled H2O, 5 μl 10× PCR buffer, 3.5 μl 25 mmol/l MgCl2, 4 μl 2.5 mmol/l dNTP, 2 μl template DNA, 2 μl primer (Invitrogen, Shanghai, China), and 0.25 μl Ex Taq polymerase (Takara BIO Inc., Tokyo, Japan) following the manufacturer's recommendations (Perkin-Elmer, Norwalk, CT). Briefly, the PCR reaction was initiated at 94 °C for 5 min, followed by 35 cycles (94 °C for 15 s, 55 °C for 30 s, and 72 °C for 45 s) and ended at 72 °C for 10 min.

Single nucleotide polymorphism genotyping by sequencing

After the PCR amplification products were purified using the PCR Cleanup Kit (Axygen Biosciences, Union City, CA), nucleotide sequence analysis was performed with an ABI 3700 sequencer (Applied Biosystems Inc., Foster City, CA). Using Clustal X (version 13.0), we compared all sequencing results to identify variations and to check the sequence map of mutable points to confirm the genotype. Positive results were further confirmed by reverse sequencing. Mutation frequencies of more than 1% were defined as the SNPs.

Analysis of SNPs

Allelic frequencies of all detected SNPs were assessed for Hardy–Weinberg equilibrium. The frequency (P) of allele A was calculated as PA=(2×NAA+NAB)/(2×N) where NAA was the number of wildtypes, NAB was the number of heterozygotes, and N was the sample size of the group. The frequency of allele B was calculated as PB=1-PA. The theoretical value of three genotypes was calculated according to the allele frequency: NAA(E)=N×(PA2),NAB(E)=2N×PA×PB,NBB(E)=N×(PB2) where (E) was the expected value and NBB was the number of homozygous mutations. The Hardy–Weinberg equilibrium for the genotype distributions was examined by the χ2 test in each group, where χ2=[NAA- NAA(E)]2/ NAA(E) + [NAB-NAB(E)]2/ NAB(E) +[NBB- NBB(E)]2/ NBB(E). The level of Hardy–Weinberg equilibrium was set at p>0.05. All samples were assumed to be from the same Mendelian population. Differences in the observed genotype and allelic frequencies between the high myopia and the control subjects were also examined by the χ2 test. Statistical analyses were performed with the SPSS software (version 13.0 for Windows; SPSS Science Inc., Chicago, IL). The power analysis for the χ2 test was run using SAS software (version 9.0; SAS Institute, Cary, NC).

Results

Clinical characteristics

The age range of high myopia onset for most cases (78.7%) was 5–17 years. The ocular characteristics of spherical power, astigmatism and axial length were similar for left and right eyes (Table 1). Of the 175 high myopia cases, almost all exhibited some level of pathological change in the fundus, and some had vision-threatening complications. Leopard-like fundus was the most common sign (45%), followed by conus (22.5%), and chorioretinal atrophy (22.5%). Other changes included complicated cataract (20%), retinal degenerations (8.3%), posterior staphyloma (7.5%), lacquer cracks (6.7%), vitreous opacities (5.8%), retinal hole (5.8%), retinal detachment (4.2%), macular hemorrhage (3.4%), and choroidal neovascularization (0.8%).

Genetic association study

We detected six exons and some flanking regions of the RARβ gene. There were five SNPs: 32574G>A, 32629G>A, 32645C>T, 32647T>G, 151973C>T (Table 3, Figure 1). Of these, only one, 32647T>G, had been reported in the National Center for Biotechnology Information (NCBI) dbSNP database. This intron SNP (NCBI notes as rs58244688 and rs2067964) was a transversion, converting pyrimidine to purine or purine to pyrimidine. The other four were transitions with 32629G>A and 32645C>T located in introns and 32574G>A and 151973C>T located in coding regions (Figure 2). Neither of these SNPs affected the amino acid sequence.

The variation frequency of the SNPs in the high myopia and the control subjects were compared using the χ2 test. There were no significant differences for the genotype and allelic frequencies for the five SNPs between the high myopia and the control subjects (Table 4). Based on the genotype frequencies, the majority of SNP genotypes were heterozygous. All detected SNPs were assessed for Hardy–Weinberg disequilibrium by the χ2 test. For all subjects, the distributions of four of the SNP RARβ genotypes were within the predictions of the Hardy–Weinberg equilibrium (p>0.05, Table 4). However, one of them, 32574G>A in the high myopia group, was not distributed accordingly (p<0.0005).

Discussion

SNPs are the most common type of human genetic variation. They can be identified using various methods, such as direct sequencing, denaturing high performance liquid chromatography, restriction enzyme digestion, and others. Direct genomic DNA sequencing yields high fidelity results with less effort. In this study, six exons and a few introns in RARβ were detected by direct sequencing. Five SNPs were found, of which four (32574G>A, 32629G>A, 32645C>T, and 151973C>T) had not been described previously. Most of the RARβ SNPs were heterozygous, and there were no differences in frequency between high myopia patients and normal subjects. As a result, there was no association between DNA sequence variation of RARβ and high myopia. Two of the SNPs were located in the coding regions of RARβ, and the distribution of one of these in the high myopia group, 32574G>A, was different from that predicted by the Hardy–Weinberg equilibrium. While neither SNP affected the amino acid sequence of the protein, these variations could affect mRNA transcription and/or the secondary structure of the mRNA. Such changes could result in altered expression levels of the RARβ protein [31].

By modulating cell proliferation, differentiation, and cell fate, RA is an important regulator of ocular embryonic and post-embryonic development [27,32]. It is a critical molecule in cultured embryonic photoreceptors as it promotes survival of neurons. It also promotes differentiation of photoreceptors but decreases differentiation of amacrine cells [33,34]. The RA content increases in the retina of chicks within 5 days after initiation of deprivation myopia, and ocular elongation is induced by oral RA in chicks and guinea pigs [23-25,35]. Retinoic acid activates RARβ [36], and it increases the expression of RARβ and downstream genes. RARβ2/RARγ2 double null mutant mice are characterized by retinal dysplasia, abnormal development of the choroid and sclera, and other eye defects [37]. Given the critical role that RA plays in scleral remodeling and the development of myopia, it is likely RA receptors, including RARβ, are important in mediating these changes.

Although RA and RARs are clearly important in the development of myopia in animal models, we did not find any association between the DNA sequence variations of RARβ and high myopia in the current study. Similarly, no variations in the nucleotide sequence of RARα were associated with myopia or hypermetropia [30], although it was differentially expressed during the development of experimental myopia in guinea pigs and chicks [24,38]. In addition, we did not find any genome-wide linkage study on high myopia with a significant log odds score within the chromosomal region containing RARβ. Nevertheless, RARβ variants could alter the responses to RA and affect downstream regulatory networks that govern the development of the ocular refractive system. While our study focused on exons, SNPs in introns that we did not detect may play a role in high myopia. In addition, expanding the sample size may help find some SNPs that influence the development of high myopia by acting as critical points in the complex pathway causing refractive errors. High myopia is a complex eye disease, and the mechanisms by which RARs contribute to myopia are worth further investigation.

Acknowledgments

We thank Britt Bromberg, PhD, Xenofile Editing, for providing editing services for this manuscript. This study was sponsored by Zhejiang Provincial Natural Science Foundation of China (R205739), the Key Research and Development Program project Grant 2008C14097 from Zhejiang Provinceand, the Development Program project grant Y2004A006 from Wenzhou, China and Zhejiang Provincial Program for the Cultivation of High-level Innovative Health talents.

References

  1. Gilmartin B. Myopia: precedents for research in the twenty-first century. Clin Experiment Ophthalmol. 2004; 32:305-24. [PMID: 15180846]
  2. McBrien NA, Cornell LM, Gentle A. Structural and ultrastructural changes to the sclera in a mammalian model of high myopia. Invest Ophthalmol Vis Sci. 2001; 42:2179-87. [PMID: 11527928]
  3. McBrien NA, Lawlor P, Gentle A. Scleral remodeling during the development of and recovery from axial myopia in the tree shrew. Invest Ophthalmol Vis Sci. 2000; 41:3713-9. [PMID: 11053267]
  4. 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]
  5. Chan WM, Ohji M, Lai TY, Liu DT, Tano Y, Lam DS. Choroidal neovascularisation in pathological myopia: an update in management. Br J Ophthalmol. 2005; 89:1522-8. [PMID: 16234465]
  6. Saw SM, Chua WH, Wu HM, Yap E, Chia KS, Stone RA. Myopia: gene-environment interaction. Ann Acad Med Singapore. 2000; 29:290-7. [PMID: 10976381]
  7. Lyhne N, Sjolie 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:1470-6. [PMID: 11734523]
  8. Teikari JM, O'Donnell J, Kaprio J, Koskenvuo M. Impact of heredity in myopia. Hum Hered. 1991; 41:151-6. [PMID: 1937488]
  9. 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]
  10. 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]
  11. Bartsocas CS, Kastrantas A. X-linked form of myopia. Hum Hered. 1981; 31:199-200. [PMID: 7262894]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. Nurnberg G, Jacobi F, Broghammer M, Becker C, Blin N, Nurnberg P, Stephani U, Pusch CM. Refinement of the MYP3 locus on human chromosome 12 in a German family with Mendelian autosomal dominant high-grade myopia by SNP array mapping. Int J Mol Med. 2008; 21:429-38. [PMID: 18360688]
  17. Lam CY, Tam POS, Fan DSP, Fan BJ, Wang DY, Lee CWS, Pang CP, Lam DSA. Genome-wide Scan Maps a Novel High Myopia Locus to 5p15. Invest Ophthalmol Vis Sci. 2008; 49:3768-78. [PMID: 18421076]
  18. Bitzer M, Feldkaemper M, Schaeffel F. Visually induced changes in components of the retinoic acid system in fundal layers of the chick. Exp Eye Res. 2000; 70:97-106. [PMID: 10644425]
  19. Stone RA, Lin T, Laties AM, Iuvone PM. Retinal dopamine and form-deprivation myopia. Proc Natl Acad Sci USA. 1989; 86:704-6. [PMID: 2911600]
  20. Seko Y, Shimokawa H, Tokoro T. In vivo and in vitro association of retinoic acid with form-deprivation myopia in the chick. Exp Eye Res. 1996; 63:443-52. [PMID: 8944551]
  21. Rohrer B, Stell WK. Basic fibroblast growth factor (bFGF) and transforming growth factor beta (TGF-beta) act as stop and go signals to modulate postnatal ocular growth in the chick. Exp Eye Res. 1994; 58:553-61. [PMID: 7925692]
  22. Honda S, Fujii S, Sekiya Y, Yamamoto M. Retinal control on the axial length mediated by transforming growth factor-beta in chick eye. Invest Ophthalmol Vis Sci. 1996; 37:2519-26. [PMID: 8933768]
  23. Seko Y, Shimizu M, Tokoro T. Retinoic acid increases in the retina of the chick with form deprivation myopia. Ophthalmic Res. 1998; 30:361-7. [PMID: 9731117]
  24. McFadden SA, Howlett MH, Mertz JR. Retinoic acid signals the direction of ocular elongation in the guinea pig eye. Vision Res. 2004; 44:643-53. [PMID: 14751549]
  25. McFadden SA, Howlett MH, Mertz JR, Wallman J. Acute effects of dietary retinoic acid on ocular components in the growing chick. Exp Eye Res. 2006; 83:949-61. [PMID: 16797531]
  26. Troilo D, Nickla DL, Mertz JR, Summers Rada JA. Change in the synthesis rates of ocular retinoic acid and scleral glycosaminoglycan during experimentally altered eye growth in marmosets. Invest Ophthalmol Vis Sci. 2006; 47:1768-77. [PMID: 16638980]
  27. Blomhoff R, Blomhoff HK. Overview of retinoid metabolism and function. J Neurobiol. 2006; 66:606-30. [PMID: 16688755]
  28. Napoli JL. Interactions of retinoid binding proteins and enzymes in retinoid metabolism. Biochim Biophys Acta. 1999; 1440:139-62. [PMID: 10521699]
  29. Yan DS, Zhou XT, Chen XY, Lu F, Wang J, Hu DN, Qu J. Expression of retinoid acid receptors in human scleral fibroblasts and regulation of growth of fibroblasts by retinoic acid. Zhonghua Yan Ke Za Zhi. 2007; 43:750-3. [PMID: 18001576]
  30. Veerappan S, Schache M, Pertile K, Islam F, Chen C, Mitchell P, Dirani M, Baird PN. The retinoic acid receptor alpha (RARA) gene is not associated with myopia, hypermetropia, and ocular biometric measures. Mol Vis. 2009; 15:1390-7. [PMID: 19626135]
  31. Shen LX, Basilion JP, Stanton VP, Jr. Single-nucleotide polymorphisms can cause different structural folds of mRNA. Proc Natl Acad Sci USA. 1999; 96:7871-6. [PMID: 10393914]
  32. Luo T, Sakai Y, Wagner E, Drager UC. Retinoids, eye development, and maturation of visual function. J Neurobiol. 2006; 66:677-86. [PMID: 16688765]
  33. Kelley MW, Turner JK, Reh TA. Retinoic acid promotes differentiation of photoreceptors in vitro. Development. 1994; 120:2091-102. [PMID: 7925013]
  34. Stenkamp DL, Gregory JK, Adler R. Retinoid effects in purified cultures of chick embryo retina neurons and photoreceptors. Invest Ophthalmol Vis Sci. 1993; 34:2425-36. [PMID: 8325750]
  35. Niederreither K, McCaffery P, Drager UC, Chambon P, Dolle P. Restricted expression and retinoic acid-induced downregulation of the retinaldehyde dehydrogenase type 2 (RALDH-2) gene during mouse development. Mech Dev. 1997; 62:67-78. [PMID: 9106168]
  36. Balmer JE, Blomhoff R. Gene expression regulation by retinoic acid. J Lipid Res. 2002; 43:1773-808. [PMID: 12401878]
  37. Grondona JM, Kastner P, Gansmuller A, Decimo D, Chambon P, Mark M. Retinal dysplasia and degeneration in RARbeta2/RARgamma2 compound mutant mice. Development. 1996; 122:2173-88. [PMID: 8681798]
  38. Morgan I, Kucharski R, Krongkaew N, Firth SI, Megaw P, Maleszka R. Screening for differential gene expression during the development of form-deprivation myopia in the chicken. Optom Vis Sci. 2004; 81:148-55. [PMID: 15127934]