Xueying Zhu
Dan Cheng
Kaiming Ruan
Meixiao Shen
Yufeng Ye- National Clinical Research Center for Ocular Diseases, Eye Hospital, Wenzhou Medical University, Wenzhou, China
Purpose: The relationship between diabetes mellitus and keratoconus remains controversial. This study aimed to assess the potential causal relationships among type 2 diabetes, glycemic traits, and the risk of keratoconus.
Methods: We used a two-sample Mendelian randomization (MR) design based on genome-wide association summary statistics. Fasting glucose, proinsulin levels, adiponectin, hemoglobin A1c (HbA1c) and type 2 diabetes with and without body mass index (BMI) adjustment were used as exposures and keratoconus was used as the outcome. MR analysis was performed using the inverse-variance weighted method, MR-Egger regression method, weighted-mode method, weighted median method and the MR-pleiotropy residual sum and outlier test (PRESSO).
Results: Results showed that genetically predicted lower fasting glucose were significantly associated with a higher risk of keratoconus [IVW: odds ratio (OR) = 0.382; 95% confidence interval (CI) = 0.261–0.560; p = 8.162 × 10−7]. Genetically predicted lower proinsulin levels were potentially linked to a higher risk of keratoconus (IVW: OR = 0.739; 95% CI = 0.568–0.963; p = 0.025). In addition, genetically predicted type 2 diabetes negatively correlated with keratoconus (IVW: BMI-unadjusted: OR = 0.869; 95% CI = 0.775–0.974, p = 0.016; BMI-adjusted: OR = 0.880, 95% CI = 0.789–0.982, p = 0.022). These associations were further corroborated by the evidence from all sensitivity analyses.
Conclusion: These findings provide genetic evidence that higher fasting glucose levels are associated with a lower risk of keratoconus. However, further studies are required to confirmed this hypothesis and to understand the mechanisms underlying this putative causative relationship.
1. Introduction
Keratoconus is the most common type of ectatic corneal disease. It is defined as progressive corneal thinning causing corneal protrusion, uneven astigmatism, and impaired vision, which can result in legal blindness if left untreated (1–4). The incidence of keratoconus varies from 1:2000 cases as documented in 1986 (5) to 1:375 cases recorded in 2016 (6). The increasing prevalence of keratoconus underscores the significance of identifying the risk factors associated with its development.
Keratoconus is a complex disease that involves a mixture of environmental and genetic factors (7, 8), but its exact etiology remains elusive. Studies on the effects of diabetes mellitus (DM) on the incidence of keratoconus differ significantly. Seiler et al. (9) was the first to demonstrated that DM is a protective factor against keratoconus. Three other studies also found similar results, with a considerably lower rate of DM in patients with keratoconus than in non-keratoconus controls (10–12). However, some conflicting studies have shown that the prevalence of DM is greater in patients with keratoconus than in the control group (13, 14). A recent meta-analysis by Hashemi et al. (15) revealed an unclear correlation between DM and keratoconus. The reasons for these contradictory study results may be due to innate biases or confounders in the observational studies, such as clinic-based case-control recruitment, small sample sizes, reverse causality, and demographic trait heterogeneity.
Mendelian randomization (MR) is a reliable approach for estimating the causal contribution of reported genetic instrumental variables to the disease outcomes of interests (16, 17). MR is less vulnerable to the impact of reverse causality or confounding factors than conventional observational research (17, 18). We used this MR study to determine the causal effect of genetically predicted type 2 diabetes and multiple glycemic traits on the risk of keratoconus in European populations.
2. Methods
2.1. Study design
To estimate the causal relationships between exposure to type 2 diabetes and glycemic traits and the risk of keratoconus, we performed a two-sample MR analysis based on summary statistical data from genome-wide association study (GWAS). We used type 2 diabetes with and without adjustment for body mass index (BMI) and four glycemic traits: fasting glucose, proinsulin levels, adiponectin, and hemoglobin A1c (HbA1c) as exposures. Keratoconus was used as an outcome measure. Applying MR analyses requires the following three essential assumptions (Figure 1): (1) genetic instrumental variables must be closely linked to exposure; (2) genetic instrumental variables are irrelevant to any confounders influencing the exposure-outcome link and (3) genetic instrumental variables affect outcomes merely through their effect on exposure.
Figure 1. Diagram of MR principles investigating the causal relationship between type 2 diabetes, glycemic traits and keratoconus. Instrumental variable assumptions: Assumption 1: genetic instrumental variables must be closely linked to exposure. Assumption 2: genetic instrumental variables are irrelevant to any confounders influencing the exposure-outcome link. Assumption 3: genetic instrumental variables affect outcomes merely through their effect on the exposure.
This study adhered to the principles of the Declaration of Helsinki. This study was approved by the Institutional Review Board of the Eye Hospital of Wenzhou Medical University as it used only publicly available data.
2.2. Genetic instruments for type 2 diabetes and glycemic traits
Intake of diets high in sugar are associated with insulin resistance, hyperglycemia and obesity (19, 20). In the present, insulin resistance (adiponectin levels), hyperglycemia (HbA1c and fasting glucose), and β-cell dysfunction (proinsulin levels) were selected as glycemic traits. GWAS summary datasets of adiponectin, fasting glucose and proinsulin levels were acquired from the MRC-IEU OpenGWAS project software1 (21). These three glycemic traits were derived from different GWASs, including adiponectin (GWASID: ieu-a-1; n = 39,883) (22), fasting glucose (GWASID: ieu-b-114; n = 133,010) (23) and proinsulin levels (GWASID: ebi-a-GCST001212; n = 10,701) (19). The selections of genetic instrumental variables for HbA1c were based on a large GWAS meta-analysis involving 123,665 participants of European ancestry without diabetes (24).2 The summary statistics for type 2 diabetes were obtained from a meta-analysis of GWAS of European ancestry (25). The studies conducted meta-analyses with and without adjustment for BMI. Single nucleotide polymorphisms (SNPs) for type 2 diabetes were obtained from the MRC-IEU OpenGWAS project (21), including type 2 diabetes adjusted for BMI (GWASID: ebi-a-GCST007516; n = 298,957) and unadjusted for BMI (GWASID: ebi-a-GCST007517; n = 298,957) (25). Detailed information of type 2 diabetes and every glycemic trait is presented in Table 1. All SNPs chosen as instrumental variants were strongly associated with the relevant exposure and reached Genome-wide significance (p < 5 × 10−8). The variants were then trimmed using linkage disequilibrium (r ≤ 0.001 within a distance of 10,000 kb for variants at the same locus). The F-statistic was used to quantify the instrumental strength for each candidate SNP and SNPs with F > 10 were retained. Effect estimates of SNP associated with type 2 diabetes were categorized as unadjusted and adjusted for BMI. The final independent SNPs determined as genetic instruments for each exposure were shown in Supplementary Tables S1, S2.
Table 1. Description of GWAS summary statistics for type 2 diabetes and glycemic traits.
2.3. GWAS summary statistics for keratoconus
Genetic variants associated with keratoconus were acquired from the first large scale GWAS. Importantly, we used the first stage meta-analysis comprising of 2,116 cases and 24,626 controls of European ancestry (26).3
2.4. Statistical analysis
For each relevant exposure, we performed two-sample MR analyses using R version 4.1.0 (R Foundation for Statistical Computing, Vienna, Austria). The methods based on the TwoSampleMR version 0.5.6 R package included the inverse-variance weighted (IVW) (27), weighted mode (28), weighted median (29), MR pleiotropy residual sum and outlier (MR-PRESSO) test (30), and MR-Egger regression methods (31). IVW is a major method for assessing the relationships among type 2 diabetes, glycemic traits, and keratoconus (32). Causal analysis of type 2 diabetes and keratoconus were categorized as unadjusted and adjusted for BMI. If the IVW approach indicated an association (p < 8.333 × 10−3 = 0.05/6, taking into account multiple testing for six exposures and one outcome) and the five MR methods had effects in a consistent direction, the results were regarded as statistically significant. p < 0.05, but greater than the significance threshold after Bonferroni correction; the five MR methods had effects in a consistent direction suggesting a potential association. For sensitivity analysis, the Egger intercept calculation (31), MR-PRESSO global test (30), the Cochran’s Q test (33), and the leave-one-out analysis (28) were used to estimate the strength of these identified associations.
3. Results
3.1. MR analysis of fasting glucose and risk of keratoconus
The results showed that a genetically predicted higher fasting glucose level was considerably associated with a lower risk of keratoconus (Table 2): IVW [odds ratio (OR) = 0.382; 95% confidence interval (CI) = 0.261–0.560; p = 8.162 × 10−7]. The uniform direction of the fasting glucose level effect indicates that it has a protective impact against keratoconus. The scatter plots of the MR analysis in Figure 2A demonstrate the effective levels of associations between fasting glucose and keratoconus. Even after Bonferroni adjustment, the IVW and MR PRESSO results remained statistically significant (p < 8.333 × 10−3). The Egger intercept indicated no horizontal pleiotropy effects (p = 0.220 > 0.05). Heterogeneity from the Cochran’s Q test was not statistically significant (p = 0.858 > 0.05). No outliers were detected in the leave-one-out analysis (Figure 2B). Furthermore, no horizontal pleiotropic outliers were identified to distort these results using MR-PRESSO (global test p = 0.862 > 0.05). Collectively, these results support the inverse causal relationship between low fasting glucose levels and keratoconus occurrence.
Table 2. Mendelian randomization estimates for associations between type 2 diabetes, glycemic traits and keratoconus.
Figure 2. MR analysis and leave-one-out analysis of the causal effect of fasting glucose on keratoconus. (A) Scatter plots for MR analyses of the causal effect of fasting glucose on keratoconus. Each line shows the slope corresponding to the estimated MR effect each method. (B) Leave-one-out analysis of the causal effect of fasting glucose on keratoconus. Every black dot represents the IVW MR method applied to estimate the causal effect of fasting glucose on keratoconus, with particular variant excluded from the analysis. The red point represents the IVW estimate using all SNPs.
3.2. MR analysis of proinsulin levels and risk of keratoconus
Genetically predicted proinsulin levels were found to be potentially inversely linked to the incidence of keratoconus (Table 2): IVW (OR = 0.739; 95% CI = 0.568–0.963; p = 0.025). The proinsulin levels effect in consistent direction indicated that it has a potential protective effect against keratoconus. The scatter plots of the MR analysis in Figure 3A demonstrate the effective level of the associations between proinsulin levels and keratoconus. The MR-Egger intercept provided no support for directional pleiotropy (p = 0.693 > 0.05). Heterogeneity from Cochran’s Q test was not statistically significant (p = 0.337 > 0.05). No outliers were detected in the leave-one-out sensitivity analysis (Figure 3B). Furthermore, no horizontal pleiotropic outliers were found in the MR-PRESSO test to distort these results (global test p = 0.335 > 0.05). In combination, these results support a potential negative relationship between low proinsulin levels and the occurrence of keratoconus.
Figure 3. MR analysis and leave-one-out analysis of the causal effect of proinsulin levels on keratoconus. (A) Scatter plots for MR analyses of the causal effect of proinsulin levels on keratoconus. Each line shows the slope corresponding to the estimated MR effect each method. (B) Leave-one-out analysis of the causal effect of proinsulin levels on keratoconus. Every black dot represents the IVW MR method applied to estimate the causal effect of proinsulin levels on keratoconus, with particular variant excluded from the analysis. The red point represents the IVW estimate using all SNPs.
3.3. MR analysis of type 2 diabetes and risk of keratoconus
IVW analyses showed that genetically predicted levels of type 2 diabetes were potentially inversely associated with keratoconus (Table 2): BMI-adjusted: OR = 0.880; 95% CI = 0.789–0.982; p = 0.035; BMI-unadjusted: OR = 0.869; 95% CI = 0.775–0.974; p = 0.016. The uniform direction of the impact of type 2 diabetes indicated that it has a potential protective effect against keratoconus. The scatter plots of the MR analyses in Figure 4A (BMI-unadjusted) and Figure 5A (BMI-adjusted) demonstrate the effective level of the associations between type 2 diabetes and keratoconus. The MR-Egger intercept test did not show any directional pleiotropy (BMI-unadjusted: p = 0.359; BMI-adjusted: p = 0.497). Heterogeneity from Cochran’s Q test was not statistically significant (BMI-unadjusted: p = 0.208; BMI-adjusted: p = 0.339). The leave-one-out analysis identified no outliers (Figure 4B: BMI-unadjusted; Figure 5B: BMI-adjusted). Furthermore, MR-PRESSO analysis did not identify horizontal pleiotropic variants, distorting these results with global test p > 0.05 (BMI-unadjusted: p = 0.225; BMI-adjusted: p = 0.370). All together, these results support type 2 diabetes is potentially inversely associated with the incidence of keratoconus.
Figure 4. MR analysis and leave-one-out analysis of the causal effect of type 2 diabetes unadjusted for BMI on keratoconus. (A) Scatter plots for MR analyses of the causal effect of type 2 diabetes unadjusted for BMI on keratoconus. The slope of each line corresponds to the estimated MR effect per method. (B) Leave-one-out analysis of the causal effect of type 2 diabetes unadjusted for BMI on keratoconus. Each black point represents the IVW MR method applied to estimate the causal effect of type 2 diabetes unadjusted for BMI on keratoconus, excluding that particular variant from the analysis. The red point represents the IVW estimate using all SNPs.
Figure 5. MR analysis and leave-one-out analysis of the causal effect of type 2 diabetes adjusted for BMI on keratoconus. (A) Scatterplots for MR analyses of the causal effect of type 2 diabetes adjusted for BMI on keratoconus. The slope of each line corresponds to the estimated MR effect per method. (B) Leave-one-out analysis of the causal effect of type 2 diabetes adjusted for BMI on keratoconus. Each black point represents the IVW MR method applied to estimate the causal effect of type 2 diabetes adjusted for BMI on keratoconus, excluding that particular variant from the analysis. The red point represents the IVW estimate using all SNPs.
4. Discussion
To the best of our knowledge, this study is the first to describe the causal associations between type 2 diabetes, glycemic traits, and the risk of keratoconus using MR analysis. This study provided evidence of a causal relationship between higher fasting glucose levels and a lower incidence of keratoconus. Meanwhile, our findings imply that increased proinsulin levels potentially decrease the risk of keratoconus. In addition, type 2 diabetes was related to a decreased incidence of keratoconus irrespective of whether type 2 diabetes was adjusted or unadjusted for BMI.
The relationship between DM and keratoconus remains controversial owing to the conflicting results from multiple studies. Several studies have shown an inverse relationship between DM and the risk of keratoconus, suggesting a protective role in the development of keratoconus (9–12). In contrast, other studies that have demonstrated a positive correlation or no association between DM and keratoconus (2, 13, 34, 35). In the present study, we discovered that higher fasting glucose levels were associated with a lower risk of keratoconus, but there was no discernible effect of HbA1c levels. HbA1c represents the average level of blood glucose during the previous two to 3 months (36). Previous studies have reported that genetically predicted fasting glucose and HbA1c levels are inconsistent and there is significant disagreement in the diagnosis of DM (37–39). Our findings indicate that adiponectin does not protect against keratoconus. Adiponectin, an endocrine hormone mainly generated and released by adipocytes, has no causal impact on glucose homeostasis and type 2 diabetes, and correlations among them in observational designs may be caused by the underlying confounding factors reported by Chen et al. (40). The precursor of insulin, proinsulin, is secreted in increased quantities when pancreatic β-cells are stressed, and prior studies have demonstrated that increased proinsulin could serve as a signal to those with prediabetes (41). In particular, our results indicate that proinsulin serves as a protective factor for keratoconus and increases proinsulin levels, thereby reducing the risk of keratoconus. An elevated proinsulin-to-insulin ratio in the blood has been postulated as a potential marker of type 2 diabetes for more than 20 years (42). Notably, our study supported an accordant trend that is elevated fasting glucose, increased proinsulin levels and type 2 diabetes are all negatively related to the incidence of keratoconus.
These interesting findings raise concerns regarding how and why fasting glucose levels affect keratoconus. These mechanisms may include modifications in corneal biomechanics and collagen crosslinking, changes in the extracellular matrix structure, oxidative stress, proteolytic activity, and increased inflammation (43). Here, we postulate on a theoretic level that DM protects keratoconus by altering the biomechanics of the cornea and increasing collagen cross-linking. A few studies have shown that the expression and activity of LOX, a copper amine oxidase that triggers the collagen cross-linking (44), increased in skin collagen (45, 46), rat retinal endothelial cells (46), and ARPE-19 cells (47) under hyperglycemic states. Moreover, the cornea has detected to contain the Lysyl oxidase (LOX) enzyme (48). Therefore, if LOX expression is upregulated in corneal cells under hyperglycemia, this could explain why individuals with DM have a lower risk of developing keratoconus (43). In addition, diabetes may protect against keratoconus through non-enzymatic approaches. Advanced glycation end product (AGE)-mediated crosslinking increases corneal stromal collagen cross-linking, thus strengthening corneal stiffness (44). Nevertheless, the underlying mechanisms of this association need to be assessed more comprehensively.
The use of MR analysis to evaluate large-scale databases using standard procedures is one of the main advantages of this study. The MR method is far less vulnerable to biases or confounders than observational studies. However, this study has some limitations. First, because individuals of European ancestry participated in this study, our findings applying to other races needs further investigation. Second, further laboratory research is required, because the effects of fasting glucose on keratoconus development have not yet been experimentally explored. Third, we only took into account BMI when selecting the instrumental variables for type 2 diabetes in this two-sample MR design. However, no discernible differences were observed between the unadjusted and adjusted BMI models, making multivariate MR a more practical method.
Our findings concluded that the relationship between low fasting levels and a high risk of keratoconus is causal. Furthermore, higher proinsulin levels decrease the risk of keratoconus. These findings, utilized in clinics, indicate that a rational sugar diet strategy may be favorable in controlling keratoconus. However, further studies are required to confirm this hypothesis and to comprehend the mechanisms underlying this putative causative relationship.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
XZ: Writing – original draft, Writing – review & editing. DC: Writing – review & editing. KR: Writing – review & editing. MS: Writing – review & editing. YY: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by research grants from the National Natural Science Foundation of China (No. 81900910) and Science and Technology Innovation Activity Plan of Zhejiang (2023R413086).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2023.1264061/full#supplementary-material
Footnotes
References
1. Godefrooij, DA, de Wit, GA, Uiterwaal, CS, Imhof, SM, and Wisse, RP. Age-specific incidence and prevalence of keratoconus: a nationwide registration study. Am J Ophthalmol. (2017) 175:169–72. doi: 10.1016/j.ajo.2016.12.015
2. Bak-Nielsen, S, Ramlau-Hansen, CH, Ivarsen, A, Plana-Ripoll, O, and Hjortdal, J. A nationwide population-based study of social demographic factors, associated diseases and mortality of keratoconus patients in Denmark from 1977 to 2015. Acta Ophthalmol. (2019) 97:497–504. doi: 10.1111/aos.13961
3. Mas Tur, V, MacGregor, C, Jayaswal, R, O’Brart, D, and Maycock, N. A review of keratoconus: diagnosis, pathophysiology, and genetics. Surv Ophthalmol. (2017) 62:770–83. doi: 10.1016/j.survophthal.2017.06.009
4. Torres Netto, EA, Al-Otaibi, WM, Hafezi, NL, Kling, S, Al-Farhan, HM, Randleman, JB, et al. Prevalence of keratoconus in paediatric patients in Riyadh, Saudi Arabia. Br J Ophthalmol. (2018) 102:1436–41. doi: 10.1136/bjophthalmol-2017-311391
5. Rabinowitz, YS. Keratoconus. Surv Ophthalmol. (1998) 42:297–319. doi: 10.1016/S0039-6257(97)00119-7
6. Godefrooij, DA, Mangen, MJ, Chan, E, O’Brart, DPS, Imhof, SM, de Wit, GA, et al. Cost-effectiveness analysis of corneal collagen crosslinking for progressive keratoconus. Ophthalmology. (2017) 124:1485–95. doi: 10.1016/j.ophtha.2017.04.011
7. Davidson, AE, Hayes, S, Hardcastle, AJ, and Tuft, SJ. The pathogenesis of keratoconus. Eye. (2014) 28:189–95. doi: 10.1038/eye.2013.278
8. Gordon-Shaag, A, Millodot, M, Shneor, E, and Liu, Y. The genetic and environmental factors for keratoconus. Biomed Res Int. (2015) 2015:795738:1–19. doi: 10.1155/2015/795738
9. Seiler, T, Huhle, S, Spoerl, E, and Kunath, H. Manifest diabetes and keratoconus: a retrospective case-control study. Graefes Arch Clin Exp Ophthalmol. (2000) 238:822–5. doi: 10.1007/s004179900111
10. Woodward, MA, Blachley, TS, and Stein, JD. The association between sociodemographic factors, common systemic diseases, and keratoconus: an analysis of a nationwide heath care claims database. Ophthalmology. (2016) 123:457–65.e2. doi: 10.1016/j.ophtha.2015.10.035
11. Naderan, M, Naderan, M, Rezagholizadeh, F, Zolfaghari, M, Pahlevani, R, and Rajabi, MT. Association between diabetes and keratoconus: a case-control study. Cornea. (2014) 33:1271–3. doi: 10.1097/ICO.0000000000000282
12. Kuo, IC, Broman, A, Pirouzmanesh, A, and Melia, M. Is there an association between diabetes and keratoconus? Ophthalmology. (2006) 113:184–90. doi: 10.1016/j.ophtha.2005.10.009
13. Kosker, M, Suri, K, Hammersmith, KM, Nassef, AH, Nagra, PK, and Rapuano, CJ. Another look at the association between diabetes and keratoconus. Cornea. (2014) 33:774–9. doi: 10.1097/ICO.0000000000000167
14. Whelchel, AE, McKay, TB, Priyadarsini, S, Rowsey, T, and Karamichos, D. Association between diabetes and keratoconus: a retrospective analysis. Sci Rep. (2019) 9:13808. doi: 10.1038/s41598-019-50095-2
15. Hashemi, H, Heydarian, S, Hooshmand, E, Saatchi, M, Yekta, A, Aghamirsalim, M, et al. The prevalence and risk factors for keratoconus: a systematic review and meta-analysis. Cornea. (2020) 39:263–70. doi: 10.1097/ICO.0000000000002150
16. Emdin, CA, Khera, AV, and Kathiresan, S. Mendelian randomization. JAMA. (2017) 318:1925–6. doi: 10.1001/jama.2017.17219
17. Davies, NM, Holmes, MV, and Davey, SG. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. (2018) 362:k601. doi: 10.1136/bmj.k601
18. Pingault, JB, O'Reilly, PF, Schoeler, T, Ploubidis, GB, Rijsdijk, F, and Dudbridge, F. Using genetic data to strengthen causal inference in observational research. Nat Rev Genet. (2018) 19:566–80. doi: 10.1038/s41576-018-0020-3
19. Strawbridge, RJ, Dupuis, J, Prokopenko, I, Barker, A, Ahlqvist, E, Rybin, D, et al. Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes. Diabetes. (2011) 60:2624–34. doi: 10.2337/db11-0415
20. Morimoto, RI, and Cuervo, AM. Proteostasis and the aging proteome in health and disease. J Gerontol A Biol Sci Med Sci. (2014) 69:S33–8. doi: 10.1093/gerona/glu049
21. Elsworth, B, Lyon, M, Alexander, T, Liu, Y, Matthews, P, Hallett, J, et al. The MRC IEU OpenGWAS data infrastructure. bioRxiv. (2020): Available at: https://doi.org/10.1101/2020.08.10.244293. [Epub ahead of preprint]
22. Dastani, Z, Hivert, MF, Timpson, N, Perry, JR, Yuan, X, and Scott, RA. Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals. PLoS Genet. (2012) 8:e1002607. doi: 10.1371/journal.pgen.1002607
23. Scott, RA, Lagou, V, Welch, RP, Wheeler, E, Montasser, ME, Luan, J, et al. Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways. Nat Genet. (2012) 44:991–1005. doi: 10.1038/ng.2385
24. Wheeler, E, Leong, A, Liu, CT, Hivert, MF, Strawbridge, RJ, and Podmore, C. Impact of common genetic determinants of hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: a transethnic genome-wide meta-analysis. PLoS Med. (2017) 14:e1002383. doi: 10.1371/journal.pmed.1002383
25. Mahajan, A, Wessel, J, Willems, SM, Zhao, W, Robertson, NR, and Chu, AY. Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes. Nat Genet. (2018) 50:559–71. doi: 10.1038/s41588-018-0084-1
26. Hardcastle, AJ, Liskova, P, Bykhovskaya, Y, McComish, BJ, Davidson, AE, Inglehearn, CF, et al. A multi-ethnic genome-wide association study implicates collagen matrix integrity and cell differentiation pathways in keratoconus. Commun Biol. (2021) 4:266. doi: 10.1038/s42003-021-01784-0
27. Burgess, S, Butterworth, A, and Thompson, SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758
28. Hemani, G, Zheng, J, Elsworth, B, Wade, KH, Haberland, V, Baird, D, et al. The MR-base platform supports systematic causal inference across the human phenome. eLife. (2018) 7:7. doi: 10.7554/eLife.34408
29. Bowden, J, Davey Smith, G, Haycock, PC, and Burgess, S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
30. Verbanck, M, Chen, CY, Neale, B, and Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
31. Bowden, J, Davey Smith, G, and Burgess, S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
32. Burgess, S, Dudbridge, F, and Thompson, SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. (2016) 35:1880–906. doi: 10.1002/sim.6835
33. Bowden, J, Del Greco, MF, Minelli, C, Zhao, Q, Lawlor, DA, Sheehan, NA, et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. (2019) 48:728–42. doi: 10.1093/ije/dyy258
34. Moon, JY, Lee, J, Park, YH, Park, EC, and Lee, SH. Incidence of keratoconus and its association with systemic comorbid conditions: a nationwide cohort study from South Korea. J Ophthalmol. (2020) 2020:1–6. doi: 10.1155/2020/3493614
35. Claessens, JLJ, Godefrooij, DA, Vink, G, Frank, LE, and Wisse, RPL. Nationwide epidemiological approach to identify associations between keratoconus and immune-mediated diseases. Br J Ophthalmol. (2022) 106:1350–4. doi: 10.1136/bjophthalmol-2021-318804
36. Goldstein, DE, Little, RR, Lorenz, RA, Malone, JI, Nathan, D, Peterson, CM, et al. Tests of glycemia in diabetes. Diabetes Care. (2004) 27:1761–73. doi: 10.2337/diacare.27.7.1761
37. Georgakis, MK, Harshfield, EL, Malik, R, Franceschini, N, Langenberg, C, Wareham, NJ, et al. Diabetes mellitus, glycemic traits, and cerebrovascular disease: a Mendelian randomization study. Neurology. (2021) 96:e1732–42. doi: 10.1212/WNL.0000000000011555
38. Abdul Murad, NA, Abdullah, N, Kamaruddin, MA, Abd Jalal, N, Ismail, N, Yusof, NAM, et al. Discordance between fasting plasma glucose (FPG) and HbA1c in diagnosing diabetes and pre-diabetes in the Malaysian cohort. J ASEAN Fed Endocr Soc. (2021) 36:127–32. doi: 10.15605/jafes.036.02.02
39. Ho-Pham, LT, Nguyen, UDT, Tran, TX, and Nguyen, TV. Discordance in the diagnosis of diabetes: comparison between HbA1c and fasting plasma glucose. PLoS One. (2017) 12:e0182192. doi: 10.1371/journal.pone.0182192
40. Chen, Z, Bai, Y, Long, X, Luo, Q, Wen, Z, Li, Y, et al. Effects of adiponectin on T2DM and glucose homeostasis: a Mendelian randomization study. Diabetes Metab Syndr Obes. (2020) 13:1771–84. doi: 10.2147/DMSO.S248352
41. Bright, DJ, Dunseath, GJ, Peter, R, and Luzio, S. Stability of proinsulin in whole blood. Clin Biochem. (2018) 52:153–5. doi: 10.1016/j.clinbiochem.2017.10.005
42. Wareham, NJ, Byrne, CD, Williams, R, Day, NE, and Hales, CN. Fasting proinsulin concentrations predict the development of type 2 diabetes. Diabetes Care. (1999) 22:262–70. doi: 10.2337/diacare.22.2.262
43. Ates, KM, Estes, AJ, and Liu, Y. Potential underlying genetic associations between keratoconus and diabetes mellitus. Adv Ophthalmol Pract Res. (2021) 1:100005. doi: 10.1016/j.aopr.2021.100005
44. McKay, TB, Priyadarsini, S, and Karamichos, D. Mechanisms of collagen crosslinking in diabetes and keratoconus. Cells. (2019) 8:1239. doi: 10.3390/cells8101239
45. Buckingham, B, and Reiser, KM. Relationship between the content of lysyl oxidase-dependent cross-links in skin collagen, nonenzymatic glycosylation, and long-term complications in type I diabetes mellitus. J Clin Invest. (1990) 86:1046–54. doi: 10.1172/JCI114807
46. Chronopoulos, A, Tang, A, Beglova, E, Trackman, PC, and Roy, S. High glucose increases lysyl oxidase expression and activity in retinal endothelial cells: mechanism for compromised extracellular matrix barrier function. Diabetes. (2010) 59:3159–66. doi: 10.2337/db10-0365
47. Coral, K, Madhavan, J, Pukhraj, R, and Angayarkanni, N. High glucose induced differential expression of lysyl oxidase and its isoform in ARPE-19 cells. Curr Eye Res. (2013) 38:194–203. doi: 10.3109/02713683.2012.720341
Keywords: Mendelian randomization, keratoconus, fasting glucose, proinsulin levels, type 2 diabetes, genome-wide association study
Citation: Zhu X, Cheng D, Ruan K, Shen M and Ye Y (2023) Causal relationships between type 2 diabetes, glycemic traits and keratoconus. Front. Med. 10:1264061. doi: 10.3389/fmed.2023.1264061
Edited by:
Hong Qi, Peking University Third Hospital, ChinaCopyright © 2023 Zhu, Cheng, Ruan, Shen and Ye. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yufeng Ye, eXlmMDU3MUBtYWlsLmV5ZS5hYy5jbg==
†These authors have contributed equally to this work