Replication of SNP associations with keratoconus in a Czech cohort

Introduction Keratoconus is a relatively frequent disease leading to severe visual impairment. Existing therapies are imperfect and clinical management may benefit from improved understanding of mechanisms leading to this disease. We aim to investigate the replication of 11 single nucleotide polymorphisms (SNPs) with keratoconus. Methods SNPs from loci previously found in association with keratoconus were genotyped in 165 keratoconus cases of Caucasian Czech origin (108 males and 57 females) and 193 population and gender-matched controls. They included rs1536482 (COL5A1), rs4839200 (KCND3), rs757219 and rs214884 (IMMP2L), rs1328083 and rs1328089 (DAOA), rs2721051 (FOXO1), rs4894535 (FNDC3B), rs4954218 (MAP3K19, RAB3GAP1), rs9938149 (ZNF469) and rs1324183 (MPDZ). A case-control association analysis was assessed using Fisher’s exact tests. Results The strongest association was found for rs1324183 (allelic test OR = 1.58; 95% CI, 1.10–2.24, p = 0.01). Statistically significant values were also obtained for rs2721051 (allelic test OR = 1.72; 95% CI, 1.07–2.77, p = 0.025) and rs4954218 (allelic test OR = 1.53; 95% CI, 1.01–2.34; p = 0.047) which showed an opposite effect direction compared to previously reported one. Conclusion Independent replication of association between two SNPs and keratoconus supports the association of these loci with the risks for the disease development, while the effect of rs4954218 warrants further investigation. Understanding the role of the genetic factors involved in keratoconus etiopathogenesis may facilitate development of novel therapies and an early detection.


Introduction
Keratoconus is a relatively frequent disease leading to severe visual impairment. Existing therapies are imperfect and clinical management may benefit from improved understanding of mechanisms leading to this disease. We aim to investigate the replication of 11 single nucleotide polymorphisms (SNPs) with keratoconus.

Conclusion
Independent replication of association between two SNPs and keratoconus supports the association of these loci with the risks for the disease development, while the effect of rs4954218 warrants further investigation. Understanding the role of the genetic factors Introduction Keratoconus (KC) is a progressive eye disoder characterized by non-inflammatory central or paracentral corneal thinning and steepening leading to irregular astigmatism that causes a decrease in visual acuity. In patients with advanced disease there is also significant corneal scarring. The disease is almost always bilateral although there may be marked asymetry. Onset occurs typically in early adolescence and progression stops before the age of 40 years. [1]. KC is a relatively common disorder, affecting approximately 1 in 2,000 people [2]. Although many cases appear sporadic, familial aggregations and twin studies suggest a strong genetic component in the KC development [3][4][5]. Specific genetic factors contributing to the risk of KC remain largely unknown. Two relatively underpowered genome-wide association studies (GWAS) using single nucleotide polymorphisms (SNPs) markers have been performed to date and identified several variants reportedly influencing susceptibility to KC [6,7]. To date, only the association of KC with rs2956540 in LOX and rs3735520 in the promoter of HGF has been independently replicated in individuals of European descent [8]. The LOX gene is also located beneath the linkage peak reported in a genome-wide linkage analysis in KC sib pair families genotyped with microsatellite markers [9].
As corneal thinning is a known genetically correlated endophenotype for KC [10], quantitative explorations of the genetic architecture of KC has attempted through associations with KC of central corneal thickness (CCT) genetic risk loci [10,11]. Li et al. (2012) [7] and Lu et al. (2013) [10] reported five loci (six independent SNPs) significantly associated with keratoconus. Of these, effects of two SNPs rs1324183 upstream of MPDZ and rs9938149 upstream of ZNF469 genes have been confirmed in relation to KC risk in an independent cohort of patients of European origin [12].
The relatively low prevalence of KC has constrained both the power of initial discovery cohorts as well as the possibilities of independent replication; indeed most variants identified did not reach in discovery stages commonly considered genome-wide significance thresholds. Yet independent replications are an integral part of the genetic association studies as, among other reasons, they contribute to the increasing credibility of GWAS, improve the accuracy of the effect size calculation and help in the elimination process of false positive findings. Given the general lack of independent replication for KC loci, additional evidence of genetic association, whether positive or negative, will improve our understanding of genetic factors and mechanisms truly influencing KC in the general population. In the current study we genotyped and analyzed in a Czech case-control panel 11 SNPs previously reported to be associated with the risk of KC development in cases of European descent.

Materials and methods
The study was approved by the Ethics Committee of the General Teaching Hospital and Charles University, Prague. Sample collection and experiments were conducted according to the principles set out in the Declaration of Helsinki. Written, informed consent was obtained from all individuals.
A total of 165 unrelated KC cases, 57 females and 108 males, participated in the study. All cases were reported for care at the Department of Ophthalmology in Prague between 2011 and 2013, and exhibited KC pattern on front sagittal curvature maps, together with localized corneal thinning in at least one eye, or had and an advanced disease with stromal scarring, Fleischer ring and Vogt striae [1]. Only patients with topographic KC index grade 1 or higher, as calculated by Pentacam build-in KC detection software (Oculus Optikgeräte GmbH, Wetzlar, Germany), were included [13]. Bilaterally grafted patients for KC were also considered as cases. 193 unrelated Czech Caucasian individuals (79 females, 114 males) self-reporting no ocular disorders impairing vision were used as controls.
To test for association, Fisher's exact test odds ratios (ORs) and p-values were calculated for each SNP, under an allelic test model which assumes co-dominance, using the 'exact2x2' package in R (https://cran.r-project.org). Secondary analyses of association under assumptions of dominant and recessive models were calculated for each SNP. An α<0.05 was considered as the threshold for statistical significance, on the grounds these variants were pre-selected because of strong prior probabilities of association. Allele frequencies in the control group were compared to European panels in the HapMap project (http://hapmap.ncbi.nlm.nih.gov/). Conventional Sanger sequencing was performed when samples had DNA concentration not reaching requirements for KASP assays and in cases of genotype calling failure. In total, 1.57% genotypes (62 out of 3960) were established using this method. In addition, for internal data quality, all 11 SNPs were Sanger-sequenced in 5 randomly selected samples (primers sequences are listed in S3 Table).

Results
The mean age of the cases was 37.18 years (standard deviation 13.26, range 15-69 years) whereas the controls were in average 39.46 years old at recruitment (standard deviation 13.68, range 20-81 years). The sex ratio was also similar between cases and controls (65% male cases compared to 60% male controls).
All SNPs in controls had allele frequencies similar to publicly available European reference panels and were in Hardy-Weinberg equilibrium (p>0.05).

Discussion
Out of 11 SNPs previously reported to influence the risk of KC in predominantly Caucasian cohorts (S1 and S2 Tables) [7,10,14], we have observed that in the Czech case-control panel rs1324183 (MPDZ locus) confers the strongest and most significant risk compared to all other SNPs included in our genotyping panel. Other variants that significantly contributed to KC genetic risk in our samples were rs2721051 (FOXO1 locus) and rs4954218 (MAP3K19/RAB3-GAP1 locus). For all variants, our results appear to support an additive genetic model since secondary analyses under assumptions of different modes of inheritance (dominant or recessive) failed to substantially improve the strength of association. All studies performed on populations of European descent (S1 and S2 Tables), including the Czech one, demonstrated that the minor allele of rs1324183 (MPZD) serves as a risk factor.
Surprisingly the effect sizes of four SNPs out of six initially associated in the original KC GWAS went in the opposite way in our samples, whereas for those initially associated with CCT (and KC); the effect size consistently went in the same direction.
Contrary to the previous studies testing individuals of European descent, where minor allele of rs4954218 (MAP3K19/RAB3GAP1) showed a protective effect, in the Czech cohort it conferred a statistically significant (p<0.05) risk for KC [7,14]. There are several possible explanations for this inconsistent direction of effect. First, it may be a case of type I or type II error in either study, especially in consideration of the modestly powered samples involved. Second, this may be another instance of the so-called "flip-flop phenomenon" [16] also seen in genetic literature for other ophthalmic diseases [17] presumably due to variations in linkage disequilibrium with the causative alleles in different populations. Future studies, using even larger number and potentially more ethnically diverse KC samples may be needed to clarify definitely the role of the rs4954218 (MAP3K19/RAB3GAP1) in KC pathogenesis.
Initially, it was speculated that association with KC of rs4954218 arises from unidentified causative variants interfering with RAB3GAP1 gene functionality, the latter being a good candidate gene based on its biological function [7]. However, according to the GRCh38 human genome assembly, this SNP is intronic to a novel protein coding gene MAP3K19 (mitogen-activated protein kinase kinase kinase 19). Studies have shown that MAP kinases regulate neuronal differentiation by activating transcription factors such as AP-1 [18,19]. Interestingly, transcriptome analysis of KC corneas revealed that 69 of 87 genes with reduced expression compared to control tissues are regulated by AP-1 [20]. The effect of rs4954218 maybe therefore associated with pathways unrelated to RAB3GAP1. However, until the exact mechanisms between KC pathogenesis and the locus tagged by rs4954218 are known, this remains hypothetical.
The study is limited by the relatively modest sample size, which may be the reason for inability to replicate some of the SNPs in the Czech cases-control panel. However, the original GWAS were based on relatively small number of cases which often leads to type I or II errors. Another limitation is the usage of population controls instead of ophthalmologically examined subjects although given the estimated prevalence of KC 1:2,000 [2] it is not very likely that an individual with an unrecognized KC phenotype would be included. In conclusion, the current replication study has positively validated association of several SNPs with KC. The results of this study raise our confidence in the relevance of association of these genetic loci with KC and should pave the way to further studies that are needed to assess the functional mechanisms and direct biological relevance of the genetic variants behind that association.
Supporting information S1