Network Integration Analysis and Immune Infiltration Analysis Reveal Potential Biomarkers for Primary Open-Angle Glaucoma

Core genes and immune dysregulation in primary open-angle glaucoma: A molecular insight

BackgroundPrimary open-angle glaucoma (POAG) is a chronic, progressive and irreversible eye disease. Currently, there is no effective way to prevent optic nerve damage.ObjectiveThis study explored POAG gene markers to identify high-risk groups at an early stage and to find new effective therapeutic targets.MethodsThe mRNA and clinical information of POAG patients and normal samples were downloaded from the Gene Expression Omnibus (GEO) database. Through Weighted correlation network analysis (WGCNA) and generalized linear models (GLM), random forests (RF), support vector machines (SVM), and extreme gradient boosting (xGB) models, key risk genes were identified and an early diagnosis model was established. Functional enrichment analysis and CIBERSORT algorithm were used to further reveal the changes in the POAG immune environment and find emerging therapeutic targets.ResultsHERPUD1, IQCK, MRPL40, SRSF7 and TMEM243 were identified as risk genes, and the prediction model and nomogram constructed based on them had good early prediction efficiency. At the mechanistic level, the heterogeneity of T cell subsets seems to be a key factor affecting the progression of POAG and has potential therapeutic value.Conclusions: HERPUD1, IQCK, MRPL40, SRSF7, and TMEM243 are of great significance for the early prediction and disease progression of POAG and have the potential value of becoming therapeutic targets.

Exploring diagnostic m6A regulators in primary open-angle glaucoma: insight from gene signature and possible mechanisms by which key genes function

PURPOSE

The purpose of this study was to interrogate the potential role of N6-methyladenosine (m6A) regulators in the process of trabecular meshwork (TM) tissue damage in patients with primary open-angle glaucoma (POAG).

METHODS

Firstly, the expression profile of m6A regulators in TM tissues of POAG patients was comprehensively analyzed by bioinformatics analysis; Plasmid transfection and siRNA gene interference were used to enhance or weaken the expression levels of YTHDC2 in human trabecular meshwork cells (HTMCs); Cell migration ability was detected by transwell chamber assay; Immunofluorescence staining assay was used to evaluate the expression of extracellular matrix (ECM) related proteins.

RESULTS

Through the analysis of GSE27276 database, 5 m6A regulators with different expression in POAG were screened out. The results of random forest model showed that these 5 m6A regulators exhibited diagnostic potential and were characteristic genes of POAG. All POAG samples could be effectively divided into two groups based on the expression levels of these 5 hub m6A regulators. Immune cell infiltration analysis indicated that the levels of activated CD8+ T cells and regulatory T cells were different in the two subtypes. HTMC oxidative stress cell model and TGF-β2 stimulation cell model were further constructed to verify the expression of the aforementioned hub m6A regulators, and it was found that YTHDC2 mRNA showed the same expression trend in both models. The silencing of YTHDC2 enhanced the migration ability of HTMCs and increased the synthesis ability of ECM. However, when YTHDC2ΔYTH, which lacks the YTH domain, is overexpressed in HTMCs, there is no significant change in the ECM synthesis ability.

CONCLUSIONS

The differentially expressed m6A regulators in TM tissues may serve as potential diagnostic biomarkers for POAG. And, in HTMCs, the expression level of YTHDC2 mRNA was changed under oxidative stress or TGF-β2 intervention, and then exerted its regulation on cell migration and ECM synthesis capability through m6A modification, which may be an important part of the disease process of POAG.

Identification of a circulating three-miRNA panel for the diagnosis of primary open angle glaucoma

PURPOSE

Conventional diagnosis of primary open angle glaucoma (POAG) needs a combination of ophthalmic examinations. An efficient assay is urgently needed for a timely POAG diagnosis. We aim to explore differential expressions of circulating microRNAs (miRNA) and provide novel miRNA biomarkers for POAG diagnosis.

METHODS

A total of 180 POAG patients and 210 age-related cataract (ARC) patients were enrolled. We collected aqueous humor (AH) and plasma samples from the recruited patients. The expressions of candidate miRNAs were measured using quantitative real time polymerase chain reaction. The diagnostic ability of candidate miRNAs was analyzed by receiver operating characteristic curve.

RESULTS

The expressions of miR-21-5p and miR-29b-3p were downregulated significantly in AH and plasma of POAG and miR-24-3p expression was significantly increased in AH and plasma of POAG, comparing with those of ARC. A three-miRNA panel was constructed by a binary logistic regression. And the panel could differentiate between POAG and ARC with an area under the curve of 0.8867 (sensitivity = 78.0%, specificity = 83.3%) in aqueous humor and 0.7547 (sensitivity = 73.8%, specificity = 81.2%) in plasma. Next, we verified the three-miRNA panel working as a potential diagnostic biomarker stable and reliable. At last, we identified related function and regulation pathways in vitro.

CONCLUSIONS

In conclusion, we built and identified a circulating three-miRNA panel as a potential diagnostic biomarker for POAG. It may be developed into an efficient assay and help improve the POAG diagnosis in the future.

Exploring Epigenetic Modifications as Potential Biomarkers and Therapeutic Targets in Glaucoma

Glaucoma, a complex and multifactorial neurodegenerative disorder, is a leading cause of irreversible blindness worldwide. Despite significant advancements in our understanding of its pathogenesis and management, early diagnosis and effective treatment of glaucoma remain major clinical challenges. Epigenetic modifications, encompassing deoxyribonucleic acid (DNA) methylation, histone modifications, and non-coding RNAs, have emerged as critical regulators of gene expression and cellular processes. The aim of this comprehensive review focuses on the emerging field of epigenetics and its role in understanding the complex genetic and molecular mechanisms underlying glaucoma. The review will provide an overview of the pathophysiology of glaucoma, emphasizing the intricacies of intraocular pressure regulation, retinal ganglion cell dysfunction, and optic nerve damage. It explores how epigenetic modifications, such as DNA methylation and histone modifications, can influence gene expression, and how these mechanisms are implicated in glaucomatous neurodegeneration and contribute to glaucoma pathogenesis. The manuscript discusses evidence from both animal models and human studies, providing insights into the epigenetic alterations associated with glaucoma onset and progression. Additionally, it discusses the potential of using epigenetic modifications as diagnostic biomarkers and therapeutic targets for more personalized and targeted glaucoma treatment.

Identification and validation of key biomarkers and potential therapeutic targets for primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is a prevalent cause of blindness worldwide, resulting in degeneration of retinal ganglion cells and permanent damage to the optic nerve. However, the underlying pathogenetic mechanisms of POAG are currently indistinct, and there has been no effective nonsurgical treatment regimen. The objective of this study is to identify novel biomarkers and potential therapeutic targets for POAG. The mRNA expression microarray datasets GSE27276 and GSE138125, as well as the single-cell high-throughput RNA sequencing (scRNA-seq) dataset GSE148371 were utilized to screen POAG-related differentially expressed genes (DEGs). Functional enrichment analyses, protein-protein interaction (PPI) analysis, and weighted gene co-expression network analysis (WGCNA) of the DEGs were performed. Subsequently, the hub genes were validated at a single-cell level, where trabecular cells were annotated, and the mRNA expression levels of target genes in different cell clusters were analyzed. Immunofluorescence and quantitative real-time PCR (qPCR) were performed for further validation. DEGs analysis identified 43 downregulated and 32 upregulated genes in POAG, which were mainly enriched in immune-related pathways, oxidative stress, and endoplasmic reticulum (ER) stress. PPI networks showed that FN1 and DUSP1 were the central hub nodes, while GPX3 and VAV3 were screened out as hub genes through WGCNA and subsequently validated by qPCR. Finally, FN1, GPX3, and VAV3 were determined to be pivotal core genes via single-cell validation. The relevant biomarkers involved in the pathogenesis of POAG, may serve as potential therapeutic targets. Further studies are necessary to unveil the mechanisms underlying the expression variations of these genes in POAG.

Construction of miRNA-mRNA regulatory network indicates potential biomarkers for primary open-angle glaucoma

BACKGROUND

Trabecular meshwork (TM) dysfunction-induced elevation of intraocular pressure has been identified as the main risk factor of irreversible optic nerve injury in Primary open‑angle glaucoma (POAG). Increasing evidences suggest that microRNA (miRNA) plays a vital role in the pathogenesis of POAG. This study aims to construct a miRNA-mRNA regulatory network and identify biomarkers for POAG.

METHODS

miRNAs and mRNAs expression profiling of TM samples from controls and POAG patients were assessed through microarray analysis. Target genes of differentially expressed miRNAs (DEmiRNAs) were predicted by miEAA and miRNet. Then GO and KEGG pathway analysis of differentially expressed mRNAs (DEmRNAs) were performed. PPI of top 30 hub genes was identified and miRNA-mRNA network was established by STRING database and Cytoscape software. GSE27276 and GSE105269 datasets were used to verify the expression of hub genes and to predict potential biomarkers in TM and aqueous humor (AH) for POAG, respectively. Finally, GSEA analysis was conducted to estimate the main signaling pathway of POAG pathogenesis.

RESULTS

A total of 29 up-regulated and 7 down-regulated miRNAs, 923 up-regulated and 887 down-regulated mRNAs were identified in TM of POAG compared with controls. Target genes and DEmRNAs were mainly enriched in nitric oxide biosynthetic process, vasopressin-regulated water reabsorption, and so on. Through miRNA-mRNA network construction, top 30 hub genes were regulated by 24 DEmiRNAs. 8 genes were aberrantly expressed in dataset GSE27276. 3 genes (CREB1, CAPZA2, SLC2A3) and 2 miRNAs (miR-106b-5p, miR-15a-5p) were identified as potential biomarkers for POAG in TM and AH, respectively. GSEA analysis revealed that these 3 genes modulated POAG through different pathways.

CONCLUSION

In this study, construction of miRNA-mRNA network and identification of biomarkers provide a novel insight into the pathogenesis, early diagnosis and treatment for POAG.

Gender-specific lncRNA-miRNA-mRNA regulatory network to reveal potential genes for primary open-angle glaucoma

BACKGROUND

Investigation of biomarkers may facilitate understanding the mechanisms of primary open-angle glaucoma (POAG) and developing therapeutic targets. This study aimed to identify potential genes based on competing endogenous RNA (ceRNA) network for POAG.

METHODS

Based on long noncoding RNAs (lncRNAs), microRNAs (miRNAs) and messenger RNAs (mRNAs) from the Gene Expression Omnibus (GEO) database, we identified differential expressed lncRNAs (DELs), differential expressed miRNAs (DEMis) and differential expressed mRNAs (DEMs) and then constructed a ceRNA network. Through weighted gene co-expression network analysis (WGCNA), we identified gender-specific genes for gender-associated ceRNA network construction, followed by the protein-protein interaction (PPI) network and functional enrichment analysis to screen hub genes and reveal their functions. The expression levels of hub genes were measured in steroid-induced ocular hypertension (SIOH) mice.

RESULTS

A total of 175 DELs, 727 DEMs and 45 DEMis were screened between control and POAG samples. Seven modules were identified through WGCNA and one module was associated with gender of POAG patients. We discovered 41 gender-specific genes for gender-associated ceRNA construction and then identified 8 genes (NAV3, C1QB, RXRB, P2RY4, ADAM15, VAV3, ZNF207 and TOP1), which were enriched in cell cycle-related pathways and immune-related pathways. C1QB, RXRB, Top1 and ZNF207 were highly interacted with other proteins. The expression levels of NAV3 and C1QB were downregulated in SIOH, while the levels of RXRB, P2RY4, ADAM15, VAV3, ZNF207 and TOP1 were upregulated in SIOH.

CONCLUSION

This study identifies hub genes associated with the pathogenesis of gender-specific POAG and provides potential biomarkers for POAG.

The Application of MicroRNAs in Glaucoma Research: A Bibliometric and Visualized Analysis

Glaucoma is similar to a neurodegenerative disorder and leads to global irreversible loss of vision. Despite extensive research, the pathophysiological mechanisms of glaucoma remain unclear, and no complete cure has yet been identified for glaucoma. Recent studies have shown that microRNAs can serve as diagnostic biomarkers or therapeutic targets for glaucoma; however, there are few bibliometric studies that focus on using microRNAs in glaucoma research. Here, we have adopted a bibliometric analysis in the field of microRNAs in glaucoma research to manifest the current tendencies and research hotspots and to present a visual map of the past and emerging tendencies in this field. In this study, we retrieved publications in the Web of Science database that centered on this field between 2007 and 2022. Next, we used VOSviewer, CiteSpace, Scimago Graphica, and Microsoft Excel to present visual representations of a co-occurrence analysis, co-citation analysis, tendencies, hotspots, and the contributions of authors, institutions, journals, and countries/regions. The United States was the main contributor. Investigative Ophthalmology and Visual Science has published the most articles in this field. Over the past 15 years, there has been exponential growth in the number of publications and citations in this field across various countries, organizations, and authors. Thus, this study illustrates the current trends, hotspots, and emerging frontiers and provides new insight and guidance for searching for new diagnostic biomarkers and clinical trials for glaucoma in the future. Furthermore, international collaborations can also be used to broaden and deepen the field of microRNAs in glaucoma research.

Integrated analysis of long non-coding RNAs and mRNAs associated with glaucoma in vitro

INTRODUCTION

In recent years, the biological functions and important roles of long non-coding RNAs (lncRNAs) have been widely reported in many diseases. Although glaucoma is the leading cause of blindness worldwide, the specific mechanisms of lncRNAs in the pathogenesis and progression of glaucoma remain unclear. Our research aims to elucidate the differentially expressed lncRNAs and mRNAs in glaucoma and to provide a basis for further exploration of the specific mechanism of action of lncRNAs in the progression of glaucoma.

METHODS

We performed RNA sequencing on samples from a pressurized model of R28 cells and performed bioinformatics analyses on the sequencing results. The expression consistency of lncRNAs in clinical samples from patients with glaucoma or cataracts was detected using real-time quantitative polymerase chain reaction (RT-qPCR).

RESULTS

RNA sequencing results showed that lncRNAs in cluster 5 were upregulated with increasing stress after typing all significantly altered lncRNAs using k-means in a cellular stress model. KEGG analysis indicated that they were associated with neurodegenerative diseases. Differentially expressed lncRNAs were verified by RT-qPCR, and the lncRNA expression levels of AC120246.2 and XLOC_006247 were significantly higher in the aqueous humor (AH) of patients with glaucoma than in those with cataracts. For LOC102551819, there was almost no expression in the AH and trabecular meshwork in patients with glaucoma but high expression was observed in the iris.

CONCLUSION

Our research proposes potential diagnostic or intervention targets for clinical applications as well as a theoretical basis for more in-depth research on the function of lncRNAs in glaucoma.

Construction of ceRNA network and identification of hub genes in aniridia-associated keratopathy using bioinformatics analysis

Aniridia-associated keratopathy (AAK) is characteristic at ocular surface of aniridia caused by haploinsufficiency of PAX6. Competing endogenous RNA (ceRNA) has been reported to play an important role in various diseases, whereas its function on AAK is unclear. The microarray data of 20 AAK patients and 20 healthy people were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed lncRNAs, miRNAs, and mRNAs were analyzed using “limma” packages and weighted gene co-expression network analysis (WGCNA). A ceRNA network was constructed by Cytoscape 3.9.1, and miR-224-5p, miR-30a-5p, and miR-204-5p were at the center of the network. CIBERSORTx algorithm and ssGSEA analyses revealed that AAK was associated with immune cell infiltration, showing that activated Mast cells increased while resting Mast cells decreased and NK cells decreased in AAK. Type II INF Response, CCR, parainflammation, T cell co-stimulation, and APC co-stimulation of AAK patients differed from healthy individuals. Additionally, the ROC curve of five genes, MITF(AUC = 0.988), RHOB(AUC = 0.973), JUN(AUC = 0.953), PLAUR (AUC = 0.925), and ARG2 (AUC = 0.915) with high confidence in predicting AAK were identified. Gene set enrichment analysis (GSEA) analysis of hub genes enriched in the IL-17 signaling pathway.

Screening of primary open-angle glaucoma diagnostic markers based on immune-related genes and immune infiltration

Purpose

Primary open-angle glaucoma (POAG) continues to be a poorly understood disease. Although there were multiple researches on the identification of POAG biomarkers, few studies systematically revealed the immune-related cells and immune infiltration of POAG. Bioinformatics analyses of optic nerve (ON) and trabecular meshwork (TM) gene expression data were performed to further elucidate the immune-related genes of POAG and identify candidate target genes for treatment.

Methods

We performed a gene analysis of publicly available microarray data, namely, the GSE27276-GPL2507, GSE2378-GPL8300, GSE9944-GPL8300, and GSE9944-GPL571 datasets from the Gene Expression Omnibus database. The obtained datasets were used as input for parallel pathway analyses. Based on random forest and support vector machine (SVM) analysis to screen the key genes, significantly changed pathways were clustered into functional categories, and the results were further investigated. CIBERSORT was used to evaluate the infiltration of immune cells in POAG tissues. A network visualizing the differences between the data in the POAG and normal groups was created. GO and KEGG enrichment analyses were performed using the Metascape database. We divided the differentially expressed mRNAs into upregulated and downregulated groups and predicted the drug targets of the differentially expressed genes through the Connectivity Map (CMap) database.

Results

A total of 49 differentially expressed genes, including 19 downregulated genes and 30 upregulated genes, were detected. Five genes ((Keratin 14) KRT14, (Hemoglobin subunit beta) HBB, (Acyl-CoA Oxidase 2) ACOX2, (Hephaestin) HEPH and Keratin 13 (KRT13)) were significantly changed. The results showed that the expression profiles of drug disturbances, including those for avrainvillamide-analysis-3, cytochalasin-D, NPI-2358, oxymethylone and vinorelbine, were negatively correlated with the expression profiles of disease disturbances. This finding indicated that these drugs may reduce or even reverse the POAG disease state.

Conclusion

This study provides an overview of the processes involved in the molecular pathogenesis of POAG in the ON and TM. The findings provide a new understanding of the molecular mechanism of POAG from the perspective of immunology.

CircXPO5 Plays a Neuroprotective Function in the Lateral Geniculate Nucleus of Glaucoma by Regulating GRIN2A

PURPOSE

Previous studies have found the neurodegeneration and atrophy of glaucomatous lateral geniculate nucleus (LGN), but the mechanism is still unknown. Circular RNA (circRNA) plays some important roles in physiological and pathological progression of the disease. In this study, we focused on the differentially expressed circRNAs and the mechanism for circXPO5 in LGN degeneration in a macaque glaucoma model.

METHODS

Using RNA-seq, we analyzed the differentially expressed circRNAs in a macaque glaucoma model. An RT-QPCR was used to check the expression of selected differentially expressed circRNAs, candidate miRNAs and mRNAs. A competing endogenous RNA (ceRNA) network analysis was performed to examine the mechanism of circXPO5 action.

RESULTS

circXPO5 significantly decreased in the glaucoma model and a ceRNA network analysis revealed that circXPO5 can bind to miR-330-5p, which also binds to GRIN2A (ionotropic receptor NMDA type subunit 2A). QPCR detection showed a decrease in GRIN2A and an increase in miR-330-5p.

CONCLUSIONS

Our earlier studies revealed that the GRIN2A gene regulates the calcium signal pathway. Decreasing of GRIN2A related with neuron apoptosis and neurodegeneration. These findings indicate that the reduction in circXPO5 may have a protective effect on neuronal apoptosis in the visual central system of glaucoma.

Identification of the potential biological target molecules related to primary open-angle glaucoma

Background

To identify the potential biological target molecules and the corresponding interaction networks in primary open-angle glaucoma (POAG) development.

Methods

The microarray datasets of GSE138125 and GSE27276 concerning lncRNA and mRNA expression profiles in trabecular meshwork of POAG were downloaded from the Gene Expression Omnibus database. The R software was applied to identify differentially expressed (DE) lncRNAs and mRNAs in POAG, and to perform GO and KEGG functional enrichment analysis. Protein–protein interaction (PPI) network and module analysis, and lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network were performed by Cytoscape software.

Results

A total of 567 DE-mRNAs were identified from GSE138125 and GSE27276, including 298 up-regulated and 269 down-regulated mRNAs, which were found enriching in biological processes of extracellular matrix organization and epidermis development, respectively. KEGG pathway enrichment analysis further revealed that module genes in PPI network were primarily involved in the AGE-PAGE, PI3K-Akt and TGF-β signaling pathways. Moreover, 897 up-regulated and 1036 down-regulated DE-lncRNAs were identified from GSE138125. Through literature review and databases searching, we obtained 712 lncRNA-miRNA and 337 miRNA-mRNA pairs based on the selected eight POAG-related miRNAs. After excluding 702 lncRNAs and 284 mRNAs that were not comprised in the DE-lncRNA and DE-mRNAs, a total of 53 lncRNA nodes, eight miRNA nodes, 10 mRNA nodes, and 78 edges were included in the final ceRNA network.

Conclusions

This study demonstrated the lncRNA and mRNA expression profiles of trabecular meshwork in POAG patients and the normal controls, and identified potentially ceRNAs and pathways which might improve the pathogenic understanding of this ocular disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12886-022-02368-0.