Low serum vitamin D is associated with axial length and risk of myopia in young children

The aim of the study was to investigate the relationship between serum 25(OH)D levels and axial length (AL) and myopia in 6-year-old children. A total of 2666 children aged 6 years participating in the birth-cohort study Generation R underwent a stepwise eye examination. First, presenting visual acuity (VA) and AL were performed. Second, automated cycloplegic refraction was measured if LogMAR VA > 0.1. Serum 25-hydroxyvitamin D [25(OH)D] was determined from blood using liquid chromatography/tandem mass spectrometry. Vitamin D related SNPs were determined with a SNP array; outdoor exposure was assessed by questionnaire. The relationships between 25(OH)D and AL or myopia were investigated using linear and logistic regression analysis. Average 25(OH)D concentration was 68.8 nmol/L (SD ± 27.5; range 4–211); average AL 22.35 mm (SD ± 0.7; range 19.2–25.3); and prevalence of myopia 2.3 % (n = 62). After adjustment for covariates, 25(OH)D concentration (per 25 nmol/L) was inversely associated with AL (β −0.043; P < 0.01), and after additional adjusting for time spent outdoors (β −0.038; P < 0.01). Associations were not different between European and non-European children (β −0.037 and β −0.039 respectively). Risk of myopia (per 25 nmol/L) was OR 0.65 (95 % CI 0.46–0.92). None of the 25(OH)D related SNPs showed an association with AL or myopia. Lower 25(OH)D concentration in serum was associated with longer AL and a higher risk of myopia in these young children. This effect appeared independent of outdoor exposure and may suggest a more direct role for 25(OH)D in myopia pathogenesis. Electronic supplementary material The online version of this article (doi:10.1007/s10654-016-0128-8) contains supplementary material, which is available to authorized users.


Introduction
In the last decades, the prevalence of myopia has increased dramatically in Asia as well as in the Western world [1][2][3]. Prevalence estimates are now around 2 % in 6-year-old children with European ethnicity, and 12 % in children of Asian descent [4,5]. These figures rise to 50 % in young European adults [6] and up to 96 % in students from South Korea [7]. Although myopic refractive error can be corrected optically by glasses, contact lenses, or refractive surgery, the longer axial length ([26 mm) increases the life-time risk of severe visual impairment and blindness due to retinal complications [8]. The basis of myopia is a developmental mismatch between the optical components of the eye [9,10], of which excessive elongation of axial length (AL) in early youth is the most important [11].
The need to reveal the etiology of myopia and develop preventive measures is urgent from a public health perspective. Associations with genetic risk variants [12,13] and environmental factors such as time spent outdoors [14][15][16] and education [4,12] have been well established [17,18]. Recent studies reported an association with serum 25-hydroxy vitamin D [25(OH)D] concentration and myopia in adolescents [19,20]. Whether this reflects the association between outdoor exposure and myopia, or whether vitamin D itself plays a role in the pathophysiology is unclear. Studies investigating the potential relation with vitamin D receptor (VDR) polymorphisms found no consistent relationships [21,22]. Serum 25(OH)D is derived from multiple sources. Cholecalciferol (vitamin D3) is formed in the skin after sunlight exposure, and also absorbed by the gut after dietary intake of e.g., fatty fish. Ergocalciferol (vitamin D2) results from intake of foods containing yeasts and fungi [23,24] Both precursors are hydroxylated in the liver into 25(OH)D. Its active metabolite 1,25(OH) 2 D is formed after transformation in the kidney [25] and is distributed to other sites of the body thereafter. In non-supplemented individuals, sunlight exposure is thought to be the main determinant of 25(OH)D [24,[26][27][28]. The main function of 1,25(OH) 2 D is regulation of calcium and phosphate metabolism in bone tissue and plasma, but it also has metabolic functions in insulin metabolism [29,30]. In neuronal disease such as cognitive decline and Parkinson disease [31,32], it can be involved in immune responses [33] and in DNA transcription and methylation [34,35]. Whether 1,25(OH) 2 D has a direct effect on eye growth is currently unclear.
The aim of this study was to investigate the association between 25(OH)D levels, AL, and the risk of myopia in children at age 6 years in a large population-based study. Additionally, influence of time spent outdoors on these relationships, and vitamin D related genotypes was studied.

Study population
This study was embedded in the Generation R Study, a population-based prospective cohort study of pregnant women and their children in Rotterdam, The Netherlands. The complete methodology has been described elsewhere [36,37]. A total of 4154 children underwent an ophthalmologic examination by trained nurses at the research center at age 6 years and underwent blood withdrawal for serum measurements. The study protocol was approved by the Medical Ethical Committee of the Erasmus Medical Center, Rotterdam (MEC 217.595/2002/20), and written informed consent was obtained from all participants. Research was conducted according to the declaration of Helsinki.

Assessment of AL and myopia
The examination included a stepwise ophthalmological examination.
Step 1 consisted of monocular visual acuity with LogMAR based LEA-charts at 3 meter distance by means of the ETDRS method, and ocular biometry including AL (mm) was measured by Zeiss IOL-master 500 (Carl Zeiss MEDITEC IOL-master, Jena, Germany) per eye; five measurements were averaged to a mean AL [38].
Step 2 was carried out in children with a LogMAR visual acuity of [0.1 in at least one eye and in children wearing prescription glasses, and included performance of automated cycloplegic refraction [Topcon auto refractor KR8900 (Topcon, Japan)] and a complete ophthalmologic work up by an ophthalmologist. Two drops (three in case of dark irises) of cyclopentolate (1 %) were administered at least 30 min before refractive error measurement. Pupil diameter was C6 mm at time of the measurement. Spherical equivalent (SE) was calculated as the sum of the full spherical value and half of the cylindrical value in accordance with standard practice, and myopia was defined as SE B -0.5D in at least one eye. Children with LogMAR visual acuity B0.1, no glasses or ophthalmic history were classified as non-myopic [39,40].

Assessment of 25(OH)D
At a median age of 6.0 y (95 % range 5.6-7.9), nonfasting blood samples were drawn by antecubital venipuncture and stored at -80°C until analysis. Serum samples were collected in all children on the examination day at the research center. The measurements of 25(OH)D (nmol/L) in the samples (110lmL serum per sample) were DEQAS certified and were conducted at the Endocrine Laboratory of the VU University Medical Center, Amsterdam, The Netherlands between July 2013 and January 2014 [41]. Serum 25(OH)D was measured with the use of isotope dilution online solid phase extraction liquid chromatography-tandem mass spectrometry, the 'gold standard' (LC-MS/MS) [42] using a deuterated internal standard [IS: 25(OH)D3-d6] (Synthetica AS, Oslo, Norway). This method is highly sensitive and has been widely used in 25(OH)D studies [43,44]. The limit of quantitation was 4.0 nmol/L; intra-assay CV was \6 %, and interassay CV was \8 % for concentrations between 25 and 180 nmol/L.

Questionnaire
Each mother completed a questionnaire regarding the daily life activities of their child. Time spent playing outdoors and time spent watching television was obtained using questions such as ''how much time does your child spend outdoors/watching television in the morning/afternoon/ evening''. Questions were asked for weekdays and weekend days separately, and answers were multiple choice (never, 0-, -1, 1-2, 2-3, 3-4 h). Total time spent in a week was summed and divided by seven to make an average h/day.

Genotyping of SNPs in vitamin D pathway
Samples were genotyped using Illumina Infinium II HumanHap610 Quad Arrays following standard manufacturer's protocols. Intensity files were analyzed using the Beadstudio Genotyping Module software v.3.2.32, and genotype calling based on default cluster files. Any sample displaying call rates below 97.5 %, excess of autosomal heterozygosity (F \ mean -4SD) and mismatch between called and phenotypic gender were excluded. Genotypes were imputed for all polymorphic SNPs from phased haplotypes in autosomal chromosomes using the 1000 Genomes GIANTv3 panel. SNPs located in genes involved in the Vitamin D metabolic pathway were studied for association with AL and presence of myopia; i.e., genes determining serum 25(OH)D levels (GC, DHCR7, CYP2R1), a gene involved in activation of serum 25(OH)D (CYP27B1), the vitamin D receptor gene (VDR), and the gene involved in deactivation of 1,25-(OH) 2 D in mitochondria (CYP24A1). A total of 33 SNPs [21,45,46] were tested, and analyses were adjusted for multiple testing using Bonferroni adjusted P value 0.05/33, P = 0.0015.

Measurement of covariates
Height and weight of children were measured by trained nurses, and BMI (weight/height 2 ) was calculated. Age was determined at the time of the visit. Income was obtained using the questionnaire and was clustered in low income (lowest tertile) and higher income. If income at the time of the visit was not available, income at birth was used. Ethnicity was obtained in the questionnaire, according to standardized criteria employed by 'Statistics Netherlands', the official national statistics agency [47], concerning the country of birth of parents and child: (1) if both parents were born in the Netherlands, the ethnicity is Dutch; (2) if one of the parents was born in another country than the Netherlands, that country was considered country of birth; (3) if both parents were born in the same country other than the Netherlands, that country was represented; (4) if the parents were born in different countries outside the Netherlands, then the country of the mother was represented; and (5) if that child and both parents were born in different countries outside the Netherlands, the country of birth of the child was represented. Ethnicity was grouped into European and non-European. To adjust for seasonality, four seasons were formed on basis of the month in which the children participated in the study (Winter: December-February, Spring: March-May, Summer: June-August, Autumn: September-November).

Statistical analysis
Separate analyses were performed for AL and myopia. Differences in covariates between myopia and children without myopia were tested using logistic regression analysis adjusting for potentially confounding effects of age and gender. The relation between 25(OH)D and AL was investigated using multivariable linear regression analysis; the relation with myopia (SE B -0.5D) was analyzed using multivariable logistic regression analysis, Covariates were only added to the model if they were significantly related with the outcome as well as with 25(OH)D. Three models were tested: model 1 only adjusted for age and gender; model 2 for age, gender, BMI, ethnicity, television watching, family income, and season visiting the research center; model 3 additionally adjusted for time spent playing outdoors. Effect estimates were determined per 25 nmol/L 25(OH)D. Beta's are presented with SE; Odds Ratios (ORs) with 95 % confidence intervals (95 % CI). Statistical analyses were performed using SPSS version 21.0 for Windows software (SPSS Inc).

Results
Demographics A flow diagram presenting the selection of children for the current analysis is shown in Supplement Figure 1. A total of 2666 children were available for analysis of serum Vitamin D and myopia; 2636 children were available for analysis of serum 25(OH)D and AL. Demographic characteristics are presented in Table 1. Children with myopia were on average somewhat older. Adjusted for age and height, girls had smaller AL than boys but not a lower frequency of myopia. Myopic children had a higher BMI, watched more television, and spent less time outdoors.
Myopia occurred more frequently in children of non-European ethnicity.

Serum 25(OH)D
The average serum 25(OH)D in the total study population was lower than the optimal level of 75 nmol/L [23]. Only 37.2 % (1023) of the children reached this optimal level; these were mostly (41.1 %) children who had been examined in summer time ( Table 2). Figure 1 shows an inverse relation between serum 25(OH)D and AL for the entire population (P \ 0.001). Most myopes had high AL and low serum 25(OH)D levels; only 18 % (11/62) of myopic children reached serum levels which corresponded to the optimal level. Table 3 shows associations between serum 25(OH)D and AL and myopia. Lower serum levels were associated with higher AL and higher risks of myopia. The estimates remained statistically significant after adjustment for covariates. The effect between serum 25(OH)D and AL remained [beta -0.033 (SE 0.012; P 0.02)] after exclusion of myopic children. The association was similar in children of European and non-European descent, but the association with AL in the relatively small non-European group failed to reach statistical significance.

Search for possible explanations
We hypothesized that our findings could be explained by outdoor exposure. Figure 2 shows the positive relation between time spent outdoors and serum 25(OH)D (Pearson, P = \ 0.001). Independent of serum 25(OH)D, time spent outdoors (hr/day) was a risk factor for AL [beta -0.034 (SE 0.012; P 0.003)]. It was not a significant risk factor for myopia (OR 0.81; 95 % CI 0.61-1.07), possibly due to the small number of myopes. The association between serum 25(OH)D and AL and myopia remained significant after adjustment for time spent outdoors (model 3). We explored possible interactions as well, but there was no significant interaction effect between 25(OH)D, ethnicity or income. Additionally, the association was tested separately in the small subgroup with missing data on time spent outdoors. The effect was similar to the effect in the group with data.  To investigate a possible genetic association between Vitamin D and eye growth, we studied genes incorporated in the Vitamin D pathway. We considered single nucleotide polymorphisms (SNPs) in genes that determine serum 25(OH)D levels, in genes involved in activation of serum 25(OH)D, in the vitamin D receptor gene (VDR), and in the gene involved in deactivation of 1,25-(OH) 2 D 3 in mitochondria (CYP24A1) (supplemental Table 1). One SNP (rs2245153) in the CYP24A1 gene showed a significant association with AL (beta 0.039; P 0.04) and myopia (OR 1.55; 95 % CI 1.04-2.31), 2 SNPs in CYP24A1 (rs4809959 beta 0.032; P 0.04 and rs3787557 beta 0.046; P 0.04) and one in the VDR (rs11568820 beta -0.042; P 0.03) only showed a significant association with axial length. P values were all insignificant after adjustment for multiple testing.

Discussion
In this cohort study of young children, we found a significant association between serum 25(OH)D levels, AL and myopia. In this study children with lower serum levels of 25(OH)D had longer AL, and those with higher 25(OH)D had a lower risk of myopia (OR 0.65; 95 % CI 0.46-0.92 per 25 nmol/L). The association remained significant after adjusting for outdoor exposure, indicating that these two closely related determinants may have some overlapping as well as separate effects on the development of myopia. Genetic variants in the vitamin D pathway genes appeared not to be related: although SNPs in the VDR and CYP24A1 genes showed some association with AL and myopia, this did not remain after adjustment for multiple testing.  Our study had strengths and weaknesses. Assets were the particularly large study sample, the inclusion of the combination of measurements of AL and myopia, and the correction for many potential confounders. The young age of our study population was a benefit as well as a potential drawback. It allowed for measurements of the determinant very close to the onset of myopia, leaving less room for confounding bias. On the other hand, it hampered the study of large effects as most children did not develop excessive eye growth yet. There were other drawbacks. We performed cycloplegia only in children with a diminished visual acuity. Reports show that our cut off value of Log-MAR VA of [0.1 had a 97.8 % sensitivity to diagnose myopia [39,40]. We therefore think that our approach did not substantially affect the number of myopes in our study, nor biased the observed associations. Finally, as the correlation between serum 25(OH)D level and time playing outdoors was relatively low in our study, our questionnaire may not have fully assessed all time spent outdoors. Not all participants filled in the questionnaire completely and data on time spent outdoors was partially missing. However, association in the sample of children without data on time spent outdoors was similar to the association in those with complete data.
A novel finding of our study was that the increase in AL in children with low 25(OH)D was already present in the physiological range of refractive error, before the onset of myopia. This implies that Vitamin D has a continuous effect on AL, and not only determines the development of myopia. We confirmed that the risk of myopia decreased with increasing 25(OH)D levels (OR 0.65) with each 25 nmol/L. The association between 25(OH)D and axial length was also significant in the European children; but failed to reach significance in the Non-European group due to low statistical power. Correction for time spent outdoors demonstrated some attenuation of the association, but did not explain it entirely. Whether this is due to residual confounding of time spent outdoors or whether Vitamin D is truly causally related with AL and myopia remains an open question. The evidence for a role of time spent outdoors in myopia is available from cross sectional studies, intervention and randomized clinical trials as well as from animal studies [15,16,48,51]. Vitamin D production is triggered by UV-exposure, not by light exposure per se. Animal studies have shown that artificial light, free of UV, can inhibit development of myopia development [48]. This may suggests that outdoor exposure and Vitamin D are independent risk factors for axial elongation and myopia. However, true causality cannot be concluded from a cross sectional study; longitudinal and functional studies are needed to provide more profound evidence.
A few previous studies have investigated the role of serum 25(OH)D in myopia. A South-Korean and an Australian study found a positive association in adolescents and young adults [19,49]. The ALSPAC study found an association with development of refractive error only for 25(OH)D 2 , not for 25(OH)D 3 in 15 years old children. A potential drawback of this study was the measurement of refraction without any cycloplegia [50]. Mutti et al. [21] found an association between SNPs in the VDR gene and myopia in a smaller study. We could not validate this association, as none of the Vitamin D related SNPs were significant after adjusting for multiple testing.
Various hypotheses underscribe a function of 25(OH)D in eye growth. One theory focusses on Vitamin D in relation to dopamine. The current view is that light exposure initiates the release of dopamine in retinal amacrine cells [51][52][53]. The released dopamine appears to influence the function of gap junctions and the size of receptive fields [54], an important determinant of eye growth. Vitamin D is known to influence dopamine metabolism in neurological disorders, such as Morbus Parkinson and restless legs syndrome [55]. In particular in Parkinson, Vitamin D protects against cell death in the substantia nigra of the dopamine secreting neuron [32,56]. Increased dopamine metabolism [57] was found in the rat brain under influence of vitamin D. In the developing rat brain, Vitamin D was found to upregulate glial derived neurotrophic factor (GDNF) which increases dopamine neurons [58]. Taken together, Vitamin D appears to strengthen the function of dopamine or dopamine secreting cells in neuronal tissues. Whether this also accounts for dopamine secreted by amacrine cells in the retina remains an intriguing question.
Another mechanism may be the regulation of DNA transcription in genes containing vitamin D response elements (VDRE, supplemental figure 2). In this case, the active intracellular 1,25(OH) 2 D binds to VDR binding protein, enters the nucleus, and forms a complex with retinoid X receptor in order to bind to VDRE and initiate transcription. VDREs are located in many genes [59]. It has been shown that retinal cells can metabolize 1,25(OH) 2 D; and this active form of vitamin D may interfere with transcription of genes that promote the myopia signaling cascade [60].
In conclusion, we found that serum levels of 25(OH)D were inversely related to AL, and that low levels increased the risk of myopia. Our data suggest that this relationship may be independent from time spent outdoors. The potential role for 25(OH)D in myopia pathogenesis should be further explored by intervention research and functional studies. Authors' contribution Willem Tideman designed and conducted the research, analyzed the data, wrote the paper and approved the final manuscript as submitted. He had primary responsibility for final content. Jan Roelof Polling designed, conducted the research, analyzed the data and critically revised all versions of the manuscript. He approved the final manuscript as submitted. Trudy Voortman provided comments and consultation regarding the analyses and manuscript and critically revised all versions of the manuscript. She approved the final manuscript as submitted. Vincent Jaddoe initiated and designed the original Generation R study, was responsible for the infrastructure in which the study is conducted, contributed to the original data collection and critically revised the manuscript. He approved the final manuscript as submitted. André Uitterlinden contributed to the analysis, provided comments and consultation regarding the analyses and manuscript. He approved the final manuscript as submitted. Albert Hofman initiated and designed the original Generation R study, was responsible for the infrastructure in which the study is conducted, contributed to the original data collection and critically revised the manuscript. He approved the final manuscript as submitted. Johannes Vingerling provided comments and consultation regarding the analyses and manuscript and critically revised all versions of the manuscript. He approved the final manuscript as submitted. Oscar Franco contributed to the analysis, provided comments and consultation regarding the analyses and manuscript. He approved the final manuscript as submitted. Caroline Klaver designed and conducted the research and wrote the paper and approved the final manuscript as submitted. She had primary responsibility for final content.

Compliance with ethical standards
Conflict of interest The authors have indicated they have no potential conflicts of interest to disclose.
Financial disclosure The authors have no financial relationships relevant to this article to disclose.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.