Generating and Using Transcriptomically Based Retinal Cell Atlases

High temporal frequency light response in mouse retina requires FAT3 signaling in bipolar cells

Vision is initiated by the reception of light by photoreceptors and subsequent processing via downstream retinal neurons. Proper cellular organization depends on the multi-functional tissue polarity protein FAT3, which is required for amacrine cell connectivity and retinal lamination. Here we investigated the retinal function of Fat3 mutant mice and found decreases in physiological and perceptual responses to high frequency flashes. These defects did not correlate with abnormal amacrine cell wiring, pointing instead to a role in bipolar cell subtypes that also express FAT3. The role of FAT3 in the response to high temporal frequency flashes depends upon its ability to transduce an intracellular signal. Mechanistically, FAT3 binds to the synaptic protein PTPσ, intracellularly, and is required to localize GRIK1 to OFF-cone bipolar cell synapses with cone photoreceptors. These findings expand the repertoire of FAT3 functions and reveal its importance in bipolar cells for high frequency light response.

Comprehensive single-cell atlas of the mouse retina

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of cellular heterogeneity by characterizing cell types across tissues and species. While several mouse retinal scRNA-seq datasets exist, each dataset is either limited in cell numbers or focused on specific cell classes, thereby hindering comprehensive gene expression analysis across all retina types. To fill the gap, we generated the largest retinal scRNA-seq dataset to date, comprising approximately 190,000 single cells from C57BL/6J mouse retinas, enriched for rare population cells via antibody-based magnetic cell sorting. Integrating this dataset with public datasets, we constructed the Mouse Retina Cell Atlas (MRCA) for wild-type mice, encompassing over 330,000 cells, characterizing 12 major classes and 138 cell types. The MRCA consolidates existing knowledge, identifies new cell types, and is publicly accessible via CELLxGENE, UCSC Cell Browser, and the Broad Single Cell Portal, providing a user-friendly resource for the mouse retina research community.

Optimizing the design of spatial genomic studies

Spatial genomic technologies characterize the relationship between the structural organization of cells and their cellular state. Despite the availability of various spatial transcriptomic and proteomic profiling platforms, these experiments remain costly and labor-intensive. Traditionally, tissue slicing for spatial sequencing involves parallel axis-aligned sections, often yielding redundant or correlated information. We propose structured batch experimental design, a method that improves the cost efficiency of spatial genomics experiments by profiling tissue slices that are maximally informative, while recognizing the destructive nature of the process. Applied to two spatial genomics studies-one to construct a spatially-resolved genomic atlas of a tissue and another to localize a region of interest in a tissue, such as a tumor-our approach collects more informative samples using fewer slices compared to traditional slicing strategies. This methodology offers a foundation for developing robust and cost-efficient design strategies, allowing spatial genomics studies to be deployed by smaller, resource-constrained labs.

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.

Identifying neuron types and circuit mechanisms in the auditory midbrain

The inferior colliculus (IC) is a critical computational hub in the central auditory pathway. From its position in the midbrain, the IC receives nearly all the ascending output from the lower auditory brainstem and provides the main source of auditory information to the thalamocortical system. In addition to being a crossroads for auditory circuits, the IC is rich with local circuits and contains more than five times as many neurons as the nuclei of the lower auditory brainstem combined. These results hint at the enormous computational power of the IC, and indeed, systems-level studies have identified numerous important transformations in sound coding that occur in the IC. However, despite decades of effort, the cellular mechanisms underlying IC computations and how these computations change following hearing loss have remained largely impenetrable. In this review, we argue that this challenge persists due to the surprisingly difficult problem of identifying the neuron types and circuit motifs that comprise the IC. After summarizing the extensive evidence pointing to a diversity of neuron types in the IC, we highlight the successes of recent efforts to parse this complexity using molecular markers to define neuron types. We conclude by arguing that the discovery of molecularly identifiable neuron types ushers in a new era for IC research marked by molecularly targeted recordings and manipulations. We propose that the ability to reproducibly investigate IC circuits at the neuronal level will lead to rapid advances in understanding the fundamental mechanisms driving IC computations and how these mechanisms shift following hearing loss.

Comprehensive single-cell atlas of the mouse retina

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of cellular heterogeneity at the single-cell resolution by classifying and characterizing cell types in multiple tissues and species. While several mouse retinal scRNA-seq reference datasets have been published, each dataset either has a relatively small number of cells or is focused on specific cell classes, and thus is suboptimal for assessing gene expression patterns across all retina types at the same time. To establish a unified and comprehensive reference for the mouse retina, we first generated the largest retinal scRNA-seq dataset to date, comprising approximately 190,000 single cells from C57BL/6J mouse whole retinas. This dataset was generated through the targeted enrichment of rare population cells via antibody-based magnetic cell sorting. By integrating this new dataset with public datasets, we conducted an integrated analysis to construct the Mouse Retina Cell Atlas (MRCA) for wild-type mice, which encompasses over 330,000 single cells. The MRCA characterizes 12 major classes and 138 cell types. It captured consensus cell type characterization from public datasets and identified additional new cell types. To facilitate the public use of the MRCA, we have deposited it in CELLxGENE, UCSC Cell Browser, and the Broad Single Cell Portal for visualization and gene expression exploration. The comprehensive MRCA serves as an easy-to-use, one-stop data resource for the mouse retina communities.

Population encoding of stimulus features along the visual hierarchy

The retina and primary visual cortex (V1) both exhibit diverse neural populations sensitive to diverse visual features. Yet it remains unclear how neural populations in each area partition stimulus space to span these features. One possibility is that neural populations are organized into discrete groups of neurons, with each group signaling a particular constellation of features. Alternatively, neurons could be continuously distributed across feature-encoding space. To distinguish these possibilities, we presented a battery of visual stimuli to the mouse retina and V1 while measuring neural responses with multi-electrode arrays. Using machine learning approaches, we developed a manifold embedding technique that captures how neural populations partition feature space and how visual responses correlate with physiological and anatomical properties of individual neurons. We show that retinal populations discretely encode features, while V1 populations provide a more continuous representation. Applying the same analysis approach to convolutional neural networks that model visual processing, we demonstrate that they partition features much more similarly to the retina, indicating they are more like big retinas than little brains.

Retinal neurons establish mosaic patterning by excluding homotypic somata from their dendritic territory

In vertebrate retina, individual neurons of the same type are distributed regularly across the tissue in a pattern known as a mosaic. Establishment of mosaics during development requires cell-cell repulsion among homotypic neurons, but the mechanisms underlying this repulsion remain unknown. Here we show that two mouse retinal cell types, OFF and ON starburst amacrine cells, establish mosaic spacing by using their dendritic arbors to repel neighboring homotypic somata. Using newly-generated transgenic tools and single cell labeling, we identify a transient developmental period when starburst somata receive extensive contacts from neighboring starburst dendrites; these serve to exclude somata from settling within the neighbor's dendritic territory. Dendrite-soma exclusion is mediated by MEGF10, a cell-surface molecule required for starburst mosaic patterning. Our results implicate dendrite-soma exclusion as a key mechanism underlying starburst mosaic spacing, and suggest that this could be a general mechanism for mosaic patterning across many cell types and species.

A simple, rapid fluorescent reporter-based method for detection of ectopic cre recombinase expression in presumed retinal cell type-targeted mouse lines

Although cell type-specific Cre recombinase-expressing mouse lines are commonly used to generate conditional knockout of genes of interest, germline recombination and ectopic "leakiness" in Cre recombinase expression in non-specific cell types has been observed in several neuronal and glial-specific Cre lines. This often leads to inadvertent loss of conditional mouse lines, requiring rederivation. It is therefore imperative to be able to monitor and validate cell type-specific Cre recombinase-mediated gene editing. Herein, we describe a simple, inexpensive, rapid ZsGreen fluor-reporter-based strategy for genotype-free identification of ectopic leakiness using a custom-designed, 3-D blue LED light box. We assessed cell type-specific expression in several allegedly specific Cre recombinase mouse lines commonly used in vision research: retinal pigment epithelium (RPE)-specific (VMD2 (Best1) Cre, RPE65 Cre); astrocyte-specific (GFAP Cre); as well as photoreceptor-bipolar progenitor cell-specific (CRX Cre). Our standardized workflow allows facile, rapid identification of ectopic and non-specific Cre recombinase expression in any presume specific Cre mouse line, without the need for genotyping and without causing animal distress.

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.

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.

Pou4f1-Tbr1 transcriptional cascade controls the formation of Jam2-expressing retinal ganglion cells

More than 40 retinal ganglion cell (RGC) subtypes have been categorized in mouse based on their morphologies, functions, and molecular features. Among these diverse subtypes, orientation-selective Jam2-expressing RGCs (J-RGCs) has two unique morphologic characteristics: the ventral-facing dendritic arbor and the OFF-sublaminae stratified terminal dendrites in the inner plexiform layer. Previously, we have discovered that T-box transcription factor T-brain 1 (Tbr1) is expressed in J-RGCs. We further found that Tbr1 is essential for the expression of Jam2, and Tbr1 regulates the formation and the dendritic morphogenesis of J-RGCs. However, Tbr1 begins to express in terminally differentiated RGCs around perinatal stage, suggesting that it is unlikely involved in the initial fate determination for J-RGC and other upstream transcription factors must control Tbr1 expression and J-RGC formation. Using the Cleavage Under Targets and Tagmentation technique, we discovered that Pou4f1 binds to Tbr1 on the evolutionary conserved exon 6 and an intergenic region downstream of the 3'UTR, and on a region flanking the promoter and the first exon of Jam2. We showed that Pou4f1 is required for the expression of Tbr1 and Jam2, indicating Pou4f1 as a direct upstream regulator of Tbr1 and Jam2. Most interestingly, the Pou4f1-bound element in exon 6 of Tbr1 possesses high-level enhancer activity, capable of directing reporter gene expression in J-RGCs. Together, these data revealed a Pou4f1-Tbr1-Jam2 genetic hierarchy as a critical pathway in the formation of J-RGC subtype.

Optimizing the design of spatial genomic studies

Spatially-resolved genomic technologies have shown promise for studying the relationship between the structural arrangement of cells and their functional behavior. While numerous sequencing and imaging platforms exist for performing spatial transcriptomics and spatial proteomics profiling, these experiments remain expensive and labor-intensive. Thus, when performing spatial genomics experiments using multiple tissue slices, there is a need to select the tissue cross sections that will be maximally informative for the purposes of the experiment. In this work, we formalize the problem of experimental design for spatial genomics experiments, which we generalize into a problem class that we call structured batch experimental design. We propose approaches for optimizing these designs in two types of spatial genomics studies: one in which the goal is to construct a spatially-resolved genomic atlas of a tissue and another in which the goal is to localize a region of interest in a tissue, such as a tumor. We demonstrate the utility of these optimal designs, where each slice is a two-dimensional plane, on several spatial genomics datasets.

Retinal TRP channels: Cell-type-specific regulators of retinal homeostasis and multimodal integration

Transient receptor potential (TRP) channels are a widely expressed family of 28 evolutionarily conserved cationic ion channels that operate as primary detectors of chemical and physical stimuli and secondary effectors of metabotropic and ionotropic receptors. In vertebrates, the channels are grouped into six related families: TRPC, TRPV, TRPM, TRPA, TRPML, and TRPP. As sensory transducers, TRP channels are ubiquitously expressed across the body and the CNS, mediating critical functions in mechanosensation, nociception, chemosensing, thermosensing, and phototransduction. This article surveys current knowledge about the expression and function of the TRP family in vertebrate retinas, which, while dedicated to transduction and transmission of visual information, are highly susceptible to non-visual stimuli. Every retinal cell expresses multiple TRP subunits, with recent evidence establishing their critical roles in paradigmatic aspects of vertebrate vision that include TRPM1-dependent transduction of ON bipolar signaling, TRPC6/7-mediated ganglion cell phototransduction, TRP/TRPL phototransduction in Drosophila and TRPV4-dependent osmoregulation, mechanotransduction, and regulation of inner and outer blood-retina barriers. TRP channels tune light-dependent and independent functions of retinal circuits by modulating the intracellular concentration of the 2nd messenger calcium, with emerging evidence implicating specific subunits in the pathogenesis of debilitating diseases such as glaucoma, ocular trauma, diabetic retinopathy, and ischemia. Elucidation of TRP channel involvement in retinal biology will yield rewards in terms of fundamental understanding of vertebrate vision and therapeutic targeting to treat diseases caused by channel dysfunction or over-activation.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

Intrinsic heterogeneity in axon regeneration

The nervous system is composed of a variety of neurons and glial cells with different morphology and functions. In the mammalian peripheral nervous system (PNS) or the lower vertebrate central nervous system (CNS), most neurons can regenerate extensively after axotomy, while the neurons in the mammalian CNS possess only limited regenerative ability. This heterogeneity is common within and across species. The studies about the transcriptomes after nerve injury in different animal models have revealed a series of molecular and cellular events that occurred in neurons after axotomy. However, responses of various types of neurons located in different positions of individuals were different remarkably. Thus, researchers aim to find the key factors that are conducive to regeneration, so as to provide the molecular basis for solving the regeneration difficulties after CNS injury. Here we review the heterogeneity of axonal regeneration among different cell subtypes in different animal models or the same organ, emphasizing the importance of comparative studies within and across species.

Homeostatic plasticity in the retina

Vision begins in the retina, whose intricate neural circuits extract salient features of the environment from the light entering our eyes. Neurodegenerative diseases of the retina (e.g., inherited retinal degenerations, age-related macular degeneration, and glaucoma) impair vision and cause blindness in a growing number of people worldwide. Increasing evidence indicates that homeostatic plasticity (i.e., the drive of a neural system to stabilize its function) can, in principle, preserve retinal function in the face of major perturbations, including neurodegeneration. Here, we review the circumstances and events that trigger homeostatic plasticity in the retina during development, sensory experience, and disease. We discuss the diverse mechanisms that cooperate to compensate and the set points and outcomes that homeostatic retinal plasticity stabilizes. Finally, we summarize the opportunities and challenges for unlocking the therapeutic potential of homeostatic plasticity. Homeostatic plasticity is fundamental to understanding retinal development and function and could be an important tool in the fight to preserve and restore vision.

Developmental errors in the common marmoset retina

Although retinal organization is remarkably conserved, morphological anomalies can be found to different extents and varieties across animal species with each presenting unique characteristics and patterns of displaced and misplaced neurons. One of the most widely used non-human primates in research, the common marmoset (Callithrix jaccus) could potentially also be of interest for visual research, but is unfortunately not well characterized in this regard. Therefore, the aim of our study was to provide a first time description of structural retinal layering including morphological differences and distinctive features in this species. Retinas from animals (n = 26) of both sexes and different ages were immunostained with cell specific antibodies to label a variety of bipolar, amacrine and ganglion cells. Misplaced ganglion cells with somata in the outermost part of the inner nuclear layer and rod bipolar cells with axon terminals projecting into the outer plexiform layer instead of the inner plexiform layer independent of age or sex of the animals were the most obvious findings, whereas misplaced amacrine cells and misplaced cone bipolar axon terminals occurred to a lesser extent. With this first time description of developmental retinal errors over a wide age range, we provide a basic characterization of the retinal system of the common marmosets, which can be taken into account for future studies in this and other animal species. The finding of misplaced ganglion cells and misplaced bipolar cell axon terminals was not reported before and displays an anatomic variation worthwhile for future analyzes of their physiological and functional impact.

Sensitivity to extracellular potassium underlies type-intrinsic differences in retinal ganglion cell excitability

Neuronal type-specific physiologic heterogeneity can be driven by both extrinsic and intrinsic mechanisms. In retinal ganglion cells (RGCs), which carry visual information from the retina to central targets, evidence suggests intrinsic properties shaping action potential (AP) generation significantly impact the responses of RGCs to visual stimuli. Here, we explored how differences in intrinsic excitability further distinguish two RCG types with distinct presynaptic circuits, alpha ON-sustained (αON-S) cells and alpha OFF-sustained (αOFF-S) cells. We found that αOFF-S RGCs are more excitable to modest depolarizing currents than αON-S RGCs but excitability plateaued earlier as depolarization increased (i.e., depolarization block). In addition to differences in depolarization block sensitivity, the two cell types also produced distinct AP shapes with increasing stimulation. αOFF-S AP width and variability increased with depolarization magnitude, which correlated with the onset of depolarization block, while αON-S AP width and variability remained stable. We then tested if differences in depolarization block observed in αON-S and αOFF-S RGCs were due to sensitivity to extracellular potassium. We found αOFF-S RGCs more sensitive to increased extracellular potassium concentration, which shifted αON-S RGC excitability to that of αOFF-S cells under baseline potassium conditions. Finally, we investigated the influence of the axon initial segment (AIS) dimensions on RGC spiking. We found that the relationship between AIS length and evoked spike rate varied not only by cell type, but also by the strength of stimulation, suggesting AIS structure alone cannot fully explain the observed differences RGC excitability. Thus, sensitivity to extracellular potassium contributes to differences in intrinsic excitability, a key factor that shapes how RGCs encode visual information.

What is a cell type and how to define it?

Cell types are the basic functional units of an organism. Cell types exhibit diverse phenotypic properties at multiple levels, making them challenging to define, categorize, and understand. This review provides an overview of the basic principles of cell types rooted in evolution and development and discusses approaches to characterize and classify cell types and investigate how they contribute to the organism's function, using the mammalian brain as a primary example. I propose a roadmap toward a conceptual framework and knowledge base of cell types that will enable a deeper understanding of the dynamic changes of cellular function under healthy and diseased conditions.

Diversification of multipotential postmitotic mouse retinal ganglion cell precursors into discrete types

The genesis of broad neuronal classes from multipotential neural progenitor cells has been extensively studied, but less is known about the diversification of a single neuronal class into multiple types. We used single-cell RNA-seq to study how newly born (postmitotic) mouse retinal ganglion cell (RGC) precursors diversify into ~45 discrete types. Computational analysis provides evidence that RGC transcriptomic type identity is not specified at mitotic exit, but acquired by gradual, asynchronous restriction of postmitotic multipotential precursors. Some types are not identifiable until a week after they are generated. Immature RGCs may be specified to project ipsilaterally or contralaterally to the rest of the brain before their type identity emerges. Optimal transport inference identifies groups of RGC precursors with largely nonoverlapping fates, distinguished by selectively expressed transcription factors that could act as fate determinants. Our study provides a framework for investigating the molecular diversification of discrete types within a neuronal class.

The role of microRNAs in glaucoma

In this review, we aim to provide a comprehensive summary of the various microRNAs (miRNAs) shown to be involved in glaucoma and intraocular pressure regulation. miRNAs are short, single-stranded, and noncoding RNAs that regulate gene expression in a number of physiological conditions and human diseases, including glaucoma. Numerous miRNAs display differential expression in glaucoma-affected tissues, such as aqueous humor, tears, trabecular meshwork, and retina analyzed from patients and animal models, suggesting their potential involvement in glaucoma pathogenesis. Several studies summarized here have also investigated the challenge of delivering intact miRNAs to target tissues in order to develop miRNA-based glaucoma therapies. We extend these reports by conducting an additional layer of analysis that integrates the interaction between glaucoma-related miRNAs and glaucoma-associated genes. We conclude with a comprehensive discussion of the therapeutic potential of miRNAs, the cellular pathways that link these miRNAs together, and the most promising miRNAs for future glaucoma research.

DeBruyn and Casagrande manuscripts on tree shrew retinal ganglion cells as a basis for cross-species retina research

The purpose of this brief communication is to make publicly available three unpublished manuscripts on the organization of retinal ganglion cells in the tree shrew. The manuscripts were authored in 1986 by Dr. Edward DeBruyn, a PhD student in the laboratory of the late Dr. Vivien Casagrande at Vanderbilt University. As diurnal animals closely related to primates, tree shrews are ideally suited for comparative analyses of visual structures including the retina. We hope that providing this basic information in a citable form inspires other groups to pursue further characterization of the tree shrew retina using modern techniques.