Identification of genetic variants in two families with Keratoconus
BMC Medical Genomics volume 16, Article number: 299 (2023)
This research investigated the genetic characteristic of two Chinese families with keratoconus (KC).
For all people in the two families with KC, their history, clinical data, and peripheral blood were collected. One hundred healthy participants without KC and 112 sporadic KC patients were recruited as the controls. Whole exome sequencing of the genomic DNA and polymerase chain reaction were conducted for all the controls and family members to verify the variants. Functional analyses of the variants was performed using the software programs.
A missense tuberous sclerosis 1 (TSC1) variant g.135797247A > G (c.622A > G, p.Ser208Gly) was detected in family 1. A single nucleotide polymorphism (SNP) rs761232139 (p.Gly235Arg) in aldehyde dehydrogenase 3 family member A1 (ALDH3A1) gene was detected in family 2. The variant c.622A > G in TSC1 and the SNP rs761232139 in ALDH3A1 were predicted as being probably damaging.
Novel variant c.622A > G in TSC1 and SNP rs761232139 in ALDH3A1 have been detected in families with KC. These two findings may play a role in the pathogenesis of KC.
Keratoconus (KC) is a corneal disorder characterized by thinning and outward protrusion of the cornea, which results in blurred vision. The etiology of KC includes genetic, behavioral, and environmental factors [1, 2]. Among KC cases, about 5–14% have a family history of KC [3, 4]. This indicates that genetic factors are associated with the pathogenesis of KC. Autosomal dominant and recessive inheritance has been reported in KC with a family history [2, 3]. Bioinformatics studies have revealed the linkage and association of causative genes; however, few mutations have been identified for this disorder [1,2,3,4].
In this study, we used bioinformatics to investigate gene mutations in two Chinese families with KC. A novel variant in tuberous sclerosis complex 1 (TSC1) and a single nucleotide polymorphism (SNP) in aldehyde dehydrogenase 3 family, member A1 (ALDH3A1) have been detected.
Participants and examinations
The subjects participating in this study were 220 people including 8 family members from two Chinese KC families, 100 healthy Chinese individuals without KC, and 112 sporadic KC patients. Among the eight family members from two two-generation Chinese KC families, three were from family 1 and five were from family 2. The 100 unrelated healthy people were normal controls, who were diagnosed with refractive errors without KC and other inherited corneal disorders. The 112 sporadic KC patients (with no significant difference in gender and age compared to the normal controls,χ2 test for gender, P = 0.93; t—test for age, P = 0.36) were also used to verify the SNP identified in the KC members from family 2. Written informed consent of the study was obtained from all subjects. All of the participants received an ophthalmological examination, which included visual acuity testing, slit-lamp and ophthalmoscope examination, and cornea evaluation using the Scheimpflug camera system (Pentacam; Oculus Optikgeräte GmbH, Wetzlar, Germany). This investigation was reviewed and approved by the Institutional Review Board of Fudan University and was performed in compliance with the Declaration of Helsinki (approval number: 2022128).
According to the Chinese Expert Consensus on the Diagnosis and Treatment of Keratoconus (2019) , the clinical diagnostic criteria of KC include progressive myopia, irregular astigmatism, and poor refractive correction with glasses; corneal topography showing corneal thinning and abnormal elevation of the anterior and posterior surfaces of the cornea; and/or decreased biomechanical indexes and corneal resistance factor in a corneal biomechanical examination. The stages of KC are as follows: 1) latent stage: the contralateral eye of a confirmed KC eye, which shows typical corneal topographic features of KC and has normal vision with uncorrected visual acuity (UCVA) ≥ 1.0; 2) incipient stage: KC is confirmed with a best spectacle-CVA (BSCVA) ≥ 0.8; 3) completion stage: KC is confirmed with a BSCVA < 0.8. Typical clinical manifestations of KC are classified into three levels: Grade 1: corneal curvature (CC) < 53.0 D, corneal thickness at the thinnest point (CTTP) > 400 μm, and BSCVA < 0.8; Grade 2: CC < 55.0 D, CTTP > 300 μm, and BSCVA < 0.3; Grade 3: CC > 55.0 D, CTTP ≤ 300 μm, and BSCVA < 0.05); and 4) scar stage: scarring of the cornea after acute KC edema has subsided.
Whole exome sequencing
Exome sequencing (ES) was performed for all participants using as previously described . Briefly, genomic DNA was extracted from peripheral blood leukocytes. Then the genomic DNA from the KC subjects was used for ES and enrichment of the exonic sequences. Afterwards, data processing and analyses were conducted. Variants have been reported (e.g., transforming growth factor beta-induced) and those presented in the KC subjects at frequencies ≤ 1% were examined in the 1000 Genomes Project. Only the variants shared by family members II:1 and I:1 in family 1 or II:2, II:3 and I:2 in family 2 were taken into account as candidate mutations.
Variant and SNP validation and cross-species conservation analysis
The variant validation and analyses were performed as previously described . Briefly, the candidate mutations were confirmed using polymerase chain reaction (PCR) and Sanger sequencing in the normal controls. The PCR primers were designed using Primer3. The validation and analyses were conducted according to the NCBI VARIANT, NCBI HomoloGene, and the 1000 Genomes Project databases. Databases such as the X1000g2015aug_eas, ExAC_EAS, and gnomAD_exome_EAS were used to check the allele frequency of the variant. Scoring according to the American College of Medical Genetics and Genomics (ACMG)  was used for interpretation of the variant. The conserved analysis was conducted using Polyphen2 (http://genetics.bwh.harvard.edu/pph2/) and multiple sequence alignment conformational alteration was conducted using the Swissmodel (https://swissmodel.expasy.org). The SNP was checked by using the Amelie database (https://amelie.stanford.edu/) and ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar/variation/2442791/). PCR and sequencing in the normal controls and sporadic patients were conducted to verify the SNP identified in the family 2.
The genograms are presented in Fig. 1 (1a: family 1; 1b: family 2). In family 1, there were two members affected with KC. The proband (II.1) was a 9-year-old boy. He presented with impaired vision, acute keratoconus edema, and rupture of the Descemet’s membrane in the right eye after rubbing the eye 10 days prior. The uncorrected distance vision acuity (UDVA) was 10/50 (right eye) and 20/50 (left eye). His corrected distance vision acuity (CDVA) was 10/50 (right eye) with -3.00 DS/-6.50 DC × 170° and 20/20 (left eye) with -2.00 DS/-1.00 DC × 165° (Fig. 2). The mean posterior elevation of the cornea (PEC) was 166 µm (right eye) and 25 µm (left eye). Another patient from family 1, subject I.1, was a 46-year-old female. Her UDVAs were 10/50 and CDVAs were 20/20 in both eyes, with -3.25 DS/-1.00 DC × 175° and -4.50 DS/-0.50 DC × 10°. The mean PEC was 15 µm in one eye and 16 µm in another. The central corneal thickness (CCT) was 493 and 491 µm. The clinical data of other family members are presented in Figure S1.
In family 2, there were three members afflicted with KC. The proband (II.3) was a 22-year-old male. His UDVAs were 10/50 in both eyes and CDVA was 20/50 (right eye) with -3.50 DS/-4.00 DC × 175° and 20/20 (left eye) with -2.50 DS/-2.00 DC × 5°. The mean PEC and CCT were 61 and 481 µm in the right eye and 20 and 530 µm in the left eye (Fig. 3). Central cornea thinning, inferior steepening, bulging outwards, and irregular astigmatism were observed in both eyes of all three subjects from family 2. The findings of other family members are presented in Figure S2.
The stages of KC subjects in the two families are shown in Table 1. In the 112 sporadic cases, the numbers of patients in different stages were 8 in stage 1, 65 in stage 2, 30 in stage 3, and 9 in stage 4.
Identification and analyses of the novel variant and SNP
A novel variant and an SNP were identified in these two families and not detected in the normal controls.
In family 1, a variant in TSC1 was detected, g.135797247A > G (c.622A > G, p.Ser208Gly), which was a missense variant (Fig. 4a). In family 2, an SNP in ALDH3A1 was detected, rs761232139 (p.Gly235Arg), and was detected in all KC members from family 2 (Fig. 4b). However, there was no significant conformational change of the protein caused by the SNP in ALDH3A1 compared with the wild-type protein (Fig. 5). The SNP rs761232139 was detected in six patients (5.4%, genotype G/A) from the 112 sporadic KC patients. Compared with the healthy controls, the rs761232139 G/A genotype and A allele frequencies (p = 0.019 and 0.018 respectively, χ2 test,) were higher in the sporadic KC patients.
The variant c.622A > G in TSC1 was predicted as probably damaging and highly conserved using software programs (Polyphen2, SIFT, PROVEAN, and Mutation Taster) (Table 2 and Fig. 4c). Moreover, the possible conformational alteration of the protein caused by the p.Ser208Gly variant in TSC1 compared with the wild-type protein was revealed by three-dimensional modeling (Fig. 6). No allele frequency of the variant c.622A > G in TSC1 was available in the databases of the X1000g2015aug_eas, ExAC_EAS, and gnomAD_exome_EAS. The score of ACMG was uncertain significance and the evidence codes were PS4 and PM2, which means that the variant is rare in the normal population, located in a critical functional domain, and pathogenic supporting via multiple in silico analyses . The allele frequency of the SNP rs761232139 in the normal population was 0.0008 (database: gnomAD_exome_EAS), very rare.
The p.Ser208Gly TSC1 variant was screened in the 112 sporadic cases and 100 normal controls. The information is shown in Table 3. Only 1 patient from the 112 sporadic patients had the variant. His clinical data were as follows: 20-year-old male, CCT (415 µm in OD and 462 µm in OS), Kmax (59.6D in OD and 44.1D in OS), PEC (72 µm in OD and 8 µm in OS), CDVA (3/20 in OD and 20/20 in OS), and stage (grade 3 in OD and grade 2 in OS). The SNP rs761232139 in ALDH3A1 was also screened in the 112 sporadic cases and 100 normal controls. Its distribution of genotype and allele frequencies were significantly different between the sporadic patients and normal controls (P < 0.05, Table 4).
KC is a complex disorder caused by multifactorial etiology, such as the interaction among behavioral, environmental, and genetic factors. Several genes associated with the familial form of KC have been studied and few mutations have been reported. Some of the familial forms of KC are related to other genetic disorders such as Down syndrome, connective tissue disorders, and Leber congenital amaurosis [1,2,3,4,5,6, 8]. Investigating the various KC subsets could shed light on its etiology. In this study, two findings that may be associated with a familial KC were identified: a missense variant in TSC1 and an SNP in ALDH3A1.
TSC1 and TSC2 genes were named because some of their mutations are related to a genetic disorder TSC [9,10,11,12]. TSC is a disorder that causes hamartomas in multiple organs (e.g., brain, eyes, skin, kidneys). The pathogenesis of TSC involves changes in cellular proliferation and differentiation. Though rare, the co-existence of TSC and KC has been reported [9,10,11] and a 4-base-pair deletion in TSC2 was identified . In the present investigation, the variant of g.135797247A > G (c.622A > G, p.Ser208Gly) located in exon 7 of TSC1 gene was detected. This variant may cause a conformational change in the protein and has not been previously reported in patients with or without TSC. The residue that this variant involves has been considered a distinct functional and/or structural sequence of the TSC1 gene. Moreover, the variant was predicted as probably damaging and highly conserved among species. In the family in which the novel variant was detected, TSC symptoms were not found in any family member.
TSC1 protein (hamartin) interacts with other proteins including TSC2 protein (tuberin), TBC17D, AKT2, MTOR, and RPS6KA1 to exert its effects [9,10,11,12,13]. The complex of TSC1, TSC2, and TBC17D proteins negatively regulates cellular proliferation and differentiation. The underlying mechanism is a canonical signaling pathway with Ras homologue enriched in the brain (Rheb) and the mammalian target of rapamycin complex 1 (mTORC1) involved [9,10,11,12,13]. AKT2, MTOR, and RPS6KA1 were screened as potential interactors of TSC1 protein using software  and were found to be expressed in ocular tissues . All of these findings indicate that TSC1 protein plays a role in the eye. The alterations in the structure and function of TSC1 protein may change its interaction with other proteins and affect cell growth and proliferation in the eye. This may cause KC, a disorder of corneal stroma hypoplasia , in patients with or without TSC.
ALDH3A1 is thought to be involved in KC etiology [14,15,16]. ALDH3A1 SNPs, rs1042183 and rs2228100, are strongly associated with KC in the Polish and Korean populations [14, 16]. The association of ALDH3A1 with intraocular pressure  and refractive astigmatism [18, 19] have also been observed. Furthermore, the relationship between astigmatism and KC has been reported, as 14.1% of patients with more than 2D astigmatism suffer from KC . Our findings further demonstrated the association of ALDH3A1 with KC. A SNP in ALDH3A1 was detected, rs761232139 (p.Gly235Arg). This SNP is located in exon 6 of ALDH3A1 gene and is predicted as highly conserved among species and probably damaging. The SNP was identified in all KC members from the family 2, and in six sporadic KC patients (5.4%) from the total 112 patients but none of the normal controls, which means the detection rate of this SNP (rs761232139) in KC patients is much higher than that in normal people. Thus, SNP (rs761232139) identified in this study was considered to increase the risk of KC occurrence and inherited in a dominant model. The residue that the SNP involves was considered a distinct structural and/or functional sequence of ALDH3A1, which is adjacent to the active site Cys-243.
ALDH3A1 protein is a very important enzyme that is highly expressed in the corneal epithelium [15, 21]. It protects the cornea from ultraviolet (UV)-induced damage in various aspects such as absorbing UV radiation directly, metabolizing toxic aldehydes produced by UV-induced lipid peroxidation, and eliminating UV-produced reactive oxygen species (ROS) . The association of accumulated ROS and toxic aldehydes with KC has been reported in previous studies [22,23,24,25,26]. Furthermore, the SNP identified in the present study is quite close to the active site Cys-243 in ALDH3A1. This cysteine residue plays an important role in the aldehydic carbonyl group attack, formation of a covalent thiohemiacetal intermediate and NAD(P) + , and interaction with Glu-209 and Glu-333 in ALDH3A1 . These functions are involved in the elimination of ROS accumulation and subsequent damage to the cells, which protects the cornea from UV radiation injury [22, 28].
In conclusion, a novel variant c.622A > G in TSC1 and a SNP rs761232139 in ALDH3A1 have been detected in two Chinese KC families. These two findings may play a role in the pathogenesis of KC. These findings expand the spectrum of identified genetic characteristics in KC. Further studies with a larger sample size are warranted to confirm the findings and reveal the mechanism by which they cause KC.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the NCBI ClinVar repository, ClinVar accession number: SCV003804500 and accessible direct link: https://www.ncbi.nlm.nih.gov/clinvar/variation/2442791/?oq=SCV003804500&m=NM_000691.5(ALDH3A1):c.703G%3EA%20(p.Gly235Arg). ClinVar accession number: SCV003804463 and accessible direct link: https://www.ncbi.nlm.nih.gov/clinvar/variation/2430291/?oq=SCV003804463&m=NM_000368.5(TSC1):c.622A%3EG%20(p.Ser208Gly).
- TSC1 :
Tuberous sclerosis 1
- ALDH3A1 :
Aldehyde dehydrogenase 3 family member A1
Polymerase chain reaction
Rong SS, Ma STU, Yu XT, Ma L, Chu WK, Chan TCY, et al. Genetic associations for keratoconus: a systematic review and meta-analysis. Sci Rep. 2017;7:4620.
Nowak DM, Gajecka M. The genetics of keratoconus. Middle East Afr J Ophthalmol. 2011;18:2–6.
Loukovitis E, Sfakianakis K, Syrmakesi P, Tsotridou E, Orfanidou M, Bakaloudi DR, et al. Genetic aspects of keratoconus: a literature review exploring potential genetic contributions and possible genetic relationships with comorbidities. Ophthalmol Ther. 2018;7:263–92.
Bykhovskaya Y, Margines B, Rabinowitz YS. Genetics in keratoconus: where are we? Eye Vis (Lond). 2016;3:16.
Chinese Medical Association Ophthalmology Branch Corneal Disease Group. Expert consensus on the diagnosis and treatment of keratoconus in China. Chin J Ophthalmol. 2019;55:891–5.
Lin Q, Zheng L, Shen Z. A novel variant in TGFBI causes keratoconus in a two-generation Chinese family. Ophthalmic Genet. 2022;43:159–63.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.
Romero-Jiménez M, Santodomingo-Rubido J, Wolffsohn JS. Keratoconus: a review. Cont Lens Anterior Eye. 2010;33:157–66 quiz 205.
Bykhovskaya Y, Fardaei M, Khaled ML, Nejabat M, Salouti R, Dastsooz H, et al. TSC1 mutations in keratoconus patients with or without tuberous sclerosis. Invest Ophthalmol Vis Sci. 2017;58:6462–9.
Wong R, Rajendram R, Poole TR. Keratoconus in tuberous sclerosis. Eye (Lond). 2007;21:287–8.
Lyall DA, Tobias ES, Srinivasan S, Willoughby C. Bilateral keratoconus in tuberous sclerosis: is there a molecular link? Can J Ophthalmol. 2012;47:e41–2.
Curatolo P, Specchio N, Aronica E. Advances in the genetics and neuropathology of tuberous sclerosis complex: edging closer to targeted therapy. Lancet Neurol. 2022;21:843–56.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.
Mok J, Joo CK. ALDH3A1 variations are strongly associated with susceptibility to keratoconus patients in Korean. Investig Ophthalmol Vis Sci. 2017;58:5658.
Khaled ML, Helwa I, Drewry M, Seremwe M, Estes A, Liu Y. Molecular and histopathological changes associated with keratoconus. Biomed Res Int. 2017;2017:7803029.
Berdyński M, Krawczyk P, Safranow K, Borzemska B, Szaflik JP, Nowakowska-Żawrocka K, et al. Common ALDH3A1 gene variant associated with keratoconus risk in the Polish population. J Clin Med. 2021;11(1):8.
Gao XR, Huang H, Nannini DR, Fan F, Kim H. Genome-wide association analyses identify new loci influencing intraocular pressure. Hum Mol Genet. 2018;27:2205–13.
Shah RL, Guggenheim JA. Genome-wide association studies for corneal and refractive astigmatism in UK Biobank demonstrate a shared role for myopia susceptibility loci. Hum Genet. 2018;137:881–96.
Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47:D1005–12.
Serdarogullari H, Tetikoglu M, Karahan H, Altin F, Elcioglu M. Prevalence of keratoconus and subclinical keratoconus in subjects with astigmatism using pentacam derived parameters. J Ophthalmic Vis Res. 2013;8:213–9.
You J, Corley SM, Wen L, Hodge C, Höllhumer R, Madigan MC, et al. RNA-Seq analysis and comparison of corneal epithelium in keratoconus and myopia patients. Sci Rep. 2018;8:389.
Stagos D, Chen Y, Cantore M, Jester JV, Vasiliou V. Corneal aldehyde dehydrogenases: multiple functions and novel nuclear localization. Brain Res Bull. 2010;81:211–8.
Chen L, Zong R, Zhou J, Ge L, Zhou T, Ma JX, et al. The oxidant role of 4-hydroxynonenal in corneal epithelium. Sci Rep. 2015;5:10630.
Gulpamuk B, Koç M, Karatepe MS, Yildiz A, Erel O, Neselioglu S, et al. Novel assay assessment of oxidative stress biomarkers in patients with keratoconus: thiol-disulfide homeostasis. Curr Eye Res. 2017;42:1215–9.
Kaur J, Kukreja S, Kaur A, Malhotra N, Kaur R. The oxidative stress in cataract patients. J Clin Diagn Res. 2012;6:1629–32.
Chen Y, Thompson DC, Koppaka V, Jester JV, Vasiliou V. Ocular aldehyde dehydrogenases: protection against ultraviolet damage and maintenance of transparency for vision. Prog Retin Eye Res. 2013;33:28–39.
Khanna M, Chen CH, Kimble-Hill A, Parajuli B, Perez-Miller S, Baskaran S, et al. Discovery of a novel class of covalent inhibitor for aldehyde dehydrogenases. J Biol Chem. 2011;286:43486–94.
Álvarez-Barrios A, Álvarez L, García M, Artime E, Pereiro R, González-Iglesias H. Antioxidant Defenses in the Human Eye: A Focus on Metallothioneins. Antioxidants (Basel). 2021;10.
The authors would like to express their gratitude to all of the study participants for their cooperation and to Medjaden Inc. for their scientific editing.
This study was supported in part by the National Natural Science Foundation of China for Young Scholars (Grant No. 82000929), the Shanghai Sailing Program (Grant No. 20YF1405000), the National Natural Science Foundation of China (Grant No. 81770955), the Project of Shanghai Science and Technology (Grant No. 20410710100), the Clinical Research Plan of SHDC (SHDC2020CR1043B), the Project of Shanghai Xuhui District Science and Technology (2020–015), the Project of Shanghai Xuhui District Science and Technology (XHLHGG202104), the Shanghai Engineering Research Center of Laser and Autostereoscopic 3D for Vision Care (20DZ2255000), and the construction of a 3D digital intelligent prevention and control platform for the whole life cycle of highly myopic patients in the Yangtze River Delta (21002411600).
Ethics approval and consent to participate
The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Fudan University (approval number: 2022128 and date of approval: 2022–09-28). Written informed consent of the study was obtained from all subjects.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
a: The corneal topography (Pentacam) reports for I.1 in family 1. b: the Corvis ST report for I.1 in family. Figure S2. a: The corneal topography (Pentacam) reports for I.2 in family 2. b: The corneal topography (Pentacam) reports for II.2 in family 2.
About this article
Cite this article
Lin, Q., Wang, X., Han, T. et al. Identification of genetic variants in two families with Keratoconus. BMC Med Genomics 16, 299 (2023). https://doi.org/10.1186/s12920-023-01738-x