Three novel variants identified within ECM-related genes in Chinese Han keratoconus patients

Identification of hub genes and molecular pathways in keratoconus by integrating bioinformatics and literature mining at the RNA level

OBJECTIVE

Keratoconus (KC) is a condition characterized by progressive corneal steepening and thinning. However, its pathophysiological mechanism remains vague. We mainly performed literature mining to extract bioinformatic and related data on KC at the RNA level. The objective of this study was to explore the potential pathological mechanisms of KC by identifying hub genes and key molecular pathways at the RNA level.

METHODS

We performed an exhaustive search of the PubMed database and identified studies that pertained to gene transcripts derived from diverse corneal layers in patients with KC. The identified differentially expressed genes were intersected, and overlapping genes were extracted for further analyses. Significantly enriched genes were screened using "Gene Ontology" (GO) and "Kyoto Encyclopedia of Genes and Genomes" (KEGG) analysis with the "Database for Annotation, Visualization, and Integrated Discovery" (DAVID) database. A protein-protein interaction (PPI) network was constructed for the significantly enriched genes using the STRING database. The PPI network was visualized using the Cytoscape software, and hub genes were screened via betweenness centrality values. Pathways that play a critical role in the pathophysiology of KC were discovered using the GO and KEGG analyses of the hub genes.

RESULTS

68 overlapping genes were obtained. Fifty genes were significantly enriched in 67 biological processes, and 16 genes were identified in 7 KEGG pathways. Moreover, 14 nodes and 32 edges were identified via the PPI network constructed using the STRING database. Multiple analyses identified 4 hub genes, 12 enriched biological processes, and 6 KEGG pathways. GO enrichment analysis showed that the hub genes are mainly involved in the positive regulation of apoptotic process, and KEGG analysis showed that the hub genes are primarily associated with the interleukin-17 (IL-17) and tumor necrosis factor (TNF) pathways. Overall, the matrix metalloproteinase 9, IL-6, estrogen receptor 1, and prostaglandin-endoperoxide synthase 2 were the potential important genes associated with KC.

CONCLUSION

Four genes, matrix metalloproteinase 9, IL-6, estrogen receptor 1, and prostaglandin endoperoxide synthase 2, as well as IL-17 and TNF pathways, are critical in the development of KC. Inflammation and apoptosis may contribute to the pathogenesis of KC.

Whole-exome sequencing screening for candidate genes and variants associated with primary sporadic keratoconus in Chinese patients

The pathogenesis of keratoconus (KC) is complex, and genetic factors play an important role. The purpose of this study was to screen and analyse candidate genes and variants in Chinese patients with primary sporadic KC. Whole-exome sequencing (WES) was performed to identify candidate genes and variants in 105 unrelated Chinese patients with primary sporadic KC. Through a series of screening processes, 54 candidate variants in 26 KC candidate genes were identified in 53 KC patients (53/105, 50.5%). These 54 candidate variants included 10 previously identified variants in 9 KC candidate genes and 44 novel variants in 20 KC candidate genes. The previously identified variants occurred in 25.7% (27/105) of patients. Of these, 4 variants (COL6A5, c.5014T > G; CAST, c.1814G > A; ZNF469, c.946G > A; and MPDZ, c.3836A > G) were identified for the first time in Chinese KC patients. The novel variants occurred in 33.3% (35/105) of patients. Of the 26 screened KC candidate genes, 11 KC candidate genes (CAT, COL12A1, FLG, HKDC1, HSPG2, PLOD1, ITGA2, TFAP2B, USH2A, WNT10A, and COL6A5) were found to be potentially pathogenic in Chinese KC patients for the first time. Gene Ontology (GO) biological process (BP) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed on the 26 KC candidate genes using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). The results showed that the KC candidate genes were significantly enriched in biological processes such as collagen fibril organization and extracellular matrix (ECM) organization and in ECM-receptor interaction and protein digestion and absorption pathways. The results further expand the spectrum of KC candidate variants and provide a basis for further KC gene studies.

Expression of Tissue Inhibitors of Metalloproteinases (TIMP-1, TIMP-2, TIMP-3, TIMP-4) in Blood Serum of Patients with Keratoconus

BACKGROUND

The etiology of keratoconus is unclear. Current evidence suggests that inflammatory and systemic mechanisms might play a role in its pathophysiology. The proper interaction of proteolytic enzymes-matrix metalloproteinases-and their specific tissue inhibitors (TIMPs) within the cornea is essential in maintaining its structure, transparency and healing processes. The aim of the study was to determine the concentration of the TIMPs TIMP-1, TIMP-2, TIMP-3, and TIMP-4 in the blood serum samples of patients with keratoconus compared to the control group.

METHODS

The study encompassed 132 patients, of which 83 people constituted the study group and 49 the control group. The concentration of selected TIMPs was determined using the Human Magnetic Luminex® Performance Assay method.

RESULTS

In the study group, the concentrations of TIMP-1 and TIMP-3 were statistically significantly reduced, and TIMP-2 and TIMP-4 increased compared to the control group. The analysis of individual TIMPs in terms of their usefulness as potential predictors of keratoconus showed high results of diagnostic sensitivity and specificity for all TIMPs, in particular for TIMP-1 and TIMP-2.

CONCLUSION

The above results may indicate systemic disturbances in the TIMPs regulation among keratoconus patients. High diagnostic sensitivity and specificity of all TIMPs, in particular TIMP-1 and TIMP-2, may confirm their participation in the etiopathogenesis of this disease.

Sequence variants contributing to dysregulated inflammatory responses across keratoconic cone surface in adolescent patients with keratoconus

Background

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Methods

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

Results

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

Conclusion

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

Novel Mutations Identified in the Chinese Han Population with Keratoconus by Next-Generation Sequencing

Aim. To identify novel mutations in keratoconus (KC) susceptibility genes in the Chinese Han population. Methods. A total of fifty-two patients with primary KC were recruited. Blood samples were collected, and genomic DNA was isolated from peripheral blood leukocytes. The entire coding region, intron-exon junctions, and promoter regions of sixteen known KC susceptibility genes were screened with next-generation sequencing technology. All identified variants were further confirmed using the Sanger sequencing technology. The Sorting Intolerant from Tolerant (SIFT), MutationTaster, and PolyPhen 2 programs were used to predict the effect of amino acid substitution on protein. Results. After removing twelve known SNPs (single nucleotide polymorphisms) and three variants predicted to be harmless, nine novel mutations were identified in eight of the fifty-two patients, including c.455C > T:p.P152L in FNDC3B; c.3636_3637del:p.R1212fs in COL4A4; c.5015G > T:p.R1672L, c.3798dupA:p.P1267fs, and c.28G > A:p.A10T in MPDZ; c.1940C > T:p.P647L in DOCK9; c.127_128insGGC:p.Q43delinsRQ in POLG; c.3019G > A:p.V1007I in IPO5; and c.624 + 7− > A in TGFBI. All nine mutations in the patients with KC were heterozygote. Conclusion. This study enlarged the gene profile of KC and should be further confirmed by well-powered, genome-wide association studies (GWAS) of Han Chinese patients.

Systematically Displaying the Pathogenesis of Keratoconus via Multi-Level Related Gene Enrichment-Based Review

Keratoconus (KC) is an etiologically heterogeneous corneal ectatic disorder. To systematically display the pathogenesis of keratoconus (KC), this study reviewed all the reported genes involved in KC, and performed an enrichment analysis of genes identified at the genome, transcription, and protein levels respectively. Combined analysis of multi-level results revealed their shared genes, gene ontology (GO), and pathway terms, to explore the possible pathogenesis of KC. After an initial search, 80 candidate genes, 2,933 transcriptional differential genes, and 947 differential proteins were collected. The candidate genes were significantly enriched in extracellular matrix (ECM) related terms, Wnt signaling pathway and cytokine activities. The enriched GO/pathway terms of transcription and protein levels highlight the importance of ECM, cell adhesion, and inflammatory once again. Combined analysis of multi-levels identified 13 genes, 43 GOs, and 12 pathways. The pathogenic relationships among these overlapping factors maybe as follows. The gene mutations/variants caused insufficient protein dosage or abnormal function, together with environmental stimulation, leading to the related functions and pathways changes in the corneal cells. These included response to the glucocorticoid and reactive oxygen species; regulation of various signaling (P13K-AKT, MAPK and NF-kappaB), apoptosis and aging; upregulation of cytokines and collagen-related enzymes; and downregulation of collagen and other ECM-related proteins. These undoubtedly lead to a reduction of extracellular components and induction of cell apoptosis, resulting in the loosening and thinning of corneal tissue structure. This study, in addition to providing information about the genes involved, also provides an integrated insight into the gene-based etiology and pathogenesis of KC.

Cascade of interactions between candidate genes reveals convergent mechanisms in keratoconus disease pathogenesis

Keratoconus is a progressive thinning, steepening and distortion of the cornea which can lead to loss of vision if left untreated. Keratoconus has a complex multifactorial etiology, with genetic and environmental components contributing to the disease pathophysiology. Studies have observed high concordance between monozygotic twins, discordance between dizygotic twins, and high familial segregation indicating the presence of a very strong genetic component in the pathogenesis of keratoconus. The use of genome-wide linkage studies on families and twins, genome-wide association studies (GWAS) on case-controls, next-generation sequencing (NGS)-based genomic screens on both familial and non-familial cohorts have led to the identification of keratoconus candidate genes with much greater success and increased resproducibility of genetic findings. This review focuses on candidate genes identified till date and attempts to understand their role in biological processes underlying keratoconus pathogenesis. In addition, using these genes I propose molecular pathways that could contribute to keratoconus pathogenesis. The pathways identified the presence of direct cross-talk between known candidate genes of keratoconus and remarkably, 28 known candidate genes have a direct relationship among themselves that involves direct protein-protein binding, regulatory activities such as activation and inhibition, chaperone, transcriptional activation/co-activation, and enzyme catalysis. This review attempts to describe these relationships and cross-talks in the context of keratoconus pathogenesis.

Update on the genetics of keratoconus

In the past few years we have seen a great acceleration of discoveries in the field of keratoconus including new treatments, diagnostic tools, genomic and molecular determinants of disease risk. Recent genome-wide association studies (GWAS) of keratoconus cases and population wide studies of variation in central corneal thickness and in corneal biomechanical properties confirmed already identified genes and found many new susceptibility variants and biological pathways. Recent findings in genetic determinants of familial keratoconus revealed functionally important variants and established first mouse model of keratoconus. Latest transcriptomic and expression studies started assessing novel non-coding RNA targets in addition to identifying tissue specific effects of coding genes. First genomic insights into better prediction of treatment outcomes are bringing the advent of genomic medicine into keratoconus clinical practice.

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.