Skip to main content

Evaluating the association between MPDZ-NF1B rs1324183 and keratoconus in an independent northwestern Chinese population

Abstract

Background

Keratoconus (KC) is a complex, non-inflammatory corneal degenerative disease. Although numerous studies have analyzed the correlation of SNP rs1324183, which located in MPDZ-NF1B gene, and KC in different populations, only few findings were repeated. In this study, to evaluate the association between rs1324183 and KC in a new independent Chinese population, we performed a replication study of the significantly associated rs1324183.

Methods

In total of 114 unrelated KC patients and 88 unrelated controls were recruited from Ningxia, China. We detected the genotypes and alleles of rs1324183 using PCR technology and Sanger sequencing and also analyzed the association between this locus and KC, its clinical parameters by statistical methods.

Results

The frequency of genotype AA (11, 9.6%) and genotypes containing allele A (47, 41.2%) of rs1324183 in KC were both higher than those of the control group. And genotype AA of rs1324183 conferred a higher risk of KC (OR > 1). Moreover, corneal parameter Belin/Ambrósio enhanced ectasia display final D value (BAD-D) had significant correlation (p = 0.002) with AA genotype of rs1324183 in KC.

Conclusions

Our replication study indicates that the results of rs1324183 associated with KC in our population is robust and further better illustrates the significance of BAD-D as a diagnostic indicator for KC. rs1324183 should be considered as the first genetic mark of KC risk in its future diagnosis.

Peer Review reports

Background

Keratoconus (KC) is a complex, non-inflammatory corneal degenerative disease characterized by corneal dilatation, cone-shaped protrusion, irregular myopia, and astigmatism, which usually affects on both eyes, but the progress of the eyes is asymmetrical [1,2,3]. KC occurs more often in adolescence, and severe visual impairment occurs with the thinning of the corneal stromal layer of cornea. It still remains one of the most frequent causes of corneal transplantation [4,5,6].

Researchers have gradually discovered that KC is no longer a rare disease. Some areas of the world, due to genetic and environmental factors, are suffering from an epidemic outbreak of KC. A review of early studies of the prevalence of KC found a range of 50–230 cases per 100,000 (0.05–0.23%) populations in Western countries [3]. However, several recent independent studies have found the prevalence to be higher. In a India study, they estimated it to be 2.3% in general population (and higher in refractive surgery cases) [7, 8]. A longitudinal cohort study from Australia found the prevalence to be as high as 1.2% among 20-year-olds [9]. The first study on population-based KC prevalence in a randomized sample in Turkey showed 2.4% (2393/100000) [10]. While a systematic Review and Meta-Analysis from 15 countries calculated a prevalence of 0.13% in the whole population [11]. Therefore, estimates for the prevalence of keratoconus have varied and it also has not been accurately determined. Different ethnicity living in the different geographic location suggests a genetic basis as well as environmental impact of disease. For example, higher prevalence than the European average has been found in Asian [12,13,14,15,16,17,18,19,20]. According to two genome-wide association studies of keratoconus, as of March 2021, approximately 66 variants from 50 genes were identified to be significantly associated with KC in different populations [21, 22]. However, the relationship between these variants and the pathogenesis of keratoconus is still unclear. Especially, whether the variants can really affect the occurrence of keratoconus and can be used as genetic markers for the early diagnosis of keratoconus remain to be further studied.

In recently years, Wang YM and Hao XD et al. both analyzed the association of SNP locus with keratoconus in Chinese population, and revealed that the rs1324183 (MPDZ-NF1B), which located between the MPDZ and NF1B genes, is a putative genetic marker for monitoring the progression of keratoconus and should be further investigated in other different Chinese populations [23, 24]. Thence, in this study, we collected a total of 114 KC patients who live in northwestern of China for the first time and further confirmed the association between the rs1324183 (MPDZ-NF1B) and KC, providing a clear basis for early clinical diagnosis in Chinese population.

Methods

Subjects

A total of 114 sporadic Chinese KC patients and 88 unrelated controls from Ningxia Hui Autonomous region were recruited from the Ningxia Eye Hospital from Nov. 2017 to Nov.2019. The diagnostic criteria of KC is based on the consensus on the diagnosis of keratoconus in China, which take in and learn from the international grading method of keratoconus [25,26,27,28] as well as grading method of keratoconus of Lixin Xie and Weiyun Shi [29]. Patients with keratoconus have a history of myopia and astigmatism, and the best corrected visual acuity (BCVA) is less than 1.0. At least one of the slit lamp microscopy will be positive, such as corneal stroma thinning, tapered lordosis, Fleischer ring, Vogt’s striae, epithelial or subepithelial scar. The corneal topography shows that the central diopter of the anterior surface of the cornea is greater than 47.00 D. Therefore, the consensus on the diagnosis of keratoconus in China is divided into 4 stages of incubation period, initial period, completion period, and scarring period. The corneal topography of the affected eye during the incubation period is normal, the corrected visual acuity is greater than or equal to 1.0. The initial patient’s best corrected visual acuity (BCVA) is greater than or equal to 0.8, but the corneal topography examination has keratoconus manifestations (central diopter of the anterior surface of the cornea is greater than 47 D). The BCVA of the affected eye during the completion period is less than 0.8, and they have at least one typical clinical symptom of KC, such as corneal stroma thinning, tapered lordosis, Fleischer ring, and Vogt striae. The affected eye in the scarring period have residual scar on the entire cornea after the acute keratoconus edema subsiding.

In this study, each KC patient has at least one eye with typical clinical characteristics from the completion period, with a BCVA < 0.8. The study protocol was approved and reviewed by the Ethics Committee on Human Research at People Hospital of Ningxia Hui Autonomous Region. Written informed consent was received from each participant or his/her legal guardians before participation, and study adhered to the tenets of the declaration of Helsinki.

All standard ophthalmic examination of participants were performed by comprehensive refractometry (VT-10, TOPCON, Japan), slit-lamp biomicroscopy (Topcon, Tokyo, Japan), intraocular pressure (IOP), anterior segment analysis system (Pentacam 70,700, Germany), and corneal biomechanics analyzer (Corvis ST 72,100, Germany). Slit lamp biomicroscopy was used to identify stromal corneal thinning, Vogt’s striae, or a Fleischer ring. Best-corrected visual acuity (BCVA) was performed by comprehensive optometry. Anterior segment analysis system measured 9 corneal parameters of the thinnest point thickness (TP), max keratometry (K max), Belin/Ambrósio enhanced ectasia display final D value (BAD-D), deviation of normality of the front elevation (Df), deviation of normality of the back elevation (Db), deviation of normality of pachymetric progression (Dp), deviation of normality of corneal thinnest point (Dt), deviation of normality of relational thickness (Da), and Ambrósio’s relational thickness (ARTH). And 6 other corneal parameters of corvis biomechanical index (CBI), deformation amplitude ratio (Da_Ratio), adjusted AP1-bIOP (SPA1), biomechanical intra occular pressure (bIOP), integrated radius, and tomographic and biomechanical index (TBI) were recorded by corneal biomechanics analyzer. All results of corneal parameters will generate the corresponding quality factor (QS). When QS > 95%, an “OK” display appears on the test instrument, indicating that the quality of the test data is acceptable. If the test quality is not good enough, they must be retested. In order to avoid detection errors, the operations mentioned above should be effectively checked by the same experienced medical technician for at least three times. The higher quality results will be selected into the test group.

The control subjects were recruited from the patients who were scheduled to undergo laser surgery for myopia and volunteers. Among them, only one eye was randomly selected for analysis in each control group. Except for astigmatism and myopia, there are no any other ophthalmological diseases. The corneal topography examination of both eyes were normal, and the best corrected visual acuity were great than or equal to 1.0.

rs1324183 (MPDZ-NF1B) detection and genotyping

rs1324183 in MPDZ-NF1B was reported to be associated with keratoconus [23, 24, 30, 31]. Genomic DNA was extracted from peripheral blood of all participants using the QIAamp DNA Mini Kit (Qiagen, Germany) according to the manufacturer’s protocol. DNA concentration was inspected by Nanodrop 2000 (Thermo Fisher, USA). The target sequence which contained the rs1324183 site (Primers, forward: 5’-TCCTACCAGCTTGTCTCCAAA-3’, reverse: 5’-ACAAGAAGCCACAAGTCTGGC-3’) were amplified by ordinary PCR (Allsheng, Hangzhou, China) and then conducted Sanger sequencing on ABI 3730 analyzer (Applied, Biosystem). Sites of variation were identified through a comparison of DNA sequences with the corresponding GenBank reference sequences using the Mutation Surveyor software, version 5.0.0. Genotyping (Heterozygous mutation, homozygous mutation and wild-type) were classified based on the results of Sanger sequencing.

Statistical analysis

All statistical analyses were performed using SPSS 22.0 (IBM Corporation, Armonk). Statistical significance was declared at a = 0.05. Kolmogorov–Smirnov test (K-S test) and the Levene test were used respectively to test whether the data conformed to the normal distribution and to evaluate the homogeneity of variance of the data between groups. The Hardy–Weinberg equilibrium (HWE) p value and minor allele frequency (MAF) were estimated using Haploview software, version 4.2 (http://www.broad.mit.edu/mpg/haploview/). We also assessed the odds risk (OR) of susceptibility locus for keratoconus carrying specific genotypes based on the binomial Logistic regression analysis. In order to investigate the most sensitive and specific indicators of Pentacam and Corvis ST in the diagnosis of keratoconus, we drew the receiver operating characteristic curve (ROC) of corneal parameters and then compared the areas under ROC (AUC) to obtain the best diagnostic parameters of KC. If the AUC is between 0.5 and 0.7, the diagnosis efficiency is low. An area between 0.7–0.9 indicates medium diagnostic efficiency. More than 0.9 indicates high diagnostic efficiency. Moreover, in keratoconus group, variance trend analysis was used to assess the difference of average values of 15 corneal parameters (K max, BAD-D, Df, Db, Dp, Dt, Da, ARTH, CBI, Da_Ratio, bIOP, TP, integrated radius, TBI, and SPA1) in different genotypes of rs1324183 (MPDZ-NF1B).

Results

Association of gender, age, and MPDZ-NF1B rs1324183 with keratoconus

The rs1324183, which locates in the intergenic region of MPDZ and NFIB, is shown in Fig. 1. Hardy–Weinberg equilibrium p values and the minor allele frequency of the studied SNP (rs1324183) are shown in Table 1. rs1324183 was in Hardy–Weinberg equilibrium in both patient and control groups, respectively. And they all had a MAF greater than 0.05 in our population. Besides, Basic characteristics of our participants were further analyzed in Table 1, and we found no significant differences in gender and age between KC cases and controls. However, the genotype distributions of rs1324183 (MPDZ-NFIB) (p = 0.032) were significantly different between the patients and controls (Table 2). Especially, the frequency of genotype AA (11, 9.6%) and genotypes containing allele A (47, 41.2%) of rs1324183 in KC were both higher than those of the control group. Genotype AA of rs1324183 (MPDZ-NFIB) had an OR of 3.475 (95% CI = 1.241–8.431, p = 0.003) under the recessive model and an OR of 2.108 (95% CI = 0.464–1.085, p = 0.001) under the dominant model (Table 2). These results confirmed that the genotype AA of rs1324183 (MPDZ-NFIB) was an important KC risk genotype in our population.

Fig. 1
figure 1

Localization of MPDZ-NFIB rs1324183 in the chromosome 9. rs1324183 locates in the intergenic region of MPDZ and NFIB in the 9p23 band of chromosome 9. It has four alleles in this gene locus, alleles: A > C,G,T. A and C are the most common alleles, while G and T are less common

Table 1 Basic characteristics of KC patients and controls
Table 2 Genotype frequencies of studied SNP and their association with susceptibility to KC

The diagnostic efficacy of corneal parameters for keratoconus

To evaluate the diagnostic efficacy of corneal parameters for keratoconus, we further analyzed the correlation between corneal parameters and keratoconus and found that, compared with the control, the KC cases had relatively higher mean db, dp, df, da, dt, BAD-D, CBI, TBI, K max, integrate radius, and Da_Ratio, while the mean of TP, bIOP, SPA1, and ARTH were relatively lower (Table 3). It suggests that all corneal parameters mentioned above may be indicators of KC diagnosis. Then, we studied the sensitivity and specificity of 15 corneal parameters in the diagnosis of KC and found that all the areas under receiver operating characteristic curves (AUC) of corneal parameters except bIOP were greater than 0.9, indicating that they have a high diagnostic efficiency for KC (Table 4). In addition, we found that the first six AUC of corneal parameters in our population were TBI (AUC = 1.000), BAD-D (AUC = 0.991), CBI (AUC = 0.986), Da (AUC = 0.978), Tp (AUC = 0.971), and Da_Ratio (0.971) (Fig. 2). Although the areas under receiver operating characteristic curves of these six corneal parameters appeared to be greater than 0.9, the youden index and specificity of TBI, BAD-D, and CBI were significantly higher than the other three parameters (Table 4). It suggests that TBI, BAD-D, and CBI are more reliable in the diagnosis of keratoconus.

Table 3 Analysis of corneal parameters in keratoconus and control group
Table 4 Analysis of sensitivity and specificity of 15 corneal parameters in the diagnosis of KC
Fig. 2
figure 2

Comparison of the area under receiver operating characteristic curve (AUC) of the first six corneal parameters. Receiver operating characteristic curve (ROC) of TBI, BAD-D, CBI, Da, Tp, and Da_Ratio. Blue line represents TBI, Green line represents BAD-D, Orange line represents CBI, Yellow line represents Da, Pink line represents Tp, and Red line represents Da_Ratio. TBI, BAD-D, and CBI have larger AUC than the other three parameters

Correlation between the corneal parameters and rs1324183 (MPDZ-NFIB) in keratoconus group

To investigate whether the risk genotype of keratoconus is related to its reliable diagnostic parameters, we analyzed the correlation between the corneal parameters and rs1324183 (MPDZ-NFIB) in our KC cases and found that the mean db, dp, da, BAD-D, K max, and integrate radius in KC patients homozygous for the risk allele A at SNP rs1324183 (26.549 ± 34.398 D, 49.091 ± 103.675 D, 3.934 ± 0.451 D, 24.948 ± 29.538 D, 70.636 ± 17.694 D, and 15.609 ± 7.676 D respectively) were significantly higher than patients with the genotype CA (9.961 ± 20.732 D, 16.975 ± 58.199 D, 2.778 ± 0.686 D, 10.022 ± 17.577 D, 54.147 ± 8.427 D, and 10.911 ± 2.718 D, respectively) and CC (7.583 ± 7.044 D, 7.481 ± 6.289D, 2.718 ± 0.794 D, 7.781 ± 5.226 D, 53.521 ± 8.753 D, and 10.479 ± 2.804 D, respectively) (Table 5). In contrast, patients with genotype AA at SNP rs1324183 had relatively lower mean Tp (312.818 ± 118.996 D), bIOP (10.500 ± 5.856 D), SPA1 (45.655 ± 19.401 D), and ARTH (63.664 ± 54.674 D) compared with the patients with CA (447.083 ± 68.893 D, 15.031 ± 3.341 D,69.686 ± 17.856 D, and 221.394 ± 122.248 D, respectively) and CC (456.418 ± 54.807 D, 15.418 ± 2.589 D, 73.057 ± 17.736 D, 242.236 ± 124.917 D) (Table 5). In addition, no other corneal parameters attained significant difference when comparing among different genotype groups. It showed that the KC risk genotype (genotype AA of rs1324183) was significantly associated with 10 specific corneal parameters of db, dp, da, BAD-D, K max, integrate radius, TP, bIOP, SPA1, and ARTH, suggesting that these 10 specific corneal parameters may play an important role in predicting the risk of KC. What’s more, one result actually echoes the findings in Tables 3 and 4, we found that only BAD-D showed a stable and strong association with KC risk in our population, indicating that the BAD-D should be confirmed in more future studies.

Table 5 Correlation of rs1324183 (MPDZ-NF1B) with corneal parameters in keratoconus

Discussion

Since the early clinical manifestations of keratoconus are not obvious, and the overlap and low sensitivity of clinically specific indicators often happens in some keratoconus and some normal corneas with abnormal morphology or thin thickness, which has been a difficulty for ophthalmologists in early diagnosis of keratoconus. Although a growing number of researchers have identified some genetic markers associated with keratoconus [21, 23, 24], the reliability of most results has yet to be widely replicated and confirmed.

In this study, we assessed the replication correlation between SNP rs1324183 and keratoconus in our northwestern population of China, which have never been confirmed in such population. SNP rs1324183 (MPDZ-NFIB, chr9:13,557,491) was first screened associated with the central corneal thickness (CCT) in an Asian cohort in 2012 [32], and then was studied in different KC population in the world [23, 24, 30, 31, 33]. They found that SNP rs1324183 was significantly associated with KC (p < 0.05), conferring a high risk towards KC, and should be used as genetic markers for the diagnosis of keratoconus in the future. In our study, similarly, we also obtained the close correlation between SNP rs1324183 and KC, in particular, the allele A and the genotype AA of rs1324183 (MPDZ-NFIB) were both confirmed as higher risks for keratoconus in our population (Tables 1 and 2). However, we did not find any differences in gender and age between our KC cases and controls (Table 1). Results concerning gender and age preponderance on keratoconus vary between previous studies. Most of them demonstrated that gender and age both associated with high prevalence of KC [34,35,36,37,38,39,40], but some people reported that there was no gender preponderance in KC and also no significant association between age and the prevalence of KC [11, 41]. In our study, no difference in gender and age between KC cases and controls could be because of the small number of subjects.

In addition, we also evaluated the diagnostic efficacy of visualized corneal biomechanics analysis system combined with Pentacam related parameters in our KC patients, 15 corneal parameters (Tp, df, db, dp, dt, da, BAD-D, K max, ARTH, CBI, TBI, integrate radius, Da_Ratio, bIOP, and SPA1) of Pentacam and Corvis ST were included and analyzed in the study. According to reports, even if no corneal morphological abnormality is observed in the early stage of keratoconus, it generally exist corneal biologic and biomechanical abnormalities [42, 43]. Therefore, we saw assessing the diagnostic efficacy of corneal biologic and biomechanical parameters for KC as vital. In the present study, we found that all corneal biologic and biomechanical parameters could differentiate the KC patient from the controls in our population (Table 3). But the sensitivity and specificity of their diagnosis for KC were different. Further study showed that the areas under the receiver operating characteristic curve (AUC) as well as the Youden index of TBI, CBI, and BAD-D were greater than those of the other 12 corneal parameters, indicating that TBI, CBI, and BAD-D all have higher diagnostic efficiency for KC (Table 4, Fig. 2). Similar findings on them were also reported by Hashemi and Salomão MQ Jr et al. [44, 45].

Notably, many studies have reported that the SNP rs1324183 (MPDZ-NFIB) is a genetic marker highly correlated with KC cases [23, 24, 30, 31, 33], and in our population, we have also further confirmed that the genotype AA of rs1324183 confers a higher risk for KC (Table 2). However, few researchers analyzed the correlation between rs1324183 (MPDZ-NFIB) and kinds of corneal biometric and biomechanical parameters in the KC cases. In this study, we explored the correlation between each corneal parameter and rs1324183 (MPDZ-NFIB) in KC patients, and found that BAD-D was the best corneal biometric parameter for diagnosing clinical KC cases. As several studies have shown BAD-D to be a strong parameter to differentiate both keratoconus and subclinical keratoconus from normal corneas [46, 47], but few studies have been conducted on its relationship with genetic markers. Therefore, this is a new evidence to support BAD-D as an effective diagnostic indicator of keratoconus. BAD-D is a comprehensive indicator of anterior segment analysis system (Pentacam), which was calculated based on the linear regression analysis of the tomographic index (the height and thickness of the anterior and posterior surface of the cornea). CBI is a comprehensive biomechanical parameter in the dynamic corneal response parameters measured by Corvis ST. And calculating TBI through BAD-D and CBI can better distinguish potential corneal dilatation. It means that the two corneal parameters, BAD-D and CBI, have a comprehensive effect on TBI. In this study, the result of TBI may be greatly affected by CBI when analyzing their correlation with the genotypes of rs1324183 in KC, and the weak correlation between CBI and KC risk genotype (genotype AA of rs1324183) is likely to weaken the correlation between TBI and KC risk genotype strongly.

All in all, although we have observed that TBI, BAD-D, and CBI, which are the top three in diagnostic efficacy, can well distinguish KC from normal people (Table 4, Fig. 2), only BAD-D is significantly associated with high-risk genotype AA of rs1324183 of KC patients. It suggests that the diagnosis criterion of KC is complex and diverse. The published results still need to be verified widely.

In this case–control study, the KC patients with typical clinical characteristics were mainly from the completion period, with a BCVA < 0.8. And we were unable to compare the correlation between the SNP rs1324183 (MPDZ-NFIB) and different degrees of keratoconus for lacking of sufficient number in patients from incubation period, initial period, and scarring period. While our results showed a potential of BAD-D in predicting the risk of keratoconus, further studies in multiple populations involving different severities of keratoconus are warranted.

Conclusions

We were able to replicate a correlation of rs1324183 (MPDZ-NF1B) with KC in our population. Our results further better illustrates the significance of BAD-D as a diagnostic indicator for KC, and the KC risk genotype, genotype AA of SNP rs1324183, would also probably be a reliable early diagnosis indicator in populations of China, and they both should be noted carefully in its future studies.

Availability of data and materials

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

Abbreviations

KC:

Keratoconus

TP:

The thinnest point thickness

K max:

Max keratometry

BAD-D:

Belin/Ambrósio enhanced ectasia display final D value

Df:

Deviation of normality of the front elevation

Db:

Deviation of normality of the back elevation

Dp:

Deviation of normality of pachymetric progression

Dt:

Deviation of normality of corneal thinnest point

Da:

Deviation of normality of relational thickness

ARTH:

Ambrósio’s relational thickness

CBI:

Corvis biomechanical index

Da_Ratio:

Deformation amplitude ratio

SPA1:

Adjusted AP1-bIOP

bIOP:

Biomechanical intra occular pressure; integrated radius

TBI:

Tomographic and biomechanical index

AUC:

Area Under Curve

ROC:

Receiver Operating Characteristics curve

OR:

Odds Risk

BCVA:

Best-corrected visual acuity

IOP:

Intraocular pressure

References

  1. Jhanji V, Sharma N, Vajpayee RB. Management of keratoconus: current scenario. Br J Ophthalmol. 2011;95(8):1044–50.

    Article  PubMed  Google Scholar 

  2. Krachmer JH, Feder RS, Belin MW. Keratoconus and related noninflammatory corneal thinning disorders. Surv Ophthalmol. 1984;28(4):293–322.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  4. Wilson SE, Klyce SD. Screening for corneal topographic abnormalities before refractive surgery. Ophthalmology. 1994;101(1):147–52.

    Article  CAS  PubMed  Google Scholar 

  5. Jun AS, Cope L, Speck C, Feng X, Lee S, Meng H, et al. Subnormal cytokine profile in the tear fluid of keratoconus patients. PLoS One. 2011;6(1):e16437.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Sarezky D, Orlin SE, Pan W, VanderBeek BL. Trends in Corneal Transplantation in Keratoconus. Cornea. 2017;36(2):131–7.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Jonas JB, Nangia V, Matin A, Kulkarni M, Bhojwani K. Prevalence and associations of keratoconus in rural maharashtra in central India: the central India eye and medical study. Am J Ophthalmol. 2009;148(5):760–5.

    Article  PubMed  Google Scholar 

  8. Shilpy N, Shah Z, Singh S, Purohit D. Prevalence of Keratoconus in Refractive Surgery Cases in Western India. Middle East Afr J Ophthalmol. 2020;27(3):156–9.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Chan E, Chong EW, Lingham G, Stevenson LJ, Sanfilippo PG, Hewitt AW, et al. Prevalence of Keratoconus Based on Scheimpflug Imaging: The Raine Study. Ophthalmology. 2021;128(4):515–21.

    Article  PubMed  Google Scholar 

  10. Özalp O, Atalay E, Yıldırım N. Prevalence and risk factors for keratoconus in a university-based population in Turkey. J Cataract Refract Surg. 2021;47(12):1524–9.

    Article  PubMed  Google Scholar 

  11. 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(2):263–70.

    Article  PubMed  Google Scholar 

  12. 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 

  13. Matthaei M, Sandhaeger H, Hermel M, Adler W, Jun AS, Cursiefen C, et al. Changing Indications in Penetrating Keratoplasty: A Systematic Review of 34 Years of Global Reporting. Transplantation. 2017;101(6):1387–99.

    Article  PubMed  Google Scholar 

  14. Waked N, Fayad AM, Fadlallah A, El Rami H. Keratoconus screening in a Lebanese students’ population. J Francais d’Ophtalmologie. 2012;35(1):23–9.

    Article  CAS  Google Scholar 

  15. Millodot M, Shneor E, Albou S, Atlani E, Gordon-Shaag A. Prevalence and associated factors of keratoconus in Jerusalem: a cross-sectional study. Ophthalmic Epidemiol. 2011;18(2):91–7.

    Article  PubMed  Google Scholar 

  16. Hashemi H, Beiranvand A, Khabazkhoob M, Asgari S, Emamian MH, Shariati M, et al. Prevalence of keratoconus in a population-based study in Shahroud. Cornea. 2013;32(11):1441–5.

    Article  PubMed  Google Scholar 

  17. Assiri AA, Yousuf BI, Quantock AJ, Murphy PJ. Incidence and severity of keratoconus in Asir province Saudi Arabia. Br J Ophthalmol. 2005;89(11):1403–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Georgiou T, Funnell CL, Cassels-Brown A, O’Conor R. Influence of ethnic origin on the incidence of keratoconus and associated atopic disease in Asians and white patients. Eye. 2004;18(4):379–83.

    Article  CAS  PubMed  Google Scholar 

  19. Pearson AR, Soneji B, Sarvananthan N, Sandford-Smith JH. Does ethnic origin influence the incidence or severity of keratoconus? Eye. 2000;4:625–8.

    Article  Google Scholar 

  20. Xu L, Wang YX, Guo Y, You QS, Jonas JB. Prevalence and associations of steep cornea/keratoconus in greater Beijing. The Beijing Eye Study. PLoS One. 2012;7(7):e39313.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  22. 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(1):266.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Hao XD, Chen P, Chen ZL, Li SX, Wang Y. Evaluating the Association between Keratoconus and Reported Genetic Loci in a Han Chinese Population. Ophthalmic Genet. 2015;36(2):132–6.

    Article  CAS  PubMed  Google Scholar 

  24. 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 

  25. Gomes JA, Tan D, Rapuano CJ, Belin MW, Ambrósio R Jr, Guell JL, et al. Group of Panelists for the Global Delphi Panel of Keratoconus and Ectatic Diseases Global consensus on keratoconus and ectatic diseases. Cornea. 2015;34(4):359–69.

    Article  PubMed  Google Scholar 

  26. Amsler M. Kératocõne classique et kératocône fruste; arguments unitaires [Classic keratocene and crude keratocene; Unitary arguments]. Ophthalmologica. 1946;111(2–3):96–101.

    Article  CAS  PubMed  Google Scholar 

  27. McMahon TT, Robin JB, Scarpulla KM, Putz JL. The spectrum of topography found in keratoconus. CLAO J. 1991;17(3):198–204.

    CAS  PubMed  Google Scholar 

  28. Rabinowitz YS. Videokeratographic indices to aid in screening for keratoconus. J Refract Surg. 1995;11(5):371–9.

    Article  CAS  PubMed  Google Scholar 

  29. Xie L, Shi W. Cornea. Beijing: People’s Medical Press China; 2007. p. 360–8.

    Google Scholar 

  30. Sahebjada S, Schache M, Richardson AJ, Snibson G, MacGregor S, Daniel M, et al. Evaluating the association between keratoconus and the corneal thickness genes in an independent Australian population. Invest Ophthalmol Vis Sci. 2013;54(13):8224–8.

    Article  CAS  PubMed  Google Scholar 

  31. Lu Y, Vitart V, Burdon KP, Khor CC, Bykhovskaya Y, Mirshahi A, 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 

  32. Cornes BK, Khor CC, Nongpiur ME, Xu L, Tay WT, Zheng Y, et al. Identification of four novel variants that influence central corneal thickness in multi-ethnic Asian populations. Hum Mol Genet. 2012;21(2):437–45.

    Article  CAS  PubMed  Google Scholar 

  33. Skorodumova LO, Belodedova AV, Sharova EI, Malyugin BE. Poisk geneticheskikh markerov dlia utochniaiushcheĭ diagnostiki keratokonusa [Search for genetic markers for precise diagnostics of keratoconus]. Biomed Khim. 2019;65(1):9–20.

    Article  CAS  PubMed  Google Scholar 

  34. 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 

  35. Ahmed AS, El-Agha MH, Khaled MO, Shousha SM. The prevalence of keratoconus in children with allergic eye disease in an Egyptian population. Eur J Ophthalmol. 2021;31(4):1571–6.

    Article  PubMed  Google Scholar 

  36. Fatima T, Acharya MC, Mathur U, Barua P. Demographic profile and visual rehabilitation of patients with keratoconus attending contact lens clinic at a tertiary eye care centre. Cont Lens Anterior Eye. 2010;33(1):19–22.

    Article  PubMed  Google Scholar 

  37. Sharma R, Titiyal JS, Prakash G, Sharma N, Tandon R, Vajpayee RB. Clinical profile and risk factors for keratoplasty and development of hydrops in North Indian patients with keratoconus. Cornea. 2009;28(4):367–70.

    Article  PubMed  Google Scholar 

  38. Saini JS, Saroha V, Singh P, Sukhija JS, Jain AK. Keratoconus in Asian eyes at a tertiary eye care facility. Clin Exp Optom. 2004;87:97–101.

    Article  PubMed  Google Scholar 

  39. Agrawal VB. Characteristics of keratoconus patients at a tertiary eye center in India. J Ophthalmic Vis Res. 2011;6(2):87–91.

    PubMed  PubMed Central  Google Scholar 

  40. Sidky MK, Hassanein DH, Eissa SA, Salah YM, Lotfy NM. Prevalence of Subclinical Keratoconus Among Pediatric Egyptian Population with Astigmatism. Clin Ophthalmol. 2020;14:905–13.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Akowuah PK, Kobia-Acquah E, Donkor R, Adjei-Anang J, Ankamah-Lomotey S. Keratoconus in Africa: A systematic review and meta-analysis. Ophthalmic Physiol Opt. 2021;41(4):736–47.

    Article  PubMed  Google Scholar 

  42. Shetty R, Arora V, Jayadev C, Nuijts RM, Kumar M, Puttaiah NK, et al. Repeatability and agreement of three Scheimpflug-based imaging systems for measuring anterior segment parameters in keratoconus. Invest Ophthalmol Vis Sci. 2014;55(8):5263–8.

    Article  PubMed  Google Scholar 

  43. Kara N, Bozkurt E, Baz O, Altinkaynak H, Dundar H, Yuksel K, et al. Corneal biomechanical properties and intraocular pressure measurement in Marfan patients. J Cataract Refract Surg. 2012;38(2):309–14.

    Article  PubMed  Google Scholar 

  44. Hashemi H, Beiranvand A, Yekta A, Maleki A, Yazdani N, Khabazkhoob M. Pentacam top indices for diagnosing subclinical and definite keratoconus. J Curr Ophthalmol. 2016;28(1):21–6.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Salomão MQ, Hofling-Lima AL, Gomes Esporcatte LP, Lopes B, Vinciguerra R, Vinciguerra P, et al. The Role of Corneal Biomechanics for the Evaluation of Ectasia Patients. Int J Environ Res Public Health. 2020;17(6):2113.

    Article  PubMed Central  Google Scholar 

  46. Sedaghat MR, Momeni-Moghaddam H, Ambro’sio R Jr, Heidari HR, Maddah N, Danesh Z, et al. Diagnostic ability of corneal shape and biomechanical parameters for detecting frank keratoconus. Cornea. 2018;37(8):1025–34.

    Article  PubMed  Google Scholar 

  47. Shetty R, Rao H, Khamar P, Sainani K, Vunnava K, Jayadev C, et al. Keratoconus screening indices and their diagnostic ability to distinguish normal from ectatic corneas. Am J Ophthalmol. 2017;181:140–8.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank all patients and controls for their participation.

Funding

This research was supported by the < National Natural Science Foundation of China > under grant < number 82060183 > , the < Natural Science Foundation of Ningxia > under grant < number 2020AAC03349 > and < number 2021AAC03302 > , and The < Key Research and Development Project of Ningxia Hui Autonomous Region > under grant < number 2021BEG02045 > .

Author information

Authors and Affiliations

Authors

Contributions

SQY drafted this manuscript, analyzed and interpreted the patient data. DL collected the patient’s information and analyzed the patient data, XLS diagnosed the patient and made contribution to supervision and final approval, JXL diagnosed the patient and made contribution to designing experiment and data analysis, MJM, SZ, and HPL collected the patient’s information, LJZ, ZM, and BYS performed the clinical examination of patient. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xunlun Sheng or Junxiu Liu.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Human Research at People Hospital of Ningxia Hui Autonomous Region (No. LS [KY] 2017010 and approved by July 3, 2017). Written informed consent was obtained from the patient.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Raw data for the control group.

Additional file 2: Table S2.

Raw data for the KC group.

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

Yuan, S., Li, D., Ma, M. et al. Evaluating the association between MPDZ-NF1B rs1324183 and keratoconus in an independent northwestern Chinese population. BMC Ophthalmol 22, 129 (2022). https://doi.org/10.1186/s12886-022-02359-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12886-022-02359-1

Keywords