Metabolic and lifestyle factors in relation to senile cataract: a Mendelian randomization study

We conducted a Mendelian randomization study to determine the associations of body mass index (BMI), type 2 diabetes (T2D), systolic blood pressure (SBP), coffee and alcohol consumption and smoking initiation with senile cataract. Independent single nucleotide polymorphisms associated with the metabolic and lifestyle factors at the p < 5 × 10–8 were selected as instrument variables. Summary-level data for senile cataract were obtained from the FinnGen consortium (20,157 cases and 154,905 non-cases) and UK Biobank study (6332 cases and 354,862 non-cases). Higher genetically predicted BMI and SBP and genetic predisposition to T2D and smoking initiation were associated with an increased risk of senile cataract. The combined odds ratios were 1.19 (95% confidence interval (CI) 1.09–1.29; p < 0.001) per one standard deviation increase in BMI (~ 4.8 kg/m2), 1.13 (95% CI 1.04–1.23; p = 0.004) per 10 mmHg increase in SBP, 1.06 (95% CI 1.03–1.09; p < 0.001) per one unit increase in log-transformed odds ratio of T2D, and 1.19 (95% CI 1.10–1.29; p < 0.001) per one standard deviation increase in prevalence of smoking initiation. Genetically predicted coffee consumption showed a suggestive association with senile cataract (odds ratio per 50% increase, 1.18, 95% CI 1.00–1.40; p = 0.050). This study suggests causal roles of obesity, T2D, SBP and smoking in senile cataract.

Data sources for senile cataract. Summary-level data (i.e., beta coefficient and corresponding standard error) for the associations of exposure-associated SNPs with senile cataract were obtained from the UK Biobank study 22 and the FinnGen consortium 23 . The senile cataract cases were defined by H25 in International Classification of Disease-10 (ICD-10), 3661 in ICD-9 and 37,402 in ICD-8 in FinnGen. We used data from the 4th release of the FinnGen consortium comprising 20,157 cases and 154,905 non-cases of Finnish descent after the exclusion of individuals with ambiguous gender, high genotype missingness (> 5%), excess heterozygosity (± 4 SDs), and non-Finnish ancestry. Data from the UK Biobank study were extracted from the results of the 2nd wave of genome-wide association analyses in UK Biobank, performed by the Neale lab (http:// www. neale lab. is/ ukbioba nk). A total of 6332 cases (defined by H25 in ICD-10) and 354,862 non-cases was included in the analysis after the removal of individuals of non-European ancestry, closely related individuals (third-degree relatives or closer), individuals with sex chromosome aneuploidies, and individuals who had withdrawn consent from the UK Biobank study.
Statistical analysis. We aligned the SNPs based on allele letter and allele frequency. The F statistic was calculated to assess the strength of genetic instruments and the F statistic for one association > 10 was deemed good and less likely to be biased by weak instrument (Table 1). We used the multiplicative random-effects inverse-variance weighted model as the main statistical method. We examined the genetic correlation on cataract between the UK Biobank and FinnGen using the ldsc software 24,25 and found a high consistency (rg = 0.92; p = 1.3 × 10 -9 ). Thus, estimates from two data sources were combined using the fixed-effect meta-analysis method. Four sensitivity analyses, including the weighted median 26 , MR-Egger 27 , MR-PRESSO 28 , and contamination mixture 29 methods, were performed to examine the robustness of the results and detect and correct for possible horizontal pleiotropy. The weighted median approach can provide consistent causal estimates at the presupposition that > 50% of weight comes from the valid instruments 26 . The MR-Egger regression can detect pleiotropic effects by its intercept (p for intercept < 0.5 indicates pleiotropy) and provide estimates after the correction for pleiotropy although it consumes statistical power 27 . The MR-PRESSO method can detect outlying SNPs and provide causal estimates after the removal of identified outliers 28 . The contamination mixture method excels at the analysis using hundreds of SNPs as instrumental variables for an exposure and allows for the presence of invalid SNPs 29 . Considering phenotypic and genetic correlations between alcohol consumption and smoking initiation 21 , we performed multivariable MR analysis 30 to minimize mutual pleiotropy. The odds ratio (OR) of senile cataract was scaled to units listed in Table 1 for exposures. We used Cochrane's Q value to assess the heterogeneity of estimates from SNPs in each analysis. Associations with a p value < 0.008 (0.05/6 outcomes) were deemed significant associations, and associations with a p value ≥ 0.008 and ≤ 0.05 were regarded as suggestive associations. All tests were two-sided and performed using the TwoSampleMR 31 , MR-PRESSO 28 and MendelianRandomization 32 packages in the R software (version 4.0.2).

Results
Among metabolic traits, higher genetically predicted BMI was associated with an increased risk of senile cataract in the FinnGen consortium, UK Biobank study and the meta-analysis ( Fig. 1). A genetically predicted one standard deviation increase in BMI corresponded to a combined OR of 1.19 (95% confidence interval (CI) 1.09, 1.29; p < 0.001). Genetic liability to type 2 diabetes and higher genetically predicted SBP were positively associated with senile cataract in FinnGen and UK Biobank, respectively, and these associations remained directionally consistent in the other data source (Fig. 1). The combined ORs of senile cataract were 1.06 (95% CI 1.03, 1.09; p < 0.001) for one unit increase in genetically predicted log-transformed OR of type 2 diabetes and 1.13 (95% CI 1.04, 1.23; p = 0.004) for genetically predicted 10 mmHg increase in SBP. Among lifestyle factors, higher genetically predicted coffee consumption was associated with an increased risk of senile cataract in the FinnGen consortium but not in the UK Biobank (Fig. 1). A suggestive association for genetically predicted coffee consumption remained in the meta-analysis (OR 1.18, 95% CI 1.00, 1.40; p = 0.050). Genetic predisposition to smoking initiation, but not genetically predicted alcohol consumption, was associated with senile cataract in both data sources (Fig. 1). The combined OR was 1.19 (95% CI 1.10, 1.29; p < 0.001) for genetically predicted one standard deviation increase in the prevalence of smoking initiation. The positive association between smoking initiation and senile cataract remained in the multivariable MR analysis with adjustment for alcohol consumption. Likewise, the null finding for alcohol consumption persisted after adjustment for smoking initiation.
Results of sensitivity analyses are presented in Table 2. The associations overall remained directionally consistent in sensitivity analyses albeit with wider CIs in the weighted median and MR-Egger methods. Most associations with the exception for the association for BMI in UK Biobank persisted in the contamination mixture model. We detected moderate heterogeneity in the analysis of type 2 diabetes and smoking initiation. However, we detected no pleiotropic effects in the MR-Egger regression (the p value of intercept > 0.1). Two, one and one outliers were observed in the analysis of BMI in FinnGen, smoking initiation in FinnGen and alcohol consumption in UK Biobank. These associations changed slightly after the removal of outliers.

Discussion
The present MR study found that higher genetically predicted BMI and SBP and genetic predisposition to type 2 diabetes and smoking initiation were associated with an elevated risk of senile cataract in a combined sample of 26,489 cases and 509,767 non-cases. There was a suggestive increased risk of senile cataract with higher genetically predicted coffee consumption, but no evidence of an association with genetically predicted moderate alcohol consumption. Our MR findings concerning the associations of metabolic syndrome components with senile cataract supported previous observational findings. In a Swedish cohort of 35,369 women followed up for a mean of 8.2 years, the metabolic syndrome and its components, abdominal adiposity, diabetes, and hypertension were associated with an increased risk for cataract extraction 5 . These associations were subsequently replicated in another cohort study of 45,049 Swedish men 6 . Meta-analyses based on observational studies also found that an elevated BMI 33 , type 2 diabetes 34 and hypertension 35 appeared to increase the risk of senile cataract. The present study based on genetic data confirmed and strengthened the causality of these associations, which indicated an increased risk of cataract among individuals with these metabolic risk factors. Therefore, maintaining a healthy BMI and blood pressure can now be more strongly recommended as prevention strategies for senile cataract. With regard to diabetes, although we observed a positive association with cataract, the magnitude of this association should be interpreted with caution given that the MR estimate was based on genetic liability to diabetes. Thus, we might have underestimated this association compared to the observational findings 34 . However, this finding still confirmed the causal detrimental role of diabetes in the development of senile cataract, which suggests that reducing the prevalence of diabetes from the population perspective and doing good diabetes management (e.g., maintaining a healthy glycemic profile) from the clinical perspective may benefit the senile cataract prevention.
Epidemiological data on the association between cigarette smoking and senile cataract are consistent and reveal an increased risk of senile cataract in both current and past smokers compared to never smokers 13,36,37 . This MR study confirmed the causal impact of smoking on senile cataract. Studies also showed smoking cessation as an efficient practice to lower the high risk of senile cataract in smokers even though it may take a long time for the increased risk to decline in heavy smokers [36][37][38] . Taken together, reducing the prevalence of smoking initiation appears to be an applicable way to decrease smoking-raised burden of senile cataract once and for all, and for current smokers smoking cessation should be promoted.
Studies on coffee consumption in relation to senile cataract are scarce and conflicting. A zone-level ecological study revealed an inverse correlation between coffee consumption and prevalence of senile cataract 9 . However, no association between coffee consumption and senile cataract was observed in a cross-sectional study 14 . Results of our study showed a positive association between genetically predicted coffee consumption and senile cataract in the Finnish population, whereas this association was not replicated in the UK Biobank study. Given these inconsistent data, more research is warranted to examine the association between coffee consumption and senile cataract.
Alcohol consumption has been inconsistently associated with senile cataract in observational studies. In a cohort of 34,713 Swedish women, the risk of senile cataract started to increase at a daily use of ≥ 1 alcoholic drinks 10 . However, a recent study with data from two large cohorts implied that low-to-moderate alcohol consumption might lower the risk of undergoing cataract surgery 39 , which is in line with a meta-analysis including observational studies published before May 2014 11 . In another meta-analysis based on 7 prospective cohort studies, the pattern of association between alcohol consumption and senile cataract risk was unclear 12 . This MR study did not detect an association between genetically predicted moderate alcohol consumption and senile cataract, which might hint that the discrepancy across above studies might be attributed to residual confounding.
There are several strengths of the present study. The major merit was the MR design, which minimized residual confounding and other biases and thus reinforced the causal reference in the associations of metabolic www.nature.com/scientificreports/ and lifestyle traits with risk of senile cataract. We confined the study population to individuals of European ancestry in all analyses with the exception for that of type 2 diabetes. Associations were examined in two independent data sources.
There are several limitations of the present study. MR studies can be biased by horizontal pleiotropy. However, we found limited evidence indicating the presence of horizontal pleiotropy based on consideration of the intercept in the MR-Egger regression, few outliers identified by the MR-PRESSO method and consistent causal estimates from sensitivity analyses. SNPs for type 2 diabetes were extracted from a genome-wide association study where less than 21% of the populations are non-European participants 18 , which might have introduced population structure bias in the analysis. However, the proportion of non-European individuals was small and genetic principal components were adjusted for in the genome-wide association studies 18 . Thus, our findings are not likely affected by this bias. There was sample overlap in certain analyses in the UK Biobank. The sample overlap might make the model overfitting and bias the causal estimate towards the observational estimate 40 . Nevertheless, this sample overlap was less likely to mislead our results since our genetic instruments were selected from large-scale genome-wide association studies and we interpreted associations based on the results of metaanalysis of FinnGen and UK Biobank where most of the weight came from FinnGen, which included 3 times as many cases as UK Biobank. In addition, the F statistics of all associations with the exception for that of alcohol consumption in FinnGen (no sample overlap) were over 10, which indicated a good strength of used genetic instruments and no weak instrument bias 40 . Given that senile cataract usually onsets at an old age and shares risk factors with other diseases, like cancer and cardiovascular disease, our MR estimates might be risked by survivor bias 41,42 . However, two populations included in this study were generally young and possibly less likely to be influenced. We could not completely rule out the possibility that selected SNPs are not associated with any confounders (i.e., the violation of the independence assumption of MR) and therefore the observed associations might be challenged by confounding. The confinement of study population to individuals of European descent might limit the generalizability of our findings to other populations. This two-sample MR study was based on several genome-wide association studies in different populations. Differences in population structures, such as female to male ratio and age distribution, might bias our findings. In addition, information on subtypes of cataract was unavailable, and thus, whether the effects of metabolic and lifestyle factors differ across cataract types needs investigation 7 . For analyses of alcohol and coffee consumption, non-linear association was unlikely to be examined in the present study based on the summary-level data.
In conclusion, this study suggests that obesity, type 2 diabetes, high SBP and smoking may be causally associated with the risk of developing senile cataract. Whether coffee consumption increases the risk of senile cataract needs verification. Moderate alcohol consumption appears not to be a causal risk factor for senile cataract.

Data availability
Data from UK Biobank (http:// www. neale lab. is/ uk-bioba nk) and FinnGen consortium (https:// r4. finng en. fi/ pheno/ H7_ CATAR ACTSE NILE) are publicly available. All data analyzed in this study are available in OSF data respiratory (https:// osf. io/ by3un/). Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.