Molecular characterization of foveal versus peripheral human retina by single-cell RNA sequencing

Differential gene expression between central and peripheral retinal regions in dogs and comparison with humans

The dog retina contains a central macula-like region, and there are reports of central retinal disorders in dogs with shared genetic etiologies with humans. Defining central/peripheral gene expression profiles may provide insight into the suitability of dogs as models for human disorders. We determined central/peripheral posterior eye gene expression profiles in dogs and interrogated inherited retinal and macular disease-associated genes for differential expression between central and peripheral regions. Bulk tissue RNA sequencing was performed on 8 mm samples of the dog central and superior peripheral regions, sampling retina and retinal pigmented epithelium/choroid separately. Reads were mapped to CanFam3.1, read counts were analyzed to determine significantly differentially expressed genes (DEGs). A similar analytic pipeline was used with a published bulk-tissue RNA sequencing human dataset. Pathways and processes involved in significantly DEGs were identified (Database for Annotation, Visualization and Integrated Discovery). Dogs and humans shared the extent and direction of central retinal differential gene expression, with multiple shared biological pathways implicated in differential expression. Many genes implicated in heritable retinal disorders in dogs and humans were differentially expressed between central and periphery. Approximately half of genes associated with human age-related macular degeneration were differentially expressed in human and dog tissues. We have identified similarities and differences in central/peripheral gene expression profiles between dogs and humans which can be applied to further define the relevance of dogs as models for human retinal disorders.

Cross-species single-cell landscapes identify the pathogenic gene characteristics of inherited retinal diseases

Introduction

Inherited retinal diseases (IRDs) affect ∼4.5 million people worldwide. Elusive pathogenic variants in over 280 genes are associated with one or more clinical forms of IRDs. It is necessary to understand the complex interaction among retinal cell types and pathogenic genes by constructing a regulatory network. In this study, we attempt to establish a panoramic expression view of the cooperative work in retinal cells to understand the clinical manifestations and pathogenic bases underlying IRDs.

Methods

Single-cell RNA sequencing (scRNA-seq) data on the retinas from 35 retina samples of 3 species (human, mouse, and zebrafish) including 259,087 cells were adopted to perform a comparative analysis across species. Bioinformatic tools were used to conduct weighted gene co-expression network analysis (WGCNA), single-cell regulatory network analysis, cell-cell communication analysis, and trajectory inference analysis.

Results

The cross-species comparison revealed shared or species-specific gene expression patterns at single-cell resolution, such as the stathmin family genes, which were highly expressed specifically in zebrafish Müller glias (MGs). Thirteen gene modules were identified, of which nine were associated with retinal cell types, and Gene Ontology (GO) enrichment of module genes was consistent with cell-specific highly expressed genes. Many IRD genes were identified as hub genes and cell-specific regulons. Most IRDs, especially the retinitis pigmentosa (RP) genes, were enriched in rod-specific regulons. Integrated expression and transcription regulatory network genes, such as congenital stationary night blindness (CSNB) genes GRK1, PDE6B, and TRPM1, showed cell-specific expression and transcription characteristics in either rods or bipolar cells (BCs). IRD genes showed evolutionary conservation (GNAT2, PDE6G, and SAG) and divergence (GNAT2, MT-ND4, and PDE6A) along the trajectory of photoreceptors (PRs) among species. In particular, the Leber congenital amaurosis (LCA) gene OTX2 showed high expression at the beginning of the trajectory of both PRs and BCs.

Conclusion

We identified molecular pathways and cell types closely connected with IRDs, bridging the gap between gene expression, genetics, and pathogenesis. The IRD genes enriched in cell-specific modules and regulons suggest that these diseases share common etiological bases. Overall, mining of interspecies transcriptome data reveals conserved transcriptomic features of retinas across species and promising applications in both normal retina anatomy and retina pathology.

Single cell transcriptomics reveals early photoreceptor states, cell-specific transcript isoforms, and cancer-predisposing features

Human cone photoreceptors differ from rods and serve as the retinoblastoma cell-of-origin. Here, we used deep full-length single-cell RNA-sequencing to distinguish post-mitotic cone and rod developmental states and cone-specific features that contribute to retinoblastomagenesis. The analyses revealed early post-mitotic cone- and rod-directed populations characterized by higher THRB or NRL regulon activities, an immature photoreceptor precursor population with concurrent cone and rod gene and regulon expression, and distinct early and late cone and rod maturation states distinguished by maturation-associated declines in RAX regulon activity. Unexpectedly, both L/M cone and rod precursors co-expressed NRL and THRB RNAs, yet they differentially expressed functionally antagonistic NRL isoforms and prematurely terminated THRB transcripts. Early L/M cone precursors exhibited successive expression of lncRNAs along with MYCN, which composed the seventh most L/M-cone-specific regulon, and SYK, which contributed to the early cone precursors' proliferative response to RB1 loss. These findings reveal previously unrecognized photoreceptor precursor states and a role for early cone-precursor-intrinsic SYK expression in retinoblastoma initiation.

VDAC in Retinal Health and Disease

The retina, a tissue of the central nervous system, is vital for vision as its photoreceptors capture light and transform it into electrical signals, which are further processed before they are sent to the brain to be interpreted as images. The retina is unique in that it is continuously exposed to light and has the highest metabolic rate and demand for energy amongst all the tissues in the body. Consequently, the retina is very susceptible to oxidative stress. VDAC, a pore in the outer membrane of mitochondria, shuttles metabolites between mitochondria and the cytosol and normally protects cells from oxidative damage, but when a cell's integrity is greatly compromised it initiates cell death. There are three isoforms of VDAC, and existing evidence indicates that all three are expressed in the retina. However, their precise localization and function in each cell type is unknown. It appears that most retinal cells express substantial amounts of VDAC2 and VDAC3, presumably to protect them from oxidative stress. Photoreceptors express VDAC2, HK2, and PKM2-key proteins in the Warburg pathway that also protect these cells. Consistent with its role in initiating cell death, VDAC is overexpressed in the retinal degenerative diseases retinitis pigmentosa, age related macular degeneration (AMD), and glaucoma. Treatment with antioxidants or inhibiting VDAC oligomerization reduced its expression and improved cell survival. Thus, VDAC may be a promising therapeutic candidate for the treatment of these diseases.

Shotgun lipidomics of human subretinal fluids under rod-dominant retina reveals cone-dominated lipids

Subretinal fluid (SRF) accumulates between photoreceptor outer segments and retinal pigment epithelium during rhegmatogenous retinal detachment. Biomolecular components such as lipids originate from cells surrounding the SRF. Knowledge of the composition of these molecules in SRF potentially provides mechanistic insight into the physiologic transfer of lipids between retinal tissue compartments. Using mass spectrometry and tandem mass spectrometry analysis on an electrospray ionization quadrupole-time-of-flight mass spectrometer, we identified a total of 115 lipid molecular species of 11 subclasses and 9 classes in two samples from two patients with rhegmatogenous retinal detachment. These included 47 glycerophosphocholines, 6 glycerophosphoethanolamines, 1 glycerophosphoinositol, 18 sphingomyelins, 9 cholesteryl esters, free cholesterol, 3 ceramides, 22 triacylglycerols and 8 free fatty acids. Glycerophosphocholines were of the highest intensity. By minimizing the formation of different adduct forms or clustering ions of different adducts, we determined the relative intensity of lipid molecular species within the same subclasses. The profiles were compared with those of retinal cells available in the published literature. The glycerophosphocholine profile of SRF was similar to that of cone outer segments, suggesting that outer segment degradation products are constitutively released into the interphotoreceptor matrix, appearing in SRF during detachment. This hypothesis was supported by the retinal distributions of corresponding lipid synthases' mRNA expression obtained from an online resource based on publicly available single-cell sequencing data. In contrast, based on lipid profiles and relevant gene expression in this study, the sources of free cholesterol and cholesteryl esters in SRF appeared more ambiguous, possibly reflecting that outer retina takes up plasma lipoproteins. Further studies to identify and quantify lipids in SRF will help better understand etiology of diseases relevant to outer retina.

Recent Achievements in the Heterogeneity of Mammalian and Human Retinal Pigment Epithelium: In Search of a Stem Cell

Retinal pigment epithelium (RPE) cells are important fundamentally for the development and function of the retina. In this regard, the study of the morphological and molecular properties of RPE cells, as well as their regenerative capabilities, is of particular importance for biomedicine. However, these studies are complicated by the fact that, despite the external morphological similarity of RPE cells, the RPE is a population of heterogeneous cells, the molecular genetic properties of which have begun to be revealed by sequencing methods only in recent years. This review carries out an analysis of the data from morphological and molecular genetic studies of the heterogeneity of RPE cells in mammals and humans, which reveals the individual differences in the subpopulations of RPE cells and the possible specificity of their functions. Particular attention is paid to discussing the properties of "stemness," proliferation, and plasticity in the RPE, which may be useful for uncovering the mechanisms of retinal diseases associated with pathologies of the RPE and finding new ways of treating them.

Expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF in the developing rhesus monkey retina

The fovea centralis (fovea) is a specialized region of the primate retina that plays crucial roles in high-resolution visual acuity and color perception. The fovea is characterized by a high density of cone photoreceptors and no rods, and unique anatomical properties that contribute to its remarkable visual capabilities. Early histological analyses identified some of the key events that contribute to foveal development, but the mechanisms that direct the specification of this area are not understood. Recently, the expression of the retinoic acid-metabolizing enzyme CYP26A1 has become a hallmark of some of the retinal specializations found in vertebrates, including the primate fovea and the high-acuity area in avian species. In chickens, the retinoic acid pathway regulates the expression of FGF8 to then direct the development of a rod-free area. Similarly, high levels of CYP26A1, CDKN1A, and NPVF expression have been observed in the primate macula using transcriptomic approaches. However, which retinal cells express these genes and their expression dynamics in the developing primate eye remain unknown. Here, we systematically characterize the expression patterns of CYP26A1, FGF8, CDKN1A, and NPVF during the development of the rhesus monkey retina, from early stages of development in the first trimester until the third trimester (near term). Our data suggest that some of the markers previously proposed to be fovea-specific are not enriched in the progenitors of the rhesus monkey fovea. In contrast, CYP26A1 is expressed at high levels in the progenitors of the fovea, while it localizes in a subpopulation of macular Müller glia cells later in development. Together these data provide invaluable insights into the expression dynamics of several molecules in the nonhuman primate retina and highlight the developmental advancement of the foveal region.

Dynamic changes in ocular shape during human development and its implications for retina fovea formation

The human fovea is known for its distinctive pit-like appearance, which results from the displacement of retinal layers superficial to the photoreceptors cells. The photoreceptors are found at high density within the foveal region but not the surrounding retina. Efforts to elucidate the mechanisms responsible for these unique features have ruled out cell death as an explanation for pit formation and changes in cell proliferation as the cause of increased photoreceptor density. These findings have led to speculation that mechanical forces acting within and on the retina during development underly the formation of foveal architecture. Here we review eye morphogenesis and retinal remodeling in human embryonic development. Our meta-analysis of the literature suggests that fovea formation is a protracted process involving dynamic changes in ocular shape that start early and continue throughout most of human embryonic development. From these observations, we propose a new model for fovea development.

The Sleep Quality- and Myopia-Linked PDE11A-Y727C Variant Impacts Neural Physiology by Reducing Catalytic Activity and Altering Subcellular Compartmentalization of the Enzyme

Recently, a Y727C variant in the dual-specific 3',5'-cyclic nucleotide phosphodiesterase 11A (PDE11A-Y727C) was linked to increased sleep quality and reduced myopia risk in humans. Given the well-established role that the PDE11 substrates cAMP and cGMP play in eye physiology and sleep, we determined if (1) PDE11A protein is expressed in the retina or other eye segments in mice, (2) PDE11A-Y7272C affects catalytic activity and/or subcellular compartmentalization more so than the nearby suicide-associated PDE11A-M878V variant, and (3) Pde11a deletion alters eye growth or sleep quality in male and female mice. Western blots show distinct protein expression of PDE11A4, but not PDE11A1-3, in eyes of Pde11a WT, but not KO mice, that vary by eye segment and age. In HT22 and COS-1 cells, PDE11A4-Y727C reduces PDE11A4 catalytic activity far more than PDE11A4-M878V, with both variants reducing PDE11A4-cAMP more so than PDE11A4-cGMP activity. Despite this, Pde11a deletion does not alter age-related changes in retinal or lens thickness or axial length, nor vitreous or anterior chamber depth. Further, Pde11a deletion only minimally changes refractive error and sleep quality. That said, both variants also dramatically alter the subcellular compartmentalization of human and mouse PDE11A4, an effect occurring independently of dephosphorylating PDE11A4-S117/S124 or phosphorylating PDE11A4-S162. Rather, re-compartmentalization of PDE11A4-Y727C is due to the loss of the tyrosine changing how PDE11A4 is packaged/repackaged via the trans-Golgi network. Therefore, the protective impact of the Y727C variant may reflect a gain-of-function (e.g., PDE11A4 displacing another PDE) that warrants further investigation in the context of reversing/preventing sleep disturbances or myopia.

The sleep quality- and myopia-linked PDE11A-Y727C variant impacts neural physiology by reducing catalytic activity and altering subcellular compartmentalization of the enzyme

Recently, a Y727C variant in the dual-specific 3',5'-cyclic nucleotide phosphodiesterase 11A (PDE11A-Y727C) was linked to increased sleep quality and reduced myopia risk in humans. Given the well-established role that the PDE11 substrates cAMP and cGMP play in eye physiology and sleep, we determined if 1) PDE11A protein is expressed in the retina or other eye segments in mouse, 2) PDE11A-Y7272C affects catalytic activity and/or subcellular compartmentalization more so than the nearby suicide-associated PDE11A-M878V variant, and 3) Pde11a deletion alters eye growth or sleep quality in male and female mice. Western blots show distinct protein expression of PDE11A4, but not PDE11A1-3, in eyes of Pde11a WT-but not KO mice-that vary by eye segment and age. In HT22 and COS-1 cells, PDE11A4-Y727C reduces PDE11A4 catalytic activity far more than PDE11A4-M878V, with both variants reducing PDE11A4-cAMP more so than PDE11A4-cGMP activity. Despite this, Pde11a deletion does not alter age-related changes in retinal or lens thickness, axial length, nor vitreous or anterior chamber depth. Further, Pde11a deletion only minimally changes refractive error and sleep quality. That said, both variants also dramatically alter the subcellular compartmentalization of human and mouse PDE11A4, an effect occurring independently of dephosphorylating PDE11A4-S117/S124 or phosphorylating PDE11A4-S162. Rather, re-compartmentalization of PDE11A4-Y727C is due to the loss of the tyrosine changing how PDE11A4 is packaged/repackaged via the trans-Golgi network. Therefore, the protective impact of the Y727C variant may reflect a gain-of-function (e.g., PDE11A4 displacing another PDE) that warrants further investigation in the context of reversing/preventing sleep disturbances or myopia.

Integrated multi-omics single cell atlas of the human retina

Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here, we present an integrated multimodal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 2.4 million cells from 55 donors, including 1.4 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) of transcriptome and chromatin accessibility, unveiling over 110 types. Engaging the retina community, we annotated each cluster, refined the Cell Ontology for the retina, identified distinct marker genes, and characterized cis-regulatory elements and gene regulatory networks (GRNs) for these cell types. Our analysis uncovered intriguing differences in transcriptome, chromatin, and GRNs across cell types. In addition, we modeled changes in gene expression and chromatin openness across gender and age. This integrated atlas also enabled the fine-mapping of GWAS and eQTL variants. Accessible through interactive browsers, this multimodal cross-donor and cross-lab HRCA, can facilitate a better understanding of retinal function and pathology.

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.

Gene Expression Within a Human Choroidal Neovascular Membrane Using Spatial Transcriptomics

Purpose

Macular neovascularization is a relatively common and potentially visually devastating complication of age-related macular degeneration. In macular neovascularization, pathologic angiogenesis can originate from either the choroid or the retina, but we have limited understanding of how different cell types become dysregulated in this dynamic process.

Methods

To study how gene expression is altered in focal areas of pathology, we performed spatial RNA sequencing on a human donor eye with macular neovascularization as well as a healthy control donor. We performed differential expression to identify genes enriched within the area of macular neovascularization and used deconvolution algorithms to predict the originating cell type of these dysregulated genes.

Results

Within the area of neovascularization, endothelial cells demonstrated increased expression of genes related to Rho family GTPase signaling and integrin signaling. Likewise, VEGF and TGFB1 were identified as potential upstream regulators that could drive the observed gene expression changes produced by endothelial and retinal pigment epithelium cells in the macular neovascularization donor. These spatial gene expression profiles were compared to previous single-cell gene expression experiments in human age-related macular degeneration as well as a model of laser-induced neovascularization in mice. As a secondary aim, we investigated regional gene expression patterns within the macular neural retina and between the macular and peripheral choroid.

Conclusions

Overall, this study spatially analyzes gene expression across the retina, retinal pigment epithelium, and choroid in health and describes a set of candidate molecules that become dysregulated in macular neovascularization.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server

MOTIVATION

Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use.

RESULTS

Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis.

AVAILABILITY AND IMPLEMENTATION

The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

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.

Selective drug delivery to the retinal cells: Biological barriers and avenues

Retinal drug delivery is a challenging, but important task, because most retinal diseases are still without any proper therapy. Drug delivery to the retina is hampered by the anatomical and physiological barriers resulting in minimal bioavailability after topical ocular and systemic administrations. Intravitreal injections are current method-of-choice in retinal delivery, but these injections show short duration of action for small molecules and low target bioavailability for many protein, gene based drugs and nanomedicines. State-of-art delivery systems are based on prolonged retention, controlled drug release and physical features (e.g. size and charge). However, drug delivery to the retina is not cell-specific and these approaches do not facilitate intracellular delivery of modern biological drugs (e.g. intracellular proteins, RNA based medicines, gene editing). In this focused review we highlight biological factors and mechanisms that form the basis for the selective retinal drug delivery systems in the future. Therefore, we are presenting current knowledge related to retinal membrane transporters, receptors and targeting ligands in relation to nanomedicines, conjugates, extracellular vesicles, and melanin binding. These issues are discussed in the light of retinal structure and cell types as well as future prospects in the field. Unlike in some other fields of targeted drug delivery (e.g. cancer research), selective delivery technologies have been rarely studied, even though cell targeted delivery may be even more feasible after local administration into the eye.

Immune cell identifier and classifier (ImmunIC) for single cell transcriptomic readouts

Single cell RNA sequencing has a central role in immune profiling, identifying specific immune cells as disease markers and suggesting therapeutic target genes of immune cells. Immune cell-type annotation from single cell transcriptomics is in high demand for dissecting complex immune signatures from multicellular blood and organ samples. However, accurate cell type assignment from single-cell RNA sequencing data alone is complicated by a high level of gene expression heterogeneity. Many computational methods have been developed to respond to this challenge, but immune cell annotation accuracy is not highly desirable. We present ImmunIC, a simple and robust tool for immune cell identification and classification by combining marker genes with a machine learning method. With over two million immune cells and half-million non-immune cells from 66 single cell RNA sequencing studies, ImmunIC shows 98% accuracy in the identification of immune cells. ImmunIC outperforms existing immune cell classifiers, categorizing into ten immune cell types with 92% accuracy. We determine peripheral blood mononuclear cell compositions of severe COVID-19 cases and healthy controls using previously published single cell transcriptomic data, permitting the identification of immune cell-type specific differential pathways. Our publicly available tool can maximize the utility of single cell RNA profiling by functioning as a stand-alone bioinformatic cell sorter, advancing cell-type specific immune profiling for the discovery of disease-specific immune signatures and therapeutic targets.

Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multisystem disease. The wide range of manifestations observed in patients with mitochondrial disease results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify nonrandom distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA-Seq, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profiled the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and from healthy control donors. Utilizing the retina as a model for complex multilineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the WT allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the nonrandom nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

Cross-species scRNA-seq reveals the cellular landscape of retina and early alterations in type 2 diabetes mice

Single-cell RNA sequencing (scRNA-seq) analysis have provided an unprecedented resolution for the studies on diabetic retinopathy (DR). However, the early changes in the retina in diabetes remain unclear. A total of 8 human and mouse scRNA-seq datasets, containing 276,402 cells were analyzed individually to comprehensively delineate the retinal cell atlas. The neural retinas were isolated from the type 2 diabetes (T2D) and control mice, and scRNA-seq analysis was conducted to evaluate the early effects of diabetes on the retina. Bipolar cell (BC) heterogeneity were identified. We found some stable BCs across multiple datasets, and explored their biological functions. A new RBC subtype (Car8_RBC) in the mouse retina was validated using the multi-color immunohistochemistry. AC149090.1 was significantly upregulated in the rod cells, ON cone BCs (CBCs), OFF CBCs, and RBCs in T2D mice. Additionally, the interneurons, especially BCs, were the most vulnerable cells to diabetes by integrating scRNA-seq and genome-wide association studies (GWAS) analyses. In conclusion, this study delineated a cross-species retinal cell atlas and uncovered the early pathological alterations in the retina of T2D mice.

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.

GENE EXPRESSION WITHIN A HUMAN CHOROIDAL NEOVASCULAR MEMBRANE USING SPATIAL TRANSCRIPTOMICS

Macular neovascularization is a relatively common and potentially visually devastating complication of age-related macular degeneration. In macular neovascularization, pathologic angiogenesis can originate from either the choroid or the retina, but we have limited understanding of how different cell types become dysregulated in this dynamic process. In this study, we performed spatial RNA sequencing on a human donor eye with macular neovascularization as well as a healthy control donor. We identified genes enriched within the area of macular neovascularization and used deconvolution algorithms to predict the originating cell type of these dysregulated genes. Within the area of neovascularization, endothelial cells were predicted to increase expression of genes related to Rho family GTPase signaling and integrin signaling. Likewise, VEGF and TGFB1 were identified as potential upstream regulators that could drive the observed gene expression changes produced by endothelial and retinal pigment epithelium cells in the macular neovascularization donor. These spatial gene expression profiles were compared to previous single-cell gene expression experiments in human age-related macular degeneration as well as a model of laser-induced neovascularization in mice. As a secondary aim, we also investigated spatial gene expression patterns within the macular neural retina and between the macular and peripheral choroid. We recapitulated previously described regional-specific gene expression patterns across both tissues. Overall, this study spatially analyzes gene expression across the retina, retinal pigment epithelium, and choroid in health and describes a set of candidate molecules that become dysregulated in macular neovascularization.

Single-cell sequencing of individual retinal organoids reveals determinants of cell fate heterogeneity

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA-seq multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10x methods produce highly concordant cell class compositions and then expand sci-Plex to analyze the cell class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to six weeks later. Our data show sci-Plex's potential to dramatically scale-up the analysis of treatment conditions on relevant human models.

Loss or gain of function? Effects of ion channel mutations on neuronal firing depend on the neuron type

Introduction

Clinically relevant mutations to voltage-gated ion channels, called channelopathies, alter ion channel function, properties of ionic currents, and neuronal firing. The effects of ion channel mutations are routinely assessed and characterized as loss of function (LOF) or gain of function (GOF) at the level of ionic currents. However, emerging personalized medicine approaches based on LOF/GOF characterization have limited therapeutic success. Potential reasons are among others that the translation from this binary characterization to neuronal firing is currently not well-understood-especially when considering different neuronal cell types. In this study, we investigate the impact of neuronal cell type on the firing outcome of ion channel mutations.

Methods

To this end, we simulated a diverse collection of single-compartment, conductance-based neuron models that differed in their composition of ionic currents. We systematically analyzed the effects of changes in ion current properties on firing in different neuronal types. Additionally, we simulated the effects of known mutations in KCNA1 gene encoding the K_V1.1 potassium channel subtype associated with episodic ataxia type 1 (EA1).

Results

These simulations revealed that the outcome of a given change in ion channel properties on neuronal excitability depends on neuron type, i.e., the properties and expression levels of the unaffected ionic currents.

Discussion

Consequently, neuron-type specific effects are vital to a full understanding of the effects of channelopathies on neuronal excitability and are an important step toward improving the efficacy and precision of personalized medicine approaches.

Overlap of Genetic Loci for Central Serous Chorioretinopathy With Age-Related Macular Degeneration

Importance

Central serous chorioretinopathy (CSC) is a serous maculopathy of unknown etiology. Two of 3 previously reported CSC genetic risk loci are also associated with AMD. Improved understanding of CSC genetics may broaden our understanding of this genetic overlap and unveil mechanisms in both diseases.

Objective

To identify novel genetic risk factors for CSC and compare genetic risk factors for CSC and AMD.

Design, Setting, and Participants

Using International Classification of Diseases, Ninth (ICD-9) and Tenth (ICD-10) Revision code-based inclusion and exclusion criteria, patients with CSC and controls were identified in both the FinnGen study and the Estonian Biobank (EstBB). Also included in a meta-analysis were previously reported patients with chronic CSC and controls. Data were analyzed from March 1 to September 31, 2022.

Main Outcomes and Measures

Genome-wide association studies (GWASs) were performed in the biobank-based cohorts followed by a meta-analysis of all cohorts. The expression of genes prioritized by the polygenic priority score and nearest-gene methods were assessed in cultured choroidal endothelial cells and public ocular single-cell RNA sequencing data sets. The predictive utility of polygenic scores (PGSs) for CSC and AMD were evaluated in the FinnGen study.

Results

A total of 1176 patients with CSC and 526 787 controls (312 162 female [59.3%]) were included in this analysis: 552 patients with CSC and 343 461 controls were identified in the FinnGen study, 103 patients with CSC and 178 573 controls were identified in the EstBB, and 521 patients with chronic CSC and 3577 controls were included in a meta-analysis. Two previously reported CSC risk loci were replicated (near CFH and GATA5) and 3 novel loci were identified (near CD34/46, NOTCH4, and PREX1). The CFH and NOTCH4 loci were associated with AMD but in the opposite direction. Prioritized genes showed increased expression in cultured choroidal endothelial cells compared with other genes in the loci (median [IQR] of log 2 [counts per million], 7.3 [0.6] vs 4.7 [3.7]; P = .004) and were differentially expressed in choroidal vascular endothelial cells in single-cell RNA sequencing data (mean [SD] fold change, 2.05 [0.38] compared with other cell types; P < 7.1 × 10-20). A PGS for AMD was predictive of reduced CSC risk (odds ratio, 0.76; 95% CI, 0.70-0.83 per +1 SD in AMD-PGS; P = 7.4 × 10-10). This association may have been mediated by loci containing complement genes.

Conclusions and Relevance

In this 3-cohort genetic association study, 5 genetic risk loci for CSC were identified, highlighting a likely role for genes involved in choroidal vascular function and complement regulation. Results suggest that polygenic AMD risk was associated with reduced risk of CSC and that this genetic overlap was largely due to loci containing complement genes.

Unique interface and dynamics of the complex of HSP90 with a specialized cochaperone AIPL1

Photoreceptor phosphodiesterase PDE6 is central for visual signal transduction. Maturation of PDE6 depends on a specialized chaperone complex of HSP90 with aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1). Disruption of PDE6 maturation underlies a severe form of retina degeneration. Here, we report a 3.9 Å cryoelectron microscopy (cryo-EM) structure of the complex of HSP90 with AIPL1. This structure reveals a unique interaction of the FK506-binding protein (FKBP)-like domain of AIPL1 with HSP90 at its dimer interface. Unusually, the N terminus AIPL1 inserts into the HSP90 lumen in a manner that was observed previously for HSP90 clients. Deletion of the 7 N-terminal residues of AIPL1 decreased its ability to cochaperone PDE6. Multi-body refinement of the cryo-EM data indicated large swing-like movements of AIPL1-FKBP. Modeling the complex of HSP90 with AIPL1 using crosslinking constraints indicated proximity of the mobile tetratricopeptide repeat (TPR) domain with the C-terminal domain of HSP90. Our study establishes a framework for future structural studies of PDE6 maturation.

Post-GWAS screening of candidate genes for refractive error in mutant zebrafish models

Genome-wide association studies (GWAS) have dissected numerous genetic factors underlying refractive errors (RE) such as myopia. Despite significant insights into understanding the genetic architecture of RE, few studies have validated and explored the functional role of candidate genes within these loci. To functionally follow-up on GWAS and characterize the potential role of candidate genes on the development of RE, we prioritized nine genes (TJP2, PDE11A, SHISA6, LAMA2, LRRC4C, KCNQ5, GNB3, RBFOX1, and GRIA4) based on biological and statistical evidence; and used CRISPR/cas9 to generate knock-out zebrafish mutants. These mutant fish were screened for abnormalities in axial length by spectral-domain optical coherence tomography and refractive status by eccentric photorefraction at the juvenile (2 months) and adult (4 months) developmental stage. We found a significantly increased axial length and myopic shift in refractive status in three of our studied mutants, indicating a potential involvement of the human orthologs (LAMA2, LRRC4C, and KCNQ5) in myopia development. Further, in-situ hybridization studies showed that all three genes are expressed throughout the zebrafish retina. Our zebrafish models provide evidence of a functional role of these three genes in refractive error development and offer opportunities to elucidate pathways driving the retina-to-sclera signaling cascade that leads to myopia.

Retinal regions shape human and murine Müller cell proteome profile and functionality

The human macula is a highly specialized retinal region with pit-like morphology and rich in cones. How Müller cells, the principal glial cell type in the retina, are adapted to this environment is still poorly understood. We compared proteomic data from cone- and rod-rich retinae from human and mice and identified different expression profiles of cone- and rod-associated Müller cells that converged on pathways representing extracellular matrix and cell adhesion. In particular, epiplakin (EPPK1), which is thought to play a role in intermediate filament organization, was highly expressed in macular Müller cells. Furthermore, EPPK1 knockout in a human Müller cell-derived cell line led to a decrease in traction forces as well as to changes in cell size, shape, and filopodia characteristics. We here identified EPPK1 as a central molecular player in the region-specific architecture of the human retina, which likely enables specific functions under the immense mechanical loads in vivo.

Large-Scale Single-Nucleus RNA Sequencing Compatible with Complex Archived Samples

Transcriptome profiling at single-cell resolution allows us to identify and assess functional cell types and cellular states, including those within degenerating ocular tissues in retinitis pigmentosa. The technology is particularly valuable when studying tissues with high cellular heterogeneity, or when specific cell types are of interest. In this chapter, we introduce a detailed protocol of a medium-throughput single-nucleus RNA sequencing technique that utilizes frozen tissue as input sample. This protocol can be executed by any researcher with basic training in molecular biology techniques. With this protocol, a single experimenter can easily process two samples per day up to cDNA amplification, and library preparations can be done in batches of 8. Routinely we can obtain ~20 K nuclei per eye from 3 to 4 library preparations.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

Immune response in retinal degenerative diseases - Time to rethink?

Retinal degeneration comprises a group of diseases whereby either the retinal neurons or the neurovascular unit degenerates leading to the loss of visual function. Although the initial cause varies in different conditions, inflammation is known to play an important role in disease pathogenesis. Recent advances in molecular and cell biology and systems biology have yielded unexpected findings, including the heterogeneity of immune cells in the degenerative retina, bidirectional neuron-microglia cross talk, and links to the gut microbiome. Here we discuss the immune response in retinal degenerative conditions, taking into account both regional (retinal) and systemic factors. We propose to classify retinal degeneration into dry and wet forms based on whether the blood-retinal barrier (BRB) is breached and fluid is accumulated in retinal parenchyma. The dry form has a relatively intact BRB and is characterised by progressive retinal thinning. Immune response to degenerative insults is dominated by the retinal defence system, which remains to be regulated by neurons. In contrast, the wet form has retinal oedema due to BRB damaged. Inflammation is executed by infiltrating immune cells as well as the retinal defence system. The gut microbiome will have easy access to the retina in wet retinal degeneration and may affect significantly retinal immune response.

Transcriptional mapping of the macaque retina and RPE-choroid reveals conserved inter-tissue transcription drivers and signaling pathways

The macula and fovea comprise a highly sensitive visual detection tissue that is susceptible to common disease processes like age-related macular degeneration (AMD). Our understanding of the molecular determinants of high acuity vision remains unclear, as few model organisms possess a human-like fovea. We explore transcription factor networks and receptor-ligand interactions to elucidate tissue interactions in the macula and peripheral retina and concomitant changes in the underlying retinal pigment epithelium (RPE)/choroid. Poly-A selected, 100 bp paired-end RNA-sequencing (RNA-seq) was performed across the macular/foveal, perimacular, and temporal peripheral regions of the neural retina and RPE/choroid tissues of four adult Rhesus macaque eyes to characterize region- and tissue-specific gene expression. RNA-seq reads were mapped to both the macaque and human genomes for maximum alignment and analyzed for differential expression and Gene Ontology (GO) enrichment. Comparison of the neural retina and RPE/choroid tissues indicated distinct, contiguously changing gene expression profiles from fovea through perimacula to periphery. Top GO enrichment of differentially expressed genes in the RPE/choroid included cell junction organization and epithelial cell development. Expression of transcriptional regulators and various disease-associated genes show distinct location-specific preference and retina-RPE/choroid tissue-tissue interactions. Regional gene expression changes in the macaque retina and RPE/choroid is greater than that found in previously published transcriptome analysis of the human retina and RPE/choroid. Further, conservation of human macula-specific transcription factor profiles and gene expression in macaque tissues suggest a conservation of programs required for retina and RPE/choroid function and disease susceptibility.

Transcriptome analysis of AAV-induced retinopathy models expressing human VEGF, TNF-α, and IL-6 in murine eyes

Retinopathies are multifactorial diseases with complex pathologies that eventually lead to vision loss. Animal models facilitate the understanding of the pathophysiology and identification of novel treatment options. However, each animal model reflects only specific disease aspects and understanding of the specific molecular changes in most disease models is limited. Here, we conducted transcriptome analysis of murine ocular tissue transduced with recombinant Adeno-associated viruses (AAVs) expressing either human VEGF-A, TNF-α, or IL-6. VEGF expression led to a distinct regulation of extracellular matrix (ECM)-associated genes. In contrast, both TNF-α and IL-6 led to more comparable gene expression changes in interleukin signaling, and the complement cascade, with TNF-α-induced changes being more pronounced. Furthermore, integration of single cell RNA-Sequencing data suggested an increase of endothelial cell-specific marker genes by VEGF, while TNF-α expression increased the expression T-cell markers. Both TNF-α and IL-6 expression led to an increase in macrophage markers. Finally, transcriptomic changes in AAV-VEGF treated mice largely overlapped with gene expression changes observed in the oxygen-induced retinopathy model, especially regarding ECM components and endothelial cell-specific gene expression. Altogether, our study represents a valuable investigation of gene expression changes induced by VEGF, TNF-α, and IL-6 and will aid researchers in selecting appropriate animal models for retinopathies based on their agreement with the human pathophysiology.

Aster la vista: Unraveling the biochemical basis of carotenoid homeostasis in the human retina

Carotenoids play pivotal roles in vision as light filters and precursor of chromophore. Many vertebrates also display the colorful pigments as ornaments in bare skin parts and feathers. Proteins involved in the transport and metabolism of these lipids have been identified including class B scavenger receptors and carotenoid cleavage dioxygenases. Recent research implicates members of the Aster protein family, also known as GRAM domain-containing (GRAMD), in carotenoid metabolism. These multi-domain proteins facilitate the intracellular movement of carotenoids from their site of cellular uptake by scavenger receptors to the site of their metabolic processing by carotenoid cleavage dioxygenases. We provide a model how the coordinated interplay of these proteins and their differential expression establishes carotenoid distribution patterns and function in tissues, with particular emphasis on the human retina.

HBEGF-TNF induce a complex outer retinal pathology with photoreceptor cell extrusion in human organoids

Human organoids could facilitate research of complex and currently incurable neuropathologies, such as age-related macular degeneration (AMD) which causes blindness. Here, we establish a human retinal organoid system reproducing several parameters of the human retina, including some within the macula, to model a complex combination of photoreceptor and glial pathologies. We show that combined application of TNF and HBEGF, factors associated with neuropathologies, is sufficient to induce photoreceptor degeneration, glial pathologies, dyslamination, and scar formation: These develop simultaneously and progressively as one complex phenotype. Histologic, transcriptome, live-imaging, and mechanistic studies reveal a previously unknown pathomechanism: Photoreceptor neurodegeneration via cell extrusion. This could be relevant for aging, AMD, and some inherited diseases. Pharmacological inhibitors of the mechanosensor PIEZO1, MAPK, and actomyosin each avert pathogenesis; a PIEZO1 activator induces photoreceptor extrusion. Our model offers mechanistic insights, hypotheses for neuropathologies, and it could be used to develop therapies to prevent vision loss or to regenerate the retina in patients suffering from AMD and other diseases.

Retinal Pigment Epithelium-Secreted VEGF-A Induces Alpha-2-Macroglobulin Expression in Endothelial Cells

Alpha-2-macroglobulin (A2M) is a protease inhibitor that regulates extracellular matrix (ECM) stability and turnover. Here, we show that A2M is expressed by endothelial cells (ECs) from human eye choroid. We demonstrate that retinal pigment epithelium (RPE)-conditioned medium induces A2M expression specifically in ECs. Experiments using chemical inhibitors, blocking antibodies, and recombinant proteins revealed a key role of VEGF-A in RPE-mediated A2M induction in ECs. Furthermore, incubation of ECs with RPE-conditioned medium reduces matrix metalloproteinase-2 gelatinase activity of culture supernatants, which is partially restored after A2M knockdown in ECs. We propose that dysfunctional RPE or choroidal blood vessels, as observed in retinal diseases such as age-related macular degeneration, may disrupt the crosstalk mechanism we describe here leading to alterations in the homeostasis of choroidal ECM, Bruch's membrane and visual function.

Boolean implication analysis of single-cell data predicts retinal cell type markers

Background

The retina is a complex tissue containing multiple cell types that are essential for vision. Understanding the gene expression patterns of various retinal cell types has potential applications in regenerative medicine. Retinal organoids (optic vesicles) derived from pluripotent stem cells have begun to yield insights into the transcriptomics of developing retinal cell types in humans through single cell RNA-sequencing studies. Previous methods of gene reporting have relied upon techniques in vivo using microarray data, or correlational and dimension reduction methods for analyzing single cell RNA-sequencing data computationally. We aimed to develop a state-of-the-art Boolean method that filtered out noise, could be applied to a wide variety of datasets and lent insight into gene expression over differentiation.

Results

Here, we present a bioinformatic approach using Boolean implication to discover genes which are retinal cell type-specific or involved in retinal cell fate. We apply this approach to previously published retina and retinal organoid datasets and improve upon previously published correlational methods. Our method improves the prediction accuracy of marker genes of retinal cell types and discovers several new high confidence cone and rod-specific genes.

Conclusions

The results of this study demonstrate the benefits of a Boolean approach that considers asymmetric relationships. We have shown a statistically significant improvement from correlational, symmetric methods in the prediction accuracy of retinal cell-type specific genes.

Furthermore, our method contains no cell or tissue-specific tuning and hence could impact other areas of gene expression analyses in cancer and other human diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04915-4.

Cone-Driven Retinal Responses Are Shaped by Rod But Not Cone HCN1

Signal integration of converging neural circuits is poorly understood. One example is in the retina where the integration of rod and cone signaling is responsible for the large dynamic range of vision. The relative contribution of rods versus cones is dictated by a complex function involving background light intensity and stimulus temporal frequency. One understudied mechanism involved in coordinating rod and cone signaling onto the shared retinal circuit is the hyperpolarization activated current (Ih) mediated by hyperpolarization-activated cyclic nucleotide-gated 1 (HCN1) channels expressed in rods and cones. Ih opposes membrane hyperpolarization driven by activation of the phototransduction cascade and modulates the strength and kinetics of the photoreceptor voltage response. We examined conditional knock-out (KO) of HCN1 from mouse rods using electroretinography (ERG). In the absence of HCN1, rod responses are prolonged in dim light which altered the response to slow modulation of light intensity both at the level of retinal signaling and behavior. Under brighter intensities, cone-driven signaling was suppressed. To our surprise, conditional KO of HCN1 from mouse cones had no effect on cone-mediated signaling. We propose that Ih is dispensable in cones because of the high level of temporal control of cone phototransduction. Thus, HCN1 is required for cone-driven retinal signaling only indirectly by modulating the voltage response of rods to limit their output.SIGNIFICANCE STATEMENT Hyperpolarization gated hyperpolarization-activated cyclic nucleotide-gated 1 (HCN1) channels carry a feedback current that helps to reset light-activated photoreceptors. Using conditional HCN1 knock-out (KO) mice we show that ablating HCN1 from rods allows rods to signal in bright light when they are normally shut down. Instead of enhancing vision this results in suppressing cone signaling. Conversely, ablating HCN1 from cones was of no consequence. This work provides novel insights into the integration of rod and cone signaling in the retina and challenges our assumptions about the role of HCN1 in cones.

Single-cell characterization of malignant phenotypes and microenvironment alteration in retinoblastoma

Retinoblastoma (RB) is the most common primary intraocular malignancy of childhood. It is known that the tumor microenvironment (TME) regulates tumorigenesis and metastasis. However, how the malignant progression in RB is determined by the heterogeneity of tumor cells and TME remains uncharacterized. Here, we conducted integrative single-cell transcriptome and whole-exome sequencing analysis of RB patients with detailed pathological and clinical measurements. By single-cell transcriptomic sequencing, we profiled around 70,000 cells from tumor samples of seven RB patients. We identified that the major cell types in RB were cone precursor-like (CP-like) and MKI67+ cone precursor (MKI67+ CP) cells. By integrating copy number variation (CNV) analysis, we found that RB samples had large clonal heterogeneity, where the malignant MKI67+ CP cells had significantly larger copy number changes. Enrichment analysis revealed that the conversion of CP-like to MKI67+ CP resulted in the loss of photoreceptor function and increased cell proliferation ability. The TME in RB was composed of tumor-associated macrophages (TAMs), astrocyte-like, and cancer-associated fibroblasts (CAFs). Particularly, during the invasion process, TAMs created an immunosuppressive environment, in which the proportion of TAMs decreased, M1-type macrophage was lost, and the TAMs-related immune functions were depressed. Finally, we identified that TAMs regulated tumor cells through GRN and MIF signaling pathways, while TAMs self-regulated through inhibition of CCL and GALECTIN signaling pathways during the invasion process. Altogether, our study creates a detailed transcriptomic map of RB with single-cell characterization of malignant phenotypes and provides novel molecular insights into the occurrence and progression of RB.

Aster proteins mediate carotenoid transport in mammalian cells

Carotenoid pigments accumulate in specific patterns in vertebrate tissues and play important roles as colorants, chromophores, and hormone precursors. However, proteins that facilitate transportation of these lipophilic pigments within cells have not been identified. We provide evidence that Aster proteins are key components for this process and show that they bind the pigments with high affinity. We observed in mice that carotenoids accumulate in tissues that express Aster-B and this accumulation can be prevented by enzymatic turnover by the BCO2 protein. Accordingly, we found opposing expression patterns of the Aster-B protein and BCO2 in the human retina that seemingly contribute to the unique carotenoid concentration in the macula lutea.

Some mammalian tissues uniquely concentrate carotenoids, but the underlying biochemical mechanism for this accumulation has not been fully elucidated. For instance, the central retina of the primate eyes displays high levels of the carotenoids, lutein, and zeaxanthin, whereas the pigments are largely absent in rodent retinas. We previously identified the scavenger receptor class B type 1 and the enzyme β-carotene-oxygenase-2 (BCO2) as key components that determine carotenoid concentration in tissues. We now provide evidence that Aster (GRAM-domain-containing) proteins, recently recognized for their role in nonvesicular cholesterol transport, engage in carotenoid metabolism. Our analyses revealed that the StART-like lipid binding domain of Aster proteins can accommodate the bulky pigments and bind them with high affinity. We further showed that carotenoids and cholesterol compete for the same binding site. We established a bacterial test system to demonstrate that the StART-like domains of mouse and human Aster proteins can extract carotenoids from biological membranes. Mice deficient for the carotenoid catabolizing enzyme BCO2 concentrated carotenoids in Aster-B protein-expressing tissues such as the adrenal glands. Remarkably, Aster-B was expressed in the human but not in the mouse retina. Within the retina, Aster-B and BCO2 showed opposite expression patterns in central versus peripheral parts. Together, our study unravels the biochemical basis for intracellular carotenoid transport and implicates Aster-B in the pathway for macula pigment concentration in the human retina.

Retinoic acid signaling mediates peripheral cone photoreceptor survival in a mouse model of retina degeneration

Retinitis Pigmentosa (RP) is a progressive, debilitating visual disorder caused by mutations in a diverse set of genes. In both humans with RP and mouse models of RP, rod photoreceptor dysfunction leads to loss of night vision, and is followed by secondary cone photoreceptor dysfunction and degeneration, leading to loss of daylight color vision. A strategy to prevent secondary cone death could provide a general RP therapy to preserve daylight color vision regardless of the underlying mutation. In mouse models of RP, cones in the peripheral retina survive long-term, despite complete rod loss. The mechanism for such peripheral cone survival had not been explored. Here, we found that active retinoic acid (RA) signaling in peripheral Muller glia is necessary for the abnormally long survival of these peripheral cones. RA depletion by conditional knockout of RA synthesis enzymes, or overexpression of an RA degradation enzyme, abrogated the extended survival of peripheral cones. Conversely, constitutive activation of RA signaling in the central retina promoted long-term cone survival. These results indicate that RA signaling mediates the prolonged peripheral cone survival in the rd1 mouse model of retinal degeneration, and provide a basis for a generic strategy for cone survival in the many diseases that lead to loss of cone-mediated vision.

In-Depth Molecular Profiling Specifies Human Retinal Microglia Identity

Microglia are the tissue-resident macrophages of the retina and brain, being critically involved in organ development, tissue homeostasis, and response to cellular damage. Until now, little is known about the molecular signature of human retinal microglia and how it differs from the one of brain microglia and peripheral monocytes. In addition, it is not yet clear to what extent murine retinal microglia resemble those of humans, which represents an important prerequisite for translational research. The present study applies fluorescence-activated cell sorting to isolate human retinal microglia from enucleated eyes and compares their transcriptional profile with the one of whole retinal tissue, human brain microglia as well as classical, intermediate and non-classical monocytes. Finally, human retinal microglia are compared to murine retinal microglia, isolated from Cx3cr1^GFP/+ mice. Whereas human retinal microglia exhibited a high grade of similarity in comparison to their counterparts in the brain, several enriched genes were identified in retinal microglia when compared to whole retinal tissue, as well as classical, intermediate, and non-classical monocytes. In relation to whole retina sequencing, several risk genes associated with age-related macular degeneration (AMD) and diabetic retinopathy (DR) were preferentially expressed in retinal microglia, indicating their potential pathophysiological involvement. Although a high degree of similarity was observed between human and murine retinal microglia, several species-specific genes were identified, which should be kept in mind when employing mouse models to investigate retinal microglia biology. In summary, this study provides detailed insights into the molecular profile of human retinal microglia, identifies a plethora of tissue-specific and species-specific genes in comparison to human brain microglia and murine retinal microglia, and thus highlights the significance of retinal microglia in human retinal diseases and for translational research approaches.

Variability in Retinal Neuron Populations and Associated Variations in Mass Transport Systems of the Retina in Health and Aging

Aging is associated with a broad range of visual impairments that can have dramatic consequences on the quality of life of those impacted. These changes are driven by a complex series of alterations affecting interactions between multiple cellular and extracellular elements. The resilience of many of these interactions may be key to minimal loss of visual function in aging; yet many of them remain poorly understood. In this review, we focus on the relation between retinal neurons and their respective mass transport systems. These metabolite delivery systems include the retinal vasculature, which lies within the inner portion of the retina, and the choroidal vasculature located externally to the retinal tissue. A framework for investigation is proposed and applied to identify the structures and processes determining retinal mass transport at the cellular and tissue levels. Spatial variability in the structure of the retina and changes observed in aging are then harnessed to explore the relation between variations in neuron populations and those seen among retinal metabolite delivery systems. Existing data demonstrate that the relation between inner retinal neurons and their mass transport systems is different in nature from that observed between the outer retina and choroid. The most prominent structural changes observed across the eye and in aging are seen in Bruch’s membrane, which forms a selective barrier to mass transfers at the interface between the choroidal vasculature and the outer retina.

Müller Cell Molecular Heterogeneity: Facts and Predictions

The retina was historically considered as an “approachable part of the brain”; advantageous, for its simplicity, to use as a model organ for deciphering cellular and molecular mechanisms underlying physiology and pathology of the nervous system. However, the most relevant discoveries arise precisely from unveiling the complexity of the retina. A complexity that partially relies on the layered organization of an extended variety of specialized neuronal and glial cellular types and subtypes. Based on functional, morphological or transcriptome data, over 40 subtypes of retinal ganglion cells or 60 subtypes of retinal amacrine cells have been described. A high degree of specialization, that may lead to segregation into functionally diverse subtypes, is also conceivable for Müller cells, a pleiotropic glial component of all vertebrate retinas. The essential role of Müller glia in retinal homeostasis maintenance involves participation in structural, metabolic and intercellular communication processes. Additionally, they are the only retinal cells that possess regenerative potential in response to injury or disease, and thus may be considered as therapeutic tools. In the assumption that functional heterogeneity might be driven by molecular heterogeneity this review aims to compile emerging evidence that could broaden our understanding of Müller cell biology and retinal physiology.

Summary statement

Müller glial cells exert multiple essential functions in retinal physiology and retinopathies reflecting perhaps the existence of distinct Müller cellular subpopulations. Harnessing Müller cell heterogeneity may serve to enhance new therapeutic approaches for retinal disease.

Carotenoid modifying enzymes in metazoans

Carotenoids constitute an essential dietary component of animals and other non-carotenogenic species which use these pigments in both their modified and unmodified forms. Animals utilize uncleaved carotenoids to mitigate light damage and oxidative stress and to signal fitness and health. Carotenoids also serve as precursors of apocarotenoids including retinol, and its retinoid metabolites, which carry out essential functions in animals by forming the visual chromophore 11-cis-retinaldehyde. Retinoids, such as all-trans-retinoic acid, can also act as ligands of nuclear hormone receptors. The fact that enzymes and biochemical pathways responsible for the metabolism of carotenoids in animals bear resemblance to the ones in plants and other carotenogenic species suggests an evolutionary relationship. We will explore some of the modes of transmission of carotenoid genes from carotenogenic species to metazoans. This apparent relationship has been successfully exploited in the past to identify and characterize new carotenoid and retinoid modifying enzymes. We will review approaches used to identify putative animal carotenoid enzymes, and we will describe methods used to functionally validate and analyze the biochemistry of carotenoid modifying enzymes encoded by animals.

Retinal ciliopathies through the lens of Bardet-Biedl Syndrome: Past, present and future

The primary cilium is a highly specialized and evolutionary conserved organelle in eukaryotes that plays a significant role in cell signaling and trafficking. Over the past few decades tremendous progress has been made in understanding the physiology of cilia and the underlying pathomechanisms of various ciliopathies. Syndromic ciliopathies consist of a group of disorders caused by ciliary dysfunction or abnormal ciliogenesis. These disorders have multiorgan involvement in addition to retinal degeneration underscoring the ubiquitous distribution of primary cilia in different cell types. Genotype-phenotype correlation is often challenging due to the allelic heterogeneity and pleiotropy of these disorders. In this review, we discuss the clinical and genetic features of syndromic ciliopathies with a focus on Bardet-Biedl syndrome (BBS) as a representative disorder. We discuss the structure and function of primary cilia and their role in retinal photoreceptors. We describe the progress made thus far in understanding the functional and genetic characterization including expression quantitative trait locus (eQTL) analysis of BBS genes. In the future directions section, we discuss the emerging technologies, such as gene therapy, as well as anticipated challenges and their implications in therapeutic development for ciliopathies.

Live imaging of retinotectal mapping reveals topographic map dynamics and a previously undescribed role for Contactin 2 in map sharpening

Organization of neuronal connections into topographic maps is essential for processing information. Yet, our understanding of topographic mapping has remained limited by our inability to observe maps forming and refining directly in vivo. Here, we used Cre-mediated recombination of a new colorswitch reporter in zebrafish to generate the first transgenic model allowing the dynamic analysis of retinotectal mapping in vivo. We found that the antero-posterior retinotopic map forms early but remains dynamic, with nasal and temporal retinal axons expanding their projection domains over time. Nasal projections initially arborize in the anterior tectum but progressively refine their projection domain to the posterior tectum, leading to the sharpening of the retinotopic map along the antero-posterior axis. Finally, using a CRISPR-mediated mutagenesis approach, we demonstrate that the refinement of nasal retinal projections requires the adhesion molecule Contactin 2. Altogether, our study provides the first analysis of a topographic map maturing in real time in a live animal and opens new strategies for dissecting the molecular mechanisms underlying precise topographic mapping in vertebrates.

Disturbed retinoid metabolism upon loss of rlbp1a impairs cone function and leads to subretinal lipid deposits and photoreceptor degeneration in the zebrafish retina

The RLBP1 gene encodes the 36 kDa cellular retinaldehyde-binding protein, CRALBP, a soluble retinoid carrier, in the visual cycle of the eyes. Mutations in RLBP1 are associated with recessively inherited clinical phenotypes, including Bothnia dystrophy, retinitis pigmentosa, retinitis punctata albescens, fundus albipunctatus, and Newfoundland rod–cone dystrophy. However, the etiology of these retinal disorders is not well understood. Here, we generated homologous zebrafish models to bridge this knowledge gap. Duplication of the rlbp1 gene in zebrafish and cell-specific expression of the paralogs rlbp1a in the retinal pigment epithelium and rlbp1b in Müller glial cells allowed us to create intrinsically cell type-specific knockout fish lines. Using rlbp1a and rlbp1b single and double mutants, we investigated the pathological effects on visual function. Our analyses revealed that rlbp1a was essential for cone photoreceptor function and chromophore metabolism in the fish eyes. rlbp1a-mutant fish displayed reduced chromophore levels and attenuated cone photoreceptor responses to light stimuli. They accumulated 11-cis and all-trans-retinyl esters which displayed as enlarged lipid droplets in the RPE reminiscent of the subretinal yellow-white lesions in patients with RLBP1 mutations. During aging, these fish developed retinal thinning and cone and rod photoreceptor dystrophy. In contrast, rlbp1b mutants did not display impaired vision. The double mutant essentially replicated the phenotype of the rlbp1a single mutant. Together, our study showed that the rlbp1a zebrafish mutant recapitulated many features of human blinding diseases caused by RLBP1 mutations and provided novel insights into the pathways for chromophore regeneration of cone photoreceptors.

Building the mega single-cell transcriptome ocular meta-atlas

BACKGROUND

The development of highly scalable single-cell transcriptome technology has resulted in the creation of thousands of datasets, >30 in the retina alone. Analyzing the transcriptomes between different projects is highly desirable because this would allow for better assessment of which biological effects are consistent across independent studies. However it is difficult to compare and contrast data across different projects because there are substantial batch effects from computational processing, single-cell technology utilized, and the natural biological variation. While many single-cell transcriptome-specific batch correction methods purport to remove the technical noise, it is difficult to ascertain which method functions best.

RESULTS

We developed a lightweight R package (scPOP, single-cell Pick Optimal Parameters) that brings in batch integration methods and uses a simple heuristic to balance batch merging and cell type/cluster purity. We use this package along with a Snakefile-based workflow system to demonstrate how to optimally merge 766,615 cells from 33 retina datsets and 3 species to create a massive ocular single-cell transcriptome meta-atlas.

CONCLUSIONS

This provides a model for how to efficiently create meta-atlases for tissues and cells of interest.

A single-cell guide to retinal development: Cell fate decisions of multipotent retinal progenitors in scRNA-seq

Recent advances in high throughput single-cell RNA sequencing (scRNA-seq) technology have enabled the simultaneous transcriptomic profiling of thousands of individual cells in a single experiment. To investigate the intrinsic process of retinal development, researchers have leveraged this technology to quantify gene expression in retinal cells across development, in multiple species, and from numerous important models of human disease. In this review, we summarize recent applications of scRNA-seq and discuss how these datasets have complemented and advanced our understanding of retinal progenitor cell competence, cell fate specification, and differentiation. Finally, we also highlight the outstanding questions in the field that advances in single-cell data generation and analysis will soon be able to answer.