Association of vitamin D receptor gene start codon (Fok1) polymorphism with high myopia.

BACKGROUND: High myopia caused primarily due to abnormal emmetropization and excessive axial ocular elongation is associated with sight-threatening ocular pathology. Muscular dysfunction of ocular ciliary muscles due to altered intracellular calcium levels can result in defective mechanotransduction of the eye and retinal defocus. The vitamin D3 receptor (VDR; a intracellular hormone receptor) is known to mediate calcium homeostasis, influencing the development of myopia.
MATERIALS AND METHODS: In the present study, a total of 206 high myopia, 98 low myopia and 250 control samples were analyzed for VDR gene Fok1 (exon 2 start codon) polymorphism using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique.
RESULTS: High myopia patients revealed decrease in the frequency of ff homozygotes (8.3%) as compared to control group (14.0%), with a corresponding increase in frequency of FF homozygotes (68.9% in high myopia vs. 62.8% in controls). The frequency of f allele carriers (Ff and ff) was increased in females of high myopia (35.6%) and low myopia cases (45.4%). Elevated frequency of f allele was found only in early age at onset cases of high myopia (0.227) and later age at onset (10-20 years) cases of low myopia (0.273) as well as in low myopia cases with parental consanguinity (0.458) (P 0.035; χ(2) = 6.692*).
CONCLUSION: The results suggest that VDR gene might not be playing a direct role in the development of myopia, but might contribute indirectly to the risk conferred by mechanical stress factors or growth/development related factors through its role in calcium homeostasis and regulation of ciliary muscle function.

Introduction

Myopia (near-sightedness), a multifactorial ocular disorder, is a common cause of visual impairment in the world.[1] Low myopia (with refractive error <6 D) is primarily physiological, whereas high-grade or pathological myopia (with refractive error >6 D) is a progressive form that is often associated with vitreo-retinal degenerative changes and carries a significant risk for vision loss.[2] The global incidence of high myopia exhibits population variation, with higher prevalence rates in Asia (as high as 70-90%), particularly in countries like China, Japan, Hong Kong and Taiwan.[34]

Myopia development is considered to be influenced by multiple genes and environmental factors, and possibly their interactions.[56] Ample molecular genetic studies indicated high heritability (60–90%) and genetic heterogeneity for non-syndromic myopia.[67] Many studies have suggested a strong genetic background for high myopia.[89] So far, 24 chromosomal loci have been identified for myopia through linkage and association studies: Xq28 (MYP1), 18p11.31 (MYP2), 12q21-31 (MYP3), 7q36 (MYP4), 17q21-22 (MYP5), 22q37.1 (MYP6), 11p13 (MYP7), 3q26 (MYP8), 4q12 (MYP9), 8p23 (MYP10), 4q22-q27 (MYP11), 2q37.1 (MYP12), Xq23 (MYP13), 1p36 (MYP14), 10q21.2 (MYP15), 5p15.33-p15.2 (MYP16), 7p15 (MYP17), 14q22.1-q24.2 (MYP18), 15q12-13, 21q22.3, 12q24, 4q21, 9q34.11 and 2q37, which provide an array of positional candidate genes, all of which are linked in some way or the other to myopia development.[10] However, attempts to identify the causative gene through mutation analysis or expression studies were not fruitful.

The vitamin D3 receptor (VDR) gene serves as a functional and positional candidate gene (located near MYP-3 locus) for myopia. Myopia is primarily the effect of abnormal emmetropization that may be caused by ciliary muscle dysfunction due to impaired calcium and mineral homeostasis or osteogenesis.[1112] VDR (intracellular hormone receptor) is a ligand inducible transcriptional factor known to mediate these functions.[13] This gene encoding VDR is mapped to chromosome 12cen-q12, comprising at least five promoter regions, eight protein-coding exons and six untranslated exons, which are alternatively spliced.[1415] Many association studies of VDR gene polymorphisms have described varied pleiotropic effects of VDR.[16] VDR gene variants are known to be associated with body mass, bone turnover and calcium absorption in several pathological conditions.[6] Of the several gene polymorphisms identified, Fok1 start codon polymorphism is the most important one that has been widely used as a genetic marker for several polygenic diseases related to calcium metabolism. It is located in the 5’ promoter region and can be considered an independent marker in the VDR gene since there is no linkage disequilibrium (LD) with any of the other VDR polymorphisms Bsm1, Apa1, Taq1, and the LD area surrounding this polymorphism seems to be very small.[17] It defines the presence of T > C transition polymorphism (ATG to ACG) in exon 2 of the VDR gene[18] that results in protein of different lengths[19] due to different translation initiation sites. The VDR gene Fok1 polymorphism, also known as start codon polymorphism, produces two different alleles (designated as f and F) distinguished on the basis of the presence or absence of Fok1 restriction site.[20] The allele f (with restriction site) initiates translation at first ATG site and encodes full-length VDR protein, while the F allele (with no restriction site) initiates translation at second ATG site and encodes shorter length protein that lacks three -NH2 terminal amino acids. The present study is an attempt to understand the association of VDR Fok1 polymorphism and its functional significance in myopia susceptibility.

Materials and Methods

Source and study design

The present case–control study was conducted during the period 2005–2008 after getting approval from the Institutional Ethical Committee. A total of 206 high myopia cases (>6 D) and around 98 low myopia (<6 D) cases (to serve as positive control group) were recruited from Maxivision Eye Hospital, Hyderabad, India, after explaining to them the purpose of the study and obtaining their informed consent. All the patients underwent a comprehensive ophthalmic examination which included cycloplegic refraction with cyclopentolate 1%, keratometry (measures central 3 mm), dilated fundoscopy and pachymetry. The mean age at onset for high myopia and low myopia cases was 24.9 ± 0.03 and 38.6 ± 0.48 years, respectively. The mean best corrected visual acuity was 20/125. Only patients who had myopia in both the eyes were included. Patients with ocular disease such as cataract or glaucoma, or a history of retinopathy or connective tissue disorders associated with myopia, such as Stickler or Marfan syndromes, were not included. Around 250 randomly selected age- and sex-matched normal healthy individuals without history of myopia and other ocular diseases were used as controls to compare with patient group. Detailed information on clinical, epidemiological and ophthalmic variables was recorded for every patient using a specified proforma. Genomic DNA from the blood samples collected from each participant in ethylenediaminetetraacetic acid (EDTA) vacutainers was isolated by using salting-out method[21] and used for polymorphic analysis using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique.

Genotyping of Fok1 polymorphism

PCR amplification (Applied Biosystem-9800 Fast Thermal Cycler) of exon-2 of VDR gene was performed with 50 ng of genomic DNA, in 100 μl PCR reaction mixture containing 10 μl 10΄ PCR buffer [100 mM Tris-HCl (pH 8.8), 500 mM KCl, 15 mM MgCl2, 0.01% w/v gelatin]; 200 μM each of dATP, dGTP, dTTP, dCTP; 5.0 pmol of each primer and 2.5 U Taq polymerase, using the following set of primer sequences: forward F5’-AGCTGGCCCTGGCACTGACTCTGCTCT-3’ and reverse R5’-ATGGAAACACCTTGCTTCTTCTCCCTC-3’ to amplify 265 bp fragment.[22] The PCR cycling conditions included initial denatuaration at 94°C for 10 min, denaturation at 94°C for 45 s, annealing at 60°C for 45 s, and extension at 72°C for 45 s for 35 cycles with final extension at 72°C for 10 min. Amplified PCR products were checked on 2% agarose gel electrophoresis. The 265-bp PCR products derived from amplification of exon 2 constituting Fok1 restriction site were genotyped based on presence (f allele) or absence (F allele) of restriction site. Upon restriction digestion, F allele produces a single band of 265 bp while f allele produces 69 and 196 bp. Genotypes of the PCR products derived from individual subjects were presented as ff, Ff or FF. Five microliter aliquots of the PCR product were incubated at 37°C for 3 hours in a 10-μl reaction mixture containing reaction buffer (50 mmol/l NaCl, 10 mmol/l Tris-HCl (pH 7.5), 10 mmol/l MgCl2, 1 mmol/l 1,4-dithiothreitol) and 4 units of the restriction endonuclease Fok1 (New England Biolabs). Aliquots of the Fok1 digested products were electrophoresed for 30 min at 100 V on 3% agarose gel and genotyped based on the band sizes by comparing with standard 100 bp ladder [Figure 1].

Figure 1

Gel picture of vitamin D3 receptor Fok1 genotypes

Statistical analysis

Statistical analyses were performed with Statistical Package for the Social Sciences (SPSS) 15.0. Allele and genotype frequencies were calculated and Chi-square test was used to test the significance of genotype distribution with myopia development. All the P values were two-sided and the level of significance was taken as P < 0.05.

Results

In the present study, the genotype distribution of Fok1 polymorphism of VDR gene in high myopia patients revealed considerable decrease in the frequency of ff homozygotes (8.3%) as compared to control group (14.0%), with a corresponding increase in frequency of FF homozygotes (68.9% in high myopia vs. 62.8% in controls) without much variation in heterozygote frequencies. When the comparison was made between high and low myopia cases (positive controls), the frequency of Ff and ff genotypes showed considerable decrease in the former compared to latter. Both high myopia cases and control cases deviated from Hardy-Weinberg equilibrium (χ2 = 15.89 for high myopia and χ2 = 38.21 for control group), while low myopia cases did not exhibit deviation from Hardy–Weinberg equilibrium. The allelic distribution revealed reduction in the f allele frequency in high myopia cases (0.196) as compared to that of control group (0.256). However, χ2 values for the distribution were insignificant for high (3.927), low (2.463) and pooled cases (3.657) [Table 1].To understand further the functional significance of this polymorphism, the genotype data were stratified with respect to epidemiological variables such as sex of the proband, age at onset, diet, presence of family history and parental consanguinity, and ophthalmic variables such as axial length (AL), anterior chamber depth (ACD), corneal thickness (CT) and corneal curvature (CC).

Table 1

Genotype distribution of VDR gene (Fok1) polymorphism with respect to the type of myopia

The frequency of f allele carriers (Ff and ff) was found to be increased in females with high myopia (35.6%) and low myopia cases (45.4%) when compared to male probands (26.7% with high myopia, 34.9% of low myopia) which reflected in the increase of f allele frequency in female patients among high myopia (0.227) and low myopia (0.29). In comparison to low myopia females, the frequency of f allele carriers was found to be reduced in high myopia females [Table 2]. The frequency of ff genotype was found to be elevated in both high and low myopia patients with early age at onset (10.6% in high myopia and 14.3% in low) as compared to patients with age at onset >10 years (7.14%: 10/140 in high myopia and 9.52%: 8/84 in low myopia). However, the allelic distribution revealed higher frequency of f allele only in early age at onset cases of high myopia (0.227) and late age at onset (10-20 years) cases of low myopia (0.273) [Table 3]. The frequency of carriers was found to be increased in vegetarian group (40.3% high myopia, 44.5% low myopia) as compared to non-vegetarian group (27.4% high myopia, 38.3% low myopia) irrespective of the type of myopia [Table 4]. When parental consanguinity and familial incidence of myopia were considered, the frequency of ff genotype did not show variation with respect to familial incidence among high and low myopia groups [Table 5], but showed an increase in consanguineous (16.7%) group compared to non-consanguineous (9.3%) group in low myopia. The heterozygote (Ff) genotype frequency was also increased in consanguineous group of low myopia, reflecting a significant increase of f allele frequency (0.458) in low myopia cases with parental consanguinity compared to cases with no consanguinity (0.226). The P value (0.035) was also found to be significant (χ 2 value = 6.692) [Table 6]. However, high myopia cases did not show much varaition in the f allele frequency with respect to consanguinity. The mean distribution of ophthalmic variables with respect to Fok1 polymorphism in high and low myopia indicated elevated mean axial length and reduced ACD in low myopia patients with ff genotype.

Table 2

Fok1 polymorphism and sex of the proband

Table 3

Fok1 polymorphism and age at onset

Table 4

Fok1 polymorphism and diet

Table 5

Fok1 polymorphism and familial incidence

Table 6

Fok1 polymorphism and parental consanguinity

Discussion

Multiple well-coordinated biological functions operate with great precision to effect visual function, and therefore, the disturbances in any of the interrelated function would disturb the normal emmetropization mechanism of the eye, contributing to ocular refraction in myopia. The ciliary muscle is one of the most important extralenticular components of the eye involved in accommodation because its force deforms the crystalline lens to produce clear retinal imagery. It was suggested that a deficit in the sympathetic inputs to the ciliary muscle may be associated with a propensity for myopia development.[2023] Excessive stress on the muscular function of ciliary muscles that can result due to mechanical factors, variations of ocular refractive components, defects in calcium homeostasis, bone mineralization, etc. may lead to defective accommodation and retinal defocus.[2425] Certain interactive genes of ocular signal transduction pathway involved in these functions, such as VDR, appear to trigger myopia progression.

Few studies reported association between VDR mediated effects and myopia development. Recent studies have revealed the functional significance of VDR gene Fok1 polymorphism and its association with lower bone mineral density (BMD) and altered calcium metabolism,[26] but there are no previous studies on the association of VDR gene Fok1 polymorphism with myopia development till date.

The results of our study on Fok1 polymorphism could not indicate any genotype with specific risk for both high and low myopia, suggesting that VDR gene might not have a direct role in the development of myopia. However, observations with respect to sex distribution indicated that f allele might show sex-specific risk to develop myopia in females, especially low myopia. Age at onset is the most important risk factor in the development of high myopia. The mean age at onset among high myopia cases in our data was found to be lower (15.1 ± 0.672 years) as compared to low myopia cases who had mean age at onset of 19.75 ± 10.78 years. The frequency of f allele carriers (Ff and ff) was found to be increased in age at onset group 10-20 years in low myopia, indicating that the presence of f allele might confound the risk conferred by growth related factors involved in low myopia development during adolescence. With regard to diet, our observations suggest that f allele might trigger myopia development in individuals who are on vegetarian diet. The presence of f allele might have an interactive role with certain dietary components, and diet taken by these individuals might have exerted cumulative risk to develop myopia. Fok1 genotype distribution in high myopia cases did not show any association with respect to familial incidence, parental consanguinity and genetically influenced ophthalmic variables including AL and ACD, indicating clearly that VDR gene may not be a direct interacting factor for myopia. Another interesting observation is the significant association of f allele with low myopia cases showing parental consanguinity, and of ff genotype with mean ACD and mean AL in low myopia only. Since the development of low myopia is more related to mechanical stress factors, it is understood that the presence of f allele might contribute further to the risk in myopia causation and progression.

Functional studies suggest that the ff genotype results in a longer protein that is a less potent transcriptional activator.[1827] The Fok1 site is located at N terminal region in exon 2 of the gene, and residues at N terminal of the first DNA binding zinc finger domain of VDR are vital for establishing contact with general transcription factor IIB. The encoded protein of F allele, which is shorter by three aminoacids, can form complex TFIIB more effectively so that elevated transcriptional regulation could be achieved compared to f allele.[28] This impaired or inefficient regulation by f allele might bring in alterations in signaling cascade affected by VDR protein. Further, it was reported that ff genotype individuals had lowered BMD[29] and are susceptible to certain diseases like breast cancer.[30] The relation between BMD and mechanical stress in ocular growth has been well documented.[31]

VDR is known to play significant role in the alteration of intracellular calcium levels which might result in impaired contraction and relaxation of ciliary muscles of the eye and defective neuromuscular transduction during myopia development. Apart from calcium homeostasis, it plays a role in bone formation due to altered mineral homeostasis that might influence ocular growth and scleral fibril architecture. It was suggested that VDR might suppress osteoblast differentiating factor, Runx2, to control osteogenesis.[32] In addition, VDR forms a heterodimer with Retinoid X receptor (RXR) which is involved in retinoid functions of the eye like accommodation, pupil responses, and aqueous humor production.[33] Hence, VDR is expected to have an indirect and functional role in the ocular growth and myopia development.

Conclusion

Although significant association was not observed between Fok1 polymorphism of VDR gene and high myopia, the increased frequency of f allele carriers observed among female probands (both high and low myopia), high myopia cases with early age at onset, low myopia cases with late age at onset and vegetarian group of both high and low myopia suggests that VDR gene might not be playing a direct role in myopia development but might contribute indirectly to the risk conferred by mechanical stress factors or growth/development related factors through its role in calcium homeostasis and regulation of ciliary muscle function. Large-scale studies, functional studies and studies of other VDR polymorphisms are necessary to understand the interaction of VDR and other ocular pathway proteins in the development of myopia.[34]