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

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. Supplementary Information The online version contains supplementary material available at 10.1186/s12886-022-02359-1.

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.

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'-TCC TAC CAG CTT GTC TCC AAA-3' , reverse: 5'-ACA AGA AGC CAC AAG TCT GGC-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/ haplo view/). 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).

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.

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.

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

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