Skip to main content

Genetic variants in the FOXO1 and ZNF469 genes are associated with keratoconus in Sweden: a case-control study

Abstract

Background

Keratoconus (KC) is characterized by pathological thinning and bulging of the cornea that may lead to visual impairment. The etiology of sporadic KC remains enigmatic despite intensive research in recent decades. The purpose of this study was to examine the relationship between previously highlighted genetic variants associated with KC and sporadic KC in a Swedish cohort.

Methods

A total of 176 patients (age 16–70 years) with sporadic KC diagnosed by Scheimpflug-topography (Pentacam) were included. The control group (n = 418; age 70 years) was a subsample originating from the Gothenburg H70 Birth Cohort Studies of ageing. Extraction of DNA from blood samples was performed according to standard procedures, and genotyping was performed using competitive allele specific PCR (KASP) technology. A total of 11 single nucleotide polymorphisms (SNPs) were selected for analysis.

Results

Statistically significant associations (p = 0.005) were found between the SNPs rs2721051 and rs9938149 and sporadic KC. These results replicate earlier research that found associations between genetic variants in the FOXO1 and BANP-ZNF469 genes and sporadic KC in other populations.

Conclusion

Genetic variations in the FOXO1 and BANP-ZNF469 genes may be involved in the pathogenesis of sporadic KC.

Peer Review reports

Background

Keratoconus (KC) is characterized by bilateral but asymmetric, pathological thinning and bulging of the cornea that may lead to visual impairment [1]. With a prevalence of approximately 1 in 2000 in the general population, KC is also a relatively common eye disease [2]. Spectacles, rigid gas-permeable contact lenses and insertion of corneal ring segments are treatment modalities commonly applied to improve visual function. Corneal collagen crosslinking (CXL) aims to halt the progression of KCs. In some cases, corneal transplantation is deemed necessary, thus replacing the thinned part of the cornea with an allogenic transplant. Eye rubbing, improper fit of contact lenses and chronic ocular irritation, have all been found to increase the risk for progression of KC [3]. In most cases, however, the etiology of KC remains enigmatic. The etiology is likely multifactorial, and several observations make a genetic component plausible. Autosomal dominant inheritance in families, higher incidence in families with one affected family member and high concordance among identical twins suggest a genetic cause of the disease [4, 5]. The results of the first genome-wide association study (GWAS) of KC by Li et al. in 2012 suggested an association of the single-nucleotide polymorphism (SNP) rs4954218 (located near the RAB3GAP1 gene) with the development of sporadic KC [6]. The findings of another GWAS by Burdon et al. indicated possible KC susceptibility through variations in the hepatocyte growth factor (HGF) gene [7]. In 2013, a comprehensive GWAS in a Caucasian population performed by Lu et al. identified 16 new loci that were associated with central corneal thickness (CCT), and two of these CCT-associated loci (FOXO1 and FNDC3B) were related to an increased risk of developing KC [8]. FOXO1 encodes the forkhead box 01 protein which is involved in the regulation of different cellular processes. These include the response of cells to oxidative stress which may contribute to corneal thinning. FOXO1 may also influence cell proliferation and differentiation which could affect corneal tissue. FNDC3B codes for a fibronectin type III protein and variations in this gene may hypothetically also affect cornea tissue leading to thinning and deformation.

The aim of this study was to examine the relationship between previously highlighted genetic variants associated with KC and central corneal thickness (CCT) and sporadic KC in a Swedish cohort.

Methods

Patients with sporadic KC

Patients with sporadic KC in this study were examined during the recruitment process of a randomized controlled clinical trial aiming to study the efficacy of CXL in progressive KC. The trial was approved by the regional ethical review board in Gothenburg, Sweden (approval number 949 − 11), and the study followed the terms of the Declaration of Helsinki. Informed written consent was obtained from each patient. Blood samples were drawn from each participant and stored for future analyses according to the Swedish Biobank Law. All study participants gave written consent regarding biobanking.

Inclusion criteria were age of 17 to 35 years and tomographic diagnosis of KC as defined by the Pentacam Scheimpflug camera topometric module (KC grade 1 or higher), the ability to stop contact lens wear (rigid and soft) at least two weeks prior to exam, and that three consecutive, topographic exams of each eye with approved quality rating at first and second visits could be obtained.

Diagnosis of KC in patients with sporadic KC

Both eyes of the patients with sporadic KC were examined with three consecutive tomography images. The measurements were performed by either of two experienced examiners (BS, UM). The tomography instrument was calibrated according to the manufacturers’ recommendations. All measurements were performed under scotopic light conditions. The patient had both eyes open during measurements. The patient fixated with the eye to be measured at the red fixation light of the Pentacam system and was instructed to blink before every measurement. The automated release mode with 25 pictures was used. The Pentacam KC staging system was used to diagnose KC in combination with individual assessment of the tomography results by a corneal surgeon (WW) with regard to the plausibility of the diagnosis. The tomographical diagnosis of KC was based on keratometry values, minimal corneal thickness, typical distribution of corneal thickness in relation to keratometry values and presence of pathological elevations of the anterior and / or poster corneal surface. Generally, keratometry values above 48 diopters and posterior corneal elevations of more than 35 μm were considered as pathological. The diagnosis of KC, however, also weighed in the overall picture as described above. All patients underwent a complete slit-lamp exam.

Controls

The Gothenburg H70 Birth Cohort Studies examine cohorts of 70-year-olds with the aim of learning more about aging. The participants in the control group in our study were chosen from a cohort of individuals born in 1944 who participated in the H70 study during 2014–2016 [9]. A subsample of 560 patients underwent extended ophthalmic examinations [10]. The H70 study was approved by the regional ethical review board, and informed, written consent, including biobankning, was obtained from each patient.

Absence of corneal ectasia in the controls

All participants in the control group underwent comprehensive ophthalmic examination, including objective refraction with an auto kerato-refractometer. Individuals with central corneal astigmatism below 1.5 diopters (D) were chosen as the control group under the assumption that the likelihood for KC was sufficiently low, especially since they were all 70 years of age and the disease rarely manifests itself or progresses after the age of 30.

SNP selection and genotyping

The selection of SNPs (i.e. FOXO1: rs2721051, COL5A1: rs7044529, RXRA-COL5A1: rs7044529, FINDC3B: rs4894535, BANP-ZNF469: rs9938149, RAB3GAP1: rs4954218, and HGF: rs2286194, rs17501108, and rs3735520) were based on findings in previous GWASs on KC and CCT [6,7,8]. Furthermore, SNPs in the LOX gene (rs10519694 and rs2956540) that have been associated with KC in several previous case-control studies were selected [11, 12]. Extraction of DNA from blood samples was performed according to standard procedures. Genotyping was performed at LGC Genomics (Hoddesdon, Herts, UK) using competitive allele specific PCR (KASP) technology. The success rate was over 95% for all genotyped SNPs, and they were in Hardy-Weinberg equilibrium. One of the selected SNPs (rs1324183 in MPDZ-NF1B) was not genotyped due to unsuccessful design of the KASP assay.

Statistical analysis

IBM SPSS Statistics software (version 28.0, IBM, Armonk, NY, USA) was used for statistical analysis. As descriptive statistics, means and standard deviations (SD) were calculated as well as proportions for categorical data. Fisher’s exact test was used for statistical analysis of the latter. Student’s t-test was used for statistical analysis of continuous data. A p-value < 0.05 was considered statistically significant.

Associations between sporadic KC and genetic variants were analyzed using logistic regression models. Two different models were tested: an additive genetic model (i.e. Rare genotypes are coded as 2, heterozygotes are coded as 1, and common genotypes are coded as 0) and a dominant genetic model (i.e. Rare genotypes and heterozygotes are coded as 1, and common genotypes are coded as 0. Only individuals of Caucasian origin (163 patients and 406 controls) were included in the analyses. The analyses were also repeated after excluding individuals born outside of Europe (including 153 patients and 399 controls in the analyses).

Post hoc power was calcualted using the G*Power software (version 3.1.9.7), including the effect sizes found in the genetic association analyses (dominant model) and a p-value level of 0.05.

Results

The sample characteristics of the KC patients and the control individuals are shown in Table 1. Significant differences were seen for sex and age between the two groups (p < 0.001). Notably, in the KC group, only 11% were women.

Table 1 Baseline characteristics of the studied individuals

The results for the genetic analysis are summarized in Table 2. Statistically significant associations were found between the SNPs rs2721051 in FOXO1 and rs9938149 in BANP-ZNF469 and sporadic KC. For rs2721051, the A allele was associated with an increased risk of disease, and for rs9938149, the C allele was associated with a decreased risk of disease (i.e. the C allele contributed with a protective effect). These associations were seen both in the additive and dominant models. None of the other investigated SNPs were associated with KC in our sample. When repeating the analyses after excluding non-European individuals, the associations with the two SNPs in FOXO1 and BANP-ZNF469 remained (results not shown). The allele frequencies of these two SNPs in the control population were similar to those reported for other European samples in the genetic data base dbSNP [13].

Table 2 Summary for the genetic analysis

Post hoc power to detect a significant finding showed a power of 81–100% for five SNPs (FOXO1: rs2721051, BANP-ZNF469: rs9938149, COL5A1: rs7044529, HGF: rs17501108, and LOX: rs2956540), a power of 50 to 71% for five SNPs (RXRA-COL5A1: rs7044529, FINDC3B: rs4894535, HGF: rs2286194, and LOX: rs10519694) and a power of 13% for one SNP (HGF: rs3735520).

Discussion

For the first time in a Swedish cohort, our study showed associations of two SNPs with sporadic KC. The SNPs were located in the genes FOXO1 and BANP-ZNF469.

As outlined earlier, the selection of SNPs was based on findings in previous GWASs on KC and CCT [6,7,8]. SNPs in the LOX gene (rs10519694 and rs2956540) that have been associated with KC in several previous case-control samples were also selected [11, 12].

FOXO1 is a member of the forkhead box (Fox) transcription factor family, and Fox proteins regulate cellular oxidative stress [14]. Increased oxidative stress may damage the trabecular meshwork and lead to elevated intraocular pressure and corneal thinning [15]. In both primary open-angle glaucoma (POAG) and KC, CCT is known to be reduced [16, 17]. In 2010, Lu et al. identified a SNP near FOXO1 (rs2721051), which accounted for approximately 1.2% of the variation in normal CCT [18]. Another study by Lu et al. found a strong association between KC risk and FOXO1 (OR = 1.62, 95% CI 1.4–1.88, P = 2.7 × 10− 10), which is comparable to the findings in our study [8].

Rare variants in the ZNF469 gene have been reported to be involved in the development of brittle cornea syndrome (BCS), which leads to thinning of the cornea and progressive visual impairment [19]. Common variants near ZNF469 were associated with CCT according to Lu et al. [18]. Furthermore, based on the finding of a KC locus close to the location of ZNF469 in the Finnish family, Tyynismaa et al. argued that mutations in the ZNF469 gene might be found in patients with KC [20]. The role of ZNF469 in the development of KC is, however, not clear. Regarding common genetic variants, our study confirmed the findings by Lu et al., showing that the CCT-reducing allele near ZNF469 (rs9938149), actually, is associated with a decreased KC risk [8].

Very recently, Hardcastle et al. reported the results of the first large-scale genome-wide association study of KC [21]. This multicenter study included 4669 cases and 116,547 controls and revealed genome-wide significant associations for 36 genomic loci, including those found to be associated with KC in our study (i.e. FOXO1 and ZNF 469). The study clearly supports a genetic component of KC. It was estimated that 12.5% of the heritability of KC in European individuals was explained by the common variants associated with KC in this study.

Among the limitations of our study is that exclusion of KC in the control group was not based on topographic examinations. The effect of this limitation should, however, be small. The prevalence of sporadic KC is rather low, and individuals with central corneal astigmatism above 1.5 diopters were excluded. The high age of the control subjects, 70 years, also precludes later development of the disease. Another limitation is that due to a relatively small sample size only five SNPs reached a power of over 80%, when post hoc power to detect associations were calculated. It is, however, important to note that for some SNPs very small effect sizes were used in the power calculations, since we included the ORs found for each specific SNP. For some of the SNPs, such as those in the LOX gene, the effect sizes (i.e. ORs) in our study are of comparable size to previous studies [22], but the p-values do not reach the level of significance. However, many of the nonsignificant results in our study show an effect size below the findings in previous studies, indicating that these genetic variants are of low importance in relation to KC in our sample [22].

Finally, it is well known that KC is more common in males [23]. This is also reflected in the significant difference between the number of males in the study group and the control group in this study. This might affect the significance of the results. It is, therefore, important to point out the that the results in this study are based on logistic regression models adjusted for sex.

Conclusion

This study demonstrates associations between two SNPs in FOXO1 and BANP-ZNF469 and KC in a Swedish sample. The results corroborate previously reported associations between these genes and sporadic KC.

Data availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. Confidential data will, however, not be shared.

References

  1. Rabinowitz YS, Rasheed K, Yang H, Elashoff J. Accuracy of ultrasonic pachymetry and videokeratography in detecting keratoconus. J Cataract Refract Surg. 1998;24(2):196–201.

    Article  CAS  PubMed  Google Scholar 

  2. Kennedy RH, Bourne WM, Dyer JA. A 48-year clinical and epidemiologic study of keratoconus. Am J Ophthalmol. 1986;101(3):267–73.

    Article  CAS  PubMed  Google Scholar 

  3. Gomes JA, Tan D, Rapuano CJ, Belin MW, Ambrosio R Jr., Guell JL, Malecaze F, Nishida K, Sangwan VS. Group of panelists for the Global Delphi Panel of K et al: global consensus on keratoconus and ectatic diseases. Cornea. 2015;34(4):359–69.

    Article  PubMed  Google Scholar 

  4. Tuft SJ, Hassan H, George S, Frazer DG, Willoughby CE, Liskova P. Keratoconus in 18 pairs of twins. Acta Ophthalmol. 2012;90(6):e482–486.

    Article  PubMed  Google Scholar 

  5. Wang Y, Rabinowitz YS, Rotter JI, Yang H. Genetic epidemiological study of keratoconus: evidence for major gene determination. Am J Med Genet. 2000;93(5):403–9.

    Article  CAS  PubMed  Google Scholar 

  6. Li X, Bykhovskaya Y, Haritunians T, Siscovick D, Aldave A, Szczotka-Flynn L, Iyengar SK, Rotter JI, Taylor KD, Rabinowitz YS. A genome-wide association study identifies a potential novel gene locus for keratoconus, one of the commonest causes for corneal transplantation in developed countries. Hum Mol Genet. 2012;21(2):421–9.

    Article  PubMed  Google Scholar 

  7. Burdon KP, Macgregor S, Bykhovskaya Y, Javadiyan S, Li X, Laurie KJ, Muszynska D, Lindsay R, Lechner J, Haritunians T, et al. Association of polymorphisms in the hepatocyte growth factor gene promoter with keratoconus. Invest Ophthalmol Vis Sci. 2011;52(11):8514–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lu Y, Vitart V, Burdon KP, Khor CC, Bykhovskaya Y, Mirshahi A, Hewitt AW, Koehn D, Hysi PG, Ramdas WD, et al. Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus. Nat Genet. 2013;45(2):155–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Rydberg Sterner T, Ahlner F, Blennow K, Dahlin-Ivanoff S, Falk H, Havstam Johansson L, Hoff M, Holm M, Horder H, Jacobsson T, et al. The Gothenburg H70 birth cohort study 2014-16: design, methods and study population. Eur J Epidemiol. 2019;34(2):191–209.

    Article  PubMed  Google Scholar 

  10. Havstam Johansson L, Skiljic D, Falk Erhag H, Ahlner F, Pernheim C, Rydberg Sterner T, Wetterberg H, Skoog I, Zetterberg M. Vision-related quality of life and visual function in a 70-year-old Swedish population. Acta Ophthalmol 2020.

  11. Bykhovskaya Y, Li X, Epifantseva I, Haritunians T, Siscovick D, Aldave A, Szczotka-Flynn L, Iyengar SK, Taylor KD, Rotter JI, et al. Variation in the lysyl oxidase (LOX) gene is associated with keratoconus in family-based and case-control studies. Invest Ophthalmol Vis Sci. 2012;53(7):4152–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Dudakova L, Jirsova K. The impairment of lysyl oxidase in keratoconus and in keratoconus-associated disorders. J Neural Transm (Vienna). 2013;120(6):977–82.

    Article  CAS  PubMed  Google Scholar 

  13. dbSNP Gd.

  14. Klotz LO, Sanchez-Ramos C, Prieto-Arroyo I, Urbanek P, Steinbrenner H, Monsalve M. Redox regulation of FoxO transcription factors. Redox Biol. 2015;6:51–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Siegfried CJ, Shui YB, Bai F, Beebe DC. Central corneal thickness correlates with oxygen levels in the human anterior chamber angle. Am J Ophthalmol. 2015;159(3):457–462e451.

    Article  PubMed  Google Scholar 

  16. Gordon MO, Beiser JA, Brandt JD, Heuer DK, Higginbotham EJ, Johnson CA, Keltner JL, Miller JP, Parrish RK 2nd, Wilson MR, et al. The ocular hypertension treatment study: baseline factors that predict the onset of primary open-angle glaucoma. Arch Ophthalmol. 2002;120(6):714–20. discussion 829– 730.

  17. Rabinowitz YS, Keratoconus. Surv Ophthalmol. 1998;42(4):297–319.

    Article  CAS  PubMed  Google Scholar 

  18. Lu Y, Dimasi DP, Hysi PG, Hewitt AW, Burdon KP, Toh T, Ruddle JB, Li YJ, Mitchell P, Healey PR, et al. Common genetic variants near the Brittle Cornea syndrome locus ZNF469 influence the blinding disease risk factor central corneal thickness. PLoS Genet. 2010;6(5):e1000947.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Abu A, Frydman M, Marek D, Pras E, Nir U, Reznik-Wolf H, Pras E. Deleterious mutations in the zinc-finger 469 gene cause brittle cornea syndrome. Am J Hum Genet. 2008;82(5):1217–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tyynismaa H, Sistonen P, Tuupanen S, Tervo T, Dammert A, Latvala T, Alitalo T. A locus for autosomal dominant keratoconus: linkage to 16q22.3-q23.1 in Finnish families. Invest Ophthalmol Vis Sci. 2002;43(10):3160–4.

    PubMed  Google Scholar 

  21. Hardcastle AJ, Liskova P, Bykhovskaya Y, McComish BJ, Davidson AE, Inglehearn CF, Li X, Choquet H, Habeeb M, Lucas SEM, et al. A multi-ethnic genome-wide association study implicates collagen matrix integrity and cell differentiation pathways in keratoconus. Commun Biol. 2021;4(1):266.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Rong SS, Ma STU, Yu XT, Ma L, Chu WK, Chan TCY, Wang YM, Young AL, Pang CP, Jhanji V, et al. Genetic associations for keratoconus: a systematic review and meta-analysis. Sci Rep. 2017;7(1):4620.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Godefrooij DA, de Wit GA, Uiterwaal CS, Imhof SM, Wisse RP. Age-specific incidence and prevalence of Keratoconus: a Nationwide Registration Study. Am J Ophthalmol. 2017;175:169–72.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors are grateful to the patients for participating in this study. The authors would also like to thank the clinical and administrative staff at the Department of Ophthalmology, Sahlgrenska University Hospital, for their support.

Funding

This work was supported by grants from the Swedish government (“Agreement concerning research and education of doctors”; ALF-GBG-441721 and ALF-GBG-725041), De Blindas Vänner, Kronprinsessan Margaretas Arbetsnämnd för Synskadade. The H70 study of individuals born in 1944 was financed by grants from the Swedish state under the agreement between the Swedish government and the county councils, the ALF agreement (ALF-GBG-716681), the Swedish Research Council (2012–5041, 2015–02830, 2019 − 01096, 2013–8717, 2017 − 00639), the Swedish Research Council for Health, Working Life and Welfare (2013 − 1202, 2018 − 00471, AGECAP 2013–2300, 2013–2496), Konung Gustaf V:s och Drottning Victorias Frimurarestiftelse, Hjärnfonden, Alzheimerfonden, Eivind och Elsa K:son Sylvans stiftelse.

Open access funding provided by University of Gothenburg.

Author information

Authors and Affiliations

Authors

Contributions

W.W. conducted the clinical research, interpreted the results and wrote the manuscript. M.Z. conceived and designed the study. A.Z. conducted the genetic analysis, interpreted the results and wrote parts of the manuscript. B.S. and U.M. participated in the clinical part of the study, examining patients and collecting research data. L.H.J. and I.S. conducted the H70 study and provided their data for the control group in this study. M.C. participated in the conception and analysis of the study. All authors contributed to writing and revision of the manuscript and all authors read and approved the final manuscript.

Corresponding author

Correspondence to Wolf Wonneberger.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the regional ethical review board in Gothenburg, Sweden (approval number 949 − 11), and the study followed the terms of the Declaration of Helsinki. Written informed consent was obtained from all participants in the study. Blood samples were stored according to the Swedish Biobank Law. All study participants gave written informed consent regarding biobanking.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wonneberger, W., Sterner, B., MacLean, U. et al. Genetic variants in the FOXO1 and ZNF469 genes are associated with keratoconus in Sweden: a case-control study. BMC Ophthalmol 24, 36 (2024). https://doi.org/10.1186/s12886-024-03299-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12886-024-03299-8

Keywords