Whole exome sequencing highlights variants in association with Keratoconus in Jordanian families
BMC Medical Genetics volume 21, Article number: 177 (2020)
Keratoconus (KC) is usually bilateral, noninflammatory progressive corneal ectasia in which the cornea becomes progressively thin and conical, resulting in myopia, irregular astigmatism, and corneal scarring.
Eight families characterized by consanguineous marriages and/or multiple keratoconic individuals were examined genetically. Whole exome sequencing was done as trio or quadro per family. The output of the filtration procedure, based on minor allele frequency (MAF) less than 0.01 for homozygous variants and MAF equals 0 for heterozygous variants, is 22 missense variants.
Based on the gene/protein function five candidate variants were highlighted in four families. Two variants were highlighted in one family within the genes MYOF and STX2, and one variant is highlighted in each of the other three families within the genes: COL6A5, ZNF676 and ZNF765.
This study is one of the very rare that highlights genetic variants in association with KC.
Keratoconus (KC) is sometimes bilateral, noninflammatory progressive corneal ectasia during which the cornea becomes progressively thin and conical, leading to myopia, irregular astigmatism, and corneal scarring. Patients with KC have cone shaped cornea (hence the name keratoconus, derived from the Greek word for cornea (‘kerato’) and cone shaped (‘conus’). It usually arises in the teenage years and progresses, eventually it stabilizes in the 3rd/4th decades . The clinical phenotypes of KC are highly variable; however, the common feature is corneal steepening which is normally detected at an early stage of the disease using computer-assisted corneal tomography. The current treatment of keratoconus is cross linking using riboflavin (vitamin B2) and ultraviolet light, which could prevent (stop) the progression of KC. In mild cases of KC, the refractive errors are usually treated by glasses or contact lenses but in more advanced cases, surgery is required in order to restore optimal visual acuity .
While the prevalence of KC in Jordan is unknown it is known that the severity of KC is related to the consanguineous marriages. For example the prevalence (per 100,000) of KC is high in societies with high consanguinity such as India (2300) , Iran (2500) , Lebanon (3330) , and Jerusalem (2340)  and low in societies with low consanguinity such as UK/Caucasian (57) , Denmark (86) , Finland (30) , Japan (17.3) , and USA (54.5) . In Jordan consanguineous marriage ranges between 25.5% in the capital Amman to 52.1% in Irbid (www.consang.net), and therefore KC is estimated to have comparable prevalence rates to those in Lebanon and Jerusalem [11, 12].
Aetiology of KC involves environmental factors such as eye rubbing, atopy, and allergies as well as genetic factors. Because of its heterogeneity, KC is either associated with genetic systemic disorders or exhibited in isolated cases (i.e non-syndromic; Online Mendelian Inheritance In Man, OMIM#148300) . Therefore, genetic factors are expected to cause non-syndromic KC especially in consanguineous families and/or families with multiple affected individuals.
The identification of genes causing non-syndromic KC has been the focus of many studies worldwide using genome wide linkage studies and genome wide association studies. Several variants have been implicated in genes coding for the primary composition of cornea including collagens and the components of the extracellular matrix. These genes include: LOX, CAST, DOCK9, IL1RN, SLC4A11, HGF, RAB3GAP1, TGFBI, ZNF469, ZEB1, VSX1, COL5A1, COL4A3, COL4A4, FNDC3B, FOXO1, MPDZ-NF1B, WNT10A, SOD1, IL1B, IL1A and MIR184. However not all analyses of these genes confirm their role in KC pathogenesis . Mutations in VSX1 and SOD1 were assigned as causative to the isolated cases of KC [15, 16]. However, contradictory results were obtained regarding the pathogenic variants within these genes. While some studies revealed missense mutations associated with KC and thus suggesting an important role of these genes in KC [15, 17,18,19], others found that the mutations do not segregate with the disease or do not prove the pathogenicity of the disease [15, 20]. Moreover, in the era of Next Generation Sequencing (NGS) mutations in these genes were revealed in the genome aggregation database (gnomAD; https://gnomad.broadinstitute.org ), suggesting that the mutations in VSX1 and SOD1 may not contribute to KC. Other studies have shown that mutations in MIR184 (the most abundant expressed microRNA in the corneal and lens epithelia) are candidates (but not consistent) causing KC and cataract [22, 23]. In Jordan, NGS was recently used to identify variants causing neurodevelopment diseases [24,25,26,27,28,29,30,31]. In one study homozygous frameshift variant in the gene GALNT14 is identified to be in association with KC .
In this study, Whole Exome Sequencing (WES) is used to identify rare variants associated with KC in eight Jordanian families.
Patients with KC were recruited in July 2017 from Sight and Insight clinic/Amman in collaboration with Dr. Ammar Hiwari and from Al Karak hospital in collaboration with Dr. Khalid Al Zubi. KC was diagnosed during a routine eye test using slit lamp examination, snellen chart vision test, refraction Corneal topography by a pentacam machine “Oculus, Germany”. KC was confirmed using Belin / Ambrósio Enhanced Ectasia Display. Because the goal of this study is the identification of genetic factors being in association with KC, the targeted study group in this study is consanguineous families and/or multiple affected individuals in the same family.
Written informed consents were obtained from all participants in this study. Written informed consents were obtained from the parents of all participants under the age of 16. All participants have provided written informed consents to publish all identifying images, personal details and clinical details anonymously. Peripheral blood samples were collected from all participants.
Molecular genetic analysis
Total genomic DNA was extracted with FlexiGene DNA kit. WES was conducted to 28 individuals (17 keratoconic and 11 healthy). Exome sequencing, bioinformatics analyses and variants validation by sanger sequencing were done as described previously .
Variants filtration procedure
We identified high-quality variants that are located in the protein coding genes including variants within the two base pair flanking the splicing sites (according to Ensembl-GRCh37.73). We maintained only the variants meeting the following quality criteria: at least 20X coverage and a mapping quality score ≥ 60. Variants meetings quality criteria were grouped into homozygous (recessive model) and heterozygous (dominant model) variants. In house controls of 200 Jordanian exomes were used in the filtration procedure.
Homozygous variants were filtered according to the following. Rare variants with minor allele frequency (MAF) less than 0.01 in the following databases: genome aggregation database (gnomAD), 1000 genome project and in the in-house sequenced controls. Shared homozygous variants between the keratconic individuals per family but heterozygous in the parents (if not keratoconic) were maintained. Filtration was then done only for the loss of function variants (LOF; including stop gain, frameshift, splice site acceptor and splice site donor) and the missense variants predicted to be possibly- or probably-damaging by Polyphen2 Humvar and predicted to be deleterious by SIFT [33,34,35]. Loss of function (LOF) variants were also excluded if the genes they are located in carry other homozygous LOF variants in gnomAD, 1000 genome project or in the in-house sequenced controls. Variants located in genes reported with non-ocular OMIM disease were excluded because we aim to identify rare variants in association with non-syndromic KC (Fig. 1).
Heterozygous variants were filtered according to three criteria. Firstly, variants that are absent in the following databases: gnomAD, 1000 genome project and in the in-house sequenced controls. Secondly, shared heterozygous variants between the keratconic individuals per family but wild type in the parents (if not affected). Filtration was then done only for the loss of function variants (LOF; including stop gain, frameshift, splice site acceptor and splice site donor) and missense variants predicted to be possibly- or probably-damaging by Polyphen2 Humvar and predicted to be deleterious by SIFT. LOF variants were also excluded if the gene they are located in carry other LOF variants in gnomAD, 1000 genome project or in the in-house sequenced controls. Thirdly, variants were excluded if they are in genes reported with non-ocular OMIM disease because we aim to identify rare variants with non-syndromic KC (Fig. 1).
Based on the filtration procedure, 22 variants were identified. All of the identified variants are missense. Four variants are homozygous in four families: KC003, KC004, KC010 and KC011 and 18 variants are heterozygous in five families: KC001(two variants), KC003 (three variants) KC005 (five variants), KC007 (one variants), KC008 (seven variants) (Table 1). None of the filtered variants is in a gene reported with OMIM disease.
The family has 3 keratoconic siblings born to non-consanguineous parents (Fig. 2a). WES was performed to the 3 siblings (II-1, II-2 and II-3) and to their mother (I-1). Two heterozygous missense variants out of 118,763 variants were identified in the genes CDC42BPA and LRRC16B. While no information about the gene LRRC16B, the gene CDC42BPA encodes a kinase protein that phosphorylates histone H1. This protein plays a role in the regulation of cytoskeleton and cell migration and abundant in the heart, brain, skeletal muscle, kidney and pancreas . Despite the fact that the identified variants (c.4586G > A:p.Arg1529Gln in the gene CDC42BPA and c.2290G > A:p.Asp764Asn LRRC16B) in these two genes are predicted to be deleterious no relationship can be established between these two variants and KC because there are no reports link any of these two genes with ocular functions.
Distantly related parents in this family have two keratoconic progeny (Fig. 2b). Performing WES to II-1, II-2 and I-1 yields 113,885 variants. According to the filtration procedure four variants are identified. One variant is homozygous in the gene POLR2M and three are heterozygous in the genes GHITM, MYOF and STX2.
The variants in the genes POLR2M and GHITM are excluded to be in association with keratoconus because the first encodes a subunit of a specific form of RNA polymerase II that act as negative regulator of transcriptional activation by the mediator complex , and the second encodes a growth hormone inducible transmembrane protein that is required for the mitochondrial tubular network and cristae organization . Usually abnormalities in such functions are in association with developmental disorders.
The variants in the genes MYOF and STX2 are more likely to be in association with keratoconus because MYOF is calcium/phospholipid-binding protein that plays a role in the plasmalemma repair mechanism of endothelial cells that permits rapid resealing of membranes disrupted by mechanical stress . While the gene STX2 encodes a protein that regulates epithelial mesenchymal interactions and epithelial cell morphogenesis and activation .
The parents in this family I-1 and I-2 belong to two different families but from the same town nonetheless they have two keratoconic progeny II-1 and II-2 (Fig. 2c). WES revealed 130,125 variants and the filtration identified one variant (c.5381G > A:p.Gly1794Glu) in the gene MLLT4 which is not likely contributing to keratoconus. This is because this gene encodes a protein involved in signalling and organization of cell junctions during embryogenesis . Dysfunction here is most likely in association with central nervous system diseases and therefore, no variant can be assigned to this family in association with keratoconus.
Distantly related parents with four keratoconic individuals including the mother were recruited in this family (Fig. 2d). The whole exome sequences in this family revealed 124,431 variants and the filtration procedure identified five variants in five different genes (Table 1). Screening the functions of the gene products of these five genes excludes any of them to be in association with keratoconus and therefore no variant can be assigned to this family.
KC007 and KC011
While the parents in these two families are first cousin with two keratoconic progeny (Fig. 2e and h), no homozygous variants were identified but rather a heterozygous variant (c.1240A > G:p.Ile414Val) in the gene ZNF676 for family KC007 and the variant c.1133C > T:p.Thr378Ile in the gene ZNF765. No detailed information is available about the proteins of these two genes except being zinc finger proteins that might be involved in transcriptional regulations. However, the identified variants in these two zinc finger proteins might be in association with KC because variants in another zinc finger protein (ZNF469) were previously reported in association with keratoconus [42,43,44,45].
In this family the mother and her daughter are keratoconic (Fig. 2f). Whole exome sequence revealed 103,750 variants and the filtered are seven variants in seven genes. By searching the functions of the proteins encoded by these seven genes, only one is prioritized COL6A5. This gene encodes a member of the collagen VI that acts as a cell-binding protein .
This is the only family in this study with one keratoconic individual who belong to nonconsanguineous parents (Fig. 2g). According to the filtration procedure only one variant was identified in the gene NPIPB5. This gene encodes a nuclear pore complex interacting protein which does not seem to be involved with keratoconus and therefore no variants are assigned to this family.
In summary, WES and the validation by sanger sequencing allowed us to highlight variants in association with KC in four families (two variants in KC003 and one variant in each of KC007, KC008, and KC011) and no variants to highlight in three families (KC001, KC004, KC005 and KC010). The families were chosen in this study with multiple keratoconic individual as a sort of evidence for the genetic contribution on KC. Nonetheless, none of the prioritized variants in the families (KC003, KC007, KC008, and KC011) is strong candidate to cause KC. This is due to the lack of concrete information connecting between the functions of the gene products (proteins) and the cornea. Therefore, further experimental assays are needed to strengthen the possible causality of the identified variants and the candidate genes on KC.
Availability of data and materials
The highlighted variants in association with keratoconus in this article are available in Leiden Open Variation Database LOVD 3.0 (https://databases.lovd.nl/) according to reference genome GRCh37.73 used in this study.
KC003 (MYOF) : https://databases.lovd.nl/shared/individuals/00307910
KC003 (STX2) :https://databases.lovd.nl/shared/individuals/00307911
KC007 (ZNF676) :https://databases.lovd.nl/shared/individuals/00307912
KC008 (COL6A5) :https://databases.lovd.nl/shared/individuals/00307914
KC011 (ZNF765) :https://databases.lovd.nl/shared/individuals/00307913
Minor Allele Frequency
Next Generation Sequencing
Whole Exome Sequencing
Online Mendelian Inheritance In Man
Genome Aggregation Database
loss of function
Genome Reference Consortium Build 37.73
Ihalainen A. Clinical and epidemiological features of keratoconus genetic and external factors in the pathogenesis of the disease. Acta Ophthalmologica Supplement. 1986;178:1–64.
Li X, Rabinowitz YS, Rasheed K, Yang H. Longitudinal study of the normal eyes in unilateral keratoconus patients. Ophthalmology. 2004;111(3):440–6.
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.
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.
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.
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.
Pearson AR, Soneji B, Sarvananthan N, Sandford-Smith JH. Does ethnic origin influence the incidence or severity of keratoconus? Eye (Lond). 2000;14(Pt 4):625–8.
Bak-Nielsen S, Ramlau-Hansen CH, Ivarsen A, Plana-Ripoll O, Hjortdal J. Incidence and prevalence of keratoconus in Denmark - an update. Acta Ophthalmol. 2019;97:752–5.
Tanabe U, Fujiki K, Ogawa A, Ueda S, Kanai A. Prevalence of keratoconus patients in Japan. Nippon Ganka Gakkai Zasshi. 1985;89(3):407–11.
Kennedy RH, Bourne WM, Dyer JA. A 48-year clinical and epidemiologic study of keratoconus. Am J Ophthalmol. 1986;101(3):267–73.
Abu Ameerh MA, Bussieres N, Hamad GI, Al Bdour MD. Topographic characteristics of keratoconus among a sample of Jordanian patients. Int J Ophthalmol. 2014;7(4):714–9.
Hanan A. Hamamy ATM, Azmy M. Al-Hadidy, Kamel M. Ajlouni. Consanguinity and genetic disorders. Profile from Jordan.pdf. Saudi Med J 2007. 2007;28(7):185–92.
Joel Sugar M, Marian S, Macasi M. What causes Keratoconus? Cornea. 2012;31(6):716–9.
Bykhovskaya Y, Margines B, Rabinowitz YS. Genetics in Keratoconus: where are we? Eye Vis (Lond). 2016;3:16.
Dash DP, George S, O'Prey D, Burns D, Nabili S, Donnelly U, et al. Mutational screening of VSX1 in keratoconus patients from the European population. Eye (Lond). 2010;24(6):1085–92.
Udar N, Atilano SR, Small K, Nesburn AB, Kenney MC. SOD1 haplotypes in familial keratoconus. Cornea. 2009;28(8):902–7.
Bisceglia L, Ciaschetti M, De Bonis P, Campo PA, Pizzicoli C, Scala C, et al. VSX1 mutational analysis in a series of Italian patients affected by keratoconus: detection of a novel mutation. Invest Ophthalmol Vis Sci. 2005;46(1):39–45.
Tang YG, Picornell Y, Su X, Li X, Yang H, Rabinowitz YS. Three VSX1 gene mutations, L159M, R166W, and H244R, are not associated with keratoconus. Cornea. 2008;27(2):189–92.
De Bonis P, Laborante A, Pizzicoli C, Stallone R, Barbano R, Longo C, et al. Mutational screening of VSX1, SPARC, SOD1, LOX, and TIMP3 in keratoconus. Mol Vis. 2011;17:2482–94.
Štabuc-Šilih M, Stražišar M, Ravnik Glavač M, Hawlina A, Glavač D. Genetics and clinical characteristics of keratoconus. Acta Dermatoven. 2010;19(2):3–10.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. bioRxiv. 2020;581:434–43.
Hughes AE, Bradley DT, Campbell M, Lechner J, Dash DP, Simpson DA, et al. Mutation altering the miR-184 seed region causes familial keratoconus with cataract. Am J Hum Genet. 2011;89(5):628–33.
Iliff BW, Riazuddin SA, Gottsch JD. A single-base substitution in the seed region of miR-184 causes EDICT syndrome. Invest Ophthalmol Vis Sci. 2012;53(1):348–53.
Froukh T. First record mutations in the genes ASPA and ARSA causing Leukodystrophy in Jordan. Biomed Res Int. 2019;2019:7.
Froukh T. Genetic study in a family with dopa-responsive dystonia revealed a novel mutation in sepiapterin reductase gene. Pak J Med Sci. 2019;35(6):1736–9.
Froukh T, Nafie O, Al Hait SAS, Laugwitz L, Sommerfeld J, Sturm M, et al. Genetic basis of neurodevelopmental disorders in 103 Jordanian families. Clin Genet. 2020;97(4):621–7.
Froukh TJ. Next generation sequencing and genome-wide genotyping identify the genetic causes of intellectual Disbailiuty in ten consanguineous families from Jordan. Tohoku J Exp Med. 2017;243(4):2970309.
Han C, Alkhater R, Froukh T, Minassian AG, Galati M, Liu RH, et al. Epileptic encephalopathy caused by mutations in the guanine nucleotide exchange factor DENND5A. Am J Hum Genet. 2016;99(6):1359–67.
Johansen A, Rosti RO, Musaev D, Sticca E, Harripaul R, Zaki M, et al. Mutations in MBOAT7, encoding Lysophosphatidylinositol Acyltransferase I, Lead to intellectual disability accompanied by epilepsy and autistic features. Am J Hum Genet. 2016;99(4):912–6.
Nafi ORB, Riess O, Buchert R, Froukh T. Two cases of variant late infantile ceroid lipofuscinosis in Jordan. World J Clin Cases. 2019;7(2):203–8.
Reuter MS, Tawamie H, Buchert R, Hosny Gebril O, Froukh T, Thiel C, et al. Diagnostic yield and novel candidate genes by exome sequencing in 152 consanguineous families with neurodevelopmental disorders. JAMA Psychiatry. 2017;74(3):293–9.
Froukh T, Hawwari A. Autosomal recessive non-syndromic Keratoconus: homozygous Frameshift variant in the candidate novel gene GALNT14. Curr Mol Med. 2019;19(9):683–7.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013;Chapter 7:Unit7 20.
Ng PC. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4.
Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC. SIFT missense predictions for genomes. Nat Protoc 2016;11(1):1–9.
Zhao Y, Loyer P, Li H, Valentine V, Kidd V, Kraft AS. Cloning and chromosomal location of a novel member of the myotonic dystrophy family of protein kinases. J Biol Chem. 1997;272(15):10013–20.
Hu X, Malik S, Negroiu CC, Hubbard K, Velalar CN, Hampton B, et al. A mediator-responsive form of metazoan RNA polymerase II. Proc Natl Acad Sci U S A. 2006;103(25):9506–11.
Oka T, Sayano T, Tamai S, Yokota S, Kato H, Fujii G, et al. Identification of a novel protein MICS1 that is involved in maintenance of mitochondrial morphology and apoptotic release of cytochrome c. Mol Biol Cell. 2008;19(6):2597–608.
Britton S, Freeman T, Vafiadaki E, Keers S, Harrison R, Bushby K, et al. The third human FER-1-like protein is highly similar to dysferlin. Genomics. 2000;68(3):313–21.
Low SH, Li X, Miura M, Kudo N, Quinones B, Weimbs T. Syntaxin 2 and endobrevin are required for the terminal step of cytokinesis in mammalian cells. Dev Cell. 2003;4(5):753–9.
Taya S, Yamamoto T, Kano K, Kawano Y, Iwamatsu A, Tsuchiya T, et al. The Ras target AF-6 is a substrate of the fam deubiquitinating enzyme. J Cell Biol. 1998;142(4):1053–62.
Loukovitis E, Tsotridou E, Vakalis T, Asteriadis S, Sousouras T, Zachariadis Z, et al. Contributions of the superoxide dismutase 1 and zinc finger protein 469 (ZNF469) genes to keratoconus. J Biol Regul Homeost Agents. 2019;33(2):457–60.
Yildiz E, Bardak H, Gunay M, Bardak Y, Imamoglu S, Ozbas H, et al. Novel zinc finger protein gene 469 (ZNF469) variants in advanced Keratoconus. Curr Eye Res. 2017;42(10):1396–400.
Yu X, Chen B, Zhang X, Shentu X. Identification of seven novel ZNF469 mutations in keratoconus patients in a Han Chinese population. Mol Vis. 2017;23:296–305.
Vincent AL, Jordan CA, Cadzow MJ, Merriman TR, McGhee CN. Mutations in the zinc finger protein gene, ZNF469, contribute to the pathogenesis of keratoconus. Invest Ophthalmol Vis Sci. 2014;55(9):5629–35.
Gara SK, Grumati P, Urciuolo A, Bonaldo P, Kobbe B, Koch M, et al. Three novel collagen VI chains with high homology to the alpha3 chain. J Biol Chem. 2008;283(16):10658–70.
The authors would like to thank all participating patients in this study, the lab technicians Aya Baraghithi and Tala Issa, and the informatician Anis Nazir.
The deanship of scientific research and graduate studies/Philadelphia University/Jordan has completely funded this study.
Ethics approval and consent to participate
This study was approved by the ethics committee of Philadelphia University in Jordan and conducted according to the principles outlined in the declaration of Helsinki in 1964. Written informed consents were obtained from all participants in this study. Written informed consents were obtained from the parents of all participants under the age of 16.
Consent for publication
All participants have provided written informed consents to publish all identifying images, personal details and clinical details anonymously. Moreover, written informed consents for publication of identifying images or other personal or clinical details were obtained from the parents or legal guardians of any participant under the age of 18.
The authors declare no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Froukh, T., Hawwari, A. & Al Zubi, K. Whole exome sequencing highlights variants in association with Keratoconus in Jordanian families. BMC Med Genet 21, 177 (2020). https://doi.org/10.1186/s12881-020-01112-z