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

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.

Differentially expressed microRNAs targeting genes in key pathways in keratoconus

Introduction: Keratoconus (KTCN) is a corneal ectasia, characterized by a progressive thinning and protrusion of the cornea, with a complex etiology involving genetic, behavioral, lifestyle, and environmental factors. Previous studies indicated that microRNAs (miRNAs) could be involved in KTCN pathogenesis. This in silico study aimed to identify precursor microRNAs (pre-miRNAs) differentially expressed in KTCN corneas and to characterize mature miRNAs and their target genes. Materials and methods: Expression levels of pre-miRNAs were retrieved from our previously obtained RNA sequencing data of 25 KTCN and 25 non-KTCN human corneas (PMID:28145428, PMID:30994860). Differential expression with FDR ≤0.01 and ≥1.5-fold changes were considered significant. Lists of target genes (target score ≥90) of mature miRNAs were obtained from miRDB. Revealed up-/downregulated miRNAs and their target genes were assessed in databases and literature. Enrichment analyses were completed applying the DAVID database. Results: From a total of 47 pre-miRNAs, six were remarkably upregulated (MIR184, MIR548I1, MIR200A, MIR6728, MIR429, MIR1299) and four downregulated (MIR6081, MIR27B, MIR23B, MIR23A) in KTCN corneas. Out of the 1,409 target genes, 220 genes with decreased and 57 genes with increased expression levels in KTCN samples vs non-KTCN samples were found. The extracellular matrix (ECM) organization, response to mechanical stimulus, regulation of cell shape, and signal transduction processes/pathways were identified as distinctive in enrichment analyses. Also, processes associated with the regulation of transcription and DNA binding were listed. Conclusion: Indicated miRNAs and their target genes might be involved in KTCN pathogenesis via disruption of crucial molecular processes, including ECM organization and signal transduction.

Genetic Correlations Among Corneal Biophysical Parameters and Anthropometric Traits

Purpose

The genetic architecture of corneal dysfunction remains poorly understood. Epidemiological and clinical evidence suggests a relationship between corneal structural features and anthropometric measures. We used global and local genetic similarity analysis to identify genomic features that may underlie structural corneal dysfunction.

Methods

We assembled genome-wide association study summary statistics for corneal features (central corneal thickness, corneal hysteresis [CH], corneal resistance factor [CRF], and the 3 mm index of keratometry) and anthropometric traits (body mass index, weight, and height) in Europeans. We calculated global genetic correlations (rg) between traits using linkage disequilibrium (LD) score regression and local genetic covariance using ρ-HESS, which partitions the genome and performs regression with LD regions. Finally, we identified genes located within regions of significant genetic covariance and analyzed patterns of tissue expression and pathway enrichment.

Results

Global LD score regression revealed significant negative correlations between height and both CH (rg = -0.12; P = 2.0 × 10-7) and CRF (rg = -0.11; P = 6.9 × 10-7). Local analysis revealed 68 genomic regions exhibiting significant local genetic covariance between CRF and height, containing 2874 unique genes. Pathway analysis of genes in regions with significant local rg revealed enrichment among signaling pathways with known keratoconus associations, including cadherin and Wnt signaling, as well as enrichment of genes modulated by copper and zinc ions.

Conclusions

Corneal biophysical parameters and height share a common genomic architecture, which may facilitate identification of disease-associated genes and therapies for corneal ectasias.

Translational Relevance

Local genetic covariance analysis enables the identification of associated genes and therapeutic targets for corneal ectatic 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.

The Impaired Wound Healing Process Is a Major Factor in Remodeling of the Corneal Epithelium in Adult and Adolescent Patients With Keratoconus

Purpose

Keratoconus (KTCN) is the most common corneal ectasia, characterized by pathological cone formation. Here, to provide an insight into the remodeling of the corneal epithelium (CE) during the course of the disease, we evaluated topographic regions of the CE of adult and adolescent patients with KTCN.

Methods

The CE samples from 17 adult and 6 adolescent patients with KTCN, and 5 control CE samples were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated toward RNA sequencing and MALDI-TOF/TOF Tandem Mass Spectrometry. Data from transcriptomic and proteomic investigations were consolidated with the morphological and clinical findings.

Results

The critical elements of the wound healing process, epithelial-mesenchymal transition, cell-cell communications, and cell-extracellular matrix interactions were altered in the particular corneal topographic regions. Abnormalities in pathways of neutrophils degranulation, extracellular matrix processing, apical junctions, IL, and IFN signaling were revealed to cooperatively disorganize the epithelial healing. Deregulation of the epithelial healing, G2M checkpoints, apoptosis, and DNA repair pathways in the middle CE topographic region in KTCN explains the presence of morphological changes in the corresponding doughnut pattern (a thin cone center surrounded by a thickened annulus). Despite similar morphological characteristics of CE samples in adolescents and adults with KTCN, their transcriptomic features were different. Values of the posterior corneal elevation differentiated adults with KTCN from adolescents with KTCN and correlated with the expression of TCHP, SPATA13, CNOT3, WNK1, TGFB2, and KRT12 genes.

Conclusions

Identified molecular, morphological, and clinical features indicate the effect of impaired wound healing on corneal remodeling in KTCN CE.

Non-allergic eye rubbing is a major behavioral risk factor for keratoconus

Since the environmental, behavioral, and socioeconomic factors in the etiology of keratoconus (KTCN) remain poorly understood, we characterized them as features influencing KTCN phenotype, and especially affecting the corneal epithelium (CE). In this case-control study, 118 KTCN patients and 73 controls were clinically examined and the Questionnaire covering the aforementioned aspects was completed and then statistically elaborated. Selected KTCN-specific findings were correlated with the outcomes of the RNA-seq assessment of the CE samples. Male sex, eye rubbing, time of using a computer after work, and dust in the working environment, were the substantial KTCN risk factors identified in multivariate analysis, with ORs of 8.66, 7.36, 2.35, and 5.25, respectively. Analyses for genes whose expression in the CE was correlated with the eye rubbing manner showed the enrichment in apoptosis (TP53, BCL2L1), chaperon-related (TLN1, CTDSP2, SRPRA), unfolded protein response (NFYA, TLN1, CTDSP2, SRPRA), cell adhesion (TGFBI, PTPN1, PDPK1), and cellular stress (TFDP1, SRPRA, CAPZB) pathways. Genes whose expression was extrapolated to the allergy status didn't contribute to IgE-related or other inflammatory pathways. Presented findings support the hypothesis of chronic mechanical corneal trauma in KTCN. Eye-rubbing causes CE damage and triggers cellular stress which through its influence on cell apoptosis, migration, and adhesion affects the KTCN phenotype.

Comprehensive Transcriptome Analysis of Patients With Keratoconus Highlights the Regulation of Immune Responses and Inflammatory Processes

Keratoconus (KTCN), characterized by the steeper curvature of the cornea and the thinner central corneal thickness, was a degenerative corneal ectasia with ambiguous etiology and mechanism. We aim to reveal underlying pathological mechanisms of KTCN by multi-level transcriptomic, integrative omics analyses. We performed RNA-sequencing on corneal epithelial samples from seven patients and seven control donors, as well as peripheral matched blood samples from three of the patients and three control donors. After RNA extraction, library construction, and sequencing, differentially expressed genes and splicing events were identified, followed by over-representation analysis. In total, 547 differential expressed genes were identified in KTCN corneal epithelial samples, which were mainly enriched in immune responses and inflammatory processes. WGCNA module analysis, the transcriptomic analysis of peripheral blood samples, multiple omics data, and the meta-analysis of GEO samples confirmed the involvement of immune and inflammatory factors. Besides, 190 and 1,163 aberrant splicing events were identified by rMATS combined with CASH methods in corneal epithelial and blood samples with KCTN. In conclusion, this comprehensive transcriptome analysis of KTCN patients based on patients’ tissue and blood samples uncovered a significant association between immune-inflammatory genes and pathways with KCTN, highlighting the contribution of the perturbed immune signaling to the pathogenesis of KCTN. Our study suggested the importance of measures to control inflammation in the treatment of KCTN.

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.

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.

Time- and Stimulus-Dependent Characteristics of Innate Immune Cells in Organ-Cultured Human Corneal Tissue

PURPOSE

The pattern of immune cells infiltrating the corneal stroma has been extensively studied in mice, but data on human tissue have been far less elaborate. To further characterize the number and differentiation state of resident immune cells in organ-cultured human corneal tissue, we employed a comprehensive bioinformatic deconvolution (xCell) of bulk RNA-sequencing (RNA-seq) data, immunohistochemistry (IHC), and flow cytometry (FC).

METHODS

A transcriptome-based analysis of immune cell types in human corneal samples was performed. The results were validated by IHC, focusing on the identification of pro-inflammatory (M1) and regulatory (M2) macrophages. A protocol was established to identify these 2 different macrophage populations in human corneal tissue by means of FC. Subsequently, corneal samples in organ culture were differentially stimulated by IL-10, IL-4 & IL-13, or LPS and macrophage populations were evaluated regarding their response to these stimuli. Furthermore, cell survival was analyzed in correlation with time in organ culture.

RESULTS

xCell-based mathematical deconvolution of bulk RNA-seq data revealed the presence of CD8 T cells, Th17 cells, dendritic cells, and macrophages as the predominant immune cell types in organ-cultured human corneal tissue. Furthermore, RNA-seq allowed the detection of different macrophage marker genes in corneal samples, including PTPRC (CD45), ITGAM (CD11b), CD14, and CD74. Our RNA-seq data showed no evidence of a relevant presence of monocytes in human corneal tissue. The presence of different macrophage subtypes was confirmed by IHC. The disintegration and subsequent FC analysis of human corneal samples showed the presence of both M1 (HLA-DR+, CD282+, CD86+, and CD284+) and M2 (CD163+ and CD206+) macrophage subtypes. Furthermore, we found that the total number of macrophages in corneal samples decreased more than the total cell count with increasing tissue culture time. Treatment with IL-10 led to higher total cell counts per cornea and to an increased expression of the M2 marker CD163 (p < 0.05) while expression levels of various M1 macrophage markers were not significantly reduced by interleukin treatment.

CONCLUSIONS

Regarding different macrophage populations, untreated human corneas showed more M1 than M2 macrophages. With increasing organ culture time, these macrophages decreased. In terms of cell dynamics, adding interleukins to the organ culture medium influenced the phenotype of macrophages within the cornea as detected by FC. Modifying the immunomodulatory properties of human grafts appears a promising approach to further reduce the risk of graft rejection in patients. In this context, treatment with interleukins was more effective in upregulating M2 macrophages than in suppressing M1 macrophages in corneal tissue.

Pathogenic alleles in microtubule, secretory granule and extracellular matrix-related genes in familial keratoconus

Keratoconus is a common corneal defect with a complex genetic basis. By whole exome sequencing of affected members from 11 multiplex families of European ancestry, we identified 23 rare, heterozygous, potentially pathogenic variants in 8 genes. These include nonsynonymous single amino acid substitutions in HSPG2, EML6 and CENPF in two families each, and in NBEAL2, LRP1B, PIK3CG and MRGPRD in three families each; ITGAX had nonsynonymous single amino acid substitutions in two families and an indel with a base substitution producing a nonsense allele in the third family. Only HSPG2, EML6 and CENPF have been associated with ocular phenotypes previously. With the exception of MRGPRD and ITGAX, we detected the transcript and encoded protein of the remaining genes in the cornea and corneal cell cultures. Cultured stromal cells showed cytoplasmic punctate staining of NBEAL2, staining of the fibrillar cytoskeletal network by EML6, while CENPF localized to the basal body of primary cilia. We inhibited the expression of HSPG2, EML6, NBEAL2 and CENPF in stromal cell cultures and assayed for the expression of COL1A1 as a readout of corneal matrix production. An upregulation in COL1A1 after siRNA inhibition indicated their functional link to stromal cell biology. For ITGAX, encoding a leukocyte integrin, we assayed its level in the sera of 3 affected families compared with 10 unrelated controls to detect an increase in all affecteds. Our study identified genes that regulate the cytoskeleton, protein trafficking and secretion, barrier tissue function and response to injury and inflammation, as being relevant to keratoconus.

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.

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.