Skip to main content

Inflammatory and genomic interactions within keratoconus susceptible patients: a nationwide registered case–control study

Abstract

Purpose

This study aimed to investigate the association between variants in the interleukin (IL)-1 gene cluster and susceptibility to keratoconus (KC) in an Iranian population.

Methods

In the case group, there were 188 KC patients diagnosed by clinical findings and corneal tomography. The control group included all 205 healthy controls with no personal or family history of eye-related, metabolic, or immune system-related disease. Using the standard salting out extraction procedure, genomic DNA was isolated from peripheral blood leukocytes. The genotypes were determined by applying agarose gel electrophoresis for the IL-1RN 86 bp VNTR and polymerase chain reaction-restriction fragment length polymorphism (PCR–RFLP) for rs16944 and rs1143634.

Results

The results showed a significant association between the IL-1β rs1143634 (rs1143634 T allele, P = 0.008) and IL-1RN 86 bp VNTR polymorphisms (LL and LS genotype, P = 0.048 and 0.012 respectively) and susceptibility to KC in the Iranian population. The genotype distributions of rs1143634 (P = 0.004) and rs2234663 (P = 0.042) significantly differed between case and control groups, with certain genotypes demonstrating a protective effect against KC. Logistic regression analysis revealed a protective effect of the IL-1RN L allele [odds ratio (OR) = 0.367, 95% confidence interval (CI): 0.240–0.562; P = 0.000] and certain haplotypes (OR = 0.628, 95% CI: 0.447–0.884; P = 0.007) against KC. However, no significant association was found for the IL-1β rs16944 polymorphism.

Conclusion

This study provides evidence for an association between variants in the IL-1 gene cluster and susceptibility to KC in an Iranian population. Further research on larger and more diverse populations is warranted to validate these findings and explore the underlying mechanisms involved.

Background

Keratoconus (KC) is a multifactorial [1], bilateral and asymmetrical eye disease that leads to progressive thinning and steepening of the cornea [2]. This leads to irregular astigmatism and reduced visual acuity [3]. The prevalence of KC is estimated to be 50–2300 cases per 100,000 people affected by geologic variation, age groups, and screening technology [4,5,6] and 0.15% of the vision service plan enrollers in the United States [7]. A recent study from Yazd province, Iran, estimated the KC rate of incidence to be 22.3 to 24.9 per 100,000 people annually [8]. In the Middle East, the prevalence of KC in the Syrian student population was estimated at 1.43% [9].

The primary environmental and behavioral risk factors of KC pathogenesis are constant eye rubbing [10] and contact lens wear [11]. Several lines of evidence also point to a genetic determinant for KC [12]; for example, familial inheritance [13] and concordance between monozygotic twins as opposed to dizygotic twins [14]. KC is a convoluted attribute that presumably involves varying penetrance, several genes [15, 16] environmental influences, and different pathways involved [17]. Elevated proteinase activity [7, 11], reduced levels of proteinase inhibitors [13, 14, 17], exacerbated oxidative damage, and keratocyte apoptosis [18, 19] are among the different mechanisms involved in the pathogenesis of the disease.

Inflammation is an important factor in the pathogenesis of KC. However, the extent of its contribution to KC is not yet known [18]. A family of 11 cytokines called interleukin (IL)-1 controls immunologic and inflammatory responses, particularly during corneal wound healing. The two primary cytokines that are released after corneal injury are IL-1α and IL-1β, both of which cause remodeling by upregulating matrix metallopeptidase (MMPs) and inducing apoptosis. IL-1 is secreted by corneal epithelial cells and enters the stroma via the tear film [18]. IL-1RN is an antagonist of IL-1 that inhibits the activities of IL-1α and IL-1β and regulates several immunological and inflammatory reactions linked to IL-1, especially during the acute stage of infection and inflammation. IL-1α, IL-1β and IL-1RN are located as clusters on chromosome 2q14 [18, 19].

This study investigated the possible association between the IL-1 gene cluster and the presence of KC in the Iranian population. We performed a genetic association study for variants in the IL-1β (rs16944 and rs1143634) and intron 2 of IL-1RN (rs2234663, 86bp VNTR). We also evaluated these variants' linkage disequilibrium (LD) and allelic interactions. We successfully demonstrated a statistically significant association between IL-1β rs1143634 and IL-1RN 86bp VNTR and susceptibility to KC in the Iranian population.

Methods

Study design, clinical examination and sample collection

This case–control genetic association study was conducted in Negah Eye Hospital, Tehran, Iran, between June 2021 and April 2023. The Committees approved the study protocol for Ethics of Tehran University of Medical Sciences (IR.SBMU.ORC.REC.1399.015). Written informed consent was obtained from the study participants in line with the Declaration of Helsinki before their inclusion in the study.

In order to elucidate the relationship for both this study and subsequent immunogenetic studies, four groups of KC cases were formed, each representing a different degree of visual limitation, as follows: Stage 4: corrected distance visual acuity (CDVA) less than 20/400, logMAR units < 0.05; Stage 3: CDVA less than 20/100, logMAR units < 0.2; Stage 2: CDVA less than 20/40, logMAR units < 0.5; Stage 1: CDVA less than 20/20, logMAR units < 1; Stage 0: CDVA = 20/20, logMAR units = 1 [3, 20,21,22]. The Pentacam HR corneal tomography (Oculus, Wetzlar, Germany) and the Corvis ST (OCULUS Optikgeräte GmbH; Wetzlar, Germany) tests were done for all suspicious cases to detect mild forms of KC (Supplementary Tables 1 and 2) [2, 23,24,25,26]. Suspect cases were defined as stage 0. All cases were bilateral and diagnoses were approved by two cornea specialists (FD, MAJ). Complete medical histories were obtained from all patients with KC. No other ocular disorders or risk factors (including contact lens use, eye rubbing, atopy, etc.) or systematic diseases (including connective tissue disorders) were uncovered. All healthy controls had no personal or familial history of eye-related, metabolic, or immune system-related disease. All participants or their legal guardians signed consent forms.

Isolation of genomic DNA and SNP genotyping

A K3EDTA tube was used to retain the 5 mL of peripheral blood drawn from the median cubital vein. Using the standard salting out extraction procedure, genomic DNA was isolated from peripheral blood leukocytes. The samples were then validated by electrophoresis on a 1% agarose gel, and their quantities were determined using ultraviolet (UV) spectroscopy and optical density ratios of 260/230 and 260/280. Polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) were used for genotyping of the targeted single nucleotide polymorphisms (SNPs). The initial stage of the PCR reactions was performed at 94 °C for 5 min before going through 35 cycles of 94 °C for 30 s, 58 °C for 30 s, 72 °C for 30 s, and an additional extension at 72 °C for 5 min. Ava1 and Tag1 restriction enzymes were used overnight to digest the target PCR product for rs16944 and rs1143634 polymorphisms, respectively. The processed PCR product underwent enzymatic digestion before being subjected to 2% agarose gel electrophoresis. By applying agarose gel electrophoresis to determine the copy number in the 86 bp-repeat, the IL-1RN 86 bp VNTR was evaluated. Supplementary Table 3 lists the forward and reverse primers, PCR products, and digested product lengths for the rs2234663, rs1143634, and rs16944 polymorphisms.

Statistical analysis

The SPSS (version 20.0, SPSS Inc., Chicago, Illinois, USA) software was used to analyze the data. Using the Hardy–Weinberg genetic balancing test, the frequency of genotype distribution in the case and control groups was assessed. The P value and χ2 value were calculated using the Chi-squared test to examine the genotype and allele frequency distribution variations and KC. The odds ratio (OR) and its 95% confidence interval (CI) were calculated using binary logistic regression analysis to examine the relative risk of illness associated with genotypes or alleles at polymorphic locations. Haplotype analysis was also performed, and LD between the SNPs was evaluated using the SNPAnalyzer software. Protective/predisposing effects of haplotypes were estimated by logistic regression analysis. Differences with P values less than 0.05 were regarded as statistically significant for all two-sided probability tests.

Results

Three hundred and ninety-three Iranian subjects were enrolled in this study, consisting of 205 healthy individuals (sex, 106 males and 99 females; age, mean ± SD: 33.26 ± 7.04 years) and 188 KC patients (sex, 97 males and 91 females; age, mean ± SD: 31.78 ± 9.18 years). The mean age distribution of the case and control groups did not differ substantially, and the demographics (age and sex) were evenly distributed across males and females (P = 0.982 and P = 0.077, respectively).

The results in Table 1 demonstrate the pertinent data of the genotype distribution of IL-1β gene rs1143634, rs16944, and IL-1RN VNTR polymorphisms. Although the genotype frequency distribution of rs2234663 for the case and control groups was not in the Hardy–Weinberg equilibrium (HWE) (P < 0.05, Table 1), rs1143634 and rs16944 genotype distribution were in the HWE (P > 0.05, Table 1). Subgroup analyses by age and sex were also performed. The genotype distributions of rs1143634 in men/women of more than 35 years of age were statistically different between case and control groups (P < 0.05 and P = 0.034, respectively). The genotype distributions of rs2234663 were significantly different only in the men's group (P = 0.016). Subgroup analysis based on sex and age for rs16944 revealed no statistically significant genotype distribution between case and control groups.

Table 1 Genotype distribution of IL-1RN 86 bp VNTR and IL-1β gene rs1143634, rs16944 polymorphisms in keratoconus cases and healthy controls

Association of rs2234663 (IL-1RN 86 bp VNTR) and IL-1β gene rs1143634 and rs16944 polymorphisms with KC

Tables 2, 3, and 4 show the logistic regression analysis and frequency distribution of the rs2234663 VNTR of IL-1RN, rs1143634, and rs16944 polymorphisms of the IL-1β gene in KC and the control group. Based on the differences of 86-bp VNTR, different alleles were expected. However, only IL-1RN-S (2-repeats) and IL-1RN-L (4 and 5 repeats) were detected in this study. We did not detect IL-1RN VNTR of 1, 3 and 6 replications in our population. The frequency distribution of LS heterozygous and LL homozygous genotypes significantly differed between the KC cases and the control group (P = 0.048 and P = 0.012, respectively). Indeed, the genotypes LS and LL were found to exert a protective effect against KC (OR = 0.431, 95% CI: 0.185–1.000; OR = 0.4, 95% CI: 0.192–0.832). The findings also revealed that, in comparison to the frequency distribution of the S allele in rs2234663, the frequency distribution of the L allele was lower in the case group than in the control group (77.9% vs. 90.6%), the difference of which was statistically significant (P = 0.000) and revealed a protective effect for the allele L of rs2234663 polymorphism against KC (OR = 0.367, 95% CI: 0.240–0.562).

Table 2 Association of IL -1RN 86 bp VNTR with keratoconus cases
Table 3 Association of IL-1β gene rs1143634 polymorphism with keratoconus cases
Table 4 Association of IL-1β gene rs16944 polymorphism with keratoconus cases

In the dominant model, the genotype distribution of LS + LL in the case–control group was statistically significant when different inheritance modes were taken into account (P = 0.011), and the regression analysis showed a protective role of this inheritance model (OR = 0.401, 95% CI: 0.194–0.828). However, in the recessive and codominant models, the frequency distribution was not significantly different (P = 0.105, P = 0.953, respectively).

Chi-squared test and logistic regression were applied to investigate the genotype/allele distribution and association of IL-1β gene rs1143634 SNP with KC risk. As indicated in Table 3, the results disclosed that the IL-1β rs1143634 allele had a significantly different distribution in cases compared to controls. However, TC and TT genotype distribution were not significantly different between case and control groups, where the allele C and genotype CC were considered the reference wild type allele/genotype. Our findings showed that individuals with allele T have an increased risk for KC disorder (OR = 1.53, 95% CI: 1.119–2.106; P = 0.008). Furthermore, the model analysis revealed that the carriers of the dominance model (TC + TT vs. CC) and codominant model (TC vs. TT + CC) had an increased risk of KC (P = 0.001). However, we found that carriers of the recessive model (TT vs. TC + CC) had no real relationship with the other KC progression risk factors (P > 0.05).

The variation between the two groups, using CC as the reference type, was not statistically significant for both the TC and CC genotypes with regards to the frequency distribution and logistic regression analysis of the IL-1β gene rs16944 (Table 4, P > 0.05). Logistic regression analysis revealed no protective/predisposing association with KC in the studied population CC and TC genotypes. Additionally, the study of several inheritance models revealed no variance in the frequency distribution of the CC and TC genotypes that could be considered statistically significant, and thus, following logistic regression analysis, there was no significant correlation with the examined disease (P > 0.05). Allele distribution was also not statistically different in cases compared to controls, and accordingly, a significant association was not observed with KC in study participants.

Haplotype analysis

The frequency of haplotypes (from left to right rs16944, rs1143634, and rs2234663) was assessed using SNPanalyzer (Table 5). We found eight haplotypes (h1–h8), where the h1 with CTL sequence was the most frequent haplotype (0.348). We further analyzed haplotypes to investigate whether any haplotype was linked to KC (Table 5). Because of the low frequency of h8 (0.009), it was omitted from downstream analysis. Considering haplotypes with frequencies of ≥ 1%, logistic regression analysis indicated a haplotype effect for h2 and h7. Our data revealed that haplotype h2 carriers had a lower risk of KC (OR = 0.628, 95% CI: 0.447–0.884; P = 0.007). In contrast, the h7 haplotype significantly increased the risk of KC (OR = 2.448, 95% CI: 1.35–4.437; P = 0.002]. We also performed LD analysis of the selected polymorphisms. Supplementary Fig. 1 shows the haplotype block formed by the LD study, which contains the two SNPs rs16944 and rs1143634 [D′ = 0.49216, r2 = 0.08612, logarithm of the odds (LOD) = 32.61988, P = 0.000).

Table 5 Distribution of different haplotypes of IL-1 cluster variants (rs16944, rs1143634, and rs2234663) and the association with keratoconus

Discussion

Apoptosis is considered an important process in KC pathogenesis. Indeed, KC corneas exhibit a hyper apoptotic phenotype [27], the apoptosis in keratocytes is thought to be triggered by IL-1. IL-1 is released in response to chronic mechanical injury to the cornea (e.g., eye rubbing) [28]. Since a genetic predisposition to the disease allows an environmental trigger to cause an inflammatory component, we examined whether variants in the IL-1 gene cluster confer susceptibility to KC in the Iranian population.

This study found a strong protective role for the IL-1RN L allele against KC. The genotypes LL and LS were associated with a reduced risk of KC (OR = 0.4 and 0.431, respectively). This association was also seen with a dominant mode of inheritance. However, no significant association between IL-1RN VNTR and the likelihood of KC was found in other experiments by Kim et al. and Palamar et al. [29, 30]. The difference, which may be explained in part by disparities in sample size, population, and ethnicity, warrants further studies on large populations of diverse ethnic groups.

The IL-1RN protein is an important negative regulator of the inflammatory response. This gene, which belongs to the IL-1 cytokine family, hinders the function of IL-1 (IL-1α and IL-1β) by interacting with the receptor IL-1R1 and keeping the coreceptor IL-1RAP (IL-1 receptor accessory protein) from associating with it for signal transduction (Fig. 1). As predicted by Signor (https://signor.uniroma2.it/), the gene is positively upregulated by NFκb- p65/p50, STAT3, and IL -10 (Fig. 1). In fact, IL-10 activates STAT3, triggering it to make the IL-1RN promoter accessible to nuclear NF-κB [31]. This highlights the value of further investigation of this pathway for future studies on KC.

Fig. 1
figure 1

Regulatory network of IL -1RN as predicted by Signor (https://signor.uniroma2.it/); IL -1RN is positively upregulated by NFκb-p65/p50, STAT3, and IL-10

IL-1β is located in the IL-1 gene cluster and contains putative functional polymorphisms associated with various human diseases [32]. The two IL -1β SNPs studied here were rs16944, C > T, and rs1143634, C > T. Our current findings on the role of rs16944 in KC is controversial; a study on an Egyptian population found the C allele to be associated with greater risk of KC [33], while in another study on the Chinese Han population, the carriers of the T allele had greater risk for KC [34]. Additionally, Palmar et al. reported no significant association between rs16944 genotypes and no susceptibility to KC was detected [30]. Although we did not detect a significant correlation between the genotype and allele frequencies of the IL-1 rs16944 variant and KC risk, it was associated with KC risk in the form of the CCL (haplotype h2) and TCS (haplotype h7) haplotypes. This observation underscores the importance of considering the role of each genetic variant, not as individual players, but as nodes in an intertwined pathway.

The variant IL-1β rs16944, C > T, is located at position −511 within the promoter region. Regulatory sequence variations can strongly affect the expression levels of their corresponding gene. Supporting evidence is available in the case of rs16944, where the T allele of IL-1β at position −511 has been reported to result in modestly increased transcriptional activity [29, 35]. Furthermore, homozygosity for the T allele at −511 and the C allele at −31 is correlated with increased IL-1β secretor phenotype [33]. The expression of IL-1β has been elevated in KC compared to normal cornea [36]. The second IL-1β polymorphism studied here is rs1143634, C > T. This SNP is located at the exon 5 position + 3954 and is considered a synonymous variant (NP_000567.1:p.Phe105 exon 5 =). This variant was previously associated with the risk of disorders like gastric cancer, endometriosis, and non-Sjogren dry eye [37,38,39]. There are a few reports on rs1143634, C > T, and susceptibility to KC. According to our results, the existence of the rs1143634 T allele could increase the risk of KC by 1.53 fold. In addition, the genotype TC could significantly increase the risk of KC in dominant and codominant modes of inheritance. However, in the Korean population, rs1143634 has not been associated with the risk of KC. The contradictory results again emphasize the significance of additional investigations in larger populations across various ethnic backgrounds.

Clinical applications including the diagnosis, treatment, and management of KC hold great promise based on the genetic results of this study. The discovery of certain IL-1 gene cluster variations linked to an increased risk of developing KC may result in the developing of genetic screening instruments that assist in identifying people who are more susceptible to the disease, which allows for earlier and more focused therapies. Furthermore, by knowing the genetic foundation of KC, we can advance customized medicine strategies that modify treatment according to a patient's genetic profile. For example, people with particular genetic markers could benefit from targeted therapy approaches meant to modify the inflammatory processes linked to KC.

It is important to recognize the limitations of this study. Firstly, corneal epithelial mapping is not included in the Pentacam equipment used in this study, despite the fact that epithelial thickness mapping is essential for KC diagnosis, especially in the early stages. Due to this restriction, occasional false-positive KC diagnoses may result, particularly in situations with anterior basement membrane dystrophy. Other than that, despite its large sample size, the study may not adequately reflect the genetic diversity of the Iranian population as a whole. Additional studies using greater and more varied cohorts are required to corroborate these results and investigate the mechanisms at work. Lastly, bias may have been introduced since environmental variables that are known to affect the development of KC, including wearing contact lenses and scratching one's eyes, were self-reported rather than objectively evaluated.

Conclusion

KC is a multifactorial disease involving both genetic and environmental components. Understanding the underlying genetic factors related to the predisposition or susceptibility can play an important role in the pathogenesis, prognostic evaluation and management. The current study demonstrates a strong association between IL-1β and IL-1RN variants and exposure to KC in Iran.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Abbreviations

KC:

Keratoconus

IL -1:

Interleukin-1

MMP:

Matrix metallopeptidase

LD:

Linkage disequilibrium

PCR:

Polymerase chain reaction

RFLP:

Restriction fragment length polymorphism

SNP:

Single nucleotide polymorphism

OR:

Odds ratio

CI:

Confidence interval

SD:

Standard deviation

HWE:

Hardy–Weinberg equilibrium

LOD:

Logarithm of the odds

References

  1. Claessens JL, Godefrooij DA, Vink G, Frank LE, Wisse RPL. Nationwide epidemiological approach to identify associations between keratoconus and immune-mediated diseases. Br J Ophthalmol. 2022;106(10):1350–4.

    Article  PubMed  Google Scholar 

  2. Ambrósio R Jr, Machado AP, Leão E, Lyra JMG, Salomão MQ, Esporcatte LGP, et al. Optimized artificial intelligence for enhanced ectasia detection using Scheimpflug-based corneal tomography and biomechanical data. Am J Ophthalmol. 2023;251:126–42.

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  4. Gokhale NS. Epidemiology of keratoconus. Indian J Ophthalmol. 2013;61(8):382.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Woodward MA, Blachley TS, Stein JD. The association between sociodemographic factors, common systemic diseases, and keratoconus: an analysis of a nationwide heath care claims database. Ophthalmology. 2016;123(3):457-65.e2.

    Article  PubMed  Google Scholar 

  6. 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(10):1436–41.

    Article  PubMed  Google Scholar 

  7. Munir SZ, Munir WM, Albrecht J. Estimated prevalence of keratoconus in the United States from a large vision insurance database. Eye Contact Lens. 2021;47(9):505–10.

    Article  PubMed  Google Scholar 

  8. Ziaei H, Jafarinasab MR, Javadi MA, Karimian F, Poorsalman H, Mahdavi M, et al. Epidemiology of keratoconus in an Iranian population. Cornea. 2012;31(9):1044–7.

    Article  PubMed  Google Scholar 

  9. Salman A, Darwish T, Ghabra M, Kailani O, Haddeh Y, Askar M, et al. Prevalence of keratoconus in a population-based study in Syria. J Ophthalmol. 2022;2022:6064533.

    Article  PubMed  PubMed Central  Google Scholar 

  10. McMonnies CW. Mechanisms of rubbing-related corneal trauma in keratoconus. Cornea. 2009;28(6):607–15.

    Article  PubMed  Google Scholar 

  11. Steahly LP. Keratoconus following contact lens wear. Ann Ophthalmol. 1978;10(9):1177–9.

    CAS  PubMed  Google Scholar 

  12. Wang YM, Ma L, Lu SY, Chan TCY, Yam JCS, Tang SM, et al. Analysis of multiple genetic loci reveals MPDZ-NF1B rs1324183 as a putative genetic marker for keratoconus. Br J Ophthalmol. 2018;102(12):1736–41.

    Article  PubMed  Google Scholar 

  13. Naderan M, Rajabi MT, Zarrinbakhsh P, Naderan M, Bakhshi A, et al. Association between family history and keratoconus severity. Curr Eye Res. 2016;41(11):1414–8.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  15. Niazi S, Moshirfar M, Alizadeh F, Doroodgar F, Baradaran-Rafii A, Filutowski O, et al. Association of 2 lysyl oxidase gene single nucleotide polymorphisms with keratoconus: a nationwide registration study. Ophthalmol Sci. 2023;3(2):100247.

    Article  PubMed  Google Scholar 

  16. Deshmukh R, Ong ZZ, Rampat R, Alió Del Barrio JL, Barua A, Ang M, et al. Management of keratoconus: an updated review. Front Med (Lausanne). 2023;10:1212314.

    Article  PubMed  Google Scholar 

  17. Bykhovskaya Y, Rabinowitz YS. Update on the genetics of keratoconus. Exp Eye Res. 2021;202:108398.

    Article  CAS  PubMed  Google Scholar 

  18. Wilson SE, He YG, Weng J, Li Q, McDowall AW, Vital M, et al. Epithelial injury induces keratocyte apoptosis: hypothesized role for the interleukin-1 system in the modulation of corneal tissue organization and wound healing. Exp Eye Res. 1996;62(4):325–38.

    Article  CAS  PubMed  Google Scholar 

  19. Behzadi P, Sameer AS, Nissar S, Banday MZ, Gajdács M, García-Perdomo HA, et al. The interleukin-1 (IL-1) superfamily cytokines and their single nucleotide polymorphisms (SNPs). J Immunol Res. 2022;2022:2054431.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Matalia H, Swarup R. Imaging modalities in keratoconus. Indian J Ophthalmol. 2013;61(8):394–400.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Kriszt Á, Losonczy G, Berta A, Takács L. Presence of Fleischer ring and prominent corneal nerves in keratoconus relatives and normal controls. Int J Ophthalmol. 2015;8(5):922–7.

    PubMed  PubMed Central  Google Scholar 

  22. Alió JL, Piñero DP, Alesón A, Teus MA, Barraquer RI, Murta J, et al. Keratoconus-integrated characterization considering anterior corneal aberrations, internal astigmatism, and corneal biomechanics. J Cataract Refract Surg. 2011;37(3):552–68.

    Article  PubMed  Google Scholar 

  23. Ambrósio R Jr, Lopes BT, Faria-Correia F, Salomão MQ, Bühren J, Roberts CJ, et al. Integration of Scheimpflug-based corneal tomography and biomechanical assessments for enhancing ectasia detection. J Refract Surg. 2017;33(7):434–43.

    Article  PubMed  Google Scholar 

  24. Mohammadpour M. Diagnostics in ocular imaging: cornea, retina, glaucoma and orbit. 1st ed. Cham: Springer; 2020.

    Google Scholar 

  25. Augustin VA, Son HS, Baur I, Zhao L, Auffarth GU, Khoramnia R. Detecting subclinical keratoconus by biomechanical analysis in tomographically regular keratoconus fellow eyes. Eur J Ophthalmol. 2022;32(2):815–22.

    Article  Google Scholar 

  26. Kataria P, Padmanabhan P, Gopalakrishnan A, Padmanaban V, Mahadik S, Ambrósio R Jr, et al. Accuracy of Scheimpflug-derived corneal biomechanical and tomographic indices for detecting subclinical and mild keratectasia in a South Asian population. J Cataract Refract Surg. 2019;45(3):328–36.

    Article  PubMed  Google Scholar 

  27. Macé M, Galiacy SD, Erraud A, Mejía JE, Etchevers H, Allouche M, et al. Comparative transcriptome and network biology analyses demonstrate antiproliferative and hyperapoptotic phenotypes in human keratoconus corneas. Invest Ophthalmol Vis Sci. 2011;52(9):6181–91.

    Article  PubMed  Google Scholar 

  28. Kim WJ, Rabinowitz YS, Meisler DM, Wilson SE. Keratocyte apoptosis associated with keratoconus. Exp Eye Res. 1999;69(5):475–81.

    Article  CAS  PubMed  Google Scholar 

  29. Kim SH, Mok JW, Kim HS, Joo CK. Association of −31T> C and −511 C> T polymorphisms in the interleukin 1 beta (IL1β) promoter in Korean keratoconus patients. Mol Vis. 2008;14:2109–16.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Palamar M, Onay H, Ozdemir TR, Arslan E, Egrilmez S, Ozkinay F, et al. Relationship between IL1β-511C> T and ILRN VNTR polymorphisms and keratoconus. Cornea. 2014;33(2):145–7.

    Article  PubMed  Google Scholar 

  31. Tamassia N, Castellucci M, Rossato M, Gasperini S, Bosisio D, Giacomelli M, et al. Uncovering an IL-10-dependent NF-KB recruitment to the IL-1ra promoter that is impaired in STAT3 functionally defective patients. FASEB J. 2010;24(5):1365–75.

    Article  CAS  PubMed  Google Scholar 

  32. Khazim K, Azulay EE, Kristal B, Cohen I. Interleukin 1 gene polymorphism and susceptibility to disease. Immunol Rev. 2018;281(1):40–56.

    Article  CAS  PubMed  Google Scholar 

  33. Hall SK, Perregaux DG, Gabel CA, Woodworth T, Durham LK, Huizinga TW, et al. Correlation of polymorphic variation in the promoter region of the interleukin-1β gene with secretion of interleukin-1β protein. Arthritis Rheum. 2004;50(6):1976–83.

    Article  CAS  PubMed  Google Scholar 

  34. Wang Y, Wei W, Zhang C, Zhang X, Liu M, Zhu X, et al. Association of interleukin-1 gene single nucleotide polymorphisms with keratoconus in Chinese Han population. Curr Eye Res. 2016;41(5):630–5.

    CAS  PubMed  Google Scholar 

  35. Chen H, Wilkins LM, Aziz N, Cannings C, Wyllie DH, Bingle C, et al. Single nucleotide polymorphisms in the human interleukin-1β gene affect transcription according to haplotype context. Hum Mol Genet. 2006;15(4):519–29.

    Article  CAS  PubMed  Google Scholar 

  36. Weng J, Mohan RR, Li Q, Wilson SE. IL-1 upregulates keratinocyte growth factor and hepatocyte growth factor mRNA and protein production by cultured stromal fibroblast cells: interleukin-1β expression in the cornea. Cornea. 1997;16(4):465–71.

    Article  CAS  PubMed  Google Scholar 

  37. Mier-Cabrera J, Cruz-Orozco O, de la Jara-Díaz J, Galicia-Castillo O, Buenrostro-Jáuregui M, Parra-Carriedo A, et al. Polymorphisms of TNF-alpha (−308), IL-1beta (+3954) and IL1-Ra (VNTR) are associated to severe stage of endometriosis in Mexican women: a case control study. BMC Womens Health. 2022;22(1):356.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Song X, Wang D, Ben B, Xiao C, Bai L, Xiao H, et al. Association between interleukin gene polymorphisms and susceptibility to gastric cancer in the Qinghai population. J Int Med Res. 2021;49(5):3000605211004755.

    Article  CAS  PubMed  Google Scholar 

  39. Na KS, Mok JW, Kim JY, Joo CK. Proinflammatory gene polymorphisms are potentially associated with Korean non-Sjogren dry eye patients. Mol Vis. 2011;17:2818–23.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors acknowledge Mr. Hamidreza Nematy, B.Sc., for technical editing and writing assistance in the preparation of the manuscript.

Funding

No external funding was received for this research.

Author information

Authors and Affiliations

Authors

Contributions

FD: conception, design, data acquisition, data interpretation, writing draft, revision; FA: conception, data acquisition, data analysis, data interpretation, writing draft, revision; SN: conception, design, data acquisition, data analysis, data interpretation, writing draft, revision; MR: data analysis; NJ: data analysis; AA: data analysis; FN: data acquisition, revision; MAJ: data interpretation, revision; JAB: data interpretation, revision; SD: data analysis; MM: data interpretation, revision; ZG: data interpretation, revision; RA: data interpretation, revision; JLA: conception, data interpretation, writing draft, revision. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jorge L. Alio.

Ethics declarations

Ethics approval and consent to participate

The study protocol for Ethics at Tehran University of Medical Sciences (IR.SBMU.ORC.REC.1399.015) was approved by the committees, and written informed consent was obtained from participants in accordance with the Declaration of Helsinki.

Consent for publication

All authors and participants have provided consent for the publication of individual person's data.

Competing interests

Jorge L. Alió is the associate editor of the journal. The other authors declare that they have no competing interests.

Supplementary Information

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

Doroodgar, F., Alizadeh, F., Niazi, S. et al. Inflammatory and genomic interactions within keratoconus susceptible patients: a nationwide registered case–control study. Eye and Vis 11, 40 (2024). https://doi.org/10.1186/s40662-024-00407-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40662-024-00407-z

Keywords