Molecular Vision 2011; 17:3128-3136 <http://www.molvis.org/molvis/v17/a337>
Received 6 September 2011 | Accepted 27 November 2011 | Published 30 November 2011

Mutation analysis of VSX1 and SOD1 in Iranian patients with keratoconus

Samira Saee-Rad,1 Hassan Hashemi,2,3 Mohammad Miraftab,3 Mohammad Reza Noori-Daloii,1 Morteza Hashemzadeh Chaleshtori,4 Reza Raoofian,1 Fatemeh Jafari,3 Wayne Greene,5 Ghasem Fakhraie,2 Farhad Rezvan,3 Mansour Heidari1

1Department of Medical Genetics, Tehran University of Medical Sciences, Tehran, Iran; 2Farabi Eye Hospital, Eye Research Center, Tehran University of Medical Sciences, Tehran, Iran; 3Noor Ophthalmology Research Center, Noor Eye Hospital, Tehran, Iran; 4Cellular and Molecular Research Center, Shahrekord University of Medical Sciences, Shahrekord, Iran; 5Biomedical Science, School of Veterinary and Biomedical Sciences, Murdoch University, South Street, Murdoch, Western Australia

Correspondence to: Dr Mansour Heidari, Department of Medical Genetics, Tehran University of Medical Sciences, Pour Sina Ave, Tehran, Iran; Phone: 98-21-8895- 3005; FAX: 98-21-8895- 3005; email: mheidari@sina.tums.ac.ir

Abstract

Purpose: To evaluate mutations in the visual system homeobox gene 1 (VSX1) and superoxide dismutase 1 (SOD1) genes with keratoconus (KTCN), direct sequencing was performed in an Iranian population.

Methods: One hundred and twelve autosomal dominant KTCN patients and fifty-two unaffected individuals from twenty-six Iranian families, as well as one hundred healthy people as controls were enrolled. Genomic DNA was extracted from whole blood sample. Then to study the possible linkage between KTCN and six known loci linkage analysis was performed using 12 short tandem repeat (STR) markers. Also, the entire coding region and intron-exon boundaries of VSX1 and SOD1 were amplified by the PCR technique in each proband. Subsequently, PCR products were subjected to direct sequencing. Co-segregation analysis of the identified mutation was conducted in the family members. An Amplification Refractory Mutation System PCR (ARMS-PCR) was additionally employed for detection of the identified mutation in healthy controls.

Results: Linkage analysis of aforementioned loci did not detect evidence for linkage to KTCN. Direct PCR sequencing revealed two single nucleotide polymorphisms (SNPs; g.1502T>G and g.9683C>T), as well as two missense mutations that have been previously reported (R166W and H244R) in VSX1. We also found three undescribed SNPs (g.4886G>A, g.4990C>G, and g.9061T>A) in SOD1. The R166W and H244R mutations were co-segregated in affected family members but not in those that were unaffected. Moreover, the ARMS-PCR strategy did not detect the identified mutations in controls.

Conclusions: Our data suggest a significant association between KTCN patients and VSX1 genetic alterations (p.R166W and p.H244R). Although our findings support VSX1 as a plausible candidate gene responsible for keratoconus, other chromosomal loci and genes could be involved in KTCN development. Taken together, our results suggest that p.R166W and p.H244R could have possible pathogenic influences on KTCN.

Introduction

Keratoconus (KTCN; OMIM 148300) is a genetically and clinically heterogeneous disease affecting the cornea that causes distortion and reduced vision. KTCN is the most common indication for corneal transplantation in developed countries [1,2]. The estimated incidence is between 1 in 500 and 1 in 2,000 individuals in the general population and the prevalence is estimated to be 54.5 per 100,000 [1]. It occurs in both genders and all ethnicities [3]. Studies suggest that the prevalence and incidence rates are higher in Asians compared to Caucasians [4,5]. A study undertaken in the Midlands area of the UK (UK), a prevalence of 4:1, and an incidence of 4.4:1 was reported in Asians compared to Caucasians while another UK study conducted in Yorkshire, found that the incidence was 7.5 times higher in Asians compared to Caucasians [6].

Several lines of evidence support the importance of genetic components in the pathogenesis of KTCN [6-8]. It has been shown that the prevalence of KTCN in first-degree relatives is significantly higher than the general population [4,5]. Approximately 6% to 23.5% of cases had familial transmission [9] which were inherited in either an X-linked or autosomal recessive or dominant trait. About 90% of pedigrees with familial KTCN exhibit an autosomal dominant inheritance with reduced penetrance, the age of onset in teenage years and variable clinical expression [1,4,6]. In addition, it has been well documented that KTCN is associated with syndromic conditions such as connective tissue disorders (osteogenesis imperfecta, Gapo syndrome, and some subtypes of Ehlers-Danlos syndrome) [1], pigmentary retinopathy, Marfan’s syndrome, Noonan’s syndrome, Apert’s syndrome, Leber congenital amaurosis, and Down syndrome [2,4,10].

KTCN appears to be a genetically heterogeneous disorder as several chromosomal regions and genes are suggested to be involved in the molecular etiology of KTCN. To date, many different chromosomal loci including 20p11-q11 (KTCN1; OMIM 148300) [11], 16q22.3-q23.1 (KTCN2; OMIM 608932) [12], 3p14-q13 (KTCN3; OMIM 608586) [13], 2p24 (KTCN4; OMIM 609271) [14], 15q22.32–24 [15,16], and 5q14.3-q21 [17], as well as genes such as VSX1 (visual system homeobox) and SOD1 (superoxide dismutase) have been implicated in KTCN pathogenesis [3,4,18-23]. In addition to genetic factors, environmental factors as well as genetic-environmental interactions could play critical roles in KTCN [9].

In this study, we conducted mutation detection of the entire VSX1 and SOD1 codon sequences in 112 patients from twenty-six Iranian families. Our results support the possible pathogenic function of VSX1 genetic variants in probands of two families with KTCN. To confirm the association of the identified VSX1 variants, we performed co-segregation analysis in affected family members and in controls. To our knowledge, this is the first report of an association between KTCN and VSX1 genetic alterations in Iranian patients.

Methods

Ophthalmological examination

This study was approved by the local Institutional Review Board (IRB) and informed consent was obtained from all affected individuals. Keratoconus was diagnosed based on the following criteria: (1) distortion of the corneal surface; (2) progressive visual acuity reduction; (3) progressive stromal thinning within central cornea using a comprehensive ophthalmic examination including, visual acuity measurement, corneal imaging including topography, and Pentacam imaging (Noor Eye Hospital, Tehran, Iran).

Molecular genetic studies

Five milliliters of peripheral blood was collected in test tubes containing 0.5 M EDTA. DNA was extracted using a QIAamp DNA mini kit (Qiagen, Hilden, Germany). Genotyping of all family members using twelve short tandem repeat (STR) markers (Table 1) from different chromosomal regions was conducted as previously reported [12-14,19,24]. Briefly, PCR amplification was typically performed in 25 μl PCR reactions, 1 U Taq DNA polymerase, 10 pmole/µl of each primer, (information on the sequences of markers, PCR primers and data on related polymorphic fragments was obtained from the UCSC genome browser,), 200 μM of dNTPs, 0.67 μl of 50 mM MgCl2, 60 ng DNA and 2.5 μl of 10× PCR buffer). The PCR conditions included an initial denaturation step for 3 min at 95 °C, 30 s at 95 °C, 45 s at 64 °C with a 1 °C decrease every second cycle down to 55 °C, then 55 °C for 14 cycles, 1 min at 72 °C for extension, and finally 10 min at 72 °C. PCR-amplified products in genotyping were separated on 12% polyacrylamide gels and the bands detected using silver staining as previously described [25]. Briefly, the gel was fixed in a solution consists of 10% acetic acid and stained with a solution containing silver nitrate and formaldehyde for 30 min. Subsequently, the gel was rinsed and developed in an ice-cold alkaline sodium carbonate solution containing formaldehyde and sodium thiosulphate.

Mutation screening was performed for the complete coding regions of VSX1 and SOD1 using specific primers (Table 2) under the same PCR conditions as above. PCR products were directly sequenced (Gene, Fanavaran, Iran). Sequence data searches were performed in non-redundant nucleic and protein databases BLAST. Also, in this study, we used different softwares such as Chromas, HaploPainter V.024 Beta, easyLinkage PLUS v5.05, and Cyrillic 2.1.

Amplification refractory mutation system PCR (ARMS-PCR)

The ARMS-PCR method was applied for the screening of identified mutations in VSX1 in unaffected and control individuals. Primer sequences are shown in Table 2. Mismatches were included to maximize discrimination of the wild-type and mutant alleles. Amplifications were performed as outlined above. All samples under investigation were analyzed simultaneously alongside positive and negative controls for the R166W and H244R mutations. To determine the genotype results, amplification products were resolved on 2% agarose gels stained with ethidium bromide.

Results

Probands from 26 Iranian families were identified with KTCN. All affected and controls were born after a normal term pregnancy. Retinoscopy, corneal topography, and Pentacam examinations of probands’ first-degree relatives discriminated KTCN from normal individuals. Affected cases were clinically examined and showed no signs and symptoms of any syndromic indication. KTCN was excluded in normal controls by retinoscopy, corneal topography, and Pentacam examination.

Linkage analysis results

To investigate the association of six known (20p11-q11, 16q22.3-q23, 3p14-q13, 2p24, 15q22.32–24, and 5q14.3-q21) genetic loci with KTCN, genomic DNA was isolated from 112 affected 52 unaffected family members. First, PCR primers were used to amplify polymorphic markers (Table 1) on chromosome 20p11-q11, 16q22.3-q23, 3p14-q13, 2p24, 15q22.32–24, and 5q14.3-q21 known to be linked to KTCN. Then, PCR products were separated by PAGE. Haplotype analysis could not define linkage between known loci with KTCN in the studied families (as shown in Appendix 1).

Mutation analysis of VSX1

Direct PCR sequencing using forward and reverse primers (Table2) was conducted to evaluate genetic alterations in the coding sequence and exon/intron boundaries of VSX1. We found two single nucleotide polymorphisms (SNPs) (g.1502T>G and g.9683C>T) and two non-synonymous mutations (H244R and R166W) in VSX1 (Figure 1 and Figure 2).

In family 1, the H244R mutation was identified in a 21-year-old female. The proband had three affected relatives (Figure 1A). To test the pathologic function of this genetic mutation, co-segregation of H244R (Figure 1B) in other affected family members were performed using PCR sequencing. We also screened the H244 R VSX1 mutation in 100 controls by ARMS-PCR (Figure 1C). Direct sequencing and ARMS-PCR results confirmed that only affected patients carried p.H244R in the heterozygous state, while Figure 1B (right) indicates the wild-type (WT) variant of VSX1 gene at codon position 244.

In family 2, the R166W mutation was detected in an 18-year-old male (Figure 2A). Figure 2B indicates R166W in II-1 (proband) and his father (I-1). Although the proband and his father presented a variable expressivity of KTCN, the R166W co-segregated among KTCN, but not in unaffected his mother.

Mutation analysis of SOD1

We additionally evaluated a possible association of KTCN with SOD1 genetic alterations. The full-length SOD1 coding sequence was screened by direct PCR sequencing. In spite of the fact that three novel SNPs (g.4886G>A, g.4990C>G, and g.9061T>A) were identified in non-coding sequences, however, sequencing the coding region did not reveal a sequence variant segregating with disease in any of the families described. The SNPs did not seem to influence the activity of SOD1 protein.

Discussion

Various genome-wide linkage analyses and mutation detection studies have reported that six loci and two genes (SOD1 and VSX1) are thought to be associated with KTCN in different ethnic groups [6]. In this study, we performed linkage analysis for six known chromosomal loci as well as a mutation detection screen in VSX1 and SOD1 in 112 affected, 52 unaffected family members and 100 normal individuals as controls in a group of Iranian patients with keratoconus. All studied families displayed an autosomal dominant pattern of KTCN with variable expressivity. According to literature KTCN is a complex disorder with multiplfactorial etiology. It has been well documented that genetic and environmental factors are associated with this disease. Different studies have shown that about 90 percent of pedigrees with familial KTCN display an autosomal dominant inheritance with incomplete penetrance [11,16]. Our pedigree analysis showed a reduced penetrance (about 58%) among Iranian families with KTCN.

Genotyping results using linkage analysis failed to detect any significant correlations between known loci and KTCN. Direct PCR sequencing did not show important genetic changes in SOD1, although three novel SNPs were found in the probands of three families with KTCN.

For VSX1, we observed two previously undescribed SNPs (g.1502T>G and g.9683C>T) in addition to two missense mutations that had been previously reported (R166W and H244R) [11].

The VSX1 homeobox gene encodes a homeodomain transcription factor that may regulate expression of the cone opsin genes early in normal development [26].Despite the fact that VSX1 remains the only major genetic element to be identified in KTCN pathogenesis, there has been debate in the literature in regard to its pathobiological function in KTCN. So far, several VSX1 protein coding changes (p.L17P, p.D144E, p.L159M, p.G160D, p.R166W, p.H244R, and p.P247R) have been reported as being potentially pathogenic in KTCN [3,27,28]. Valleix et al. [29] reported that H244R VSX1 is associated with selective cone ON bipolar cell dysfunction and macular degeneration in a posterior polymorphous corneal dystrophy (PPCD) family [29]. However, different studies in several ethnic groups [3,30,31] did not detect any mutations in VSX1 [21,32,33].

To evaluate the possible function of p.H244R genetic alteration we undertook co-segregation analysis and mutation detection in the subject families and 100 unaffected individuals, respectively. The mutation was not observed in 100 unaffected people who were examined by corneal topography followed by Pentacam evaluation. Therefore, for first time, our results support the notion that p.H244R could represent a pathogenic change in KTCN, although it has been previously identified [11]. Héon et al. [11] showed that the VSX1 mutation p.H244R co-segregated with the disease in two family members. However, they also observed this genetic change in two controls (n=277). For that reason, they suggested that this change may or may not be disease-causing [11]. These results could be due to several factors, such as the mode of inheritance, gene-gene interaction, gene-environment interaction, and genetic background. About 90% of pedigrees with familial KTCN are predominately transmitted in an autosomal-dominant fashion with an incomplete penetrance and variable expressivity [9]. A comparison of the findings obtained in our study and from other studies raises a critical question: what could be the biologic significance of the observation of VSX1 H244R alteration in KTCN? With respect to this event, there are several possibilities that may account for the presence of the VSX1 H244R change observed in our family including:

random appearance in the affected members without any biologic significance,

the H244R mutation may not directly cause the disease but its pathogenic functions could be due to gene-gene interactions, gene-environment interactions, or genetic background, or

the presence of this change in the unaffected individual suggests that this could be a causative mutation with incomplete penetrance.

Our results suggest that the H244R VSX1 change may be a pathogenic variant with incomplete penetrance.

The R166W VSX1 mutation was also initially found in an isolated case of keratoconus with visual impairment for whom a corneal graft was required in adulthood [11]. We observed this mutation in the proband of a family with two affected patients with KTCN. Our results showed that R166W was co-segregated in two affected family members, but not in unaffected individuals and controls. The R166W alters the highly conserved third amino acid of the DNA binding homeodomain (HD). Dorval et al. [34] suggested that this mutation causes keratoconus in humans by impairing VSX1 DNA binding [34].

Three SNPs were identified in non-coding sequences of the SOD1 gene. This gene is located on chromosome 21 and functions to destroy free superoxide radicals in the body [22]. The association of SOD1 mutation with KTCN was originally reported by Udar et al. [22]. SOD1 mutations were screened in 15 unrelated individuals, each with a family history of KTCN. Results from this study determined an IVS2+50del7 change within intron 2 in two families. Also, they observed that this 7-base deletion segregated with the KTCN subjects in a studied pedigree [22]. However, an independent study failed to define any association between SOD1 mutations and KTCN [34]. Stabuc-Silih and colleagues [33] studied the association of KTCN with VSX1 and SOD1 gene mutations in 113 Slovenian patients with sporadic and familial KTCN by direct sequencing. They found no causative disease mutation in the SOD1 gene but a significant association was detected between a VSX1 polymorphism (627+23G>A) and KTCN [33]. Our findings are in agreement with published data suggesting that other genetic and non-genetic factors are involved the pathogenesis of KTCN.

We conclude that the R166W and H244R VSX1 variants might play critical roles in the pathogenesis of KTCN. To test the potential pathogenic relevance of these variants we performed co-segregation analysis in all affected and unaffected family members. To rule out the incomplete penetrance of autosomal dominant KTCN in unaffected family members, we conducted precise clinical examinations of unaffected individuals; a crucial step in an association study involving a molecular genetic analysis of an autosomal dominant disorder. An accurate diagnosis and exclusion of keratoconus in patients and controls (n=100), respectively, was made using corneal topography and Pentacam evaluation. Our findings suggest that the R166W and H2446R mutations might be involved in KTCN pathogenesis.

Appendix 1. Twenty six pedigrees of KTCN.

Acknowledgments

We thank all the patients and healthy participants in this project. This work was supported by a grant from the Tehran University of Medical Sciences (TUMS).

References

  1. 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]
  2. Stabuc-Silih M, Ravnik-Glavac M, Glavac D, Hawlina M, Strazisar M. Polymorphisms in COL4A3 and COL4A4 genes associated with keratoconus. Mol Vis. 2009; 15:2848-60. [PMID: 20029656]
  3. Abu-Amero KK, Kalantan H, Al-Muammar AM. Analysis of the VSX1 gene in keratoconus patients from Saudi Arabia. Mol Vis. 2011; 17:667-72. [PMID: 21403853]
  4. 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]
  5. Rabinowitz YS. The genetics of keratoconus. Ophthalmol Clin North Am. 2003; 16:607-20. [PMID: 14741001]
  6. Romero-Jiménez M, Santodomingo-Rubido J, Wolffsohn JS. Keratoconus: a review. Cont Lens Anterior Eye. 2010; 33:157-66. [PMID: 20537579]
  7. Newsome DAFJ, Hassell JR, Krachmer JH, Rodrigues MM, Katz SI. Detection of specific collagen types in normal and keratoconus corneas. Invest Ophthalmol Vis Sci. 1981; 20:738-50. [PMID: 7016805]
  8. Wang Y, Rotter JI, Yang H. Genetic epidemiological study of keratoconus: evidence for major gene determination. Am J Med Genet. 2000; 93:403-9. [PMID: 10951465]
  9. Bisceglia L, De Bonis P, Pizzicoli C, Fischetti L, Laborante A, Di Perna M, Giuliani F, Delle Noci N, Buzzonetti L, Zelante L. Linkage analysis in keratoconus: replication of locus 5q21.2 and identification of other suggestive Loci. Invest Ophthalmol Vis Sci. 2009; 50:1081-6. [PMID: 18978346]
  10. Cullen JF. Mongolism (Down's Syndrome) and keratoconus. Br J Ophthalmol. 1963; 47:321-30. [PMID: 14189698]
  11. 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]
  12. 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]
  13. Brancati F, Valente EM, Sarkozy A, Feher 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]
  14. Hutchings H, Ginisty H, Le Gallo M, Levy D, Stoesser F, Rouland JF, Arne JL, Lalaux MH, Calvas P, Roth MP, Hovnanian A, Malecaze F. Identification of a new locus for isolated familial keratoconus at 2p24. J Med Genet. 2005; 42:88-94. [PMID: 15635082]
  15. Dash DP, Silvestri G, Hughes AE. Fine mapping of the keratoconus with cataract locus on chromosome 15q and candidate gene analysis. Mol Vis. 2006; 12:499-505. [PMID: 16735990]
  16. 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]
  17. Tang YG, Rabinowitz YS, Taylor KD, Li X, Hu M, Picornell Y, Yang H. Genomewide linkage scan in a multigeneration Caucasian pedigree identifies a novel locus for keratoconus on chromosome 5q14.3-q21.1. Genet Med. 2005; 7:397-405. [PMID: 16024971]
  18. Bisceglia L, Ciaschetti M, De Bonis P, Campo PAP, 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]
  19. 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 the European population. Eye (Lond). 2010; 24:1085-92. [PMID: 19763142]
  20. Mok J-W, Baek S-J, Joo C-K. VSX1 gene variants are associated with keratoconus in unrelated Korean patients. J Hum Genet. 2008; 53:842-9. [PMID: 18626569]
  21. 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]
  22. Udar N, Atilano SR, Brown DJ, Holguin B, Small K, Nesburn AB, Kenney MC. SOD1: a candidate gene for keratoconus. Invest Ophthalmol Vis Sci. 2006; 47:3345-51. [PMID: 16877401]
  23. Udar N, Atilano SR, Small K, Nesburn AB, Kenney MC. SOD1 haplotypes in familial keratoconus. Cornea. 2009; 28:902-7. [PMID: 19654524]
  24. Héon E, Mathers WD, Alward WL, Weisenthal RW, Sunden SL, Fishbaugh JA, Taylor CM, Krachmer JH, Sheffield VC, Stone EM. Linkage of posterior polymorphous corneal dystrophy to 20q11. Hum Mol Genet. 1995; 4:485-8. [PMID: 7795607]
  25. Saffari M, Dinehkabodi OS, Ghaffari SH, Modarressi MH, Mansouri F, Heidari M. Identification of novel p53 target genes by cDNA AFLP in glioblastoma cells. Cancer Lett. 2009; 273:316-22. [PMID: 18814959]
  26. 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]
  27. Grünauer-Kloevekorn C, Duncker GI. Keratoconus: epidemiology, risk factors and diagnosis. Klin Monatsbl Augenheilkd. 2006; 223:493-502. [PMID: 16804819]
  28. Tanwar M, Kumar M, Nayak B, Pathak D, Sharma N, Titiyal JS, Dada R. VSX1 gene analysis in keratoconus. Mol Vis. 2010; 16:2395-401. [PMID: 21139977]
  29. Valleix S, Nedelec B, Rigaudiere F, Dighiero P, Pouliquen Y, Renard G, Le Gargasson JF, Delpech M. H244R VSX1 is associated with selective cone ON bipolar cell dysfunction and macular degeneration in a PPCD family. Invest Ophthalmol Vis Sci. 2006; 47:48-54. [PMID: 16384943]
  30. Cozma I, Atherley C, James NJ. Influence of ethnic origin on the incidence of keratoconus and associated atopic disease in Asian and white patients. Eye (Lond). 2005; 19:924-5. [PMID: 15389278]
  31. Pearson AR, Soneji B, Sarvananthan N, Sandford-Smith JH. Does ethnic origin influence the incidence or severity of keratoconus? Eye (Lond). 2000; 14:625-8. [PMID: 11040911]
  32. 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]
  33. Stabuc-Silih M, Strazisar M, Hawlina M, Glavac D. Absence of pathogenic mutations in VSX1 and SOD1 genes in patients with keratoconus. Cornea. 2010; 29:172-6. [PMID: 20023586]
  34. Dorval KM, Bobechko BP, Ahmad KF, Bremner R. Transcriptional activity of the paired-like homeodomain proteins CHX10 and VSX1. J Biol Chem. 2005; 280:10100-8. [PMID: 15647262]