The transcriptional state of a cell reflects a variety of biological factors, from cell-type-specific features to transient processes such as the cell cycle, all of which may be of interest. However, identifying such aspects from noisy single-cell RNA-seq data remains challenging. We developed pathway and gene set overdispersion analysis (PAGODA) to resolve multiple, potentially overlapping aspects of transcriptional heterogeneity by testing gene sets for coordinated variability among measured cells.

N-methyl-D-aspartate receptors (NMDARs) in all hippocampal areas play an essential role in distinct processes of memory formation as well as in sustaining cell survival of postnatally generated neurons in the dentate gyrus (DG). In contrast to the beneficial effects, over-activation of NMDARs has been implicated in many acute and chronic neurological diseases, reason why therapeutic approaches and clinical trials involving receptor blockade have been envisaged for decades. Here we employed genetically engineered mice to study the long-term effect of NMDAR ablation on selective hippocampal neuronal populations. Ablation of either GluN1 or GluN2B causes degeneration of the DG. The neuronal demise affects mature neurons specifically in the dorsal DG and is NMDAR subunit-dependent. Most importantly, the degenerative process exacerbates with increasing age of the animals. These results lead us to conclude that mature granule cells in the dorsal DG undergo neurodegeneration following NMDAR ablation in aged mouse. Thus, caution needs to be exerted when considering long-term administration of NMDAR antagonists for therapeutic purposes.

Alternative splicing is prevalent in the mammalian brain. To interrogate the functional role of alternative splicing in neural development, we analyzed purified neural progenitor cells (NPCs) and neurons from developing cerebral cortices, revealing hundreds of differentially spliced exons that preferentially alter key protein domains-especially in cytoskeletal proteins-and can harbor disease-causing mutations. We show that Ptbp1 and Rbfox proteins antagonistically govern the NPC-to-neuron transition by regulating neuron-specific exons. Whereas Ptbp1 maintains apical progenitors partly through suppressing a poison exon of Flna in NPCs, Rbfox proteins promote neuronal differentiation by switching Ninein from a centrosomal splice form in NPCs to a non-centrosomal isoform in neurons. We further uncover an intronic human mutation within a PTBP1-binding site that disrupts normal skipping of the FLNA poison exon in NPCs and causes a brain-specific malformation. Our study indicates that dynamic control of alternative splicing governs cell fate in cerebral cortical development.

Neurogenesis in the adult brain, a powerful mechanism for neuronal plasticity and brain repair, is altered by aging and pathological conditions, including metabolic disorders. The search for mechanisms and therapeutic solutions to alter neurogenesis requires understanding of cell kinetics within neurogenic niches using a high-throughput quantitative approach. The challenge is in the dynamic nature of the process and multiple cell types involved, each having several potential modes of division or cell fate. Here we show that cell kinetics can be revealed through a combination of the BrdU/EdU pulse-chase, based on the circadian pattern of DNA replication, and a differential equations model that describes time-dependent cell densities. The model is validated through the analysis of cell kinetics in the cerebellar neurogenic niche of normal young adult male zebrafish, with cells quantified in 2D (sections), and with neuronal fate and reactivation of stem cells confirmed in 3D whole-brain images (CLARITY). We then reveal complex alterations in cell kinetics associated with accelerated aging due to chronic high caloric intake. Low activity of neuronal stem cells in this condition persists 2 months after reverting to normal diet, and is accompanied by overproduction of transient amplifying cells, their accelerated cell death, and slow migration of postmitotic progeny. This combined experimental and mathematical approach should allow for relatively high-throughput analysis of early signs of pathological and age-related changes in neurogenesis, evaluation of specific therapeutic targets, and drug efficacy.SIGNIFICANCE STATEMENT Understanding normal cell kinetics of adult neurogenesis and the type of cells affected by a pathological process is needed to develop effective prophylactic and therapeutic measures directed at specific cell targets. Complex time-dependent mechanisms involved in the kinetics of multiple cell types require a combination of experimental and mathematical modeling approaches. This study demonstrates such a combined approach by comparing normal neurogenesis with that altered by diet-induced accelerated aging in adult zebrafish.

Recent single-cell RNA-seq protocols based on droplet microfluidics use massively multiplexed barcoding to enable simultaneous measurements of transcriptomes for thousands of individual cells. The increasing complexity of such data creates challenges for subsequent computational processing and troubleshooting of these experiments, with few software options currently available. Here, we describe a flexible pipeline for processing droplet-based transcriptome data that implements barcode corrections, classification of cell quality, and diagnostic information about the droplet libraries. We introduce advanced methods for correcting composition bias and sequencing errors affecting cellular and molecular barcodes to provide more accurate estimates of molecular counts in individual cells.

Engineered nanoparticles are smaller than 100 nm and designed to improve or creating even new physico-chemical properties. Consequently, toxicological properties of materials may change as size reaches the nm size-range. We examined outcomes related to the central nervous system in the offspring following maternal inhalation exposure to nanosized carbon black particles (Printex 90).

Chromatin immunoprecipitation (ChIP) followed by high-throughput DNA sequencing (ChIP-seq) has become a valuable and widely used approach for mapping the genomic location of transcription-factor binding and histone modifications in living cells. Despite its widespread use, there are considerable differences in how these experiments are conducted, how the results are scored and evaluated for quality, and how the data and metadata are archived for public use. These practices affect the quality and utility of any global ChIP experiment. Through our experience in performing ChIP-seq experiments, the ENCODE and modENCODE consortia have developed a set of working standards and guidelines for ChIP experiments that are updated routinely. The current guidelines address antibody validation, experimental replication, sequencing depth, data and metadata reporting, and data quality assessment. We discuss how ChIP quality, assessed in these ways, affects different uses of ChIP-seq data. All data sets used in the analysis have been deposited for public viewing and downloading at the ENCODE (http://encodeproject.org/ENCODE/) and modENCODE (http://www.modencode.org/) portals.

Chromatin insulator elements and associated proteins have been proposed to partition eukaryotic genomes into sets of independently regulated domains. Here we test this hypothesis by quantitative genome-wide analysis of insulator protein binding to Drosophila chromatin. We find distinct combinatorial binding of insulator proteins to different classes of sites and uncover a novel type of insulator element that binds CP190 but not any other known insulator proteins. Functional characterization of different classes of binding sites indicates that only a small fraction act as robust insulators in standard enhancer-blocking assays. We show that insulators restrict the spreading of the H3K27me3 mark but only at a small number of Polycomb target regions and only to prevent repressive histone methylation within adjacent genes that are already transcriptionally inactive. RNAi knockdown of insulator proteins in cultured cells does not lead to major alterations in genome expression. Taken together, these observations argue against the concept of a genome partitioned by specialized boundary elements and suggest that insulators are reserved for specific regulation of selected genes.

Gating properties and surface trafficking of AMPA receptors (AMPARs) are modulated by auxiliary subunits. Here we studied the function of coexpressed auxiliary subunits belonging to two different classes. We focused on TARP γ-8 and CKAMP44 in dentate gyrus (DG) granule cells, since both subunits are highly expressed in this cell type. TARP γ-8 and CKAMP44 decrease the rate of deactivation but have an opposing influence on receptor desensitization, which accounts for their differential modulation of synaptic short-term plasticity. Furthermore, long-term plasticity (LTP) requires TARP γ-8 but not CKAMP44. The coexpression of both auxiliary subunits is necessary for the efficient targeting of AMPARs to the cell surface of DG granule cells. Finally, electrophysiological and biochemical evidence support the notion that CKAMP44 and TARP γ-8 can be contained in the same AMPAR complex.

Using 40 known human-specific LTR sequences, we have derived a consensus sequence for an evolutionary young HERV-K (HML-2) LTR family, which was named the HS family. In the human genome the HS family is represented by approximately 150-160 LTR sequences, 90% of them being human-specific (hs). The family can be subdivided into two subfamilies differing in five linked nucleotide substitutions: HS-a and HS-b of 5.8 and 10.3 Myr evolutionary ages, respectively. The HS-b subfamily members were transpositionally active both before the divergence of the human and chimpanzee ancestor lineages and after it in both lineages. The HS-a subfamily comprises only hs LTRs. These and other data strongly suggest that at least three "master genes" of HERV-K (HML-2) LTRs were active in the human ancestor lineage after the human-chimpanzee divergence. We also found hs HERV-K (HML-2) LTRs integrations in introns of 12 human genes and identified 13 new hs HERV-K (HML-2) LTRs.

Ca(2+)-permeable AMPA receptors are densely expressed in the spinal dorsal horn, but their functional significance in pain processing is not understood. By disrupting the genes encoding GluR-A or GluR-B, we generated mice exhibiting increased or decreased numbers of Ca(2+)-permeable AMPA receptors, respectively. Here, we demonstrate that AMPA receptors are critical determinants of nociceptive plasticity and inflammatory pain. A reduction in the number of Ca(2+)-permeable AMPA receptors and density of AMPA channel currents in spinal neurons of GluR-A-deficient mice is accompanied by a loss of nociceptive plasticity in vitro and a reduction in acute inflammatory hyperalgesia in vivo. In contrast, an increase in spinal Ca(2+)-permeable AMPA receptors in GluR-B-deficient mice facilitated nociceptive plasticity and enhanced long-lasting inflammatory hyperalgesia. Thus, AMPA receptors are not mere determinants of fast synaptic transmission underlying basal pain sensitivity as previously thought, but are critically involved in activity-dependent changes in synaptic processing of nociceptive inputs.

In many vertebrates, postnatally generated neurons often migrate long distances to reach their final destination, where they help shape local circuit activity. Concerted action of extrinsic stimuli is required to regulate long-distance migration. Some migratory principles are evolutionarily conserved, whereas others are species and cell type specific. Here we identified a serotonergic mechanism that governs migration of postnatally generated neurons in the mouse brain. Serotonergic axons originating from the raphe nuclei exhibit a conspicuous alignment with subventricular zone-derived neuroblasts. Optogenetic axonal activation provides functional evidence for serotonergic modulation of neuroblast migration. Furthermore, we show that the underlying mechanism involves serotonin receptor 3A (5HT3A)-mediated calcium influx. Thus, 5HT3A receptor deletion in neuroblasts impaired speed and directionality of migration and abolished calcium spikes. We speculate that serotonergic modulation of postnatally generated neuroblast migration is evolutionarily conserved as indicated by the presence of serotonergic axons in migratory paths in other vertebrates.

To understand the function of cortical circuits, it is necessary to catalog their cellular diversity. Past attempts to do so using anatomical, physiological or molecular features of cortical cells have not resulted in a unified taxonomy of neuronal or glial cell types, partly due to limited data. Single-cell transcriptomics is enabling, for the first time, systematic high-throughput measurements of cortical cells and generation of datasets that hold the promise of being complete, accurate and permanent. Statistical analyses of these data reveal clusters that often correspond to cell types previously defined by morphological or physiological criteria and that appear conserved across cortical areas and species. To capitalize on these new methods, we propose the adoption of a transcriptome-based taxonomy of cell types for mammalian neocortex. This classification should be hierarchical and use a standardized nomenclature. It should be based on a probabilistic definition of a cell type and incorporate data from different approaches, developmental stages and species. A community-based classification and data aggregation model, such as a knowledge graph, could provide a common foundation for the study of cortical circuits. This community-based classification, nomenclature and data aggregation could serve as an example for cell type atlases in other parts of the body.

The peripheral nervous system responds to a wide variety of sensory stimuli, a process that requires great neuronal diversity. These diverse neurons are closely associated with glial cells originating from the neural crest. However, the molecular nature and diversity among peripheral glia are not understood. Here, we used single-cell RNA sequencing to profile developing and mature glia from somatosensory dorsal root ganglia and auditory spiral ganglia. We found that glial precursors (GPs) in these two systems differ in their transcriptional profiles. Despite their unique features, somatosensory and auditory GPs undergo convergent differentiation to generate molecularly uniform myelinating and non-myelinating Schwann cells. By contrast, somatosensory and auditory satellite glial cells retain system-specific features. Lastly, we identified a glial signature gene set, providing new insights into commonalities among glia across the nervous system. This survey of gene expression in peripheral glia constitutes a resource for understanding functions of glia across different sensory modalities.

Stellate cells are principal neurons in the entorhinal cortex that contribute to spatial processing. They also play a role in the context of Alzheimer's disease as they accumulate Amyloid beta early in the disease. Producing human stellate cells from pluripotent stem cells would allow researchers to study early mechanisms of Alzheimer's disease, however, no protocols currently exist for producing such cells. In order to develop novel stem cell protocols, we characterize at high resolution the development of the porcine medial entorhinal cortex by tracing neuronal and glial subtypes from mid-gestation to the adult brain to identify the transcriptomic profile of progenitor and adult stellate cells. Importantly, we could confirm the robustness of our data by extracting developmental factors from the identified intermediate stellate cell cluster and implemented these factors to generate putative intermediate stellate cells from human induced pluripotent stem cells. Six transcription factors identified from the stellate cell cluster including RUNX1T1, SOX5, FOXP1, MEF2C, TCF4, EYA2 were overexpressed using a forward programming approach to produce neurons expressing a unique combination of RELN, SATB2, LEF1 and BCL11B observed in stellate cells. Further analyses of the individual transcription factors led to the discovery that FOXP1 is critical in the reprogramming process and omission of RUNX1T1 and EYA2 enhances neuron conversion. Our findings contribute not only to the profiling of cell types within the developing and adult brain's medial entorhinal cortex but also provides proof-of-concept for using scRNAseq data to produce entorhinal intermediate stellate cells from human pluripotent stem cells in-vitro.

Single-cell RNA sequencing is often applied in study designs that include multiple individuals, conditions or tissues. To identify recurrent cell subpopulations in such heterogeneous collections, we developed Conos, an approach that relies on multiple plausible inter-sample mappings to construct a global graph connecting all measured cells. The graph enables identification of recurrent cell clusters and propagation of information between datasets in multi-sample or atlas-scale collections.

Neuroserpin is a serine-protease inhibitor mainly expressed in the CNS and involved in the inhibition of the proteolytic cascade. Animal models confirmed its neuroprotective role in perinatal hypoxia-ischaemia and adult stroke. Although neuroserpin may be a potential therapeutic target in the treatment of the aforementioned conditions, there is still no information in the literature on its distribution during human brain development. The present study provides a detailed description of the changing spatiotemporal patterns of neuroserpin focusing on physiological human brain development. Five stages were distinguished within our examined age range which spanned from the 7th gestational week until adulthood. In particular, subplate and deep cortical plate neurons were identified as the main sources of neuroserpin production between the 25th gestational week and the first postnatal month. Our immunohistochemical findings were substantiated by single cell RNA sequencing data showing specific neuronal and glial cell types expressing neuroserpin. The characterization of neuroserpin expression during physiological human brain development is essential for forthcoming studies which will explore its involvement in pathological conditions, such as perinatal hypoxia-ischaemia and adult stroke in human.

NUT midline carcinoma (NMC), a subtype of squamous cell cancer, is one of the most aggressive human solid malignancies known. NMC is driven by the creation of a translocation oncoprotein, BRD4-NUT, which blocks differentiation and drives growth of NMC cells. BRD4-NUT forms distinctive nuclear foci in patient tumors, which we found correlate with ∼100 unprecedented, hyperacetylated expanses of chromatin that reach up to 2 Mb in size. These "megadomains" appear to be the result of aberrant, feed-forward loops of acetylation and binding of acetylated histones that drive transcription of underlying DNA in NMC patient cells and naïve cells induced to express BRD4-NUT. Megadomain locations are typically cell lineage-specific; however, the cMYC and TP63 regions are targeted in all NMCs tested and play functional roles in tumor growth. Megadomains appear to originate from select pre-existing enhancers that progressively broaden but are ultimately delimited by topologically associating domain (TAD) boundaries. Therefore, our findings establish a basis for understanding the powerful role played by large-scale chromatin organization in normal and aberrant lineage-specific gene transcription.

Pain alters opioid reinforcement, presumably via neuroadaptations within ascending pain pathways interacting with the limbic system. Nerve injury increases expression of glutamate receptors and their associated Homer scaffolding proteins throughout the pain processing pathway. Homer proteins, and their associated glutamate receptors, regulate behavioral sensitivity to various addictive drugs. Thus, we investigated a potential role for Homers in the interactions between pain and drug reward in mice. Chronic constriction injury (CCI) of the sciatic nerve elevated Homer1b/c and/or Homer2a/b expression within all mesolimbic structures examined and for the most part, the Homer increases coincided with elevated mGluR5, GluN2A/B, and the activational state of various down-stream kinases. Behaviorally, CCI mice showed pain hypersensitivity and a conditioned place-aversion (CPA) at a low heroin dose that supported conditioned place-preference (CPP) in naïve controls. Null mutations of Homer1a, Homer1, and Homer2, as well as transgenic disruption of mGluR5-Homer interactions, either attenuated or completely blocked low-dose heroin CPP, and none of the CCI mutant strains exhibited heroin-induced CPA. However, heroin CPP did not depend upon full Homer1c expression within the nucleus accumbens (NAC), as CPP occurred in controls infused locally with small hairpin RNA-Homer1c, although intra-NAC and/or intrathecal cDNA-Homer1c, -Homer1a, and -Homer2b infusions (to best mimic CCI's effects) were sufficient to blunt heroin CPP in uninjured mice. However, arguing against a simple role for CCI-induced increases in either spinal or NAC Homer expression for heroin CPA, cDNA infusion of our various cDNA constructs either did not affect (intrathecal) or attenuated (NAC) heroin CPA. Together, these data implicate increases in glutamate receptor/Homer/kinase activity within limbic structures, perhaps outside the NAC, as possibly critical for switching the incentive motivational properties of heroin following nerve injury, which has relevance for opioid psychopharmacology in individuals suffering from neuropathic pain.

Significant heterogeneities in gene expression among individual cells are typically interrogated using single whole cell approaches. However, tissues that have highly interconnected processes, such as in the brain, present unique challenges. Single-nucleus RNA sequencing (SNS) has emerged as an alternative method of assessing a cell's transcriptome through the use of isolated nuclei. However, studies directly comparing expression data between nuclei and whole cells are lacking. Here, we have characterized nuclear and whole cell transcriptomes in mouse single neurons and provided a normalization strategy to reduce method-specific differences related to the length of genic regions. We confirmed a high concordance between nuclear and whole cell transcriptomes in the expression of cell type and metabolic modeling markers, but less so for a subset of genes associated with mitochondrial respiration. Therefore, our results indicate that single-nucleus transcriptome sequencing provides an effective means to profile cell type expression dynamics in previously inaccessible tissues.