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.

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.

Molecular Classification and Comparative Taxonomics of Foveal and Peripheral Cells in Primate Retina.

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.

T Helper Cell Cytokines Modulate Intestinal Stem Cell Renewal and Differentiation.

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.

Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics.

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.

Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets.

Cells, the basic units of biological structure and function, vary broadly in type and state. Single-cell genomics can characterize cell identity and function, but limitations of ease and scale have prevented its broad application. Here we describe Drop-seq, a strategy for quickly profiling thousands of individual cells by separating them into nanoliter-sized aqueous droplets, associating a different barcode with each cell's RNAs, and sequencing them all together. Drop-seq analyzes mRNA transcripts from thousands of individual cells simultaneously while remembering transcripts' cell of origin. We analyzed transcriptomes from 44,808 mouse retinal cells and identified 39 transcriptionally distinct cell populations, creating a molecular atlas of gene expression for known retinal cell classes and novel candidate cell subtypes. Drop-seq will accelerate biological discovery by enabling routine transcriptional profiling at single-cell resolution. VIDEO ABSTRACT.

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.

Hydrophobic CDR3 residues promote the development of self-reactive T cells.

Studies of individual T cell antigen receptors (TCRs) have shed some light on structural features that underlie self-reactivity. However, the general rules that can be used to predict whether TCRs are self-reactive have not been fully elucidated. Here we found that the interfacial hydrophobicity of amino acids at positions 6 and 7 of the complementarity-determining region CDR3ß robustly promoted the development of self-reactive TCRs. This property was found irrespective of the member of the ß-chain variable region (Vß) family present in the TCR or the length of the CDR3ß. An index based on these findings distinguished Vß2(+), Vß6(+) and Vß8.2(+) regulatory T cells from conventional T cells and also distinguished CD4(+) T cells selected by the major histocompatibility complex (MHC) class II molecule I-A(g7) (associated with the development of type 1 diabetes in NOD mice) from those selected by a non-autoimmunity-promoting MHC class II molecule I-A(b). Our results provide a means for distinguishing normal T cell repertoires versus autoimmunity-prone T cell repertoires.

Magnitude and Kinetics of CD8+ T Cell Activation during Hyperacute HIV Infection Impact Viral Set Point.

CD8(+) T cells contribute to the control of HIV, but it is not clear whether initial immune responses modulate the viral set point. We screened high-risk uninfected women twice a week for plasma HIV RNA and identified 12 hyperacute infections. Onset of viremia elicited a massive HIV-specific CD8(+) T cell response, with limited bystander activation of non-HIV memory CD8(+) T cells. HIV-specific CD8(+) T cells secreted little interferon-γ, underwent rapid apoptosis, and failed to upregulate the interleukin-7 receptor, known to be important for T cell survival. The rapidity to peak CD8(+) T cell activation and the absolute magnitude of activation induced by the exponential rise in viremia were inversely correlated with set point viremia. These data indicate that rapid, high magnitude HIV-induced CD8(+) T cell responses are crucial for subsequent immune control of acute infection, which has important implications for HIV vaccine design.

A single-cell survey of the small intestinal epithelium.

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.

Cell atlas of aqueous humor outflow pathways in eyes of humans and four model species provides insight into glaucoma pathogenesis.

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.

Massively parallel single-nucleus RNA-seq with DroNc-seq.

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.

Continuous immunotypes describe human immune variation and predict diverse responses.

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.

Therapeutic efficacy of potent neutralizing HIV-1-specific monoclonal antibodies in SHIV-infected rhesus monkeys.

Human immunodeficiency virus type 1 (HIV-1)-specific monoclonal antibodies with extraordinary potency and breadth have recently been described. In humanized mice, combinations of monoclonal antibodies have been shown to suppress viraemia, but the therapeutic potential of these monoclonal antibodies has not yet been evaluated in primates with an intact immune system. Here we show that administration of a cocktail of HIV-1-specific monoclonal antibodies, as well as the single glycan-dependent monoclonal antibody PGT121, resulted in a rapid and precipitous decline of plasma viraemia to undetectable levels in rhesus monkeys chronically infected with the pathogenic simian-human immunodeficiency virus SHIV-SF162P3. A single monoclonal antibody infusion afforded up to a 3.1 log decline of plasma viral RNA in 7 days and also reduced proviral DNA in peripheral blood, gastrointestinal mucosa and lymph nodes without the development of viral resistance. Moreover, after monoclonal antibody administration, host Gag-specific T-lymphocyte responses showed improved functionality. Virus rebounded in most animals after a median of 56 days when serum monoclonal antibody titres had declined to undetectable levels, although, notably, a subset of animals maintained long-term virological control in the absence of further monoclonal antibody infusions. These data demonstrate a profound therapeutic effect of potent neutralizing HIV-1-specific monoclonal antibodies in SHIV-infected rhesus monkeys as well as an impact on host immune responses. Our findings strongly encourage the investigation of monoclonal antibody therapy for HIV-1 in humans.

Automatic Classification of Cellular Expression by Nonlinear Stochastic Embedding (ACCENSE).

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.

Statistical linkage analysis of substitutions in patient-derived sequences of genotype 1a hepatitis C virus nonstructural protein 3 exposes targets for immunogen design.

UNLABELLED: Chronic hepatitis C virus (HCV) infection is one of the leading causes of liver failure and liver cancer, affecting around 3% of the world's population. The extreme sequence variability of the virus resulting from error-prone replication has thwarted the discovery of a universal prophylactic vaccine. It is known that vigorous and multispecific cellular immune responses, involving both helper CD4(+) and cytotoxic CD8(+) T cells, are associated with the spontaneous clearance of acute HCV infection. Escape mutations in viral epitopes can, however, abrogate protective T-cell responses, leading to viral persistence and associated pathologies. Despite the propensity of the virus to mutate, there might still exist substitutions that incur a fitness cost. In this paper, we identify groups of coevolving residues within HCV nonstructural protein 3 (NS3) by analyzing diverse sequences of this protein using ideas from random matrix theory and associated methods. Our analyses indicate that one of these groups comprises a large percentage of residues for which HCV appears to resist multiple simultaneous substitutions. Targeting multiple residues in this group through vaccine-induced immune responses should either lead to viral recognition or elicit escape substitutions that compromise viral fitness. Our predictions are supported by published clinical data, which suggested that immune genotypes associated with spontaneous clearance of HCV preferentially recognized and targeted this vulnerable group of residues. Moreover, mapping the sites of this group onto the available protein structure provided insight into its functional significance. An epitope-based immunogen is proposed as an alternative to the NS3 epitopes in the peptide-based vaccine IC41. IMPORTANCE: Despite much experimental work on HCV, a thorough statistical study of the HCV sequences for the purpose of immunogen design was missing in the literature. Such a study is vital to identify epistatic couplings among residues that can provide useful insights for designing a potent vaccine. In this work, ideas from random matrix theory were applied to characterize the statistics of substitutions within the diverse publicly available sequences of the genotype 1a HCV NS3 protein, leading to a group of sites for which HCV appears to resist simultaneous substitutions possibly due to deleterious effect on viral fitness. Our analysis leads to completely novel immunogen designs for HCV. In addition, the NS3 epitopes used in the recently proposed peptide-based vaccine IC41 were analyzed in the context of our framework. Our analysis predicts that alternative NS3 epitopes may be worth exploring as they might be more efficacious.

Coordinate linkage of HIV evolution reveals regions of immunological vulnerability.

Cellular immune control of HIV is mediated, in part, by induction of single amino acid mutations that reduce viral fitness, but compensatory mutations limit this effect. Here, we sought to determine if higher order constraints on viral evolution exist, because some coordinately linked combinations of mutations may hurt viability. Immune targeting of multiple sites in such a multidimensionally conserved region might render the virus particularly vulnerable, because viable escape pathways would be greatly restricted. We analyzed available HIV sequences using a method from physics to reveal distinct groups of amino acids whose mutations are collectively coordinated ("HIV sectors"). From the standpoint of mutations at individual sites, one such group in Gag is as conserved as other collectively coevolving groups of sites in Gag. However, it exhibits higher order conservation indicating constraints on the viability of viral strains with multiple mutations. Mapping amino acids from this group onto protein structures shows that combined mutations likely destabilize multiprotein structural interactions critical for viral function. Persons who durably control HIV without medications preferentially target the sector in Gag predicted to be most vulnerable. By sequencing circulating viruses from these individuals, we find that individual mutations occur with similar frequency in this sector as in other targeted Gag sectors. However, multiple mutations within this sector are very rare, indicating previously unrecognized multidimensional constraints on HIV evolution. Targeting such regions with higher order evolutionary constraints provides a novel approach to immunogen design for a vaccine against HIV and other rapidly mutating viruses.

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.

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.