Search

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

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

The 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.

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.

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 Müller 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.

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.

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.

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.

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.

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.

© 2020 National Academy of Sciences. All rights reserved. Increased intraocular pressure (IOP) represents a major risk factor for glaucoma, a prevalent eye disease characterized by death of retinal ganglion cells; lowering IOP is the only proven treatment strategy to delay disease progression. The main determinant of IOP is the equilibrium between production and drainage of aqueous humor, with compromised drainage generally viewed as the primary contributor to dangerous IOP elevations. Drainage occurs through two pathways in the anterior segment of the eye called conventional and uveoscleral. To gain insights into the cell types that comprise these pathways, we used high-throughput single-cell RNA sequencing (scRNAseq). From ∼24,000 single-cell transcriptomes, we identified 19 cell types with molecular markers for each and used histological methods to localize each type. We then performed similar analyses on four organisms used for experimental studies of IOP dynamics and glaucoma: cynomolgus macaque (Macaca fascicularis), rhesus macaque (Macaca mulatta), pig (Sus scrofa), and mouse (Mus musculus). Many human cell types had counterparts in these models, but differences in cell types and gene expression were evident. Finally, we identified the cell types that express genes implicated in glaucoma in all five species. Together, our results provide foundations for investigating the pathogenesis of glaucoma and for using model systems to assess mechanisms and potential interventions.

Unprecedented technological advances in single-cell RNA-sequencing (scRNA-seq) technology have now made it possible to profile genome-wide expression in single cells at low cost and high throughput. There is substantial ongoing effort to use scRNA-seq measurements to identify the "cell types" that form components of a complex tissue, akin to taxonomizing species in ecology. Cell type classification from scRNA-seq data involves the application of computational tools rooted in dimensionality reduction and clustering, and statistical analysis to identify molecular signatures that are unique to each type. As datasets continue to grow in size and complexity, computational challenges abound, requiring analytical methods to be scalable, flexible, and robust. Moreover, careful consideration needs to be paid to experimental biases and statistical challenges that are unique to these measurements to avoid artifacts. This chapter introduces these topics in the context of cell-type identification, and outlines an instructive step-by-step example bioinformatic pipeline for researchers entering this field.

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.

© 2020, The Author(s). 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.

© 2020 National Academy of Sciences. All rights reserved. Increased intraocular pressure (IOP) represents a major risk factor for glaucoma, a prevalent eye disease characterized by death of retinal ganglion cells; lowering IOP is the only proven treatment strategy to delay disease progression. The main determinant of IOP is the equilibrium between production and drainage of aqueous humor, with compromised drainage generally viewed as the primary contributor to dangerous IOP elevations. Drainage occurs through two pathways in the anterior segment of the eye called conventional and uveoscleral. To gain insights into the cell types that comprise these pathways, we used high-throughput single-cell RNA sequencing (scRNAseq). From ∼24,000 single-cell transcriptomes, we identified 19 cell types with molecular markers for each and used histological methods to localize each type. We then performed similar analyses on four organisms used for experimental studies of IOP dynamics and glaucoma: cynomolgus macaque (Macaca fascicularis), rhesus macaque (Macaca mulatta), pig (Sus scrofa), and mouse (Mus musculus). Many human cell types had counterparts in these models, but differences in cell types and gene expression were evident. Finally, we identified the cell types that express genes implicated in glaucoma in all five species. Together, our results provide foundations for investigating the pathogenesis of glaucoma and for using model systems to assess mechanisms and potential interventions.

© 2019 Elsevier Inc. Neuronal types in the central nervous system differ dramatically in their resilience to injury or other insults. Here we studied the selective resilience of mouse retinal ganglion cells (RGCs) following optic nerve crush (ONC), which severs their axons and leads to death of ∼80% of RGCs within 2 weeks. To identify expression programs associated with differential resilience, we first used single-cell RNA-seq (scRNA-seq) to generate a comprehensive molecular atlas of 46 RGC types in adult retina. We then tracked their survival after ONC; characterized transcriptomic, physiological, and morphological changes that preceded degeneration; and identified genes selectively expressed by each type. Finally, using loss- and gain-of-function assays in vivo, we showed that manipulating some of these genes improved neuronal survival and axon regeneration following ONC. This study provides a systematic framework for parsing type-specific responses to injury and demonstrates that differential gene expression can be used to reveal molecular targets for intervention.

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.

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.

Increased intraocular pressure (IOP) represents a major risk factor for glaucoma, a prevalent eye disease characterized by death of retinal ganglion cells; lowering IOP is the only proven treatment strategy to delay disease progression. The main determinant of IOP is the equilibrium between production and drainage of aqueous humor, with compromised drainage generally viewed as the primary contributor to dangerous IOP elevations. Drainage occurs through two pathways in the anterior segment of the eye called conventional and uveoscleral. To gain insights into the cell types that comprise these pathways, we used high-throughput single-cell RNA sequencing (scRNAseq). From ∼24,000 single-cell transcriptomes, we identified 19 cell types with molecular markers for each and used histological methods to localize each type. We then performed similar analyses on four organisms used for experimental studies of IOP dynamics and glaucoma: cynomolgus macaque (Macaca fascicularis), rhesus macaque (Macaca mulatta), pig (Sus scrofa), and mouse (Mus musculus). Many human cell types had counterparts in these models, but differences in cell types and gene expression were evident. Finally, we identified the cell types that express genes implicated in glaucoma in all five species. Together, our results provide foundations for investigating the pathogenesis of glaucoma and for using model systems to assess mechanisms and potential interventions.

In the small intestine, a niche of accessory cell types supports the generation of mature epithelial cell types from intestinal stem cells (ISCs). It is unclear, however, if and how immune cells in the niche affect ISC fate or the balance between self-renewal and differentiation. Here, we use single-cell RNA sequencing (scRNA-seq) to identify MHC class II (MHCII) machinery enrichment in two subsets of Lgr5⁺ ISCs. We show that MHCII⁺ Lgr5⁺ ISCs are non-conventional antigen-presenting cells in co-cultures with CD4⁺ T helper (Th) cells. Stimulation of intestinal organoids with key Th cytokines affects Lgr5⁺ ISC renewal and differentiation in opposing ways: pro-inflammatory signals promote differentiation, while regulatory cells and cytokines reduce it. In vivo genetic perturbation of Th cells or MHCII expression on Lgr5⁺ ISCs impacts epithelial cell differentiation and IEC fate during infection. These interactions between Th cells and Lgr5⁺ ISCs, thus, orchestrate tissue-wide responses to external signals. Intestinal stem cells act as non-conventional antigen presenting cells, and these interactions with T helper cells modulate ISC renewal and differentiation to shape the intestine. Keywords: gut biology; intestinal stem cells; ISCs; T helper; Th; mucosal immunity; MHC class II; MHCII; tuft cells; T regulatory; Treg; single cell RNA-seq; scRNA-seq; epithelial differentiation; stem cell renewal, National Institute of Health (U.S.) (Award 1DP2OD020839), National Institute of Health (U.S.) (Grant CA211184), National Institute of Health (U.S.) (Grant AG045144)

Increased intraocular pressure (IOP) represents a major risk factor for glaucoma, a prevalent eye disease characterized by death of retinal ganglion cells; lowering IOP is the only proven treatment strategy to delay disease progression. The main determinant of IOP is the equilibrium between production and drainage of aqueous humor, with compromised drainage generally viewed as the primary contributor to dangerous IOP elevations. Drainage occurs through two pathways in the anterior segment of the eye called conventional and uveoscleral. To gain insights into the cell types that comprise these pathways, we used high-throughput single-cell RNA sequencing (scRNAseq). From ∼24,000 single-cell transcriptomes, we identified 19 cell types with molecular markers for each and used histological methods to localize each type. We then performed similar analyses on four organisms used for experimental studies of IOP dynamics and glaucoma: cynomolgus macaque (Macaca fascicularis), rhesus macaque (Macaca mulatta), pig (Sus scrofa), and mouse (Mus musculus). Many human cell types had counterparts in these models, but differences in cell types and gene expression were evident. Finally, we identified the cell types that express genes implicated in glaucoma in all five species. Together, our results provide foundations for investigating the pathogenesis of glaucoma and for using model systems to assess mechanisms and potential interventions.

Neuronal types in the central nervous system differ dramatically in their resilience to injury or other insults. Here we studied the selective resilience of mouse retinal ganglion cells (RGCs) following optic nerve crush (ONC), which severs their axons and leads to death of ∼80% of RGCs within 2 weeks. To identify expression programs associated with differential resilience, we first used single-cell RNA-seq (scRNA-seq) to generate a comprehensive molecular atlas of 46 RGC types in adult retina. We then tracked their survival after ONC; characterized transcriptomic, physiological, and morphological changes that preceded degeneration; and identified genes selectively expressed by each type. Finally, using loss- and gain-of-function assays in vivo, we showed that manipulating some of these genes improved neuronal survival and axon regeneration following ONC. This study provides a systematic framework for parsing type-specific responses to injury and demonstrates that differential gene expression can be used to reveal molecular targets for intervention.

High-acuity vision in primates, including humans, is mediated by a small central retinal region called the fovea. As more accessible organisms lack a fovea, its specialized function and its dysfunction in ocular diseases remain poorly understood. We used 165,000 single-cell RNA-seq profiles to generate comprehensive cellular taxonomies of macaque fovea and peripheral retina. More than 80% of >60 cell types match between the two regions but exhibit substantial differences in proportions and gene expression, some of which we relate to functional differences. Comparison of macaque retinal types with those of mice reveals that interneuron types are tightly conserved. In contrast, projection neuron types and programs diverge, despite exhibiting conserved transcription factor codes. Key macaque types are conserved in humans, allowing mapping of cell-type and region-specific expression of >190 genes associated with 7 human retinal diseases. Our work provides a framework for comparative single-cell analysis across tissue regions and species.

High-acuity vision in primates, including humans, is mediated by a small central retinal region called the fovea. As more accessible organisms lack a fovea, its specialized function and its dysfunction in ocular diseases remain poorly understood. We used 165,000 single-cell RNA-seq profiles to generate comprehensive cellular taxonomies of macaque fovea and peripheral retina. More than 80% of >60 cell types match between the two regions but exhibit substantial differences in proportions and gene expression, some of which we relate to functional differences. Comparison of macaque retinal types with those of mice reveals that interneuron types are tightly conserved. In contrast, projection neuron types and programs diverge, despite exhibiting conserved transcription factor codes. Key macaque types are conserved in humans, allowing mapping of cell-type and region-specific expression of >190 genes associated with 7 human retinal diseases. Our work provides a framework for comparative single-cell analysis across tissue regions and species.

Neuronal types in the central nervous system differ dramatically in their resilience to injury or other insults. Here we studied the selective resilience of mouse retinal ganglion cells (RGCs) following optic nerve crush (ONC), which severs their axons and leads to death of ∼80% of RGCs within 2 weeks. To identify expression programs associated with differential resilience, we first used single-cell RNA-seq (scRNA-seq) to generate a comprehensive molecular atlas of 46 RGC types in adult retina. We then tracked their survival after ONC; characterized transcriptomic, physiological, and morphological changes that preceded degeneration; and identified genes selectively expressed by each type. Finally, using loss- and gain-of-function assays in vivo, we showed that manipulating some of these genes improved neuronal survival and axon regeneration following ONC. This study provides a systematic framework for parsing type-specific responses to injury and demonstrates that differential gene expression can be used to reveal molecular targets for intervention.

Patterns of gene expression can be used to characterize and classify neuronal types. It is challenging, however, to generate taxonomies that fulfill the essential criteria of being comprehensive, harmonizing with conventional classification schemes, and lacking superfluous subdivisions of genuine types. To address these challenges, we used massively parallel single-cell RNA profiling and optimized computational methods on a heterogeneous class of neurons, mouse retinal bipolar cells (BCs). From a population of ∼25,000 BCs, we derived a molecular classification that identified 15 types, including all types observed previously and two novel types, one of which has a non-canonical morphology and position. We validated the classification scheme and identified dozens of novel markers using methods that match molecular expression to cell morphology. This work provides a systematic methodology for achieving comprehensive molecular classification of neurons, identifies novel neuronal types, and uncovers transcriptional differences that distinguish types within a class.

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

The identification of cell types and marker genes is critical for dissecting neural development and function, but the size and complexity of the brain has hindered the comprehensive discovery of cell types. We combined single-cell RNA-seq (scRNA-seq) with anatomical brain registration to create a comprehensive map of the zebrafish habenula, a conserved forebrain hub involved in pain processing and learning. Single-cell transcriptomes of ∼13,000 habenular cells with 4× cellular coverage identified 18 neuronal types and dozens of marker genes. Registration of marker genes onto a reference atlas created a resource for anatomical and functional studies and enabled the mapping of active neurons onto neuronal types following aversive stimuli. Strikingly, despite brain growth and functional maturation, cell types were retained between the larval and adult habenula. This study provides a gene expression atlas to dissect habenular development and function and offers a general framework for the comprehensive characterization of other brain regions.

During embryogenesis, cells acquire distinct fates by transitioning through transcriptional states. To uncover these transcriptional trajectories during zebrafish embryogenesis, we sequenced 38,731 cells and developed URD, a simulated diffusion-based computational reconstruction method. URD identified the trajectories of 25 cell types through early somitogenesis, gene expression along them, and their spatial origin in the blastula. Analysis of Nodal signaling mutants revealed that their transcriptomes were canalized into a subset of wild-type transcriptional trajectories. Some wild-type developmental branch points contained cells that express genes characteristic of multiple fates. These cells appeared to trans-specify from one fate to another. These findings reconstruct the transcriptional trajectories of a vertebrate embryo, highlight the concurrent canalization and plasticity of embryonic specification, and provide a framework with which to reconstruct complex developmental trees from single-cell transcriptomes.

In the small intestine, a niche of accessory cell types supports the generation of mature epithelial cell types from intestinal stem cells (ISCs). It is unclear, however, if and how immune cells in the niche affect ISC fate or the balance between self-renewal and differentiation. Here, we use single-cell RNA sequencing (scRNA-seq) to identify MHC class II (MHCII) machinery enrichment in two subsets of Lgr5+ ISCs. We show that MHCII+ Lgr5+ ISCs are non-conventional antigen-presenting cells in co-cultures with CD4+ T helper (Th) cells. Stimulation of intestinal organoids with key Th cytokines affects Lgr5+ ISC renewal and differentiation in opposing ways: pro-inflammatory signals promote differentiation, while regulatory cells and cytokines reduce it. In vivo genetic perturbation of Th cells or MHCII expression on Lgr5+ ISCs impacts epithelial cell differentiation and IEC fate during infection. These interactions between Th cells and Lgr5+ ISCs, thus, orchestrate tissue-wide responses to external signals.

Dendritic cells (DCs) and monocytes play a central role in pathogen sensing, phagocytosis, and antigen presentation and consist of multiple specialized subtypes. However, their identities and interrelationships are not fully understood. Using unbiased single-cell RNA sequencing (RNA-seq) of ~2400 cells, we identified six human DCs and four monocyte subtypes in human blood. Our study reveals a new DC subset that shares properties with plasmacytoid DCs (pDCs) but potently activates T cells, thus redefining pDCs; a new subdivision within the CD1C+ subset of DCs; the relationship between blastic plasmacytoid DC neoplasia cells and healthy DCs; and circulating progenitor of conventional DCs (cDCs). Our revised taxonomy will enable more accurate functional and developmental analyses as well as immune monitoring in health and disease.

The identification of cell types and marker genes is critical for dissecting neural development and function, but the size and complexity of the brain has hindered the comprehensive discovery of cell types. We combined single-cell RNA-seq (scRNA-seq) with anatomical brain registration to create a comprehensive map of the zebrafish habenula, a conserved forebrain hub involved in pain processing and learning. Single-cell transcriptomes of ∼13,000 habenular cells with 4× cellular coverage identified 18 neuronal types and dozens of marker genes. Registration of marker genes onto a reference atlas created a resource for anatomical and functional studies and enabled the mapping of active neurons onto neuronal types following aversive stimuli. Strikingly, despite brain growth and functional maturation, cell types were retained between the larval and adult habenula. This study provides a gene expression atlas to dissect habenular development and function and offers a general framework for the comprehensive characterization of other brain regions.

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

In the small intestine, a niche of accessory cell types supports the generation of mature epithelial cell types from intestinal stem cells (ISCs). It is unclear, however, if and how immune cells in the niche affect ISC fate or the balance between self-renewal and differentiation. Here, we use single-cell RNA sequencing (scRNA-seq) to identify MHC class II (MHCII) machinery enrichment in two subsets of Lgr5+ ISCs. We show that MHCII+ Lgr5+ ISCs are non-conventional antigen-presenting cells in co-cultures with CD4+ T helper (Th) cells. Stimulation of intestinal organoids with key Th cytokines affects Lgr5+ ISC renewal and differentiation in opposing ways: pro-inflammatory signals promote differentiation, while regulatory cells and cytokines reduce it. In vivo genetic perturbation of Th cells or MHCII expression on Lgr5+ ISCs impacts epithelial cell differentiation and IEC fate during infection. These interactions between Th cells and Lgr5+ ISCs, thus, orchestrate tissue-wide responses to external signals.

During embryogenesis, cells acquire distinct fates by transitioning through transcriptional states. To uncover these transcriptional trajectories during zebrafish embryogenesis, we sequenced 38,731 cells and developed URD, a simulated diffusion-based computational reconstruction method. URD identified the trajectories of 25 cell types through early somitogenesis, gene expression along them, and their spatial origin in the blastula. Analysis of Nodal signaling mutants revealed that their transcriptomes were canalized into a subset of wild-type transcriptional trajectories. Some wild-type developmental branch points contained cells that express genes characteristic of multiple fates. These cells appeared to trans-specify from one fate to another. These findings reconstruct the transcriptional trajectories of a vertebrate embryo, highlight the concurrent canalization and plasticity of embryonic specification, and provide a framework with which to reconstruct complex developmental trees from single-cell transcriptomes.

Patterns of gene expression can be used to characterize and classify neuronal types. It is challenging, however, to generate taxonomies that fulfill the essential criteria of being comprehensive, harmonizing with conventional classification schemes, and lacking superfluous subdivisions of genuine types. To address these challenges, we used massively parallel single-cell RNA profiling and optimized computational methods on a heterogeneous class of neurons, mouse retinal bipolar cells (BCs). From a population of ∼25,000 BCs, we derived a molecular classification that identified 15 types, including all types observed previously and two novel types, one of which has a non-canonical morphology and position. We validated the classification scheme and identified dozens of novel markers using methods that match molecular expression to cell morphology. This work provides a systematic methodology for achieving comprehensive molecular classification of neurons, identifies novel neuronal types, and uncovers transcriptional differences that distinguish types within a class.

Intestinal epithelial cells absorb nutrients, respond to microbes, function as a barrier and help to coordinate immune responses. Here we report profiling of 53,193 individual epithelial cells from the small intestine and organoids of mice, which enabled the identification and characterization of previously unknown subtypes of intestinal epithelial cell and their gene signatures. We found unexpected diversity in hormone-secreting enteroendocrine cells and constructed the taxonomy of newly identified subtypes, and distinguished between two subtypes of tuft cell, one of which expresses the epithelial cytokine Tslp and the pan-immune marker CD45, which was not previously associated with non-haematopoietic cells. We also characterized the ways in which cell-intrinsic states and the proportions of different cell types respond to bacterial and helminth infections: Salmonella infection caused an increase in the abundance of Paneth cells and enterocytes, and broad activation of an antimicrobial program; Heligmosomoides polygyrus caused an increase in the abundance of goblet and tuft cells. Our survey highlights previously unidentified markers and programs, associates sensory molecules with cell types, and uncovers principles of gut homeostasis and response to pathogens.

Single-nucleus RNA sequencing (sNuc-seq) profiles RNA from tissues that are preserved or cannot be dissociated, but it does not provide high throughput. Here, we develop DroNc-seq: massively parallel sNuc-seq with droplet technology. We profile 39,111 nuclei from mouse and human archived brain samples to demonstrate sensitive, efficient, and unbiased classification of cell types, paving the way for systematic charting of cell atlases. Keywords: Cellular neuroscience; Gene expression; Gene expression analysis; RNA sequencing

The immune system consists of many specialized cell populations that communicate with each other to achieve systemic immune responses. Our analyses of various measured immune cell population frequencies in healthy humans and their responses to diverse stimuli show that human immune variation is continuous in nature, rather than characterized by discrete groups of similar individuals. We show that the same three key combinations of immune cell population frequencies can define an individual’s immunotype and predict a diverse set of functional responses to cytokine stimulation. We find that, even though interindividual variations in specific cell population frequencies can be large, unrelated individuals of younger age have more homogeneous immunotypes than older individuals. Across age groups, cytomegalovirus seropositive individuals displayed immunotypes characteristic of older individuals. The conceptual framework for defining immunotypes suggested by our results could guide the development of better therapies that appropriately modulate collective immunotypes, rather than individual immune components. Keywords: human immune variation; immune cell composition; systems immunology; aging, National Institutes of Health (U.S.) (Grant R01 HL120724)

Mass cytometry enables an unprecedented number of parameters to be measured in individual cells at a high throughput, but the large dimensionality of the resulting data severely limits approaches relying on manual "gating." Clustering cells based on phenotypic similarity comes at a loss of single-cell resolution and often the number of subpopulations is unknown a priori. Here we describe ACCENSE, a tool that combines nonlinear dimensionality reduction with density-based partitioning, and displays multivariate cellular phenotypes on a 2D plot. We apply ACCENSE to 35-parameter mass cytometry data from CD8(+) T cells derived from specific pathogen-free and germ-free mice, and stratify cells into phenotypic subpopulations. Our results show significant heterogeneity within the known CD8(+) T-cell subpopulations, and of particular note is that we find a large novel subpopulation in both specific pathogen-free and germ-free mice that has not been described previously. This subpopulation possesses a phenotypic signature that is distinct from conventional naive and memory subpopulations when analyzed by ACCENSE, but is not distinguishable on a biaxial plot of standard markers. We are able to automatically identify cellular subpopulations based on all proteins analyzed, thus aiding the full utilization of powerful new single-cell technologies such as mass cytometry.

Mass cytometry enables an unprecedented number of parameters to be measured in individual cells at a high throughput, but the large dimensionality of the resulting data severely limits approaches relying on manual "gating." Clustering cells based on phenotypic similarity comes at a loss of single-cell resolution and often the number of subpopulations is unknown a priori. Here we describe ACCENSE, a tool that combines nonlinear dimensionality reduction with density-based partitioning, and displays multivariate cellular phenotypes on a 2D plot. We apply ACCENSE to 35-parameter mass cytometry data from CD8(+) T cells derived from specific pathogen-free and germ-free mice, and stratify cells into phenotypic subpopulations. Our results show significant heterogeneity within the known CD8(+) T-cell subpopulations, and of particular note is that we find a large novel subpopulation in both specific pathogen-free and germ-free mice that has not been described previously. This subpopulation possesses a phenotypic signature that is distinct from conventional naive and memory subpopulations when analyzed by ACCENSE, but is not distinguishable on a biaxial plot of standard markers. We are able to automatically identify cellular subpopulations based on all proteins analyzed, thus aiding the full utilization of powerful new single-cell technologies such as mass cytometry.

Single-nucleus RNA sequencing (sNuc-seq) profiles RNA from tissues that are preserved or cannot be dissociated, but it does not provide high throughput. Here, we develop DroNc-seq: massively parallel sNuc-seq with droplet technology. We profile 39,111 nuclei from mouse and human archived brain samples to demonstrate sensitive, efficient, and unbiased classification of cell types, paving the way for systematic charting of cell atlases.

Intestinal epithelial cells absorb nutrients, respond to microbes, function as a barrier and help to coordinate immune responses. Here we report profiling of 53,193 individual epithelial cells from the small intestine and organoids of mice, which enabled the identification and characterization of previously unknown subtypes of intestinal epithelial cell and their gene signatures. We found unexpected diversity in hormone-secreting enteroendocrine cells and constructed the taxonomy of newly identified subtypes, and distinguished between two subtypes of tuft cell, one of which expresses the epithelial cytokine Tslp and the pan-immune marker CD45, which was not previously associated with non-haematopoietic cells. We also characterized the ways in which cell-intrinsic states and the proportions of different cell types respond to bacterial and helminth infections: Salmonella infection caused an increase in the abundance of Paneth cells and enterocytes, and broad activation of an antimicrobial program; Heligmosomoides polygyrus caused an increase in the abundance of goblet and tuft cells. Our survey highlights previously unidentified markers and programs, associates sensory molecules with cell types, and uncovers principles of gut homeostasis and response to pathogens.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are a promising platform for cardiac studies in vitro, and possibly for tissue repair in humans. However, hiPSC-CM cells tend to retain morphology, metabolism, patterns of gene expression, and electrophysiology similar to that of embryonic cardiomyocytes. We grew hiPSC-CM in patterned islands of different sizes and shapes, and measured the effect of island geometry on action potential waveform and calcium dynamics using optical recordings of voltage and calcium from 970 islands of different sizes. hiPSC-CM in larger islands showed electrical and calcium dynamics indicative of greater functional maturity. We then compared transcriptional signatures of the small and large islands against a developmental time course of cardiac differentiation. Although island size had little effect on expression of most genes whose levels differed between hiPSC-CM and adult primary CM, we identified a subset of genes for which island size drove the majority (58%) of the changes associated with functional maturation. Finally, we patterned hiPSC-CM on islands with a variety of shapes to probe the relative contributions of soluble factors, electrical coupling, and direct cell-cell contacts to the functional maturation. Collectively, our data show that optical electrophysiology is a powerful tool for assaying hiPSC-CM maturation, and that island size powerfully drives activation of a subset of genes involved in cardiac maturation.

Intestinal epithelial cells absorb nutrients, respond to microbes, function as a barrier and help to coordinate immune responses. Here we report profiling of 53,193 individual epithelial cells from the small intestine and organoids of mice, which enabled the identification and characterization of previously unknown subtypes of intestinal epithelial cell and their gene signatures. We found unexpected diversity in hormone-secreting enteroendocrine cells and constructed the taxonomy of newly identified subtypes, and distinguished between two subtypes of tuft cell, one of which expresses the epithelial cytokine Tslp and the pan-immune marker CD45, which was not previously associated with non-haematopoietic cells. We also characterized the ways in which cell-intrinsic states and the proportions of different cell types respond to bacterial and helminth infections: Salmonella infection caused an increase in the abundance of Paneth cells and enterocytes, and broad activation of an antimicrobial program; Heligmosomoides polygyrus caused an increase in the abundance of goblet and tuft cells. Our survey highlights previously unidentified markers and programs, associates sensory molecules with cell types, and uncovers principles of gut homeostasis and response to pathogens.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are a promising platform for cardiac studies in vitro, and possibly for tissue repair in humans. However, hiPSC-CM cells tend to retain morphology, metabolism, patterns of gene expression, and electrophysiology similar to that of embryonic cardiomyocytes. We grew hiPSC-CM in patterned islands of different sizes and shapes, and measured the effect of island geometry on action potential waveform and calcium dynamics using optical recordings of voltage and calcium from 970 islands of different sizes. hiPSC-CM in larger islands showed electrical and calcium dynamics indicative of greater functional maturity. We then compared transcriptional signatures of the small and large islands against a developmental time course of cardiac differentiation. Although island size had little effect on expression of most genes whose levels differed between hiPSC-CM and adult primary CM, we identified a subset of genes for which island size drove the majority (58%) of the changes associated with functional maturation. Finally, we patterned hiPSC-CM on islands with a variety of shapes to probe the relative contributions of soluble factors, electrical coupling, and direct cell-cell contacts to the functional maturation. Collectively, our data show that optical electrophysiology is a powerful tool for assaying hiPSC-CM maturation, and that island size powerfully drives activation of a subset of genes involved in cardiac maturation.

Dendritic cells (DCs) and monocytes play a central role in pathogen sensing, phagocytosis, and antigen presentation and consist of multiple specialized subtypes. However, their identities and interrelationships are not fully understood. Using unbiased single-cell RNA sequencing (RNA-seq) of ~2400 cells, we identified six human DCs and four monocyte subtypes in human blood. Our study reveals a new DC subset that shares properties with plasmacytoid DCs (pDCs) but potently activates T cells, thus redefining pDCs; a new subdivision within the CD1C+ subset of DCs; the relationship between blastic plasmacytoid DC neoplasia cells and healthy DCs; and circulating progenitor of conventional DCs (cDCs). Our revised taxonomy will enable more accurate functional and developmental analyses as well as immune monitoring in health and disease.

The immune system consists of many specialized cell populations that communicate with each other to achieve systemic immune responses. Our analyses of various measured immune cell population frequencies in healthy humans and their responses to diverse stimuli show that human immune variation is continuous in nature, rather than characterized by discrete groups of similar individuals. We show that the same three key combinations of immune cell population frequencies can define an individual's immunotype and predict a diverse set of functional responses to cytokine stimulation. We find that, even though interindividual variations in specific cell population frequencies can be large, unrelated individuals of younger age have more homogeneous immunotypes than older individuals. Across age groups, cytomegalovirus seropositive individuals displayed immunotypes characteristic of older individuals. The conceptual framework for defining immunotypes suggested by our results could guide the development of better therapies that appropriately modulate collective immunotypes, rather than individual immune components.

Understanding how populations of human T cells leverage cellular heterogeneity, plasticity, and diversity to achieve a wide range of functional flexibility, particularly during dynamic processes such as development, differentiation, and antigenic response, is a core challenge that is well suited for single-cell analysis. Hypothesis-free evaluation of cellular states and subpopulations by transcriptional profiling of single T cells can identify relationships that may be obscured by targeted approaches such as FACS sorting on cell-surface antigens, or bulk expression analysis. While this approach is relevant to all cell types, it is of particular interest in the study of T cells for which classical phenotypic criteria are now viewed as insufficient for distinguishing different T cell subtypes and transitional states, and defining the changes associated with dysfunctional T cell states in autoimmunity and tumor-related exhaustion. This unit describes a protocol to generate single-cell transcriptomic libraries of human blood CD4+ and CD8+ T cells, and also introduces the basic bioinformatic steps to process the resulting sequence data for further computational analysis. We show how cellular subpopulations can be identified from transcriptional data, and derive characteristic gene expression signatures that distinguish these states. We believe single-cell RNA-seq is a powerful technique to study the cellular heterogeneity in complex tissues, a paradigm that will be of great value for the immune system.