Molecular Vision 2006;
12:499-505 <http://www.molvis.org/molvis/v12/a58/> Received 30 September 2005 | Accepted 26 April 2006 | Published 12 May 2006 |
Download Reprint |
Fine mapping of the keratoconus with cataract locus on chromosome 15q and candidate gene analysis
Durga Prasad Dash,1
Giuliana Silvestri,2
Anne E. Hughes1
Departments of 1Medical Genetics and 2Ophthalmology, Queen's University of Belfast, Belfast, United Kingdom
Correspondence to: Durga Prasad Dash, Department of Ophthalmology, Institute of Clinical Science, Block A, Queen's University of Belfast, Royal Victoria Hospital, Grosvenor Road, Belfast, BT12 6BA, United Kingdom; Phone: 0044-7754404969; FAX: 0044-2890632699; email: d.dash@queens-belfast.ac.uk
Abstract
Purpose: To report the fine mapping of the keratoconus with cataract locus on chromosome 15q and the mutational analysis of positional candidate genes.
Methods: Genotyping of two novel microsatellite markers and a single nucleotide polymorphism (SNP) in the critical region of linkage for keratoconus with cataract on 15q was performed. Positional candidate genes (MORF4L1, KIAA1055, ETFA, AWP1, REC14, KIAA1199, RCN2, FAH, IDH3A, MTHFS, ADAMTS7, MAN2C1, PTPN9, KIAA1024, ARNT2, BCL2A1, ISL2, C15ORF22 (P24B), DNAJA4, FLJ14594, CIB2 (KIP2), C15ORF5, and PSMA4) prioritized on the basis of ocular expression and probable function were screened by PCR-based DNA sequencing methods.
Results: We report the refinement of the linkage region for keratoconus with cataract to an interval of approximately 5.5 Mb flanked by the MAN2C1 gene and the D15S211 marker on chromosome 15q. Mutational analysis of positional candidate genes detected many sequence variations and single nucleotide polymorphisms. None of the sequence variants were considered pathogenic as they were also found in unaffected family members and normal control DNA samples.
Conclusions: Fine mapping of the keratoconus with cataract locus on 15q has reduced the linked region to 5.5 Mb, thereby excluding 28 candidate genes. A further 23 candidate genes were excluded by direct sequencing methods, although a pathogenic genomic rearrangement or exonic deletion would not have been detected.
Introduction
Keratoconus (KC; OMIM 148300), the most common corneal dystrophy, is a bilateral, noninflammatory progressive corneal ectasia. Clinically, the cornea becomes progressively thin and conical which leads to myopia, irregular astigmatism, and corneal scarring. The transparency and refractive state of the cornea is a prerequisite for normal vision. The disease usually arises in the teenage years, eventually stabilizing in the third and fourth decades [1]. It occurs with no ethnic or gender preponderance and causes significant visual impairment in young adults. No specific treatment exists except to replace the corneal tissue by surgery (corneal transplantation) when the visual acuity can no longer be corrected by contact lenses. In the Western world, KC is the most common indication for corneal transplantation [2]. The underlying biochemical processes and pathobiology of keratoconus remain poorly understood.
The incidence of KC is 1 in 2,000 in the general population [1,3]. A family history is present only in a minority of cases, however, one of the major etiological factors is genetic [1,4-6]. KC is believed to be inherited autosomally, because of familial occurrence [7], a higher concordance rate of the trait in monozygotic twins then dizygotic twins [8] and its prevalence in first degree relatives is 15-67 times higher than the general population [6]. Most studies describe autosomal dominant inheritance, with incomplete penetrance or variable expressivity [1,9,10].
Mutations in the VSX-1 transcription factor (OMIM 605020) were identified in 4.7% of patients with keratoconus (KTCN1; OMIM 148300) and also in posterior polymorphous corneal dystrophy (PPCD1; OMIM 122000) [11]. However, traditional mapping methods to identify the genetic basis of KC have been limited by the lack of large multigeneration families for study. Genome-wide scans to localize the KC gene or genes had not been reported until our group identified a large three-generation family with KC and anterior polar cataract and reported linkage to chromosome 15q22.32-24.2 [12]. Recently loci for autosomal dominantly inherited keratoconus have been mapped to chromosomes 21 [13], 16q22.3-q23.1 (KTCN2; OMIM 608932) [10], 3p14-q13 (KTCN3; OMIM 608586) [14], 2p24 (KTCN4; OMIM 609271) [15] and 5q14.1-q21.3 (NCBI).
The keratoconus with cataract locus on chromosome 15q was identified in a large Northern Irish family of three generations affected by combined autosomal dominant early onset anterior polar cataract and clinically severe keratoconus [12]. The disease gene in this family was successfully mapped to chromosome 15q22.32-24.2 within a 6.5 Mb region flanked by markers CYP11A and D15S211 [12]. CTSH, CRABP1, IREB2, and RASGRF1 were excluded previously as the causative gene for keratoconus with cataract [12]. The purpose of the present study was to narrow the critical region of linkage of keratoconus with cataract by genotyping novel microsatellite markers and single nucleotide polymorphisms identified during sequencing of candidate genes. Fine mapping of a keratoconus with cataract locus on 15q has reduced the linked region to 5.5 Mb, thereby excluding 28 potential candidate genes. A further 23 positional candidate genes were prioritized for analysis based on expression studies, molecular and biochemical evidence (Table 1) and excluded by PCR-based DNA sequencing.
Methods
Genomic DNA was extracted from the peripheral blood leukocytes of all available family members using Puregene DNA purification kit (Gentra Systems, Inc., Minneapolis, MN).
Single nucleotide polymorphism typing
The single nucleotide polymorphism (SNP) rs1128933 in exon 20 of MAN2C1 (GenBank NM_006715) was typed in key recombinants using the forward primer 5'-CTG GAG ACA CGG TAT AGG CTG G-3' and the reverse primer 5'-AGT GTA CCT GGG AGT GGG AAG G-3'. Cycle sequencing was performed with the BigDye Terminator Cycle Sequencing kit (Applied Biosystems, Foster City, CA) on PCR products using the forward PCR primer and analyzed on the ABI Prism 3100 DNA Sequencing System (Applied Biosystems).
Identification and typing of microsatellite markers
Two novel microsatellite markers on chromosome 15 located at 73.8 Mb between MAN2C1 and ETFA and at 78.4 Mb (within ARNT2) were identified from the human genome database (HGD). PCR was performed on genomic DNA from family members using the following primers: 5'-Fam-CTG TAT GCA GCG ATC AGT GAG C-3' and 5'-GTA GGA GAC AGA GAC ACT CAG TCC-3' for the marker at 73.8 Mb and 5'-Fam-TTA CAT CAC TTG CAT TGC TTC C-3' and 5'-AAT CAC AGC CAA ACA TAG ATG C-3' for the ARNT2 intragenic marker. Diluted PCR products were electrophoresed on an ABI 3100 genetic analyzer. Alleles were sized by Genotyper software (ABI) and scored using in-house software, C-allele (C-ALLELE program, Dr. Shane Mckee, Department of Medical Genetics, Queen's University of Belfast, UK). Linkage analysis was performed using the Fastlink linkage program version 5.1.
PCR-based DNA sequencing
Analysis of positional candidate genes was carried out by PCR-based DNA sequencing using dye terminator chemistry (ABI). Primers were designed using Primer Detective, version 1.01 (Clontech Labs Inc., CA) for all individual exons plus at least 50 bp of flanking intronic sequence, and 5' and 3' untranslated regions. For ADAMTS7, long-range PCR using Elongase (Invitrogen, Paisley, UK) was carried out as a preliminary step, using gene-specific primers that did not amplify related pseudogenes. After dilution of product (1 in 1000) nested PCR was performed before sequencing as already described. All primer pairs and PCR conditions are available on request.
Results
One SNP and two novel microsatellite markers were identified toward the extremities of the previously linked interval on chromosome 15q [12] in a family with autosomal dominant keratoconus and cataract. Utilizing these markers, family members were genotyped to clarify the position of the recombination in key individuals. Typing of rs1128933 (Figure 1) within exon 20 of MAN2C1 showed III23 is recombinant at the MAN2C1 locus and more proximal markers, whereas III21, previously shown to be recombinant at CYP11A, is not recombinant at MAN2C1. The novel microsatellite markers located at 73.8 and 78.4 Mb (ARNT2) were fully informative for tracking haplotypes inherited by the crucial individuals in this family; however, none showed recombination at these loci. The further haplotype analysis described here has reduced the candidate interval of the keratoconus with cataract locus by 1.0 Mb, to the region flanked by rs1128933 (within MAN2C1) and D15S211.
Candidate genes positioned within this 5.5 Mb region (Figure 2) were identified and prioritized for analysis based on expression studies, molecular, and biochemical evidence (Table 1). Numerous sequence variations were found (Table 2) in positional candidate genes (MORF4L1, KIAA1055, ETFA, AWP1, REC14, KIAA1199, RCN2, FAH, IDH3A, MTHFS, ADAMTS7, MAN2C1, PTPN9, KIAA1024, ARNT2, BCL2A1, ISL2, C15ORF22 (P24B), DNAJA4, FLJ14594, CIB2 (KIP2), C15ORF5, and PSMA4). None of the sequence variants were considered pathogenic as they were also found in unaffected family members and unrelated, ethnically matched normal control DNA samples.
Discussion
The refinement of the linkage region for the keratoconus with cataract locus on chromosome 15q [15] to 5.5 Mb between MAN2C1 and D15S211 has excluded 28 positional candidate genes. CTSH, CRABP1, IREB2, and RASGRF1 which remained within the refined disease interval were excluded previously by direct sequencing [15]. The reduced disease gene interval is exceedingly gene rich containing a total of 75 positional candidate genes, of which 27 have now been screened without identifying a causative mutation. There is limited scope for further refinement of the linkage region, which cannot be reduced below a core interval of 4.6 Mb with the present family members. One family has been reported in the literature with corneal pathology and congenital cataract which maps to chr15q [16]. The authors of this study have proposed the term EDICT syndrome for the clinical findings of endothelial dystrophy, iris hypoplasia, congenital cataract, and stromal thinning (EDICT) which were inherited as an autosomal dominant trait in their family. The locus for EDICT syndrome maps to a 26 Mb region of chromosome 15q which overlaps our region, and may therefore be allelic, however the small size of this family restricts its potential for improved mapping. Although the clinical phenotype in this family shows considerable differences from the EDICT syndrome [16], it is possible that both families share a common genetic basis. The EDICT syndrome includes features of significant anterior dysgenesis such as microcornea, iris hypoplasia, endothelial abnormalities and importantly the snydrome been reported to have congenital anterior polar cataract. Our family differs in that those affected have been documented as having no abnormalities or signs of anterior segment dysgenesis at birth, with the corneal and lens changes developing after the age of 5 years. Sequential examination of our family has indicated that the anterior polar cataract and keratoconus are developmental and not congenital. The ocular phenotype in our family is also unusual as the occurrence of keratoconus and developmental anterior polar cataract has not to our knowledge been previously described. Isolated congenital cataracts are usually inherited in an autosomal dominant fashion with mutations predominantly reported in the lens crystallins and connexins [17]. Anterior polar cataracts are developmental lens opacities and have been mapped to 14q24-qter (CTAA1; OMIM 115650) [18] and 17p13 (CTAA2; OMIM 601202) [19]. Corneal guttata (focal thickenings of the Descemet's membrane of the corneal endothelium) have been reported in association with anterior polar cataract (OMIM 121390) [20,21]. Developmentally the human lens is derived from surface ectoderm and the cornea from surface ectoderm and neuroectoderm. The ocular phenotype in this family may reflect a single developmental anomaly. Pathogenic mutations have been reported in a number of transcription factors which result in ocular phenotypes combining cataract, anterior segment dysgenesis and specific iris defects (aniridia, coloboma) [22-24].
Human VSX-1 (OMIM 605020) is a member of the Vsx1 group of vertebrate paired-like homeodomain transcription factors [25,26]. Mutations in the VSX-1 transcription factor were identified in 7% [11] of keratoconus patients [27]. To date, five pathogenic heterozygous missense mutations have been reported in VSX-1 associated with keratoconus: L17P, R166W, L159M, D144E, and P247R [11,27]. The aryl hydrocarbon receptor nuclear translocator (ARNT2) in our candidate interval is also member of a novel transcription factor family consisting of a conserved basic helix-loop-helix (bHLH) structural motif contiguous with a PAS domain. A W278X mutation in the aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1) gene on 17p was found in a Pakistani family with Leber's congenital amaurosis and anterior keratoconus [28]. We have not found any pathogenic mutation in ARNT2 or in ISL2 or MORF4L1 (also called MRG15) which also code for transcription factors.
ADAMTS7 was a strong candidate gene as the ADAMTS (a disintegrin-like and metalloprotease domain with thrombospondin type 1 motif) families are secreted enzymes metalloproteases with a prominent role in extracellular matrix proteolysis [29] The proteolytic activity of this gene may be related to development of keratoconus in which there is evidence implicating matrix metalloproteinases in the pathogenesis [30]. Other ADAMTS7 family members such as ADAMTS18 (16q22.3-q23.1) and ADAMTS9 (3p14-q13) have been mapped within other reported keratoconus loci [10,13]. ADAMTS7 presented sequencing difficulties because of its extremely close homology to several genes located on chromosome 15q and elsewhere, necessitating the use of long-range PCR to ensure specificity. Although no causative mutation was found in ADAMTS7, several novel SNPs were identified.
The most interesting of the remaining positional candidate genes are RNUT1, IMP3 (C15orf12), CSPG4, FBXO22, PSTPIP1, HMG20A, LRRN6A, ACSBG1 (BG1), CHRNA5, CHRNA3, CHRNB4, transcription factor genes ZNF291 and SIN3A1 and will be studied in the future.
The pathogenic mutation for keratoconus with cataract in this family may be in one of the lower priority genes that have not yet been screened, however, the possibility of a mutation in a high priority gene that is not amenable to detection by PCR-based sequencing must also be considered. Exonic deletions, duplications and other rearrangements would fall into this category. We have applied multiplex ligation-dependent probe amplification (MLPA) [31-34] to three genes (data not shown) and have found no evidence of copy number change. The mutation may also be located in a conserved intergenic region of regulatory importance. Chromosome 15q shows a high frequency of neocentromeres, which harbour a high density of clinically important duplicons, indicating a significant evolutionary history of this region [35]. Rearrangements in this region have previously been associated with diseases [36-41].
Quantitative analysis of SNP genotypes at numerous loci spanning the critical region of the keratoconus with cataract locus is included in future research plans. This will enable more precise mapping of the linkage interval and also provide evidence of any anomalous dosage effects that may be of relevance to the disease process. Identification of genes implicated in the pathogenesis of keratoconus is vital for the development of novel treatments and early preventative strategies for this common corneal disorder.
Acknowledgements
This work was supported by grants from The Guide Dogs for the Blind Association, The Northern Ireland HPSS Research and Development Office, and The Royal College of Surgeons of Edinburgh. We are most grateful to the family who participated in this study. We are thankful to Mr. C. E. Willoughby, Department of Ophthalmology and Vision Sciences, Queen's University of Belfast, UK for discussions and suggestions regarding the project.
References
1. Rabinowitz YS. Keratoconus. Surv Ophthalmol 1998; 42:297-319.
2. Dobbins KR, Price FW Jr, Whitson WE. Trends in the indications for penetrating keratoplasty in the midwestern United States. Cornea 2000; 19:813-6.
3. Kennedy RH, Bourne WM, Dyer JA. A 48-year clinical and epidemiologic study of keratoconus. Am J Ophthalmol 1986; 101:267-73.
4. Ihalainen A. Clinical and epidemiological features of keratoconus genetic and external factors in the pathogenesis of the disease. Acta Ophthalmol Suppl 1986; 178:1-64.
5. Bron AJ. Keratoconus. Cornea 1988; 7:163-9.
6. 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.
7. Falls HF, Allen AW. Dominantly inherited keratoconus. J Genet Hum 1969; 17:317-24.
8. Parker J, Ko WW, Pavlopoulos G, Wolfe PJ, Rabinowitz YS, Feldman ST. Videokeratography of keratoconus in monozygotic twins. J Refract Surg 1996; 12:180-3.
9. Edwards M, McGhee CN, Dean S. The genetics of keratoconus. Clin Experiment Ophthalmol 2001; 29:345-51.
10. 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.
11. Heon 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.
12. 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.
13. 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.
14. 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.
15. 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.
16. Jun AS, Broman KW, Do DV, Akpek EK, Stark WJ, Gottsch JD. Endothelial dystrophy, iris hypoplasia, congenital cataract, and stromal thinning (edict) syndrome maps to chromosome 15q22.1-q25.3. Am J Ophthalmol 2002; 134:172-6.
17. Graw J. Congenital hereditary cataracts. Int J Dev Biol 2004; 48:1031-44.
18. Moross T, Vaithilingam SS, Styles S, Gardner HA. Autosomal dominant anterior polar cataracts associated with a familial 2;14 translocation. J Med Genet 1984; 21:52-3.
19. Berry V, Ionides AC, Moore AT, Plant C, Bhattacharya SS, Shiels A. A locus for autosomal dominant anterior polar cataract on chromosome 17p. Hum Mol Genet 1996; 5:415-9.
20. Dohlman CH. Familial congenital cornea guttata in association with anterior polar cataract. Acta Ophthalmol Suppl 1951; 29:445-73.
21. Traboulsi EI, Weinberg RJ. Familial congenital cornea guttata with anterior polar cataracts. Am J Ophthalmol 1989; 108:123-5.
22. Jamieson RV, Munier F, Balmer A, Farrar N, Perveen R, Black GC. Pulverulent cataract with variably associated microcornea and iris coloboma in a MAF mutation family. Br J Ophthalmol 2003; 87:411-2.
23. Hanson I, Churchill A, Love J, Axton R, Moore T, Clarke M, Meire F, van Heyningen V. Missense mutations in the most ancient residues of the PAX6 paired domain underlie a spectrum of human congenital eye malformations. Hum Mol Genet 1999; 8:165-72.
24. Azuma N, Yamaguchi Y, Handa H, Hayakawa M, Kanai A, Yamada M. Missense mutation in the alternative splice region of the PAX6 gene in eye anomalies. Am J Hum Genet 1999; 65:656-63.
25. 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.
26. 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.
27. 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.
28. Damji KF, Sohocki MM, Khan R, Gupta SK, Rahim M, Loyer M, Hussein N, Karim N, Ladak SS, Jamal A, Bulman D, Koenekoop RK. Leber's congenital amaurosis with anterior keratoconus in Pakistani families is caused by the Trp278X mutation in the AIPL1 gene on 17p. Can J Ophthalmol 2001; 36:252-9.
29. Somerville RP, Longpre JM, Jungers KA, Engle JM, Ross M, Evanko S, Wight TN, Leduc R, Apte SS. Characterization of ADAMTS-9 and ADAMTS-20 as a distinct ADAMTS subfamily related to Caenorhabditis elegans GON-1. J Biol Chem 2003; 278:9503-13.
30. Ye HQ, Maeda M, Yu FS, Azar DT. Differential expression of MT1-MMP (MMP-14) and collagenase III (MMP-13) genes in normal and wounded rat corneas. Invest Ophthalmol Vis Sci 2000; 41:2894-9.
31. Schouten JP, McElgunn CJ, Waaijer R, Zwijnenburg D, Diepvens F, Pals G. Relative quantification of 40 nucleic acid sequences by multiplex ligation-dependent probe amplification. Nucleic Acids Res 2002; 30:e57.
32. Stern RF, Roberts RG, Mann K, Yau SC, Berg J, Ogilvie CM. Multiplex ligation-dependent probe amplification using a completely synthetic probe set. Biotechniques 2004; 37:399-405.
33. White SJ, Vink GR, Kriek M, Wuyts W, Schouten J, Bakker B, Breuning MH, den Dunnen JT. Two-color multiplex ligation-dependent probe amplification: detecting genomic rearrangements in hereditary multiple exostoses. Hum Mutat 2004; 24:86-92.
34. Sellner LN, Taylor GR. MLPA and MAPH: new techniques for detection of gene deletions. Hum Mutat 2004; 23:413-9.
35. Ventura M, Mudge JM, Palumbo V, Burn S, Blennow E, Pierluigi M, Giorda R, Zuffardi O, Archidiacono N, Jackson MS, Rocchi M. Neocentromeres in 15q24-26 map to duplicons which flanked an ancestral centromere in 15q25. Genome Res 2003; 13:2059-68.
36. Zhou JY, Taguchi T, Siegfried JM, Jhanwar SC, Resau J, Testa JR. Characterization of 9q;15q whole-arm translocation derivatives in non-small cell lung carcinomas by fluorescence in situ hybridization. Cancer Genet Cytogenet 1993; 69:1-6.
37. Shohat M, Shohat T, Rimoin DL, Mohandas T, Heckenlively J, Magenis RE, Davidson MB, Korenberg JR. Rearrangement of chromosome 15 in the region q11.2----q12 in an individual with obesity syndrome and her normal mother. Am J Med Genet 1990; 37:173-7.
38. Guanti G. A new case of rearrangement of chromosome 15 associated with Prader Willi syndrome. Clin Genet 1980; 17:423-7.
39. Schluth C, Mattei MG, Mignon-Ravix C, Salman S, Alembik Y, Willig J, Ginglinger E, Jeandidier E. Intrachromosomal triplication for the distal part of chromosome 15q. Am J Med Genet A 2005; 136:179-84.
40. Pujana MA, Nadal M, Guitart M, Armengol L, Gratacos M, Estivill X. Human chromosome 15q11-q14 regions of rearrangements contain clusters of LCR15 duplicons. Eur J Hum Genet 2002; 10:26-35.
41. Freeman SB, May KM, Pettay D, Fernhoff PM, Hassold TJ. Paternal uniparental disomy in a child with a balanced 15;15 translocation and Angelman syndrome. Am J Med Genet 1993; 45:625-30.