A Genome-Wide Scan for MicroRNA-Related Genetic Variants Associated With Primary Open-Angle Glaucoma

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank

BACKGROUND

Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown.

METHODS

Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis.

RESULTS

We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia.

CONCLUSIONS

Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Genetic Regulation of Human isomiR Biogenesis

MicroRNAs play a critical role in regulating gene expression post-transcriptionally. Variations in mature microRNA sequences, known as isomiRs, arise from imprecise cleavage and nucleotide substitution or addition. These isomiRs can target different mRNAs or compete with their canonical counterparts, thereby expanding the scope of miRNA post-transcriptional regulation. Our study investigated the relationship between cis-acting single-nucleotide polymorphisms (SNPs) in precursor miRNA regions and isomiR composition, represented by the ratio of a specific 5'-isomiR subtype to all isomiRs identified for a particular mature miRNA. Significant associations between 95 SNP-isomiR pairs were identified. Of note, rs6505162 was significantly associated with both the 5'-extension of hsa-miR-423-3p and the 5'-trimming of hsa-miR-423-5p. Comparison of breast cancer and normal samples revealed that the expression of both isomiRs was significantly higher in tumors than in normal tissues. This study sheds light on the genetic regulation of isomiR maturation and advances our understanding of post-transcriptional regulation by microRNAs.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

The 3' UTR polymorphisms rs3742330 in DICER1 and rs10719 in DROSHA genes are not associated with primary open-angle and angle-closure glaucoma: As case-control study

AIM

In a retrospective and exploratory case-control study, we examined the genetic association of two common polymorphisms in the 3' untranslated region (UTR) of DICER1 (rs3742330) and DROSHA (rs10719) genes in primary open-angle glaucoma (POAG) and primary angle-closure glaucoma (PACG), and its related clinical phenotypes in a Saudi cohort.

METHODS

DNA genotyping was performed using TaqMan real-time PCR assays in 500 participants, including 152 POAG, 102 PACG, and 246 non-glaucomatous controls. Statistical analyses were performed to examine the association(s).

RESULTS

Allele and genotype frequency of rs3742330 and rs10719 did not vary significantly in POAG and PACG compared to controls. No significant deviation was observed from Hardy-Weinberg Equilibrium (p > 0.05). Gender stratification revealed no significant allelic/genotype association with glaucoma types. Also, these polymorphisms showed no significant genotype effect on clinical markers such as intraocular pressure, cup/disc ratio, and the number of antiglaucoma medications. Logistic regression showed no effect of age, sex, rs3742330, and rs10719 genotypes on the risk of disease outcome. We also examined a combined allelic effect of rs3742330 (A>G) and rs10719 (A>G). However, none of the allelic combinations significantly affected POAG and PACG.

CONCLUSIONS

The 3' UTR polymorphisms rs3742330 and rs10719 of DICER1 and DROSHA genes are not associated with POAG and PACG or its related glaucoma indices in this Middle-Eastern cohort of Saudi Arab ethnicity. However, there is a need to validate the results on a broader population and other ethnicities.

In silico prioritisation of microRNA-associated common variants in multiple sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls.

METHODS

We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3'UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3'UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE.

RESULTS

We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3'UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites.

CONCLUSIONS

We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3'UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3'UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Signatures of genetic variation in human microRNAs point to processes of positive selection and population-specific disease risks

The occurrence of natural variation in human microRNAs has been the focus of numerous studies during the last 20 years. Most of them have been focused on the role of specific mutations in disease, while a minor proportion seek to analyse microRNA diversity in the genomes of human populations. We analyse the latest human microRNA annotations in the light of the most updated catalogue of genetic variation provided by the 1000 Genomes Project. By means of the in silico analysis of microRNA genetic variation we show that the level of evolutionary constraint of these sequences is governed by the interplay of different factors, like their evolutionary age or genomic location. The role of mutations in the shaping of microRNA-driven regulatory interactions is emphasized with the acknowledgement that, while the whole microRNA sequence is highly conserved, the seed region shows a pattern of higher genetic diversity that appears to be caused by the dramatic frequency shifts of a fraction of human microRNAs. We highlight the participation of these microRNAs in population-specific processes by identifying that not only the seed, but also the loop, are particularly differentiated regions among human populations. The quantitative computational comparison of signatures of population differentiation showed that candidate microRNAs with the largest differences are enriched in variants implicated in gene expression levels (eQTLs), selective sweeps and pathological processes. We explore the implication of these evolutionary-driven microRNAs and their SNPs in human diseases, such as different types of cancer, and discuss their role in population-specific disease risk.

Mechanobiological responses of astrocytes in optic nerve head due to biaxial stretch

Background

Elevated intraocular pressure (IOP) is the main risk factor for glaucoma, which might cause the activation of astrocytes in optic nerve head. To determine the effect of mechanical stretch on the astrocytes, we investigated the changes in cell phenotype, proteins of interest and signaling pathways under biaxial stretch.

Method

The cultured astrocytes in rat optic nerve head were stretched biaxially by 10 and 17% for 24 h, respectively. Then, we detected the morphology, proliferation and apoptosis of the stretched cells, and performed proteomics analysis. Protein expression was analyzed by Isobaric tags for relative and absolute quantification (iTRAQ) mass spectrometry. Proteins of interest and signaling pathways were screened using Gene Ontology enrichment analysis and pathway enrichment analysis, and the results were verified by western blot and the gene-chip data from Gene Expression Omnibus (GEO) database.

Result

The results showed that rearrangement of the actin cytoskeleton in response to stimulation by mechanical stress and proliferation rate of astrocytes decreased under 10 and 17% stretch condition, while there was no significant difference on the apoptosis rate of astrocytes in both groups. In the iTRAQ quantitative experiment, there were 141 differential proteins in the 10% stretch group and 140 differential proteins in the 17% stretch group. These proteins include low-density lipoprotein receptor-related protein (LRP6), caspase recruitment domain family, member 10 (CARD10), thrombospondin 1 (THBS1) and tetraspanin (CD81). The western blot results of LRP6, THBS1 and CD81 were consistent with that of iTRAQ experiment. ANTXR2 and CARD10 were both differentially expressed in the mass spectrometry results and GEO database. We also screened out the signaling pathways associated with astrocyte activation, including Wnt/β–catenin pathway, NF-κB signaling pathway, PI3K-Akt signaling pathway, MAPK signaling pathway, Jak-STAT signaling pathway, ECM-receptor interaction, and transforming growth factor-β (TGF-β) signaling pathway.

Conclusion

Mechanical stimulation can induce changes in cell phenotype, some proteins and signaling pathways, which might be associated with astrocyte activation. These proteins and signaling pathways may help us have a better understanding on the activation of astrocytes and the role astrocyte activation played in glaucomatous optic neuropathy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12886-022-02592-8.

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.

MicroRNAs in the aqueous humor of patients with different types of glaucoma

Purpose

The aim of the study was to compare the frequency and the level of expression of selected miRNAs in the aqueous humor of patients with various types of glaucoma.

Methods

The studied group consisted of 42 patients with glaucoma: 19 with primary open-angle glaucoma (POAG), 14 with pseudoexfoliation glaucoma (PEXG), 9 with primary angle closure glaucoma (PACG), and the control group of 36 patients with senile cataract without glaucoma. The real-time polymerase chain reaction method was used to analyze the expression of miRNAs.

Results

There were no significant differences in the frequency and the level of miRNA expression between various types of glaucoma. There was a tendency for hsa-miR-6722-3p and hsa-miR-184 to be expressed more frequently in PEXG and hsa-miR-1260b in POAG. The expression levels of hsa-miR-1260b and hsa-miR-6515-3p were correlated with age in POAG. Target annotation and functional analyses showed that genes targeted by the most frequently expressed miRNAs (hsa-miR-1202, -1260b, -184, -187-5p, -6515-3p, -6722-3p, and hsa-mir-4634) are involved mainly in response to hypoxia, cardiovascular system development, and apoptosis.

Conclusion

Hsa-miR-1260b was the most abundantly expressed among studied miRNAs and may be a potential biomarker of clinical status in PEXG and PACG.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00417-021-05214-z.

Crosstalk between MicroRNA and Oxidative Stress in Primary Open-Angle Glaucoma

Reactive oxygen species (ROS) plays a key role in the pathogenesis of primary open-angle glaucoma (POAG), a chronic neurodegenerative disease that damages the trabecular meshwork (TM) cells, inducing apoptosis of the retinal ganglion cells (RGC), deteriorating the optic nerve head, and leading to blindness. Aqueous humor (AH) outflow resistance and intraocular pressure (IOP) elevation contribute to disease progression. Nevertheless, despite the existence of pharmacological and surgical treatments, there is room for the development of additional treatment approaches. The following review is aimed at investigating the role of different microRNAs (miRNAs) in the expression of genes and proteins involved in the regulation of inflammatory and degenerative processes, focusing on the delicate balance of synthesis and deposition of extracellular matrix (ECM) regulated by chronic oxidative stress in POAG related tissues. The neutralizing activity of a couple of miRNAs was described, suggesting effective downregulation of pro-inflammatory and pro-fibrotic signaling pathways, including nuclear factor kappa-light-chain-enhancer of activated B cells (NF-kB), transforming growth factor-beta 2 (TGF-β2), Wnt/β-Catenin, and PI3K/AKT. In addition, with regards to the elevated IOP in many POAG patients due to increased outflow resistance, Collagen type I degradation was stimulated by some miRNAs and prevented ECM deposition in TM cells. Mitochondrial dysfunction as a consequence of oxidative stress was suppressed following exposure to different miRNAs. In contrast, increased oxidative damage by inhibiting the mTOR signaling pathway was described as part of the action of selected miRNAs. Summarizing, specific miRNAs may be promising therapeutic targets for lowering or preventing oxidative stress injury in POAG patients.

Role of Non-Coding Variants in Brugada Syndrome

Brugada syndrome (BrS) is an inherited electrical heart disease associated with a high risk of sudden cardiac death (SCD). The genetic characterization of BrS has always been challenging. Although several cardiac ion channel genes have been associated with BrS, SCN5A is the only gene that presents definitive evidence for causality to be used for clinical diagnosis of BrS. However, more than 65% of diagnosed cases cannot be explained by variants in SCN5A or other genes. Therefore, in an important number of BrS cases, the underlying mechanisms are still elusive. Common variants, mostly located in non-coding regions, have emerged as potential modulators of the disease by affecting different regulatory mechanisms, including transcription factors (TFs), three-dimensional organization of the genome, or non-coding RNAs (ncRNAs). These common variants have been hypothesized to modulate the interindividual susceptibility of the disease, which could explain incomplete penetrance of BrS observed within families. Altogether, the study of both common and rare variants in parallel is becoming increasingly important to better understand the genetic basis underlying BrS. In this review, we aim to describe the challenges of studying non-coding variants associated with disease, re-examine the studies that have linked non-coding variants with BrS, and provide further evidence for the relevance of regulatory elements in understanding this cardiac disorder.

An update: mechanisms of microRNA in primary open-angle glaucoma

Glaucoma is a disease with characteristic optic neuropathy and loss of vision, leading to blindness, and primary open-angle glaucoma (POAG) is the most common glaucoma type throughout the world. Genetic susceptibility is the main factor in POAG, and most susceptibility genes cause changes in microRNA expression and function, thereby leading to POAG occurrence and development. Increasing evidence indicates that many microRNAs are involved in the regulation of intraocular pressure (IOP) and play an important role in the increase in IOP in POAG. Additionally, microRNA is closely related to optic nerve damage factors (mechanical stress, hypoxia and inflammation). This review discusses the effect of single-nucleotide polymorphisms in POAG-related genes on microRNA and the value of microRNA in the diagnosis and treatment of POAG.

Differential miRNA Expression: Signature for Glaucoma in Pseudoexfoliation

Purpose

To investigate the microRNA (miRNA) profile in patients with different stages of pseudoexfoliation (PXF).

Methods

Peripheral blood of patients with PXF (naïve to medical therapy and with no systemic disease/drugs) with ocular hypertension (OHT) and pseudoexfoliation glaucoma (PXG) was evaluated in triplicate for miRNA profiling using polymerase chain reaction (PCR) arrays. Those identified in the discovery stage were validated with evaluation of serum transforming growth factor-β1 (TGF-β1) levels by ELISA. The downstream targets of TGF-β1 and unfolded protein response (UPR) were analyzed using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Predicted targets of the identified miRNA and KEGG pathway analysis were done using miRbase and DIANA tools mirPathv3.1.

Results

We found hsa-miR-122-5p, hsa-miR-124-3p and hsa-miR-424-5p to be upregulated in PXG targeting 3 specific pathways namely TGF-β1, fibrosis/ECM and proteoglycan metabolism with common effectors like SMAD/3/2. The unfolded protein response (UPR) genes were significantly downregulated in all stages of PXF suggesting this as the key mechanism for protein aggregates in PXF syndrome. Serum TGF-β1 was significantly upregulated as disease progressed to later stages in PXG. This elevation in advanced stages was associated with significantly differential expression of downstream pathways and fibrotic genes in OHT compared to PXG predominantly through the SMAD3, a canonical pathway marker.

Conclusion

Circulatory miRNA differentially regulating TGF-β1 and downstream targets including UPR genes may be the key mechanisms for glaucoma onset in PXF.

Polymorphism analysis of miR182 and CDKN2B genes in Greek patients with primary open angle glaucoma

Glaucoma is a progressive optic neuropathy resulting from retinal ganglion cells death; it represents one of the leading causes of irreversible blindness worldwide. Although, primary open angle glaucoma (POAG) is the most common type of the disease, the pathogenesis of POAG and the genetic factors contributing to disease development remain poorly understood. The aim of this study was to investigate whether the polymorphisms rs76481776 in miR182 gene and rs3217992 in cyclin-dependent kinase inhibitor-2B (CDKN2B) gene are risk factors for POAG in a series of patients of Greek origin. A case-control study was conducted including 120 patients with POAG and 113 unaffected healthy controls of Greek origin, surveyed for polymorphisms with potential correlation to POAG. DNA from each individual was tested for the miR182 rs76481776 and CDKN2B rs3217992 polymorphisms. Regarding the miR182 rs76481776 polymorphism, the T allele occurred with significantly higher frequency in POAG patients compared to controls (OR: 2.62, 95% CI: 1.56-4.39; p = 0.0002). The CDKN2B rs3217992 A allele frequency was found significantly increased in POAG patients compared to healthy individuals (OR: 1.72, 95% CI: 1.18-2.49; p = 0.005). Therefore, both rs76481776 polymorphism in miR182 gene and rs3217992 polymorphism in CDKN2B gene seem to be associated with the development of POAG in a Greek population. The carriers of the T allele of rs76481776 in miR182 and the carriers of the A allele of rs3217992 in CDKN2B have an increased risk of developing POAG.

Objectives, design and main findings until 2020 from the Rotterdam Study

The Rotterdam Study is an ongoing prospective cohort study that started in 1990 in the city of Rotterdam, The Netherlands. The study aims to unravel etiology, preclinical course, natural history and potential targets for intervention for chronic diseases in mid-life and late-life. The study focuses on cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1700 research articles and reports. This article provides an update on the rationale and design of the study. It also presents a summary of the major findings from the preceding 3 years and outlines developments for the coming period.

Update on the genetics of primary open-angle glaucoma

Affecting nearly 80 million individuals, glaucoma is the number one cause of irreversible blindness in the world. This ocular disease describes a set of optic neuropathies of which primary open angle glaucoma (POAG) is the most common. POAG is associated with progressive visual field deterioration resulting from damage to the optic nerve and loss of retinal ganglion cells. Risk factors for POAG include elevated intraocular pressure, aging, African and Hispanic ancestry, and a positive family history of POAG. Multiple genes have been found to contribute to POAG. Much of POAG genetics and pathology has yet to be explained. Recent genome-wide association studies have identified a large number of novel loci associated with POAG and its endophenotypes. Genomic and proteomic profiling of biofluids has contributed to our knowledge of differential gene expression in POAG. Functional studies both in cell culture and animal models have confirmed the effects of variants and differential gene expression on ocular physiology while in silico analyses have increased our understanding of disease risk and progression so that we might better diagnose and treat this complex genetic illness.

[Molecular and genetic aspects of age-related macular degeneration and glaucoma]

In most cases, age-related macular degeneration (AMD) and glaucoma are considered multi-factor diseases that lead to irreversible blindness in senior population of developed countries. Among different types of these diseases, around 5% are monogenic. Studying their molecular and genetic aspects can lay the basis for improvement of diagnostic methods, prognosis of the risks of development, manner of progression and treatment outcomes, as well as creation of new therapy methods. The article reviews modern understanding of the etiopathogenesis of AMD and glaucoma and describes their interrelations.

microRNA Profiling in Glaucoma Eyes With Varying Degrees of Optic Neuropathy by Using Next-Generation Sequencing

Purpose

To explore the microRNA (miRNA) profile and its putative role in glaucomatous optic neuropathy by using next-generation sequencing.

Methods

Aqueous humor (AH) samples were collected from 19 primary open-angle glaucoma (POAG) eyes and 17 cataract eyes before surgery. Next-generation sequencing was performed for RNA samples extracted from 18 AH samples, and the bioinformatics approach was applied for samples with adequate clean data output. The other 18 samples were used for quantitative PCR validation of sequencing results.

Results

In total, 12 (six POAG and six cataract controls) samples with sufficient clean data output after sequencing were used for further data analysis. Four hundred sixty-six and 480 mature miRNAs were detected in the POAG and cataract control groups, respectively. Among them, 164 miRNAs were detected in all POAG samples, and 96 miRNAs were detected in all cataract control samples. Furthermore, 88 miRNAs were identified as differently expressed between POAG and cataract control eyes. In addition, 16 miRNAs were differently expressed between POAG eyes with severe visual field damage and eyes with moderate visual field damage. This differential expression was predicted to regulate thiamine metabolism, purine metabolism, and transcriptional misregulation. Relative expression patterns of hsa-miR-184, hsa-miR-486-5p, and hsa-miR-93-5p were confirmed by quantitative PCR.

Conclusions

This study comprehensively demonstrated the miRNA expression profile in the AH of POAG eyes, especially the differential expression of miRNA in eyes with varying degrees of visual field damage, which, together with the underlying miRNA-related pathways, indicate new targets for the pathogenesis and progression of POAG.