Association study of IGF1 polymorphisms with susceptibility to high myopia in a Japanese population

Abstract

Purpose

Polymorphisms in the insulin-like growth factor 1 (IGF1) gene were previously associated with high or extreme myopia in Caucasian and Chinese populations. In the present study, we investigated whether IGF1 polymorphisms are associated with high myopia in a Japanese population.

Methods

A total of 446 Japanese patients with high myopia (≤−9.00 diopters) and 481 Japanese healthy controls (+1.50 diopters to −1.50 diopters) were recruited. We genotyped seven tagging single-nucleotide polymorphisms (SNPs) in IGF1 and assessed allelic and haplotypic diversity in cases and controls.

Results

There were no statistically significant differences in the allele frequencies of IGF1 SNPs and genotypes between cases and controls (P>0.05). However, the A allele of rs5742629 and the G allele of rs12423791 were associated with a moderately increased risk of high myopia (odds ratio [OR] =1.20 and OR =1.21, respectively) with borderline statistical significance (P=0.0502, corrected P (Pc) =0.21 and P=0.064, Pc=0.29, respectively). The haplotype consisting of the A allele of rs5742629 and the G allele of rs12423791 was marginally associated with the risk of high myopia (P=0.041; OR =1.21); this association was not significant after correction (Pc=0.19).

Conclusion

We found that the IGF1 SNPs are not significantly associated with high myopia in our Japanese population. Our results are in contrast to a previous study in which extreme myopia cases had significantly higher frequencies of the G allele of rs5742629 and the C allele of rs12423791 than controls. Therefore, the IGF1 SNPs may not be important factors for susceptibility to high myopia in all populations. Further genetic studies are needed to elucidate the possible contributions of the IGF1 region to the development of high myopia.

Keywords:

• high myopia
• IGF1
• association study
• polymorphism

Introduction

Myopia is a very common refractive error that significantly impacts public health and economics around the world. High myopia is generally defined by an axial length greater than 26 mm or a spherical refractive error less than or equal to −6 diopters (D) that can cause blindness associated with an increased risk of various ocular and systemic diseases, including retinal detachment, submacular hemorrhage, glaucoma, and macular degeneration.^1,2 High myopia imposes enormous economic and social burdens due to its higher prevalence in the populations of the People’s Republic of China, Singapore, Japan, and other countries when compared to the global average.^3–6 Therefore, it is very important to identify the risk factors associated with high myopia and to establish preventive strategies for high myopia.

Although the cause of myopia is unclear, family-based studies and twin studies have shown that genetic factors are involved in the development and progression of myopia.^7–14 Familial linkage studies have attempted to identify candidate genes for myopia, and significant linkages have been reported at 18 loci (MYP1 to MYP18).¹⁵ Many genome-wide association studies for myopia or high myopia have recently been performed, and many candidate genes and loci for myopia or high myopia have been identified.^16–26

The gene encoding insulin-like growth factor 1 (IGF1) is highly conserved between species, and it encodes a polypeptide that plays an important role in regulating cell proliferation, differentiation, and apoptosis.^27–29 Previous animal studies demonstrated that IGF1 also contributes to eye growth and myopia development in several species.^30–32 In addition, recent genetic studies indicated that single-nucleotide polymorphisms (SNPs) in IGF1 were significantly associated with high or extreme myopia in Caucasian and Chinese populations,^33–35 suggesting that IGF1 may play an important role in the development of high myopia through genetic polymorphisms. However, replication studies by Rydzanicz et al³⁶ with a Polish family cohort and by Miyake et al³⁷ with a Japanese population reported a lack of association of IGF1 SNPs with high or extreme myopia. Therefore, further genetic studies are needed to clarify the contribution of IGF1 SNPs to the development of myopia.

The aim of this study was to investigate whether genetic variants in IGF1 are associated with the risk of developing high myopia in our Japanese population.

Methods

Subjects

We recruited 446 unrelated Japanese individuals with high myopia (refractive error ≤−9.00 D in at least one eye) and 481 unrelated (not related to each other or to the patients) healthy Japanese controls (+1.50 D to −1.50 D in both eyes) at Yokohama City University and Okada Eye Clinic in Japan. The controls were sex-matched (47.7% male) to the patients, with an age range of 24 years to 75 years (mean: 39.3 ± 11.0 years) (Table 1). All participants had similar social backgrounds and resided in the same urban area. All participants were diagnosed by comprehensive ophthalmologic tests, including axial length, fundus examination, spherical power, and corneal curvature (autorefractors: NIDEK ARK-730A, ARK-700A TOPCON KP-8100P, and Biometer/Pachymeter AL-2000; Tomey Corporation, Nagoya, Japan). The individuals with high myopia had no known genetic diseases associated with myopia and/or high myopia, including glaucoma, keratoconus, and Marfan syndrome. 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 each participating institute.

Table 1 Characteristics of the study participants

Parameters	Patients	Controls
Individuals, number	446	481
Male, %	44.6	47.7
Age, years (mean ± SD)	37.9 ± 11.9	39.3 ± 11.0
Spherical refractive errors*, diopters
Right eye (mean ± SD)	−11.7 ± 2.24	Na
Left eye (mean ± SD)	−11.7 ± 2.27	Na
Axial length, mm
Right eye (mean ± SD)	28.0 ± 1.16	Na
Left eye (mean ± SD)	28.0 ± 1.22	Na
Corneal refraction, mm
Right eye (mean ± SD)	43.8 ± 1.45	Na
Left eye (mean ± SD)	43.8 ± 1.49	Na

Notes:

* Patients, spherical refractive errors ≤−9.00 diopters in at least one eye; controls, spherical refractive errors within ± 1.50 diopters in both eyes.

Abbreviations: SD, standard deviation; Na, not applicable.

Patients’ ages ranged from 12 years to 76 years (mean: 37.9 ± 11.9 years), and 44.6% of the patients were male. The average spherical refractive errors were −11.7 ± 2.24 D (range: −6.75 D to −22.75 D) in the right eye and −11.7 ± 2.27 D (range: −8.50 D to −22.5 D) in the left eye. The average axial length was 28.0 ± 1.16 mm (range: 26.0–33.1 mm) for the right eye and 28.0 ± 1.22 mm (range: 26.0–34.7 mm) for the left eye. The average corneal refraction was 43.8 ± 1.45 D (range: 39.5–47.8 D) for the right eye and 43.8 ± 1.49 D (range 39.8–53.0 D) for the left eye (Table 1).

DNA extraction and IGF1 genotyping

The QIAamp DNA Blood Maxi Kit (Qiagen GmbH, Hilden, Germany) was used to collect peripheral blood lymphocytes and to extract genomic deoxyribonucleic acid (DNA) from peripheral blood cells. Procedures were performed under standardized conditions to prevent variation in DNA quality.

We evaluated seven tagging SNPs covering IGF1 and its flanking regions (3 kb upstream and 3 kb downstream): rs6214; rs11111262; rs972936; rs5742629; rs12423791; rs2162679; and rs5742612 (Table 2). These SNPs were selected from HapMap Japanese data (minor allele frequency ≥10%; pairwise r² ≥0.8) (http://hapmap.ncbi.nlm.nih.gov/) to capture 100% of common IGF1 variants in HapMap. Of these SNPs, rs6214 and rs12423791 were associated with high myopia in previous studies.^33–35 Genotyping was performed using the TaqMan 5′ exonuclease assay with validated TaqMan primer-probe sets supplied by Applied Biosystems, Inc. (Foster City, CA, USA). Polymerase chain reaction (PCR) was performed using a reaction mixture with a total volume of 10 μL containing 1X TaqMan Universal PCR Master Mix (Applied Biosystems, Inc.), 24 nM of each primer-probe set, and 3 ng of 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’s fluorescence signal was detected using the StepOnePlus Real-Time PCR System (Applied Biosystems, Inc.).

Table 2 Allele frequencies of seven SNPs in IGF1

SNPs	Chr	Position (build 37.1)	Gene location	Alleles (1 >.2)	Risk allele	Risk allele frequency, %		P	Pc	OR (95% CI)
						Cases (n=446)	Controls (n=481)
rs6214	12	102,793,569	3′-UTr	A >G	G	35.9	33.8	0.35		1.10 (0.91–1.33)
rs11111262	12	102,798,177	Intron 3	G >A	G	80.8	80.6	0.88		1.02 (0.81–1.28)
rs972936	12	102,824,921	Intron 2	A >G	A	54.3	50.5	0.11		1.16 (0.97–1.39)
rs5742629	12	102,857,263	Intron 2	A >G	A	60.3	55.8	0.0502	0.21	1.20 (1.00–1.45)
rs12423791	12	102,858,828	Intron 2	G >C	G	73.4	69.5	0.064	0.29	1.21 (0.99–1.48)
rs2162679	12	102,871,259	Intron 1	A >G	A	68.4	66.3	0.34		1.10 (0.90–1.33)
rs5742612	12	102,874,864	5′-UTR	T >C	T	69.7	67.6	0.32		1.11 (0.91–1.35)

Note: P-values less than 0.1 were corrected for multiple hypothesis testing (Pc) by Haploview using 10,000 permutations.

Abbreviations: SNPs, single-nucleotide polymorphisms; IGF1, insulin-like growth factor 1 gene; Chr, chromosome; n, number; Pc, corrected P-value; OR, odds ratio; Cl, confidence interval; UTR, untranslated region.

Statistical analysis

The Hardy-Weinberg equilibrium, allele frequencies, haplotype frequencies, and linkage disequilibrium were assessed using Haploview 4.1 (Broad Institute, Cambridge, MA, USA).³⁸ Differences in allele and haplotype frequencies between cases and controls were assessed with the x² text. The obtained P-values were corrected (Pc) for multiple hypothesis testing using a permutation test (10,000 iterations) in Haploview.

Results

The genotype frequencies of all seven tagging SNPs were in Hardy-Weinberg equilibrium in controls. Figure 1 shows the overall linkage disequilibrium patterns for the seven SNPs in 927 individuals. The seven SNPs showed a similar pattern of linkage disequilibrium between cases and controls (data not shown). Strong linkage disequilibrium was observed throughout the IGF1 region, and three haplotype blocks were suggested by our analysis; rs6214 and rs11111262 were located in block 1 (D’=1.00; r²=0.44), rs5742629 and rs12423791 were in block 2 (D’=0.99; r²=0.54), and rs2162679 and rs5742612 were in block 3 (D’=1.00; r²=0.94).

Figure 1 Linkage disequilibrium plot of seven SNPs in IGF1.

Notes: The D’ value and r² value (in parentheses) corresponding to each SNP pair are expressed as a percentage and are shown within the respective square. Higher D’ values are indicated by a brighter red.
Abbreviations: SNPs, single-nucleotide polymorphisms; IGF1, insulin-like growth factor 1 gene.

Allele frequencies of seven SNPs were determined in cases and controls (Table 2). No statistically significant association was detected for any of the SNPs between cases and controls (P>0.05). However, the A allele of rs5742629 and the G allele of rs12423791 were associated with a moderately increased risk of high myopia (odds ratio [OR] =1.20; 95% confidence interval [CI] =1.00–1.45 and OR =1.21, 95% CI =0.99–1.48, respectively) with borderline statistical significance (P=0.0502, Pc=0.21 and P=0.064, Pc=0.29, respectively).

Table 3 contains the haplotype analysis of rs5742629 and rs12423791 in block 2. The haplotype, AG, created by the risk alleles was marginally associated with the risk of high myopia (P=0.041; OR =1.21; 95% CI =1.01–1.46), but the association was not significant after multiple-testing correction (Pc=0.19). In block 1 and block 3, there were no differences in the frequencies between cases and controls for any of the haplotypes examined in this study (data not shown).

Table 3 Haplotype frequencies of rs5742629 and rs12423791

Haplotypes	Frequency, %		P	Pc	OR (95% CI)
	Cases (n=446)	Controls (n=481)
rs5742629-rs12423791
AG	60.3	55.6	0.041	0.19	1.21 (1.01–1.46)
GC	26.6	30.2	0.079	0.34	0.83 (0.68–1.02)
GG	13.1	13.9	0.61		0.93 (0.71–1.22)

Note: P-values less than 0.1 were corrected for multiple hypothesis testing (Pc) by Haploview using 10,000 permutations.

Abbreviations: n, number; Pc, corrected P-value; OR, odds ratio; CI, confidence interval.

Discussion

The aim of the present study was to assess whether polymorphisms in the IGF1 region affect the development of high myopia in our Japanese population. To this end, we genotyped seven tagging SNPs in IGF1. Here, we report a lack of associations between IGF1 SNPs and high myopia in our Japanese patients.

Metlapally et al³³ reported that rs6214 in IGF1 exhibited a significant association with several types of myopia, including high myopia (≤−5.00 D in at least one eye) in Caucasians. On the other hand, Rydzanicz et al³⁶ did not find any significant association between IGF1 SNPs, (including rs6214) and high myopia (≤−6.00 D in at least one eye and ≤−5.00 D in the second eye) in a Polish family cohort. Zhuang et al³⁴ demonstrated that IGF1 rs12423791, but not rs6214, was significantly associated with extreme myopia (axial length ≥26.00 mm and <−10.00 D in both eyes) in a Chinese population.³⁴ In addition, a haplotype consisting of rs12423791 and rs5742629 was also associated with extreme myopia.³⁴ Mak et al³⁵ reported that a haplotype including rs12423791 was associated with high myopia (≤−8.00 D in both eyes) in another Chinese population. However, the replication study by Miyake et al³⁷ indicated that none of the major tagging SNPs in IGF1 were associated with high (axial length ≥26.00 mm in both eyes) or extreme (≥28.00 mm in both eyes) myopia in a Japanese population.

In the present study of another Japanese population, we found the lack of an association between rs6214 and high myopia (≤−9.00 D in at least one eye), which is in agreement with previous replication studies.^34–37 On the other hand, two SNPs (rs5742629 and rs12423791) and one haplotype (rs5742629-rs12423791) showed a borderline but not significant association, and patients harboring the A allele of rs5742629 and the G allele of rs12423791 had a moderately increased risk of high myopia. Our results are in contrast to those of a previous study,³⁴ in which extreme myopia cases had significantly higher frequencies of the G allele of rs5742629 and the C allele of rs12423791 than controls; these alleles constituted a risk factor for extreme myopia cases in this previous study.³⁴ This disparity suggests that the alleles of rs5742629 and rs12423791 are not a risk factor for susceptibility to high/extreme myopia.

At least three possible explanations exist for the differing results observed among these association studies. First, differences in the sample sizes of high or extreme myopia cases (628 from Caucasians,³³ 127 from Poland,³⁶ 302 from the People’s Republic of China,³⁴ another 300 from the People’s Republic of China,³⁵ 1,339 from Japan,³⁷ and another 446 from Japan in this study) may be important, because limited sample sizes can sometimes lead to false-positive or negative results in genetic studies. For example, the powers to detect the SNPs with an OR =1.3–1.5 and a minor allele frequency of 40% in the Polish³⁶ and two Chinese studies^34,35 are 15%–34% and 32%–73%, respectively. With less frequent SNPs with a minor allele frequency of 20%, the Polish³⁶ and two Chinese studies^34,35 showed lower power (11%–25% and 23%–59%, respectively) at the same significance level. Second, the differences in the recruitment criteria employed for high/extreme cases may lead to different results among association studies. However, we did not find any relationship between recruitment criteria and association levels. To assess this hypothesis in detail, it will be necessary to perform a meta-analysis of the present and previous studies. Finally, it is possible that environmental factor(s) are required for the association of IGF1 SNPs with high/extreme myopia. In addition to genetic factors, environmental factors such as near work (reading, studying, and computer use) and a lack of outdoor activity are clearly involved in myopia development, and it has been suggested that genetic and environmental interactions are important in establishing one’s predisposition to myopia.^39,40IGF1 SNPs may not be able to contribute to myopia development alone; an interaction between IGF1 SNPs and particular environmental factor(s) may be important for this development. Since the present and previous studies^33–37 did not assess joint contributions to disease risk from IGF1 SNPs and environmental factors, it remains to be seen whether the association of IGF1 SNPs with high/extreme myopia depends on particular environmental factor(s).

The present and previous studies^33–37 used common IGF1 SNPs. Since common variants have only a limited capacity to tag rare variants, these studies have low power to detect associations with rare variants. Both common and rare variants are believed to contribute to the genetic susceptibility to complex diseases.⁴¹ Although previous studies^33–35 that demonstrated positive results have shown significant associations of SNPs with minor allele frequencies of at least 20% with high myopia, there is a possibility that rare IGF1 variants can be causal variants of high myopia, suggesting the need to assess the role of rare IGF1 variants in high myopia.

Conclusion

In conclusion, we found that the IGF1 SNPs are not significantly associated with high myopia in our Japanese population. Our results do not support the positive findings of previous studies, suggesting that the IGF1 SNPs are not important risk factors for susceptibility to high myopia in all populations. However, as IGF1 variants may still affect the risk of high/extreme myopia, further genetic studies with larger sample sizes, other ethnic populations, rare variants, and consideration of gene–environment interactions are needed to clarify the contribution of IGF1 variants to disease development.

Acknowledgments

This study was supported by JSPS KAKENHI grant number 23590382. We sincerely thank all of the participants for their participation in this study, as well as all of the medical staff involved in sample collection and diagnosis.

Disclosure

The authors report no conflicts of interest in this work.