Purpose Keratoconus (KTCN) is a noninflammatory thinning and anterior protrusion of the cornea that results in steepening and distortion of the cornea, altered refractive powers, and reduced visual acuity. Several loci responsible for a familial form of KTCN have been mapped, however; no mutations in any genes have been identified for any of these loci. There is also evidence that VSX1 and SOD1 may be involved in the etiology of KTCN. The purpose of this study was to verify the available data and to identify a new keratoconus susceptibility locus. Methods KTCN without other ocular or systemic features was diagnosed in 18 families. VSX1 and SOD1 sequencing was performed on affected individuals and control subjects. Genomewide linkage analysis was then performed in all families using polymorphic microsatellite markers with an average spacing of 5 cM. Next, single-nucleotide polymorphism (SNP) arrays, fluorescence in situ hybridization (FISH) analysis, and a comparative genomic hybridization array were used in one family to assess a candidate region on 13q32. Results All previously reported KTCN loci were excluded. VSX1 and SOD1 were sequenced, and no potentially functional variants were found. One KTCN family yielded a maximum multipoint parametric LOD score of 4.1 and multipoint non-parametric linkage (NPL) LOD score of 3.2. Multipoint linkage and haplotype analysis narrowed the locus to a 5.6-Mb region between the SNPs rs9516572 and rs3825523 on 13q32. Conclusions The results exclude VSX1 and SOD1 as potential disease-causing genes in these families and localize a novel gene for keratoconus to a 5.6-Mb interval on 13q32.

Keratoconus (KTCN), a non-inflammatory corneal disorder characterized by stromal thinning, represents a major cause of corneal transplantations. Genetic and environmental factors have a role in the etiology of this complex disease. Previously reported linkage analysis revealed that chromosomal region 13q32 is likely to contain causative gene(s) for familial KTCN. Consequently, we have chosen eight positional candidate genes in this region: MBNL1, IPO5, FARP1, RNF113B, STK24, DOCK9, ZIC5 and ZIC2, and sequenced all of them in 51 individuals from Ecuadorian KTCN families and 105 matching controls. The mutation screening identified one mutation and three sequence variants showing 100% segregation under a dominant model with KTCN phenotype in one large Ecuadorian family. These substitutions were found in three different genes: c.2262A4C (p.Gln754His) and c.720+43A4G in DOCK9; c.2377-132A4C in IPO5 and c.1053+29G4C in STK24. PolyPhen analyses predicted that c.2262A4C (Gln754His) is possibly damaging for the protein function and structure. Our results suggest that c.2262A4C (p.Gln754His) mutation in DOCK9 may contribute to the KTCN phenotype in the large KTCN-014 family.

Although copy number changes of 5q31 have been rarely reported, deletions have been associated with some common characteristics, such as short stature, failure to thrive, developmental delay (DD)/intellectual disability (ID), club feet, dislocated hips, and dysmorphic features. We report on three individuals with deletions and two individuals with duplications at 5q31, ranging from 3.6 Mb to 8.1 Mb and 830 kb to 3.4 Mb in size, respectively. All five copy number changes are apparently de novo and involve several genes that are important in developmental pathways, including PITX1, SMAD5, and WNT8A. The individuals with deletions have characteristic features including DD, short stature, club feet, cleft or high palate, dysmorphic features, and skeletal anomalies. Haploinsufficiency of PITX1, a transcription factor important for limb development, is likely the cause for the club feet, skeletal anomalies, and cleft/high palate, while additional genes, including SMAD5 and WNT8A, may also contribute to additional phenotypic features. Two patients with deletions also presented with corneal anomalies. To identify a causative gene for the corneal anomalies, we sequenced candidate genes in a family with apparent autosomal dominant keratoconus with suggestive linkage to 5q31, but no mutations in candidate genes were found. The duplications are smaller than the deletions, and the patients with duplications have nonspecific features. Although development is likely affected by increased dosage of the genes in the region, the developmental disruption appears less severe than that seen with deletion. © 2011 Wiley‐Liss, Inc.

The analyses of selected genes have led to identification of numerous sequence variants in the examined Ecuadorian family. Both substitution c.214+242C > T in IL1RN and novel deletion c.2558+149_2558+203del54 in SLC4A11 were observed significantly more frequently in family members with KTCN (P = 0.004525 and P = 0.00761, respectively), suggesting involvement of these two genes in KTCN etiology in the studied family.

Keratoconus (KTCN) is a protrusion and thinning of the cornea, resulting in impairment of visual function. The extreme genetic heterogeneity makes it difficult to discover factors unambiguously influencing the KTCN phenotype. In this study, we used whole-exome sequencing (WES) and Sanger sequencing to reduce the number of candidate genes at the 5q31.1-q35.3 locus and to prioritize other potentially relevant variants in an Ecuadorian family with KTCN. We applied WES in two affected KTCN individuals from the Ecuadorian family that showed a suggestive linkage between the KTCN phenotype and the 5q31.1-q35.3 locus. Putative variants identified by WES were further evaluated in this family using Sanger sequencing. Exome capture discovered a total of 173 rare (minor allele frequency o0.001 in control population) nonsynonymous variants in both affected individuals. Among them, 16 SNVs were selected for further evaluation. Segregation analysis revealed that variants c.475T4G in SKP1, c.671G4A in PROB1, and c.527G4A in IL17B in the 5q31.1-q35.3 linkage region, and c.850G4A in HKDC1 in the 10q22 locus completely segregated with the phenotype in the studied KTCN family. We demonstrate that a combination of various techniques significantly narrowed the studied genomic region and reduced the list of the putative exonic variants. Moreover, since this locus overlapped two other chromosomal regions previously recognized in distinct KTCN studies, our findings suggest that this 5q31.1-q35.3 locus might be linked with KTCN.