Article Text
Statistics from Altmetric.com
Keratoconus (OMIM148300) is a bilateral, non-inflammatory, slowly progressive, corneal ectasia that is a major cause of corneal transplant. Characteristically, the cornea becomes thin and conical, with myopia and irregular astigmatism that leads to vision impairment. The incidence of keratoconus is between 50 and 230 per 100 000, with remarkable differences between ethnic groups.1 Although the pathogenesis of keratoconus still is unknown, little doubt exists about an underlying genetic background. A positive family history is found in 6–8% of patients with keratoconus, and concordance is high among monozygotic twins.1,2 In rare instances, keratoconus is inherited either as a mendelian trait or is associated with a genetic disorder (for example, Down syndrome, Leber’s amaurosis, or connective tissue diseases). For most patients, however, the disease is sporadic.3 The lack of multiple familial cases is a major obstacle to linkage analysis. To date, three loci have been associated with keratoconus, none of them through a genomewide search in a single informative family. A putative 6.8 cM locus was identified in a family affected by keratoconus by direct scan of chromosome 21 only, while an association at 20q12 was found in seven related Tasmanian patients.4,5 More recently, a non-parametric linkage analysis in 20 Finnish small pedigrees with autosomal dominant keratoconus identified a third locus to chromosome 16q.6 Finally, two distinct heterozygous mutations in the VSX1 homeobox gene were identified with a candidate gene approach in two families affected by keratoconus.7
Recent advances in computerised topographic diagnostic techniques enable higher accuracy in the diagnosis of keratoconus, including forme fruste, which eases the identification of extensive families affected by keratoconus. We describe an Italian family with autosomal dominant pure keratoconus and the localisation of a novel keratoconus locus to chromosome 3p14–q13.
METHODS
Clinical studies
We identified an Italian pedigree with autosomal dominant keratoconus spanning two generations. Fourteen family members underwent ophthalmological and slit lamp examination, and all but one (II:9) had corneal topography (Keratron; Opticon 2000, Rome, Italy). A diagnosis of definite keratoconus was established when central or paracentral corneal thinning or conical protrusion of the cornea with or without stromal scar tissue, or both, was found, in association with a videokeratographic picture typical of keratoconus. Patients who previously had undergone penetrating corneal transplant also received a diagnosis of definite keratoconus. A diagnosis of latent or forme fruste keratoconus was established when no signs of keratoconus were found on slit lamp examination but when corneal topography showed abnormal findings, such as central steepening, inferior–superior asymmetry, differences between right and left central corneal power, and skewing of the steepest radial axes above and below the horizontal meridian.8
Key points
-
An Italian pedigree with 11 members in two generations with autosomal dominant keratoconus (OMIM 148300) was investigated.
-
Linkage to all known keratoconus loci on chromosomes 16, 20, and 21 was excluded. A genomewide search allowed mapping of a novel keratoconus locus to chromosome 3p14–q13, with a maximum logarithm of odds score of 3.09. Analysis of key recombinant members placed the keratoconus locus to a 43.7 cM region flanked by markers D3S1600 and D3S1278. This spans about 53 Mb on the physical map and includes about 4.5 Mb centromeric DNA.
-
Mutation analysis of COL8A1, an interesting candidate gene located within the linked region, did not show any pathogenic mutations.
-
This is the first locus for autosomal dominant pure keratoconus identified by means of a genomewide search performed in a single informative family.
Linkage analysis
After informed consent, 15 family members and six spouses were blood sampled. Genomic DNA was extracted with standard procedures. Linkage of the disease with known keratoconus loci was excluded by genotyping microsatellite markers that spanned chromosome 21, 20q12, 16q, and 20p11–q11 loci. For the genomewide search, 380 markers that covered the 22 autosomes with an average distance of 10 cM were amplified by polymerase chain reaction in a Perkin Elmer 9700 thermal cycler (Perkin Elmer, Wellesley, MA, USA) with panels of fluorescently labelled primers (ABI Prism Linkage Mapping Set version 2; Applied Biosystems, Foster City, CA, USA). Polymerase chain reactions were performed in microtitre plates in a final volume of 10 μl containing 50 ng genomic DNA, 5 pmol of each primer, 2 mM each of dATP, dCTP, dGTP, and dTTP, 1 μl GeneAmp 10× magnesium free PCR Gold buffer (Applied Biosystems, Foster City, CA, USA), 2.5 mM magnesium chloride, and 0.5 units AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA). Polymerase chain reaction conditions for these markers already have been optimised by the manufacturers, and all markers could be amplified under the same conditions. Reaction mixes were first heated to 95°C for 11 minutes; subsequent cycling conditions were 94°C for 15 seconds, 55°C for 15 seconds, and 72°C for 30 seconds (repeated for 10 cycles); 89°C for 15 seconds, 55°C for 15 seconds, and 72°C for 30 seconds (repeated for 20 cycles); and a final extension of 72°C for 10 minutes. After amplification, PCR products from each panel were pooled in a tube with Genescan 500HD size standard (Applied Biosystems, Foster City, CA, USA) and loaded on a 3100 automated DNA sequencer (Applied Biosystems, Foster City, CA, USA). Analysis of the amplified fragments was performed with Genescan and Genotyper software (Applied Biosystems, Foster City, CA, USA). Two point logarithm of odds (LOD) scores were generated with the Fastlink version of the Mlink program;9 we assumed an equal man–woman recombination rate, autosomal dominant inheritance, a gene frequency of 0.0001, and reduced penetrance (0.80). Allele frequencies were estimated based on a random sample of 16 chromosomes (12 from the six healthy spouses and four from the two people in generation I). Multipoint LOD scores were obtained with Simwalk2 software.10
We constructed haplotypes manually and assigned phase on the basis of the smallest number of recombinants. To saturate the linked region, we selected five novel polymorphic markers from the Human Genome Working Draft with Tandem Repeat Finder software.11 We took genetic distances between markers from the Marshfield genetic map (www.marshfieldclinic.org/research/genetics/default.htm), while physical distances were taken from the University of California Santa Cruz draft of the human genome (www.genome.ucsc.edu) (table 1). Polymerase chain reactions were prepared in a volume of 10 μl as described above. After heating for 11 minutes at 95°C, subsequent cycling conditions were 94°C for 30 seconds, 52–60°C for 30 seconds, and 72°C for 30 seconds (repeated for 28–33 cycles). We optimised the annealing temperature and number of cycles to produce an approximately constant PCR yield. Primer pairs and PCR conditions for these markers are available on request.
Pairwise LOD scores between keratoconus and markers on chromosome 3p14–q13
Candidate gene analysis
The five exons of the alpha1 type VIII collagen (COL8A1; RefSeq NM_001850) were amplified from genomic DNA under the following conditions. For PCR, 250 ng of genomic DNA was used as a template with 10× PCR buffer, 2.5 mM magnesium chloride, 10 mM of nucleotides (dNTP), 50 pmol of each primer, and 1.25 units of Taq Promega Polymerase (Promega, Madison, WI, USA) in a total volume of 50 μl. Primer pairs are available on request. The coding regions and exon–intron boundaries were analysed by direct sequencing of purified PCR products with an ABI BigDye terminator Sequencing Kit and ABI 3100 Capillary Array Sequencer (Applied Biosystem, Foster City, CA). Sequences were analysed with Sequencing Analysis and SeqScape software (Applied Biosystem, Foster City, CA, USA).
RESULTS
Clinical examination
The mode of inheritance was autosomal dominant, with male to male transmission. Seven people received a diagnosis of definite keratoconus, while four were diagnosed with forme fruste keratoconus. A simplified family tree is presented in fig 1; table 2 summarises the clinical features of the 11 affected family members.
Clinical data of family members affected by keratoconus
Pedigree of family and haplotypes of markers spanning the novel keratoconus locus. ▪ affected man; □ unaffected man; • affected woman; ○ unaffected woman; / deceased; ? unknown phenotype;1 ▒ haplotype that segregated with disease in the family.
Individual III:6 underwent bilateral corneal transplant at age 30 years for advanced keratoconus. She had complained of progressive visual loss since the age of 18 years, with photophobia and excessive tearing. In six people (II:5, II:9, II:11, III:1, III:5, and III:8), a diagnosis of keratoconus was established on clinical and slit lamp examination. Videokeratography was subsequently performed in all but individual II:9 and confirmed the diagnosis (fig 2). All these patients presented with photophobia or excessive tearing, sometimes associated with vision impairment, except individual II:5, who did not have any ocular symptoms. The age at onset varied among affected individuals, ranging from 18 to 42 years (table 2). Patient II:5 complained of visus reduction that could be corrected with lenses, but could not report the age of onset. Four individuals received a diagnosis of forme fruste keratoconus on the basis of videokeratography (III:2, III:4, III:9, and III:10) (fig 2). One of them (III:2) reported mild photophobia. Individual II:1 was not available for slit lamp examination or videokeratography examination. His clinical history was characterised by very mild ocular symptoms (that is, mild photophobia) since the age of 60 years, but he could not be examined, as he was in a wheelchair because of limb girdle muscular dystrophy. This man was considered as “unknown phenotype” for linkage purposes. Three family members (II:7, III:3, and III:7) did not complain of any ocular symptom, and all ophthalmologic examinations were normal (fig 2).
Videokeratography showing: normal corneal presentation in individual III:7 (top panel), forme fruste keratoconus in patient III:9 (middle panel); and keratoconus in patient III:8 (bottom panel).
Linkage and molecular genetic analysis
Linkage to known keratoconus loci (on chromosomes 21, 20q12, and 16q) and to the VSX1 gene on 20p11–q11 was excluded (data not shown). The genomewide search produced negative or non-significant LOD scores for all but 13 microsatellite markers on chromosomes 1, 2, 3, 6, 7, 11, and 13. The regions that surrounded these loci and regions that surrounded non-informative markers were saturated with more markers, and haplotypes were manually constructed. The segregation of different haplotypes in affected individuals and the negative LOD scores obtained allowed exclusion of all loci except the pericentromeric region of chromosome 3, in which a common haplotype was shared by the 11 affected family members. The upper and lower boundaries were determined by recombination events in individuals II:1 and III:2, which identified a 43.7 cM linked region flanked by markers D3S1600 and D3S1278. this region spanned about 53 Mb and included about 4.5 Mb of centromeric DNA. All genotyped markers in the region generated positive LOD scores, with a maximum value of 3.09 for marker Chr3:77.08 (θ = 0; penetrance = 0.80). Logarithm of odds scores and intermarker distances are shown in table 1; fig 3 shows the multipoint graph.
Multipoint LOD scores between the disease and markers on chromosome 3p14–q13.
In silico analysis of the linked region identified more than 100 genes. An interesting candidate gene was COL8A1, which codes for the human alpha1(VIII) chain of type VIII collagen, highly expressed in the cornea. Sequencing of the whole COL8A1 coding region and exon–intron boundaries did not show any pathogenic mutation. A single nucleotide, G to T substitution, not previously described, was found at position −59 in the untranslated exon 3 of the gene. Nine affected people and three spouses (II:4, II:8, and II:10) were heterozygous for the G to T change, but individuals III:4 and III:9 were homozygous. This substitution therefore was considered as a polymorphic variant.
DISCUSSION
We report a clinical and molecular investigation of a large central Italian kindred affected by pure keratoconus and the mapping of a novel keratoconus locus to the pericentromeric region of chromosome 3.
The phenotype mainly was characterised by slowly progressive visual impairment and ocular symptoms, such as photophobia and tearing. The age of onset, clinical presentation, and disease progression varied widely among the 11 affected family members (table 2). For instance, individual III:6 underwent corneal transplant because of severe keratoconus at age 30 years, while patient III:4 was still asymptomatic at age 46 and received a diagnosis of forme fruste keratoconus on videokeratography. A true lack of penetrance was also observed, as individual III:3, aged 39 years, shared the disease haplotype but the diagnosis of keratoconus was excluded by slit lamp examination and videokeratography.
Reduced penetrance and variable expression of autosomal dominant keratoconus have been reported in several pedigrees.4–7,12 However, it must be noted that the detection rate for keratoconus depends strictly on the sensitivity of the diagnostic techniques used.1,3 In this family, 10 patients were evaluated by both slit lamp examination and videokeratography: four (40%) received a diagnosis of forme fruste keratoconus only by videokeratography, while slit lamp examination was normal. These data confirm that abortive forms of keratoconus are common and can be missed without a thorough examination.
The high rate of forme fruste keratoconus, as well as the incomplete penetrance, can partly explain the paucity of large keratoconus kindred in the literature. In fact, although a number of putative keratoconus loci have been identified so far,1,4–6 this is the first locus mapped by a genomewide search in a single informative family.
We obtained a maximum two point LOD score of 3.09 for marker Chr3-77083 (at θ = 0), while multipoint linkage analysis generated an LOD score >3.00 in a 43.7 cM interval between markers D3S1600 and D3S1278. Although these LOD scores represent only tentative evidence of linkage, a common haplotype shared by all keratoconus family members was observed only in this region, while the rest of the genome was excluded thoroughly by linkage and haplotype analyses. The gene frequency of autosomal dominant keratoconus is not known, but families with clear mendelian inheritance are much rarer than sporadic cases, which have an incidence of approximately 1 per 2000 in the general population.1 For this reason, a gene frequency of 1 in 10 000 was assumed for linkage calculation in this family.
An interesting gene within the mapped region is COL8A1, which encodes for the human alpha1(VIII) chain of type VIII collagen. A number of collagen genes (that is, collagen types I and III–VIII) have been considered excellent candidates for keratoconus on the basis of their expression pattern in different layers of the cornea.1 In particular, studies of the Descemet’s membrane have shown that type VIII collagen is a major component of the hexagonal lattice of this tissue.13,14 Recently, missense mutations in COL8A2, which encodes the alpha2 chain of type VIII collagen, have been identified in two forms of corneal endothelial dystrophy, underlining its role in influencing the terminal differentiation of corneal endothelial cells.15 We could not detect any pathogenic mutations in the coding sequences and exon–intron boundaries of COL8A1. As collagen VIII is known to exist in vivo as heterotrimers or homotrimers, deletions of the entire gene are unlikely to play a pathogenic role in keratoconus. Intronic and regulative regions of the COL8A1 gene were not analysed, however, so the occurrence of aberrant splicing or expression could not be ruled out.
The identification of additional autosomal dominant families linked to this locus will help evaluate its role among familial keratoconus, refine the linked region, and ease the identification of positional candidate genes.
REFERENCES
Footnotes
-
Funding: This work was supported by grants from the Italian Ministry of Health, “Ricerca Corrente 2003” and the Italian Ministry of Education, University and Research, “Progetto Facoltà 2002”.
-
Conflicts of interest: None declared.