- Case report
- Open access
- Published:
Case report: a novel KERA mutation associated with cornea plana and its predicted effect on protein function
BMC Medical Genetics volume 16, Article number: 40 (2015)
Abstract
Background
Cornea plana (CNA) is a hereditary congenital abnormality of the cornea characterized by reduced corneal curvature, extreme hypermetropia, corneal clouding and hazy corneal limbus. The recessive form, CNA2, is associated with homozygous or compound heterozygous mutations of the keratocan gene (KERA) on chromosome 12q22. To date, only nine different disease-associated KERA mutations, including four missense mutations, have been described.
Case presentation
In this report, we present clinical data from a Turkish family with autosomal recessive cornea plana. In some of the affected individuals, hypotrichosis was found. KERA was screened for mutations using Sanger sequencing. We detected a novel KERA variant, p.(Ile225Thr), that segregates with the disease in the homozygous form. The three-dimensional structure of keratocan protein was modelled, and we showed that this missense variation is predicted to destabilize the structure of keratocan, leading to the classical ocular phenotype in the affected individuals. All the four known missense mutations, including the variation found in this family, affect the conserved residues of the leucine rich repeat domain of keratocan. These mutations are predicted to result in destabilization of the protein.
Conclusion
We present the 10th pathogenic KERA mutation identified so far. Protein modelling is a useful tool in predicting the effect of missense mutations. This case underline the importance of the leucin rich repeat domain for the protein function, and this knowledge will ease the interpretation of future findings of mutations in these areas in other families with cornea plana.
Background
Cornea plana is a rare hereditary congenital abnormality of the cornea characterized by reduced corneal curvature, extreme hypermetropia, corneal clouding and hazy corneal limbus. There is an autosomal dominant (CNA1; OMIM 121400) and a more severe autosomal recessive form (CNA2; OMIM 217300) of the disorder. The recessive form is more frequently associated with additional ocular manifestations, such as iris malformation and adhesions between cornea and iris [1–4]. The gene defective in CNA1 has not been identified yet, while homozygous or compound heterozygous mutations in the keratocan gene (KERA) on chromosome 12q22 are shown to cause CNA2 [5]. The protein encoded by KERA, keratocan, is a keratan sulfate proteoglycan with a core of leucine-rich repeats (LRR), flanked by clusters of cystein, essential to maintain the three-dimensional structure of the protein. Keratocan is an extracellular matrix molecule of the corneal stroma that is important for normal cornea morphogenesis and involved in maintaining corneal transparency [6]. Studies with mice show expression of this gene during eye development and in adult mice, it is exclusively expressed in the cornea [7]. To date, only nine different disease-associated KERA mutations have been reported in various populations [5, 8–11]: Three nonsense mutations, four missense mutations, one splice-site mutation, and one frameshift (single-nucleotide deletion) mutation. All the reported missense mutations are within the LRR domain, which is considered to be essential for the binding between keratocan and the collagen fibrils in the corneal matrix during development. A disruption of the protein affects both the structure and the transparency of the cornea [12].
In this report, we present the clinical data and describe a new KERA mutation in a Turkish family with autosomal recessive cornea plana. We suggest the genetic variation to result in destabilization of the protein structure. We furthermore apply protein modelling to predict the effect of the missense mutations on the structure of keratocan.
Case presentation
Clinical studies
A pedigree of the family is shown in Fig. 1a. The index individuals (IV:4 and IV:5) are dizygotic male twins, who were 11 years old at the time of clinical examination. They were born at gestational age 36 to healthy and unrelated parents. Birth weights were 2300 g (IV:4) and 2352 g (IV:5). Shortly after birth, IV:4 and IV:5 were found to be excessive hypermetropic (+11,75D/+11,75D and 12,75D/+13,75D) and both had corneal clouding and flat corneas. The boys have been wearing glasses since the age of 3 months, and have been monitored frequently. The clouding of the corneas decreased over the first 8 months of life. Clinical and eye examinations were carried out at the age of 10 years showing hazy corneal limbi with vessels were noticed (Fig. 1c). IV:4 has a small temporal iridocorneal synechi. Axial lengths were within normal range (A-scan using Hiscan opticon). ERG, motility, pupillary reaction, funduscopy and intraocular pressure were normal. Keratometry (Zeiss keratometer) showed bilateral flat corneas and high astigmatism in both boys (Additional file 1: Table S1).
The boys’ general development has been in the lower range of normal. They attend a normal school at normal age level, but perform in the low end of their class. They both have slight behavioural problems and motor tics, which are also found in their unaffected younger brother. From the age of 9 years IV:4 has experienced progressive frontal hair loss on vertex, the hair being lighter in colour and lusterless (Fig. 1d). Hair microscopy was normal without any signs of hair shaft disorder. IV:5 has slight general hypermobility and pes planus.
Family history
There were several affected individuals in the family, and the pedigree was suggestive of autosomal recessive inheritance (Fig. 1a). Clinical information of the individuals III:10, III:11, IV:3, IV:9, V:1 and V:2 was extracted from charts of previous examinations. For individuals III:12, IV:1, IV:2 and IV:8 clinical information was obtained from family members. Clinical details on family members are found in Table 1, Fig. 1b and Additional file 1: Table S1.
Molecular studies
The clinical features and family history prompted us to sequence KERA. Genomic DNA was isolated from peripheral blood lymphocytes using Chemagnic Magnetic Separation Modile 1 (Chemagen, Baesweiler, Germany). Coding exons (exon 2 and 3) of KERA were amplified from genomic DNA of the investigated family members using PCR. Purified PCR products were sequenced using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). Sequencing products were purified and run on an ABI 3130 XL genetic analyzer (AppliedBiosystems). Sequences of the coding exons and 10 bp of the flanking introns were analyzed using SeqScape software (AppliedBiosystems). Primer sequences and PCR conditions are available on request.
A novel missense mutation was detected in exon 2 of KERA, c.674T>C, which is predicted to result in substitution of isoleucine at position 225 with threonine, p.(Ile225Thr). All affected family members investigated for the mutation (IV:3, IV:4 and IV:5) were homozygous, while the clinically unaffected parents (III: 3 and III:4), the brother (IV:6), a maternal aunt (III:7) and a maternal uncle (III:5) were all heterozygous for the mutation. . The mutation is absent from the 1000 Human Genome [13] and NHLBI Exome Sequencing [14] databases. The Ile225 lies within the LRR domain and is a highly conserved amino acid. The in silico prediction programs SIFT [15] and Polyphen2 [16] predict the substitution to be pathogenic (SIFT score 0.00 and Polyphen2 score 0.999).
Protein modelling
We modelled the three-dimensional structure of the keratocan protein to predict the effect of the present mutation and the three previously published missense mutations on protein structure/function.
The LRR protein core of decorin (PDB entry 1xku_A) [17] was used as template. The template was chosen using BLASTp[18] to search against the protein data bank. ClustalW [19] was used to align the two sequences and SWISS-MODEL [20] was used to change the amino acids. Finally, the structure was investigated using Pymol [21]. The number of LRR repeats were predicted using LRRfinder [22]. Keratocan has a LRR structure, where the leucine rich motif is repeated several times. It has the shape of a curved solenoid, where the concave side of the structure is made up of a parallel beta-sheet and the convex side is made up of a more diverse number of structural elements (Fig. 2d and e). The leucine-rich motif LXXLXLXXNXL at the concave side of the tandem is highly conserved: The L may be leucine or a hydrophobic amino acid, N may be asparagine or cysteine, and X may be any amino acid (Fig. 2c). The interior of the LRR domain is made up of hydrophobic residues – with the exception of the conserved asparagines (Asn, N) which create an Asn-ladder stabilizing the structure as seen in Fig. 2e [23–25]. Previously, 10 LRR repeats were suggested for keratocan [5], but the LRRfinder [22] suggests presence of an 11th LRR repeat starting at amino acid position 322 (Fig. 2c).
We found that the p.(Asn131Asp) and p.(Asn247Ser) mutations affect the conserved Asn131 and Asn247 residues of the 3rd and 7th LRR-repeat. The substitution of an Asn with an Aspartic acid or a Serine will change a hydrogen bond donor to a hydrogen bond acceptor or to neutral. The important stabilizing structure of the Asn-ladder, which is seen in Fig. 2e, is likely to be destroyed.
The p.(Thr215Lys) mutation affects the Thr215 in the beginning of the 7th repeat where the threonine is pointing into the structure. Substitution of this residue with lysine is likely to destroy the structure in this area due to steric hindrance of the large lysine side chain.
The p.(Ile107Arg) and the current mutation, p.(Ile225Thr) affects the 11th residue within the 2nd and the 7th repeat (the end of the beta-strain) where the isoleucine is pointing into the structure, stabilizing the hydrophobic core. Substitution of the hydrophobic isoleucine with the hydrophilic and bulky arginine or the hydrophilic threonine will most likely destabilize the structure. It is interesting that the two mutations p.(Ile107Arg) and p.(Asn131Asp) have the same special orientation with respect to each other as the two mutations p.(Ile225Thr) and p.(Asn247Ser).
Discussion
This family, with a novel missense mutation KERA [p.(Ile225Thr)], exhibits the classical ocular phenotype of autosomal recessive cornea plana, with excessive hypermetropia, hazy corneal limbus, and (although variable) iris malformation and iridocorneal adhesions.
Only nine different disease-associated mutations have been previously reported in KERA [4;5;8–11] (Fig. 2a, b, c). There are several lines of evidence suggesting that p.(Ile225Thr) is causative: The mutation segregates with the disease in the family and it affects a conserved residue located within the LRR domain which is predicted to be important for protein function. In silico analysis predicts the mutation to be damaging and protein modelling suggests that substitution of isoleucine at position 225 of the 7th LRR repeat will result in destabilization of the protein structure. We therefore consider this mutation the most likely cause of cornea plana in this family.
Three affected family members were homozygous for the mutation. The parents of patient IV:3 are consanguineous. The parents of the index patients IV:4 and IV:5 are not, to their own knowledge, consanguineous, but descend from the same area. Since they are both heterozygous for the same mutation, we speculate that they are distantly related. If the mutation is found in other Turkish families, however, a founder effect could be considered.
Several family members have sparse scalp hair, and one of the twins exhibits severe frontal hypotrichosis. Although keratocan has not been linked to abnormalities of the hair, another protein member of the small LRR proteoglycan family, called Tsukushi (TSK), is found to be expressed at the restricted areas of hair follicles during the morphogenesis, and targeted disruption of TSK causes the hair cycle to be delayed in mice [17]. However, the expression pattern of keratocan in the hair follicle has not been examined. Hair loss has not previously been described in connection with cornea plana. Pedigree analysis does not show a clear connection between hair loss and the mutation - family members not homozygous for the mutation also show some thinning of the hair, which might be androgenetic. If more cases of hair loss in cornea plana patients emerge, expression studies of keratocan in the hair follicle should be considered.
Conclusions
KERA mutations are very rare and we present the 10th pathogenic mutation identified so far. Protein modelling is a useful tool in predicting the effect of missense mutations. Our studies suggest that the five amino acid substitutions associated with cornea plana are likely to exercise their effect by destabilizing the structure of the keratocan protein, underlining the importance of this domain for the protein function. This is valuable knowledge when interpreting the clinical significance of future KERA mutations.
Consent
The project conformed to the tenets of the Declaration of Helsinki and was approved by the Regional Scientific Ethical Committee for Copenhagen, Denmark. Written informed consent was obtained from the patient for publication of this Case report and any accompanying images. A copy of the written consent is available for review by the Series Editor of this journal.
References
Khan AO, Aldahmesh M, Meyer B. Corneal ectasia and hydrops in a patient with autosomal recessive cornea plana. Ophthalmic Genet. 2006;27(3):99–101.
Khan AO, Aldahmesh MA, Al-Gehedan S, Meyer BF, Alkuraya FS. Corneal decompensation in recessive cornea plana. Ophthalmic Genet. 2009;30(3):142–5.
Tahvanainen E, Forsius H, Kolehmainen J, Damsten M, Fellman J, de la Chapelle A. The genetics of cornea plana congenita. J Med Genet. 1996;33(2):116–9.
Khan AO, Aldahmesh M, Meyer B. Recessive cornea plana in the Kingdom of Saudi Arabia. Ophthalmology. 2006;113(10):1773–178.
Pellegata NS, Dieguez-Lucena JL, Joensuu T, Lau S, Montgomery KT, Krahe R, et al. Mutations in KERA, encoding keratocan, cause cornea plana. Nat Genet. 2000;25(1):91–5.
Tasheva ES, Funderburgh JL, Funderburgh ML, Corpuz LM, Conrad GW. Structure and sequence of the gene encoding human keratocan. DNA Seq. 1999;10(1):67–74.
Liu CY, Shiraishi A, Kao CW, Converse RL, Funderburgh JL, Corpuz LM, et al. The cloning of mouse keratocan cDNA and genomic DNA and the characterization of its expression during eye development. J Biol Chem. 1998;273(35):22584–8.
Ebenezer ND, Patel CB, Hariprasad SM, Chen LL, Patel RJ, Hardcastle AJ, et al. Clinical and molecular characterization of a family with autosomal recessive cornea plana. Arch Ophthalmol. 2005;123(9):1248–53.
Khan A, Al-Saif A, Kambouris M. A novel KERA mutation associated with autosomal recessive cornea plana. Ophthalmic Genet. 2004;25(2):147–52.
Liskova P, Hysi PG, Williams D, Ainsworth JR, Shah S, de la Chapelle A, et al. Study of p.N247S KERA mutation in a British family with cornea plana. Mol Vis. 2007;13:1339–47.
Dudakova L, Palos M, Hardcastle AJ, Liskova P. Corneal endothelial findings in a Czech patient with compound heterozygous mutations in KERA. Ophthalmic Genet. 2014;35(4):252–4.
Hassell JR, Birk DE. The molecular basis of corneal transparency. Exp Eye Res. 2010;91(3):326–35.
1000 Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, Depristo MA, Durbin RM, et al. An integreated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.
Tennessen JA, Bigham AW, O’connor TD, Fu W, Kenny EE, Gravel S, et al. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2010;337(6090):64–9.
Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11(5):863–74.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Scott PG, McEwan PA, Dodd CM, Bergmann EM, Bishop PN, Bella J. Crystal structure of the dimeric protein core of decorin, the archetypal small leucine-rich repeat proteoglycan. Proc Natl Acad Sci U S A. 2004;101(44):15633–8.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;1;25(17):3389–402.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.
Arnold K, Bordoli L, Kopp J, Schwede T. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;15;22(2):195–201.
The PyMOL Molecular Graphics System [computer program]. Version Version 1.5.0.4 Schrödinger, New York, USA: LLC; 2014.
Offord V, Coffey TJ, Werling D. LRRfinder: a web application for the identification of leucine-rich repeats and an integrative Toll-like receptor database. Dev Comp Immunol. 2010;34(10):1035–41.
Bella J, Hindle KL, McEwan PA, Lovell SC. The leucine-rich repeat structure. Cell Mol Life Sci. 2008;65(15):2307–33.
Kobe B, Deisenhofer J. The leucine-rich repeat: a versatile binding motif. Trends Biochem Sci. 1994;19(10):415–21.
Kobe B, Deisenhofer J. Proteins with leucine-rich repeats. Curr Opin Struct Biol. 1995;5(3):409–16.
Acknowledgements
We would like to thank the family for their participation in this study. BB is currently supported by a scholarship from the University of Copenhagen.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LR conceived of the study, performed clinical examination of the patients and drafted the manuscript. HJ performed eye examination on the patients. AB performed dermatological examination and hair microscopy of the patients. BB carried out the molecular genetic studies, and helped to draft the manuscript. PH modelled and analysed the protein structure. ZT and KG participated in the design of the study, and helped to draft the manuscript. All authors read and approved the final manuscript.
Additional file
Additional file 1: Table S1.
Best corrected visual acuity, biometric characteristics and refractive error of affected family members. Topography of patient IV:4.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Roos, L., Bertelsen, B., Harris, P. et al. Case report: a novel KERA mutation associated with cornea plana and its predicted effect on protein function. BMC Med Genet 16, 40 (2015). https://doi.org/10.1186/s12881-015-0179-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12881-015-0179-9