The Human and Mouse Enteric Nervous System at Single-Cell Resolution.

The enteric nervous system (ENS) coordinates diverse functions in the intestine but has eluded comprehensive molecular characterization because of the rarity and diversity of cells. Here we develop two methods to profile the ENS of adult mice and humans at single-cell resolution: RAISIN RNA-seq for profiling intact nuclei with ribosome-bound mRNA and MIRACL-seq for label-free enrichment of rare cell types by droplet-based profiling. The 1,187,535 nuclei in our mouse atlas include 5,068 neurons from the ileum and colon, revealing extraordinary neuron diversity. We highlight circadian expression changes in enteric neurons, show that disease-related genes are dysregulated with aging, and identify differences between the ileum and proximal/distal colon. In humans, we profile 436,202 nuclei, recovering 1,445 neurons, and identify conserved and species-specific transcriptional programs and putative neuro-epithelial, neuro-stromal, and neuro-immune interactions. The human ENS expresses risk genes for neuropathic, inflammatory, and extra-intestinal diseases, suggesting neuronal contributions to disease.

snRNA-seq reveals a subpopulation of adipocytes that regulates thermogenesis.

Adipose tissue is usually classified on the basis of its function as white, brown or beige (brite)1. It is an important regulator of systemic metabolism, as shown by the fact that dysfunctional adipose tissue in obesity leads to a variety of secondary metabolic complications2,3. In addition, adipose tissue functions as a signalling hub that regulates systemic metabolism through paracrine and endocrine signals4. Here we use single-nucleus RNA-sequencing (snRNA-seq) analysis in mice and humans to characterize adipocyte heterogeneity. We identify a rare subpopulation of adipocytes in mice that increases in abundance at higher temperatures, and we show that this subpopulation regulates the activity of neighbouring adipocytes through acetate-mediated modulation of their thermogenic capacity. Human adipose tissue contains higher numbers of cells of this subpopulation, which could explain the lower thermogenic activity of human compared to mouse adipose tissue and suggests that targeting this pathway could be used to restore thermogenic activity.

The brain parenchyma has a type I interferon response that can limit virus spread.

The brain has a tightly regulated environment that protects neurons and limits inflammation, designated "immune privilege." However, there is not an absolute lack of an immune response. We tested the ability of the brain to initiate an innate immune response to a virus, which was directly injected into the brain parenchyma, and to determine whether this response could limit viral spread. We injected vesicular stomatitis virus (VSV), a transsynaptic tracer, or naturally occurring VSV-derived defective interfering particles (DIPs), into the caudate-putamen (CP) and scored for an innate immune response and inhibition of virus spread. We found that the brain parenchyma has a functional type I interferon (IFN) response that can limit VSV spread at both the inoculation site and among synaptically connected neurons. Furthermore, we characterized the response of microglia to VSV infection and found that infected microglia produced type I IFN and uninfected microglia induced an innate immune response following virus injection.

Preferential Budding of Vesicular Stomatitis Virus from the Basolateral Surface of Polarized Epithelial Cells Is Not Solely Directed by Matrix Protein or Glycoprotein.

Vesicular stomatitis virus has been shown to bud basolaterally, and the matrix protein, but not glycoprotein, was proposed to mediate this asymmetry. Using polarized T84 monolayers, we demonstrate that no single viral protein is sufficient for polarized budding. Particles are released from the apical and basolateral surfaces and are indistinguishable, indicating that there is no apical assembly defect. We propose that aspects of host cell polarity create a more efficient budding process at the basolateral surface.

Lassa-vesicular stomatitis chimeric virus safely destroys brain tumors.

UNLABELLED: High-grade tumors in the brain are among the deadliest of cancers. Here, we took a promising oncolytic virus, vesicular stomatitis virus (VSV), and tested the hypothesis that the neurotoxicity associated with the virus could be eliminated without blocking its oncolytic potential in the brain by replacing the neurotropic VSV glycoprotein with the glycoprotein from one of five different viruses, including Ebola virus, Marburg virus, lymphocytic choriomeningitis virus (LCMV), rabies virus, and Lassa virus. Based on in vitro infections of normal and tumor cells, we selected two viruses to test in vivo. Wild-type VSV was lethal when injected directly into the brain. In contrast, a novel chimeric virus (VSV-LASV-GPC) containing genes from both the Lassa virus glycoprotein precursor (GPC) and VSV showed no adverse actions within or outside the brain and targeted and completely destroyed brain cancer, including high-grade glioblastoma and melanoma, even in metastatic cancer models. When mice had two brain tumors, intratumoral VSV-LASV-GPC injection in one tumor (glioma or melanoma) led to complete tumor destruction; importantly, the virus moved contralaterally within the brain to selectively infect the second noninjected tumor. A chimeric virus combining VSV genes with the gene coding for the Ebola virus glycoprotein was safe in the brain and also selectively targeted brain tumors but was substantially less effective in destroying brain tumors and prolonging survival of tumor-bearing mice. A tropism for multiple cancer types combined with an exquisite tumor specificity opens a new door to widespread application of VSV-LASV-GPC as a safe and efficacious oncolytic chimeric virus within the brain. IMPORTANCE: Many viruses have been tested for their ability to target and kill cancer cells. Vesicular stomatitis virus (VSV) has shown substantial promise, but a key problem is that if it enters the brain, it can generate adverse neurologic consequences, including death. We tested a series of chimeric viruses containing genes coding for VSV, together with a gene coding for the glycoprotein from other viruses, including Ebola virus, Lassa virus, LCMV, rabies virus, and Marburg virus, which was substituted for the VSV glycoprotein gene. Ebola and Lassa chimeric viruses were safe in the brain and targeted brain tumors. Lassa-VSV was particularly effective, showed no adverse side effects even when injected directly into the brain, and targeted and destroyed two different types of deadly brain cancer, including glioblastoma and melanoma.

Tissue and cellular spatiotemporal dynamics in colon aging.

Tissue structure and molecular circuitry in the colon can be profoundly impacted by systemic age-related effects, but many of the underlying molecular cues remain unclear. Here, we built a cellular and spatial atlas of the colon across three anatomical regions and 11 age groups, encompassing ~1,500 mouse gut tissues profiled by spatial transcriptomics and ~400,000 single nucleus RNA-seq profiles. We developed a new computational framework, cSplotch, which learns a hierarchical Bayesian model of spatially resolved cellular expression associated with age, tissue region, and sex, by leveraging histological features to share information across tissue samples and data modalities. Using this model, we identified cellular and molecular gradients along the adult colonic tract and across the main crypt axis, and multicellular programs associated with aging in the large intestine. Our multi-modal framework for the investigation of cell and tissue organization can aid in the understanding of cellular roles in tissue-level pathology.

Assessment of spatial transcriptomics for oncology discovery.

Tumor heterogeneity is a major challenge for oncology drug discovery and development. Understanding of the spatial tumor landscape is key to identifying new targets and impactful model systems. Here, we test the utility of spatial transcriptomics (ST) for oncology discovery by profiling 40 tissue sections and 80,024 capture spots across a diverse set of tissue types, sample formats, and RNA capture chemistries. We verify the accuracy and fidelity of ST by leveraging matched pathology analysis, which provides a ground truth for tissue section composition. We then use spatial data to demonstrate the capture of key tumor depth features, identifying hypoxia, necrosis, vasculature, and extracellular matrix variation. We also leverage spatial context to identify relative cell-type locations showing the anti-correlation of tumor and immune cells in syngeneic cancer models. Lastly, we demonstrate target identification approaches in clinical pancreatic adenocarcinoma samples, highlighting tumor intrinsic biomarkers and paracrine signaling.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Author Correction: Nuclei multiplexing with barcoded antibodies for single-nucleus genomics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors.

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

Nuclei multiplexing with barcoded antibodies for single-nucleus genomics.

Single-nucleus RNA-seq (snRNA-seq) enables the interrogation of cellular states in complex tissues that are challenging to dissociate or are frozen, and opens the way to human genetics studies, clinical trials, and precise cell atlases of large organs. However, such applications are currently limited by batch effects, processing, and costs. Here, we present an approach for multiplexing snRNA-seq, using sample-barcoded antibodies to uniquely label nuclei from distinct samples. Comparing human brain cortex samples profiled with or without hashing antibodies, we demonstrate that nucleus hashing does not significantly alter recovered profiles. We develop DemuxEM, a computational tool that detects inter-sample multiplets and assigns singlets to their sample of origin, and validate its accuracy using sex-specific gene expression, species-mixing and natural genetic variation. Our approach will facilitate tissue atlases of isogenic model organisms or from multiple biopsies or longitudinal samples of one donor, and large-scale perturbation screens.

Glycolytic reliance promotes anabolism in photoreceptors.

Vertebrate photoreceptors are among the most metabolically active cells, exhibiting a high rate of ATP consumption. This is coupled with a high anabolic demand, necessitated by the diurnal turnover of a specialized membrane-rich organelle, the outer segment, which is the primary site of phototransduction. How photoreceptors balance their catabolic and anabolic demands is poorly understood. Here, we show that rod photoreceptors in mice rely on glycolysis for their outer segment biogenesis. Genetic perturbations targeting allostery or key regulatory nodes in the glycolytic pathway impacted the size of the outer segments. Fibroblast growth factor signaling was found to regulate glycolysis, with antagonism of this pathway resulting in anabolic deficits. These data demonstrate the cell autonomous role of the glycolytic pathway in outer segment maintenance and provide evidence that aerobic glycolysis is part of a metabolic program that supports the biosynthetic needs of a normal neuronal cell type.

Detection and manipulation of live antigen-expressing cells using conditionally stable nanobodies.

The ability to detect and/or manipulate specific cell populations based upon the presence of intracellular protein epitopes would enable many types of studies and applications. Protein binders such as nanobodies (Nbs) can target untagged proteins (antigens) in the intracellular environment. However, genetically expressed protein binders are stable regardless of antigen expression, complicating their use for applications that require cell-specificity. Here, we created a conditional system in which the stability of an Nb depends upon an antigen of interest. We identified Nb framework mutations that can be used to rapidly create destabilized Nbs. Fusion of destabilized Nbs to various proteins enabled applications in living cells, such as optogenetic control of neural activity in specific cell types in the mouse brain, and detection of HIV-infected human cells by flow cytometry. These approaches are generalizable to other protein binders, and enable the rapid generation of single-polypeptide sensors and effectors active in cells expressing specific intracellular epitopes.

Cell type-specific manipulation with GFP-dependent Cre recombinase.

There are many transgenic GFP reporter lines that allow the visualization of specific populations of cells. Using such lines for functional studies requires a method that transforms GFP into a molecule that enables genetic manipulation. We developed a method that exploits GFP for gene manipulation, Cre recombinase dependent on GFP (CRE-DOG), a split component system that uses GFP and its derivatives to directly induce Cre/loxP recombination. Using plasmid electroporation and AAV viral vectors, we delivered CRE-DOG to multiple GFP mouse lines, which led to effective recombination selectively in GFP-labeled cells. Furthermore, CRE-DOG enabled optogenetic control of these neurons. Beyond providing a new set of tools for manipulation of gene expression selectively in GFP(+) cells, we found that GFP can be used to reconstitute the activity of a protein not known to have a modular structure, suggesting that this strategy might be applicable to a wide range of proteins.

Vesicular stomatitis virus enables gene transfer and transsynaptic tracing in a wide range of organisms.

Current limitations in technology have prevented an extensive analysis of the connections among neurons, particularly within nonmammalian organisms. We developed a transsynaptic viral tracer originally for use in mice, and then tested its utility in a broader range of organisms. By engineering the vesicular stomatitis virus (VSV) to encode a fluorophore and either the rabies virus glycoprotein (RABV-G) or its own glycoprotein (VSV-G), we created viruses that can transsynaptically label neuronal circuits in either the retrograde or anterograde direction, respectively. The vectors were investigated for their utility as polysynaptic tracers of chicken and zebrafish visual pathways. They showed patterns of connectivity consistent with previously characterized visual system connections, and revealed several potentially novel connections. Further, these vectors were shown to infect neurons in several other vertebrates, including Old and New World monkeys, seahorses, axolotls, and Xenopus. They were also shown to infect two invertebrates, Drosophila melanogaster, and the box jellyfish, Tripedalia cystophora, a species previously intractable for gene transfer, although no clear evidence of transsynaptic spread was observed in these species. These vectors provide a starting point for transsynaptic tracing in most vertebrates, and are also excellent candidates for gene transfer in organisms that have been refractory to other methods.

Human CFEOM1 mutations attenuate KIF21A autoinhibition and cause oculomotor axon stalling.

The ocular motility disorder "Congenital fibrosis of the extraocular muscles type 1" (CFEOM1) results from heterozygous mutations altering the motor and third coiled-coil stalk of the anterograde kinesin, KIF21A. We demonstrate that Kif21a knockin mice harboring the most common human mutation develop CFEOM. The developing axons of the oculomotor nerve's superior division stall in the proximal nerve; the growth cones enlarge, extend excessive filopodia, and assume random trajectories. Inferior division axons reach the orbit but branch ectopically. We establish a gain-of-function mechanism and find that human motor or stalk mutations attenuate Kif21a autoinhibition, providing in vivo evidence for mammalian kinesin autoregulation. We identify Map1b as a Kif21a-interacting protein and report that Map1b⁻/⁻ mice develop CFEOM. The interaction between Kif21a and Map1b is likely to play a critical role in the pathogenesis of CFEOM1 and highlights a selective vulnerability of the developing oculomotor nerve to perturbations of the axon cytoskeleton.