Molecular states underlying neuronal cell type development and plasticity in the whisker cortex.

Mouse whisker somatosensory cortex (wS1) is a major model system to study the experience-dependent plasticity of cortical neuron physiology, morphology, and sensory coding. However, the role of sensory experience in regulating neuronal cell type development and gene expression in wS1 remains poorly understood. We assembled and annotated a transcriptomic atlas of wS1 during postnatal development comprising 45 molecularly distinct neuronal types that can be grouped into eight excitatory and four inhibitory neuron subclasses. Using this atlas, we examined the influence of whisker experience from postnatal day (P) 12, the onset of active whisking, to P22, on the maturation of molecularly distinct cell types. During this developmental period, when whisker experience was normal, ~250 genes were regulated in a neuronal subclass-specific fashion. At the resolution of neuronal types, we found that only the composition of layer (L) 2/3 glutamatergic neuronal types, but not other neuronal types, changed substantially between P12 and P22. These compositional changes resemble those observed previously in the primary visual cortex (V1), and the temporal gene expression changes were also highly conserved between the two regions. In contrast to V1, however, cell type maturation in wS1 is not substantially dependent on sensory experience, as 10-day full-face whisker deprivation did not influence the transcriptomic identity and composition of L2/3 neuronal types. A one-day competitive whisker deprivation protocol also did not affect cell type identity but induced moderate changes in plasticity-related gene expression. Thus, developmental maturation of cell types is similar in V1 and wS1, but sensory deprivation minimally affects cell type development in wS1.

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.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions. The specialized downstream circuit for rod signalling, called the primary rod pathway, is well characterized in mammals, but circuitry for rod signalling in non-mammals is largely unknown. Here we demonstrate that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA sequencing, we identified two bipolar cell types in zebrafish that are related to mammalian rod bipolar cell (RBCs), the only bipolar type that directly carries rod signals from the outer to the inner retina in the primary rod pathway. By combining electrophysiology, histology and ultrastructural reconstruction of the zebrafish RBCs, we found that, similar to mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells postsynaptic to one RBC type is strikingly similar to that of mammalian RBCs and their amacrine partners, suggesting that the cell types and circuit design of the primary rod pathway emerged before the divergence of teleost fish and mammals. The second RBC type, which forms separate pathways, was either lost in mammals or emerged in fish.

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex.

How early sensory experience during "critical periods" of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent (Tan et al., Neuron, 108(4), 2020; Cheng et al., Cell, 185(2), 2022). Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial-ventricular axis in V1. By applying multi-tasking theory (Adler et al., Cell Systems, 8, 2019), we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types ("A", "B", and "C"). By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell type-specific gene expression programs. By contrast, the second program impacts L2/3 cell type identity, regulating a subset of cell type-specific genes and shifting the distribution of cells within the L2/3 manifold, with a depression of the B-type and C-type and a gain of the A-type. Through an integrated analysis of spatial transcriptomic measurements with single-nucleus RNA-seq data from our previous study, we describe how vision patterns L2/3 cortical cell types during the postnatal critical period. Significance statement: Layer 2/3 (L2/3) glutamatergic neurons are important sites of experience-dependent plasticity and learning in the mammalian cortex. Their properties vary continuously with cortical depth and depend upon experience. Here, by applying spatial transcriptomics and different computational approaches in the mouse primary visual cortex, we show that vision regulates orthogonal gene expression programs underlying cell states and cell types. Visual deprivation not only induces an activity-dependent cell state, but also alters the composition of L2/3 cell types, which are appropriately described as a transcriptomic continuum. Our results provide insights into how experience shapes transcriptomes that may, in turn, sculpt brain wiring, function, and behavior.

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.

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.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage. Thus, it has been long assumed that the primary rod pathway evolved in mammals. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs, both zebrafish RBC types connect with all rods in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs, suggesting that the cell types and circuit design of the primary rod pathway have emerged before the divergence of teleost fish and amniotes. The second RBC type, which forms separate pathways, is either lost in mammals or emerged in fish.

Ancient origin of the rod bipolar cell pathway in the vertebrate retina.

Vertebrates rely on rod photoreceptors for vision in low-light conditions1. Mammals have a specialized downstream circuit for rod signaling called the primary rod pathway, which comprises specific cell types and wiring patterns that are thought to be unique to this lineage2-6. Thus, it has been long assumed that the primary rod pathway evolved in mammals3,5-7. Here, we challenge this view by demonstrating that the mammalian primary rod pathway is conserved in zebrafish, which diverged from extant mammals ~400 million years ago. Using single-cell RNA-sequencing, we identified two bipolar cell (BC) types in zebrafish that are related to mammalian rod BCs (RBCs) of the primary rod pathway. By combining electrophysiology, histology, and ultrastructural reconstruction of the zebrafish RBCs, we found that, like mammalian RBCs8, both zebrafish RBC types connect with all rods and red-cones in their dendritic territory, and provide output largely onto amacrine cells. The wiring pattern of the amacrine cells post-synaptic to one RBC type is strikingly similar to that of mammalian RBCs. This suggests that the cell types and circuit design of the primary rod pathway may have emerged before the divergence of teleost fish and amniotes (mammals, bird, reptiles). The second RBC type in zebrafish, which forms separate pathways from the first RBC type, is either lost in mammals or emerged in fish to serve yet unknown roles.

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 needs (Baden et al., 2020). One might expect that retinal cell types 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 teleost fish, a bird, a reptile and a lamprey. Molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells [RGCs] and Muller glia) is striking, with transcriptomic differences across species correlated with evolutionary distance. Major subclasses are also conserved, whereas variation among types within classes or subclasses is more pronounced. However, an integrative analysis revealed that numerous types are shared across species based on conserved gene expression programs that likely trace back to the common ancestor of jawed vertebrates. The degree of variation among types increases from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified mammalian orthologs of midget RGCs, which comprise >80% of RGCs in the human retina, subserve high-acuity vision, and were believed to be primate-specific (Berson, 2008); in contrast, the mouse orthologs comprise <2% of mouse RGCs. Projections 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 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.

Defining Selective Neuronal Resilience and Identifying Targets of Neuroprotection and Axon Regeneration Using Single-Cell RNA Sequencing: Computational Approaches.

We describe a computational workflow to analyze single-cell RNA-sequencing (scRNA-seq) profiles of axotomized retinal ganglion cells (RGCs) in mice. Our goal is to identify differences in the dynamics of survival among 46 molecularly defined RGC types together with molecular signatures that correlate with these differences. The data consists of scRNA-seq profiles of RGCs collected at six time points following optic nerve crush (ONC) (see companion chapter by Jacobi and Tran). We use a supervised classification-based approach to map injured RGCs to type identities and quantify type-specific differences in survival at 2 weeks post crush. As injury-related changes in gene expression confound the inference of type identity in surviving cells, the approach deconvolves type-specific gene signatures from injury responses by using an iterative strategy that leverages measurements along the time course. We use these classifications to compare expression differences between resilient and susceptible subpopulations, identifying potential mediators of resilience. The conceptual framework underlying the method is sufficiently general for analysis of selective vulnerability in other neuronal systems.

Temporal single-cell atlas of non-neuronal retinal cells reveals dynamic, coordinated multicellular responses to central nervous system injury.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2+ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Vision-Dependent and -Independent Molecular Maturation of Mouse Retinal Ganglion Cells.

The development and connectivity of retinal ganglion cells (RGCs), the retina's sole output neurons, are patterned by activity-independent transcriptional programs and activity-dependent remodeling. To inventory the molecular correlates of these influences, we applied high-throughput single-cell RNA sequencing (scRNA-seq) to mouse RGCs at six embryonic and postnatal ages. We identified temporally regulated modules of genes that correlate with, and likely regulate, multiple phases of RGC development, ranging from differentiation and axon guidance to synaptic recognition and refinement. Some of these genes are expressed broadly while others, including key transcription factors and recognition molecules, are selectively expressed by one or a few of the 45 transcriptomically distinct types defined previously in adult mice. Next, we used these results as a foundation to analyze the transcriptomes of RGCs in mice lacking visual experience due to dark rearing from birth or to mutations that ablate either bipolar or photoreceptor cells. 98.5% of visually deprived (VD) RGCs could be unequivocally assigned to a single RGC type based on their transcriptional profiles, demonstrating that visual activity is dispensable for acquisition and maintenance of RGC type identity. However, visual deprivation significantly reduced the transcriptomic distinctions among RGC types, implying that activity is required for complete RGC maturation or maintenance. Consistent with this notion, transcriptomic alternations in VD RGCs significantly overlapped with gene modules found in developing RGCs. Our results provide a resource for mechanistic analyses of RGC differentiation and maturation, and for investigating the role of activity in these processes.

Unified classification of mouse retinal ganglion cells using function, morphology, and gene expression.

Classification and characterization of neuronal types are critical for understanding their function and dysfunction. Neuronal classification schemes typically rely on measurements of electrophysiological, morphological, and molecular features, but aligning such datasets has been challenging. Here, we present a unified classification of mouse retinal ganglion cells (RGCs), the sole retinal output neurons. We use visually evoked responses to classify 1,859 mouse RGCs into 42 types. We also obtain morphological or transcriptomic data from subsets and use these measurements to align the functional classification to publicly available morphological and transcriptomic datasets. We create an online database that allows users to browse or download the data and to classify RGCs from their light responses using a machine learning algorithm. This work provides a resource for studies of RGCs, their upstream circuits in the retina, and their projections in the brain, and establishes a framework for future efforts in neuronal classification and open data distribution.

Solving neurodegeneration: common mechanisms and strategies for new treatments.

Across neurodegenerative diseases, common mechanisms may reveal novel therapeutic targets based on neuronal protection, repair, or regeneration, independent of etiology or site of disease pathology. To address these mechanisms and discuss emerging treatments, in April, 2021, Glaucoma Research Foundation, BrightFocus Foundation, and the Melza M. and Frank Theodore Barr Foundation collaborated to bring together key opinion leaders and experts in the field of neurodegenerative disease for a virtual meeting titled "Solving Neurodegeneration". This "think-tank" style meeting focused on uncovering common mechanistic roots of neurodegenerative disease and promising targets for new treatments, catalyzed by the goal of finding new treatments for glaucoma, the world's leading cause of irreversible blindness and the common interest of the three hosting foundations. Glaucoma, which causes vision loss through degeneration of the optic nerve, likely shares early cellular and molecular events with other neurodegenerative diseases of the central nervous system. Here we discuss major areas of mechanistic overlap between neurodegenerative diseases of the central nervous system: neuroinflammation, bioenergetics and metabolism, genetic contributions, and neurovascular interactions. We summarize important discussion points with emphasis on the research areas that are most innovative and promising in the treatment of neurodegeneration yet require further development. The research that is highlighted provides unique opportunities for collaboration that will lead to efforts in preventing neurodegeneration and ultimately vision loss.

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.

Vision-dependent specification of cell types and function in the developing cortex.

The role of postnatal experience in sculpting cortical circuitry, while long appreciated, is poorly understood at the level of cell types. We explore this in the mouse primary visual cortex (V1) using single-nucleus RNA sequencing, visual deprivation, genetics, and functional imaging. We find that vision selectively drives the specification of glutamatergic cell types in upper layers (L) (L2/3/4), while deeper-layer glutamatergic, GABAergic, and non-neuronal cell types are established prior to eye opening. L2/3 cell types form an experience-dependent spatial continuum defined by the graded expression of ∼200 genes, including regulators of cell adhesion and synapse formation. One of these genes, Igsf9b, a vision-dependent gene encoding an inhibitory synaptic cell adhesion molecule, is required for the normal development of binocular responses in L2/3. In summary, vision preferentially regulates the development of upper-layer glutamatergic cell types through the regulation of cell-type-specific gene expression programs.

Dietary suppression of MHC class II expression in intestinal epithelial cells enhances intestinal tumorigenesis.

Little is known about how interactions of diet, intestinal stem cells (ISCs), and immune cells affect early-stage intestinal tumorigenesis. We show that a high-fat diet (HFD) reduces the expression of the major histocompatibility complex class II (MHC class II) genes in intestinal epithelial cells, including ISCs. This decline in epithelial MHC class II expression in a HFD correlates with reduced intestinal microbiome diversity. Microbial community transfer experiments suggest that epithelial MHC class II expression is regulated by intestinal flora. Mechanistically, pattern recognition receptor (PRR) and interferon-gamma (IFNγ) signaling regulates epithelial MHC class II expression. MHC class II-negative (MHC-II-) ISCs exhibit greater tumor-initiating capacity than their MHC class II-positive (MHC-II+) counterparts upon loss of the tumor suppressor Apc coupled with a HFD, suggesting a role for epithelial MHC class II-mediated immune surveillance in suppressing tumorigenesis. ISC-specific genetic ablation of MHC class II increases tumor burden cell autonomously. Thus, HFD perturbs a microbiome-stem cell-immune cell interaction that contributes to tumor initiation in the intestine.

Generating and Using Transcriptomically Based Retinal Cell Atlases.

It has been known for over a century that the basic organization of the retina is conserved across vertebrates. It has been equally clear that retinal cells can be classified into numerous types, but only recently have methods been devised to explore this diversity in unbiased, scalable, and comprehensive ways. Advances in high-throughput single-cell RNA sequencing (scRNA-seq) have played a pivotal role in this effort. In this article, we outline the experimental and computational components of scRNA-seq and review studies that have used them to generate retinal atlases of cell types in several vertebrate species. These atlases have enabled studies of retinal development, responses of retinal cells to injury, expression patterns of genes implicated in retinal disease, and the evolution of cell types. Recently, the inquiry has expanded to include the entire eye and visual centers in the brain. These studies have enhanced our understanding of retinal function and dysfunction and provided tools and insights for exploring neural diversity throughout the brain.

Molecular classification of zebrafish retinal ganglion cells links genes to cell types to behavior.

Retinal ganglion cells (RGCs) form an array of feature detectors, which convey visual information to central brain regions. Characterizing RGC diversity is required to understand the logic of the underlying functional segregation. Using single-cell transcriptomics, we systematically classified RGCs in adult and larval zebrafish, thereby identifying marker genes for >30 mature types and several developmental intermediates. We used this dataset to engineer transgenic driver lines, enabling specific experimental access to a subset of RGC types. Expression of one or few transcription factors often predicts dendrite morphologies and axonal projections to specific tectal layers and extratectal targets. In vivo calcium imaging revealed that molecularly defined RGCs exhibit specific functional tuning. Finally, chemogenetic ablation of eomesa+ RGCs, which comprise melanopsin-expressing types with projections to a small subset of central targets, selectively impaired phototaxis. Together, our study establishes a framework for systematically studying the functional architecture of the visual system.

Cell Atlas of The Human Fovea and Peripheral Retina.

Most irreversible blindness results from retinal disease. To advance our understanding of the etiology of blinding diseases, we used single-cell RNA-sequencing (scRNA-seq) to analyze the transcriptomes of ~85,000 cells from the fovea and peripheral retina of seven adult human donors. Utilizing computational methods, we identified 58 cell types within 6 classes: photoreceptor, horizontal, bipolar, amacrine, retinal ganglion and non-neuronal cells. Nearly all types are shared between the two retinal regions, but there are notable differences in gene expression and proportions between foveal and peripheral cohorts of shared types. We then used the human retinal atlas to map expression of 636 genes implicated as causes of or risk factors for blinding diseases. Many are expressed in striking cell class-, type-, or region-specific patterns. Finally, we compared gene expression signatures of cell types between human and the cynomolgus macaque monkey, Macaca fascicularis. We show that over 90% of human types correspond transcriptomically to those previously identified in macaque, and that expression of disease-related genes is largely conserved between the two species. These results validate the use of the macaque for modeling blinding disease, and provide a foundation for investigating molecular mechanisms underlying visual processing.