Single-Cell Analysis of Human Retina Identifies Evolutionarily Conserved and Species-Specific Mechanisms Controlling Development

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.

Non-canonical Wnt pathway expression in the developing mouse and human retina

The non-canonical Wnt pathway is an evolutionarily conserved pathway essential for tissue patterning and development across species and tissues. In mammals, this pathway plays a role in neuronal migration, dendritogenesis, axon growth, and synapse formation. However, its role in development and synaptogenesis of the human retina remains less established. In order to address this knowledge gap, we analyzed publicly available single-cell RNA sequencing (scRNAseq) datasets for mouse retina, human retina, and human retinal organoids over multiple developmental time points during outer retinal maturation. We identified ligands, receptors, and mediator genes with a putative role in retinal development, including those with novel or species-specific expression, and validated this expression using fluorescence in situ hybridization (FISH). By quantifying outer nuclear layer (ONL) versus inner nuclear layer (INL) expression, we provide evidence for the differential expression of certain non-canonical Wnt signaling components in the developing mouse and human retina during outer plexiform layer (OPL) development. Importantly, we identified distinct expression patterns of mouse and human FZD3 and WNT10A, as well as previously undescribed expression, such as for mouse Wnt2b in Chat+ starburst amacrine cells. Human retinal organoids largely recapitulated the human non-canonical Wnt pathway expression. Together, this work provides the basis for further study of non-canonical Wnt signaling in mouse and human retinal development and synaptogenesis.

Temporal patterning of the vertebrate developing neural tube

The chronologically ordered generation of distinct cell types is essential for the establishment of neuronal diversity and the formation of neuronal circuits. Recently, single-cell transcriptomic analyses of various areas of the developing vertebrate nervous system have provided evidence for the existence of a shared temporal patterning program that partitions neurons based on the timing of neurogenesis. In this review, I summarize the findings that lead to the proposal of this shared temporal program before focusing on the developing spinal cord to discuss how temporal patterning in general and this program specifically contributes to the ordered formation of neuronal circuits.

Perivascular neurons instruct 3D vascular lattice formation via neurovascular contact

The vasculature of the central nervous system is a 3D lattice composed of laminar vascular beds interconnected by penetrating vessels. The mechanisms controlling 3D lattice network formation remain largely unknown. Combining viral labeling, genetic marking, and single-cell profiling in the mouse retina, we discovered a perivascular neuronal subset, annotated as Fam19a4/Nts-positive retinal ganglion cells (Fam19a4/Nts-RGCs), directly contacting the vasculature with perisomatic endfeet. Developmental ablation of Fam19a4/Nts-RGCs led to disoriented growth of penetrating vessels near the ganglion cell layer (GCL), leading to a disorganized 3D vascular lattice. We identified enriched PIEZO2 expression in Fam19a4/Nts-RGCs. Piezo2 loss from all retinal neurons or Fam19a4/Nts-RGCs abolished the direct neurovascular contacts and phenocopied the Fam19a4/Nts-RGC ablation deficits. The defective vascular structure led to reduced capillary perfusion and sensitized the retina to ischemic insults. Furthermore, we uncovered a Piezo2-dependent perivascular granule cell subset for cerebellar vascular patterning, indicating neuronal Piezo2-dependent 3D vascular patterning in the brain.

Retinoid Synthesis Regulation by Retinal Cells in Health and Disease

Vision starts in retinal photoreceptors when specialized proteins (opsins) sense photons via their covalently bonded vitamin A derivative 11cis retinaldehyde (11cis-RAL). The reaction of non-enzymatic aldehydes with amino groups lacks specificity, and the reaction products may trigger cell damage. However, the reduced synthesis of 11cis-RAL results in photoreceptor demise and suggests the need for careful control over 11cis-RAL handling by retinal cells. This perspective focuses on retinoid(s) synthesis, their control in the adult retina, and their role during retina development. It also explores the potential importance of 9cis vitamin A derivatives in regulating retinoid synthesis and their impact on photoreceptor development and survival. Additionally, recent advancements suggesting the pivotal nature of retinoid synthesis regulation for cone cell viability are discussed.

Age-related dysregulation of the retinal transcriptome in African turquoise killifish

Age-related vision loss caused by retinal neurodegenerative pathologies is becoming more prevalent in our ageing society. To understand the physiological and molecular impact of ageing on retinal homeostasis, we used the short-lived African turquoise killifish, a model known to naturally develop central nervous system (CNS) ageing hallmarks and vision loss. Bulk and single-cell RNA-sequencing (scRNAseq) of three age groups (6-, 12-, and 18-week-old) identified transcriptional ageing fingerprints in the killifish retina, unveiling pathways also identified in the aged brain, including oxidative stress, gliosis, and inflammageing. These findings were comparable to observations in the ageing mouse retina. Additionally, transcriptional changes in genes related to retinal diseases, such as glaucoma and age-related macular degeneration, were observed. The cellular heterogeneity in the killifish retina was characterized, confirming the presence of all typical vertebrate retinal cell types. Data integration from age-matched samples between the bulk and scRNAseq experiments revealed a loss of cellular specificity in gene expression upon ageing, suggesting potential disruption in transcriptional homeostasis. Differential expression analysis within the identified cell types highlighted the role of glial/immune cells as important stress regulators during ageing. Our work emphasizes the value of the fast-ageing killifish in elucidating molecular signatures in age-associated retinal disease and vision decline. This study contributes to the understanding of how age-related changes in molecular pathways may impact CNS health, providing insights that may inform future therapeutic strategies for age-related pathologies.

Pluripotent stem cell-derived models of retinal disease: Elucidating pathogenesis, evaluating novel treatments, and estimating toxicity

Blindness poses a growing global challenge, with approximately 26% of cases attributed to degenerative retinal diseases. While gene therapy, optogenetic tools, photosensitive switches, and retinal prostheses offer hope for vision restoration, these high-cost therapies will benefit few patients. Understanding retinal diseases is therefore key to advance effective treatments, requiring in vitro models replicating pathology and allowing quantitative assessments for drug discovery. Pluripotent stem cells (PSCs) provide a unique solution given their limitless supply and ability to differentiate into light-responsive retinal tissues encompassing all cell types. This review focuses on the history and current state of photoreceptor and retinal pigment epithelium (RPE) cell generation from PSCs. We explore the applications of this technology in disease modelling, experimental therapy testing, biomarker identification, and toxicity studies. We consider challenges in scalability, standardisation, and reproducibility, and stress the importance of incorporating vasculature and immune cells into retinal organoids. We advocate for high-throughput automation in data acquisition and analyses and underscore the value of advanced micro-physiological systems that fully capture the interactions between the neural retina, RPE, and choriocapillaris.

Ciliopathy patient variants reveal organelle-specific functions for TUBB4B in axonemal microtubules

Tubulin, one of the most abundant cytoskeletal building blocks, has numerous isotypes in metazoans encoded by different conserved genes. Whether these distinct isotypes form cell type- and context-specific microtubule structures is poorly understood. Based on a cohort of 12 patients with primary ciliary dyskinesia as well as mouse mutants, we identified and characterized variants in the TUBB4B isotype that specifically perturbed centriole and cilium biogenesis. Distinct TUBB4B variants differentially affected microtubule dynamics and cilia formation in a dominant-negative manner. Structure-function studies revealed that different TUBB4B variants disrupted distinct tubulin interfaces, thereby enabling stratification of patients into three classes of ciliopathic diseases. These findings show that specific tubulin isotypes have distinct and nonredundant subcellular functions and establish a link between tubulinopathies and ciliopathies.

Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

Deciphering the spatiotemporal transcriptional and chromatin accessibility of human retinal organoid development at the single-cell level

Molecular information on the early stages of human retinal development remains scarce due to limitations in obtaining early human eye samples. Pluripotent stem cell-derived retinal organoids (ROs) provide an unprecedented opportunity for studying early retinogenesis. Using a combination of single cell RNA-seq and spatial transcriptomics we present for the first-time a single cell spatiotemporal transcriptome of RO development. Our data demonstrate that ROs recapitulate key events of retinogenesis including optic vesicle/cup formation, presence of a putative ciliary margin zone, emergence of retinal progenitor cells and their orderly differentiation to retinal neurons. Combining the scRNA- with scATAC-seq data, we were able to reveal cell-type specific transcription factor binding motifs on accessible chromatin at each stage of organoid development, and to show that chromatin accessibility is highly correlated to the developing human retina, but with some differences in the temporal emergence and abundance of some of the retinal neurons.

Single-Cell RNA Sequencing in Organ and Cell Transplantation

Single-cell RNA sequencing is a high-throughput novel method that provides transcriptional profiling of individual cells within biological samples. This method typically uses microfluidics systems to uncover the complex intercellular communication networks and biological pathways buried within highly heterogeneous cell populations in tissues. One important application of this technology sits in the fields of organ and stem cell transplantation, where complications such as graft rejection and other post-transplantation life-threatening issues may occur. In this review, we first focus on research in which single-cell RNA sequencing is used to study the transcriptional profile of transplanted tissues. This technology enables the analysis of the donor and recipient cells and identifies cell types and states associated with transplant complications and pathologies. We also review the use of single-cell RNA sequencing in stem cell implantation. This method enables studying the heterogeneity of normal and pathological stem cells and the heterogeneity in cell populations. With their remarkably rapid pace, the single-cell RNA sequencing methodologies will potentially result in breakthroughs in clinical transplantation in the coming years.

Age-related dysregulation of the retinal transcriptome in African turquoise killifish

Age-related vision loss caused by retinal neurodegenerative pathologies is becoming more prevalent in our ageing society. To understand the physiological and molecular impact of ageing on retinal homeostasis, we used the short-lived African turquoise killifish, a model known to naturally develop central nervous system (CNS) ageing hallmarks and vision loss. Bulk and single-cell RNA-sequencing (scRNA-seq) of three age groups (6-, 12-, and 18-week-old) identified transcriptional ageing fingerprints in the killifish retina, unveiling pathways also identified in the aged brain, including oxidative stress, gliosis, and inflammageing. These findings were comparable to observations in ageing mouse retina. Additionally, transcriptional changes in genes related to retinal diseases, such as glaucoma and age-related macular degeneration, were observed. The cellular heterogeneity in the killifish retina was characterised, confirming the presence of all typical vertebrate retinal cell types. Data integration from age-matched samples between the bulk and scRNA-seq experiments revealed a loss of cellular specificity in gene expression upon ageing, suggesting potential disruption in transcriptional homeostasis. Differential expression analysis within the identified cell types highlighted the role of glial/immune cells as important stress regulators during ageing. Our work emphasises the value of the fast-ageing killifish in elucidating molecular signatures in age-associated retinal disease and vision decline. This study contributes to the understanding of how age-related changes in molecular pathways may impact CNS health, providing insights that may inform future therapeutic strategies for age-related pathologies.

Retinylidene chromophore hydrolysis from mammalian visual and non-visual opsins

Rhodopsin (Rho) and cone opsins are essential for detection of light. They respond via photoisomerization, converting their Schiff-base-adducted 11-cis-retinylidene chromophores to the all-trans configuration, eliciting conformational changes to activate opsin signaling. Subsequent Schiff-base hydrolysis releases all-trans-retinal, initiating two important cycles that maintain continuous vision-the Rho photocycle and visual cycle pathway. Schiff-base hydrolysis has been thoroughly studied with photoactivated Rho but not with cone opsins. Using established methodology, we directly measured the formation of Schiff-base between retinal chromophores with mammalian visual and nonvisual opsins of the eye. Next, we determined the rate of light-induced chromophore hydrolysis. We found that retinal hydrolysis from photoactivated cone opsins was markedly faster than from photoactivated Rho. Bovine retinal G protein-coupled receptor (bRGR) displayed rapid hydrolysis of its 11-cis-retinylidene photoproduct to quickly supply 11-cis-retinal and re-bind all-trans-retinal. Hydrolysis within bRGR in native retinal pigment epithelium microsomal membranes was >6-times faster than that of bRGR purified in detergent micelles. N-terminal-targeted antibodies significantly slowed bRGR hydrolysis, while C-terminal antibodies had no effect. Our study highlights the much faster photocycle of cone opsins relative to Rho and the crucial role of RGR in chromophore recycling in daylight. By contrast, in our experimental conditions, bovine peropsin did not form pigment in the presence of all-trans-retinal nor with any mono-cis retinal isomers, leaving uncertain the role of this opsin as a light sensor.

Toward Retinal Organoids in High-Throughput

Human retinal organoids recapitulate the cellular diversity, arrangement, gene expression, and functional aspects of the human retina. Protocols to generate human retinal organoids from pluripotent stem cells are typically labor intensive, include many manual handling steps, and the organoids need to be maintained for several months until they mature. To generate large numbers of human retinal organoids for therapy development and screening purposes, scaling up retinal organoid production, maintenance, and analysis is of utmost importance. In this review, we discuss strategies to increase the number of high-quality retinal organoids while reducing manual handling steps. We further review different approaches to analyze thousands of retinal organoids with currently available technologies and point to challenges that still await to be overcome both in culture and analysis of retinal organoids.

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.

Aberrant gene expression yet undiminished retinal ganglion cell genesis in iPSC-derived models of optic nerve hypoplasia

BACKGROUND

Optic nerve hypoplasia (ONH), the leading congenital cause of permanent blindness, is characterized by a retinal ganglion cell (RGC) deficit at birth. Multifactorial developmental events are hypothesized to underlie ONH and its frequently associated neurologic and endocrine abnormalities; however, environmental influences are unclear and genetic underpinnings are unexplored. This work investigates the genetic contribution to ONH RGC production and gene expression using patient induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs).

MATERIALS AND METHODS

iPSCs produced from ONH patients and controls were differentiated to ROs. RGC genesis was assessed using immunofluorescence and flow cytometry. Flow-sorted BRN3+ cells were collected for RNA extraction for RNA-Sequencing. Differential gene expression was assessed using DESeq2 and edgeR. PANTHER was employed to identify statistically over-represented ontologies among the differentially expressed genes (DEGs). DEGs of high interest to ONH were distinguished by assessing function, mutational constraint, and prior identification in ONH, autism and neurodevelopmental disorder (NDD) studies.

RESULTS

RGC genesis and survival were similar in ONH and control ROs. Differential expression of 70 genes was identified in both DESeq2 and edgeR analyses, representing a ~ 4-fold higher percentage of DEGs than in randomized study participants. DEGs showed trends towards over-representation of validated NDD genes and ONH exome variant genes. Among the DEGs, RAPGEF4 and DMD had the greatest number of disease-relevant features.

CONCLUSIONS

ONH genetic background was not associated with impaired RGC genesis but was associated with DEGs exhibiting disease contribution potential. This constitutes some of the first evidence of a genetic contribution to ONH.

RNA fusion in human retinal development

Chimeric RNAs have been found in both cancerous and healthy human cells. They have regulatory effects on human stem/progenitor cell differentiation, stemness maintenance, and central nervous system development. However, whether they are present in human retinal cells and their physiological functions in the retinal development remain unknown. Based on the human embryonic stem cell-derived retinal organoids (ROs) spanning from days 0 to 120, we present the expression atlas of chimeric RNAs throughout the developing ROs. We confirmed the existence of some common chimeric RNAs and also discovered many novel chimeric RNAs during retinal development. We focused on CTNNBIP1-CLSTN1 (CTCL) whose downregulation caused precocious neuronal differentiation and a marked reduction of neural progenitors in human cerebral organoids. CTCL is universally present in human retinas, ROs, and retinal cell lines, and its loss-of-function biases the progenitor cells toward retinal pigment epithelial cell fate at the expense of retinal cells. Together, this work provides a landscape of chimeric RNAs and reveals evidence for their critical role in human retinal development.

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.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes through Müller glia (MG) reprogramming and asymmetric cell division that produces a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, do MG reprogram to a developmental retinal progenitor cell (RPC) state? Second, to what extent does regeneration recapitulate retinal development? And finally, does loss of different retinal cell subtypes induce unique MG regeneration responses? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. Here we show that injury induces MG to reprogram to a state similar to late-stage RPCs. However, there are major transcriptional differences between MGPCs and RPCs, as well as major transcriptional differences between activated MG and MGPCs when different retinal cell subtypes are damaged. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes.

Pathway centric analysis for single-cell RNA-seq and spatial transcriptomics data with GSDensity

Advances in single-cell technology have enabled molecular dissection of heterogeneous biospecimens at unprecedented scales and resolutions. Cluster-centric approaches are widely applied in analyzing single-cell data, however they have limited power in dissecting and interpreting highly heterogenous, dynamically evolving data. Here, we present GSDensity, a graph-modeling approach that allows users to obtain pathway-centric interpretation and dissection of single-cell and spatial transcriptomics (ST) data without performing clustering. Using pathway gene sets, we show that GSDensity can accurately detect biologically distinct cells and reveal novel cell-pathway associations ignored by existing methods. Moreover, GSDensity, combined with trajectory analysis can identify curated pathways that are active at various stages of mouse brain development. Finally, GSDensity can identify spatially relevant pathways in mouse brains and human tumors including those following high-order organizational patterns in the ST data. Particularly, we create a pan-cancer ST map revealing spatially relevant and recurrently active pathways across six different tumor types.

ASCL1 induces neurogenesis in human Müller glia

In mammals, loss of retinal cells due to disease or trauma is an irreversible process that can lead to blindness. Interestingly, regeneration of retinal neurons is a well established process in some non-mammalian vertebrates and is driven by the Müller glia (MG), which are able to re-enter the cell cycle and reprogram into neurogenic progenitors upon retinal injury or disease. Progress has been made to restore this mechanism in mammals to promote retinal regeneration: MG can be stimulated to generate new neurons in vivo in the adult mouse retina after the over-expression of the pro-neural transcription factor Ascl1. In this study, we applied the same strategy to reprogram human MG derived from fetal retina and retinal organoids into neurons. Combining single cell RNA sequencing, single cell ATAC sequencing, immunofluorescence, and electrophysiology we demonstrate that human MG can be reprogrammed into neurogenic cells in vitro.

Comparative Analysis of Single-cell and Single-nucleus RNA-sequencing in a Rabbit Model of Retinal Detachment-related Proliferative Vitreoretinopathy

Purpose

Proliferative vitreoretinopathy (PVR) is the most common cause of failure of retinal reattachment surgery, and the molecular changes leading to this aberrant wound healing process are currently unknown. Our ultimate goal is to study PVR pathogenesis by employing single-cell transcriptomics to dissect cellular heterogeneity.

Design

Here we aimed to compare single-cell RNA sequencing (scRNA-seq)  and single-nucleus RNA-sequencing (snRNA-seq) of retinal PVR samples in the rabbit model.

Participants

Unilateral induction of PVR lesions in rabbit eyes with contralateral eyes serving as controls.

Methods

Proliferative vitreoretinopathy was induced unilaterally in Dutch Belted rabbits. At different timepoints after PVR induction, retinas were dissociated into either cells or nuclei suspension and processed for scRNA-seq or snRNA-seq.

Main Outcome Measures

Single cell and nuclei transcriptomic profiles of retinas after PVR induction.

Results

Single-cell RNA sequencing and snRNA-seq were conducted on retinas at 4 hours and 14 days after disease induction. Although the capture rate of unique molecular identifiers and genes were greater in scRNA-seq samples, overall gene expression profiles of individual cell types were highly correlated between scRNA-seq and snRNA-seq. A major disparity between the 2 sequencing modalities was the cell type capture rate, however, with glial cell types overrepresented in scRNA-seq, and inner retinal neurons were enriched by snRNA-seq. Furthermore, fibrotic Müller glia were overrepresented in snRNA-seq samples, whereas reactive Müller glia were overrepresented in scRNA-seq samples. Trajectory analyses were similar between the 2 methods, allowing for the combined analysis of the scRNA-seq and snRNA-seq data sets.

Conclusions

These findings highlight limitations of both scRNA-seq and snRNA-seq analysis and imply that use of both techniques together can more accurately identify transcriptional networks critical for aberrant fibrogenesis in PVR than using either in isolation.

Financial Disclosures

Proprietary or commercial disclosure may be found after the references.

Multimodal spatiotemporal phenotyping of human retinal organoid development

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

Centrins control chicken cone cell lipid droplet dynamics through lipid-droplet-localized SPDL1

As evolutionarily conserved organelles, lipid droplets (LDs) carry out numerous functions and have various subcellular localizations in different cell types and species. In avian cone cells, there is a single apically localized LD. We demonstrated that CIDEA (cell death inducing DFFA like effector a) and microtubules promote the formation of the single LD in chicken cone cells. Centrins, which are well-known centriole proteins, target to the cone cell LD via their C-terminal calcium-binding domains. Centrins localize on cone cell LDs with the help of SPDL1-L (spindle apparatus coiled-coil protein 1-L), a previously uncharacterized isoform of the kinetochore-associated dynein adaptor SPDL1. The loss of CETN3 or overexpression of a truncated CETN1 abrogates the apical localization of the cone cell LD. Simulation analysis showed that multiple LDs or a single mispositioned LD reduces the light sensitivity. Collectively, our findings identify a role of centrins in the regulation of cone cell LD localization, which is important for the light sensitivity of cone cells.

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.

Robust reprogramming of glia into neurons by inhibition of Notch signaling and NFI factors in adult mammalian retina

Generation of neurons through direct reprogramming has emerged as a promising therapeutic approach for neurodegenerative diseases. Despite successful applications in vitro , in vivo implementation has been hampered by low efficiency. In this study, we present a highly efficient strategy for reprogramming retinal glial cells into neurons by simultaneously inhibiting key negative regulators. By suppressing Notch signaling through the removal of its central mediator Rbpj, we induced mature Müller glial cells to reprogram into bipolar and amacrine neurons in uninjured adult mouse retinas, and observed that this effect was further enhanced by retinal injury. We found that specific loss of function of Notch1 and Notch2 receptors in Müller glia mimicked the effect of Rbpj deletion on Müller glia-derived neurogenesis. Integrated analysis of multiome (scRNA- and scATAC-seq) and CUT&Tag data revealed that Rbpj directly activates Notch effector genes and genes specific to mature Müller glia while also indirectly represses the expression of neurogenic bHLH factors. Furthermore, we found that combined loss of function of Rbpj and Nfia/b/x resulted in a robust conversion of nearly all Müller glia to neurons. Finally, we demonstrated that inducing Müller glial proliferation by AAV (adeno-associated virus)-mediated overexpression of dominant- active Yap supports efficient levels of Müller glia-derived neurogenesis in both Rbpj - and Nfia/b/x/Rbpj - deficient Müller glia. These findings demonstrate that, much like in zebrafish, Notch signaling actively represses neurogenic competence in mammalian Müller glia, and suggest that inhibition of Notch signaling and Nfia/b/x in combination with overexpression of activated Yap could serve as an effective component of regenerative therapies for degenerative retinal diseases.

Dynamic human retinal pigment epithelium (RPE) and choroid architecture based on single-cell transcriptomic landscape analysis

The retinal pigment epithelium (RPE) and choroid are located behind the human retina and have multiple functions in the human visual system. Knowledge of the RPE and choroid cells and their gene expression profiles are fundamental for understanding retinal disease mechanisms and therapeutic strategies. Here, we sequenced the RNA of about 0.3 million single cells from human RPE and choroids across two regions and seven ages, revealing regional and age differences within the human RPE and choroid. Cell-cell interactions highlight the broad connectivity networks between the RPE and different choroid cell types. Moreover, the transcription factors and their target genes change during aging. The coding of somatic variations increases during aging in the human RPE and choroid at the single-cell level. Moreover, we identified ELN as a candidate for improving RPE degeneration and choroidal structure during aging. The mapping of the molecular architecture of the human RPE and choroid improves our understanding of the human vision support system and offers potential insights into the intervention targets for retinal diseases.

Episodic live imaging of cone photoreceptor maturation in GNAT2-EGFP retinal organoids

Fluorescent reporter pluripotent stem cell-derived retinal organoids are powerful tools to investigate cell type-specific development and disease phenotypes. When combined with live imaging, they enable direct and repeated observation of cell behaviors within a developing retinal tissue. Here, we generated a human cone photoreceptor reporter line by CRISPR/Cas9 genome editing of WTC11-mTagRFPT-LMNB1 human induced pluripotent stem cells (iPSCs) by inserting enhanced green fluorescent protein (EGFP) coding sequences and a 2A self-cleaving peptide at the N-terminus of guanine nucleotide-binding protein subunit alpha transducin 2 (GNAT2). In retinal organoids generated from these iPSCs, the GNAT2-EGFP alleles robustly and exclusively labeled immature and mature cones. Episodic confocal live imaging of hydrogel immobilized retinal organoids allowed tracking of the morphological maturation of individual cones for >18 weeks and revealed inner segment accumulation of mitochondria and growth at 12.2 μm3 per day from day 126 to day 153. Immobilized GNAT2-EGFP cone reporter organoids provide a valuable tool for investigating human cone development and disease.

Identification of the immune-associated characteristics and predictive biomarkers of keratoconus based on single-cell RNA-sequencing and bulk RNA-sequencing

Background

Whether keratoconus (KC) is an inflammatory disease is currently debated. Hence, we aimed to investigate the immune-related features of KC based on single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (bulk RNA-seq) data.

Methods

scRNA-seq data were obtained from the Genome Sequence Archive (GSA), bulk RNA-seq data were obtained from the Gene Expression Omnibus (GEO), and immune-associated genes(IAGs) were obtained from the ImmPort database. Cell clusters of KC were annotated, and different cell clusters were then selected. The IAG score of each cell was calculated using the AUCell package. Three bulk RNA-seq datasets were merged and used to identify the differentially expressed genes (DEGs), biological functions, and immune characteristics. Weighted gene coexpression network analysis (WGCNA) was used to select the IAG score-related hub genes. Based on scRNA-seq and bulk RNA-seq analyses, three machine learning algorithms, including random forest (RF), support vector machine (SVM), and least absolute shrinkage and selection operator (LASSO) regression analysis, were used to identify potential prognostic markers for KC. A predictive nomogram was developed based on prognostic markers.

Results

Six cell clusters were identified in KC, and decreased corneal stromal cell-5 (CSC-5) and increased CSC-6 were found in KC. CSC and immune cell clusters had the highest IAG scores. The bulk RNA-seq analysis identified 1362 DEGs (553 upregulated and 809 downregulated) in KC. We found different immune cell populations and differentially expressed cytokines in KC. More than three key IAG score-related modules and 367 genes were identified. By integrating the scRNA-seq and bulk RNA-seq analyses, 250 IAGs were selected and then incorporated into three machine learning models, and 10 IAGs (CEP112, FYN, IFITM1, IGFBP5, LPIN2, MAP1B, RNASE1, RUNX3, SMIM10, and SRGN) were identified as potential prognostic genes that were significantly associated with cytokine and matrix metalloproteinase(MMP)1-14 expression. Finally, a predictive nomogram was constructed and validated.

Conclusion

Taken together, our results identified CSCs and immune cell clusters that may play a key role during KC progression by regulating immunological features and maintaining cell stability.

Molecular mechanisms controlling vertebrate retinal patterning, neurogenesis, and cell fate specification

This review covers recent advances in understanding the molecular mechanisms controlling neurogenesis and specification of the developing retina, with a focus on insights obtained from comparative single cell multiomic analysis. We discuss recent advances in understanding the mechanisms by which extrinsic factors trigger transcriptional changes that spatially pattern the optic cup (OC) and control the initiation and progression of retinal neurogenesis. We also discuss progress in unraveling the core evolutionarily conserved gene regulatory networks (GRNs) that specify early- and late-state retinal progenitor cells (RPCs) and neurogenic progenitors and that control the final steps in determining cell identity. Finally, we discuss findings that provide insight into regulation of species-specific aspects of retinal patterning and neurogenesis, including consideration of key outstanding questions in the field.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

Human retinal ganglion cell neurons generated by synchronous BMP inhibition and transcription factor mediated reprogramming

In optic neuropathies, including glaucoma, retinal ganglion cells (RGCs) die. Cell transplantation and endogenous regeneration offer strategies for retinal repair, however, developmental programs required for this to succeed are incompletely understood. To address this, we explored cellular reprogramming with transcription factor (TF) regulators of RGC development which were integrated into human pluripotent stem cells (PSCs) as inducible gene cassettes. When the pioneer factor NEUROG2 was combined with RGC-expressed TFs (ATOH7, ISL1, and POU4F2) some conversion was observed and when pre-patterned by BMP inhibition, RGC-like induced neurons (RGC-iNs) were generated with high efficiency in just under a week. These exhibited transcriptional profiles that were reminiscent of RGCs and exhibited electrophysiological properties, including AMPA-mediated synaptic transmission. Additionally, we demonstrated that small molecule inhibitors of DLK/LZK and GCK-IV can block neuronal death in two pharmacological axon injury models. Combining developmental patterning with RGC-specific TFs thus provided valuable insight into strategies for cell replacement and neuroprotection.

Osteopontin drives retinal ganglion cell resiliency in glaucomatous optic neuropathy

Chronic neurodegeneration and acute injuries lead to neuron losses via diverse processes. We compared retinal ganglion cell (RGC) responses between chronic glaucomatous conditions and the acute injury model. Among major RGC subclasses, αRGCs and intrinsically photosensitive RGCs (ipRGCs) preferentially survive glaucomatous conditions, similar to findings in the retina subject to axotomy. Focusing on an αRGC intrinsic factor, Osteopontin (secreted phosphoprotein 1 [Spp1]), we found an ectopic neuronal expression of Osteopontin (Spp1) in other RGCs subject to glaucomatous conditions. This contrasted with the Spp1 downregulation subject to axotomy. αRGC-specific Spp1 elimination led to significant αRGC loss, diminishing their resiliency. Spp1 overexpression led to robust neuroprotection of susceptible RGC subclasses under glaucomatous conditions. In contrast, Spp1 overexpression did not significantly protect RGCs subject to axotomy. Additionally, SPP1 marked adult human RGC subsets with large somata and SPP1 expression in the aqueous humor correlated with glaucoma severity. Our study reveals Spp1's role in mediating neuronal resiliency in glaucoma.

Dopamine signaling from ganglion cells directs layer-specific angiogenesis in the retina

During central nervous system (CNS) development, a precisely patterned vasculature emerges to support CNS function. How neurons control angiogenesis is not well understood. Here, we show that the neuromodulator dopamine restricts vascular development in the retina via temporally limited production by an unexpected neuron subset. Our genetic and pharmacological experiments demonstrate that elevating dopamine levels inhibits tip-cell sprouting and vessel growth, whereas reducing dopamine production by all retina neurons increases growth. Dopamine production by canonical dopaminergic amacrine interneurons is dispensable for these events. Instead, we found that temporally restricted dopamine production by retinal ganglion cells (RGCs) modulates vascular development. RGCs produce dopamine precisely during angiogenic periods. Genetically limiting dopamine production by ganglion cells, but not amacrines, decreases angiogenesis. Conversely, elevating ganglion-cell-derived dopamine production inhibits early vessel growth. These vasculature outcomes occur downstream of vascular endothelial growth factor receptor (VEGFR) activation and Notch-Jagged1 signaling. Jagged1 is increased and subsequently inhibits Notch signaling when ganglion cell dopamine production is reduced. Our findings demonstrate that dopaminergic neural activity from a small neuron subset functions upstream of VEGFR to serve as developmental timing cue that regulates vessel growth.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes. This regeneration requires Müller glia (MG) to reprogram and divide asymmetrically to produce a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, does loss of different retinal cell subtypes induce unique MG regeneration responses? Second, do MG reprogram to a developmental retinal progenitor cell state? And finally, to what extent does regeneration recapitulate retinal development? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. While MG reprogram to a state similar to late-stage retinal progenitors in developing retinas, there are transcriptional differences between reprogrammed MG/MGPCs and late progenitors, as well as reprogrammed MG in outer and inner retinal damage models. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes. This work identifies major differences between gene regulatory networks activated following the selective loss of different subtypes of retina neurons, as well as between retinal regeneration and development.

Self-formation of concentric zones of telencephalic and ocular tissues and directional retinal ganglion cell axons

The telencephalon and eye in mammals are originated from adjacent fields at the anterior neural plate. Morphogenesis of these fields generates telencephalon, optic-stalk, optic-disc, and neuroretina along a spatial axis. How these telencephalic and ocular tissues are specified coordinately to ensure directional retinal ganglion cell (RGC) axon growth is unclear. Here, we report self-formation of human telencephalon-eye organoids comprising concentric zones of telencephalic, optic-stalk, optic-disc, and neuroretinal tissues along the center-periphery axis. Initially-differentiated RGCs grew axons towards and then along a path defined by adjacent PAX2+ VSX2+ optic-disc cells. Single-cell RNA sequencing of these organoids not only confirmed telencephalic and ocular identities but also identified expression signatures of early optic-disc, optic-stalk, and RGCs. These signatures were similar to those in human fetal retinas. Optic-disc cells in these organoids differentially expressed FGF8 and FGF9; FGFR inhibitions drastically decreased early RGC differentiation and directional axon growth. Through the RGC-specific cell-surface marker CNTN2 identified here, electrophysiologically excitable RGCs were isolated under a native condition. Our findings provide insight into the coordinated specification of early telencephalic and ocular tissues in humans and establish resources for studying RGC-related diseases such as glaucoma.

Establishing Functional Retina in a Dish: Progress and Promises of Induced Pluripotent Stem Cell-Based Retinal Neuron Differentiation

The human eye plays a critical role in vision perception, but various retinal degenerative diseases such as retinitis pigmentosa (RP), glaucoma, and age-related macular degeneration (AMD) can lead to vision loss or blindness. Although progress has been made in understanding retinal development and in clinical research, current treatments remain inadequate for curing or reversing these degenerative conditions. Animal models have limited relevance to humans, and obtaining human eye tissue samples is challenging due to ethical and legal considerations. Consequently, researchers have turned to stem cell-based approaches, specifically induced pluripotent stem cells (iPSCs), to generate distinct retinal cell populations and develop cell replacement therapies. iPSCs offer a novel platform for studying the key stages of human retinogenesis and disease-specific mechanisms. Stem cell technology has facilitated the production of diverse retinal cell types, including retinal ganglion cells (RGCs) and photoreceptors, and the development of retinal organoids has emerged as a valuable in vitro tool for investigating retinal neuron differentiation and modeling retinal diseases. This review focuses on the protocols, culture conditions, and techniques employed in differentiating retinal neurons from iPSCs. Furthermore, it emphasizes the significance of molecular and functional validation of the differentiated cells.

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.

Atonal homolog 7 (ATOH7) confers neuroprotection for photoreceptor cells in glaucoma via inhibition of the notch pathway

Single-cell RNA sequencing (scRNA-seq) technologies enable the profiling and analysis of the transcriptomes of single cells and hold promise for clarifying gene mechanisms at single-cell resolution. We based this study on scRNA-seq data to reveal glaucoma-related genes and downstream pathways with neuroprotection effects. The scRNA-seq datasets related to glaucoma of retinal tissue samples of human beings and Atonal Homolog 7 (ATOH7)-null mice were obtained from the GEO database. The 74 top marker genes and 20 cell clusters were obtained in human retinal tissue samples. The key gene ATOH7 was found after the intersection with genes from GeneCards data. In the ATOH7-null mouse retinal tissue samples, pseudotime inference demonstrated significant changes in cell differentiation. Moreover, mouse retinal photoreceptor cells (PRCs) were cultured and treated with lentivirus carrying oe-ATOH7 alone or in combination with Notch signaling pathway activator Jagged-1/FC, after which cell biological functions were determined. The involvement of ATOH7 in glaucoma was identified through regulating PRCs. Furthermore, ATOH7 conferred neuroprotection in PRCs in glaucoma by mediating the Notch signaling pathway. In vitro data confirmed that ATOH7 overexpression promoted the differentiation of PRCs and inhibited their apoptosis by suppressing the Notch signaling pathway. The evidence provided by our study highlighted the involvement of ATOH7 in the blockade of the Notch signaling pathway, resulting in the neuroprotection for PRCs in glaucoma.

Developmental adaptation of rod photoreceptor number via photoreception in melanopsin (OPN4) retinal ganglion cells

Photoreception, a form of sensory experience, is essential for normal development of the mammalian visual system. Detecting photons during development is a prerequisite for visual system function - from vision's first synapse at the cone pedicle and maturation of retinal vascular networks, to transcriptional establishment and maturation of cell types within the visual cortex. Consistent with this theme, we find that the lighting environment regulates developmental rod photoreceptor apoptosis via OPN4-expressing intrinsically photosensitive retinal ganglion cells (ipRGCs). Using a combination of genetics, sensory environment manipulations, and computational approaches, we establish a molecular pathway in which light-dependent glutamate release from ipRGCs is detected via a transiently expressed kainate receptor (GRIK3) in immature rods localized to the inner retina. Communication between ipRGCs and nascent inner retinal rods appears to be mediated by unusual hybrid neurites projecting from ipRGCs that sense light before eye-opening. These structures, previously referred to as outer retinal dendrites (ORDs), span the ipRGC-immature rod distance over the first postnatal week and contain the machinery for sensory detection (melanopsin, OPN4) and axonal/anterograde neurotransmitter release (Synaptophysin, and VGLUT2). Histological and computational assessment of human mid-gestation development reveal conservation of several hallmarks of an ipRGC-to-immature rod pathway, including displaced immature rods, transient GRIK3 expression in the rod lineage, and the presence of ipRGCs with putative neurites projecting deep into the developing retina. Thus, this analysis defines a retinal retrograde signaling pathway that links the sensory environment to immature rods via ipRGC photoreceptors, allowing the visual system to adapt to distinct lighting environments priory to eye-opening.

Integration of human stem cell-derived in vitro systems and mouse preclinical models identifies complex pathophysiologic mechanisms in retinal dystrophy

Rare DRAM2 coding variants cause retinal dystrophy with early macular involvement via unknown mechanisms. We found that DRAM2 is ubiquitously expressed in the human eye and expression changes were observed in eyes with more common maculopathy such as Age-related Macular Degeneration (AMD). To gain insights into pathogenicity of DRAM2-related retinopathy, we used a combination of in vitro and in vivo models. We found that DRAM2 loss in human pluripotent stem cell (hPSC)-derived retinal organoids caused the presence of additional mesenchymal cells. Interestingly, Dram2 loss in mice also caused increased proliferation of cells from the choroid in vitro and exacerbated choroidal neovascular lesions in vivo. Furthermore, we observed that DRAM2 loss in human retinal pigment epithelial (RPE) cells resulted in increased susceptibility to stress-induced cell death in vitro and that Dram2 loss in mice caused age-related photoreceptor degeneration. This highlights the complexity of DRAM2 function, as its loss in choroidal cells provided a proliferative advantage, whereas its loss in post-mitotic cells, such as photoreceptor and RPE cells, increased degeneration susceptibility. Different models such as human pluripotent stem cell-derived systems and mice can be leveraged to study and model human retinal dystrophies; however, cell type and species-specific expression must be taken into account when selecting relevant systems.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes. This regeneration requires Müller glia (MG) to reprogram and divide asymmetrically to produce a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, does loss of different retinal cell subtypes induce unique MG regeneration responses? Second, do MG reprogram to a developmental retinal progenitor cell state? And finally, to what extent does regeneration recapitulate retinal development? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. While MG reprogram to a state similar to late-stage retinal progenitors in developing retinas, there are transcriptional differences between reprogrammed MG/MGPCs and late progenitors, as well as reprogrammed MG in outer and inner retinal damage models. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes. This work identifies major differences between gene regulatory networks activated following the selective loss of different subtypes of retina neurons, as well as between retinal regeneration and development.

The efficient induction of human retinal ganglion-like cells provides a platform for studying optic neuropathies

Retinal ganglion cells (RGCs) are essential for vision perception. In glaucoma and other optic neuropathies, RGCs and their optic axons undergo degenerative change and cell death; this can result in irreversible vision loss. Here we developed a rapid protocol for directly inducing RGC differentiation from human induced pluripotent stem cells (hiPSCs) by the overexpression of ATOH7, BRN3B, and SOX4. The hiPSC-derived RGC-like cells (iRGCs) show robust expression of various RGC-specific markers by whole transcriptome profiling. A functional assessment was also carried out and this demonstrated that these iRGCs display stimulus-induced neuronal activity, as well as spontaneous neuronal activity. Ethambutol (EMB), an effective first-line anti-tuberculosis agent, is known to cause serious visual impairment and irreversible vision loss due to the RGC degeneration in a significant number of treated patients. Using our iRGCs, EMB was found to induce significant dose-dependent and time-dependent increases in cell death and neurite degeneration. Western blot analysis revealed that the expression levels of p62 and LC3-II were upregulated, and further investigations revealed that EMB caused a blockade of lysosome-autophagosome fusion; this indicates that impairment of autophagic flux is one of the adverse effects of that EMB has on iRGCs. In addition, EMB was found to elevate intracellular reactive oxygen species (ROS) levels increasing apoptotic cell death. This could be partially rescued by the co-treatment with the ROS scavenger NAC. Taken together, our findings suggest that this iRGC model, which achieves both high yield and high purity, is suitable for investigating optic neuropathies, as well as being useful when searching for potential drugs for therapeutic treatment and/or disease prevention.

Cell-type specification in the retina: Recent discoveries from transcriptomic approaches

Understanding the formation of the complex nervous system hinges on decoding the mechanism that specifies a vast array of neuronal types, each endowed with a unique morphology, physiology, and connectivity. As a pivotal step towards addressing this problem, seminal work has been devoted to characterizing distinct neuronal types. In recent years, high-throughput, single-cell transcriptomic methods have enabled a rapid inventory of cell types in various regions of the nervous system, with the retina exhibiting complete molecular characterization across many vertebrate species. This invaluable resource has furnished a fresh perspective for investigating the molecular principles of cell-type specification, thereby advancing our understanding of retinal development. Accordingly, this review focuses on the most recent transcriptomic characterizations of retinal cells, with a particular focus on amacrine cells and retinal ganglion cells. These investigations have unearthed new insights into their cell-type specification.

Inducible Rbpms-CreERT2 Mouse Line for Studying Gene Function in Retinal Ganglion Cell Physiology and Disease

Retinal ganglion cells (RGCs) are the sole output neurons conveying visual stimuli from the retina to the brain, and dysfunction or loss of RGCs is the primary determinant of visual loss in traumatic and degenerative ocular conditions. Currently, there is a lack of RGC-specific Cre mouse lines that serve as invaluable tools for manipulating genes in RGCs and studying the genetic basis of RGC diseases. The RNA-binding protein with multiple splicing (RBPMS) is identified as the specific marker of all RGCs. Here, we report the generation and characterization of a knock-in mouse line in which a P2A-CreERT2 coding sequence is fused in-frame to the C-terminus of endogenous RBPMS, allowing for the co-expression of RBPMS and CreERT2. The inducible Rbpms-CreERT2 mice exhibited a high recombination efficiency in activating the expression of the tdTomato reporter gene in nearly all adult RGCs as well as in differentiated RGCs starting at E13.5. Additionally, both heterozygous and homozygous Rbpms-CreERT2 knock-in mice showed no detectable defect in the retinal structure, visual function, and transcriptome. Together, these results demonstrated that the Rbpms-CreERT2 knock-in mouse can serve as a powerful and highly desired genetic tool for lineage tracing, genetic manipulation, retinal physiology study, and ocular disease modeling in RGCs.

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.

Application of Human Stem Cell Derived Retinal Organoids in the Exploration of the Mechanisms of Early Retinal Development

The intricate neural circuit of retina extracts salient features of the natural world and forms bioelectric impulse as the origin of vision. The early development of retina is a highly complex and coordinated process in morphogenesis and neurogenesis. Increasing evidence indicates that stem cells derived human retinal organoids (hROs) in vitro faithfully recapitulates the embryonic developmental process of human retina no matter in the transcriptome, cellular biology and histomorphology. The emergence of hROs greatly deepens on the understanding of early development of human retina. Here, we reviewed the events of early retinal development both in animal embryos and hROs studies, which mainly comprises the formation of optic vesicle and optic cup shape, differentiation of retinal ganglion cells (RGCs), photoreceptor cells (PRs) and its supportive retinal pigment epithelium cells (RPE). We also discussed the classic and frontier molecular pathways up to date to decipher the underlying mechanisms of early development of human retina and hROs. Finally, we summarized the application prospect, challenges and cutting-edge techniques of hROs for uncovering the principles and mechanisms of retinal development and related developmental disorder. hROs is a priori selection for studying human retinal development and function and may be a fundamental tool for unlocking the unknown insight into retinal development and disease.