The sizes of different neuronal populations within the CNS are precisely controlled, but whether neuronal number is coordinated between cell types is unknown. We examined the covariance structure of 12 different retinal cell types across 30 genetically distinct lines of mice, finding minimal covariation when comparing synaptically connected or developmentally related cell types. Variation mapped to one or more genomic loci for each cell type, but rarely were these shared, indicating minimal genetic coregulation of final number. Multiple genes, therefore, participate in the specification of the size of every population of retinal neuron, yet genetic variants work largely independent of one another during development to modulate those numbers, yielding substantial variability in the convergence ratios between pre- and postsynaptic populations. Density-dependent cellular interactions in the outer plexiform layer overcome this variability to ensure the formation of neuronal circuits that maintain constant retinal coverage and complete afferent sampling.

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 present study interrogated a quantitative trait locus (QTL) on Chr 4 associated with the population sizes of two types of bipolar cell in the mouse retina. This locus was identified by quantifying the number of rod bipolar cells and Type 2 cone bipolar cells across a panel of recombinant inbred (RI) strains of mice derived from two inbred laboratory strains, C57BL/6J (B6/J) and A/J, and mapping a proportion of that variation in cell number, for each cell type, to this shared locus. There, we identified the candidate gene X Kell blood group precursor related family member 8 homolog (Xkr8). While Xkr8 has no documented role in the retina, we localize robust expression in the mature retina via in situ hybridization, confirm its developmental presence via immunolabeling, and show that it is differentially regulated during the postnatal period between the B6/J and A/J strains using qPCR. Microarray analysis, derived from whole eye mRNA from the entire RI strain set, demonstrates significant negative correlation of Xkr8 expression with the number of each of these two types of bipolar cells, and the variation in Xkr8 expression across the strains maps a cis-eQTL, implicating a regulatory variant discriminating the parental genomes. Xkr8 plasmid electroporation during development yielded a reduction in the number of bipolar cells in the retina, while sequence analysis of Xkr8 in the two parental strain genomes identified a structural variant in the 3' UTR that may disrupt mRNA stability, and two SNPs in the promoter that create transcription factor binding sites. We propose that Xkr8, via its participation in mediating cell death, plays a role in the specification of bipolar cell number in the retina.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

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

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

We previously reported that vision specifies Layer 2/3 (L2/3) glutamatergic cell-type identity in the primary visual cortex (V1). Using unsupervised clustering of single-nucleus RNA-sequencing data, we identified molecularly distinct L2/3 cell types in normal-reared (NR) and dark-reared (DR) mice, but the two sets exhibited poor correspondence. Here, we show that classification of cell types was confounded in DR by vision-dependent gene programs that are orthogonal to gene programs underlying cell-type identity. A focused clustering analysis successfully matches cell types between DR and NR, suggesting that cell identity-defining gene programs persist under vision deprivation but are overshadowed by vision-dependent transcriptomic variation. Using multi-tasking theory we show that L2/3 cell types form a continuum between three cell-archetypes. Visual deprivation markedly shifts this distribution along the continuum. Thus, dark-rearing markedly influences cell states thereby masking cell-type-identities and changes the distribution of L2/3 types along a transcriptomic continuum.

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.

Genome-wide association studies (GWAS) have revealed risk alleles for ulcerative colitis (UC). To understand their cell type specificities and pathways of action, we generate an atlas of 366,650 cells from the colon mucosa of 18 UC patients and 12 healthy individuals, revealing 51 epithelial, stromal, and immune cell subsets, including BEST4+ enterocytes, microfold-like cells, and IL13RA2+IL11+ inflammatory fibroblasts, which we associate with resistance to anti-TNF treatment. Inflammatory fibroblasts, inflammatory monocytes, microfold-like cells, and T cells that co-express CD8 and IL-17 expand with disease, forming intercellular interaction hubs. Many UC risk genes are cell type specific and co-regulated within relatively few gene modules, suggesting convergence onto limited sets of cell types and pathways. Using this observation, we nominate and infer functions for specific risk genes across GWAS loci. Our work provides a framework for interrogating complex human diseases and mapping risk variants to cell types and pathways.

Islet-1 (Isl1) is a LIM-homeodomain (LIM-HD) transcription factor that functions in a combinatorial manner with other LIM-HD proteins to direct the differentiation of distinct cell types within the central nervous system and many other tissues. A study of pancreatic cell lines showed that Isl1 is alternatively spliced generating a second isoform, Isl1β, which is missing 23 amino acids within the C-terminal region. This study examines the expression of the canonical and alternative Isl1 transcripts across other tissues, in particular, within the retina, where Isl1 is required for the differentiation of multiple neuronal cell types. The alternative splicing of Isl1 is shown to occur in multiple tissues, but the relative abundance of Isl1α and Isl1β expression varies greatly across them. In most tissues, Isl1α is the more abundant transcript, but in others the transcripts are expressed equally, or the alternative splice variant is dominant. Within the retina, differential expression of the two Isl1 transcripts increases as a function of development, with dynamic changes in expression peaking at E16.5 and again at P10. At the cellular level, individual retinal ganglion cells vary in their expression, with a subset of small-to-medium sized cells expressing only the alternative isoform. The functional significance of the difference in protein sequence between the two Isl1 isoforms was also assessed using a luciferase assay, demonstrating that the alternative isoform forms a less effective transcriptional complex for activating gene expression. These results demonstrate the differential presence of the canonical and alternative isoforms of Isl1 amongst retinal ganglion cell classes. As Isl1 participates in the differentiation of multiple cell types within the CNS, the present results support a role for alternative splicing in the establishment of cellular diversity in the developing nervous system.

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.

Speciation leads to adaptive changes in organ cellular physiology and creates challenges for studying rare cell-type functions that diverge between humans and mice. Rare cystic fibrosis transmembrane conductance regulator (CFTR)-rich pulmonary ionocytes exist throughout the cartilaginous airways of humans1,2, but limited presence and divergent biology in the proximal trachea of mice has prevented the use of traditional transgenic models to elucidate ionocyte functions in the airway. Here we describe the creation and use of conditional genetic ferret models to dissect pulmonary ionocyte biology and function by enabling ionocyte lineage tracing (FOXI1-CreERT2::ROSA-TG), ionocyte ablation (FOXI1-KO) and ionocyte-specific deletion of CFTR (FOXI1-CreERT2::CFTRL/L). By comparing these models with cystic fibrosis ferrets3,4, we demonstrate that ionocytes control airway surface liquid absorption, secretion, pH and mucus viscosity-leading to reduced airway surface liquid volume and impaired mucociliary clearance in cystic fibrosis, FOXI1-KO and FOXI1-CreERT2::CFTRL/L ferrets. These processes are regulated by CFTR-dependent ionocyte transport of Cl- and HCO3-. Single-cell transcriptomics and in vivo lineage tracing revealed three subtypes of pulmonary ionocytes and a FOXI1-lineage common rare cell progenitor for ionocytes, tuft cells and neuroendocrine cells during airway development. Thus, rare pulmonary ionocytes perform critical CFTR-dependent functions in the proximal airway that are hallmark features of cystic fibrosis airway disease. These studies provide a road map for using conditional genetics in the first non-rodent mammal to address gene function, cell biology and disease processes that have greater evolutionary conservation between humans and ferrets.

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.

Despite advances in high-dimensional cellular analysis, the molecular profiling of dynamic behaviors of cells in their native environment remains a major challenge. We present a method that allows us to couple the physiological behaviors of cells in an intact murine tissue to deep molecular profiling of individual cells. This method enabled us to establish a novel molecular signature for a striking migratory cellular behavior following injury in murine airways.

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.

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.

The two populations of cholinergic amacrine cells in the inner nuclear layer (INL) and the ganglion cell layer (GCL) differ in their spatial organization in the mouse retina, but the basis for this difference is not understood. The present investigation examined this issue in six strains of mice that differ in their number of cholinergic cells, addressing how the regularity, packing, and spacing of these cells varies as a function of strain, layer, and density. The number of cholinergic cells was lower in the GCL than in the INL in all six strains. The nearest neighbor and Voronoi domain regularity indexes as well as the packing factor were each consistently lower for the GCL. While these regularity indexes and the packing factor were largely stable across variation in density, the effective radius was inversely related to density for both the GCL and INL, being smaller and more variable in the GCL. Consequently, despite the lower densities in the GCL, neighboring cells were more likely to be positioned closer to one another than in the higher-density INL, thereby reducing regularity and packing. This difference in the spatial organization of cholinergic cells may be due to the cells in the GCL having been passively displaced by fascicles of optic axons and an expanding retinal vasculature during development. In support of this interpretation, we show such displacement of cholinergic somata relative to their dendritic stalks and a decline in packing efficiency and regularity during postnatal development that is more severe for the GCL.

Little is known about how metabolites couple tissue-specific stem cell function with physiology. Here we show that, in the mammalian small intestine, the expression of Hmgcs2 (3-hydroxy-3-methylglutaryl-CoA synthetase 2), the gene encoding the rate-limiting enzyme in the production of ketone bodies, including beta-hydroxybutyrate (βOHB), distinguishes self-renewing Lgr5+ stem cells (ISCs) from differentiated cell types. Hmgcs2 loss depletes βOHB levels in Lgr5+ ISCs and skews their differentiation toward secretory cell fates, which can be rescued by exogenous βOHB and class I histone deacetylase (HDAC) inhibitor treatment. Mechanistically, βOHB acts by inhibiting HDACs to reinforce Notch signaling, instructing ISC self-renewal and lineage decisions. Notably, although a high-fat ketogenic diet elevates ISC function and post-injury regeneration through βOHB-mediated Notch signaling, a glucose-supplemented diet has the opposite effects. These findings reveal how control of βOHB-activated signaling in ISCs by diet helps to fine-tune stem cell adaptation in homeostasis and injury.

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.

Transgenic mouse lines are routinely employed to label and manipulate distinct cell types. The transgene generally comprises cell-type specific regulatory elements linked to a cDNA encoding a reporter or other protein. However, off-target expression seemingly unrelated to the regulatory elements in the transgene is often observed, it is sometimes suspected to reflect influences related to the site of transgene integration in the genome. To test this hypothesis, we used a proximity ligation-based method, Targeted Locus Amplification (TLA), to map the insertion sites of three well-characterized transgenes that appeared to exhibit insertion site-dependent expression in retina. The nearest endogenous genes to transgenes HB9-GFP, Mito-P, and TYW3 are Cdh6, Fat4 and Khdrbs2, respectively. For two lines, we demonstrate that expression reflects that of the closest endogenous gene (Fat4 and Cdh6), even though the distance between transgene and endogenous gene is 550 and 680 kb, respectively. In all three lines, the transgenes decrease expression of the neighboring endogenous genes. In each case, the affected endogenous gene was expressed in at least some of the cell types that the transgenic line has been used to mark and study. These results provide insights into the effects of transgenes and endogenous genes on each other's expression, demonstrate that mapping insertion site is valuable for interpreting results obtained with transgenic lines, and indicate that TLA is a reliable method for integration site discovery.