Mouse Retinal Cell Atlas: Molecular Identification of over Sixty Amacrine Cell Types

Experimental Framework for Assessing Mouse Retinal Regeneration Through Single-Cell RNA-Sequencing

Retinal degenerative diseases including age-related macular degeneration and glaucoma are estimated to currently affect more than 14 million people in the United States, with an increased prevalence of retinal degenerations in aged individuals. An expanding aged population who are living longer forecasts an increased prevalence and economic burden of visual impairments. Improvements to visual health and treatment paradigms for progressive retinal degenerations slow vision loss. However, current treatments fail to remedy the root cause of visual impairments caused by retinal degenerations-loss of retinal neurons. Stimulation of retinal regeneration from endogenous cellular sources presents an exciting treatment avenue for replacement of lost retinal cells. In multiple species including zebrafish and Xenopus, Müller glial cells maintain a highly efficient regenerative ability to reconstitute lost cells throughout the organism's lifespan, highlighting potential therapeutic avenues for stimulation of retinal regeneration in humans. Here, we describe how the application of single-cell RNA-sequencing (scRNA-seq) has enhanced our understanding of Müller glial cell-derived retinal regeneration, including the characterization of gene regulatory networks that facilitate/inhibit regenerative responses. Additionally, we provide a validated experimental framework for cellular preparation of mouse retinal cells as input into scRNA-seq experiments, including insights into experimental design and analyses of resulting data.

Applying Super-Resolution and Tomography Concepts to Identify Receptive Field Subunits in the Retina

Spatially nonlinear stimulus integration by retinal ganglion cells lies at the heart of various computations performed by the retina. It arises from the nonlinear transmission of signals that ganglion cells receive from bipolar cells, which thereby constitute functional subunits within a ganglion cell's receptive field. Inferring these subunits from recorded ganglion cell activity promises a new avenue for studying the functional architecture of the retina. This calls for efficient methods, which leave sufficient experimental time to leverage the acquired knowledge for further investigating identified subunits. Here, we combine concepts from super-resolution microscopy and computed tomography and introduce super-resolved tomographic reconstruction (STR) as a technique to efficiently stimulate and locate receptive field subunits. Simulations demonstrate that this approach can reliably identify subunits across a wide range of model variations, and application in recordings of primate parasol ganglion cells validates the experimental feasibility. STR can potentially reveal comprehensive subunit layouts within only a few tens of minutes of recording time, making it ideal for online analysis and closed-loop investigations of receptive field substructure in retina recordings.

Top-down modulation of the retinal code via histaminergic neurons of the hypothalamus

The mammalian retina is considered an autonomous circuit, yet work dating back to Ramon y Cajal indicates that it receives inputs from the brain. How such inputs affect retinal processing has remained unknown. We confirmed brain-to-retina projections of histaminergic neurons from the mouse hypothalamus. Histamine application ex vivo altered the activity of various retinal ganglion cells (RGCs), including direction-selective RGCs that gained responses to high motion velocities. These results were reproduced in vivo with optic tract recordings where histaminergic retinopetal axons were activated chemogenetically. Such changes could improve vision of fast-moving objects (e.g., while running), which fits with the known increased activity of histaminergic neurons during arousal. An antihistamine drug reduced optomotor responses to high-speed moving stimuli in freely moving mice. In humans, the same antihistamine nonuniformly modulated visual sensitivity across the visual field, indicating an evolutionary conserved function of the histaminergic system. Our findings expose a previously unappreciated role for brain-to-retina projections in modulating retinal function.

Retinal neuroanatomy of two emerging model organisms, the spiny mouse (Acomys dimidiatus) and the Mongolian gerbil (Meriones unguiculatus)

Current research using animal models to investigate retinal cell biology and model retinal degenerative diseases largely utilize small mammals that are nocturnal and lack the ability to restore lost vision. In contrast, the Mongolian gerbil (Meriones) is a diurnal rodent with good photopic vision, and the spiny mouse (Acomys) is a small desert-dwelling rodent with remarkable regenerative capabilities. The goal of this study was to identify antibodies that detect retinal cell classes in Meriones and Acomys, and to describe the retinal anatomy of these two species in comparison to outbred laboratory mice (Mus musculus). Immunohistochemistry was performed on retinal sections with antibodies for various retinal cell types. Sections were imaged by light, fluorescence, and confocal microscopy. Cell density, morphology, and placement were compared between species qualitatively and quantitatively. Our analyses revealed a classic assembly of retinal cells in Meriones and Acomys, with a few deviations compared to Mus. Meriones displayed the highest density of cones and Acomys the lowest. A higher density of bipolar cell bodies in the proximal portion of the inner nuclear layer was observed in both Acomys and Meriones compared to Mus, and both species exhibited an increase in amacrine cell density compared to Mus. Our results provide a foundation for future research into the visual system adaptations of these interesting species.

ON and OFF starburst amacrine cells are controlled by distinct cholinergic pathways

Cholinergic signaling in the retina is mediated by acetylcholine (ACh) released from starburst amacrine cells (SACs), which are key neurons for motion detection. SACs comprise ON and OFF subtypes, which morphologically show mirror symmetry to each other. Although many physiological studies on SACs have targeted ON cells only, the synaptic computation of ON and OFF SACs is assumed to be similar. Recent studies demonstrated that gene expression patterns and receptor types differed between ON and OFF SACs, suggesting differences in their functions. Here, we compared cholinergic signaling pathways between ON and OFF SACs in the mouse retina using the patch clamp technique. The application of ACh increased GABAergic feedback, observed as postsynaptic currents to SACs, in both ON and OFF SACs; however, the mode of GABAergic feedback differed. Nicotinic receptors mediated GABAergic feedback in both ON and OFF SACs, while muscarinic receptors mediated GABAergic feedback in ON SACs only in adults. Neither tetrodotoxin, which blocked action potentials, nor LY354740, which blocked neurotransmitter release from SACs, eliminated ACh-induced GABAergic feedback in SACs. These results suggest that ACh-induced GABAergic feedback in ON and OFF SACs is regulated by different feedback mechanisms in adults and mediated by non-spiking amacrine cells other than SACs.

Advances in the study of the influence of photoreceptors on the development of myopia

This review examines the pivotal role of photoreceptor cells in ocular refraction development, focusing on dopamine (DA) as a key neurotransmitter. Contrary to the earlier view favoring cone cells, recent studies have highlighted the substantial contributions of both rod and cone cells to the visual signaling pathways that influence ocular refractive development. Notably, rod cells appeared to play a central role. Photoreceptor cells interact intricately with circadian rhythms, color vision pathways, and other neurotransmitters, all of which are crucial for the complex mechanisms driving the development of myopia. This review emphasizes that ocular refractive development results from a coordinated interplay between diverse cell types, signaling pathways, and neurotransmitters. This perspective has significant implications for unraveling the complex mechanisms underlying myopia and aiding in the development of more effective prevention and treatment strategies.

Dysregulated energy and protein homeostasis and the loss of GABAergic amacrine cells in aging retina

Aging is a major risk factor for the development or the worsening of retinal degenerative conditions. The intricate network of the neural retina determined that the retinal aging is a complicated process. The aim of this study is to delineate the transcriptomic changes of major retinal neurons during aging in C57BL/6 mice at single-cell level. We analyzed the transcriptional profiles of the photoreceptor, bipolar, amacrine, and Müller glial cells of 1.5-2 and 24-30 months old mice using single-cell RNA sequencing technique. We selectively confirmed the differences in gene expression using immunofluorescence staining and RNA in situ hybridization analysis. We found that each retinal cell type had unique changes upon aging. However, they all showed signs of dysregulated glucose and energy metabolism, and perturbed proteostasis. In particular, old Müller glia exhibited the most profound changes, including the upregulation of cell metabolism, stress-responses, antigen-presentation and immune responses and metal ion homeostasis. The dysregulated gliogenesis and differentiation was confirmed by the presence of Müller glia expressing rod-specific genes in the inner nuclear layer and the outer plexiform layer of the old retina. We further pinpointed the specific loss of GABAergic amacrine cells in old retina. Our study emphasized changes of amacrine and Müller glia during retinal aging, provided resources for further research on the molecular and cellular regulatory mechanisms underlying aging-associated retinal deterioration.

Metabolic Deficits in the Retina of a Familial Dysautonomia Mouse Model

Neurodegenerative retinal diseases such as glaucoma, diabetic retinopathy, Leber's hereditary optic neuropathy (LHON), and dominant optic atrophy (DOA) are marked by progressive death of retinal ganglion cells (RGC). This decline is promoted by structural and functional mitochondrial deficits, including electron transport chain (ETC) impairments, increased oxidative stress, and reduced energy (ATP) production. These cellular mechanisms associated with progressive optic nerve atrophy have been similarly observed in familial dysautonomia (FD) patients, who experience gradual loss of visual acuity due to the degeneration of RGCs, which is thought to be caused by a breakdown of mitochondrial structures, and a disruption in ETC function. Retinal metabolism plays a crucial role in meeting the elevated energetic demands of this tissue, and recent characterizations of FD patients' serum and stool metabolomes have indicated alterations in central metabolic processes and potential systemic deficits of taurine, a small molecule essential for retina and overall eye health. The present study sought to elucidate metabolic alterations that contribute to the progressive degeneration of RGCs observed in FD. Additionally, a critical subpopulation of retinal interneurons, the dopaminergic amacrine cells, mediate the integration and modulation of visual information in a time-dependent manner to RGCs. As these cells have been associated with RGC loss in the neurodegenerative disease Parkinson's, which shares hallmarks with FD, a targeted analysis of the dopaminergic amacrine cells and their product, dopamine, was also undertaken. One dimensional (1D) proton (1H) nuclear magnetic resonance (NMR) spectroscopy, mass spectrometry, and retinal histology methods were employed to characterize retinae from the retina-specific Elp1 conditional knockout (CKO) FD mouse model (Pax6-Cre; Elp1LoxP/LoxP). Metabolite alterations correlated temporally with progressive RGC degeneration and were associated with reduced mitochondrial function, alterations in ATP production through the Cahill and mini-Krebs cycles, and phospholipid metabolism. Dopaminergic amacrine cell populations were reduced at timepoints P30-P90, and dopamine levels were 25-35% lower in CKO retinae compared to control retinae at P60. Overall, this study has expanded upon our current understanding of retina pathology in FD. This knowledge may apply to other retinal diseases that share hallmark features with FD and may help guide new avenues for novel non-invasive therapeutics to mitigate the progressive optic neuropathy in FD.

Reconfiguration of the visual code and retinal cell type complement in closely related diurnal and nocturnal mice

How does evolution act on neuronal populations to match computational characteristics to functional demands? We address this problem by comparing visual code and retinal cell composition in closely related murid species with different behaviours. Rhabdomys pumilio are diurnal and have substantially thicker inner retina and larger visual thalamus than nocturnal Mus musculus. High-density electrophysiological recordings of visual response features in the dorsal lateral geniculate nucleus (dLGN) reveals that Rhabdomys attains higher spatiotemporal acuity both by denser coverage of the visual scene and a selective expansion of elements of the code characterised by non-linear spatiotemporal summation. Comparative analysis of single cell transcriptomic cell atlases reveals that realignment of the visual code is associated with increased relative abundance of bipolar and ganglion cell types supporting OFF and ON-OFF responses. These findings demonstrate how changes in retinal cell complement can reconfigure the coding of visual information to match changes in visual needs.

Morphological and electrophysiological characterization of a novel displaced astrocyte in the mouse retina

Astrocytes throughout the central nervous system are heterogeneous in both structure and function. This diversity leads to tissue-specific specialization where morphology is adapted to the surrounding neuronal circuitry, as seen in Bergman glia of the cerebellum and Müller glia of the retina. Because morphology can be a differentiating factor for cellular classification, we recently developed a mouse where glial-fibrillary acidic protein (GFAP)-expressing cells stochastically label for full membranous morphology. Here we utilize this tool to investigate whether morphological and electrophysiological features separate types of mouse retinal astrocytes. In this work, we report on a novel glial population found in the inner plexiform layer and ganglion cell layer which expresses the canonical astrocyte markers GFAP, S100β, connexin-43, Sox2 and Sox9. Apart from their retinal layer localization, these cells are unique in their radial distribution. They are notably absent from the mid-retina but are heavily concentrated near the optic nerve head, and to a lesser degree the peripheral retina. Additionally, their morphology is distinct from both nerve fiber layer astrocytes and Müller glia, appearing more similar to amacrine cells. Despite this structural similarity, these cells lack protein expression of common neuronal markers. Additionally, they do not exhibit action potentials, but rather resemble astrocytes and Müller glia in their small amplitude, graded depolarization to both light onset and offset. Their structure, protein expression, physiology, and intercellular connections suggest that these cells are astrocytic, displaced from their counterparts in the nerve fiber layer. As such, we refer to these cells as displaced retinal astrocytes.

Cholecystokinin receptor type A are involved in the circadian rhythm of the mouse retina

The retina is the only organ projecting external light to the suprachiasmatic nucleus. Cholecystokinin receptor type A (Cckar/Cckar) is one of the essential factors for light reception in retinal cells. As there was a lack of literature on the matter, we aimed to elucidate the cause of the time-dependent phase change in clock gene expression. We found that Cckar mRNA expression in retinal cells exhibited diurnal variations. The rhythm of expression of the clock gene Per1/Per2 in retinal cells was altered in Cckar -/- mice. The light sensitivity of retinal cells was evaluated in wild-type mice, which showed c-Fos was activated in the ganglion cell layer more than in the inner granular layer. This increase in the number of c-Fos-positive cells was suppressed by lorglumide, a Cckar antagonist. Treatment of rat retina primary cells with lorglumide suppressed Per2 transcription, which was altered in a time-dependent manner relative to the Per2 expression. Light irradiation studies in Cckar -/- mice did not exhibit an increase in Period expression in the suprachiasmatic nucleus. These results indicate that Cckar is among the factors that regulate the cycle of clock genes on the retina. Cckar knockout attenuates the light responsiveness of suprachiasmatic nucleus and reduces the expression amplitude of Period genes in the retina. Thus, Cckar may contribute to entrainment of the light environment and maintenance of the expression cycle of Period gene, which is one of the core clock genes.

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's functions and reveal its importance in bipolar cells for high frequency light response.

Fluorescent identification of axons, dendrites and soma of neuronal retinal ganglion cells with a genetic marker as a tool for facilitating the study of neurodegeneration

This study characterizes a fluorescent Slc17a6 -tdTomato neuronal reporter mouse line offering strong labeling in axons throughout the optic nerve, dendrites and soma in 99% of retinal ganglion cells (RGCs). The model facilitates neuronal assessment ex vivo with wholemounts quantified to show neurodegeneration following optic nerve crush or elevated IOP as related to glaucoma, in vitro with robust Ca 2+ responses to P2X7 receptor stimulation in neuronal cultures, and in vivo using a confocal scanning laser ophthalmoscope (cSLO). While the tdTomato signal showed strong overlap with RGC markers, BRN3A and RBPMS, there was no cross-labeling of displaced amacrine cells in the ganglion cell layer. Controls indicated no impact of Slc17a6 -tdTomato expression on light-dependent neuronal function, as determined with a microelectrode array (MEA), or on structure, as measured with optical coherence tomography (OCT). In summary, this novel neuronal reporter mouse model offers an effective means to increase the efficiency for real-time, specific visualization of retinal ganglion cells. It holds substantial promise for enhancing our understanding of RGC pathology in glaucoma and other diseases of the optic nerve, and could facilitate the screening of targeted therapeutic interventions for neurodegeneration.

Abstract Figure

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.

A specific olfactory bulb interneuron subtype Tpbg/5T4 generated at embryonic and neonatal stages

Various mammals have shown that sensory stimulation plays a crucial role in regulating the development of diverse structures, such as the olfactory bulb (OB), cerebral cortex, hippocampus, and retina. In the OB, the dendritic development of excitatory projection neurons like mitral/tufted cells is influenced by olfactory experiences. Odor stimulation is also essential for the dendritic development of inhibitory OB interneurons, such as granule and periglomerular cells, which are continuously produced in the ventricular-subventricular zone throughout life. Based on the morphological and molecular features, OB interneurons are classified into several subtypes. The role for each interneuron subtype in the control of olfactory behavior remains poorly understood due to lack of each specific marker. Among the several OB interneuron subtypes, a specific granule cell subtype, which expresses the oncofetal trophoblast glycoprotein (Tpbg or 5T4) gene, has been reported to be required for odor detection and discrimination behavior. This review will primarily focus on elucidating the contribution of different granule cell subtypes, including the Tpbg/5T4 subtype, to olfactory processing and behavior during the embryonic and adult stages.

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.

The vertebrate retina: a window into the evolution of computation in the brain

Animal brains are probably the most complex computational machines on our planet, and like everything in biology, they are the product of evolution. Advances in developmental and palaeobiology have been expanding our general understanding of how nervous systems can change at a molecular and structural level. However, how these changes translate into altered function - that is, into 'computation' - remains comparatively sparsely explored. What, concretely, does it mean for neuronal computation when neurons change their morphology and connectivity, when new neurons appear or old ones disappear, or when transmitter systems are slowly modified over many generations? And how does evolution use these many possible knobs and dials to constantly tune computation to give rise to the amazing diversity in animal behaviours we see today? Addressing these major gaps of understanding benefits from choosing a suitable model system. Here, I present the vertebrate retina as one perhaps unusually promising candidate. The retina is ancient and displays highly conserved core organisational principles across the entire vertebrate lineage, alongside a myriad of adjustments across extant species that were shaped by the history of their visual ecology. Moreover, the computational logic of the retina is readily interrogated experimentally, and our existing understanding of retinal circuits in a handful of species can serve as an anchor when exploring the visual circuit adaptations across the entire vertebrate tree of life, from fish deep in the aphotic zone of the oceans to eagles soaring high up in the sky.

Retinal ganglion cell type-specific expression of synuclein family members revealed by scRNA-sequencing

Synuclein family members (Snca, Sncb, and Scng) are expressed in the retina, but their precise locations and roles are poorly understood. We performed an extensive analysis of the single-cell transcriptome in healthy and injured retinas to investigate their expression patterns and roles. We observed the expression of all synuclein family members in retinal ganglion cells (RGCs), which remained consistent across species (human, mouse, and chicken). We unveiled differential expression of Snca across distinct clusters (highly expressed in most), while Sncb and Sncg displayed uniform expression across all clusters. Further, we observed a decreased expression in RGCs following traumatic axonal injury. However, the proportion of α-Syn-positive RGCs in all RGCs and α-Syn-positive intrinsically photosensitive retinal ganglion cells (ipRGCs) in all ipRGCs remained unaltered. Lastly, we identified changes in communication patterns preceding cell death, with particular significance in the pleiotrophin-nucleolin (Ptn-Ncl) and neural cell adhesion molecule signaling pathways, where communication differences were pronounced between cells with varying expression levels of Snca. Our study employs an innovative approach using scRNA-seq to characterize synuclein expression in health retinal cells, specifically focusing on RGC subtypes, advances our knowledge of retinal physiology and pathology.

Classification and analysis of retinal interneurons by computational structure under natural scenes

Inhibitory neurons are diverse across the brain, but for the visual system we lack the ability to functionally classify these neurons under complex natural stimuli. Here we take the approach of classifying retinal amacrine cell responses to natural scenes using optical recording and an interpretable neural network model. We fit mouse amacrine cell responses to a two-layer convolutional neural network model of a class shown previously to accurately capture salamander ganglion cell responses to natural scenes. Using an approach from interpretable machine learning, we determined for each stimulus the model interneurons that generated each amacrine response, analogous to the set of bipolar cells that target the amacrine population. From this analysis we clustered amacrine cells not by their natural scene responses, but by the model presynaptic neurons that constructed those responses, conservatively finding approximately seven groups by this approach. By analyzing the set of model presynaptic input neurons for each amacrine cluster, we find that distributed rather than dedicated inputs generate natural scene responses for different amacrine cell types. Additional analyses revealed distinct transient and sustained modes exhibited by the network during the response to simple flashes. These results give insight into the computational structure of how the diverse amacrine cell population responds to natural scenes, and generate multiple quantitative hypotheses for how synaptic inputs generate those responses.

FOXC1 regulates endothelial CD98 (LAT1/4F2hc) expression in retinal angiogenesis and blood-retina barrier formation

Angiogenesis, the growth of new blood vessels from pre-existing vasculature, is essential for the development of new organ systems, but transcriptional control of angiogenesis remains incompletely understood. Here we show that FOXC1 is essential for retinal angiogenesis. Endothelial cell (EC)-specific loss of Foxc1 impairs retinal vascular growth and expression of Slc3a2 and Slc7a5, which encode the heterodimeric CD98 (LAT1/4F2hc) amino acid transporter and regulate the intracellular transport of essential amino acids and activation of the mammalian target of rapamycin (mTOR). EC-Foxc1 deficiency diminishes mTOR activity, while administration of the mTOR agonist MHY-1485 rescues perturbed retinal angiogenesis. EC-Foxc1 expression is required for retinal revascularization and resolution of neovascular tufts in a model of oxygen-induced retinopathy. Foxc1 is also indispensable for pericytes, a critical component of the blood-retina barrier during retinal angiogenesis. Our findings establish FOXC1 as a crucial regulator of retinal vessels and identify therapeutic targets for treating retinal vascular disease.

Spatial patterning controls neuron numbers in the Drosophila visual system

Neurons must be made in the correct proportions to communicate with the appropriate synaptic partners and form functional circuits. In the Drosophila visual system, multiple subtypes of distal medulla (Dm) inhibitory interneurons are made in distinct, reproducible numbers-from 5 to 800 per optic lobe. These neurons are born from a crescent-shaped neuroepithelium called the outer proliferation center (OPC), which can be subdivided into specific domains based on transcription factor and growth factor expression. We fate mapped Dm neurons and found that more abundant neural types are born from larger neuroepithelial subdomains, while less abundant subtypes are born from smaller ones. Additionally, morphogenetic Dpp/BMP signaling provides a second layer of patterning that subdivides the neuroepithelium into smaller domains to provide more granular control of cell proportions. Apoptosis appears to play a minor role in regulating Dm neuron abundance. This work describes an underappreciated mechanism for the regulation of neuronal stoichiometry.

Differential Expression Analysis Identifies Candidate Synaptogenic Molecules for Wiring Direction-Selective Circuits in the Retina

An organizational feature of neural circuits is the specificity of synaptic connections. A striking example is the direction-selective (DS) circuit of the retina. There are multiple subtypes of DS retinal ganglion cells (DSGCs) that prefer motion along one of four preferred directions. This computation is mediated by selective wiring of a single inhibitory interneuron, the starburst amacrine cell (SAC), with each DSGC subtype preferentially receiving input from a subset of SAC processes. We hypothesize that the molecular basis of this wiring is mediated in part by unique expression profiles of DSGC subtypes. To test this, we first performed paired recordings from isolated mouse retinas of both sexes to determine that postnatal day 10 (P10) represents the age at which asymmetric synapses form. Second, we performed RNA sequencing and differential expression analysis on isolated P10 ON-OFF DSGCs tuned for either nasal or ventral motion and identified candidates which may promote direction-specific wiring. We then used a conditional knock-out strategy to test the role of one candidate, the secreted synaptic organizer cerebellin-4 (Cbln4), in the development of DS tuning. Using two-photon calcium imaging, we observed a small deficit in directional tuning among ventral-preferring DSGCs lacking Cbln4, though whole-cell voltage-clamp recordings did not identify a significant change in inhibitory inputs. This suggests that Cbln4 does not function primarily via a cell-autonomous mechanism to instruct wiring of DS circuits. Nevertheless, our transcriptomic analysis identified unique candidate factors for gaining insights into the molecular mechanisms that instruct wiring specificity in the DS circuit.

Asymmetric Activation of ON and OFF Pathways in the Degenerated Retina

Retinal prosthetics are one of the leading therapeutic strategies to restore lost vision in patients with retinitis pigmentosa and age-related macular degeneration. Much work has described patterns of spiking in retinal ganglion cells (RGCs) in response to electrical stimulation, but less work has examined the underlying retinal circuitry that is activated by electrical stimulation to drive these responses. Surprisingly, little is known about the role of inhibition in generating electrical responses or how inhibition might be altered during degeneration. Using whole-cell voltage-clamp recordings during subretinal electrical stimulation in the rd10 and wild-type (wt) retina, we found electrically evoked synaptic inputs differed between ON and OFF RGC populations, with ON cells receiving mostly excitation and OFF cells receiving mostly inhibition and very little excitation. We found that the inhibition of OFF bipolar cells limits excitation in OFF RGCs, and a majority of both pre- and postsynaptic inhibition in the OFF pathway arises from glycinergic amacrine cells, and the stimulation of the ON pathway contributes to inhibitory inputs to the RGC. We also show that this presynaptic inhibition in the OFF pathway is greater in the rd10 retina, compared with that in the wt retina.

Abnormal α-Synuclein Aggregates Cause Synaptic- and Microcircuit-Specific Deficits in the Retinal Rod Pathway

α-Synuclein (α-Syn) is a key determinator of Parkinson disease (PD) pathology, but synapse and microcircuit pathologies in the retina underlying visual dysfunction are poorly understood. Herein, histochemical and ultrastructural analyses and ophthalmologic measurements in old transgenic M83 PD model (mice aged 16 to 18 months) indicated that abnormal α-Syn aggregation in the outer plexiform layer (OPL) was associated with degeneration in the C-terminal binding protein 2 (CtBP2)+ ribbon synapses of photoreceptor terminals and protein kinase C alpha (PKCα)+ rod bipolar cell terminals, whereas α-Syn aggregates in the inner retina correlated with the reduction and degeneration of tyrosine hydroxylase- and parvalbumin-positive amacrine cells. Phosphorylated Ser129 α-synuclein expression was strikingly restricted in the OPL, with the most severe degenerations in the entire retina, including mitochondrial degeneration and loss of ribbon synapses in 16- to 18-month-old mice. These synapse- and microcircuit-specific deficits of the rod pathway at the CtBP2+ rod terminals and PKCα+ rod bipolar and amacrine cells were associated with attenuated a- and b-wave amplitudes and oscillatory potentials on the electroretinogram. They were also associated with the impairment of visual functions, including reduced contrast sensitivity and impairment of the middle range of spatial frequencies. Collectively, these findings demonstrate that α-Syn aggregates cause the synapse- and microcircuit-specific deficits of the rod pathway and the most severe damage to the OPL, providing the retinal synaptic and microcircuit basis for visual dysfunctions in PD.

Pten regulates endocytic trafficking of cell adhesion and Wnt signaling molecules to pattern the retina

The retina is exquisitely patterned, with neuronal somata positioned at regular intervals to completely sample the visual field. Here, we show that phosphatase and tensin homolog (Pten) controls starburst amacrine cell spacing by modulating vesicular trafficking of cell adhesion molecules and Wnt proteins. Single-cell transcriptomics and double-mutant analyses revealed that Pten and Down syndrome cell adhesion molecule Dscam) are co-expressed and function additively to pattern starburst amacrine cell mosaics. Mechanistically, Pten loss accelerates the endocytic trafficking of DSCAM, FAT3, and MEGF10 off the cell membrane and into endocytic vesicles in amacrine cells. Accordingly, the vesicular proteome, a molecular signature of the cell of origin, is enriched in exocytosis, vesicle-mediated transport, and receptor internalization proteins in Pten conditional knockout (PtencKO) retinas. Wnt signaling molecules are also enriched in PtencKO retinal vesicles, and the genetic or pharmacological disruption of Wnt signaling phenocopies amacrine cell patterning defects. Pten thus controls vesicular trafficking of cell adhesion and signaling molecules to establish retinal amacrine cell mosaics.

Evolutionary and developmental specialization of foveal cell types in the marmoset

In primates, high-acuity vision is mediated by the fovea, a small specialized central region of the retina. The fovea, unique to the anthropoid lineage among mammals, undergoes notable neuronal morphological changes during postnatal maturation. However, the extent of cellular similarity across anthropoid foveas and the molecular underpinnings of foveal maturation remain unclear. Here, we used high-throughput single-cell RNA sequencing to profile retinal cells of the common marmoset (Callithrix jacchus), an early divergent in anthropoid evolution from humans, apes, and macaques. We generated atlases of the marmoset fovea and peripheral retina for both neonates and adults. Our comparative analysis revealed that marmosets share almost all their foveal types with both humans and macaques, highlighting a conserved cellular structure among primate foveas. Furthermore, by tracing the developmental trajectory of cell types in the foveal and peripheral retina, we found distinct maturation paths for each. In-depth analysis of gene expression differences demonstrated that cone photoreceptors and Müller glia (MG), among others, show the greatest molecular divergence between these two regions. Utilizing single-cell ATAC-seq and gene-regulatory network inference, we uncovered distinct transcriptional regulations differentiating foveal cones from their peripheral counterparts. Further analysis of predicted ligand-receptor interactions suggested a potential role for MG in supporting the maturation of foveal cones. Together, these results provide valuable insights into foveal development, structure, and evolution.

Retinal VIP-amacrine cells: their development, structure, and function

Amacrine cells (ACs) are the most structurally and functionally diverse neuron type in the retina. Different ACs have distinct functions, such as neuropeptide secretion and inhibitory connection. Vasoactive intestinal peptide (VIP) -ergic -ACs are retina gamma-aminobutyric acid (GABA) -ergic -ACs that were discovered long ago. They secrete VIP and form connections with bipolar cells (BCs), other ACs, and retinal ganglion cells (RGCs). They have a specific structure, density, distribution, and function. They play an important role in myopia, light stimulated responses, retinal vascular disease and other ocular diseases. Their significance in the study of refractive development and disease is increasing daily. However, a systematic review of the structure and function of retinal VIP-ACs is lacking. We discussed the detailed characteristics of VIP-ACs from every aspect across species and providing systematic knowledge base for future studies. Our review led to the main conclusion that retinal VIP-ACs develop early, and although their morphology and distribution across species are not the same, they have similar functions in a wide range of ocular diseases based on their function of secreting neuropeptides and forming inhibitory connections with other cells.

Regulation of Diseases-Associated Microglia in the Optic Nerve by Lipoxin B4 and Ocular Hypertension

Background

The resident astrocyte-retinal ganglion cell (RGC) lipoxin circuit is impaired during retinal stress, which includes ocular hypertension-induced neuropathy. Lipoxin B4 produced by homeostatic astrocytes directly acts on RGCs to increase survival and function in ocular hypertension-induced neuropathy. RGC death in the retina and axonal degeneration in the optic nerve are driven by the complex interactions between microglia and macroglia. Whether LXB4 neuroprotective actions include regulation of other cell types in the retina and/or optic nerve is an important knowledge gap.

Methods

Cellular targets and signaling of LXB4 in the retina were defined by single-cell RNA sequencing. Retinal neurodegeneration was induced by injecting silicone oil into the anterior chamber of the mouse eyes, which induced sustained and stable ocular hypertension. Morphological characterization of microglia populations in the retina and optic nerve was established by MorphOMICs and pseudotime trajectory analyses. The pathways and mechanisms of action of LXB4 in the optic nerve were investigated using bulk RNA sequencing. Transcriptomics data was validated by qPCR and immunohistochemistry. Differences between experimental groups was assessed by Student's t-test and one-way ANOVA.

Results

Single-cell transcriptomics identified microglia as a primary target for LXB4 in the healthy retina. LXB4 downregulated genes that drive microglia environmental sensing and reactivity responses. Analysis of microglial function revealed that ocular hypertension induced distinct, temporally defined, and dynamic phenotypes in the retina and, unexpectedly, in the distal myelinated optic nerve. Microglial expression of CD74, a marker of disease-associated microglia in the brain, was only induced in a unique population of optic nerve microglia, but not in the retina. Genetic deletion of lipoxin formation correlated with the presence of a CD74 optic nerve microglia population in normotensive eyes, while LXB4 treatment during ocular hypertension shifted optic nerve microglia toward a homeostatic morphology and non-reactive state and downregulated the expression of CD74. Furthermore, we identified a correlation between CD74 and phospho-phosphoinositide 3-kinases (p-PI3K) expression levels in the optic nerve, which was reduced by LXB4 treatment.

Conclusion

We identified early and dynamic changes in the microglia functional phenotype, reactivity, and induction of a unique CD74 microglia population in the distal optic nerve as key features of ocular hypertension-induced neurodegeneration. Our findings establish microglia regulation as a novel LXB4 target in the retina and optic nerve. LXB4 maintenance of a homeostatic optic nerve microglia phenotype and inhibition of a disease-associated phenotype are potential neuroprotective mechanisms for the resident LXB4 pathway.

Atypical retinal function in a mouse model of Fragile X syndrome

Altered function of peripheral sensory neurons is an emerging mechanism for symptoms of autism spectrum disorders. Visual sensitivities are common in autism, but whether differences in the retina might underlie these sensitivities is not well-understood. We explored retinal function in the Fmr1 knockout model of Fragile X syndrome, focusing on a specific type of retinal neuron, the "sustained On alpha" retinal ganglion cell. We found that these cells exhibit changes in dendritic structure and dampened responses to light in the Fmr1 knockout. We show that decreased light sensitivity is due to increased inhibitory input and reduced E-I balance. The change in E-I balance supports maintenance of circuit excitability similar to what has been observed in cortex. These results show that loss of Fmr1 in the mouse retina affects sensory function of one retinal neuron type. Our findings suggest that the retina may be relevant for understanding visual function in Fragile X syndrome.

A Non-canonical Excitatory PV RGC-PV SC Visual Pathway for Mediating the Looming-evoked Innate Defensive Response

Parvalbumin-positive retinal ganglion cells (PV+ RGCs) are an essential subset of RGCs found in various species. However, their role in transmitting visual information remains unclear. Here, we characterized PV+ RGCs in the retina and explored the functions of the PV+ RGC-mediated visual pathway. By applying multiple viral tracing strategies, we investigated the downstream of PV+ RGCs across the whole brain. Interestingly, we found that the PV+ RGCs provided direct monosynaptic input to PV+ excitatory neurons in the superficial layers of the superior colliculus (SC). Ablation or suppression of SC-projecting PV+ RGCs abolished or severely impaired the flight response to looming visual stimuli in mice without affecting visual acuity. Furthermore, using transcriptome expression profiling of individual cells and immunofluorescence colocalization for RGCs, we found that PV+ RGCs are predominant glutamatergic neurons. Thus, our findings indicate the critical role of PV+ RGCs in an innate defensive response and suggest a non-canonical subcortical visual pathway from excitatory PV+ RGCs to PV+ SC neurons that regulates looming visual stimuli. These results provide a potential target for intervening and treating diseases related to this circuit, such as schizophrenia and autism.

Role of PKN1 in Retinal Cell Type Formation

We recently identified PKN1 as a developmentally active gatekeeper of the transcription factor neuronal differentiation-2 (NeuroD2) in several brain areas. Since NeuroD2 plays an important role in amacrine cell (AC) and retinal ganglion cell (RGC) type formation, we aimed to study the expression of NeuroD2 in the postnatal retina of WT and Pkn1-/- animals, with a particular focus on these two cell types. We show that PKN1 is broadly expressed in the retina and that the gross retinal structure is not different between both genotypes. Postnatal retinal NeuroD2 levels were elevated upon Pkn1 knockout, with Pkn1-/- retinae showing more NeuroD2+ cells in the lower portion of the inner nuclear layer. Accordingly, immunohistochemical analysis revealed an increased amount of AC in postnatal and adult Pkn1-/- retinae. There were no differences in horizontal cell, bipolar cell, glial cell and RGC numbers, nor defective axon guidance to the optic chiasm or tract upon Pkn1 knockout. Interestingly, we did, however, see a specific reduction in SMI-32+ α-RGC in Pkn1-/- retinae. These results suggest that PKN1 is important for retinal cell type formation and validate PKN1 for future studies focusing on AC and α-RGC specification and development.

Cone Synaptic function is modulated by the leucine rich repeat (LRR) adhesion molecule LRFN2

Daylight vision is mediated by cone photoreceptors in vertebrates, which synapse with bipolar cells (BCs) and horizontal (HCs) cells. This cone synapse is functionally and anatomically complex, connecting to 8 types of depolarizing BCs (DBCs) and 5 types of hyperpolarizing BCs (HBCs) in mice. The dendrites of DBCs and HCs cells make invaginating ribbon synapses with the cone axon terminal, while HBCs form flat synapses with the cone pedicles. The molecular architecture that underpins this organization is relatively poorly understood. To identify new proteins involved in synapse formation and function we used an unbiased proteomic approach and identified LRFN2 (leucine-rich repeat and fibronectin III domain-containing 2) as a component of the DBC signaling complex. LRFN2 is selectively expressed at cone terminals and co-localizes with PNA, and other DBC signalplex members. In LRFN2 deficient mice, the synaptic markers: LRIT3, ELFN2, mGluR6, TRPM1 and GPR179 are properly localized. Similarly, LRFN2 expression and localization is not dependent on these synaptic proteins. In the absence of LRFN2 the cone-mediated photopic electroretinogram b-wave amplitude is reduced at the brightest flash intensities. These data demonstrate that LRFN2 absence compromises normal synaptic transmission between cones and cone DBCs.Significance Statement Signaling between cone photoreceptors and the downstream bipolar cells is critical to normal vision. Cones synapse with 13 different types of bipolar cells forming an invaginating ribbon synapses with 8 types, and flat synapses with 5 types, to form one of the most complex synapses in the brain. In this report a new protein, LRFN2 (leucine-rich repeat and fibronectin III domain-containing 2), was identified that is expressed at the cone synapse. Using Lrfn2 knockout mice we show LRFN2 is required for the normal cone signaling.

Single-cell RNA sequencing reveals roles of unique retinal microglia types in early diabetic retinopathy

BACKGROUND

The pathophysiological mechanisms of diabetic retinopathy (DR), a blinding disease, are intricate. DR was thought to be a microvascular disease previously. However, growing studies have indicated that the retinal microglia-induced inflammation precedes microangiopathy. The binary concept of microglial M1/M2 polarization paradigms during inflammatory activation has been debated. In this study, we confirmed microglia had the most significant changes in early DR using single-cell RNA sequencing.

METHODS

A total of five retinal specimens were collected from donor SD rats. Changes in various cells of the retina at the early stage of DR were analyzed using single-cell sequencing technology.

RESULTS

We defined three new microglial subtypes at cellular level, including two M1 types (Egr2+ M1 and Egr2- M1) and one M2 type. We also revealed the anatomical location between these subtypes, the dynamic changes of polarization phenotypes, and the possible activation sequence and mutual activation regulatory mechanism of different cells. Furthermore, we constructed an inflammatory network involving microglia, blood-derived macrophages and other retinal nonneuronal cells. The targeted study of new disease-specific microglial subtypes can shorten the time for drug screening and clinical application, which provided insight for the early control and reversal of DR.

CONCLUSIONS

We found that microglia show the most obvious differential expression changes in early DR and reveal the changes in microglia in a high-glucose microenvironment at the single-cell level. Our comprehensive analysis will help achieve early reversal and control the occurrence and progression of DR.

Touch receptor end-organ innervation and function require sensory neuron expression of the transcription factor Meis2

Touch sensation is primarily encoded by mechanoreceptors, called low-threshold mechanoreceptors (LTMRs), with their cell bodies in the dorsal root ganglia. Because of their great diversity in terms of molecular signature, terminal endings morphology, and electrophysiological properties, mirroring the complexity of tactile experience, LTMRs are a model of choice to study the molecular cues differentially controlling neuronal diversification. While the transcriptional codes that define different LTMR subtypes have been extensively studied, the molecular players that participate in their late maturation and in particular in the striking diversity of their end-organ morphological specialization are largely unknown. Here we identified the TALE homeodomain transcription factor Meis2 as a key regulator of LTMRs target-field innervation in mice. Meis2 is specifically expressed in cutaneous LTMRs, and its expression depends on target-derived signals. While LTMRs lacking Meis2 survived and are normally specified, their end-organ innervations, electrophysiological properties, and transcriptome are differentially and markedly affected, resulting in impaired sensory-evoked behavioral responses. These data establish Meis2 as a major transcriptional regulator controlling the orderly formation of sensory neurons innervating peripheral end organs required for light touch.

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.

Sustained Extracellular Electrical Stimulation Modulates the Permeability of Gap Junctions in rd1 Mouse Retina with Photoreceptor Degeneration

Neurons build vast gap junction-coupled networks (GJ-nets) that are permeable to ions or small molecules, enabling lateral signaling. Herein, we investigate (1) the effect of blinding diseases on GJ-nets in mouse retinas and (2) the impact of electrical stimulation on GJ permeability. GJ permeability was traced in the acute retinal explants of blind retinal degeneration 1 (rd1) mice using the GJ tracer neurobiotin. The tracer was introduced via the edge cut method into the GJ-net, and its spread was visualized in histological preparations (fluorescent tagged) using microscopy. Sustained stimulation was applied to modulate GJ permeability using a single large electrode. Our findings are: (1) The blind rd1 retinas displayed extensive intercellular coupling via open GJs. Three GJ-nets were identified: horizontal, amacrine, and ganglion cell networks. (2) Sustained stimulation significantly diminished the tracer spread through the GJs in all the cell layers, as occurs with pharmaceutical inhibition with carbenoxolone. We concluded that the GJ-nets of rd1 retinas remain coupled and functional after blinding disease and that their permeability is regulatable by sustained stimulation. These findings are essential for understanding molecular signaling in diseases over coupled networks and therapeutic approaches using electrical implants, such as eliciting visual sensations or suppressing cortical seizures.

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.

A presynaptic source drives differing levels of surround suppression in two mouse retinal ganglion cell types

In early sensory systems, cell-type diversity generally increases from the periphery into the brain, resulting in a greater heterogeneity of responses to the same stimuli. Surround suppression is a canonical visual computation that begins within the retina and is found at varying levels across retinal ganglion cell types. Our results show that heterogeneity in the level of surround suppression occurs subcellularly at bipolar cell synapses. Using single-cell electrophysiology and serial block-face scanning electron microscopy, we show that two retinal ganglion cell types exhibit very different levels of surround suppression even though they receive input from the same bipolar cell types. This divergence of the bipolar cell signal occurs through synapse-specific regulation by amacrine cells at the scale of tens of microns. These findings indicate that each synapse of a single bipolar cell can carry a unique visual signal, expanding the number of possible functional channels at the earliest stages of visual processing.

Biotin-cGMP and -cAMP are able to permeate through the gap junctions of some amacrine cells in the mouse retina despite their large size

Gap junctions transmit electrical signals in neurons and serve metabolic coupling and chemical communication. Gap junctions are made of intercellular channels with large pores, allowing ions and small molecules to permeate. In the mammalian retina, intercellular coupling fulfills many vital functions in visual signal processing but is also implicated in promoting cell death after insults, such as excitotoxicity or hypoxia. Conversely, some studies also suggested a role for retinal gap junctions in neuroprotection. Recently, gap junctions were also advocated as conduits for therapeutic drug delivery in neurodegenerative disorders. This requires the permeation of rather large molecules through retinal gap junctions. However, the permeability of retinal networks for molecules >0.6 kDa has not been tested systematically. Here, we used the cut-loading method and probed gap junctional networks in the mouse retina for their permeability to cGMP and cAMP coupled to Biotin, using the well-characterized tracer Neurobiotin as control. Biotin-cGMP and -cAMP have a molecular weight of >0.8 kDa. We show that they cannot pass the gap junctions of horizontal cells but can permeate through the gap junctions of specific amacrine cells in the inner retina. These amacrine cells do not comprise AII amacrine cells and nitric oxide-releasing amacrine cells but some unknown type. In summary, we show that some retinal gap junctions are large enough to let molecules >0.8 kDa pass, making the intercellular delivery of therapeutic agents - already successfully exploited, for example, in cancer - also feasible in neurodegenerative diseases.

MarsGT: Multi-omics analysis for rare population inference using single-cell graph transformer

Rare cell populations are key in neoplastic progression and therapeutic response, offering potential intervention targets. However, their computational identification and analysis often lag behind major cell types. To fill this gap, we introduce MarsGT: Multi-omics Analysis for Rare population inference using a Single-cell Graph Transformer. It identifies rare cell populations using a probability-based heterogeneous graph transformer on single-cell multi-omics data. MarsGT outperforms existing tools in identifying rare cells across 550 simulated and four real human datasets. In mouse retina data, it reveals unique subpopulations of rare bipolar cells and a Müller glia cell subpopulation. In human lymph node data, MarsGT detects an intermediate B cell population potentially acting as lymphoma precursors. In human melanoma data, it identifies a rare MAIT-like population impacted by a high IFN-I response and reveals the mechanism of immunotherapy. Hence, MarsGT offers biological insights and suggests potential strategies for early detection and therapeutic intervention of disease.

Synaptic Origins of the Complex Receptive Field Structure in Primate Smooth Monostratified Retinal Ganglion Cells

Considerable progress has been made in studying the receptive fields of the most common primate retinal ganglion cell (RGC) types, such as parasol RGCs. Much less is known about the rarer primate RGC types and the circuitry that gives rise to noncanonical receptive field structures. The goal of this study was to analyze synaptic inputs to smooth monostratified RGCs to determine the origins of their complex spatial receptive fields, which contain isolated regions of high sensitivity called "hotspots." Interestingly, smooth monostratified RGCs co-stratify with the well-studied parasol RGCs and are thus constrained to receiving input from bipolar and amacrine cells with processes sharing the same layer, raising the question of how their functional differences originate. Through 3D reconstructions of circuitry and synapses onto ON smooth monostratified and ON parasol RGCs from central macaque retina, we identified four distinct sampling strategies employed by smooth and parasol RGCs to extract diverse response properties from co-stratifying bipolar and amacrine cells. The two RGC types differed in the proportion of amacrine cell input, relative contributions of co-stratifying bipolar cell types, amount of synaptic input per bipolar cell, and spatial distribution of bipolar cell synapses. Our results indicate that the smooth RGC's complex receptive field structure arises through spatial asymmetries in excitatory bipolar cell input which formed several discrete clusters comparable with physiologically measured hotspots. Taken together, our results demonstrate how the striking differences between ON parasol and ON smooth monostratified RGCs arise from distinct strategies for sampling a common set of synaptic inputs.

Molecular identification of wide-field amacrine cells in mouse retina that encode stimulus orientation

Visual information processing is sculpted by a diverse group of inhibitory interneurons in the retina called amacrine cells. Yet, for most of the >60 amacrine cell types, molecular identities and specialized functional attributes remain elusive. Here, we developed an intersectional genetic strategy to target a group of wide-field amacrine cells (WACs) in mouse retina that co-express the transcription factor Bhlhe22 and the Kappa Opioid Receptor (KOR; B/K WACs). B/K WACs feature straight, unbranched dendrites spanning over 0.5 mm (∼15° visual angle) and produce non-spiking responses to either light increments or decrements. Two-photon dendritic population imaging reveals Ca 2+ signals tuned to the physical orientations of B/K WAC dendrites, signifying a robust structure-function alignment. B/K WACs establish divergent connections with multiple retinal neurons, including unexpected connections with non-orientation-tuned ganglion cells and bipolar cells. Our work sets the stage for future comprehensive investigations of the most enigmatic group of retinal neurons: WACs.

Heterogeneous presynaptic receptive fields contribute to directional tuning in starburst amacrine cells

The processing of visual information by retinal starburst amacrine cells (SACs) involves transforming excitatory input from bipolar cells (BCs) into directional calcium output. While previous studies have suggested that an asymmetry in the kinetic properties of BCs along the soma-dendritic axes of the postsynaptic cell could enhance directional tuning at the level of individual branches, it remains unclear whether biologically relevant presynaptic kinetics contribute to direction selectivity (DS) when visual stimulation engages the entire dendritic tree. To address this question, we built multicompartmental models of the bipolar-SAC circuit and trained them to boost directional tuning. We report that despite significant dendritic crosstalk and dissimilar directional preferences along the dendrites that occur during whole-cell stimulation, the rules that guide BC kinetics leading to optimal DS are similar to the single-dendrite condition. To correlate model predictions to empirical findings, we utilized two-photon glutamate imaging to study the dynamics of bipolar release onto ON- and OFF-starburst dendrites in the murine retina. We reveal diverse presynaptic dynamics in response to motion in both BC populations; algorithms trained on the experimental data suggested that the differences in the temporal release kinetics are likely to correspond to heterogeneous receptive field properties among the different BC types, including the spatial extent of the center and surround components. In addition, we demonstrate that circuit architecture composed of presynaptic units with experimentally recorded dynamics could enhance directional drive but not to levels that replicate empirical findings, suggesting other DS mechanisms are required to explain SAC function. Our study provides new insights into the complex mechanisms underlying DS in retinal processing and highlights the potential contribution of presynaptic kinetics to the computation of visual information by SACs.

Evolutionary and Developmental Specialization of Foveal Cell Types in the Marmoset

In primates, high-acuity vision is mediated by the fovea, a small specialized central region of the retina. The fovea, unique to the anthropoid lineage among mammals, undergoes notable neuronal morphological changes during postnatal maturation. However, the extent of cellular similarity across anthropoid foveas and the molecular underpinnings of foveal maturation remain unclear. Here, we used high throughput single cell RNA sequencing to profile retinal cells of the common marmoset ( Callithrix jacchus ), an early divergent in anthropoid evolution from humans, apes, and macaques. We generated atlases of the marmoset fovea and peripheral retina for both neonates and adults. Our comparative analysis revealed that marmosets share almost all its foveal types with both humans and macaques, highlighting a conserved cellular structure among primate foveas. Furthermore, by tracing the developmental trajectory of cell types in the foveal and peripheral retina, we found distinct maturation paths for each. In-depth analysis of gene expression differences demonstrated that cone photoreceptors and Müller glia, among others, show the greatest molecular divergence between these two regions. Utilizing single-cell ATAC-seq and gene-regulatory network inference, we uncovered distinct transcriptional regulations differentiating foveal cones from their peripheral counterparts. Further analysis of predicted ligand-receptor interactions suggested a potential role for Müller glia in supporting the maturation of foveal cones. Together, these results provide valuable insights into foveal development, structure, and evolution.

Significance statement

The sharpness of our eyesight hinges on a tiny retinal region known as the fovea. The fovea is pivotal for primate vision and is susceptible to diseases like age-related macular degeneration. We studied the fovea in the marmoset-a primate with ancient evolutionary ties. Our data illustrated the cellular and molecular composition of its fovea across different developmental ages. Our findings highlighted a profound cellular consistency among marmosets, humans, and macaques, emphasizing the value of marmosets in visual research and the study of visual diseases.

Nfia Is Critical for AII Amacrine Cell Production: Selective Bipolar Cell Dependencies and Diminished ERG

The nuclear factor one (NFI) transcription factor genes Nfia, Nfib, and Nfix are all enriched in late-stage retinal progenitor cells, and their loss has been shown to retain these progenitors at the expense of later-generated retinal cell types. Whether they play any role in the specification of those later-generated fates is unknown, but the expression of one of these, Nfia, in a specific amacrine cell type may intimate such a role. Here, Nfia conditional knockout (Nfia-CKO) mice (both sexes) were assessed, finding a massive and largely selective absence of AII amacrine cells. There was, however, a partial reduction in type 2 cone bipolar cells (CBCs), being richly interconnected to AII cells. Counts of dying cells showed a significant increase in Nfia-CKO retinas at postnatal day (P)7, after AII cell numbers were already reduced but in advance of the loss of type 2 CBCs detected by P10. Those results suggest a role for Nfia in the specification of the AII amacrine cell fate and a dependency of the type 2 CBCs on them. Delaying the conditional loss of Nfia to the first postnatal week did not alter AII cell number nor differentiation, further suggesting that its role in AII cells is solely associated with their production. The physiological consequences of their loss were assessed using the ERG, finding the oscillatory potentials to be profoundly diminished. A slight reduction in the b-wave was also detected, attributed to an altered distribution of the terminals of rod bipolar cells, implicating a role of the AII amacrine cells in constraining their stratification.SIGNIFICANCE STATEMENT The transcription factor NFIA is shown to play a critical role in the specification of a single type of retinal amacrine cell, the AII cell. Using an Nfia-conditional knockout mouse to eliminate this population of retinal neurons, we demonstrate two selective bipolar cell dependencies on the AII cells; the terminals of rod bipolar cells become mis-stratified in the inner plexiform layer, and one type of cone bipolar cell undergoes enhanced cell death. The physiological consequence of this loss of the AII cells was also assessed, finding the cells to be a major contributor to the oscillatory potentials in the electroretinogram.

Evolution of neuronal cell classes and types in the vertebrate retina

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

Modular interneuron circuits control motion sensitivity in the mouse retina

Neural computations arise from highly precise connections between specific types of neurons. Retinal ganglion cells (RGCs) with similar stratification patterns are positioned to receive similar inputs but often display different response properties. In this study, we used intersectional mouse genetics to achieve single-cell type labeling and identified an object motion sensitive (OMS) AC type, COMS-AC(counter-OMS AC). Optogenetic stimulation revealed that COMS-AC makes glycinergic synapses with the OMS-insensitive HD2p-RGC, while chemogenetic inactivation showed that COMS-AC provides inhibitory control to HD2p-RGC during local motion. This local inhibition, combined with the inhibitory drive from TH2-AC during global motion, explains the OMS-insensitive feature of HD2p-RGC. In contrast, COMS-AC fails to make synapses with W3(UHD)-RGC, allowing it to exhibit OMS under the control of VGlut3-AC and TH2-AC. These findings reveal modular interneuron circuits that endow structurally similar RGC types with different responses and present a mechanism for redundancy-reduction in the retina to expand coding capacity.

Suppressing DNMT3a Alleviates the Intrinsic Epigenetic Barrier for Optic Nerve Regeneration and Restores Vision in Adult Mice

The limited regenerative potential of the optic nerve in adult mammals presents a major challenge for restoring vision after optic nerve trauma or disease. The mechanisms of this regenerative failure are not fully understood1,2. Here, through small-molecule and genetic screening for epigenetic modulators3, we identify DNA methyltransferase 3a (DNMT3a) as a potent inhibitor of axon regeneration in mouse and human retinal explants. Selective suppression of DNMT3a in retinal ganglion cells (RGCs) by gene targeting or delivery of shRNA leads to robust, full-length regeneration of RGC axons through the optic nerve and restoration of vision in adult mice after nerve crush injury. Genome-wide bisulfite and transcriptome profiling in combination with single nucleus RNA-sequencing of RGCs revealed selective DNA demethylation and reactivation of genetic programs supporting neuronal survival and axonal growth/regeneration by DNMT3a deficiency. This was accompanied by the suppression of gene networks associated with apoptosis and inflammation. Our results identify DNMT3a as the central orchestrator of an RGC-intrinsic mechanism that limits optic nerve regeneration. Suppressing DNMT3a expression in RGCs unlocks the epigenetic switch for optic nerve regeneration and presents a promising therapeutic avenue for effectively reversing vision loss resulted from optic nerve trauma or diseases.

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.

Awake responses suggest inefficient dense coding in the mouse retina

The structure and function of the vertebrate retina have been extensively studied across species with an isolated, ex vivo preparation. Retinal function in vivo, however, remains elusive, especially in awake animals. Here, we performed single-unit extracellular recordings in the optic tract of head-fixed mice to compare the output of awake, anesthetized, and ex vivo retinas. While the visual response properties were overall similar across conditions, we found that awake retinal output had in general (1) faster kinetics with less variability in the response latencies; (2) a larger dynamic range; and (3) higher firing activity, by ~20 Hz on average, for both baseline and visually evoked responses. Our modeling analyses further showed that such awake response patterns convey comparable total information but less efficiently, and allow for a linear population decoder to perform significantly better than the anesthetized or ex vivo responses. These results highlight distinct retinal behavior in awake states, in particular suggesting that the retina employs dense coding in vivo, rather than sparse efficient coding as has been often assumed from ex vivo studies.