Changes in Nuclear Gene Expression Related to Mitochondrial Function Affect Extracellular Matrix, Collagens, and Focal Adhesion in Keratoconus.

Purpose: Mitochondrial DNA (mtDNA) abnormalities were previously found to be causative in the pathogenesis of various diseases. Here, comprehensive mitochondrial and nuclear sequence and transcript analyses, along with analyses of the methylation aspects of nuclear genes related to mitochondrial function, were performed in patients with keratoconus (KTCN) to evaluate their contribution to the KTCN pathogenesis. Methods: Blood mtDNA of 42 KTCN and 51 non-KTCN individuals was Sanger sequenced and analyzed along with the previously obtained corneal RNA-sequencing data of 20 KTCN and 21 non-KTCN individuals. In addition, the expression and methylation of mtDNA genes and 1223 mitochondria-related nuclear genes were evaluated. Results: The mtDNA sequence alterations detected in blood coincided with variants identified in transcripts of the matched corneal tissues. In KTCN corneas, 97 mitochondria-related genes were deregulated, including TGFB1, P4HB, and BCL2, which are involved in the extracellular matrix (ECM) organization, collagen formation, and focal adhesion pathways. No changes in the expression of mtDNA transcripts and no differentially methylated genes among the assessed mitochondrial-nuclear gene sets were found. Conclusions: The absence of corneal-specific mtDNA variants indicates that there is no direct relationship between mitochondrial sequence variability and KTCN phenotype in the studied individuals. However, the identified KTCN-specific transcriptomic alterations of the nuclear genes directly related to the mitochondria functioning point to their possible involvement in the ECM organization, collagen formation, and focal adhesion. Translational Relevance: The identification of abnormalities within nuclear genes regulating ECM formation, collagen synthesis, and/or focal adhesion may form the basis of future treatment strategies or predict the progression of corneal changes in KTCN.

Further evaluation of differential expression of keratoconus candidate genes in human corneas.

BACKGROUND: Keratoconus (KTCN) is a progressive eye disease, characterized by changes in the shape and thickness of the cornea that results in loss of visual acuity. While numerous KTCN candidate genes have been identified, the genetic etiology of the disease remains undetermined. To further investigate and verify the contribution of particular genetic factors to KTCN, we assessed 45 candidate genes previously indicated as involved in KTCN etiology based on transcriptomic and genomic data. METHODS: The RealTime ready Custom Panel, covering 45 KTCN candidate genes and two reference transcripts, has been designed. Then, the expression profiles have been assessed using the RT-qPCR assay in six KTCN and six non-KTCN human corneas, obtained from individuals undergoing a penetrating keratoplasty procedure. RESULTS: In total, 35 genes exhibiting differential expression between KTCN and non-KTCN corneas have been identified. Among these genes were ones linked to the extracellular matrix formation, including collagen synthesis or the TGF-ß, Hippo, and Wnt signaling pathways. The most downregulated transcripts in KTCN corneas were CTGF, TGFB3, ZNF469, COL5A2, SMAD7, and SPARC, while TGFBI and SLC4A11 were the most upregulated ones. Hierarchical clustering of expression profiles demonstrated almost clear separation between KTCN and non-KTCN corneas. The gene expression levels determined using RT-qPCR showed a strong correlation with previous RNA sequencing (RNA-Seq) results. CONCLUSIONS: A strong correlation between RT-qPCR and earlier RNA-Seq data confirms the possible involvement of genes from collagen synthesis and the TGF-ß, Hippo, and Wnt signaling pathways in KTCN etiology. Our data also revealed altered expression of several genes, such as LOX, SPARC, and ZNF469, in which single nucleotide variants have been frequently identified in KTCN. These findings further highlight the heterogeneous nature of KTCN.

Accumulation of sequence variants in genes of Wnt signaling and focal adhesion pathways in human corneas further explains their involvement in keratoconus.

BACKGROUND: Keratoconus (KTCN) is a protrusion and thinning of the cornea, resulting in loss of visual acuity. The etiology of KTCN remains unclear. The purpose of this study was to assess the potential involvement of new genetic variants in KTCN etiology based on both the genomic and transcriptomic findings recognized in the same corneal tissues. METHODS: Corneal tissues derived from five unrelated Polish individuals with KTCN were examined using exome sequencing (ES), followed by enrichment analyses. For comparison purposes, the datasets comprising ES data of five randomly selected Polish individuals without ocular abnormalities and five Polish patients with high myopia were used. Expression levels of selected genes from the overrepresented pathways were obtained from the previous RNA-Seq study. RESULTS: Exome capture discovered 117 potentially relevant variants that were further narrowed by gene overrepresentation analyses. In each of five patients, the assessment of functional interactions revealed rare (MAF ≤ 0.01) DNA variants in at least one gene from Wnt signaling (VANGL1, WNT1, PPP3CC, LRP6, FZD2) and focal adhesion (BIRC2, PAK6, COL4A4, PPP1R12A, PTK6) pathways. No genes involved in pathways enriched in KTCN corneas were overrepresented in our control sample sets. CONCLUSIONS: The results of this first pilot ES profiling of human KTCN corneas emphasized that accumulation of sequence variants in several genes from Wnt signaling and/or focal adhesion pathways might cause the phenotypic effect and further points to a complex etiology of KTCN.

Collagen synthesis disruption and downregulation of core elements of TGF-ß, Hippo, and Wnt pathways in keratoconus corneas.

To understand better the factors contributing to keratoconus (KTCN), we performed comprehensive transcriptome profiling of human KTCN corneas for the first time using an RNA-Seq approach. Twenty-five KTCN and 25 non-KTCN corneas were enrolled in this study. After RNA extraction, total RNA libraries were prepared and sequenced. The discovery RNA-Seq analysis (in eight KTCN and eight non-KTCN corneas) was conducted first, after which the replication RNA-Seq experiment was performed on a second set of samples (17 KTCN and 17 non-KTCN corneas). Over 82% of the genes and almost 75% of the transcripts detected as differentially expressed in KTCN and non-KTCN corneas were confirmed in the replication study using another set of samples. We used these differentially expressed genes to generate a network of KTCN-deregulated genes. We found an extensive disruption of collagen synthesis and maturation pathways, as well as downregulation of the core elements of the TGF-ß, Hippo, and Wnt signaling pathways influencing corneal organization. This first comprehensive transcriptome profiling of human KTCN corneas points further to a complex etiology of KTCN.

KTCNlncDB-a first platform to investigate lncRNAs expressed in human keratoconus and non-keratoconus corneas.

Keratoconus (KTCN, OMIM 148300) is a degenerative eye disorder characterized by progressive stromal thinning that leads to a conical shape of the cornea, resulting in optical aberrations and even loss of visual function. The biochemical background of the disease is poorly understood, which motivated us to perform RNA-Seq experiment, aimed at better characterizing the KTCN transcriptome and identification of long non-coding RNAs (lncRNAs) that might be involved in KTCN etiology. The in silico functional studies based on predicted lncRNA:RNA base-pairings led us to recognition of a number of lncRNAs possibly regulating genes with known or plausible links to KTCN. The lncRNA sequences and data regarding their predicted functions in controlling the RNA processing and stability are available for browse, search and download in KTCNlncDB (http://rhesus.amu.edu.pl/KTCNlncDB/), the first online platform devoted to KTCN transcriptome.Database URL: http://rhesus.amu.edu.pl/KTCNlncDB/.

Evidence against ZNF469 being causative for keratoconus in Polish patients.

PURPOSE: Keratoconus (KTCN) is a degenerative disorder characterized by stromal thinning and protrusion of the cornea, resulting in severe impairment of visual function. A recent study proposed that rare heterozygous mutations in ZNF469 determine KTCN aetiology. METHODS: To investigate the contribution of ZNF469 to KTCN, we Sanger sequenced ZNF469 in 42 unrelated Polish patients with KTCN and 49 Polish individuals with high myopia (HM) and compared the results with whole-exome sequencing (WES) data performed in 268 Polish individuals without ocular abnormalities. RESULTS: The average number of ZNF469 non-synonymous variants was 16.31 and 16.0 for individuals with KTCN and HM, respectively (p = 0.3724). All identified variants were previously reported. Alternative allele frequency (AAF) was determined based on the WES results. Among missense variants, only one (rs528085780) has AAF ≤ 0.001 and was identified in one patient with sporadic KTCN. However, the resulting Arg1864Lys substitution was not predicted to be deleterious. CONCLUSION: In summary, we have not found a significant enrichment of sequence variants in ZNF469 in Polish patients with KTCN. High prevalence of ZNF469 variants identified in our KTCN group is typical for a common genetic variation observed in general population. Our findings indicate that variation in ZNF469 is not responsible for KTCN and other genetic variants are involved in the development and progression of this disease in Polish patients.

Molecular Screening of Keratoconus Susceptibility Sequence Variants in VSX1, TGFBI, DOCK9, STK24, and IPO5 Genes in Polish Patients and Novel TGFBI Variant Identification.

PURPOSE: Keratoconus (KTCN) is a degenerative disorder of the eye that results in the conical shape and thinning of the cornea and is a leading cause for corneal transplantations. A number of studies suggest that genetic factors play a role in KTCN etiology. Some candidate gene variants have recently been shown to be associated with KTCN. The purpose of our study was to verify the role of VSX1, TGFBI, DOCK9, IPO5, and STK24 sequence variants in Polish KTCN patients. METHODS: Forty-two Polish patients with sporadic KTCN and 50 control individuals were enrolled into this study. Both affected and unaffected individuals underwent detailed ophthalmic examination. The mutations screening in the candidate genes was performed by the direct sequencing method. RESULTS: Analysis of VSX1, TGFBI, DOCK9, IPO5, and STK24 genes identified numerous sequence variants. Variants c.-264_-255delGGGGTGGGGT, c.627 + 23G > A, c.809-6_809-5insT, and c.*200G > T in the VSX1 gene, and heterozygous c.1598G > A mutation (Arg533Gln) in exon 12 of TGFBI were detected for the first time in KTCN patients. Two known sequence variants of TGFBI c.1620T > C (Phe540Phe) and c.1678 + 23G > A were observed in KTCN patients and control individuals. The newly reported c.717 + 43A > G substitution in intron 7 of DOCK9 was identified in both KTCN patients and healthy individuals. CONCLUSIONS: Our investigation showed that KTCN-related sequence variants of analyzed genes were found in a very small proportion of the studied patients indicating that genes other than VSX1, TGFBI, DOCK9, IPO5, and STK24 are involved in the development and progression of KTCN in Polish patients. Our results support the hypothesis about the genetic heterogeneity of KTCN.

Lack of association between polymorphisms of the DNA base excision repair genes MUTYH and hOGG1 and keratoconus in a Polish subpopulation.

INTRODUCTION: Keratoconus (KC) is a non-inflammatory thinning of the cornea and a leading indication for corneal transplantation. Oxidative stress plays a role in the pathogenesis of this disease. The products of the hOGG1 and MUTYH genes play an important role in the repair of oxidatively modified DNA in the base excision repair pathway. We hypothesized that variability in these genes may change susceptibility to oxidative stress and predispose individuals to the development of KC. We investigated the possible association between the c.977C>G polymorphism of the hOGG1 gene (rs1052133) and the c.972G>C polymorphism of the MUTYH gene (rs3219489) and KC occurrence as well as the modulation of this association by some KC risk factors. MATERIAL AND METHODS: A total of 205 patients with KC and 220 controls were included in this study. The polymorphisms were genotyped with polymerase chain reaction (PCR) restriction fragment length polymorphism and PCR-confronting two-pair primer techniques. Differences in genotype and allele frequency distributions were evaluated using the χ(2) test, and KC risk was estimated with an unconditional multiple logistic regression with and without adjustment for co-occurrence of visual impairment, allergies, sex and family history for KC. RESULTS: We did not find any association between the genotypes and combined genotypes of the c.977C>G polymorphism of the hOGG1 gene and the c.972G>C polymorphism of the MUTYH gene and the occurrence of KC. CONCLUSIONS: Our findings suggest that the c.977C>G-hOGG1 polymorphism and the c.972G>C-MUTYH polymorphism may not be linked with KC occurrence in this Polish subpopulation.

Variation in DNA Base Excision Repair Genes in Fuchs Endothelial Corneal Dystrophy.

BACKGROUND Fuchs endothelial corneal dystrophy (FECD) is a corneal disease characterized by abnormalities in the Descemet membrane and the corneal endothelium. The etiology of this disease is poorly understood. An increased level of oxidative DNA damage reported in FECD corneas suggests a role of DNA base excision repair (BER) genes in its pathogenesis. In this work, we searched for the association between variation of the PARP-1, NEIL1, POLG, and XRCC1 genes and FECD occurrence. MATERIAL AND METHODS This study was conducted on 250 FECD patients and 353 controls using polymerase chain reaction-restriction fragment length polymorphism, high-resolution melting analysis, and the TaqMan® SNP Genotyping Assay. RESULTS We observed that the A/A genotype and the A allele of the c.1196A>G polymorphism of the XRCC1 gene were positively correlated with an increased FECD occurrence, whereas the G allele had the opposite effect. A weak association between the C/G genotype of the g.46438521G>C polymorphism of the NEIL1 gene and an increased incidence of FECD was also detected. Haplotypes of both polymorphisms of the XRCC1 were associated with FECD occurrence. No association of the c.2285T>C, c.-1370T>A and c.580C>T polymorphisms of the PARP-1, POLG and XRCC1 genes, respectively, with FECD occurrence was observed. CONCLUSIONS Our results suggest that the c.1196A>G polymorphism in the XRCC1 gene may be an independent genetic risk factor for FECD.

Polymorphism of the APEX nuclease 1 gene in keratoconus and Fuchs endothelial corneal dystrophy.

Human APEX nuclease 1 (APEX1) plays an important role in the repair of oxidative DNA lesions through base excision repair. It may influence the development of oxidative stress-related diseases. The aim of this study was to determine the relationship between the genotypes of the c.444 T>G (rs1130409) and c.-468 T>G (rs1760944) polymorphisms in the APEX1 gene and the occurrence of two oxidative stress-related eye diseases: keratoconus (KC) and Fuchs endothelial corneal dystrophy (FECD). The study involved 250 patients with KC, 209 patients with FECD, and 350 control subjects. All of the patients and control subjects underwent a detailed ophthalmic examination. The polymorphisms were genotyped by mismatch polymerase chain reaction restriction fragment length polymorphism (mismatch PCR-RFLP). We observed that the G/T and T/T genotypes of the c.-468 T>G polymorphism were respectively associated with a decreased occurrence of KC (OR 0.54, 95% CI 0.37-0.95; p = 0.030) and an increased occurrence of KC (OR 1.87, 95% CI 1.06-3.32; p = 0.032). None of these polymorphisms showed any association with FECD. Furthermore, no other association was observed, including haplotypes of the two polymorphisms. Our findings suggest that the c.-468 T>G polymorphism of the APEX1 gene may play a role in the pathogenesis of KC.

Polymorphism of the flap endonuclease 1 gene in keratoconus and Fuchs endothelial corneal dystrophy.

Oxidative stress is implicated in the pathogenesis of many diseases, including serious ocular diseases, keratoconus (KC) and Fuchs endothelial corneal dystrophy (FECD). Flap endonuclease 1 (FEN1) plays an important role in the repair of oxidative DNA damage in the base excision repair pathway. We determined the association between two single nucleotide polymorphisms (SNPs), c.-441G>A (rs174538) and g.61564299G>T (rs4246215), in the FEN1 gene and the occurrence of KC and FECD. This study involved 279 patients with KC, 225 patients with FECD and 322 control individuals. Polymerase chain reaction (PCR) and length polymorphism restriction fragment analysis (RFLP) were applied. The T/T genotype of the g.61564299G>T polymorphism was associated with an increased occurrence of KC and FECD. There was no association between the c.-441G>A polymorphism and either disease. However, the GG haplotype of both polymorphisms was observed more frequently and the GT haplotype less frequently in the KC group than the control. The AG haplotype was associated with increased FECD occurrence. Our findings suggest that the g.61564299G>T and c.-441G>A polymorphisms in the FEN1 gene may modulate the risk of keratoconus and Fuchs endothelial corneal dystrophy.

Polymorphism of the transferrin gene in eye diseases: keratoconus and Fuchs endothelial corneal dystrophy.

Oxidative stress may play a role in the pathogenesis of keratoconus (KC) and Fuchs endothelial corneal dystrophy (FECD). Iron may promote the stress by the Fenton reaction, so its homeostasis should be strictly controlled. Transferrin is essential for iron homeostasis because it transports iron from plasma into cells. The malfunction of transferrin, which may be caused by variation in its gene (TF) variation, may contribute to oxidative stress and change KC and FECD risk. To verify this hypothesis we investigated the association between three polymorphisms of the TF gene, g.3296G>A (rs8177178), g.3481A>G (rs8177179), and c.-2G>A (rs1130459), and KC and FECD occurrence. Genotyping was performed in blood lymphocytes in 216 patients with KC, 130 patients with FECD and 228 controls by PCR-RFLP. We studied also the influence of other risk factors. The A/A genotype and the A allele of the g.3296G>A polymorphism were associated with KC occurrence, while the G allele was negatively correlated with it. We observed a decrease in KC occurrence associated with the A/G genotype of the g.3481A>G polymorphism. We did not find any association between the c.-2G>A polymorphism and KC. No association was found between all three polymorphisms and FECD occurrence.

Novel mutation and three other sequence variants segregating with phenotype at keratoconus 13q32 susceptibility locus.

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.2262A>C (p.Gln754His) and c.720+43A>G in DOCK9; c.2377-132A>C in IPO5 and c.1053+29G>C in STK24. PolyPhen analyses predicted that c.2262A>C (Gln754His) is possibly damaging for the protein function and structure. Our results suggest that c.2262A>C (p.Gln754His) mutation in DOCK9 may contribute to the KTCN phenotype in the large KTCN-014 family.