Cell atlas of aqueous humor outflow pathways in eyes of humans and four model species provides insight into glaucoma pathogenesis

Single-cell RNA sequencing analysis of the retina under acute high intraocular pressure

JOURNAL/nrgr/04.03/01300535-202419110-00032/figure1/v/2024-03-08T184507Z/r/image-tiff High intraocular pressure causes retinal ganglion cell injury in primary and secondary glaucoma diseases, yet the molecular landscape characteristics of retinal cells under high intraocular pressure remain unknown. Rat models of acute hypertension ocular pressure were established by injection of cross-linked hyaluronic acid hydrogel (Healaflow®). Single-cell RNA sequencing was then used to describe the cellular composition and molecular profile of the retina following high intraocular pressure. Our results identified a total of 12 cell types, namely retinal pigment epithelial cells, rod-photoreceptor cells, bipolar cells, Müller cells, microglia, cone-photoreceptor cells, retinal ganglion cells, endothelial cells, retinal progenitor cells, oligodendrocytes, pericytes, and fibroblasts. The single-cell RNA sequencing analysis of the retina under acute high intraocular pressure revealed obvious changes in the proportions of various retinal cells, with ganglion cells decreased by 23%. Hematoxylin and eosin staining and TUNEL staining confirmed the damage to retinal ganglion cells under high intraocular pressure. We extracted data from retinal ganglion cells and analyzed the retinal ganglion cell cluster with the most distinct expression. We found upregulation of the B3gat2 gene, which is associated with neuronal migration and adhesion, and downregulation of the Tsc22d gene, which participates in inhibition of inflammation. This study is the first to reveal molecular changes and intercellular interactions in the retina under high intraocular pressure. These data contribute to understanding of the molecular mechanism of retinal injury induced by high intraocular pressure and will benefit the development of novel therapies.

Single-cell transcriptomic Atlas of aging macaque ocular outflow tissues

The progressive degradation in the trabecular meshwork (TM) is related to age-related ocular diseases like primary open-angle glaucoma. However, the molecular basis and biological significance of the aging process in TM have not been fully elucidated. Here, we established a dynamic single-cell transcriptomic landscape of aged macaque TM, wherein we classified the outflow tissue into 12 cell subtypes and identified mitochondrial dysfunction as a prominent feature of TM aging. Furthermore, we divided TM cells into 13 clusters and performed an in-depth analysis on cluster 0, which had the highest aging score and the most significant changes in cell proportions between the two groups. Ultimately, we found that the APOE gene was an important differentially expressed gene in cluster 0 during the aging process, highlighting the close relationship between cell migration and extracellular matrix regulation, and TM function. Our work further demonstrated that silencing the APOE gene could increase migration and reduce apoptosis by releasing the inhibition on the PI3K-AKT pathway and downregulating the expression of extracellular matrix components, thereby increasing the aqueous outflow rate and maintaining intraocular pressure within the normal range. Our work provides valuable insights for future clinical diagnosis and treatment of glaucoma.

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.

Single-cell RNA sequencing in exploring the pathogenesis of diabetic retinopathy

Diabetic retinopathy (DR) is a leading cause of irreversible blindness in the working-age populations. Despite decades of research on the pathogenesis of DR for clinical care, a comprehensive understanding of the condition is still lacking due to the intricate cellular diversity and molecular heterogeneity involved. Single-cell RNA sequencing (scRNA-seq) has made the high-throughput molecular profiling of cells across modalities possible which has provided valuable insights into complex biological systems. In this review, we summarise the application of scRNA-seq in investigating the pathogenesis of DR, focusing on four aspects. These include the identification of differentially expressed genes, characterisation of key cell subpopulations and reconstruction of developmental 'trajectories' to unveil their state transition, exploration of complex cell‒cell communication in DR and integration of scRNA-seq with genome-wide association studies to identify cell types that are most closely related to DR risk genetic loci. Finally, we discuss the future challenges and expectations associated with studying DR using scRNA-seq. We anticipate that scRNA-seq will facilitate the discovery of mechanisms and new treatment targets in the clinical care landscape for patients with DR. KEY POINTS: Progress in scRNA-seq for diabetic retinopathy (DR) research includes studies on DR patients, non-human primates, and the prevalent mouse models. scRNA-seq facilitates the identification of differentially expressed genes, pivotal cell subpopulations, and complex cell-cell interactions in DR at single-cell level. Future scRNA-seq applications in DR should target specific patient subsets and integrate with single-cell and spatial multi-omics approaches.

iPSC-derived cells stimulate ABCG2+/NES+ endogenous trabecular meshwork cell proliferation and tissue regeneration

A major risk factor for glaucoma, the first leading cause of irreversible blindness worldwide, is the decellularisation of the trabecular meshwork (TM) in the conventional outflow pathway. Stem cell-based therapy, particularly the utilisation of induced pluripotent stem cells (iPSCs), presents an enticing potential for tissue regeneration and intraocular pressure (IOP) maintenance in glaucoma. We have previously observed that differentiated iPSCs can stimulate endogenous cell proliferation in the TM, a pivotal factor in TM regeneration and aqueous humour outflow restoration. In this study, we investigated the response of TM cells in vivo after interacting with iPSC-derived cells and identified two subpopulations responsible for this relatively long-term tissue regeneration: ATP Binding Cassette Subfamily G Member 2 (ABCG2)-positive cells and Nestin (NES)-positive cells. We further uncovered that alterations of these responsive cells are linked to ageing and different glaucoma etiologies, suggesting that ABCG2+ subpopulation decellularization could serve as a potential risk factor for TM decellularization in glaucoma. Taken together, our findings illustrated the proliferative subpopulations in the conventional outflow pathway when stimulated with iPSC-derived cells and defined them as TM precursors, which may be applied to develop novel therapeutic approaches for glaucoma.

The future of rapid and automated single-cell data analysis using reference mapping

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

Addressing neurodegeneration in glaucoma: Mechanisms, challenges, and treatments

Glaucoma is the leading cause of irreversible blindness globally. The disease causes vision loss due to neurodegeneration of the retinal ganglion cell (RGC) projection to the brain through the optic nerve. Glaucoma is associated with sensitivity to intraocular pressure (IOP). Thus, mainstay treatments seek to manage IOP, though many patients continue to lose vision. To address neurodegeneration directly, numerous preclinical studies seek to develop protective or reparative therapies that act independently of IOP. These include growth factors, compounds targeting metabolism, anti-inflammatory and antioxidant agents, and neuromodulators. Despite success in experimental models, many of these approaches fail to translate into clinical benefits. Several factors contribute to this challenge. Firstly, the anatomic structure of the optic nerve head differs between rodents, nonhuman primates, and humans. Additionally, animal models do not replicate the complex glaucoma pathophysiology in humans. Therefore, to enhance the success of translating these findings, we propose two approaches. First, thorough evaluation of experimental targets in multiple animal models, including nonhuman primates, should precede clinical trials. Second, we advocate for combination therapy, which involves using multiple agents simultaneously, especially in the early and potentially reversible stages of the disease. These strategies aim to increase the chances of successful neuroprotective treatment for glaucoma.

To be or not to be - Decoding the Trabecular Meshwork Cell Identity

The trabecular meshwork within the conventional outflow apparatus is critical in maintaining intraocular pressure homeostasis. In vitro studies employing primary cell cultures of the human trabecular meshwork (hTM) have conventionally served as surrogates for investigating the pathobiology of TM dysfunction. Despite its abundant use, translation of outcomes from in vitro studies to ex vivo and/or in vivo studies remains a challenge. Given the cell heterogeneity, performing single-cell RNA sequencing comparing primary hTM cell cultures to hTM tissue may provide important insights on cellular identity and translatability, as such an approach has not been reported before. In this study, we assembled a total of 14 primary hTM in vitro samples across passages 1-4, including 4 samples from individuals diagnosed with glaucoma. This dataset offers a comprehensive transcriptomic resource of primary hTM in vitro scRNA-seq data to study global changes in gene expression in comparison to cells in tissue in situ. We have performed extensive preprocessing and quality control, allowing the research community to access and utilize this public resource.

High throughput functional profiling of genes at intraocular pressure loci reveals distinct networks for glaucoma

INTRODUCTION

Primary open angle glaucoma (POAG) is a leading cause of blindness globally. Characterized by progressive retinal ganglion cell degeneration, the precise pathogenesis remains unknown. Genome-wide association studies (GWAS) have uncovered many genetic variants associated with elevated intraocular pressure (IOP), one of the key risk factors for POAG. We aimed to identify genetic and morphological variation that can be attributed to trabecular meshwork cell (TMC) dysfunction and raised IOP in POAG.

METHODS

62 genes across 55 loci were knocked-out in a primary human TMC line. Each knockout group, including five non-targeting control groups, underwent single-cell RNA-sequencing (scRNA-seq) for differentially-expressed gene (DEG) analysis. Multiplexed fluorescence coupled with CellProfiler image analysis allowed for single-cell morphological profiling.

RESULTS

Many gene knockouts invoked DEGs relating to matrix metalloproteinases and interferon-induced proteins. We have prioritized genes at four loci of interest to identify gene knockouts that may contribute to the pathogenesis of POAG, including ANGPTL2, LMX1B, CAV1, and KREMEN1. Three genetic networks of gene knockouts with similar transcriptomic profiles were identified, suggesting a synergistic function in trabecular meshwork cell physiology. TEK knockout caused significant upregulation of nuclear granularity on morphological analysis, while knockout of TRIOBP, TMCO1 and PLEKHA7 increased granularity and intensity of actin and the cell-membrane.

CONCLUSION

High-throughput analysis of cellular structure and function through multiplex fluorescent single-cell analysis and scRNA-seq assays enabled the direct study of genetic perturbations at the single-cell resolution. This work provides a framework for investigating the role of genes in the pathogenesis of glaucoma and heterogenous diseases with a strong genetic basis.

Exploring the Genetic Landscape of Childhood Glaucoma

Childhood glaucoma, a significant cause of global blindness, represents a heterogeneous group of disorders categorized into primary or secondary forms. Primary childhood glaucoma stands as the most prevalent subtype, comprising primary congenital glaucoma (PCG) and juvenile open-angle glaucoma (JOAG). Presently, multiple genes are implicated in inherited forms of primary childhood glaucoma. This comprehensive review delves into genetic investigations into primary childhood glaucoma, with a focus on identifying causative genes, understanding their inheritance patterns, exploring essential biological pathways in disease pathogenesis, and utilizing animal models to study these mechanisms. Specifically, attention is directed towards genes such as CYP1B1 (cytochrome P450 family 1 subfamily B member 1), LTBP2 (latent transforming growth factor beta binding protein 2), TEK (TEK receptor tyrosine kinase), ANGPT1 (angiopoietin 1), and FOXC1 (forkhead box C1), all associated with PCG; and MYOC (myocilin), associated with JOAG. Through exploring these genetic factors, this review aims to deepen our understanding of the intricate pathogenesis of primary childhood glaucoma, thereby facilitating the development of enhanced diagnostic and therapeutic strategies.

Circadian rhythm of intraocular pressure

Intraocular pressure (IOP) plays a crucial role in glaucoma development, involving the dynamics of aqueous humor (AH). AH flows in from the ciliary body and exits through the trabecular meshwork (TM). IOP follows a circadian rhythm synchronized with the suprachiasmatic nucleus (SCN), the circadian pacemaker. The SCN resets peripheral clocks through sympathetic nerves or adrenal glucocorticoids (GCs). IOP's circadian rhythm is governed by circadian time signals, sympathetic noradrenaline (NE), and GCs, rather than the local clock. The activity of Na+/K+-ATPase in non-pigmented epithelial cells in the ciliary body can influence the nocturnal increase in IOP by enhancing AH inflow. Conversely, NE, not GCs, can regulate the IOP rhythm by suppressing TM macrophage phagocytosis and AH outflow. The activation of the β1-adrenergic receptor (AR)-mediated EPAC-SHIP1 signal through the ablation of phosphatidylinositol triphosphate may govern phagocytic cup formation. These findings could offer insights for better glaucoma management, such as chronotherapy.

Adult stem cells in the eye: Identification, characterisation, and therapeutic application in ocular regeneration - A review

Adult stem cells, present in various parts of the human body, are undifferentiated cells that can proliferate and differentiate to replace dying cells within tissues. Stem cells have specifically been identified in the cornea, trabecular meshwork, crystalline lens, iris, ciliary body, retina, choroid, sclera, conjunctiva, eyelid, lacrimal gland, and orbital fat. The identification of ocular stem cells broadens the potential therapeutic strategies for untreatable eye diseases. Currently, stem cell transplantation for corneal and conjunctival diseases remains the most common stem cell-based therapy in ocular clinical management. Lens epithelial stem cells have been applied in the treatment of paediatric cataracts. Several early-phase clinical trials for corneal and retinal regeneration using ocular stem cells are also underway. Extensive preclinical studies using ocular stem cells have been conducted, showing encouraging outcomes. Ocular stem cells currently demonstrate great promise in potential treatments of eye diseases. In this review, we focus on the identification, characterisation, and therapeutic application of adult stem cells in the eye.

Role of mechanically-sensitive cation channels Piezo1 and TRPV4 in trabecular meshwork cell mechanotransduction

Glaucoma is one of the leading causes of irreversible blindness in developed countries, and intraocular pressure (IOP) is primary and only treatable risk factor, suggesting that to a significant extent, glaucoma is a disease of IOP disorder and pathological mechanotransduction. IOP-lowering ways are limited to decreaseing aqueous humour (AH) production or increasing the uveoscleral outflow pathway. Still, therapeutic approaches have been lacking to control IOP by enhancing the trabecular meshwork (TM) pathway. Trabecular meshwork cells (TMCs) have endothelial and myofibroblast properties and are responsible for the renewal of the extracellular matrix (ECM). Mechanosensitive cation channels, including Piezo1 and TRPV4, are abundantly expressed in primary TMCs and trigger mechanostress-dependent ECM and cytoskeletal remodelling. However, prolonged mechanical stimulation severely affects cellular biosynthesis through TMC mechanotransduction, including signaling, gene expression, ECM remodelling, and cytoskeletal structural changes, involving outflow facilities and elevating IOP. As for the functional coupling relationship between Piezo1 and TRPV4 channels, inspired by VECs and osteoblasts, we hypothesized that Piezo1 may also act upstream of TRPV4 in glaucomatous TM tissue, mediating the activation of TRPV4 via Ca2+ inflow or Ca2+ binding to phospholipase A2(PLA2), and thus be involved in increasing TM outflow resistance and elevated IOP. Therefore, this review aims to help identify new potential targets for IOP stabilization in ocular hypertension and primary open-angle glaucoma by understanding the mechanical transduction mechanisms associated with the development of glaucoma and may provide ideas into novel treatments for preventing the progression of glaucoma by targeting mechanotransduction.

Toward universal cell embeddings: integrating single-cell RNA-seq datasets across species with SATURN

Analysis of single-cell datasets generated from diverse organisms offers unprecedented opportunities to unravel fundamental evolutionary processes of conservation and diversification of cell types. However, interspecies genomic differences limit the joint analysis of cross-species datasets to homologous genes. Here we present SATURN, a deep learning method for learning universal cell embeddings that encodes genes' biological properties using protein language models. By coupling protein embeddings from language models with RNA expression, SATURN integrates datasets profiled from different species regardless of their genomic similarity. SATURN can detect functionally related genes coexpressed across species, redefining differential expression for cross-species analysis. Applying SATURN to three species whole-organism atlases and frog and zebrafish embryogenesis datasets, we show that SATURN can effectively transfer annotations across species, even when they are evolutionarily remote. We also demonstrate that SATURN can be used to find potentially divergent gene functions between glaucoma-associated genes in humans and four other species.

Cell type evolution reconstruction across species through cell phylogenies of single-cell RNA sequencing data

The origin and evolution of cell types has emerged as a key topic in evolutionary biology. Driven by rapidly accumulating single-cell datasets, recent attempts to infer cell type evolution have largely been limited to pairwise comparisons because we lack approaches to build cell phylogenies using model-based approaches. Here we approach the challenges of applying explicit phylogenetic methods to single-cell data by using principal components as phylogenetic characters. We infer a cell phylogeny from a large, comparative single-cell dataset of eye cells from five distantly related mammals. Robust cell type clades enable us to provide a phylogenetic, rather than phenetic, definition of cell type, allowing us to forgo marker genes and phylogenetically classify cells by topology. We further observe evolutionary relationships between diverse vessel endothelia and identify the myelinating and non-myelinating Schwann cells as sister cell types. Finally, we examine principal component loadings and describe the gene expression dynamics underlying the function and identity of cell type clades that have been conserved across the five species. A cell phylogeny provides a rigorous framework towards investigating the evolutionary history of cells and will be critical to interpret comparative single-cell datasets that aim to ask fundamental evolutionary questions.

Transcriptomic profiling of Schlemm's canal cells reveals a lymphatic-biased identity and three major cell states

Schlemm's canal (SC) is central in intraocular pressure regulation but requires much characterization. It has distinct inner and outer walls, each composed of Schlemm's canal endothelial cells (SECs) with different morphologies and functions. Recent transcriptomic studies of the anterior segment added important knowledge, but were limited in power by SEC numbers or did not focus on SC. To gain a more comprehensive understanding of SC biology, we performed bulk RNA sequencing on C57BL/6J SC, blood vessel, and lymphatic endothelial cells from limbal tissue (~4500 SECs). We also analyzed mouse limbal tissues by single-cell and single-nucleus RNA sequencing (C57BL/6J and 129/Sj strains), successfully sequencing 903 individual SECs. Together, these datasets confirm that SC has molecular characteristics of both blood and lymphatic endothelia with a lymphatic phenotype predominating. SECs are enriched in pathways that regulate cell-cell junction formation pointing to the importance of junctions in determining SC fluid permeability. Importantly, and for the first time, our analyses characterize 3 molecular classes of SECs, molecularly distinguishing inner wall from outer wall SECs and discovering two inner wall cell states that likely result from local environmental differences. Further, and based on ligand and receptor expression patterns, we document key interactions between SECs and cells of the adjacent trabecular meshwork (TM) drainage tissue. Also, we present cell type expression for a collection of human glaucoma genes. These data provide a new molecular foundation that will enable the functional dissection of key homeostatic processes mediated by SECs as well as the development of new glaucoma therapeutics.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

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.

TIE1 and TEK signalling, intraocular pressure, and primary open-angle glaucoma: a Mendelian randomization study

BACKGROUND

In primary open-angle glaucoma (POAG), lowering intraocular pressure (IOP) is the only proven way of slowing vision loss. Schlemm's canal (SC) is a hybrid vascular and lymphatic vessel that mediates aqueous humour drainage from the anterior ocular chamber. Animal studies support the importance of SC endothelial angiopoietin-TEK signalling, and more recently TIE1 signalling, in maintaining normal IOP. However, human genetic support for a causal role of TIE1 and TEK signalling in lowering IOP is currently lacking.

METHODS

GWAS summary statistics were obtained for plasma soluble TIE1 (sTIE1) protein levels (N = 35,559), soluble TEK (sTEK) protein levels (N = 35,559), IOP (N = 139,555) and POAG (Ncases = 16,677, Ncontrols = 199,580). Mendelian randomization (MR) was performed to estimate the association of genetically proxied TIE1 and TEK protein levels with IOP and POAG liability. Where significant MR estimates were obtained, genetic colocalization was performed to assess the probability of a shared causal variant (PPshared) versus distinct (PPdistinct) causal variants underlying TIE1/TEK signalling and the outcome. Publicly available single-nucleus RNA-sequencing data were leveraged to investigate differential expression of TIE1 and TEK in the human ocular anterior segment.

RESULTS

Increased genetically proxied TIE1 signalling and TEK signalling associated with a reduction in IOP (- 0.21 mmHg per SD increase in sTIE1, 95% CI = - 0.09 to - 0.33 mmHg, P = 6.57 × 10-4, and - 0.14 mmHg per SD decrease in sTEK, 95% CI = - 0.03 to - 0.25 mmHg, P = 0.011), but not with POAG liability. Colocalization analysis found that the probability of a shared causal variant was greater for TIE1 and IOP than for TEK and IOP (PPshared/(PPdistinct + PPshared) = 0.98 for TIE1 and 0.30 for TEK). In the anterior segment, TIE1 and TEK were preferentially expressed in SC, lymphatic, and vascular endothelium.

CONCLUSIONS

This study provides novel human genetic support for a causal role of both TIE1 and TEK signalling in regulating IOP. Here, combined evidence from cis-MR and colocalization analyses provide stronger support for TIE1 than TEK as a potential IOP-lowering therapeutic target.

Cell atlas of trabecular meshwork in glaucomatous non-human primates and DEGs related to tissue contract based on single-cell transcriptomics

As the major channel of aqueous humor outflow, dysfunction of trabecular meshwork (TM) can lead to intraocular pressure elevating, which can trigger primary open-angle glaucoma (POAG). In this study, we use single-cell RNA sequencing (scRNA-seq) technique to build an atlas and further explore the spontaneous POAG and healthy macaques cellular heterogeneity associated with the dysfunction of TM contraction. We built the TM atlas, which identified 14 different cell types. In Beam A, Beam B, Beam C, and smooth muscle cell (SMC) cell types, we first found multiple genes associated with TM contraction (e.g., TPM1, ACTC1, TNNT1), determining their differential expression in the POAG and healthy groups. In addition, the microstructural alterations in TM of POAG non-human primates were observed, which was compact and collapsed. Thus, our study indicated that TPM1 may be a key target for regulating TM structure, contraction function, and resistance of aqueous humor outflow.

Glaucoma as a Tauopathy-Is It the Missing Piece in the Glaucoma Puzzle?

Glaucoma is a chronic neurodegenerative disorder affecting the visual system which can result in vision loss and blindness. The pathogenetic mechanisms underlying glaucomatous optic neuropathy are ultimately enigmatic, prompting ongoing investigations into its potential shared pathogenesis with other neurodegenerative neurological disorders. Tauopathies represent a subclass of neurodegenerative diseases characterized by the abnormal deposition of tau protein within the brain and consequent microtubule destabilization. The extended spectrum of tauopathies includes conditions such as frontotemporal dementias, progressive supranuclear palsy, chronic traumatic encephalopathy, and Alzheimer's disease. Notably, recent decades have witnessed emerging documentation of tau inclusion among glaucoma patients, providing substantiation that this ocular disease may similarly manifest features of tauopathies. These studies found that: (i) aggregated tau inclusions are present in the somatodendritic compartment of RGCs in glaucoma patients; (ii) the etiology of the disease may affect tau splicing, phosphorylation, oligomerization, and subcellular localization; and (iii) short interfering RNA against tau, administered intraocularly, significantly decreased retinal tau accumulation and enhanced RGC somas and axon survival, demonstrating a crucial role for tau modifications in ocular hypertension-induced neuronal injury. Here, we examine the most recent evidence surrounding the interplay between tau protein dysregulation and glaucomatous neurodegeneration. We explore the novel perspective of glaucoma as a tau-associated disorder and open avenues for cross-disciplinary collaboration and new treatment strategies.

Single-cell transcriptomics of the ocular anterior segment: a comprehensive review

Elucidating the cellular and genetic composition of ocular tissues is essential for uncovering the pathophysiology of ocular diseases. Since the introduction of single-cell RNA sequencing (scRNA-seq) in 2009, vision researchers have performed extensive single-cell analyses to better understand transcriptome complexity and heterogeneity of ocular structures. This technology has revolutionized our ability to identify rare cell populations and to make cross-species comparisons of gene expression in both steady state and disease conditions. Importantly, single-cell transcriptomic analyses have enabled the identification of cell-type specific gene markers and signalling pathways between ocular cell populations. While most scRNA-seq studies have been conducted on retinal tissues, large-scale transcriptomic atlases pertaining to the ocular anterior segment have also been constructed in the past three years. This timely review provides vision researchers with an overview of scRNA-seq experimental design, technical limitations, and clinical applications in a variety of anterior segment-related ocular pathologies. We review open-access anterior segment-related scRNA-seq datasets and illustrate how scRNA-seq can be an indispensable tool for the development of targeted therapeutics.

Differential expression of PIEZO1 and PIEZO2 mechanosensitive channels in ocular tissues implicates diverse functional roles

PIEZO1 and PIEZO2 are mechanosensitive ion channels that regulate many important physiological processes including vascular blood flow, touch, and proprioception. As the eye is subject to mechanical stress and is highly perfused, these channels may play important roles in ocular function and intraocular pressure regulation. PIEZO channel expression in the eye has not been well defined, in part due to difficulties in validating available antibodies against PIEZO1 and PIEZO2 in ocular tissues. It is also unclear if PIEZO1 and PIEZO2 are differentially expressed. To address these questions, we used single-molecule fluorescence in situ hybridization (smFISH) together with transgenic reporter mice expressing PIEZO fusion proteins under the control of their endogenous promoters to compare the expression and localization of PIEZO1 and PIEZO2 in mouse ocular tissues relevant to glaucoma. We detected both PIEZO1 and PIEZO2 expression in the trabecular meshwork, ciliary body, and in the ganglion cell layer (GCL) of the retina. Piezo1 mRNA was more abundantly expressed than Piezo2 mRNA in these ocular tissues. Piezo1 but not Piezo2 mRNA was detected in the inner nuclear layer and outer nuclear layer of the retina. Our results suggest that PIEZO1 and PIEZO2 are differentially expressed and may have distinct roles as mechanosensors in glaucoma-relevant ocular tissues.

Liquid-biopsy proteomics combined with AI identifies cellular drivers of eye aging and disease in vivo

Single-cell analysis in living humans is essential for understanding disease mechanisms, but it is impractical in non-regenerative organs, such as the eye and brain, because tissue biopsies would cause serious damage. We resolve this problem by integrating proteomics of liquid biopsies with single-cell transcriptomics from all known ocular cell types to trace the cellular origin of 5,953 proteins detected in the aqueous humor. We identified hundreds of cell-specific protein markers, including for individual retinal cell types. Surprisingly, our results reveal that retinal degeneration occurs in Parkinson's disease, and the cells driving diabetic retinopathy switch with disease stage. Finally, we developed artificial intelligence (AI) models to assess individual cellular aging and found that many eye diseases not associated with chronological age undergo accelerated molecular aging of disease-specific cell types. Our approach, which can be applied to other organ systems, has the potential to transform molecular diagnostics and prognostics while uncovering new cellular disease and aging mechanisms.

Role of integrins in the development of fibrosis in the trabecular meshwork

Primary open angle glaucoma (POAG) is a progressive and chronic disease exhibiting many of the features of fibrosis. The extracellular matrix (ECM) in the trabecular meshwork (TM) undergoes extensive remodeling and enhanced rigidity, resembling fibrotic changes. In addition, there are changes associated with myofibroblast activation and cell contractility that further drives tissue fibrosis and stiffening. This review discusses what is known about the integrins in the TM and their involvement in fibrotic processes.

miRNA Studies in Glaucoma: A Comprehensive Review of Current Knowledge and Future Perspectives

Glaucoma, a neurodegenerative disorder that leads to irreversible blindness, remains a challenge because of its complex nature. MicroRNAs (miRNAs) are crucial regulators of gene expression and are associated with glaucoma and other diseases. We aimed to review and discuss the advantages and disadvantages of miRNA-focused molecular studies in glaucoma through discussing their potential as biomarkers for early detection and diagnosis; offering insights into molecular pathways and mechanisms; and discussing their potential utility with respect to personalized medicine, their therapeutic potential, and non-invasive monitoring. Limitations, such as variability, small sample sizes, sample specificity, and limited accessibility to ocular tissues, are also addressed, underscoring the need for robust protocols and collaboration. Reproducibility and validation are crucial to establish the credibility of miRNA research findings, and the integration of bioinformatics tools for miRNA database creation is a valuable component of a comprehensive approach to investigate miRNA aberrations in patients with glaucoma. Overall, miRNA research in glaucoma has provided significant insights into the molecular mechanisms of the disease, offering potential biomarkers, diagnostic tools, and therapeutic targets. However, addressing challenges such as variability and limited tissue accessibility is essential, and further investigations and validation will contribute to a deeper understanding of the functional significance of miRNAs in glaucoma.

Towards Universal Cell Embeddings: Integrating Single-cell RNA-seq Datasets across Species with SATURN

Analysis of single-cell datasets generated from diverse organisms offers unprecedented opportunities to unravel fundamental evolutionary processes of conservation and diversification of cell types. However, inter-species genomic differences limit the joint analysis of cross-species datasets to homologous genes. Here, we present SATURN, a deep learning method for learning universal cell embeddings that encodes genes' biological properties using protein language models. By coupling protein embeddings from language models with RNA expression, SATURN integrates datasets profiled from different species regardless of their genomic similarity. SATURN has a unique ability to detect functionally related genes co-expressed across species, redefining differential expression for cross-species analysis. We apply SATURN to three species whole-organism atlases and frog and zebrafish embryogenesis datasets. We show that cell embeddings learnt in SATURN can be effectively used to transfer annotations across species and identify both homologous and species-specific cell types, even across evolutionarily remote species. Finally, we use SATURN to reannotate the five species Cell Atlas of Human Trabecular Meshwork and Aqueous Outflow Structures and find evidence of potentially divergent functions between glaucoma associated genes in humans and other species.

Single-cell transcriptome sequencing reveals tumor heterogeneity in family neuroblastoma

Neuroblastoma(NB) is the most common extracranial solid tumor in childhood, and it is now believed that some patients with NB have an underlying genetic susceptibility, which may be one of the reasons for the multiplicity of NB patients within a family line. Even within the same family, the samples show great variation and can present as ganglioneuroblastoma or even benign ganglioneuroma. The genomics of NB is still unclear and more in-depth studies are needed to reveal its key components. We first performed single-cell RNA sequencing(sc-RNAseq) analysis on clinical specimens of two family neuroblastoma(FNB) and four sporadic NB cases. A complete transcriptional profile of FNB was constructed from 18,394 cells from FNB, and we found that SDHD may be genetically associated with FNB and identified a prognostic related CAF subtype in FNB: Fib-4. Single-cell flux estimation analysis (scFEA) results showed that malignant cells were associated with arginine spermine, oxaloacetate and hypoxanthine, and that malignant cells metabolize lactate at lower levels than T cells. Our study provides new resources and ideas for the development of the genomics of family NB, and the mechanisms of cell-to-cell interactions and communication and the metabolic landscape will provide new therapeutic targets.

FYN regulates aqueous humor outflow and IOP through the phosphorylation of VE-cadherin

The exact sites and molecules that determine resistance to aqueous humor drainage and control intraocular pressure (IOP) need further elaboration. Proposed sites include the inner wall of Schlemms's canal and the juxtacanalicular trabecular meshwork ocular drainage tissues. The adherens junctions (AJs) of Schlemm's canal endothelial cells (SECs) must both preserve the blood-aqueous humor (AQH) barrier and be conducive to AQH drainage. How homeostatic control of AJ permeability in SC occurs and how such control impacts IOP is unclear. We hypothesized that mechano-responsive phosphorylation of the junctional molecule VE-CADHERIN (VEC) by SRC family kinases (SFKs) regulates the permeability of SEC AJs. We tested this by clamping IOP at either 16 mmHg, 25 mmHg, or 45 mmHg in mice and then measuring AJ permeability and VEC phosphorylation. We found that with increasing IOP: 1) SEC AJ permeability increased, 2) VEC phosphorylation was increased at tyrosine-658, and 3) SFKs were activated at the AJ. Among the two SFKs known to phosphorylate VEC, FYN, but not SRC, localizes to the SC. Furthermore, FYN mutant mice had decreased phosphorylation of VEC at SEC AJs, dysregulated IOP, and reduced AQH outflow. Together, our data demonstrate that increased IOP activates FYN in the inner wall of SC, leading to increased phosphorylation of AJ VEC and, thus, decreased resistance to AQH outflow. These findings support a crucial role of mechanotransduction signaling in IOP homeostasis within SC in response to IOP. These data strongly suggest that the inner wall of SC partially contributes to outflow resistance.

Artificial Intelligence Models for Cell Type and Subtype Identification Based on Single-Cell RNA Sequencing Data in Vision Science

Single-cell RNA sequencing (scRNA-seq) provides a high throughput, quantitative and unbiased framework for scientists in many research fields to identify and characterize cell types within heterogeneous cell populations from various tissues. However, scRNA-seq based identification of discrete cell-types is still labor intensive and depends on prior molecular knowledge. Artificial intelligence has provided faster, more accurate, and user-friendly approaches for cell-type identification. In this review, we discuss recent advances in cell-type identification methods using artificial intelligence techniques based on single-cell and single-nucleus RNA sequencing data in vision science. The main purpose of this review paper is to assist vision scientists not only to select suitable datasets for their problems, but also to be aware of the appropriate computational tools to perform their analysis. Developing novel methods for scRNA-seq data analysis remains to be addressed in future studies.

PLAE Web App Enables Powerful Searching and Multiple Visualizations Across One Million Unified Single-Cell Ocular Transcriptomes

Purpose

To create a high-performance reactive web application to query single-cell gene expression data across cell type, species, study, and other factors.

Methods

We updated the content and structure of the underlying data (single cell Eye in a Disk [scEiaD]) and wrote the web application PLAE (https://plae.nei.nih.gov) to visualize and explore the data.

Results

The new portal provides quick visualization of over a million individual cells from vertebrate eye and body transcriptomes encompassing four species, 60 cell types, six ocular tissues, and 23 body tissues across 35 publications. To demonstrate the value of this unified pan-eye dataset, we replicated known neurogenic and cone macula markers in addition to proposing six new cone human region markers.

Conclusions

The PLAE web application offers the eye community a powerful and quick means to test hypotheses related to gene expression across a highly diverse, community-derived database.

Translational Relevance

The PLAE resource enables any researcher or clinician to study and research gene expression patterning across a wide variety of curated ocular cell types with a responsive web app.

Characterization of extracellular matrix deposited by segmental trabecular meshwork cells

PURPOSE

Biophysical and biochemical attributes of the extracellular matrix are major determinants of cell fate in homeostasis and disease. Ocular hypertension and glaucoma are diseases where the trabecular meshwork tissue responsible for aqueous humor egress becomes stiffer accompanied by changes in its matrisome in a segmental manner with regions of high or low flow. Prior studies demonstrate these alterations in the matrix are dynamic in response to age and pressure changes. The underlying reason for segmentation or differential response to pressure and stiffening are unknown. This is largely due to a lack of appropriate models (in vitro or ex vivo) to study this phenomena.

METHODS

Primary trabecular meshwork cells were isolated from segmental flow regions, and cells were cultured for 4 weeks in the presence or absence or dexamethasone to obtain cell derived matrices (CDM). The biomechanical attributes of the CDM, composition of the matrisome, and incidence of crosslinks were determined by atomic force microscopy and mass spectrometry.

RESULTS

Data demonstrate that matrix deposited by cells from low flow regions are stiffer and exhibit a greater number of immature and mature crosslinks, and that these are exacerbated in the presence of steroid. We also show a differential response of high or low flow cells to steroid via changes observed in the matrix composition. However, no correlations were observed between elastic moduli and presence or absence of mature and immature crosslinks in the CDMs.

CONCLUSION

Regardless of a direct correlation between matrix stiffness and crosslinks, we observed distinct differences in the composition and mechanics of the matrices deposited by segmental flow cells. These results suggest distinct differences in cellular identify and likely a basis for mechanical memory post isolation and culture. Nevertheless, we conclude that although a mechanistic basis for matrix stiffness was undetermined in this study, it is a viable tool to study cell-matrix interactions and further our understanding of trabecular meshwork pathobiology.

Differential roles of FOXC2 in the trabecular meshwork and Schlemm's canal in glaucomatous pathology

Impaired development and maintenance of Schlemm's canal (SC) are associated with perturbed aqueous humor outflow and intraocular pressure. The angiopoietin (ANGPT)/TIE2 signaling pathway regulates SC development and maintenance, whereas the molecular mechanisms of crosstalk between SC and the neural crest (NC)-derived neighboring tissue, the trabecular meshwork (TM), are poorly understood. Here, we show NC-specific forkhead box (Fox)c2 deletion in mice results in impaired SC morphogenesis, loss of SC identity, and elevated intraocular pressure. Visible-light optical coherence tomography analysis further demonstrated functional impairment of the SC in response to changes in intraocular pressure in NC-Foxc2 -/- mice, suggesting altered TM biomechanics. Single-cell RNA-sequencing analysis identified that this phenotype is predominately characterized by transcriptional changes associated with extracellular matrix organization and stiffness in TM cell clusters, including increased matrix metalloproteinase expression, which can cleave the TIE2 ectodomain to produce soluble TIE2. Moreover, endothelial-specific Foxc2 deletion impaired SC morphogenesis because of reduced TIE2 expression, which was rescued by deleting the TIE2 phosphatase VE-PTP. Thus, Foxc2 is critical in maintaining SC identity and morphogenesis via TM-SC crosstalk.

Posterior corneoscleral limbus: Architecture, stem cells, and clinical implications

The limbus is a transition from the cornea to conjunctiva and sclera. In human eyes, this thin strip has a rich variation of tissue structures and composition, typifying a change from scleral irregularity and opacity to corneal regularity and transparency; a variation from richly vascularized conjunctiva and sclera to avascular cornea; the neural passage and drainage of aqueous humor. The limbal stroma is enriched with circular fibres running parallel to the corneal circumference, giving its unique role in absorbing small pressure changes to maintain corneal curvature and refractivity. It contains specific niches housing different types of stem cells for the corneal epithelium, stromal keratocytes, corneal endothelium, and trabecular meshwork. This truly reflects the important roles of the limbus in ocular physiology, and the limbal functionality is crucial for corneal health and the entire visual system. Since the anterior limbus containing epithelial structures and limbal epithelial stem cells has been extensively reviewed, this article is focused on the posterior limbus. We have discussed the structural organization and cellular components of the region beneath the limbal epithelium, the characteristics of stem cell types: namely corneal stromal stem cells, endothelial progenitors and trabecular meshwork stem cells, and recent advances leading to the emergence of potential cell therapy options to replenish their respective mature cell types and to correct defects causing corneal abnormalities. We have reviewed different clinical disorders associated with defects of the posterior limbus and summarized the available preclinical and clinical evidence about the developing topic of cell-based therapy for corneal disorders.

The Association between Vascular Abnormalities and Glaucoma-What Comes First?

Glaucoma is a leading cause of irreversible blindness worldwide. While intraocular pressure (IOP) presents a major risk factor, the underlying pathophysiology still remains largely unclear. The correlation between vascular abnormalities and glaucoma has been deliberated for decades. Evidence for a role played by vascular factors in the pathogenesis of glaucomatous neurodegeneration has already been postulated. In addition, the fact that glaucoma causes both structural and functional changes to retinal blood vessels has been described. This review aims to investigate the published evidence concerning the relationship between vascular abnormalities and glaucoma, and to provide an overview of the "chicken or egg" dilemma in glaucoma. In this study, several biomarkers of glaucoma progression from a vascular perspective, including endothelin-1 (ET-1), nitric oxide, vascular endothelial growth factor (VEGF), and matrix metalloproteinases (MMPs), were identified and subsequently assessed for their potential as pharmacological intervention targets.

Transcriptomic analysis of the ocular posterior segment completes a cell atlas of the human eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here, we used high-throughput single-nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in six extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases-for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ~151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the data presented here complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as the neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Evaluation of trabecular meshwork-specific promoters in vitro and in vivo using scAAV2 vectors expressing C3 transferase

AIM

To evaluate the potential of two trabecular meshwork (TM)-specific promoters, Chitinase 3-like 1 (Ch3L1) and matrix gla protein (MGP), for improving specificity and safety in glaucoma gene therapy based on self-complementary AAV2 (scAAV2) vector technologies.

METHODS

An scAAV2 vector with C3 transferase (C3) as the reporter gene (scAAV2-C3) was selected. The scAAV2-C3 vectors were driven by Ch3L1 (scAAV2-Ch3L1-C3), MGP (scAAV2-MGP-C3), enhanced MGP (scAAV2-eMGP-C3) and cytomegalovirus (scAAV2-CMV-C3), respectively. The cultured primary human TM cells were treated with each vector at different multiplicities of infections. Changes in cell morphology were observed by phase contrast microscopy. Actin stress fibers and Rho GTPases/Rho-associated protein kinase pathway-related molecules were assessed by immunofluorescence staining, real-time quantitative polymerase chain reaction and Western blot. Each vector was injected intracamerally into the one eye of each rat at low and high doses respectively. In vivo green fluorescence was visualized by a Micron III Retinal Imaging Microscope. Intraocular pressure (IOP) was monitored using a rebound tonometer. Ocular responses were evaluated by slit-lamp microscopy. Ocular histopathology analysis was examined by hematoxylin and eosin staining.

RESULTS

In TM cell culture studies, the vector-mediated C3 expression induced morphologic changes, disruption of actin cytoskeleton and reduction of fibronectin expression in TM cells by inhibiting the Rho GTPases/Rho-associated protein kinase signaling pathway. At the same dose, these changes were significant in TM cells treated with scAAV2-CMV-C3 or scAAV2-Ch3L1-C3, but not in cells treated with scAAV2-eMGP-C3 or scAAV2-MGP-C3. At low-injected dose, the IOP was significantly decreased in the scAAV2-Ch3L1-C3-injected eyes but not in scAAV2-MGP-C3-injected and scAAV2-eMGP-C3-injected eyes. At high-injected dose, significant IOP reduction was observed in the scAAV2-eMGP-C3-injected eyes but not in scAAV2-MGP-C3-injected eyes. Similar to scAAV2-CMV-C3, scAAV2-Ch3L1-C3 vector showed efficient transduction both in the TM and corneal endothelium. In anterior segment tissues of scAAV2-eMGP-C3-injected eyes, no obvious morphological changes were found except for the TM. Inflammation was absent.

CONCLUSION

In scAAV2-transduced TM cells, the promoter-driven efficiency of Ch3L1 is close to that of cytomegalovirus, but obviously higher than that of MGP. In the anterior chamber of rat eye, the transgene expression pattern of scAAV2 vector is presumably affected by MGP promoter, but not by Ch3L1 promoter. These findings would provide a useful reference for improvement of specificity and safety in glaucoma gene therapy using scAAV2 vector.

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.

Rare protective variants and glaucoma-relevant cell stressors modulate Angiopoietin-like 7 expression

Rare missense and nonsense variants in the Angiopoietin-like 7 (ANGPTL7) gene confer protection from primary open-angle glaucoma (POAG), though the functional mechanism remains uncharacterized. Interestingly, a larger variant effect size strongly correlates with in silico predictions of increased protein instability (r = -0.98), suggesting that protective variants lower ANGPTL7 protein levels. Here, we show that missense and nonsense variants cause aggregation of mutant ANGPTL7 protein in the endoplasmic reticulum (ER) and decreased levels of secreted protein in human trabecular meshwork (TM) cells; a lower secreted:intracellular protein ratio strongly correlates with variant effects on intraocular pressure (r = 0.81). Importantly, accumulation of mutant protein in the ER does not increase expression of ER stress proteins in TM cells (P > 0.05 for all variants tested). Cyclic mechanical stress, a glaucoma-relevant physiologic stressor, also significantly lowers ANGPTL7 expression in primary cultures of human Schlemm's canal (SC) cells (-2.4-fold-change, P = 0.01). Collectively, these data suggest that the protective effects of ANGPTL7 variants in POAG stem from lower levels of secreted protein, which may modulate responses to physiologic and pathologic ocular cell stressors. Downregulation of ANGPTL7 expression may therefore serve as a viable preventative and therapeutic strategy for this common, blinding disease.

Neurofibromatosis Type 1-Associated Optic Pathway Gliomas: Current Challenges and Future Prospects

Optic pathway glioma (OPG) occurs in as many as one-fifth of individuals with the neurofibromatosis type 1 (NF1) cancer predisposition syndrome. Generally considered low-grade and slow growing, many children with NF1-OPGs remain asymptomatic. However, due to their location within the optic pathway, ~20-30% of those harboring NF1-OPGs will experience symptoms, including progressive vision loss, proptosis, diplopia, and precocious puberty. While treatment with conventional chemotherapy is largely effective at attenuating tumor growth, it is not clear whether there is much long-term recovery of visual function. Additionally, because these tumors predominantly affect young children, there are unique challenges to NF1-OPG diagnosis, monitoring, and longitudinal management. Over the past two decades, the employment of authenticated genetically engineered Nf1-OPG mouse models have provided key insights into the function of the NF1 protein, neurofibromin, as well as the molecular and cellular pathways that contribute to optic gliomagenesis. Findings from these studies have resulted in the identification of new molecular targets whose inhibition blocks murine Nf1-OPG growth in preclinical studies. Some of these promising compounds have now entered into early clinical trials. Future research focused on defining the determinants that underlie optic glioma initiation, expansion, and tumor-induced optic nerve injury will pave the way to personalized risk assessment strategies, improved tumor monitoring, and optimized treatment plans for children with NF1-OPG.

Aqueous Humor TGF-β2 and Its Association With Intraocular Pressure in a Naturally Occurring Large Animal Model of Glaucoma

Purpose

Transforming growth factor (TGF)-β2 has been widely implicated in human glaucoma pathology. The purpose of this study was to determine the source of TGF-β2 in aqueous humor (AH) and its relationship with intraocular pressure (IOP) in an inherited large animal model of glaucoma.

Methods

Sixty-six glaucomatous cats homozygous for LTBP2 mutation, and 42 normal cats were studied. IOP was measured weekly by rebound tonometry. AH was collected by anterior chamber paracentesis from each eye under general anesthesia, and serum samples collected from venous blood concurrently. Concentrations of total, active and latent TGF-β2 in AH and serum samples were measured by quantitative sandwich immunoassay. For comparisons between groups, unpaired t-test or Mann Whitney test were used, with P < 0.05 considered significant. The relationships between TGF-β2 concentrations and IOP values were examined by Pearson's correlation coefficient and generalized estimating equation.

Results

IOP and AH TGF-β2 concentrations were significantly higher in glaucomatous than in normal cats. AH TGF-β2 showed a significant, robust positive correlation with IOP in glaucomatous cats (r = 0.83, R2 = 0.70, P < 0.0001). Serum TGF-β2 did not correlate with AH TGF-β2 and was not significantly different between groups. TGF-β2 mRNA and protein expression were significantly increased in local ocular tissues in glaucomatous cats.

Conclusions

Enhanced, local ocular production of TGF-β2 with a robust positive association with IOP was identified in this spontaneous feline glaucoma model, providing a foundation for preclinical testing of novel therapeutics to limit disease-associated AH TGF-β2 elevation and signaling in glaucoma.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

A new gene-scoring method for uncovering novel glaucoma-related genes using non-negative matrix factorization based on RNA-seq data

Early diagnosis and treatment of glaucoma are challenging. The discovery of glaucoma biomarkers based on gene expression data could potentially provide new insights for early diagnosis, monitoring, and treatment options of glaucoma. Non-negative Matrix Factorization (NMF) has been widely used in numerous transcriptome data analyses in order to identify subtypes and biomarkers of different diseases; however, its application in glaucoma biomarker discovery has not been previously reported. Our study applied NMF to extract latent representations of RNA-seq data from BXD mouse strains and sorted the genes based on a novel gene scoring method. The enrichment ratio of the glaucoma-reference genes, extracted from multiple relevant resources, was compared using both the classical differentially expressed gene (DEG) analysis and NMF methods. The complete pipeline was validated using an independent RNA-seq dataset. Findings showed our NMF method significantly improved the enrichment detection of glaucoma genes. The application of NMF with the scoring method showed great promise in the identification of marker genes for glaucoma.

Degradomic Identification of Membrane Type 1-Matrix Metalloproteinase as an ADAMTS9 and ADAMTS20 Substrate

The secreted metalloproteases ADAMTS9 and ADAMTS20 are implicated in extracellular matrix proteolysis and primary cilium biogenesis. Here, we show that clonal gene-edited RPE-1 cells in which ADAMTS9 was inactivated, and which constitutively lack ADAMTS20 expression, have morphologic characteristics distinct from parental RPE-1 cells. To investigate underlying proteolytic mechanisms, a quantitative terminomics method, terminal amine isotopic labeling of substrates was used to compare the parental and gene-edited RPE-1 cells and their medium to identify ADAMTS9 substrates. Among differentially abundant neo-amino (N) terminal peptides arising from secreted and transmembrane proteins, a peptide with lower abundance in the medium of gene-edited cells suggested cleavage at the Tyr314-Gly315 bond in the ectodomain of the transmembrane metalloprotease membrane type 1-matrix metalloproteinase (MT1-MMP), whose mRNA was also reduced in gene-edited cells. This cleavage, occurring in the MT1-MMP hinge, that is, between the catalytic and hemopexin domains, was orthogonally validated both by lack of an MT1-MMP catalytic domain fragment in the medium of gene-edited cells and restoration of its release from the cell surface by reexpression of ADAMTS9 and ADAMTS20 and was dependent on hinge O-glycosylation. A C-terminally semitryptic MT1-MMP peptide with greater abundance in WT RPE-1 medium identified a second ADAMTS9 cleavage site in the MT1-MMP hemopexin domain. Consistent with greater retention of MT1-MMP on the surface of gene-edited cells, pro-MMP2 activation, which requires cell surface MT1-MMP, was increased. MT1-MMP knockdown in gene-edited ADAMTS9/20-deficient cells restored focal adhesions but not ciliogenesis. The findings expand the web of interacting proteases at the cell surface, suggest a role for ADAMTS9 and ADAMTS20 in regulating cell surface activity of MT1-MMP, and indicate that MT1-MMP shedding does not underlie their observed requirement in ciliogenesis.

Recently Approved Drugs for Lowering and Controlling Intraocular Pressure to Reduce Vision Loss in Ocular Hypertensive and Glaucoma Patients

Serious vision loss occurs in patients affected by chronically raised intraocular pressure (IOP), a characteristic of many forms of glaucoma where damage to the optic nerve components causes progressive degeneration of retinal and brain neurons involved in visual perception. While many risk factors abound and have been validated for this glaucomatous optic neuropathy (GON), the major one is ocular hypertension (OHT), which results from the accumulation of excess aqueous humor (AQH) fluid in the anterior chamber of the eye. Millions around the world suffer from this asymptomatic and progressive degenerative eye disease. Since clinical evidence has revealed a strong correlation between the reduction in elevated IOP/OHT and GON progression, many drugs, devices, and surgical techniques have been developed to lower and control IOP. The constant quest for new pharmaceuticals and other modalities with superior therapeutic indices has recently yielded health authority-approved novel drugs with unique pharmacological signatures and mechanism(s) of action and AQH drainage microdevices for effectively and durably treating OHT. A unique nitric oxide-donating conjugate of latanoprost, an FP-receptor prostaglandin (PG; latanoprostene bunod), new rho kinase inhibitors (ripasudil; netarsudil), a novel non-PG EP2-receptor-selective agonist (omidenepag isopropyl), and a form of FP-receptor PG in a slow-release intracameral implant (Durysta) represent the additions to the pharmaceutical toolchest to mitigate the ravages of OHT. Despite these advances, early diagnosis of OHT and glaucoma still lags behind and would benefit from further concerted effort and attention.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

A temporal single cell transcriptome atlas of zebrafish anterior segment development

Anterior segment dysgenesis (ASD), resulting in vision impairment, stems from maldevelopment of anterior segment (AS) tissues. Incidence of ASD has been linked to malfunction of periocular mesenchyme cells (POM). POM cells specify into anterior segment mesenchyme (ASM) cells which colonize and produce AS tissues. In this study we uncover ASM developmental trajectories associated with formation of the AS. Using a transgenic line of zebrafish that fluorescently labels the ASM throughout development, Tg[foxc1b:GFP], we isolated GFP+ ASM cells at several developmental timepoints (48-144 hpf) and performed single cell RNA sequencing. Clustering analysis indicates subdifferentiation of ASM as early as 48 hpf and subsequent diversification into corneal epithelium/endothelium/stroma, or annular ligament (AL) lineages. Tracking individual clusters reveals common developmental pathways, up to 72 hpf, for the AL and corneal endothelium/stroma and distinct pathways for corneal epithelium starting at 48 hpf. Spatiotemporal validation of over 80 genes found associated with AS development demonstrates a high degree of conservation with mammalian trabecular meshwork and corneal tissues. In addition, we characterize thirteen novel genes associated with annular ligament and seven with corneal development. Overall, the data provide a molecular verification of the long-standing hypothesis that POM derived ASM give rise to AS tissues and highlight the high degree of conservation between zebrafish and mammals.

Application of a Magnetic Platform in α6 Integrin-Positive iPSC-TM Purification

The emergence of induced pluripotent stem cell (iPSC) technology has provided a new approach to regenerating decellularized trabecular meshwork (TM) in glaucoma. We have previously generated iPSC-derived TM (iPSC-TM) using a medium conditioned by TM cells and verified its function in tissue regeneration. Because of the heterogeneity of iPSCs and the isolated TM cells, iPSC-TM cells appear to be heterogeneous, which impedes our understanding of how the decellularized TM may be regenerated. Herein, we developed a protocol based on a magnetic-activated cell sorting (MACS) system or an immunopanning (IP) method for sorting integrin subunit alpha 6 (ITGA6)-positive iPSC-TM, an example of the iPSC-TM subpopulation. We first analyzed the purification efficiency of these two approaches by flow cytometry. In addition, we also determined cell viability by analyzing the morphologies of the purified cells. To conclude, the MACS-based purification could yield a higher ratio of ITGA6-positive iPSC-TM and maintain a relatively higher cell viability than the IP-based method, allowing for the preparation of any iPSC-TM subpopulation of interest and facilitating a better understanding of the regenerative mechanism of iPSC-based therapy.

Genome-wide analysis of genetic pleiotropy and causal genes across three age-related ocular disorders

Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.

Characterization of extracellular matrix deposited by segmental trabecular meshwork cells

Biophysical and biochemical attributes of the extracellular matrix are major determinants of cell fate in homeostasis and disease. Ocular hypertension and glaucoma are diseases where the trabecular meshwork tissue responsible for aqueous humor egress becomes stiffer accompanied by changes in its matrisome in a segmental manner with regions of high or low flow. Prior studies demonstrate these alterations in the matrix are dynamic in response to age and pressure changes. The underlying reason for segmentation or differential response to pressure and stiffening are unknown. This is largely due to a lack of appropriate models ( in vitro or ex vivo ) to study this phenomena. In this study, we characterize the biomechanical attributes, matrisome, and incidence of crosslinks in the matrix deposited by primary cells isolated from segmental flow regions and when treated with glucocorticosteroid. Data demonstrate that matrix deposited by cells from low flow regions are stiffer and exhibit a greater number of immature and mature crosslinks, and that these are exacerbated in the presence of steroid. We also show a differential response of high or low flow cells to steroid via changes observed in the matrix composition. We conclude that although a mechanistic basis for matrix stiffness was undetermined in this study, it is a viable tool to study cell-matrix interactions and further our understanding of trabecular meshwork pathobiology.