Molecular Vision 2010; 16:2395-2401 <>
Received 8 September 2010 | Accepted 9 November 2010 | Published 16 November 2010

VSX1 gene analysis in keratoconus

Mukesh Tanwar,1 Manoj Kumar,1 Bhagabat Nayak,2 Dhananjay Pathak,1 Namrata Sharma,2 Jeewan S. Titiyal,2 Rima Dada1

The first two authors contributed equally to the manuscript

1Laboratory For Molecular Reproduction and Genetics, Department of Anatomy, All India Institute of Medical Sciences, Ansari Nagar, New Delhi, India; 2Dr. R.P. Centre for Ophthalmic Sciences, All India Institute of Medical Sciences, Ansari Nagar, New Delhi, India

Correspondence to: Dr. Rima Dada, Associate Professor, Lab. For Molecular Reproduction and Genetics, Department of Anatomy, All India Institute of Medical Sciences, Ansari Nagar, New Delhi, India-110029; Phone: +91-11-26546716; FAX: +91-11-26588663; email:


Purpose: To screen the visual system homeobox 1 (VSX1) gene in keratoconus patients.

Methods: The enntire coding region of VSX1, including intron-exon boundaries were amplified in keratoconus cases (n=50) and controls (n=50). All sequences were analyzed against the ensemble sequence (ENSG00000100987) for VXS1.

Results: Sequencing analysis showed four alterations (p.A182A, p.R217H, p.P237P, and g.25059612C>T) in VSX1 of which g.25059612C>T (in intron 2) was found to be novel. Of these four, p.A182A and p.P237P were present in both cases as well as controls while p.R217H and g.25059612C>T were limited to cases only. All these changes were non-pathogenic.

Conclusions: In our study no pathogenic VSX1 mutation was identified. The role of VSX1 in the pathogenesis of keratoconus is still controversial. VSX1 mutations are responsible for a very small fraction of all observed keratoconus cases. The absence of pathogenic mutations in VSX1 in our patients indicates that other genetic loci like 13q32 as suggested by a recent study may be involved in the pathogenesis of this disorder.


Keratoconus (KTCN; OMIM 148300) is a bilateral, non-inflammatory, and progredient corneal ectasia [1]. There is no specific treatment for keratoconus except corneal transplantation. It has an estimated incidence between 1/500 to 1/2,000 persons throughout the world. The disease usually arises in the teens and stabilizes in the third and fourth decade of life [2]. It occurs with no ethnic or gender preponderance and causes significant visual impairment. Most cases of keratoconus are sporadic but some (5%–10%) have a positive family history [2,3]. In such cases both autosomal dominant and recessive patterns of inheritance have been documented [4-6]. The exact pathogenesis of keratoconus is still unknown. Genome-wide linkage analyses has identified several chromosomal loci and genes that may be associated with keratoconus [6-9]; however, some were eventually excluded [10,11] while for others a conclusive association with the disease is yet to be established. Mutations in the visual system homeobox 1 (VSX1) gene in keratoconus have been reported in different studies [12-15].

VSX1 is a member of the Vsx1 group of vertebrate paired-like homeodomain transcription factors. It has been localized to human chromosome 20p11-q11. Initially VSX1 was chosen for screening mutations in posterior polymorphous corneal dystrophy (PPCD) and keratoconus [8]. VSX1 is considered important in ocular development and is particularly involved in the developing cornea. Expression in human was demonstrated in embryonic craniofacial, adult corneal, and adult retinal cDNA libraries [16]. VSX1 mRNA has been found in the outer tier of the inner nuclear layer of the human retina and the cornea [17]. Mutations in this gene are also associated with posterior polymorphous dystrophy. VSX1 is highly conserved across many species [18-20]. Several mutations, such as p.D144E, p.G160D, p.P247R, p.L159M, p.R166W, and p.H244R have been reported by various groups [14,17,21] but a definite pathogenic role of these mutations in keratoconus is not yet established. In this study we present the results of VSX1 gene analysis in 50 keratoconus patients and controls from north India.


Patient selection and DNA isolation

The research followed the tenets of the Declaration of Helsinki in the treatment of the subject reported herein. The study was approved by institutional review board (IRB # IRB00006862) of the All India Institute of Medical Sciences (AIIMS) and all participants gave their written informed consent. A total of fifty keratoconus patients (Table 1) presented (during April 2009 to April 2010) at the Dr. R. P. Centre for Ophthalmic Sciences (AIIMS, New Delhi, India) were enrolled in this study. Clinical evaluation involved Ultrasonic Pachymetry, videokeratography (VKG), Orbscan, visual testing, fundoscopy, slitlamp-biomicroscopy, and retinoscopy. Of these patients, 29 were males and 21 were females. The mean age of presentation was 18.2 years. Diagnosis of keratoconus involved the presence of characteristic topographic features, such as inferior or central corneal steepening, or an asymmetric bowtie pattern with skewing of the radial axes, and the presence of one or more of the following characteristic clinical features in one or both eyes: conical corneal deformation, munsen sign, corneal stromal thinning, a Fleischer ring or Vogt striae. All cases were sporadic without any family history.

All keratoconus cases secondary to causes like trauma, surgery, Ehlers Danlos syndrome, Osteogenesis Imperfecta, and pellucid marginal degeneration were excluded from the study.

After informed consent, detailed personal, medical and occupational history was collected and a family tree up to three generations was drawn. Fifty ethnically matched normal individuals without any ocular disorder were enrolled as controls. Health information was obtained from controls through the questionnaire; all underwent ophthalmological examination and blood sample (5 ml) was collected in EDTA (EDTA) vacutainers (Greiner Bio-One GmbH, Frickenhausen, Germany) from patients and controls for DNA extraction. DNA was extracted from whole blood samples of all patients and controls using the phenol-chloroform method.

PCR and DNA sequencing

All the coding regions of VSX1 including exon-intron junctions were amplified using a set of five oligonucleotide primers (Table 2). Each reaction was performed in a 25 µl mixture containing 0.2 µM each primer, 0.5 U Taq DNA polymerase (Biogene, New Delhi, India), 2.5 µl 10× PCR buffer (Biogene) with 2.5 mM MgCl2, and approximately 100 ng genomic DNA. Thermal cycling was performed in a thermal cycler (My Cycler; Biorad, Gurgaon, India) under the following conditions: initial denaturation for 3 min at 95 °C; 35 cycles of 94 °C for 30 s, 55 0-60 °C for 45 s, 72 °C for 60 s; and a final extension for 10 min at 72 °C.

All PCR products were analyzed on 1.8% agarose gel stained with ethidium-bromide (EtBr; 10 mg/ml). Agarose gels were analyzed using gel documentation system (Applied Biosystems, Carlsbad, CA). Amplified PCR products were purified using gel/PCR DNA fragments extraction kit (DF100; Geneaid Biotech Ltd., Sijhih City Taiwan). The purified PCR products were sent for sequencing at MCLAB (Molecular Cloning Laboratories, South San Francisco, CA)

Nucleotide sequences were compared with the VSX1ensembl reference sequence.

Insilico analysis of missense mutations

Two homology based programs PolyPhen (Polymorphism Phenotyping) and SIFT (Sorting Intolerant From Tolerant) analysis tool were used to predict the functional impact of missense changes identified in this study.

PolyPhen structurally analyzes an amino acid polymorphism and predicts whether that amino acid change is likely to be deleterious to protein function [22,23]. The prediction is based on the position-specific independent counts (PSIC) score derived from multiple sequence alignments of observations. PolyPhen scores of >2.0 indicate the polymorphism is probably damaging to protein function. Scores of 1.5–2.0 are possibly damaging, and scores of <1.5 are likely benign.

SIFT is a sequence homology-based tool that sorts intolerant from tolerant amino acid substitutions and predicts whether an amino acid substitution in a protein will have a phenotypic effect [24,25]. SIFT is based on the premise that protein evolution is correlated with protein function. Positions important for function should be conserved in an alignment of the protein family, whereas unimportant positions should appear diverse in an alignment. Positions with normalized probabilities less than 0.05 are predicted to be deleterious and, those greater than or equal to 0.05 are predicted to be tolerated.

We have also used improved Splice Site predictor tool [26] for prediction the effect of an intronic nucleotide change on splicing.


DNA sequencing analysis of 50 patients and 50 controls revealed a total of four nucleotide changes (Table 3) of which one was novel and 3 have been previously reported. Details of these changes are given below.

Alanine182Alanine (p.A182A)

In this mutation a single nucleotide adenine (A) was replaced by guanine (G) at g.25059546 (rs12480307); c.546; codon 182 resulted in codon change GCA>GCG resulting in synonymous change p.ala182ala (p.A182A; Figure 1). This change was present as homozygous change in 25 cases and 29 controls.

Arginine217Histidine (p.R217H)

In this mutation a single nucleotide thymine (T) was replaced by adenine at position g.25059442 (rs6138482); cDNA position c.650; codon 217. This change resulted in a codon change from CGC>CAC resulting in non-synonymous change p.arg217his (p.R217H; Figure 2) in protein. This change was present in only one case and was homozygous but was absent in controls.

Proline237Proline (p.P237P)

In this mutation a single nucleotide T was replaced by A at position g.25059381 (rs56157240); c.711; codon 237 resulted in a codon change CCT>CCA which predicts a synonymous change p.pro237pro (p.P237P; Figure 3). This change was present in 18/50 cases (7 were homozygous and 11 were heterozygous) being also present in 14/50 controls (9 were homozygous and 5 were heterozygous).

Cytosine to Thymine in intron 2

A novel single nucleotide change C>T at g.25059612 (Figure 4) was present in three cases but absent in controls. Alteration is located in 2nd intron (IVS3–24C>T). This change was registered at GenBank with accession number GU471016.

PolyPhen and SIFT analysis of p.R217H showed that it is non-pathogenic (SIFT score is >0.05 and PSIC score is <1.5).

Improved splice site predictor tool analysis of g.25059612C>T showed that this location (g.25059612) is not present at splice site and may not create splicing error in VSX1 mRNA. Conclusively, no pathogenic change was identified in our patients.


In this study we analyzed VSX1 in 50 unrelated keratoconus patients and controls from north India. In our patient’s, males were affected more than females. Mutations in VSX1 gene have been identified in association with keratoconus [12,14,15,17,21]. Human VSX1 is a member of the CVC domain containing paired-like class of homeoproteins. VSX1 expression in humans was detected in embryonic craniofacial, adult retinal, and adult corneal tissues [17,27]. The role of VSX1 in keratoconus is still ambiguous. Previous studies have shown that the pathogenesis of keratoconus is very complex and several genes and gene environmental interactions play a critical role in disease prognosis. In fact, VSX1 may have a pleiotropic action among the tissues of the cornea leading to keratoconus in some cases as observed for the transforming growth factor (TGFBI) gene; which causes four distinct autosomal dominant corneal diseases [28].

In this study, four sequence variations were detected of which 3 have been previously reported [15,20] and one was novel. SIFT and PolyPhen analysis of p.217H showed that it to be a non-pathogenic change. Similar findings have already been reported in European population [15].

The VSXI variants reported in various studies include p.L17P, p.D144E, p.N151S, p.L159M, p.G160V, p.G160D, p.R166W, p.Q175H, p.H244R, and p.P247R [12-15,17,21]. Some of these were initially reported to be pathogenic but their pathogenicity could not be confirmed as some of these variants were observed in unaffected individuals also. He’on and associates [14] identified a compound heterozygous change with p.P247R and p.G160D and reported p.G160D to be pathogenic and p.P247R to be nonpathogenic. Another study reported p.P247R as a pathogenic change because it was found to be co-segregating with keratoconus. Similarly, p.D144E mutation was initially reported as pathogenic [20,21,] but subsequent studies identified its presence in unaffected individuals and suggested this to be a non-pathogenic polymorphism [29,30]. The variants p.R166W, p.H244R, and p.L159M have been identified in keratoconus patients but these changes did not segregate with the disease phenotype in their family members and hence were not considered sufficiently significant to support a pathogenic role in keratoconus [31]. Similarly, variants p.G160V and p.N151S have been identified in patients from the Korean population [12] but these changes have not been reported in other populations.

In a recent study (2009) from India, VSX1 was screened in 66 keratoconus cases and a potentially pathogenic change (p.Q175H) was identified in one case only [13]. In this study,  the second from India, no pathogenic change  was indentified in VSX1. Similar results have been published recently [29,32,33]. So lack of possibly pathogenic changes in VSX1 gene in keratoconus patients suggests that mutations of VSX1 could only be responsible for a very small fraction of all observed cases and need to be investigated in different populations. This also suggests that other genetic loci like 13q32 as suggested by Gajecka et al. [34] may be involved in the pathogenesis of keratoconus.


The authors would like to thank the families of the patients for their co-operation. The author Mukesh Tanwar thanks University Grants Commission (UGC) Govt. of India and Manoj Kumar thanks Indian Council of Medical Research (ICMR) for providing Senior Research Fellowship (SRF).


  1. Krachmer JH, Feder RS, Belin MW. Keratoconus and related noninflammatory corneal thinning disorders. Surv Ophthalmol. 1984; 28:293-322. [PMID: 6230745]
  2. Rabinowitz YS. Keratoconus. Surv Ophthalmol. 1998; 42:297-319. [PMID: 9493273]
  3. Kennedy RH, Bourne WM, Dyer JA. A 48-year clinical and epidemiologic study of keratoconus. Am J Ophthalmol. 1986; 101:267-73. [PMID: 3513592]
  4. Hughes AE, Dash DP, Jackson AJ, Frazer DG, Silvestri G. Familial keratoconus with cataract: linkage to the long arm of chromosome 15 and exclusion of candidate genes. Invest Ophthalmol Vis Sci. 2003; 44:5063-6. [PMID: 14638698]
  5. Wang Y, Rabinowitz YS, Rotter JI, Yang H. Genetic epidemiological study of keratoconus: evidence for major gene determination. Am J Med Genet. 2000; 93:403-9. [PMID: 10951465]
  6. Tyynismaa H, Sistonen P, Tuupanen S, Tervo T, Dammert A, Latvala T, Alitalo T. A locus for autosomal dominant keratoconus: linkage to 16q22.3-q23.1 in Finnish families. Invest Ophthalmol Vis Sci. 2002; 43:3160-4. [PMID: 12356819]
  7. Brancati F, Valente EM, Sarkozy A, Fehèr J, Castori M, Del Duca P, Mingarelli R, Pizzuti A, Dallapiccola B. A locus for autosomal dominant keratoconus maps to human chromosome 3p14–q13. J Med Genet. 2004; 41:188-92. [PMID: 14985379]
  8. Héon E, Greenberg A, Kopp KK, Rootman D, Vincent AL, Billingsley G, Priston M, Dorval KM, Chow RL, McInnes RR, Heathcote G, Westall C, Sutphin JE, Semina E, Bremner R, Stone EM. VSX1: a gene for posterior polymorphous dystrophy and keratoconus. Hum Mol Genet. 2002; 11:1029-36. [PMID: 11978762]
  9. Hughes AE, Dash DP, Jackson AJ, Frazer DG, Silvestri G. Familial keratoconus with cataract: linkage to the long arm of chromosome 15 and exclusion of candidate genes. Invest Ophthalmol Vis Sci. 2003; 44:5063-6. [PMID: 14638698]
  10. Rabinowitz YS, Maumenee IH, Lundergan MK, Puffenberger E, Zhu D, Antonarakis S, Francomano CA. Molecular genetic analysis in autosomal dominant keratoconus. Cornea. 1992; 11:302-8. [PMID: 1358551]
  11. Fullerton J, Paprocki P, Foote S, Mackey DA, Williamson R, Forrest S. Identity-by-descent approach to gene localization in eight individuals affected by keratoconus from north-west Tasmania, Australia. Hum Genet. 2002; 110:462-70. [PMID: 12073017]
  12. Mok JW, Baek SJ, Joo CK. VSX1 gene variants are associated with keratoconus in unrelated Korean patients. J Hum Genet. 2008; 53:842-9. [PMID: 18626569]
  13. Paliwal P, Singh A, Tandon R, Titiyal JS, Sharma A. A novel VSX1 mutation identified in an individual with keratoconus in India. Mol Vis. 2009; 15:2475-9. [PMID: 19956409]
  14. Héon E, Greenberg A, Kopp KK, Rootman D, Vincent AL, Billingsley G, Priston M, Dorval KM, Chow RL, McInnes RR, Heathcote G, Westall C, Sutphin JE, Semina E, Bremner R, Stone EM. VSX1: a gene for posterior polymorphous dystrophy and keratoconus. Hum Mol Genet. 2002; 11:1029-36. [PMID: 11978762]
  15. Dash DP, George S, O’Prey D, Burns D, Nabili S, Donnelly U, Hughes AE, Silvestri G, Jackson J, Frazer D, Heon E, Willoughby CE. Mutational screening of VSX1 in keratoconus patients from European population. Eye. 2010; 24:1085-92. [PMID: 19763142]
  16. Semina EV, Mintz-Hittner HA, Murray JC. Isolation and characterization of a novel human paired-like homeodomain-containing transcription factor gene, VSX1, expressed in ocular tissues. Genomics. 2000; 63:289-93. [PMID: 10673340]
  17. Hayashi T, Huang J, Deeb SS. RINX (VSX1), a novel homeobox gene expressed in the inner nuclear layer of the adult retina. Genomics. 2000; 67:128-39. [PMID: 10903837]
  18. Burmeister M, Novak J, Liang MY, Basu S, Ploder L, Hawes NL, Vidgen D, Hoover F, Goldman D, Kalnins VI, Roderick TH, Taylor BA, Hankin MH, McInnes RR. Ocular retardation mouse caused by Chx10 homeobox null allele: impaired retinal progenitor proliferation and bipolar cell differentiation. Nat Genet. 1996; 12:376-84. [PMID: 8630490]
  19. Chow RL, Snow B, Novak J, Looser J, Freund C, Vidgen D, Ploder L, McInnes RR. Vsx1, a rapidly evolving paired-like homeobox gene expressed in cone bipolar cells. Mech Dev. 2001; 109:315-22. [PMID: 11731243]
  20. Bisceglia L, Ciaschetti M, De Bonis P, Campo PA, Pizzicoli C, Scala C, Grifa M, Ciavarella P, Delle Noci N, Vaira F, Macaluso C, Zelante L. VSX1 mutational analysis in a series of Italian patients affected by keratoconus: detection of a novel mutation. Invest Ophthalmol Vis Sci. 2005; 46:39-45. [PMID: 15623752]
  21. Eran P, Almogit A, David Z, Wolf HR, Hana G, Yaniv B, Elon P, Isaac A. The D144E substitution in the VSX1 gene: a nonpathogenic variant or a disease causing mutation. Ophthalmic Genet. 2008; 29:53-9. [PMID: 18484309]
  22. Ramensky V, Bork P, Sunyaev S. Human non-synonymous SNPs: server and survey. Nucleic Acids Res. 2002; 30:3894-900. [PMID: 12202775]
  23. Sunyaev S, Ramensky V, Koch I, Lathe W, III. Kondrashiv, Bork P. Prediction of deleterious human alleles. Hum Mol Genet. 2001; 10:591-7. [PMID: 11230178]
  24. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009; 4:1073-81. [PMID: 19561590]
  25. Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003; 31:3812-4. [PMID: 12824425]
  26. Reese MG, Eeckman FH, Kulp D, Haussler D. Improved Splice Site Detection in Genie. J Comput Biol. 1997; 4:311-23. [PMID: 9278062]
  27. Rozen S, Skaletsky HJ. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Totowa, NJ: Humana Press; 2000. p. 365–386.
  28. Korvatska E, Munier FL, Djemal A, Wang MX, Frueh B, Chiou AG, Uffer S, Ballestrazzi E, Braunstein RE, Forster RK, Culbertson WW, Boman H, Zografos L, Schorderet DF. Mutation hot-spot in 5q31-linked corneal dystrophies. Am J Hum Genet. 1998; 62:320-4. [PMID: 9463327]
  29. Aldave AJ, Yellore VS, Salem AK, Yoo GL, Rayner SA, Yang H, Tang GY, Piconell Y, Rabinowitz YS. No VSX1 gene mutations associated with keratoconus. Invest Ophthalmol Vis Sci. 2006; 47:2820-2. [PMID: 16799019]
  30. Liskova P, Ebenezer ND, Hysi PG, Gwilliam R, El-Ashry MF, Moodaley LC, Hau S, Twa M, Tuft SJ, Bhatacharya SS. Molecular analysis of the VSX1 gene in familial keratoconus. Mol Vis. 2007; 13:1887-91. [PMID: 17960127]
  31. 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:189-92. [PMID: 18216574]
  32. Stabuc-Silih M, Strazisar M, Hawlina M, Glavac D. Absence of pathogenic mutations in VSX1 and SOD genes in patients with keratoconus. Cornea. 2010; 29:172-6. [PMID: 20023586]
  33. Stabuc-Silih M, Strazisar M, Ravnik-Glavak M, Hawlina M, Glavac D. Genetics and clinical characterstics of keratoconus. Acta Dermatovenerol Alp Panonica Adriat. 2010; 19:3-10. [PMID: 20664914]
  34. Gajecka M, Radhakrishna U, Winters D, Nath SK, Rydzanicz M, Ratnamala U, Ewing K, Molinari A, Pitarque JA, Lee K, Leal SM, Bejjani BA. Localization of a gene for keratoconus to a 5.6-Mb interval on 13q32. Invest Ophthalmol Vis Sci. 2009; 50:1531-9. [PMID: 19011015]