Mapping mRNA Expression of Glaucoma Genes in the Healthy Mouse Eye

Exploring single-cell RNA sequencing as a decision-making tool in the clinical management of Fuchs' endothelial corneal dystrophy

Single-cell RNA sequencing (scRNA-seq) has enabled the identification of novel gene signatures and cell heterogeneity in numerous tissues and diseases. Here we review the use of this technology for Fuchs' Endothelial Corneal Dystrophy (FECD). FECD is the most common indication for corneal endothelial transplantation worldwide. FECD is challenging to manage because it is genetically heterogenous, can be autosomal dominant or sporadic, and progress at different rates. Single-cell RNA sequencing has enabled the discovery of several FECD subtypes, each with associated gene signatures, and cell heterogeneity. Current FECD treatments are mainly surgical, with various Rho kinase (ROCK) inhibitors used to promote endothelial cell metabolism and proliferation following surgery. A range of emerging therapies for FECD including cell therapies, gene therapies, tissue engineered scaffolds, and pharmaceuticals are in preclinical and clinical trials. Unlike conventional disease management methods based on clinical presentations and family history, targeting FECD using scRNA-seq based precision-medicine has the potential to pinpoint the disease subtypes, mechanisms, stages, severities, and help clinicians in making the best decision for surgeries and the applications of therapeutics. In this review, we first discuss the feasibility and potential of using scRNA-seq in clinical diagnostics for FECD, highlight advances from the latest clinical treatments and emerging therapies for FECD, integrate scRNA-seq results and clinical notes from our FECD patients and discuss the potential of applying alternative therapies to manage these cases clinically.

Mitochondrial DNA D-loop variants correlate with a primary open-angle glaucoma subgroup

Introduction

Primary open-angle glaucoma (POAG) is a characteristic optic neuropathy, caused by degeneration of the optic nerve-forming neurons, the retinal ganglion cells (RGCs). High intraocular pressure (IOP) and aging have been identified as major risk factors; yet the POAG pathophysiology is not fully understood. Since RGCs have high energy requirements, mitochondrial dysfunction may put the survivability of RGCs at risk. We explored in buffy coat DNA whether mtDNA variants and their distribution throughout the mtDNA could be risk factors for POAG.

Methods

The mtDNA was sequenced from age- and sex-matched study groups, being high tension glaucoma (HTG, n=71), normal tension glaucoma patients (NTG, n=33), ocular hypertensive subjects (OH, n=7), and cataract controls (without glaucoma; n=30), all without remarkable comorbidities.

Results

No association was found between the number of mtDNA variants in genes encoding proteins, tRNAs, rRNAs, and in non-coding regions in the different study groups. Next, variants that controls shared with the other groups were discarded. A significantly higher number of exclusive variants was observed in the D-loop region for the HTG group (~1.23 variants/subject), in contrast to controls (~0.35 variants/subject). In the D-loop, specifically in the 7S DNA sub-region within the Hypervariable region 1 (HV1), we found that 42% of the HTG and 27% of the NTG subjects presented variants, while this was only 14% for the controls and OH subjects. As we have previously reported a reduction in mtDNA copy number in HTG, we analysed if specific D-loop variants could explain this. While the majority of glaucoma patients with the exclusive D-loop variants m.72T>C, m.16163 A>G, m.16186C>T, m.16298T>C, and m.16390G>A presented a mtDNA copy number below controls median, no significant association between these variants and low copy number was found and their possible negative role in mtDNA replication remains uncertain. Approximately 38% of the HTG patients with reduced copy number did not carry any exclusive D-loop or other mtDNA variants, which indicates that variants in nuclear-encoded mitochondrial genes, environmental factors, or aging might be involved in those cases.

Conclusion

In conclusion, we found that variants in the D-loop region may be a risk factor in a subgroup of POAG, possibly by affecting mtDNA replication.

Low mitochondrial DNA copy number in buffy coat DNA of primary open-angle glaucoma patients

Primary open-angle glaucoma (POAG) is characterized by optic nerve degeneration and irreversible loss of retinal ganglion cells (RGCs). The pathophysiology is not fully understood. Since RGCs have a high energy demand, suboptimal mitochondrial function may put the survival of these neurons at risk. In the present study, we explored whether mtDNA copy number or mtDNA deletions could reveal a mitochondrial component in POAG pathophysiology. Buffy coat DNA was isolated from EDTA blood of age- and sex-matched study groups, namely POAG patients with high intraocular pressure (IOP) at diagnosis (high tension glaucoma: HTG; n = 97), normal tension glaucoma patients (NTG, n = 37), ocular hypertensive controls (n = 9), and cataract controls (without glaucoma; n = 32), all without remarkable comorbidities. The number of mtDNA copies was assessed through qPCR quantification of the mitochondrial D-loop and nuclear B2M gene. Presence of the common 4977 base pair mtDNA deletion was assessed by a highly sensitive breakpoint PCR. Analysis showed that HTG patients had a lower number of mtDNA copies per nuclear DNA than NTG patients (p-value <0.01, Dunn test) and controls (p-value <0.001, Dunn test). The common 4977 base pair mtDNA deletion was not detected in any of the participants. A lower mtDNA copy number in blood of HTG patients suggests a role for a genetically defined, deficient mtDNA replication in the pathology of HTG. This may cause a low number of mtDNA copies in RGCs, which together with aging and high IOP, may lead to mitochondrial dysfunction, and contribute to glaucoma pathology.

Bioinformatic Prioritization and Functional Annotation of GWAS-Based Candidate Genes for Primary Open-Angle Glaucoma

Background: Primary open-angle glaucoma (POAG) is the most prevalent glaucoma subtype, but its exact etiology is still unknown. In this study, we aimed to prioritize the most likely ‘causal’ genes and identify functional characteristics and underlying biological pathways of POAG candidate genes. Methods: We used the results of a large POAG genome-wide association analysis study from GERA and UK Biobank cohorts. First, we performed systematic gene-prioritization analyses based on: (i) nearest genes; (ii) nonsynonymous single-nucleotide polymorphisms; (iii) co-regulation analysis; (iv) transcriptome-wide association studies; and (v) epigenomic data. Next, we performed functional enrichment analyses to find overrepresented functional pathways and tissues. Results: We identified 142 prioritized genes, of which 64 were novel for POAG. BICC1, AFAP1, and ABCA1 were the most highly prioritized genes based on four or more lines of evidence. The most significant pathways were related to extracellular matrix turnover, transforming growth factor-β, blood vessel development, and retinoic acid receptor signaling. Ocular tissues such as sclera and trabecular meshwork showed enrichment in prioritized gene expression (>1.5 fold). We found pleiotropy of POAG with intraocular pressure and optic-disc parameters, as well as genetic correlation with hypertension and diabetes-related eye disease. Conclusions: Our findings contribute to a better understanding of the molecular mechanisms underlying glaucoma pathogenesis and have prioritized many novel candidate genes for functional follow-up studies.

Effect of NK-5962 on Gene Expression Profiling of Retina in a Rat Model of Retinitis Pigmentosa

Purpose: NK-5962 is a key component of photoelectric dye-coupled polyethylene film, designated Okayama University type-retinal prosthesis (OUReP™). Previously, we found that NK-5962 solution could reduce the number of apoptotic photoreceptors in the eyes of the Royal College of Surgeons (RCS) rats by intravitreal injection under a 12 h light/dark cycle. This study aimed to explore possible molecular mechanisms underlying the anti-apoptotic effect of NK-5962 in the retina of RCS rats. Methods: RCS rats received intravitreal injections of NK-5962 solution in the left eye at the age of 3 and 4 weeks, before the age of 5 weeks when the speed in the apoptotic degeneration of photoreceptors reaches its peak. The vehicle-treated right eyes served as controls. All rats were housed under a 12 h light/dark cycle, and the retinas were dissected out at the age of 5 weeks for RNA sequence (RNA-seq) analysis. For the functional annotation of differentially expressed genes (DEGs), the Metascape and DAVID databases were used. Results: In total, 55 up-regulated DEGs, and one down-regulated gene (LYVE1) were found to be common among samples treated with NK-5962. These DEGs were analyzed using Gene Ontology (GO) term enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome pathway analyses. We focused on the up-regulated DEGs that were enriched in extracellular matrix organization, extracellular exosome, and PI3K–Akt signaling pathways. These terms and pathways may relate to mechanisms to protect photoreceptor cells. Moreover, our analyses suggest that SERPINF1, which encodes pigment epithelium-derived factor (PEDF), is one of the key regulatory genes involved in the anti-apoptotic effect of NK-5962 in RCS rat retinas. Conclusions: Our findings suggest that photoelectric dye NK-5962 may delay apoptotic death of photoreceptor cells in RCS rats by up-regulating genes related to extracellular matrix organization, extracellular exosome, and PI3K–Akt signaling pathways. Overall, our RNA-seq and bioinformatics analyses provide insights in the transcriptome responses in the dystrophic RCS rat retinas that were induced by NK-5962 intravitreal injection and offer potential target genes for developing new therapeutic strategies for patients with retinitis pigmentosa.

Plasma GDF-15 concentration is not elevated in open-angle glaucoma

Aim

Recently, the level of growth differentiation factor 15 (GDF-15) in blood, was proposed as biomarker to detect mitochondrial dysfunction. In the current study, we evaluate this biomarker in open-angle glaucoma (OAG), as there is increasing evidence that mitochondrial dysfunction plays a role in the pathophysiology of this disease.

Methods

Plasma GDF-15 concentrations were measured with ELISA in 200 OAG patients and 61 age-matched controls (cataract without glaucoma). The OAG patient group consisted of high tension glaucoma (HTG; n = 162) and normal tension glaucoma (NTG; n = 38). Groups were compared using the Kruskal-Wallis nonparametric test with Dunn’s multiple comparison post-hoc correction. GDF-15 concentration was corrected for confounders identified with forward linear regression models.

Results

Before correcting for confounders, median plasma GDF-15 levels was significantly lower in the combined OAG group (p = 0.04), but not when analysing HTG and NTG patients separately. Forward linear regression analysis showed that age, gender, smoking and systemic hypertension were significant confounders affecting GDF-15 levels. After correction for these confounders, GDF-15 levels in OAG patients were no longer significantly different from controls. Subgroup analysis of the glaucoma patients did not show a correlation between disease severity and plasma GDF-15, but did reveal that for NTG patients, intake of dietary supplements, which potentially improve mitochondrial function, correlated with lower plasma GDF-15.

Conclusion

The present study suggests that plasma GDF-15 is not suited as biomarker of mitochondrial dysfunction in OAG patients.