The Mendelian Randomization (MR) analyses utilized data from a publicly accessible GWAS (https://gwas.mrcieu.ac.uk/), which had already obtained the necessary ethical and informed consent approvals during the initial data collection. Consequently, no further ethical or informed consent was required for the study. The software and the databases used are listed in the Table of Materials.
1. Exposure data
GWAS summary information on relative dietary macronutrient intake of fat, protein, carbohydrates, and sugar was obtained from over 235,000 individuals of European ancestries23. Dietary intake of the previous day (UKB) or habitual (all other cohorts) was captured using comprehensive food item questionnaires. The relative intake of each macronutrient was corrected by the total energy intake, allowing for nonlinear effects. Following stringent quality control procedures, the final aggregated sample size is 235,391 for sugar and 268,922 for fat, protein, and carbohydrate23. The basic characteristics of the participants included in the GWAS of relative dietary macronutrient intake are summarized in Table 1.
Genetic association data for micronutrients were sourced from a comprehensive review of the existing literature and the availability of robust GWASs in European populations. Briefly, GWAS summary statistics were included for seven essential minerals, serum calcium, copper, iron, magnesium, phosphorus, selenium (measured in blood or toenail samples), and zinc24,25,26,27,28,29,30,31, as well as seven vitamins, including 25-hydroxyvitamin D, vitamin A1 (retinol), vitamin B6, vitamin B9 (folic acid), vitamin B12, vitamin C, and vitamin E32,33,34,35,36,37. The basic characteristics of the participants included in the GWAS of micronutrients are summarized in Table 1.
2. Outcome data
The outcome data were obtained from the IEU GWAS Project (https://gwas.mrcieu.ac.uk/), developed at the MRC Integrative Epidemiology Unit at the University of Bristol. Specifically, myopia was examined as the outcome using data from source: ukb-b-6353 (n = 460,536), cataracts from source: ebi-a-GCST90018814 (n = 491,877), glaucoma from source: ebi-a-GCST90013865 (n = 406,927), and early age-related macular degeneration from source: ebi-a-GCST010723 (n = 105,248). The basic characteristics of the participants included in the GWAS of myopia and other major ocular disorders are summarized in Table 2. The definitions of ocular disease phenotypes were based on the original GWAS datasets. Myopia (ukb-b-6353) was derived from the UK Biobank and defined based on self-reported use of glasses or contact lenses mainly for distance viewing (short-sightedness). Cataract and glaucoma outcomes were obtained from large GWAS meta-analyses based on clinically diagnosed cases identified through medical records or hospital episode statistics. Early age-related macular degeneration (AMD) was defined according to standardized diagnostic criteria used in the original GWAS consortium studies. Differences in phenotype definitions across studies may introduce some heterogeneity and should be considered when interpreting the results.
3. Selection of instrumental variables
Single-nucleotide polymorphisms (SNPs) associated with exposures at genome-wide significance (P < 5 × 10⁻8) were selected as instrumental variables. Independent SNPs were selected as instrumental variables with an r2 threshold of < 0.001, based on their distance within 10,000 kb in the 1000G European data, using the “TwoSampleMR” R package. The strength of the selected instruments was assessed using the F statistic38. An F statistic greater than 10 is generally considered indicative of strong instrumental variables39. All analyses were harmonized to ensure alignment of effect alleles between exposure and outcome datasets.
4. Mendelian randomization estimates
To estimate the causal relationship between dietary macronutrients and micronutrients and myopia and other major ocular disorders, the “TwoSampleMR” R package was applied for MR analysis. Genetic variants associated with each of the dietary macronutrients and micronutrients were used as instrumental variables to examine their causal relationship with myopia and other major ocular disorders. MR analysis relies on three key assumptions: (1) the relevance assumption, whereby the genetic variants are strongly associated with the exposure; (2) the independence assumption, whereby the variants are independent of confounders; and (3) the exclusion restriction assumption, whereby the variants influence the outcome only through the exposure of interest. Multiple MR methods were implemented, including inverse variance weighted (IVW), weighted median, weighted mode, and MR-Egger, with IVW serving as the primary approach and the other methods used to assess the robustness of the findings. The IVW method provides the most precise causal estimate when all instrumental variables are valid or when horizontal pleiotropy is balanced across variants. The weighted median estimator can generate consistent estimates even if up to 50% of the genetic instruments are invalid. The MR-Egger regression method allows for the detection of directional pleiotropy by estimating an intercept term that reflects the average pleiotropic effect across genetic variants. Initially, causal estimates were calculated using fixed effects IVW models, and random effects IVW models were employed when significant heterogeneity was detected (p < 0.05). To account for multiple comparisons across multiple exposures and outcomes, false discovery rate (FDR) correction using the Benjamini–Hochberg method was applied. Statistical significance was defined as FDR-adjusted P < 0.05, while P < 0.05 was considered nominally significant. All analyses were performed in R (version 4.1.2). To clarify the analytical framework, a schematic representation of the MR design and assumptions is shown in Figure 1.
5. Heterogeneity and sensitivity analyses
To evaluate the robustness of the MR results, heterogeneity tests were conducted to detect potential inconsistencies among the instrumental variables. Heterogeneity was assessed using the Cochran Q statistic, where a P value of < 0.05 indicated significant heterogeneity40. Sensitivity analyses were further performed to assess potential pleiotropy, where genetic variants may influence the outcome through biological pathways other than the exposure of interest. Directional pleiotropy was evaluated using the MR-Egger regression intercept test41. In addition, the Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method was applied to identify and remove potential pleiotropic outliers (p < 0.05) and to provide corrected causal estimates when necessary41.